|
|
本帖最后由 haijiayuan 于 2023-2-8 14:04 编辑
如图所示,想使ab之间产生特定方向(横向)的压力进行md模拟,采用的策略是对ab进行恒力牵引,请问如何将(kJ mol-1 nm-1)换算成ab之间的压力?
;pul com NVT
pull = yes
pull-pbc-ref-prev-step-com = yes
pull-ngroups = 2
pull-ncoords = 1
pull-group1-name = a
pull-group1-pbcatom = 111
pull-group2-name = b
pull-group2-pbcatom = 222
pull-coord1-type = constant-force
pull-coord1-geometry = direction-periodic
pull-coord1-groups = 1 2
pull-coord1-vec = 0.0 0.0 -1
pull-coord1-k = 0.9477367 ;kJ mol-1 nm-1
|
-
22.png
(690.62 KB, 下载次数 Times of downloads: 8)
|