计算化学公社

标题: 如何用过mc-pdft优化过交叉点？ [打印本页]
作者
Author: Freeman    时间: 2022-5-10 15:25
标题: 如何用过mc-pdft优化过交叉点？
大家好。我想用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

  12. &SEWARD

  14. &RASSCF
  15. Charge = 1
  16. NACTEL = 2
  17. RAS2 = 2
  18. CIRoot = 2 2 1
  19. fileorb=5ci+.mc.RasOrb

  21. >>> DoWhile

  23. &SEWARD

  25. &RASSCF
  26.   Charge = 1
  27.   NACTEL = 2
  28.   RAS2 = 2
  29.   CIRoot = 2 2 1
  30.   CMSI

  32. &MCPDFT
  33.   KSDFT=TPBE
  34.   Grad
  35.   MSPDft

  37. &SLAPAF

  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得到解答，就把解答复制到另一个平台去。


作者
Author: hebrewsnabla    时间: 2022-5-10 15:58
耐心一点，这不才一天没人回复么，会有人回的。或者你可以在issue里面艾特负责这块的开发者Jie Bao，或者发邮件给他。
作者
Author: zjxitcc    时间: 2022-5-10 16:02
我估计大家可能习惯了计算化学公社上的快速回复我看世界上大部分量化论坛，回复都没这里快。我在Dalton那里提的问题，几天能有人回复我就烧高香了，多的是几个月没回复的
作者
Author: Freeman    时间: 2022-5-10 16:05
zjxitcc 发表于 2022-5-10 16:02
我估计大家可能习惯了计算化学公社上的快速回复我看世界上大部分量化论坛，回复都没这里快。我在Dalto ...

是这么回事。我前几天收到了一个issue的回复，已经忘了是什么问题了，点进去一看才发现，哦，去年这个时候提的。
作者
Author: lijingbai2009    时间: 2022-5-11 05:03
通常CI优化都是要分别计算两个态的梯度，然后计算非绝热耦合常数。你这明显还在计算S0的梯度，后面的计算还没开始。

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

建议楼主测试一下并行与非并行的计算效率，你这个体系500多基函数，CAS（2，2），还只有两个态，40核并行有点杀鸡焉用牛刀了。。。
作者
Author: Freeman    时间: 2022-5-14 11:16
给大家贴一个模板，用最新的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

  12. &SEWARD

  14. &RASSCF
  15. Charge = 1
  16. NACTEL = 2
  17. RAS2 = 2
  18. CIRoot = 2 2 1
  19. fileorb=5ci+.mc.RasOrb

  21. >>> DoWhile

  23. &SEWARD
  24.   DOANalytical

  26. &RASSCF
  27.   Charge = 1
  28.   NACTEL = 2
  29.   RAS2 = 2
  30.   CIRoot = 2 2 1
  31.   CMSI
  32.   rlxr = 1

  34. &MCPDFT
  35.   KSDFT=T:PBE
  36.   Grad
  37.   MSPDft

  39. &ALASKA

  41. >>COPY $WorkDir/$Project.RunFile $WorkDir/RUNFILE2

  43. &RASSCF
  44.   Charge = 1
  45.   NACTEL = 2
  46.   RAS2 = 2
  47.   CIRoot = 2 2 1
  48.   CMSI
  49.   rlxr = 2

  51. &MCPDFT
  52.   KSDFT=T:PBE
  53.   Grad
  54.   MSPDft

  56. &ALASKA

  58. &SLAPAF

  60. >>> EndDo
复制代码








欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3