计算化学公社

 找回密码 Forget password
 注册 Register

求助:Gromacs模拟钴酸锂固液界面问题

查看数: 2664 | 评论数: 16 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
xjw
发布时间: 2023-11-11 14:04

正文摘要:

老师好,我想做一个关于钴酸锂正极和电解液的固液界面的模拟,但是好像没看到钴酸锂的力场(结构如下图一),我在谷歌查到一篇关于钴酸锂的文章,他们用的UFF力场,但是好像只有原子对的力场参数(图二),没有看到 ...

回复 Reply

chuxuedexiaobai 发表于 Post on 2024-6-28 10:53:51
Jill-W 发表于 2024-6-26 22:29
请问你解决了嘛?我目前也是固液界面,然后不固定固体的话,固体石墨烯界面就会出现滑移,请问怎么解决的 ...

目前找到的问题是我使用的gromacs2020版本有bug导致我冻结时会报错,要时平面固定就只能冻结然后只能跑nvt,所以的换一个gromacs版本
Jill-W 发表于 Post on 2024-6-26 22:29:29
chuxuedexiaobai 发表于 2024-4-18 14:18
我看冻结不能和压力耦合一起用,因为后面还要跑npt,所以没有用冻结

请问你解决了嘛?我目前也是固液界面,然后不固定固体的话,固体石墨烯界面就会出现滑移,请问怎么解决的呢?
chuxuedexiaobai 发表于 Post on 2024-4-18 14:18:30
slxc920113 发表于 2024-4-18 11:47
直接冻结组不就得了。

我看冻结不能和压力耦合一起用,因为后面还要跑npt,所以没有用冻结
slxc920113 发表于 Post on 2024-4-18 11:47:19
chuxuedexiaobai 发表于 2024-4-18 10:29
想问一下,gmx里面对于金属基面的限制文件应该怎么写才能固定金属面中所有的原子呢,我铁面19200个铁原子 ...

直接冻结组不就得了。
slxc920113 发表于 Post on 2024-4-18 11:46:56
gjs 发表于 2024-4-17 17:38
请问使用AutoFF产生一个Li板的力场文件应该怎么做呢,我使用金属/共轭有机框架选择的OPLS力场,但是显示 ...

瞎用力场,Fe金属你得用FCC-Metal
chuxuedexiaobai 发表于 Post on 2024-4-18 10:29:49
gjs 发表于 2024-4-17 17:38
请问使用AutoFF产生一个Li板的力场文件应该怎么做呢,我使用金属/共轭有机框架选择的OPLS力场,但是显示 ...

想问一下,gmx里面对于金属基面的限制文件应该怎么写才能固定金属面中所有的原子呢,我铁面19200个铁原子,我在itp,posre.itp均只写了一个,然后位置限制后跑完nvt铁面的原子就散乱了,对于金属面的固定要怎么设置才是正确的呀,可以解答一下嘛,非常感谢
gjs 发表于 Post on 2024-4-17 17:38:56
本帖最后由 gjs 于 2024-4-17 17:41 编辑
slxc920113 发表于 2024-1-9 17:06
可以将就用。虽然混合规则不同,但是定性看一些现象没什么问题。

请问使用AutoFF产生一个Li板的力场文件应该怎么做呢,我使用金属/共轭有机框架选择的OPLS力场,但是显示的Li板中的Li都是带电荷的,按道理来说应该不带电荷的呀。

锂板-OPLS.png (385.55 KB, 下载次数 Times of downloads: 41)

锂板-OPLS.png
xjw 发表于 Post on 2024-1-9 19:50:32
slxc920113 发表于 2024-1-9 17:06
可以将就用。虽然混合规则不同,但是定性看一些现象没什么问题。

好的,感谢
slxc920113 发表于 Post on 2024-1-9 17:06:37
xjw 发表于 2024-1-9 16:14
您好,请问下我用Autoff网站生成的LCO的UFF力场,但是我电解液想用opls力场,这中方案可行吗?

可以将就用。虽然混合规则不同,但是定性看一些现象没什么问题。
xjw 发表于 Post on 2024-1-9 16:14:27
slxc920113 发表于 2023-11-14 11:27
你需要自己先切好晶面,设置好真空层。
如果你想要使用COMPASS力场,可以在AuToFF的晶体模块那儿上传cif结 ...

您好,请问下我用Autoff网站生成的LCO的UFF力场,但是我电解液想用opls力场,这中方案可行吗?
xjw 发表于 Post on 2023-11-15 21:30:42
喵星大佬 发表于 2023-11-14 06:40
对于这种体系,键伸缩+键角已经足够了,实际上你需要的只有Co-O的键伸缩参数,Co-O-Co和O-Co-O的键角弯曲参 ...

收到,感谢
xjw 发表于 Post on 2023-11-15 20:17:21
slxc920113 发表于 2023-11-14 11:27
你需要自己先切好晶面,设置好真空层。
如果你想要使用COMPASS力场,可以在AuToFF的晶体模块那儿上传cif结 ...

收到,感谢,我试试
slxc920113 发表于 Post on 2023-11-14 11:27:34
你需要自己先切好晶面,设置好真空层。
如果你想要使用COMPASS力场,可以在AuToFF的晶体模块那儿上传cif结构,选择COMPASS力场,依次点击下一步,最后下载Lammps的拓扑文件。

如果你想要使用UFF力场,可以在MOF模块那儿上传cif,选择UFF力场,不需要额外选择电荷模型(QEq电荷对于这类体系误差很大,我之前做锂电体系的时候对Li-Co-Ni-Mn-O调了一套BCI电荷规则,默认在UFF力场中)也是根据提示依次点击下一步,最后下载Gromacs或者Lammps的拓扑文件,软件会自动给二面角加上谐振的约束,建议手动删除减小计算量。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
xjw + 3 好物!

查看全部评分 View all ratings

喵星大佬 发表于 Post on 2023-11-14 06:40:13
对于这种体系,键伸缩+键角已经足够了,实际上你需要的只有Co-O的键伸缩参数,Co-O-Co和O-Co-O的键角弯曲参数(UFF的八面体键角项形式转换到gmx支持的GROMOS的cos形式键角可以参考OBGMX的原文,除了Co-O-Co之外在UFF里本来就都有定义,那个在UFF2MOF定义了,相互之间的vdW项用的Buckingham形式势函数你上面也写了,自己换算一下然后用x2top就可以把拓扑构建好了,反正你也只会去做高温相钴酸锂,这种结构里二面角没啥意义
不放心可以用AutoFF看看

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

GMT+8, 2025-8-16 07:33 , Processed in 0.395282 second(s), 34 queries , Gzip On.

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