LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Free energy
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Free energy


From: Lamm Gro <lammgro@...24...>
Date: Wed, 26 Jul 2017 15:16:46 +0900

Dear Lammps users ,

I try to calculate the free energy of water with FEP method .
( free energy of one water molecule inside of a box of it ) .
Unfortunately I just get 0.0 amount for free energy ! 
Can you please let me know what might be problem ?

This is a part of my input file : 
=================================================
pair_style  hybrid lj/cut/tip4p/long/soft 2 1 1 1 0.1546 2.0 0.5 10.0 12.0 &
          lj/cut/tip4p/long 4 3 1 1 0.1546 12.0
kspace_style pppm/tip4p 1.0e-4

pair_modify tail no

pair_coeff    1    1  lj/cut/tip4p/long/soft 0.0000  1.0000 1.0  # HW-Hw
pair_coeff    1    2  lj/cut/tip4p/long/soft 0.0000  1.0000 1.0  # HW-OW 
pair_coeff    2    2  lj/cut/tip4p/long/soft 0.1852  3.1589 1.0  # Ow Ow

pair_coeff    1    3  lj/cut/tip4p/long/soft 0.0000  1.0000 0.0  # HW-HW 
pair_coeff    1    4  lj/cut/tip4p/long/soft 0.0000  1.0000 0.0  # HW-OW 
pair_coeff    2    3  lj/cut/tip4p/long/soft 0.0000  1.0000 0.0  # OW-HW 
pair_coeff    2    4  lj/cut/tip4p/long/soft 0.1852  3.1589 0.0  # OW-OW 

pair_coeff    3    3  lj/cut/tip4p/long      0.0000  1.0000   # HW-HW 
pair_coeff    3    4  lj/cut/tip4p/long      0.0000  1.0000   # HW-OW 
pair_coeff    4    4  lj/cut/tip4p/long      0.1852  3.1589   # OW-OW 

bond_coeff    1       450    0.9572
angle_coeff   1       55     104.52

set type 1 charge   0.5564
set type 2 charge  -1.1128
set type 3 charge   0.5564
set type 4 charge  -1.1128

fix fshake all shake 0.0001 20 0 b 1 a 1
velocity   all create 300 432567 dist uniform

neighbor        2.0 bin
neigh_modify    delay 0 every 1 check yes

set type 3*4 charge 0.0

timestep 1.0

fix P1 all npt temp 300.0 300.0 100 iso 1.0 1.0 1000
run 20000000
 
reset_timestep 0

variable lambda equal ramp(0.0,1.0)
variable q3  equal 0.5564*v_lambda
variable q4  equal -1.1128*v_lambda

fix adapt all adapt/fep 5000000 &
  pair lj/cut/tip4p/long/soft lambda 3*4 1*2 v_lambda &
  atom charge 3 v_q3 &
  atom charge 4 v_q4 &
  after yes

fix pnt all print 5000000 "adapt  lambda = ${lambda} q3 = ${q3}  q4 = ${q4}"

variable dlambda equal 0.01
variable dq3  equal 0.5564*v_dlambda
variable dq4  equal -1.1128*v_dlambda

compute FEP all fep 300.0 &
  pair lj/cut/tip4p/long/soft lambda 4 2 v_dlambda &
  atom charge 3 v_dq3 &
  atom charge 4 v_dq4

fix FEP all ave/time 1000 5000 5000000 c_FEP[1] c_FEP[2] file FEP01.txt

run 100000000
=======================================================

Best Regards,
Saeed.

Attachment: data.water
Description: Binary data

Attachment: in.fep
Description: Binary data