计算化学公社

标题: 求助 Ag电极体系力场选择 [打印本页]

作者
Author: elpa 时间: 2026-1-20 21:17
标题: 求助 Ag电极体系力场选择
本帖最后由 elpa 于 2026-1-21 16:05 编辑

请教下Ag金属电极，和电解液体系进行模拟，对金属Ag 的力场参数(主要是LJ参数)该怎么选。目前查到的Hendrik Heinz的这个Interface力场是否适合模拟表面体系？https://pubs.acs.org/doi/10.1021/jp801931d

放一个类似结构图

体系主要是由两边的Ag电极，水，CO2，以及Cl离子等组成。

UFF力场的Ag参数适合吗？ Ag- Ag ϵ = 18.11591 K, σ = 2.80455 Å

作者
Author: elpa 时间: 2026-1-21 11:47
本帖最后由 elpa 于 2026-1-23 22:22 编辑

发现一个文章，好像可以做，不知是否合理。https://pubs.acs.org/doi/10.1021/la402839q

另一个新模型2-Point Polarizable model of Silver and parameters - Mendeley Data，参数是
3.5 31 Agc 5540282.0000 8355.00000
3.5 31 Age 237441.0000 358.50000

作者
Author: sobereva 时间: 2026-1-23 03:05
Heinz的参数可以用。UFF无法描述金属单质

作者
Author: elpa 时间: 2026-1-23 21:36
本帖最后由 elpa 于 2026-1-25 10:38 编辑

sobereva 发表于 2026-1-23 03:05
Heinz的参数可以用。UFF无法描述金属单质

Heinz的参数看起来比较大(描述Ag-Ag作用肯定合理)，导致Ag与其他原子比如C, O之间的作用会不会太强？
10号是Ag，前1-9分别是 C N O H Cl 之类，这些元素ε比较小。以下是参数：pair_coeff 1 1 0.066000 3.500000 pair_coeff 2 2 0.030000 2.500000 pair_coeff 3 3 0.170000 3.250000 pair_coeff 4 4 0.066000 3.500000 pair_coeff 5 5 0.0536 2.80 pair_coeff 6 6 0.1570 3.05 pair_coeff 7 7 0.71 4.02 pair_coeff 8 8 0.1553 3.166 pair_coeff 9 9 0.0 0.1 pair_coeff 10 10 4.56 2.6326
以下是混合参数：

(, 下载次数 Times of downloads: 5)

作者
Author: sobereva 时间: 2026-1-24 13:29

elpa 发表于 2026-1-23 21:36
Heinz的参数看起来比较大(描述Ag-Ag作用肯定合理)，导致Ag与其他原子比如C, O之间的作用会不会太强？

搞清楚LJ参数的混合规则

作者
Author: elpa 时间: 2026-1-25 10:35

sobereva 发表于 2026-1-24 13:29
搞清楚LJ参数的混合规则

用的lammps pair_modify mix arithmetic，ϵ对两者开平方根。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)