Constrained MD using andersen thermostat at 300K
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 2
- Joined: Wed Feb 11, 2015 6:47 am
- License Nr.: #5-1160
Constrained MD using andersen thermostat at 300K
I have repeated my constrained Md calculation of a single water molecule using andersen thermostat and keeping bond length and bond angle fixed at 300K using smaller time step. Here also after initial movement O atom moves in a strange way. I am attaching my input files. Please help me to solve the problem.
POSCAR
H2O
0.52918 ! scaling parameter
20 0 0
0 20 0
0 0 20
O H
1 2
cart
10.00 10.00 10.00
11.10 8.57 10.00
11.10 11.43 10.00
INCAR
ALGO = V
MAXMIX = 40
IBRION = 0
NSW = 10000
NBLOCK = 10
POTIM = 0.3
TEBEG = 300; TEEND = 300
ISYM = 0
MDALGO = 1
ANDERSEN PROB = 0.0
LREAL =.TRUE.
LCHARG =.FALSE.
NELMIN = 4
PREC = LOW
ISTART = 0
ISMEAR = 0; SIGMA=0.1
Displacement of O atom:
step x(ang) y(ang) z(ang) r(ang)
10 5.2921 5.2929 5.2855 9.1629
5000 5.3447 5.3065 5.2299 9.1694
10000 5.3453 5.2286 5.2840 9.1559
POSCAR
H2O
0.52918 ! scaling parameter
20 0 0
0 20 0
0 0 20
O H
1 2
cart
10.00 10.00 10.00
11.10 8.57 10.00
11.10 11.43 10.00
INCAR
ALGO = V
MAXMIX = 40
IBRION = 0
NSW = 10000
NBLOCK = 10
POTIM = 0.3
TEBEG = 300; TEEND = 300
ISYM = 0
MDALGO = 1
ANDERSEN PROB = 0.0
LREAL =.TRUE.
LCHARG =.FALSE.
NELMIN = 4
PREC = LOW
ISTART = 0
ISMEAR = 0; SIGMA=0.1
Displacement of O atom:
step x(ang) y(ang) z(ang) r(ang)
10 5.2921 5.2929 5.2855 9.1629
5000 5.3447 5.3065 5.2299 9.1694
10000 5.3453 5.2286 5.2840 9.1559
-
- Global Moderator
- Posts: 1817
- Joined: Mon Nov 18, 2019 11:00 am
Re: Constrained MD using andersen thermostat at 300K
Hi,
We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.
Best wishes,
VASP