求助lammps固液相双相模型NVT弛豫问题

3046083274p

Lammps中，我构建了一个固液相的双相模型，如图所示，左边是Ca(OH)2溶液（cvff力场描述），右边是凝胶模型（clayFF力场描述），中间是两层石墨烯层（目的是阻止两侧原子穿越）。目前我正在进行该模型的NVT弛豫。发现弛豫过程中，能量最小化进行几十步就结束了，并且溶液模型的原子会从左侧进入固相模型中，并且填充了固相中的真空层（使得我之后的液相扩散无法进行）。请问为什么会出现这样的问题，怎样解决呢？谢谢！下附lammps的NVT弛豫in文件

1. echo screen
   
2. units          real                # 单位制：real（kcal/mol, Å, fs）
   
3. atom_style     full              # 原子类型定义[1,3](@ref)
   
4. boundary       p p p               # 周期性边界条件[3,10](@ref)
   
5. timestep       0.5                 # 时间步长0.5 fs
   
6. angle_style harmonic
   
7. bond_style harmonic
   
8. neighbor 2.0 bin
   
9. neigh_modify every 1 delay 0 check yes
   
10. 
    
11. # ------------------- 读取文件 -------------------
    
12. pair_style     lj/cut/coul/long 12 12  # LJ截断12Å，库仑截断12Å[6,8](@ref)
    
13. kspace_style pppm 1e-4
    
14. read_data combined_x_wall2.data
    
15. 
    
16. group graphene type 13
    
17. fix freeze graphene setforce 0 0 0
    
18. 
    
19. # ------------------- 阶段1：能量最小化-------------------
    
20. min_style cg
    
21. minimize 1e-6 1e-8 10000 100000
    
22. 
    
23. # ------------------- 阶段2：NVT -------------------
    
24. fix            1 all nvt temp 300 300 100 # Nose-Hoover恒温器，Tdamp=0.1 ps[4,5](@ref)
    
25. thermo         1000                # 每1000步输出热力学信息
    
26. thermo_style   custom step temp pe etotal press vol
    
27. dump 3 all xyz 1000 dump3.xyz
    
28. run            200000              # 100 ps（200,000 steps * 0.5 fs/step）
    
29. write_data nvt_Layer.data  # 输出弛豫后状态

复制代码

lmch

如下槽点：
echo screen          无用
timestep       0.5      步长过小
neigh_modify every 1 delay 0 check yes   把默认值又写了一遍
另外最重要的，你石墨烯设置pair_coeff了吗？

aw6we6

如果不希望左侧溶液通过左侧边界进入到固相中，可以把左右侧周期性边界变成f，并增加上fix zwalls all wall/reflect zhi EDGE zlo EDGE（x，y，z根据自己的模型设置）