| 查看: 1277 | 回复: 0 | |||
[求助]
关于利用gamess中的spin flip tddft (sf-tddft)寻找conical intersection的问题
|
|
小弟初学gamess，希望利用其中的sf-tddft寻找势能曲线的圆锥交叉，参考的文献为J. Phys. Chem. A, 2009, 113 (46), pp 12749–12753,该文利用sf-tddft研究了乙烯S1和S0的CI,可是我没办法重复出文献中的结果。以下是我的输入文件： $CONTRL SCFTYP=UHF DFTTYP=BHHLYP TDDFT=SPNFLP MULT=3 RUNTYP=CONICAL EXETYP=RUN MAXIT=200 $END $SYSTEM TIMLIM=525600 MEMORY=6200000 $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 $END $guess guess=huckel $end $tddft nstate=5 iroot=1 tammd=.true. $end $CONICL OPTTYP=PENALTY IXROOT(1)=1,4 $END $STATPT NSTEP=1000 $END $scf dirscf=.t. $end $DATA Title C1 C 6.0 0.0000000000 0.0000000000 -0.6635189607 C 6.0 0.0000000000 0.0000000000 0.6635189607 H 1.0 0.9162439386 0.0000000000 -1.2297557119 H 1.0 -0.9162439386 0.0000000000 -1.2297557119 H 1.0 0.9162439386 0.0000000000 1.2297557119 H 1.0 -0.9162439386 0.0000000000 1.2297557119 $END 其中乙烯的结构是文献给的基态结构，IXROOT(1)=1,4是因为1，4分别对应S0 和S1，可是我的计算结果和文献对不上，只能得到一个两个CC键不断增长CH2相互远离的结构。 文献中计算细节的描述如下： The BHHLYP hybrid functional (50% Hartree-Fock plus 50% Becke exchange34 with Lee-Yang-Parr correlation35) was employed in this work because SFDFT benchmark calculations suggest better performance with a larger fraction of Hartree-Fock exchange.24 The basis sets employed were 6-31G(d,p) and aug-cc-pVTZ.36,37 The former was used in the MS-CASPT2 study of ref 21 and the latter in the MR-CI calculations in ref 15. The BroydenFletcher-Goldfarb-Shannon (BFGS) quasi-Newton scheme was adopted to minimize the objective function in eq 3. The same optimization criteria, thresholds, and parameters were used as discussed in ref 21, where three criteria were considered to achieve convergence, the change in the objective function f and the parallel and perpendicular components of the gradient of f with respect to the direction of the gradient vector of the penalty (the second term in eq 3). They were minimized simultaneously to be lower than a given threshold. No symmetry constraint was applied during the CI optimizations. 请各位前辈帮忙看一下我的问题究竟出在哪里，另外想问大家IXROOT后边括号里的数字有什么意义，手册里没有找到。还有就是通过观察输出文件我发现在中间步骤（NSEARCH）中第一和第四激发态已经不再对应S0和S1，并且出现了很多单重态和三重态混合的激发态（S SQUARE=1）请问这种情况怎样解释和解决。 另外附上了输出文件，请各位斧正。 问题很多，麻烦大家了，万分感谢！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : eth-spnflp-ex-conical.log
2015-06-19 19:22:38, 1.44 M
» 猜你喜欢
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有62人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有0人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
