24小时热门版块排行榜    

查看: 1856  |  回复: 0

youminguo

金虫 (小有名气)

[求助] 关于LAMMPS 如果定义charge

我想模拟碳原子在晶格铁中扩散,现在学要外加电场 ,把碳原子看成离子,请问如何定义input文件中碳原子,
#------------LAMMPS carbon diffusion in FCC-iron input-----------



#------------initialization-----------
dimension                  3
units                        metal
boundary                 p p p
atom_style                 atomic


#----------simulation box
variable                  x index 1100
variable                  L equal 5
lattice                  bcc 2.8553
region                   box block 0 $L 0 $L 0 $L
create_box                2 box
create_atoms              1 box
create_atoms             2 single 1.4276 1.4276 2.885

# ----------------------------- Settings atoms

group                 Fe type 1
group              C type 2
mass                                1 55.84500
mass                                2 12.01070
atom_modify         C charge
set                type 2 charge 0.5

neighbor                        0.3 bin
neigh_modify                every 5 delay 0 check yes


# ----------------------------- Potential
pair_style                        eam/alloy
pair_coeff                        * *  Fe-Cbecquart2011.eam Fe C

# -----------------------------
velocity        all create $x 825577 dist gaussian
timestep      0.001

compute 2 C msd
fix             2 all npt temp $x $x 0.2 iso 0 0 10
fix EFIELD all efield 1.5 0.0 0.0
run 10000
unfix EFIELD
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 youminguo 的主题更新
信息提示
请填处理意见