计算化学公社

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

[Lammps] 如何编写martini3在lammps中的in文件?

[复制链接 Copy URL]

7

帖子

0

威望

33

eV
积分
40

Level 2 能力者

想请问各位老师,如何编写martini3在lammps中的in文件?或者说,如何转化可以使martini3力场可以在lammps中运行?有没有比较好的教程或者方法的,谢谢各位大佬

48

帖子

0

威望

451

eV
积分
499

Level 3 能力者

2#
发表于 Post on 2025-3-14 08:35:33 | 只看该作者 Only view this author
直接写就好,Martini官网有提供martini3的itp格式文件,里面有详细的每种珠子的12-6势的参数,不过用的单位是nm和KJ/mol,需要转化为埃米和Kcal/mol,

至少我没找到可以直接使用Martini3力场生成lammps拓扑文件的软件,不过可以试着用 PolyPly生成gmx格式的文件,然后用GRO2LAM转化为lammps的data文件

不过这个流程我没实际做过,我是写python脚本转化的

GRO2LAM: https://github.com/hernanchavezthielemann/GRO2LAM

PolyPly:https://github.com/marrink-lab/polyply_1.0

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
含光君 + 5 谢谢分享

查看全部评分 View all ratings

7

帖子

0

威望

33

eV
积分
40

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-3-19 00:44:50 | 只看该作者 Only view this author
jrfjrf123 发表于 2025-3-14 08:35
直接写就好,Martini官网有提供martini3的itp格式文件,里面有详细的每种珠子的12-6势的参数,不过用的单位 ...

好的,谢谢指教。另外我想问一下,像gro2lam和intermol这种软件,是只能生成data文件吗?因为我看网上有的说能生成in文件,有的说不能,有点迷糊

48

帖子

0

威望

451

eV
积分
499

Level 3 能力者

4#
发表于 Post on 2025-3-19 22:33:29 | 只看该作者 Only view this author
本帖最后由 jrfjrf123 于 2025-3-19 22:41 编辑
230820071 发表于 2025-3-19 00:44
好的,谢谢指教。另外我想问一下,像gro2lam和intermol这种软件,是只能生成data文件吗?因为我看网上有 ...

gro2lam可以看这篇文章:


From GROMACS to LAMMPS-GRO2LAM.pdf (3.29 MB, 下载次数 Times of downloads: 6)


gromacs里有些md运行参数和lammps是不同的,不能直接转化,

比如gromacs里的v-rescale热浴并不等同于LAMMPS中的fix temp/rescale,详细的讨论可以看这个帖子:http://bbs.keinsci.com/thread-51538-1-1.html


还有压浴的C-rescale也是没有对应的方法的,但fix npt 是结合了Parrinello and Rahman提出Parrinello-Rahman 压浴方法(这可以看lammps手册中fix npt这一部分),可以直接转化


所以我觉得还是转化结构文件,in文件自己写好一点

7

帖子

0

威望

33

eV
积分
40

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-3-20 13:33:42 | 只看该作者 Only view this author
jrfjrf123 发表于 2025-3-19 22:33
gro2lam可以看这篇文章:

好的,谢谢指教

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 21:35 , Processed in 0.158218 second(s), 24 queries , Gzip On.

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