| 查看: 4084 | 回复: 22 | |||
[交流]
求助：lammps非周期边界已有6人参与
|
|
向大家请教一下，lammps里的非周期边界f,s,m是什么意思？物理意义是什么？手册上只是形象的描述，不算解释吧。 lammps里其他的命令是怎么利用非周期边界条件的？只是用于原子丢失的判定吗？ 假如我设置一个很大的盒子，原子只占里面的很小部分，肯定不会丢失，这时f,s,m下是一样的结果吗？ |
回复此楼» 猜你喜欢
审稿
已经有6人回复
-
江南大学到瑞士招聘，称取消非升即走，改预聘+长聘
已经有21人回复
-
基金得中
已经有4人回复
-
F口信息学部拿面上，大概需要什么样的成果
已经有16人回复
-
34岁读博士晚吗
已经有38人回复
-
国家公派访问学者申请结果出了吗？
已经有4人回复
-
关于读博感觉自己很抓马
已经有9人回复
-
面青地会评时间？？？
已经有4人回复
-
2024国社科通讯评审
已经有39人回复
-
Nature 11日发文，中国著名学者们称造假迫不得已
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一下衰减的周期函数是不是周期函数
已经有6人回复
- lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
已经有8人回复
- fluent里面将进出口设置为周期性边界条件
已经有17人回复
- gambit中周期性边界网格的画法
已经有5人回复
- fluent中如何设置周期边界？
已经有4人回复
- read_restart非周期性边界条件
已经有4人回复
- fluent周期边界计算的困惑
已经有9人回复
- 用lammps模拟填充水的碳管的热导率，能设置碳管为固定边界，水为周期性边界条件吗？
已经有6人回复
- 3维周期性边界条件设置
已经有9人回复
- 石墨中加磁场，周期性边界条件怎么定义
已经有5人回复
- 分子动力学模拟时如何利用周期性边界条件？？
已经有5人回复
- fix deform 指令拉伸求教
已经有10人回复
- MS分子模拟的截断半径该怎么设？
已经有4人回复
- 周期性边界条件相关问题
已经有6人回复
- 如何设置周期边界条件
已经有9人回复
- 周期性边界造成的RMSD异常值
已经有12人回复
- s s p边界 原子数增多的问题
已经有5人回复
- 【求助】vasp为什么采用周期性边界条件（或超原胞模型）来处理原子分子表面体系等，有
已经有8人回复
- 【讨论】分子动力学初始速度的选取
已经有16人回复
- 【讨论】周期性边界条件什么时候用
已经有11人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
2楼2013-12-26 10:04:39
3楼2013-12-27 12:48:56
4楼2013-12-27 15:46:58
lhllookatme
木虫 (初入文坛)
- 模拟EPI: 1
- 应助: 4 (幼儿园)
- 金币: 3109.7
- 红花: 1
- 帖子: 25
- 在线: 19.3小时
- 虫号: 1232860
- 注册: 2011-03-14
- 性别: GG
- 专业: 动力学与控制
5楼2013-12-28 14:28:34
6楼2013-12-30 09:44:01
|
in文件： dimension 3 units metal boundary p p f atom_style full neighbor 3.0 bin neigh_modify delay 5 every 1 check yes read_data gufinal2.data # force field pair_style hybrid airebo 2.5 0 0 lj/cut 10.0 pair_coeff * * airebo CH.airebo C C C C pair_coeff 1 2*4 lj/cut 0.00284 3.40 pair_coeff 2 3*4 lj/cut 0.00284 3.40 pair_coeff 3 4 lj/cut 0.00284 3.40 # define region before define group region rfix block INF INF INF INF INF -8 units box region r1 block INF INF INF INF -8 -4.5 units box region r2 block INF INF INF INF -4.5 -1 units box region r3 block INF INF INF INF -1 INF units box region th1 block INF INF INF 2 -8 INF units box region th2 block INF INF 31 INF -8 INF units box group gfix region rfix group g1 region r1 group g2 region r2 group g3 region r3 group gth1 region th1 group gth2 region th2 group gth union gth1 gth2 group gthg12 subtract gth g3 group gthg3 subtract gth gthg12 group g12 union g1 g2 group gf12 union gfix g12 # compute and velocity compute temp6 gthg3 temp/com variable vg3x equal vcm(g3,x)#x component of velocity of the center of mass of group g3 variable pg3x equal xcm(g3,x)#x component of position of the center of mass of group g3 variable pg3y equal xcm(g3,y)#y component of position of the center of mass of group g3 variable pg3z equal xcm(g3,z)#z component of position of the center of mass of group g3 velocity gfix set 0.0 0.0 0.0 units box velocity g12 create 10 4928459 rot yes dist gaussian units box velocity g3 create 10 4928459 rot yes dist gaussian units box # fixes fix 1 all nve fix 2 gthg12 langevin 10 10 0.3 699483 zero yes fix 3 gthg3 langevin 10 10 0.1 699483 zero yes fix_modify 3 temp temp6 fix 4 gfix setforce 0.0 0.0 0.0 #relax min_modify dmax 0.01 minimize 1.0e-6 1.0e-8 1000 100000 neigh_modify exclude group gfix gfix timestep 0.001 thermo 100 thermo_style custom step v_pg3x v_pg3y v_pg3z v_vg3x dump 1 all xyz 200 gufinal2.xyz dump 2 all atom 200 gufinal2.lammpstrj run 50000 velocity g3 set 1.00 NULL NULL units box sum yes run 100000 盒子尺寸： 0.000000000 85.2168 xlo xhi 0.000000000 34.4332 ylo yhi -8.500000000 18.257102924 zlo zhi 模拟的是碳纳米管在石墨烯上运动，z方向原子不会丢失。 只改 boundary p p f 这句， v_pg3x v_pg3y v_pg3z v_vg3x 结果不一样。 修改data文件里的 -8.500000000 18.257102924 zlo zhi 使盒子尺寸稍微变化一点，结果也不一样。 |
7楼2013-12-30 09:59:17
8楼2013-12-31 12:00:34
9楼2014-01-02 17:04:16
10楼2014-01-07 21:00:32
20