| 查看: 824 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者wsws1234将赠送您 10 个金币 | ||
[求助]
能量最小化和驰豫问题
|
||
小弟在做烃类在石英表面的吸附，初始模型构建如图所示，其中石英基底在模拟中固定，下面是我的in文件，我想先能量最小化，然后在npt系综下驰豫得到一个常温常压下合理的状态，但是烃类分子老往盒子上边的真空区域跑，是什么原因导致的呢？？各位大侠帮忙看看我in文件中关于能量最小化和驰豫部分的参数和命令是否合理 # 3d text simulation # Initialization dimension 3 boundary p p p units real atom_style full neighbor 2.0 bin neigh_modify every 1 delay 5 check yes pair_style lj/cut/coul/long 10.0 pair_modify mix arithmetic kspace_style pppm 1.0e-4 bond_style harmonic angle_style harmonic dihedral_style opls read_data data.text region 1 block INF INF INF INF 0 16.7 group jidi region 1 region 2 block INF INF INF INF 16.7 40 group move region 2 fix 1 jidi setforce 0.0 0.0 0.0 min_style sd minimize 1.0e-8 1.0e-8 10000 10000 velocity move create 298.0 4928459 dist gaussian velocity jidi set 0 0 0 fix 2 jidi nve fix 3 move npt temp 298.0 298.0 100.0 iso 1.0 1.0 100.0 timestep 1 thermo 10 thermo_style one dump 1 all atom 2000 dump.text run 1000000  微信截图_20200808161948.png |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
40