计算化学公社

 找回密码 Forget password
 注册 Register
Views: 255|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] 求助:怎么将packmol的pdb文件转化成正确的lammpsdata文件

[复制链接 Copy URL]

8

帖子

0

威望

35

eV
积分
43

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 WUIDHAU 于 2025-11-21 22:46 编辑

我是一名初学者,我使用packmol构建了一份过氧化氢水溶液,其中的h2o.pdb h2o2.pdb是在ms里画的导出的,packmol的inp文件内容如下
tolerance 2.0filetype pdb
output model_box.pdb
structure H2O.pdb
    number 200
    inside box 3.0   3.0   3.0  22.0 22.0 27.0
end structure
structure H2O2.pdb
    number 10
    inside box 3.0   3.0   3.0  22.0 22.0 27.0   
end structure
之后我使用vmd的topo命令生成了一份data,文件的表头如下
LAMMPS data file. CGCMM style. atom_style full generated by VMD/TopoTools v1.7 on Fri Nov 21 19:11:22 +0800 2025
640 atoms
430 bonds
0 angles
0 dihedrals
0 impropers
2 atom types
1 bond types
0 angle types
0 dihedral types
0 improper types
-0.500000 25  xlo xhi
-0.500000 25  ylo yhi
-0.500000 30  zlo zhi
# Pair Coeffs
#
# 1  H
# 2  O
# Bond Coeffs
#
# 1  
Masses
1 1.007940 # H
2 15.999400 # O
Atoms # full

这样的话我感觉很奇怪,1是键角信息全无,2是过氧化氢和水的氧氢原子被识别为同一种,我感觉会出问题,求助以下该怎么解决这个问题
另外我的我是使用此data进行模拟时,总是qeq无法平衡,即使放大范围也不行,以下是我的势函数和能量最小化代码及势函数文件,文件取自一篇文献,文献做的是过氧化氢水溶液抛光cu的
pair_style reaxff NULL
pair_coeff * * cu.txt  C H O   Cu   
fix 101 all qeq/reax 1 0.0 10.0 1.0e-6 reaxff
min_style           cg
minimize            1e-10 1e-12 10000 10000
能请大佬看看问题吗

cu.txt

18.52 KB, 下载次数 Times of downloads: 2

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 03:34 , Processed in 0.198305 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list