计算化学公社

标题: 求助:pull中pull-coord1-k计算问题(已解决) [打印本页]

作者
Author:
haijiayuan    时间: 2023-2-8 10:24
标题: 求助:pull中pull-coord1-k计算问题(已解决)
本帖最后由 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


作者
Author:
haijiayuan    时间: 2023-2-8 14:22
jerkwin的在线计算兼单位换算器 可以换算
作者
Author:
zyp1999    时间: 2023-3-1 11:28
本帖最后由 zyp1999 于 2023-3-1 15:55 编辑
haijiayuan 发表于 2023-2-8 14:22
jerkwin的在线计算兼单位换算器 可以换算

这个在线单位换算的话,假如我需要加压150mpa,盒子面积30nm^2,是换算完是270.996吗C:\Users\YT-PC592\Desktop\微信图片_20230301112623.png
作者
Author:
haijiayuan    时间: 2023-4-11 15:41
zyp1999 发表于 2023-3-1 11:28
这个在线单位换算的话,假如我需要加压150mpa,盒子面积30nm^2,是换算完是270.996吗

应该是这样
作者
Author:
姜酱将    时间: 2023-12-25 20:33
zyp1999 发表于 2023-3-1 11:28
这个在线单位换算的话,假如我需要加压150mpa,盒子面积30nm^2,是换算完是270.996吗

请问一下按照这个换算出来的力你能实现这个压缩的效果吗
作者
Author:
zyp1999    时间: 2024-1-2 16:42
姜酱将 发表于 2023-12-25 20:33
请问一下按照这个换算出来的力你能实现这个压缩的效果吗

我上面说错了,换出来应该是2709.96,然后这样换出来的力是有压缩效果的




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3