Problems of IR and dipole calculation
Ross, Sun
sunshou... at yahoo.com.cn
Thu Oct 28 00:36:22 UTC 2010
Hi Mathieu and Fawzi,
I tested the dipoles again and again. The problem still happened
(dipole jumps and bad wannier center), then I took Mr. Teodoro Laino's
advise to set the inputfile as follows, and it works. The wannier
centers are at reasonable position and the values of dipoles seem
right. But I still get a question, as CPMD mannual goes "TUCKERMAN
Method: - box size : molecule + 3 Ang. border AND 2*size of charge
distribution". If center molecule at 0. 0. 0., that means "molecule +
3 Ang. border" fails, right? And it still puzzles me that why setting
"CENTER_POINT 0. 0. 0." could give the right result. Any advise?
1) CENTER your system in 0. 0. 0.
&TOPOLOGY
&CENTER_COORDINATES
CENTER_POINT 0. 0. 0.
&END
&END
2) Use MT (which is translational invariant) instead of WAVELET (which
forces the system to be between 0 and L) as POISSON_SOLVER.
Regards,
Thank very much for your kind help. ^_^
More information about the CP2K-user
mailing list