LAMMPS模拟（二）碳纳米管表面冰块脱附模拟

最新推荐文章于 2023-11-18 20:17:37 发布

转载最新推荐文章于 2023-11-18 20:17:37 发布 · 840 阅读

CC 4.0 BY-SA版权

原文链接：https://mp.weixin.qq.com/s?__biz=Mzg5Njc1MzgwMg==&mid=2247485092&idx=2&sn=6e87aabbbe5fd13dd2d4c7c157c0fc41&chksm=c07d0f41f70a8657b3adb155e6628f98e33730a6ee27ace3875e9ba578f20688eeb977127150&scene=126&&sessionid=0

文章标签：

#java #小程序 #python #编程语言 #css

本文详细介绍了使用LAMMPS进行冰晶附着力模拟的设置过程，包括参数设定、单元系统、原子类型定义、设置、运行步骤，以及关键变量如ε_WG和T的设置。特别提到了'元素魔方'这个科研工具的推荐。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

以下为in文件内容

# ----------------- Init Section-----------------


 


  variable eps_WG equal 0.21084


  variable T equal 200


 


# ----------------- Init Section-----------------


   units        real


   atom_style   full


       neigh_modifyevery 1 delay 5


       newtonon on


   pair_style   lj/cut/tip4p/long 1 21 1 0.1577 12.0 8.5


   bond_style   harmonic


   angle_style  harmonic


   kspace_style pppm/tip4p 0.0001


 


# ----------------- Atom Definition Section-----------------


 


read_data "ice_adhesion.data"


 


# ----------------- Settings Section-----------------


 


   bond_coeff   1     1000.0 0.9572


   angle_coeff  1   1000.0 104.52


   pair_coeff   1 1  0.21084 3.1668 #OO


       pair_coeff   1 2 0.0     0.0  #OH


       pair_coeff   1 3 ${eps_WG}     3.19  #OC 2.39 3.55


   pair_coeff   2 2  0.0    0.0  #HH


       pair_coeff   2 3 0.0     0.0  #HC


       pair_coeff   3 3 0.0     0.0  #CC


      


   group spce type  1  2


       groupgraphene type 3


 


# ----------------- Run Section-----------------


 


 timestep 2.0


 thermo 20


 fixfSHAKE spce shake 0.0001 10 0 b 1 a 1


 compute spceTemp spce temp


 velocity spce create $T 4928459 rot yes distgaussian temp spceTemp


 


 dumpmyDump all atom 100 dump.lammpstrj


 


 fix  fxnvt spce nvt temp $T $T 100.0


 fix_modify fxnvt temp spceTemp


 


 compute force graphene property/atom fz


 compute shear graphene reduce sum c_force


 thermo_style custom step temp c_shear


 thermo_modify temp spceTemp


 


 run  500000


 


 variable acc_g equal elapsed*5.734e-9


 fix1 spce gravity v_acc_g vector 0 0 1


 


 run  800000

了解更多代码详情，请扫描下方二维码！

PS：推荐一个巨好用的科研小程序“元素魔方”，专注【材料测试+模拟计算】科研服务，承接第一性原理、分子动力学、量子化学计算等服务，有需要的朋友可以点击文章开头小程序，了解详情哦~

如需投稿或合作

发送邮箱YlmmEditor@ylmm-tech.com

获取更多科研硬货可关注公众号

获取更多福利、经验交流可加群