| 查看: 1873 | 回复: 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 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有229人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有21人回复
-
欢迎加入课题组
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
看到同一作者 1 种数据，2 篇文章，间隔 4 年
已经有26人回复
- lammps里log文件里一些数据的单位？
已经有10人回复
- matlab调试问题
已经有3人回复
- 关于势函数的求助
已经有23人回复
- 【原创】MCM-41的孔径分布的计算
已经有13人回复
- 【讨论】势函数深刻的意义
已经有3人回复
- 【求助】DL_POLY的config文件如何生成??
已经有98人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
