请教一下3个中性的有机分子（DOL）+一个Li原子高斯计算的自旋量子数设置

yuda_li

输入文件结构是3个中性的有机分子（DOL）+一个Li原子

想计算的是体系在整体为0价、正1、负1三种情况下的吉布斯自由能，以下是输入文件：

%chk=myfile.chk
# gfinput b3lyp/6-311G(d,p) Freq

pvdf-hfp

0   1
O   18.669249   35.297791   39.718655
O   18.589733   34.840206   37.606609
C   19.429352   36.357082   39.080540
C   19.060286   36.166855   37.604740
C   17.801764   34.700108   38.738522
H   20.517099   36.271900   39.352039
H   18.957983   37.214581   39.563713
H   18.324253   36.881969   37.238945
H   19.882317   36.157974   36.835403
H   17.630344   33.679489   38.979702
H   16.869658   35.374676   38.629623
O   19.499519   33.631874   42.307442
O   20.331427   31.765799   41.218410
C   19.165506   32.516785   43.121120
C   19.426548   31.298868   42.254803
C   20.673910   33.084965   41.609737
H   18.121519   32.631180   43.459805
H   19.723637   32.547459   44.051411
H   19.992214   30.490099   42.698536
H   18.573257   30.935469   41.723316
H   20.911814   33.660210   40.662735
H   21.547939   33.148064   42.318134
O   18.620348   37.831078   44.919422
O   18.286861   35.956467   43.750118
C   19.024178   36.763309   45.796612
C   18.579334   35.539959   45.085232
C   17.766298   37.223698   43.963306
H   18.579159   36.783058   46.810627
H   20.086744   36.682159   45.884830
H   19.452723   34.852772   44.979309
H   17.649275   35.123100   45.436657
H   16.653372   37.237076   44.311810
H   17.617668   37.748989   43.038330
O   22.377625   37.449097   41.990139
O   20.159716   37.233685   41.856403
C   21.936403   38.788185   41.751125
C   20.413395   38.630497   41.846825
C   21.264406   36.799183   42.553932
H   22.260193   39.060501   40.734657
H   22.394920   39.348068   42.521450
H   20.051514   39.162277   42.759823
H   19.911413   39.117989   41.056339
H   21.321312   35.697475   42.478741
H   21.019491   37.013775   43.627544
Li   18.399912   35.318062   41.873756


这里我设置了价态为0，自旋量子数为1
但是最终结果报错为：
The combination of multiplicity 1 and   163 electrons is impossible.
Error termination via Lnk1e in /home/admin/opt/gaussian/g16/l301.exe at Thu Feb 29 12:19:10 2024.
Job cpu time:       0 days  0 hours  0 minutes  0.5 seconds.
Elapsed time:       0 days  0 hours  0 minutes  0.5 seconds.
File lengths (MBytes):  RWF=     17 Int=      0 D2E=      0 Chk=      2 Scr=      2


看到其他老师解释过这个问题，可能是建模的时候少了一个氢原子，但是我检查过了，输入文件是没错的，请教一下各位老师应该怎么办？

zjxitcc

离合理的计算还有很大距离。
（1）你还没做结构优化（即使做了结构优化，还要确认没有虚频）。直接freq算出来的热力学量没有意义。
（2）你想说的是自旋多重度为1，不是自旋量子数为1。
（3）Li原子是二重态的，含有奇数个电子，其他分子都是偶数个电子，所以整体是奇数个电子，当前体系自旋多重度应为2。
（4）对于绝大多数普通用户而言，关键词gfinput对他们没有任何用处，没必要写。
（5）没有加色散校正，急需仔细阅读
《谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的》http://sobereva.com/557
《谈谈“计算时是否需要加DFT-D3色散校正？”》http://sobereva.com/413

yuda_li

zjxitcc 发表于 2024-3-1 15:51
离合理的计算还有很大距离。
（1）你还没做结构优化（即使做了结构优化，还要确认没有虚频）。直接freq算 ...

感谢老师回复，再请教一下老师：（3）中，如果整体是+1价和-1价的话，自旋多重度应该怎么设置呢？

zjxitcc

yuda_li 发表于 2024-3-1 16:09
感谢老师回复，再请教一下老师：（3）中，如果整体是+1价和-1价的话，自旋多重度应该怎么设置呢？

（1）若当前体系整体为中性，则当前体系自旋多重度为2，那么整体是+1价和-1价时自旋多重度都是1。
（2）实际上当前体系整体不太可能为中性，一般来说算这种体系的人都有实际应用背景，他们一般碰到的都是Li+阳离子，所以整体带+1正电荷，自旋多重度为1。而整体得、失一个电子的话，对应的是中性和带+2电荷，自旋多重度都是2。
我不知道你的实际研究问题背景，所以给了两种答案，希望你自己知道自己在算什么。

yuda_li

zjxitcc 发表于 2024-3-1 16:14
（1）若当前体系整体为中性，则当前体系自旋多重度为2，那么整体是+1价和-1价时自旋多重度都是1。
（2）实 ...

谢谢老师的（1）的回答，（2）的话，我是看到别人论文里的公式以及描述，他也是3个中性分子包裹一个Li原子，实际上是在模拟锂电池溶液里的溶剂化体系，您可以看一下我给你发的这个原文，如果我没理解错的话，文章中的GN描述的应该就是中性体系，希望老师指点！原文我放在评论下了，回复好像没办法发图片

yuda_li

@zjxitcc