计算化学公社

标题: 关于Dalton软件里空腔半径的设置问题 [打印本页]

作者
Author: seantan521 时间: 2022-6-9 09:27
标题: 关于Dalton软件里空腔半径的设置问题
本帖最后由 seantan521 于 2022-6-9 09:39 编辑

各位老师和同行们好，我近期尝试用Dalton计算一个Pt磷光配合物的磷光辐射跃迁速率和激子衰减寿命，因为分子量比较大，必须先自定义空腔半径（不定义或者按照手册上默认方法写都会报错），但是几次下来都出现大量warnings，感觉应该是空腔半径定义的不合适（我给所有的非氢原子都分别定义了空腔半径，其中C:1.9, N:1.7, O:1.6, Pt:2.0）。想请各位比较熟悉Dalton操作的大佬们给些建议，先谢过各位了！
主要参考的资料如下：
Dalton 官方手册：https://daltonprogram.org/manuals/dalton2020manual.pdf
Dalton 官方论坛：http://forum.daltonprogram.org/
计算化学公社的帖子：http://bbs.keinsci.com/thread-2455-1-1.html
计算化学公社的帖子：http://bbs.keinsci.com/thread-13848-1-1.html
其他网上资料：https://cloud.tencent.com/developer/article/1668799

附上warnings提示如下：
200 WARNINGS have been issued by DALTON,
output from 'grep -n -i WARNING'  (max 10 warnings):
1244: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1248: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1253: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1257: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1264: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1267: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1273: ** WARNING ** A VERY POOR TESSELATION HAS BEEN CHOSEN
1280: WARNING: The distance between center of tessera  119and  311 is  0.000320, less than  0.000400 A
1375: Warning, element          46       872 of SI too big: set to zero
1376: Warning, element       119       308 of SI too big: set to zero

附上mol文件主要内容如下：
ATOMBASIS
test molecule
Generated by Multiwfn
Atomtypes=5 Angstrom Nosymmetry charge=0
Charge=1.0 Atoms=45 Basis=6-31G*
H1 4.81963025 2.01397617 0.02365220
往后省略若干
Charge=6.0 Atoms=54 Basis=6-31G*
C1 4.64519300 -0.12280483 0.12337778
往后省略若干
Charge=7.0 Atoms=3 Basis=6-31G*
N1 1.86498614 0.40194913 -0.03559965
N2 -2.05777264 0.00934545 -0.73407419
N3 -7.34391578 0.01804382 0.36807535
Charge=8.0 Atoms=1 Basis=6-31G*
O1 -0.01860170 2.57685031 -0.69703662
Charge=78.0 Atoms=1 Basis=ecp-sdd-DZ ECP=ecp-sdd-DZ
Pt1 -0.08465529 0.50023450 -0.42196961

附上dal文件主要内容如下：
**DALTON INPUT
.RUN RESPONS
*PCM
.SOLVNT
METHYLCL
.ICESPH /给特定的非H原子自定义半径
2
.NESFP    /确定非H原子的数目
59
.NEQRSP /溶剂对响应计算的非平衡贡献
*PCMCAV
.INA    /非H原子的.mol的顺序排列
1       /在.mol文件中的第一个原子
此处省略若干
X       /在.mol文件中的第X个原子
.RIN
1.9       /第一个原子的自定义半径
此处省略若干
2.0    /第X个原子的自定义半径
.AREATS /H原子的空腔半径为0.3
0.3
**WAVE FUNCTIONS
.DFT
B3LYP /在二次响应计算中目前只能使用B3LYP,CAM-B3LYP
*SCF INPUT
.MAX DIIS ITERATIONS
300
**RESPONSE
*QUADRATIC
.ECPHOS /使用赝势
.APROP /计算偶极矩A,即偶极矩算符
XDIPLEN
.APROP
YDIPLEN
.APROP
ZDIPLEN
.BPROP /即有效电荷自旋轨道算符
X1SPNSCA
.BPROP
Y1SPNSCA
.BPROP
Z1SPNSCA
.ISPABC
0 1 1
.ROOTS /计算的态的数目
1
**END OF INPUT

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