| 查看: 2273 | 回复: 6 | ||
[求助]
LAMMPS模拟黑磷，板材弯曲，寻求原因 已有3人参与
|
|
我在黑磷弛豫的过程中，黑磷板竟然弯曲了，这是一个不合理的现象，怎么解决呢？ 现贴出脚本： # temperature:300K #---------参数设置----------------------------- units metal #原子类型 dimension 3 #结构类型 boundary p p p #边界性 atom_style atomic #原子风格 read_data 123.data #读取文档 #----------------------------------------------- #--------建立盒子------------------------------- region type_2 block INF INF INF INF -100 100 units box #创建模拟盒子 group 2_atom region type_2 #定义原子组 set group 2_atom type 2 #----------------------------------------------- #-------力场参数-------------------------------- pair_style hybrid sw zbl 3.0 4.0 #原子间作用势 pair_coeff 1 1 zbl 15.0 15.0 # pair_coeff 2 2 zbl 15.0 15.0 pair_coeff 3 3 zbl 1.0 1.0 pair_coeff 1 3 zbl 1.0 1.0 pair_coeff 2 3 zbl 1.0 1.0 pair_coeff * * sw black_p.sw P1 P2 NULL neighbor 1.0 bin #截断半径 neigh_modify every 1 delay 5 check yes #--------------------------------------------- #---------输出辐照前原子位置------------------- timestep 0.1e-3 #步长 thermo 100000 #打印周期 velocity all create 300 5287287 mom yes rot yes dist gaussian #创建模拟材料速度 dump before_relax all xyz 100000 before_relax.xyz #打印材料位置 #------------------------------------------ #---------输出辐照后原子位置--------------- fix 1 all npt temp 300.0 300.0 0.05 iso 0 0 0.5 dump after_frelax all xyz 100000 after_relax.xyz run 2000000 #------------------------------------------ #---------开始退火------------------------ fix cool all nve dump cool all xyz 100000 cool.xyz run 2000000 #----------------------------------------- #---------输出文件----------------------- write_restart relaxation #------------------------------------------  9[W`~4$OD1PSF}IKFY3R_{8.png@oxox6085 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
chinaboy798
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 299.3
- 散金: 20
- 红花: 8
- 帖子: 177
- 在线: 69.6小时
- 虫号: 2007979
- 注册: 2012-09-17
- 性别: GG
- 专业: 金属材料的力学行为
2楼2017-07-01 16:57:16
3楼2017-07-02 22:16:50
4楼2017-07-02 23:18:56
5楼2017-10-02 02:01:55
6楼2017-10-09 21:07:58
z841888077
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 730.9
- 散金: 21
- 红花: 1
- 帖子: 93
- 在线: 76.5小时
- 虫号: 2330591
- 注册: 2013-03-08
- 专业: 无机纳米化学
7楼2017-10-22 09:11:54
