计算化学公社

标题: Molpro中equation of motion coupled cluster的input文件设置求助 [打印本页]

作者
Author: Rainy_2020 时间: 2020-7-17 01:39
标题: Molpro中equation of motion coupled cluster的input文件设置求助
本帖最后由 Rainy_2020 于 2020-7-17 05:55 编辑

各位老师，我是新手，想学习一下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
}

作者
Author: sobereva 时间: 2020-7-17 09:19
手册里搜EOM-CCSD，直接就有现成的例子

作者
Author: Rainy_2020 时间: 2020-7-17 12:14

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
}

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

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

作者
Author: Rainy_2020 时间: 2020-7-17 22:52

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, 一直都没有尝试成功。还请老师指教呀。

作者
Author: Rainy_2020 时间: 2020-7-20 09:18

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

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

作者
Author: beefly 时间: 2020-7-21 17:09

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

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

作者
Author: Rainy_2020 时间: 2020-7-22 03:32

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

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

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