| 查看: 1899 | 回复: 0 | ||
[求助]
关于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 |
回复此楼» 猜你喜欢
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有147人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有31人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有2人回复
-
PVA溶解
已经有1人回复
-
基金申请
已经有40人回复
找到一些相关的精华帖子，希望有用哦~
看到同一作者 1 种数据，2 篇文章，间隔 4 年
已经有26人回复
- lammps里log文件里一些数据的单位？
已经有10人回复
- matlab调试问题
已经有3人回复
- 关于势函数的求助
已经有23人回复
- 【原创】MCM-41的孔径分布的计算
已经有13人回复
- 【讨论】势函数深刻的意义
已经有3人回复
- 【求助】DL_POLY的config文件如何生成??
已经有98人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人