如何计算·CH2-CH2-CH2-CH2·的能量？

糖糖糖豆9988

一般的分子，都是先用便宜方法结构优化，然后在计算较高精度单点；对于双自由基，自旋极化单重态，按照sob老师“谈谈片段组合波函数与自旋极化单重态”、“CASSCF计算双自由基以及双自由基特征的计算”等文的介绍，我的理解是：第一步：# ub3lyp/6-31g(d) guess(fragment=2)得到自旋极化初猜波函数；第二步：# ub3lyp/6-31g(d) guess=(naturalorbitals,read,save,only)，即基于UNO确定双自由基特征；第三步：#P CAS(2,2)/6-31G* guess=read，得到双自由基的CASSCF能量。1）请问老师们，这样理解正确吗？2）上述第一至三步骤有没有体现双自由基结构的优化呢？
3）对称破缺波函数可以通过#P UHF/6-31G* guess=mix或者# ub3lyp/6-31g(d) guess(fragment=2)片段组合得到，但是 a）guess=mix方法使用的最初输入文件里·CH2-CH2-CH2-CH2·双自由基的结构、b)抽取两个片段的·CH2-CH2-CH2-CH2·结构分别是怎么得到呢？是直接从丁烷删掉俩H原子吗？
4）按照“详谈使用Gaussian做势能面扫描”一文中的“5.6 柔性扫描实例4：把环丁烷拉成丁烷双自由基”的介绍，使用# UM062X/6-31G* guess(mix,always) opt=modredundant pop=always nosymm，然后优化端点结构，也是计算单点能的一种方法，只不过是从环丁烷出发？

zjxitcc

（1）你的理解是正确的。但你描述的步骤并不是唯一做法。举个例子，例子见底下
（2）没有体现结构优化。
（3）你想要什么结构，直接在GaussView里编辑即可，鼠标点几点就行。
（4）既然叫“柔性扫描实例”，那就是柔性扫描，跟算单点能没啥关系。

例子在这里：写个gjf文件，用MOKIT中的automr提交

1. %mem=8GB
   
2. %nprocshared=4
   
3. #p CASSCF/cc-pVDZ
   
4. 
   
5. mokit{}
   
6. 
   
7. 0 1
   
8. C           -1.61451807    0.61952440    0.00000000
   
9. H           -1.25786364   -0.38928561    0.00000000
   
10. H           -2.68451807    0.61953758    0.00000000
    
11. C           -1.10117585    1.34548067    1.25740497
    
12. H           -1.45623220    2.35485376    1.25642745
    
13. H           -0.03117766    1.34377398    1.25838372
    
14. C           -1.61681550    0.62115040    2.51480761
    
15. H           -1.26206655   -0.38833093    2.51559568
    
16. H           -2.68681328    0.62318320    2.51401717
    
17. C           -1.10303147    1.34679310    3.77221313
    
18. H           -1.45771312    2.35629805    3.77138289
    
19. H           -0.03303386    1.34468902    3.77304506
    

复制代码

算之前不用指定活性空间，不用看轨道，算完后自然可以发现是CAS(2,2)，还有*_CASSCF_NO.fch文件，可以用GaussView或Multiwfn+VMD观看自然轨道和轨道占据数。