24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

菜鸟他哥

银虫 (初入文坛)

应助: 1 (幼儿园)
金币: 308.9
散金: 10
红花: 1
帖子: 40
在线: 18.7小时
虫号: 5236302
注册: 2016-11-14
专业: 半导体电子器件

[求助] 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

回复此楼