计算化学公社

标题: 关于在Gaussian中使用nosymm关键词后体系朝向仍然发生变化 [打印本页]

作者
Author: gingaorb 时间: 2024-6-30 23:07
标题: 关于在Gaussian中使用nosymm关键词后体系朝向仍然发生变化
各位小哥哥，各位大佬们，我最近又遇到了一个奇怪的问题，想请教一下大家呀~
是这样的，最近我想做过渡态+电场方面的计算，然后我在没有电场时，读取.fchk文件，进行偶极矩方向的确认，如图1所示的那样（并且也通过输出的.log文件得到了确认）：
Dipole moment (field-independent basis, Debye):
X= -1.6635 Y= 1.4643 Z= -2.3394 Tot= 3.2225

然后呢~为了保险起见，我使用如下关键词，读取这个.fchk对应的chk文件进行了又一遍的计算：
%oldchk=La-1H2O-TS.chk
%nprocshared=64
%mem=60GB
%chk=TS0
# opt=(ts, calcfc, noeigen) freq b3lyp/def2tzvp empiricaldispersion=gd3bj pop=(full,nbo) nosymm geom=allcheck guess=read field=read

0.0000 0.0000 0.0000

计算得到过渡态结构后，诡异的事情发生了！！！得到其偶极矩变为：（1.3289，-1.1805，-2.6911）和3.2251 Debye。Gaussian它居然又把结构旋转了！！！虽然偶极矩的大小没有明显变化，但X、Y的方向居然反向了！！！为什么鸭？明明我都已经用了“nosymm”关键词了~nosymm怎么在这里突然失效了？到底是哪里出了问题鸭？

真是让我感到好困惑呀

~
希望大家能够不吝赐教啊！谢谢大家~

作者
Author: sobereva 时间: 2024-7-1 18:21
直接把坐标写在输入文件里，此时用了nosymm就绝对不可能有坐标整体的变化

作者
Author: 王OldMan 时间: 2024-9-10 16:24

sobereva 发表于 2024-7-1 18:21
直接把坐标写在输入文件里，此时用了nosymm就绝对不可能有坐标整体的变化

您好，老师。我没理解您说的意思，按我的理解输入文件中(.gjf)不就是有包含原子的坐标了吗？您说的“直接把坐标写在输入文件”是什么意思？该如何操作？
目前还是小白，如果问的比较脑瘫请不吝赐教

作者
Author: sobereva 时间: 2024-9-10 22:55

王OldMan 发表于 2024-9-10 16:24
您好，老师。我没理解您说的意思，按我的理解输入文件中(.gjf)不就是有包含原子的坐标了吗？您说的“直 ...

楼主用了geom=allcheck，输入文件里没原子坐标

作者
Author: 王OldMan 时间: 2024-9-11 09:32

sobereva 发表于 2024-9-10 22:55
楼主用了geom=allcheck，输入文件里没原子坐标

好的，谢谢老师。我遇到的问题是：我先用关键词（# opt freq b3lyp/6-31g(d) em=gd3bj nosymm field=z-40）对分子结构进行优化，后面我再用关键词（# m062x/def2TZVP field=z-40 int=ultrafine）计算结构优化后分子的单点能。计算单点能后发现，分子的坐标发生了翻转（即各原子的坐标关于x-y平面发生翻转），通过查计算单点能的输出文件，发现分子的（Standard orientation）坐标发生了变化。请问老师坐标变化的问题该如何解决？是否计算分子在电场条件下的单点能时也应加上关键词（nosymm）？

作者
Author: sobereva 时间: 2024-9-11 11:44

王OldMan 发表于 2024-9-11 09:32
好的，谢谢老师。我遇到的问题是：我先用关键词（# opt freq b3lyp/6-31g(d) em=gd3bj nosymm field=z-40 ...

nosymm

作者
Author: 王OldMan 时间: 2024-9-11 15:10

sobereva 发表于 2024-9-11 11:44
nosymm

感谢sob老师

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