计算双光子吸收截面时，基组如何输入

weixue

sob老师，您好！我在用dalton计算双光子吸收截面时，一直报错，其中里面有Ru(钌），C、N、O元素，应该怎么输入基组比较合适？

小范范1989

给个我们组做过例子：
ATOMBASIS                                                                     
using the Mixed basis                                                         
spin-orbit and phosphorescence                                                
Atomtypes=6 Nosymmetry Angstrom                                                
Charge=29.0 Atoms=2 Basis=lanl2dz_ecp  ECP=lanl2dz_ecp
Cu      -1.207406000      0.828055000      0.522245000
Cu       1.207515000     -0.827826000     -0.522247000
Charge=17.0 Atoms=2 Basis=6-31G*   
Cl      -3.454008000      0.582470000      0.939558000
Cl       3.454049000     -0.581893000     -0.939472000
Charge=15.0 Atoms=2 Basis=6-31G*   
P        0.344005000      0.616113000      2.176023000
P       -0.343877000     -0.616024000     -2.176046000
Charge=7.0 Atoms=2 Basis=6-31G*   
N        1.127027000     -1.807524000      1.307358000
N       -1.126920000      1.807660000     -1.307403000
Charge=6.0 Atoms=36 Basis=6-31G*   
C        0.776575000     -1.154932000      2.430044000
C        3.143103000      0.965786000      2.577830000
C        1.524004000     -3.089989000      1.376862000
C        4.315398000      1.711028000      2.498488000
C        1.959981000      1.456821000      2.015403000
C       -0.260920000      1.117290000      3.826534000

zjxitcc

基组那是最简单的了，除了自己手写，还可以用程序自动生成。例如用Gaussian算个几秒的任务

1. %chk=a.chk
   
2. #p HF genecp nosymm int=nobasistransform guess(only,save)
   
3. 
   
4. title
   
5. 
   
6. 0 1
   
7. [坐标]
   
8. 
   
9. Ru 0
   
10. SDD
    
11. ****
    
12. H C N O 0
    
13. def2SVP
    
14. ****
    
15. 
    

复制代码

交给高斯算几秒，完成后转化得到fch文件
formchk a.chk a.fch

从fch文件生成Dalton输入文件（.dal和.mol两个文件）
fch2dal a.fch

你可以在.dal文件中增加你要的关键词、删掉此时没用的轨道信息（因为高斯写的是only），基组信息就在.mol文件里。fch2dal小程序免费开源，可在https://gitlab.com/jxzou/mokit下载。

weixue

小范范1989 发表于 2021-11-3 08:49
给个我们组做过例子：
ATOMBASIS                                                                      ...

谢谢您的回复，我想再问一下，Cu的那个基组那里写ecp是什么？我写的Ru的 Basis=ecp-sdd-DZ可以吗

weixue

小范范1989 发表于 2021-11-3 08:49
给个我们组做过例子：
ATOMBASIS                                                                      ...

您好，老师。请问为什么在lanl2dz后面加ecp

小范范1989

weixue 发表于 2021-11-4 07:52
您好，老师。请问为什么在lanl2dz后面加ecp

赝势