计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 如何用CP2K计算水空气界面处的电场大小

如何用CP2K计算水空气界面处的电场大小

查看数: 2848 | 评论数: 10 | 收藏 Add to favorites 8

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

gingaorb

发布时间: 2022-8-25 22:46

正文摘要:

各位小哥哥，各位大佬们，请教一个问题呀~ 是这样的，我的研究课题是关于水空气界面的分子动力学模拟，我花了两个月时间的摸索，才终于用256个水分子建立的水层用CP2K跑出来一个NVT的平衡构型，像这篇文献中所展 ...

回复 Reply

sobereva 发表于 Post on 2023-3-17 22:02:55

chen0201 发表于 2023-3-17 20:27
您好，我想请问下，水溶液中，不需要外加电势，也有电场吗？

体系里总是有电场，因为静电势并不是处处相同。和加不加外势没关系

chen0201 发表于 Post on 2023-3-17 20:27:56

您好，我想请问下，水溶液中，不需要外加电势，也有电场吗？

sobereva 发表于 Post on 2022-9-27 06:05:34

gingaorb 发表于 2022-9-26 17:57
我来回帖啦！之前关于CP2K输出的电场强度的单位问题解决了~确实是Hartree/bohr/e~参考站长大大的帖子[/ba ...

&MD里设COMVEL_TOL 0消除平动试试

sobereva 发表于 Post on 2022-9-14 03:20:42

gingaorb 发表于 2022-9-13 21:56
我来回帖啦！之前在站长大大的指导下，我成功用Multiwfn绘制出了电场在z方向的平面平均值的图像，虽然花 ...

这种问题，拿个简单体系比如氟化氢，用CP2K和Gaussian（ORCA等亦可）都计算，看CP2K的值在什么单位时大致能对上Multiwfn基于Gaussian的fch文件（或ORCA的molden文件，等等）算的电场便知。Multiwfn载入波函数文件后在主功能3里输入坐标就能得到a.u.为单位的静电势，手动做有限差分求导数就能得到电场。

gingaorb 发表于 Post on 2022-8-26 18:47:11

sobereva 发表于 2022-8-26 14:47
CP2K直接就能输出各个位置电场的格点数据：https://manual.cp2k.org/cp2k-202 ... NT/EFIELD_CUBE.html

...

太感谢站长大大了！我先按您说的试试看！

sobereva 发表于 Post on 2022-8-26 14:47:12

CP2K直接就能输出各个位置电场的格点数据：https://manual.cp2k.org/cp2k-202 ... NT/EFIELD_CUBE.html

Multiwfn创建CP2K输入文件的功能（http://sobereva.com/587）里选-3 Set exporting cube file也可以直接设置要求产生电场的cube文件

CP2K给的是三维空间中各个点的，可以再用Multiwfn获得不同z坐标的平面平均的值，参考
使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线
http://sobereva.com/638（http://bbs.keinsci.com/thread-28225-1-1.html）

Jack 发表于 Post on 2022-8-26 07:10:57

文章中给了公式了。对电势微分就是电场了

gingaorb 发表于 Post on 2022-8-25 22:48:16

这里附上我所参考的这篇文章的全文~

acs.jpclett.7b01257.pdf

3.11 MB, 下载次数 Times of downloads: 70

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register