计算化学公社

标题: 请教一下3个中性的有机分子（DOL）+一个Li原子高斯计算的自旋量子数设置 [打印本页]

作者
Author: yuda_li 时间: 2024-3-1 15:27
标题: 请教一下3个中性的有机分子（DOL）+一个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

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

作者
Author: zjxitcc 时间: 2024-3-1 15:51
本帖最后由 zjxitcc 于 2024-3-1 15:54 编辑

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

作者
Author: yuda_li 时间: 2024-3-1 16:09

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

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

作者
Author: zjxitcc 时间: 2024-3-1 16:14
本帖最后由 zjxitcc 于 2024-3-1 16:15 编辑

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

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

作者
Author: yuda_li 时间: 2024-3-1 16:20

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

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

作者
Author: yuda_li 时间: 2024-3-1 16:21
@zjxitcc

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