计算化学公社

标题: 求助lammps固液相双相模型NVT弛豫问题 [打印本页]
作者
Author: 3046083274p    时间: 2025-5-9 12:41
标题: 求助lammps固液相双相模型NVT弛豫问题
本帖最后由 3046083274p 于 2025-5-9 20:25 编辑

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

  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

  16. group graphene type 13
  17. fix freeze graphene setforce 0 0 0

  19. # ------------------- 阶段1：能量最小化-------------------
  20. min_style cg
  21. minimize 1e-6 1e-8 10000 100000

  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  # 输出弛豫后状态
复制代码

作者
Author: lmch    时间: 2025-5-10 01:17
如下槽点：
echo screen          无用
timestep       0.5      步长过小
neigh_modify every 1 delay 0 check yes   把默认值又写了一遍
另外最重要的，你石墨烯设置pair_coeff了吗？
作者
Author: aw6we6    时间: 2025-7-1 09:50
如果不希望左侧溶液通过左侧边界进入到固相中，可以把左右侧周期性边界变成f，并增加上fix zwalls all wall/reflect zhi EDGE zlo EDGE（x，y，z根据自己的模型设置）



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