计算化学公社

标题: LAMMPS做GCMC的in文件怎么写？ [打印本页]
作者
Author: Divine亦倾    时间: 2017-9-13 19:36
标题: LAMMPS做GCMC的in文件怎么写？
请教下，用LAMMPS里的fix gcmc 命令做MC模拟，那个In文件怎么写啊？一直出错，求告知。。。
作者
Author: winterzen    时间: 2017-9-13 22:12
LAMMPS源码包里examples/gcmc目录下自带了两个例子，以下是自带的LJ简单流体的in：

  1. # GCMC for LJ simple fluid, no dynamics
  2. # T = 2.0
  3. # rho ~ 0.5
  4. # p ~ 1.5
  5. # mu_ex ~ 0.0
  6. # comparable to Frenkel and Smit GCMC Case Study, Figure 5.8

  8. # variables modifiable using -var command line switch

  10. variable        mu index -1.25
  11. variable        temp index 2.0
  12. variable        disp index 1.0
  13. variable        lbox index 5.0

  15. # global model settings

  17. units           lj
  18. atom_style      atomic
  19. pair_style      lj/cut 3.0
  20. pair_modify        tail no # turn of to avoid triggering full_energy

  22. # box

  24. region                box block 0 ${lbox} 0 ${lbox} 0 ${lbox}
  25. create_box        1 box

  27. # lj parameters

  29. pair_coeff        * * 1.0 1.0
  30. mass                * 1.0

  32. # gcmc

  34. fix             mygcmc all gcmc 1 100 100 1 29494 ${temp} ${mu} ${disp}

  36. # averaging

  38. variable        rho equal density
  39. variable        p equal press
  40. variable        nugget equal 1.0e-8
  41. variable        lambda equal 1.0
  42. variable             muex equal ${mu}-${temp}*ln(density*${lambda}+${nugget})
  43. fix                 ave all ave/time 10 100 1000 v_rho v_p v_muex ave one file rho_vs_p.dat
  44. variable        rhoav equal f_ave[1]
  45. variable        pav equal f_ave[2]
  46. variable        muexav equal f_ave[3]

  48. # output

  50. variable        tacc equal f_mygcmc[2]/(f_mygcmc[1]+${nugget})
  51. variable        iacc equal f_mygcmc[4]/(f_mygcmc[3]+${nugget})
  52. variable        dacc equal f_mygcmc[6]/(f_mygcmc[5]+${nugget})
  53. compute_modify  thermo_temp dynamic yes
  54. thermo_style    custom step temp press pe ke density atoms v_iacc v_dacc v_tacc v_rhoav v_pav v_muexav
  55. thermo          1000

  57. # run

  59. run             10000
复制代码
作者
Author: zjuhjx    时间: 2017-9-14 21:29
看你做什么体系了，其实可以考虑用一些其他的软件来做MC，比如TOWHEE或者RASPA
作者
Author: Divine亦倾    时间: 2017-9-18 10:08
winterzen 发表于 2017-9-13 22:12
LAMMPS源码包里examples/gcmc目录下自带了两个例子，以下是自带的LJ简单流体的in：

谢谢哦，我去多了解下
作者
Author: Divine亦倾    时间: 2017-9-18 10:09
zjuhjx 发表于 2017-9-14 21:29
看你做什么体系了，其实可以考虑用一些其他的软件来做MC，比如TOWHEE或者RASPA

目前就会MS ，LAMMPS还在学，in文件还不太会写
作者
Author: zjuhjx    时间: 2017-9-20 18:00
Divine亦倾 发表于 2017-9-18 10:09
目前就会MS ，LAMMPS还在学，in文件还不太会写

那我说实话，你MC的基础理论要先看懂的比较多再做比较好，我也是在学，感觉MC里面有的概念不是很好理解。哈哈，加油
作者
Author: Divine亦倾    时间: 2017-9-24 15:40
zjuhjx 发表于 2017-9-20 18:00
那我说实话，你MC的基础理论要先看懂的比较多再做比较好，我也是在学，感觉MC里面有的概念不是很好理解。 ...

好的呀
作者
Author: luogaoyang123    时间: 2023-3-12 18:42
本帖最后由 luogaoyang123 于 2023-8-10 20:10 编辑
winterzen 发表于 2017-9-13 22:12
LAMMPS源码包里examples/gcmc目录下自带了两个例子，以下是自带的LJ简单流体的in：


作者
Author: lltll    时间: 2023-9-22 16:45
本帖最后由 lltll 于 2023-9-22 16:46 编辑

结果如何？in文件跑成功了吗？可以分享一下in文件吗？我已经error和warning了一年多了，还没从坑里爬出来
作者
Author: lltll    时间: 2023-9-25 23:43
lltll 发表于 2023-9-22 16:45
结果如何？in文件跑成功了吗？可以分享一下in文件吗？我已经error和warning了一年多了，还没从坑里爬出来{: ...

+1，一年多了，各种各样的报错，可以拉个群抱团吗？
作者
Author: Cattttttttt    时间: 2023-9-26 21:35
lltll 发表于 2023-9-25 23:43
+1，一年多了，各种各样的报错，可以拉个群抱团吗？

+1+1
作者
Author: 刘小胖    时间: 2023-10-24 09:43
Cattttttttt 发表于 2023-9-26 21:35
+1+1

+1+1
作者
Author: NWH    时间: 2025-2-23 13:32
请问楼上有群了吗，想抱团＋1
作者
Author: CHECHE    时间: 2025-7-26 16:36
NWH 发表于 2025-2-23 13:32
请问楼上有群了吗，想抱团＋1

大佬有群吗 能一起交流下不



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3