求助gaussian加电场计算电偶极矩的问题

小蜗 · 发表于 Post on 2020-8-17 16:31:58

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

doublezhang · 发表于 Post on 2020-8-17 17:29:19

需要写field=read吗？有点怪怪的

wzkchem5 · 发表于 Post on 2020-8-17 18:18:03

算偶极矩，尤其是加了电场情况下的偶极矩，基组应该加弥散函数。如果基组不够弥散，基组就会限制电子云变形的程度，导致低估偶极矩的变化。

sobereva · 发表于 Post on 2020-8-17 20:24:21

field=read和field=x+300完全是冲突的，只能写一个，如果用前者，坐标末尾后头明确写明外场矢量。

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

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

小蜗 · 发表于 Post on 2020-8-17 21:50:40

谢谢sob老师~

小蜗 · 发表于 Post on 2020-8-17 21:51:29

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

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

小蜗 · 发表于 Post on 2020-8-17 21:52:22

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

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

hcy002 · 发表于 Post on 2020-8-17 21:52:38

请问，field=x+300是分子x方向加电场吗？ 0.003au相当于多少V?

小蜗 · 发表于 Post on 2020-8-18 09:39:26

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老师有讲哦~

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

[Gaussian/gview] 求助gaussian加电场计算电偶极矩的问题