计算化学公社

标题: 计算双光子吸收截面时，基组如何输入 [打印本页]
作者
Author: weixue    时间: 2021-11-3 08:14
标题: 计算双光子吸收截面时，基组如何输入
sob老师，您好！我在用dalton计算双光子吸收截面时，一直报错，其中里面有Ru(钌），C、N、O元素，应该怎么输入基组比较合适？
作者
Author: 小范范1989    时间: 2021-11-3 08:49
给个我们组做过例子：
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
作者
Author: zjxitcc    时间: 2021-11-3 09:33
本帖最后由 zjxitcc 于 2021-11-3 09:36 编辑

基组那是最简单的了，除了自己手写，还可以用程序自动生成。例如用Gaussian算个几秒的任务
  1. %chk=a.chk
  2. #p HF genecp nosymm int=nobasistransform guess(only,save)

  4. title

  6. 0 1
  7. [坐标]

  9. Ru 0
  10. SDD
  11. ****
  12. H C N O 0
  13. def2SVP
  14. ****

复制代码

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

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

你可以在.dal文件中增加你要的关键词、删掉此时没用的轨道信息（因为高斯写的是only），基组信息就在.mol文件里。fch2dal小程序免费开源，可在https://gitlab.com/jxzou/mokit下载。
作者
Author: weixue    时间: 2021-11-3 11:05
小范范1989 发表于 2021-11-3 08:49
给个我们组做过例子：
ATOMBASIS                                                                      ...

谢谢您的回复，我想再问一下，Cu的那个基组那里写ecp是什么？我写的Ru的 Basis=ecp-sdd-DZ可以吗
作者
Author: weixue    时间: 2021-11-4 07:52
小范范1989 发表于 2021-11-3 08:49
给个我们组做过例子：
ATOMBASIS                                                                      ...

您好，老师。请问为什么在lanl2dz后面加ecp
作者
Author: 小范范1989    时间: 2021-11-4 08:25
weixue 发表于 2021-11-4 07:52
您好，老师。请问为什么在lanl2dz后面加ecp

赝势



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