关于在Gaussian中使用nosymm关键词后体系朝向仍然发生变化

gingaorb

各位小哥哥，各位大佬们，我最近又遇到了一个奇怪的问题，想请教一下大家呀~
是这样的，最近我想做过渡态+电场方面的计算，然后我在没有电场时，读取.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怎么在这里突然失效了？到底是哪里出了问题鸭？



真是让我感到好困惑呀~
希望大家能够不吝赐教啊！谢谢大家~

sobereva

直接把坐标写在输入文件里，此时用了nosymm就绝对不可能有坐标整体的变化

王OldMan

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

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

sobereva

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

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

王OldMan

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）？

sobereva

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

nosymm

王OldMan

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

感谢sob老师