Dear siesta users

For 20 atoms in unit cell we want to constrain all atoms moves in all
directions except for atoms number 11 and 16 that are fixes just in x and y
directions and are free to move in z direction.

W use this block GeometryConstraints in fdf input file


%block GeometryConstraints

position from -1 to -10
position from -12 to -15
position from -17 to -20
center 11 16         1.0  0.0  0.0
center 11 16         0.0  1.0  0.0

%endblock GeometryConstraints

is it true?


*Thanks in advance*

Responder a