Molpro中equation of motion coupled cluster的input文件设置求助

Rainy_2020

各位老师，我是新手，想学习一下molpro中equation of motion coupled cluster(EOM-CC)的用法。现在准备从一个简单的 CuF6 cluster入手，Cu是2+离子，d轨道有9个电子，一共只有5个组态。在EOM-CC中，想算这个5个组态，不太清楚在molpro的input文件中具体应该怎么设置？试着用了下UCCSD(T), 并没有成功。
各位有对molpro中equation of motion coupled cluster的用法比较熟的吗？ 希望能够详细指教一下。

***, CuF6*****
MEMORY,1300,M;
DKROLL=2
Symmetry, NOSYM
ANGSTROM
GEOMETRY={
   7
1Cu-6F
Cu1    2.07073   -2.07073     0.00000
F2     2.07073   -2.07073    -1.96218
F2     2.07073   -2.07073     1.96218
F2    -0.181974  -2.07073     0.00000
F2     4.32343   -2.07073     0.00000
F2     2.07073   -0.181974    0.00000
F2     2.07073   -3.95948     0.00000
}

BASIS={
spdf, Cu1,aug-cc-pwCVTZ-DK; c;
spd, F2, aug-cc-pVDZ-DK; c;
}

TEXT, ~~~~~~~~~~~~~~~~
TEXT, Get the orbitals
TEXT, ~~~~~~~~~~~~~~~~
{ MATROP
  READ,SV,ORBITALS,,,FILE=o_hf_loc
  SAVE,SV,2101.2,ORBITALS
}
{ MERGE
  ORBITAL,2101.2
  
  MOVE,  1.1,   1.1,  1.1  ! Cu: 1s
  MOVE,  2.1,   2.1,  2.1  ! Cu: 2s
  MOVE,  3.1,   5.1,  3.1  ! Cu: 2p
  MOVE, 12.1,  12.1,  6.1  ! Cu: 3s
  MOVE, 13.1,  15.1,  7.1  ! Cu: 3p
  MOVE,  6.1,  11.1, 10.1  ! F: 1s [6]
  MOVE, 16.1,  21.1, 16.1  ! F: 2s [6]
  MOVE, 27.1,  44.1, 22.1  ! F: 2p [18]
  MOVE, 22.1,  26.1, 40.1  ! Cu: 3d [5]

  MOVE, 45.1, 292.1, 45.1  ! ...
  SAVE,2141.2
}

TEXT,~~~~~~~~~~~~~~~~~
TEXT,1.    5 CAS
TEXT,~~~~~~~~~~~~~~~~~
{ MCSCF
  OCC, 44; CLOSED, 39; frozen, 21,2141.2
  START,2141.2
  WF, 87,1, 1                        !! 2S=1
  state,    5
  MAXITER, 30
  PRINT, REF,REF1,CIVECTOR
  natorb,,ci;
  ORBITAL,2142.2
}
{ POP;
  INDIVIDUAL
}
{ LOCALI,PIPEK
  ORBITAL,2142.2
  LOCAO
  GROUP, 22.1, - 39.1
  GROUP, 40.1, - 44.1
}
{ MATROP
  LOAD, RO,ORB,2142.2,SET=3
  WRITE,RO,o_cas_5r new
}

PUT,    molden,o_cas_5rplot; NOSORT;
ORBITAL,2142.2,SET=3

{ UCCSD(T), thrden=1.00d-06, thrvar=1.00d-05
  ORBITAL, 2142.2, SET=3
  core
  direct
!  eom, -4.1, trans=1
  SAVE, 5130.2
}

sobereva

手册里搜EOM-CCSD，直接就有现成的例子

Rainy_2020

sobereva 发表于 2020-7-17 09:19
手册里搜EOM-CCSD，直接就有现成的例子

谢谢老师提醒！
我在open-shell cc计算中碰到这样一个报错信息：
1PROGRAM * CCSD (Unrestricted open-shell coupled cluster)   

NUMBER OF OCCUPIED ORBITALS NOT CONSISTENT WITH NUMBER OF ELECTRONS FOR CCSD. NELEC=  87  NCORE= 15  NCLOS= 24  NOCC= 44  MS2= 1

相应的UCCDS(T) input文件是：
{ uccsd(t), thrden=1.00d-06, thrvar=1.00d-05
  ORBITAL, 2142.2, SET=3
  WF, 87,1, 1
!!  eom, -5.1 !!trans=1
  SAVE, 5130.2
}

  请老师帮忙看看是哪里出了问题，谢谢啦                                 

beefly

molpro的开壳层cc不支持eom方法。如果是2018或以上版本，加一个电子变闭壳层，然后通过ip-ccsd算电离基态、激发态。
orca、cfour、qchem、gamess、nwchem在这方面做得更好

Rainy_2020

beefly 发表于 2020-7-17 18:14
molpro的开壳层cc不支持eom方法。如果是2018或以上版本，加一个电子变闭壳层，然后通过ip-ccsd算电离基态、 ...

谢谢老师！如果用ORCA的话，这里我还有一个问题是，在CASSCF中怎么优化inactive orbitals, 比如我在这个帖子里的molpro的input文件：
{ MCSCF
OCC, 44; CLOSED, 39; frozen, 21,2141.2
  START,2141.2
  WF, 87,1, 1                        
  state,    5
  MAXITER, 30
  PRINT, REF,REF1,CIVECTOR
  natorb,,ci;
  ORBITAL,2142.2
}
如果想在ORCA中优化inactive space中22-39的orbitals，在ORCA中具体应该怎么设置。我查看了ORCA的tuturial, 一直都没有尝试成功。还请老师指教呀。

Rainy_2020

beefly 发表于 2020-7-17 18:14
molpro的开壳层cc不支持eom方法。如果是2018或以上版本，加一个电子变闭壳层，然后通过ip-ccsd算电离基态、 ...

老师，想问一下，你知道ip-ccsd在molpro的input里面应该怎么设置吗？ 谢谢啦！

beefly

Rainy_2020 发表于 2020-7-20 09:18
老师，想问一下，你知道ip-ccsd在molpro的input里面应该怎么设置吗？ 谢谢啦！

手册里没写，需要问作者。建议用orca，更省事，功能也更多

Rainy_2020

beefly 发表于 2020-7-21 17:09
手册里没写，需要问作者。建议用orca，更省事，功能也更多

谢谢您，老师！ 我还有另外一个问题，在ORCA计算中优化orbitals, 计算中是不是默认优化所有的core, inactive和active orbitals？ 还是也能自己选择所要优化的obritals？