| 查看: 2640 | 回复: 14 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hcmzj668 的 2 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | |||
[交流]
怎么基于STEOM-DLPNO-CCSD计算rISC
|
|||
|
各位同行们，正常的td-dft计算risc (或isc)已按照orca手册成功，现在想计算基于steom-dlpno-ccsd的risc，不知道步骤和输入文件怎么编写，遂来求助。 手册中写道："theinput is similar to those discussed above. here esd(isc) should be used on themain input to indicate an intersystem crossing calculation, and the hessians should be provided byiscishessian and iscfshessian for the initial and final states, respectively. please note that the geometryused in the input file should correspond to that of the final state, specified through the iscfshessian flag. the relevantmatrix elements can be computed using any method available in orca.as a simple example, one could compute theexcited singlet and ground triplet geometries and hessians for anthracene usingtd-dft. then, compute the spin-orbit coupling (soc) matrix elements for aspecific triplet spin-sublevel using the same method (see the details below),potentially employing methods like casscf, mrci,steom-ccsd" 即：先得到gs，es和三重态的结构和hessians，然后基于gs结构做steom-dlpno-ccsd下的单线态和三重态的能量以用于dele。 但是，手册中没有给出具体的输入例子，就有点无从下手了，于是产生了如下问题和想法，麻烦各位老师帮忙建议。 1. “gs，es和三重态的结构和hessians”可以用td-dft方法得到，后直接使用？ 2. 基于上述得到的gs结构，做steom-dlpno-ccsd计算得到s1和t1的能量，是如下结果不呢？ steom-ccsdresults (singlets) iroot= 1: 0.109392 au 2.977 ev 24008.9 cm**-1 steom-ccsd results (triplets) iroot= 1: 0.089216 au 2.428 ev 19580.7 cm**-1 所以dele for risc be : 19580.7 - 24008.9 = -4,428.2 cm-1? 3. 最后，最重要的步骤计算risc: a)是用手册中提供的例子作为输入文件（td-dft），将上述td-dft得到的s1-t1结构和hessian， steom-dlpno-ccsd得到的能量，写入如下输入文件中来计算？ %tddft nroots 5 sroot 1 troot 1 trootssl 0 1 -1 #t1的三个亚能级 dosoc %esd iscishess "t1.hess" iscfshess "s1.hess" dele xxx # 用steom-dlpno-ccsd计算得到的 （应该是3个job，此处省略了） b)还是直接用td-dft得到的s1-t1结构和hessian，按照steom-dlpno-ccsd的输入文件计算？（这里不会写，尝试着写了如下输入文件，不知道对还是错） !steom-dlpno-ccsddef2-tzvp def2-tzvp/c def2/j tightscf tightpno esd(isc) %maxcore 5000 %scf maxiter 7000 end %mdci nroots 3 #如何控制t1的3个亚能级呢？（直接按照手册中的td-dft方式写3个job？） dorootwise true othresh 0.001 vthresh 0.001 tcutpnosingles 1e-11 ndav 400 dostoresteom true dosimpledens false addl2term true dtol 1e-5 dotriplet true steomsoc true end %esd iscishess "t1-opt.hess" iscfshess "s1-opt.hess" approxaden true # 这个 keyword会直接去自动获取dele不呢？ usej true doht false temp 298.15 %pal nprocs 24 end * xyzfile s0-opt.xyz 0 1 最后，麻烦各位行家了，祝大家生活愉快，科研顺利。 |
回复此楼» 猜你喜欢
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
钯炭、铂炭、钌炭等加氢催化剂
已经有27人回复
-
物理化学论文润色/翻译怎么收费?
已经有282人回复
-
AgI的FT-IR谱图
已经有1人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
【急招】 北京工业大学 第二批博士申请-杰青课题组1-2名额，25日截止！！
已经有12人回复
-
1，6-己二胺甲基化合成四甲基己二胺
已经有2人回复
-
江苏科技大学2026年招聘启事
已经有0人回复
-
求助电子版书籍《半导体材料计算》
已经有0人回复
-
气相色谱出峰异常求助
已经有1人回复
» 抢金币啦！回帖就可以得到:
-
江西理工大学2026年博士研究生招生报名公告（第二批次），可以联系马胜灿老师
+2/900
-
山东征女友，坐标济南
+1/168
-
半导体封测的“最后一公里”，为什么越来越需要“与世隔绝”的环境？
+1/89
-
测试█TEM/ EPR/ XPS/PY-GCMS/TG-IR/XRF/BET/MIP/核磁/EA/ICP，VX: 761711562。
+1/88
-
北京航空航天大学-仿生界面材料科学全国重点实验室郭林院士团队诚聘博士后
+1/83
-
江苏大学优青 胡兴好教授课题组博士后招聘
+1/79
-
西交利物浦大学招收26年12月入学奖学金博士生1名（人机协作康复的交互式数字孪生）
+1/34
-
大连工业大学杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
+1/33
-
【通知】北京信息科技大学仪器科学与光电工程学院招收调剂生（2026），快来报名吧！
+1/32
-
湖南师范大学（211，双一流）-有机化学-申请审核制博士招生
+1/31
-
【急招】“双一流”高校-新能源材料课题组招收2026年秋季入学博士生1名(湘潭大学)
+1/29
-
【通知】北京信息科技大学仪器科学与光电工程学院招收博士研究生（2026）
+2/22
-
【20260528】湖南工业大学许利剑教授招收2026级生物医学工程博士研究生
+1/12
-
“双一流”南京邮电大学--海优团队--招收2026博士生1-2名（6月15日前有效）
+1/12
-
【26.9月入学】211大学补招学博1人（AI或控制方向）
+1/9
-
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/7
-
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/5
-
上海技物所多维探测课题组博士后招聘
+1/4
-
深圳大学材料学院“智能分子材料团队”招聘博士后（连续流微反应器方向优先）
+1/3
-
北理工国家杰青团队招博士后
+1/2
hcmzj66813楼2024-11-26 19:33:33
2楼2024-11-25 13:53:13
3楼2024-11-25 14:27:07
11楼2024-11-26 11:53:04
12楼2024-11-26 11:57:07
14楼2024-11-26 22:40:43
简单回复
nono20094楼
2024-11-25 15:43
回复
hcmzj668(金币+1): 谢谢参与
·
2024-11-25 16:39
回复
hcmzj668(金币+1): 谢谢参与
庚辰辛丑6楼
2024-11-25 21:22
回复
hcmzj668(金币+1): 谢谢参与
psylhh7楼
2024-11-25 22:07
回复
hcmzj668(金币+1): 谢谢参与
tzynew8楼
2024-11-25 23:10
回复
hcmzj668(金币+1): 谢谢参与
lymjjh9楼
2024-11-26 08:34
回复
hcmzj668(金币+1): 谢谢参与
lymjjh10楼
2024-11-26 08:35
回复
月夜风波15楼
2024-11-29 15:13
回复
hcmzj668(金币+1): 谢谢参与
