XUXINGYA 发表于 2025-5-21 22:45 您好,目前关于Na-F和Na-O最大峰高不同的问题已经解决,就是通过调节所有原子电荷,都乘以0.8。https://www.sciencedirect.com/sc ... i/S2451929424004728这篇文献有提到为什么要乘以0.8 |
spinel 发表于 2025-4-26 13:42 您好,电荷数*0.8之后您解决了吗,我也遇到和你一样的问题,但是*0.8之后体系从阴阳离子聚集变成了出现一大块空洞区域,不知道你后续如何了 |
slxc920113 发表于 2025-4-26 12:27 好的,谢谢老师 |
spinel 发表于 2025-4-25 15:45 所有的原子电荷都乘以系数,目的是削弱静电作用,本质上是因为采用的力场无法描述阴离子的极化作用,采用缩放电荷的形式进行修正。通常scale系数是1/折光率。 |
slxc920113 发表于 2025-4-25 15:33 老师,您可以简单解释一下为什么这样操作之后可以改变趋势吗?以及乘以0.8应该修改哪个位置的电荷呢? |
spinel 发表于 2025-4-25 15:25 趋势上应该是的。 |
slxc920113 发表于 2025-4-25 14:23 老师,请问您的意思是把Li+和PF6-的电荷都乘以0.8,就可以实现和opls-AA力场一样的结果了吗(即金属离子和溶剂之间的RDF更占优势)? |
BUAANirvana 发表于 2025-3-5 16:19 把Li和PF6的电荷乘以0.8,溶剂分子不变。 |
BUAANirvana 发表于 2025-3-5 16:19 同学,请问你的问题解决了吗,我后续重新模拟了文献(https://advanced.onlinelibrary.w ... 1002/aenm.202405363)中的RDF,该文献中用的也是GAFF力场,我使用和文献中一样的参数,粒子数量和盒子的构建也和文献中一样,依然是Na与F之间的g(r)占优势。 |
yxdd98 发表于 2025-3-5 12:38 现在有个新体系用sobtop产生力场文件可能更适合,OPLS-AA感觉乱糟糟的,LigParGen产生力场文件老出问题。另外就是好奇为什么GAFF会产生这样的结果,担心有什么没注意到的地方。 |
既然采用OPLS-AA力场得到和文献类似的结果了,为什么还要改用GAFF力场呢?是为了测试还是说必须要用GAFF呢? |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-13 21:21 , Processed in 0.156270 second(s), 25 queries , Gzip On.