LAMMPS做GCMC的in文件怎么写？

Divine亦倾

请教下，用LAMMPS里的fix gcmc 命令做MC模拟，那个In文件怎么写啊？一直出错，求告知。。。

winterzen

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
   
7. 
   
8. # variables modifiable using -var command line switch
   
9. 
   
10. variable        mu index -1.25
    
11. variable        temp index 2.0
    
12. variable        disp index 1.0
    
13. variable        lbox index 5.0
    
14. 
    
15. # global model settings
    
16. 
    
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
    
21. 
    
22. # box
    
23. 
    
24. region                box block 0 ${lbox} 0 ${lbox} 0 ${lbox}
    
25. create_box        1 box
    
26. 
    
27. # lj parameters
    
28. 
    
29. pair_coeff        * * 1.0 1.0
    
30. mass                * 1.0
    
31. 
    
32. # gcmc
    
33. 
    
34. fix             mygcmc all gcmc 1 100 100 1 29494 ${temp} ${mu} ${disp}
    
35. 
    
36. # averaging
    
37. 
    
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]
    
47. 
    
48. # output
    
49. 
    
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
    
56. 
    
57. # run
    
58. 
    
59. run             10000
    

复制代码

zjuhjx

看你做什么体系了，其实可以考虑用一些其他的软件来做MC，比如TOWHEE或者RASPA

Divine亦倾

winterzen 发表于 2017-9-13 22:12
LAMMPS源码包里examples/gcmc目录下自带了两个例子，以下是自带的LJ简单流体的in：

谢谢哦，我去多了解下

Divine亦倾

zjuhjx 发表于 2017-9-14 21:29
看你做什么体系了，其实可以考虑用一些其他的软件来做MC，比如TOWHEE或者RASPA

目前就会MS ，LAMMPS还在学，in文件还不太会写

zjuhjx

Divine亦倾 发表于 2017-9-18 10:09
目前就会MS ，LAMMPS还在学，in文件还不太会写

那我说实话，你MC的基础理论要先看懂的比较多再做比较好，我也是在学，感觉MC里面有的概念不是很好理解。哈哈，加油

Divine亦倾

zjuhjx 发表于 2017-9-20 18:00
那我说实话，你MC的基础理论要先看懂的比较多再做比较好，我也是在学，感觉MC里面有的概念不是很好理解。 ...

好的呀

luogaoyang123

winterzen 发表于 2017-9-13 22:12
LAMMPS源码包里examples/gcmc目录下自带了两个例子，以下是自带的LJ简单流体的in：

lltll

结果如何？in文件跑成功了吗？可以分享一下in文件吗？我已经error和warning了一年多了，还没从坑里爬出来

lltll

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

+1，一年多了，各种各样的报错，可以拉个群抱团吗？

Cattttttttt

lltll 发表于 2023-9-25 23:43
+1，一年多了，各种各样的报错，可以拉个群抱团吗？

+1+1

刘小胖

Cattttttttt 发表于 2023-9-26 21:35
+1+1

+1+1

NWH

请问楼上有群了吗，想抱团＋1

CHECHE

NWH 发表于 2025-2-23 13:32
请问楼上有群了吗，想抱团＋1

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