如何用过mc-pdft优化过交叉点？

Freeman

大家好。我想用molcas的mc-pdft优化一个交叉点。molcas支持找交叉点，也支持mc-pdft的解析梯度，但是官网手册上只有用casscf优化交叉点的输入文件，没有mc-pdft的，我只能自己乱试。试过的输入文件如下：

1. &GATEWAY
   
2.   Coord = 5ci+.orca.xyz
   
3.   Basis
   
4.    def2-SVP
   
5.    ricd
   
6.   Constraints
   
7.     a = Ediff 1 2
   
8.   Values
   
9.     a = 0.0
   
10.   End of constraints
    
11. 
    
12. &SEWARD
    
13. 
    
14. &RASSCF
    
15. Charge = 1
    
16. NACTEL = 2
    
17. RAS2 = 2
    
18. CIRoot = 2 2 1
    
19. fileorb=5ci+.mc.RasOrb
    
20. 
    
21. >>> DoWhile
    
22. 
    
23. &SEWARD
    
24. 
    
25. &RASSCF
    
26.   Charge = 1
    
27.   NACTEL = 2
    
28.   RAS2 = 2
    
29.   CIRoot = 2 2 1
    
30.   CMSI
    
31. 
    
32. &MCPDFT
    
33.   KSDFT=TPBE
    
34.   Grad
    
35.   MSPDft
    
36. 
    
37. &SLAPAF
    
38. 
    
39. >>> EndDo

复制代码

不奏效。mc-pdft的解析梯度正常算完，但是&ALASKA后面，自动启动了对S0态的数值梯度计算，看起来没有在找交叉点，于是我就把它停了。输出文件：https://gitlab.com/Molcas/OpenMo ... 6e08ed/5ci++.mc.out
请问大家，我该怎么改写输入文件呢？
PS：我在openmolcas的gitlab主页issue区也提问了（https://gitlab.com/Molcas/OpenMolcas/-/issues/344），但至今都没得到回复。要是我在公社论坛或gitlab得到解答，就把解答复制到另一个平台去。

hebrewsnabla

耐心一点，这不才一天没人回复么，会有人回的。或者你可以在issue里面艾特负责这块的开发者Jie Bao，或者发邮件给他。

zjxitcc

我估计大家可能习惯了计算化学公社上的快速回复我看世界上大部分量化论坛，回复都没这里快。我在Dalton那里提的问题，几天能有人回复我就烧高香了，多的是几个月没回复的

Freeman

zjxitcc 发表于 2022-5-10 16:02
我估计大家可能习惯了计算化学公社上的快速回复我看世界上大部分量化论坛，回复都没这里快。我在Dalto ...

是这么回事。我前几天收到了一个issue的回复，已经忘了是什么问题了，点进去一看才发现，哦，去年这个时候提的。

lijingbai2009

通常CI优化都是要分别计算两个态的梯度，然后计算非绝热耦合常数。你这明显还在计算S0的梯度，后面的计算还没开始。

我注意到你使用了40个MPI线程，所以你的梯度计算是数值方法。Molcas的解析梯度不支持并行。说明书里说MCPDFT的Grad使用的是解析梯度，不知道会不会和数值梯度计算冲突。

建议楼主测试一下并行与非并行的计算效率，你这个体系500多基函数，CAS（2，2），还只有两个态，40核并行有点杀鸡焉用牛刀了。。。

Freeman

给大家贴一个模板，用最新的molcas运行是可以优化CI的（必须是最新，因为最近molcas的mc-pdft部分可能更新得很频繁）。但是非常非常慢：50来个原子，cas(2,2)，八小时才算2个点。不知道是咋回事儿，可能还需要测试一下不同进程、线程的并行方式。

1. &GATEWAY
   
2.   Coord = 5ci+.orca.xyz
   
3.   Basis
   
4.    def2-SVP
   
5.    ricd
   
6.   Constraints
   
7.     a = Ediff 1 2
   
8.   Values
   
9.     a = 0.0
   
10.   End of constraints
    
11. 
    
12. &SEWARD
    
13. 
    
14. &RASSCF
    
15. Charge = 1
    
16. NACTEL = 2
    
17. RAS2 = 2
    
18. CIRoot = 2 2 1
    
19. fileorb=5ci+.mc.RasOrb
    
20. 
    
21. >>> DoWhile
    
22. 
    
23. &SEWARD
    
24.   DOANalytical
    
25. 
    
26. &RASSCF
    
27.   Charge = 1
    
28.   NACTEL = 2
    
29.   RAS2 = 2
    
30.   CIRoot = 2 2 1
    
31.   CMSI
    
32.   rlxr = 1
    
33. 
    
34. &MCPDFT
    
35.   KSDFT=T:PBE
    
36.   Grad
    
37.   MSPDft
    
38. 
    
39. &ALASKA
    
40. 
    
41. >>COPY $WorkDir/$Project.RunFile $WorkDir/RUNFILE2
    
42. 
    
43. &RASSCF
    
44.   Charge = 1
    
45.   NACTEL = 2
    
46.   RAS2 = 2
    
47.   CIRoot = 2 2 1
    
48.   CMSI
    
49.   rlxr = 2
    
50. 
    
51. &MCPDFT
    
52.   KSDFT=T:PBE
    
53.   Grad
    
54.   MSPDft
    
55. 
    
56. &ALASKA
    
57. 
    
58. &SLAPAF
    
59. 
    
60. >>> EndDo
    

复制代码