[ase-users] Atoms' Constraints
陳佑任
b00504051 at ntu.edu.tw
Wed May 10 10:00:19 CEST 2017
Dear ASE users,
I am using ase.io.read() to read in following POSCAR in VASP format to do a normal mode analysis:
---
H O
1.00000000000000
10.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 10.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 10.0000000000000000
H O
2 1
Selective dynamics
Cartesian
5.1707077079 4.0137770615 5.2080185115 F F F
5.0954999575 5.5514905157 5.1493322874 T T T
5.5902191105 4.8193934419 5.5566036301 T T T
0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
---
Although it does not make too much sense here, I have fixed one of the H atoms hoping to obtain the eigenvalues of a partial Hessian instead of the full Hessian. After using the ase.io.read function, it shows: ('atoms:', Atoms(symbols='H2O', pbc=True, cell=[10.0, 10.0, 10.0], constraint=FixAtoms(indices=[0]))). If I understand this correctly, ASE has recognized that the first atom should be fixed. However, when I run the vibrational calculation using vib.run() and vib.summary(), it shows the following.
---
---------------------
# meV cm^-1
---------------------
0 136.1i 1097.5i
1 61.2i 493.9i
2 37.2i 299.6i
3 24.1i 194.1i
4 3.9 31.2
5 27.6 222.6
6 156.8 1264.8
7 172.2 1389.0
8 447.4 3608.8
---------------------
Zero-point energy: 0.404 eV
---
There is a problem that I cannot figure out. Since I fix one H atom, I expect that it should only show 6 values, but it gives me 9 values. In other words, the degrees of freedom of the first atoms has not been fixed. In line with that, the POSCAR files generated by ASE do feature structures that have the first atom moved from the initial position. Can anyone help me with fixing the atoms for the normal mode analysis?
Yu-Ren Chen
Institute of Atomic and Molecular Sciences, Academia Sinica
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listserv.fysik.dtu.dk/pipermail/ase-users/attachments/20170510/6d3909e5/attachment-0001.html>
More information about the ase-users
mailing list