计算化学公社

标题: 求助gaussian加电场计算电偶极矩的问题 [打印本页]

作者
Author: 小蜗 时间: 2020-8-17 16:31
标题: 求助gaussian加电场计算电偶极矩的问题
各位老师好，最近刚接触了gaussian，想用电场计算一个分子的偶极矩变化情况。电场关键词是#p opt b3lyp/6-31g(d) nosymm field=read field=x+300，预期目标是分子在电场下的偶极矩会随电场增大或缩小，甚至会出现偶极矩方向翻转的现象（文献中有类似情况）但是，自己计算发现，加0.003au的电场，偶极矩确实会改变，之后再加多大的电场偶极矩都为恒定值。想问问各位老师，是哪里出了问题，非常感谢！

作者
Author: doublezhang 时间: 2020-8-17 17:29
需要写field=read吗？有点怪怪的

作者
Author: wzkchem5 时间: 2020-8-17 18:18
算偶极矩，尤其是加了电场情况下的偶极矩，基组应该加弥散函数。如果基组不够弥散，基组就会限制电子云变形的程度，导致低估偶极矩的变化。

作者
Author: sobereva 时间: 2020-8-17 20:24
field=read和field=x+300完全是冲突的，只能写一个，如果用前者，坐标末尾后头明确写明外场矢量。

不可能继续加大电场完全不影响偶极矩，只是程度问题

如果你对偶极矩定量数值有要求，当前基组完全没法用，最起码得def2-SVPD，仔细看此文
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

作者
Author: 小蜗 时间: 2020-8-17 21:50
谢谢sob老师~

作者
Author: 小蜗 时间: 2020-8-17 21:51

doublezhang 发表于 2020-8-17 17:29
需要写field=read吗？有点怪怪的

啊好像是我写的有问题耶刚接触还是要再学习的感谢你哇

作者
Author: 小蜗 时间: 2020-8-17 21:52

wzkchem5 发表于 2020-8-17 18:18
算偶极矩，尤其是加了电场情况下的偶极矩，基组应该加弥散函数。如果基组不够弥散，基组就会限制电子云变形 ...

好的好的我会记住你的建议谢谢哇

作者
Author: hcy002 时间: 2020-8-17 21:52
请问，field=x+300是分子x方向加电场吗？ 0.003au相当于多少V?

作者
Author: 小蜗 时间: 2020-8-18 09:39

hcy002 发表于 2020-8-17 21:52
请问，field=x+300是分子x方向加电场吗？ 0.003au相当于多少V?

1a.u.电场
=1 Hartree/bohr/e
=(2625500J/mol)/0.529177Angstrom/(1.602176*10^-19C)
=2625500/6.022^23/0.529177/1.602176*10^19 J/C/Angstrom
=51.423 J/C/Angstrom
=51.423 V/Angstrom
sob老师有讲哦~

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