Hello everyone,

I'm having some problems with an optimization of the dimensions of the
unit cell under certain constraints. I want to perform a relaxation of
an orthorhombic unit cell where the product of two of the unit cell
vectors (a·b) remain constant, so the area along the XY plane remain
constant.

I tried to use the constr.f file to do this, but the cell vector matrix
can be only used as an input, so I can not change the values of the unit
cell vectors from here.

If I use the "MD.ConstantVolume true" option and constrain the stress
tensor along the zz, xz, yz and xy components, so I obtain an
orthorhombic cell with fixed Z, and consequently if it has fixed
volume,  I guessed that the area of the XY plane would be also constant.
But it also doesn't keep the XY area fixed, as the constant volume
constraint is applied after the stress tensor constraints, so the Z
component is altered when the cell is expanded or contracted to keep the
volume constant.

I don't know if I'm missing something and there is actually a simple way
to keep constant the area of one of the faces of the orthorhombic cell.

Thank you

Óscar



AVISO LEGAL. Este mensaje puede contener información reservada y confidencial. 
Si usted no es el destinatario no está autorizado a copiar, reproducir o 
distribuir este mensaje ni su contenido. Si ha recibido este mensaje por error, 
le rogamos que lo notifique al remitente.
Le informamos de que sus datos personales, que puedan constar en este mensaje, serán 
tratados en calidad de responsable de tratamiento por la UNIVERSIDAD NACIONAL DE EDUCACIÓN 
A DISTANCIA (UNED) c/ Bravo Murillo, 38, 28015-MADRID-, con la finalidad de mantener el 
contacto con usted. La base jurídica que legitima este tratamiento, será su consentimiento, 
el interés legítimo o la necesidad para gestionar una relación contractual o similar. En 
cualquier momento podrá ejercer sus derechos de acceso, rectificación, supresión, 
oposición, limitación al tratamiento o portabilidad de los datos, ante la UNED, 
Departamento de Política Jurídica de Seguridad de la 
Información<https://www.uned.es/dpj>, o a través de la Sede 
electrónica<https://sede.uned.es/> de la Universidad.
Para más información visite nuestra Política de 
Privacidad<https://descargas.uned.es/publico/pdf/Politica_privacidad_UNED.pdf>.
-- 
SIESTA is supported by the Spanish Research Agency (AEI) and by the European 
H2020 MaX Centre of Excellence (http://www.max-centre.eu/)

Responder a