“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 388|回复 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-4-13 11:42 , Processed in 0.229206 second(s), 23 queries , Gzip On.

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