“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

[VASP] 求助:AIMD中使用slow growth方法时,CV取配位数时,如何原子间平衡j键长?

[复制链接 Copy URL]

257

帖子

0

威望

1151

eV
积分
1408

Level 4 (黑子)

如题,我想用隐式溶剂模型计算一下金属锌(002)面一个锌原子的溶出能,打算采用slow growth的方法跑一条10 ps的AIMD轨迹,配位数打算定义为选取锌原子与其他所有锌原子的配位数,但是对于平衡键长的计算,我没什么概念,应该如何计算或者如何定义?请各位大神不吝赐教,给些许点播,不胜感谢!
努力挖掘科学的本质

529

帖子

1

威望

2463

eV
积分
3012

Level 5 (御坂)

2#
发表于 Post on 2021-11-15 08:05:06 | 只看该作者 Only view this author
我不清楚你说的平衡键长具体指什么。如果是switching function中指数项里面的分子上的那个长度的话,我曾经在 https://aip.scitation.org/doi/10.1063/5.0020240 的Chemical reaction in vacuum中看到调整这个东西的操作。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
Penson + 5 谢谢分享

查看全部评分 View all ratings

257

帖子

0

威望

1151

eV
积分
1408

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-11-15 23:58:15 | 只看该作者 Only view this author
Daniel_Arndt 发表于 2021-11-15 08:05
我不清楚你说的平衡键长具体指什么。如果是switching function中指数项里面的分子上的那个长度的话,我曾经 ...

非常感谢,应该是指数项里面的分母吧?以下公式中的d?那确实是我想要的东西
努力挖掘科学的本质

257

帖子

0

威望

1151

eV
积分
1408

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2021-11-17 17:30:43 | 只看该作者 Only view this author
接着这个帖子追问一下,SG计算跑起来后,无论怎样调整配位数设置的平衡键长,都会出现下面的错误,请问有大佬知道原因吗?附上我的输入文件以供参考

POTCAR

576.85 KB, 下载次数 Times of downloads: 9

POSCAR

5.31 KB, 下载次数 Times of downloads: 27

KPOINTS

77 Bytes, 下载次数 Times of downloads: 15

INCAR

2.17 KB, 下载次数 Times of downloads: 63

ICONST

522 Bytes, 下载次数 Times of downloads: 91

努力挖掘科学的本质

97

帖子

0

威望

1167

eV
积分
1264

Level 4 (黑子)

5#
发表于 Post on 2021-12-17 10:29:03 | 只看该作者 Only view this author
Penson 发表于 2021-11-17 17:30
接着这个帖子追问一下,SG计算跑起来后,无论怎样调整配位数设置的平衡键长,都会出现下面的错误,请问有大 ...

函数作图,一清二楚呀https://zuotu.91maths.com/#W3sid ... WQiOlsiMSIsIjEiXX1d
man

257

帖子

0

威望

1151

eV
积分
1408

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-12-19 11:09:12 | 只看该作者 Only view this author
elpa 发表于 2021-12-17 10:29
函数作图,一清二楚呀https://zuotu.91maths.com/#W3sidHlwZSI6MCwiZXEiOiIoMS0oeC8xLjUpXjgpLygxLSh4LzE ...

请问您给的函数是对应我这个体系的?如何确认d, n, m这些参数的呢?
努力挖掘科学的本质

239

帖子

0

威望

5297

eV
积分
5536

Level 6 (一方通行)

7#
发表于 Post on 2021-12-20 16:51:16 | 只看该作者 Only view this author
截断半径可以取一个比两种原子的共价半径之和稍大的值。还可以看看相关文献是怎么设置的。
我自己的做法是先用距离为cv跑一下,看看自由能势垒的位置,把它作为截断半径

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
Penson + 5 谢谢分享

查看全部评分 View all ratings

256

帖子

3

威望

1922

eV
积分
2238

Level 5 (御坂)

8#
发表于 Post on 2022-6-18 02:02:56 | 只看该作者 Only view this author
我有个问题请教您下,VASP跑的slow-growth的时候ICONST中配位的写法是怎样的,我看您写完平衡键长后最后的status值设置为5,我看下官网,status值设置为5不是跑Metadynamics的时候的设置吗?

257

帖子

0

威望

1151

eV
积分
1408

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-6-19 00:26:00 | 只看该作者 Only view this author
wangyueda 发表于 2022-6-18 02:02
我有个问题请教您下,VASP跑的slow-growth的时候ICONST中配位的写法是怎样的,我看您写完平衡键长后最后的s ...

应该按照官网是比较可靠的
努力挖掘科学的本质

本版积分规则 Credits rule

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

GMT+8, 2026-4-21 06:00 , Processed in 0.588633 second(s), 27 queries , Gzip On.

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