计算化学公社

标题: Born-Mayer-Huggins势函数相关 [打印本页]

作者
Author:
chemdumpling    时间: 2023-12-4 19:50
标题: Born-Mayer-Huggins势函数相关
各位大佬,本人在使用lammps模拟混合熔盐体系时对势函数系数有疑问,在有的文献中Born-Mayer-Huggins势的描述为:
(, 下载次数 Times of downloads: 32)
此式可使用混合计算势函数系数,但相比如下:
(, 下载次数 Times of downloads: 28)
此为lammps说明书中对BMH的描述,其存在首项以及尾项的差别,请问此分歧应如何解释?在计算时是等同于D=qi*qj还是其他算法?还望大佬们百忙之中解惑答疑,不胜感激。

作者
Author:
chemdumpling    时间: 2023-12-5 10:30
因为很多文献中给出的为第一个式,而在使用lammps中BMH式的pair_style born/coul/long时所需要参数为A,σ,ρ,C,D五个pair_coeff,无法得到D值让我很疑惑
作者
Author:
shinkon    时间: 2023-12-6 00:22
D=qi*qj?D下面的分母不是r的8次方吗,这能=?
作者
Author:
chemdumpling    时间: 2023-12-6 08:24
shinkon 发表于 2023-12-6 00:22
D=qi*qj?D下面的分母不是r的8次方吗,这能=?

感觉不能等于,但我也很疑惑,不是这样的话那没给出的D怎么计算?文献里似乎也没说忽略一部分的范德华力啊
作者
Author:
Chengen    时间: 2023-12-14 20:59
你好,你可以贴一下你提到的Born势函数公式的那篇文献吗?
你说的式(1)中的公式与Matsui开发的势函数参数是一致的,可以通过转化计算得到设置Buckingham势函数的参数。
我所看到的文献中对Born势函数的数学描述是这样的:
(, 下载次数 Times of downloads: 28)
DOI: 10.1016/j.jnoncrysol.2014.12.030
或者是这样的:
(, 下载次数 Times of downloads: 25)
DOI: 10.1016/j.ces.2020.116285
作者
Author:
chemdumpling    时间: 2023-12-19 16:24
本帖最后由 chemdumpling 于 2023-12-19 16:38 编辑
Chengen 发表于 2023-12-14 20:59
你好,你可以贴一下你提到的Born势函数公式的那篇文献吗?
你说的式(1)中的公式与Matsui开发的势函数参 ...

我看到比较新的文献是:
DOI:10.1016/j.nanoen.2017.07.020
还有较久远的:
DOI:10.1080/00268977300102671
他这边对后两项的描述是范德华力中偶极子-偶极子和偶极子-四极子色散能,但是具体这两式中差的项是忽略了还是在特定体系下不予考虑这边完全没有给出解释,此外在90年代的文章中对Aij的计算也是分两种方法,现在就是很迷。

作者
Author:
chemdumpling    时间: 2023-12-19 16:26
Chengen 发表于 2023-12-14 20:59
你好,你可以贴一下你提到的Born势函数公式的那篇文献吗?
你说的式(1)中的公式与Matsui开发的势函数参 ...

首项的差别我现在大概清楚可能是文献中特定体系他不考虑库伦相互作用
作者
Author:
Gassifflen    时间: 2024-4-24 09:06
chemdumpling 发表于 2023-12-5 10:30
因为很多文献中给出的为第一个式,而在使用lammps中BMH式的pair_style born/coul/long时所需要参数为A,σ ...

您好,请问这个问题您解决了吗,这个势函数的pair_coeff应该怎么设置,C,D值应该怎么取,请指教
作者
Author:
thor    时间: 2024-7-22 16:57
chemdumpling 发表于 2023-12-19 16:24
我看到比较新的文献是:
DOI:10.1016/j.nanoen.2017.07.020
还有较久远的:

大部分文献应该是这个形式吧,你贴的文献里面的形式并不是你给的那样。不一样的话应该就是忽略了偶极子和四极子之间的相互作用
(, 下载次数 Times of downloads: 24)




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3