| 查看: 1070 | 回复: 2 | ||
[求助]
LAMMPS能量计算值不正确
|
|
大家好，关于LAMMPS模拟中，能量相关的部分，有一些问题想要请教各位。 我所采用的in文件及data文件，是moltemplate软件中的OPLS-AA力场正丁烷范例，条件为260K，1atm。在按照范例运行之后，发现内能/焓的计算值相比于NIST的数据小50%左右。我尝试了其他的温度和压力，发现能量（U及H）的计算值要比NIST数据小很多，而且集中在4-5kcal/mol附近，变化很小。 我尝试修改了special_bonds的参数，虽然值发生变化但这应该是不正确的。因此想要请教各位是否出现过同样的问题，以及对应的解决方法。 谢谢大家！ |
回复此楼» 猜你喜欢
酰胺脱乙酰基
已经有9人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有3人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有6人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
|
这里是In文件中的运行部分 # ------------------------------- Run Section ------------------------------- dump dumpeq1 all custom 50 traj_eq1_min.lammpstrj id mol type x y z ix iy iz thermo_style custom step temp press ke etotal enthalpy thermo 50 thermo_modify norm yes velocity all create 900.0 4928459 dist gaussian # -- Equilibration: part 1: initial minimization -- minimize 1.0e-5 1.0e-7 100000 400000 undump dumpeq1 write_data system_after_eq1_min.data # -- Equilibration part 2: reorienting the molecules (NVT) -- timestep 1.0 dump dumpeq2 all custom 200 traj_eq2_reorient.lammpstrj id mol type x y z ix iy iz fix fxlan all langevin 900.0 900.0 120 48279 fix_modify fxlan energy yes fix fxnve all nve run 4000 velocity all scale 900.0 sum yes unfix fxlan unfix fxnve fix fxnvt all nvt temp 900.0 900.0 100.0 fix_modify fxnvt energy yes run 10000 undump dumpeq2 write_data system_after_eq2_reorient.data unfix fxnvt # -- equilibration part 3: Equilibrating the density (NPT) -- dump dumpeq3 all custom 200 traj_eq3_npt.lammpstrj id mol type x y z ix iy iz fix fxnpt all npt temp 900.0 260.0 100.0 iso 500.0 500.0 1000.0 drag 2.0 timestep 1.0 run 20000 unfix fxnpt fix fxnpt all npt temp 260.0 260.0 100.0 iso 500.0 1.0 1000.0 drag 2.0 timestep 1.0 run 100000 write_data system_after_eq3_npt.data |
2楼2019-05-05 10:07:04
3楼2019-05-05 10:08:31