计算化学公社

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

[Lammps] Born-Mayer-Huggins势函数相关

[复制链接 Copy URL]

16

帖子

0

威望

233

eV
积分
249

Level 3 能力者

各位大佬,本人在使用lammps模拟混合熔盐体系时对势函数系数有疑问,在有的文献中Born-Mayer-Huggins势的描述为:

此式可使用混合计算势函数系数,但相比如下:

此为lammps说明书中对BMH的描述,其存在首项以及尾项的差别,请问此分歧应如何解释?在计算时是等同于D=qi*qj还是其他算法?还望大佬们百忙之中解惑答疑,不胜感激。

16

帖子

0

威望

233

eV
积分
249

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2023-12-5 10:30:23 | 只看该作者 Only view this author
因为很多文献中给出的为第一个式,而在使用lammps中BMH式的pair_style born/coul/long时所需要参数为A,σ,ρ,C,D五个pair_coeff,无法得到D值让我很疑惑

161

帖子

0

威望

560

eV
积分
721

Level 4 (黑子)

极天的流星雨

3#
发表于 Post on 2023-12-6 00:22:33 | 只看该作者 Only view this author
D=qi*qj?D下面的分母不是r的8次方吗,这能=?
主要使用软件:VASP
目前从事稀土发光相关的第一性原理计算。
以前曾从事:蛋白质水溶液的分子动力学计算(NAMD、Gromacs);钙钛矿相关的第一性原理计算。

16

帖子

0

威望

233

eV
积分
249

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-12-6 08:24:00 | 只看该作者 Only view this author
shinkon 发表于 2023-12-6 00:22
D=qi*qj?D下面的分母不是r的8次方吗,这能=?

感觉不能等于,但我也很疑惑,不是这样的话那没给出的D怎么计算?文献里似乎也没说忽略一部分的范德华力啊

3

帖子

0

威望

191

eV
积分
194

Level 3 能力者

5#
发表于 Post on 2023-12-14 20:59:17 | 只看该作者 Only view this author
你好,你可以贴一下你提到的Born势函数公式的那篇文献吗?
你说的式(1)中的公式与Matsui开发的势函数参数是一致的,可以通过转化计算得到设置Buckingham势函数的参数。
我所看到的文献中对Born势函数的数学描述是这样的:

DOI: 10.1016/j.jnoncrysol.2014.12.030
或者是这样的:

DOI: 10.1016/j.ces.2020.116285

16

帖子

0

威望

233

eV
积分
249

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-12-19 16:24:47 | 只看该作者 Only view this author
本帖最后由 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的计算也是分两种方法,现在就是很迷。

16

帖子

0

威望

233

eV
积分
249

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-12-19 16:26:53 | 只看该作者 Only view this author
Chengen 发表于 2023-12-14 20:59
你好,你可以贴一下你提到的Born势函数公式的那篇文献吗?
你说的式(1)中的公式与Matsui开发的势函数参 ...

首项的差别我现在大概清楚可能是文献中特定体系他不考虑库伦相互作用

1

帖子

0

威望

17

eV
积分
18

Level 1 能力者

8#
发表于 Post on 2024-4-24 09:06:26 | 只看该作者 Only view this author
chemdumpling 发表于 2023-12-5 10:30
因为很多文献中给出的为第一个式,而在使用lammps中BMH式的pair_style born/coul/long时所需要参数为A,σ ...

您好,请问这个问题您解决了吗,这个势函数的pair_coeff应该怎么设置,C,D值应该怎么取,请指教

134

帖子

0

威望

1664

eV
积分
1798

Level 5 (御坂)

9#
发表于 Post on 2024-7-22 16:57:24 | 只看该作者 Only view this author
chemdumpling 发表于 2023-12-19 16:24
我看到比较新的文献是:
DOI:10.1016/j.nanoen.2017.07.020
还有较久远的:

大部分文献应该是这个形式吧,你贴的文献里面的形式并不是你给的那样。不一样的话应该就是忽略了偶极子和四极子之间的相互作用

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 01:27 , Processed in 0.301228 second(s), 29 queries , Gzip On.

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