24小时热门版块排行榜    

查看: 1645  |  回复: 3
【悬赏金币】回答本帖问题,作者Han_ping将赠送您 5 个金币

Han_ping

新虫 (初入文坛)

[求助] 用TDDFT如何进行激发态结构优化以及相关参数设置

大家好!我最近在尝试使用Turbomole优化结构的S1激发态,遇到优化失败的问题。我不确定是我control文件中没有设置$exopt 1(看manual说的是$soes 1的时候,无需设置$exopt 1)还是我的提交脚本设置错误?下面是我的control文件和提交脚本,希望各位大侠不吝赐教! 以及我不知道如何设置参数$exopt
$title
$symmetry c1
$redundant    file=coord
$user-defined bonds    file=coord
$coord    file=coord
$optimize
internal   on
redundant  on
cartesian  off
global     off
basis      off
$atoms
cd 1-13                                                                        \
   basis =cd def2-TZVP                                                         \
   ecp   =cd def2-ecp                                                          \
   jbas  =cd def2-TZVP
se 14-23                                                                       \
   basis =se def2-TZVP                                                         \
   jbas  =se def2-TZVP
h  24-42                                                                       \
   basis =h def2-TZVP                                                          \
   jbas  =h def2-TZVP                                                          \
   charge=  1.50000000
h  43-51                                                                       \
   basis =h def2-TZVP                                                          \
   jbas  =h def2-TZVP                                                          \
   charge=  0.50000000
cl 52                                                                          \
   basis =cl def2-TZVP                                                         \
   jbas  =cl def2-TZVP
o  53-54                                                                       \
   basis =o def2-TZVP                                                          \
   jbas  =o def2-TZVP
c  55-57                                                                       \
   basis =c def2-TZVP                                                          \
   jbas  =c def2-TZVP
h  58-61  
   basis =h def2-TZVP                                                          \
   jbas  =h def2-TZVP
$basis    file=basis
$ecp    file=basis
$scfmo   file=mos
$scfiterlimit      500
$thize     0.10000000E-04
$thime        5
$scfdamp   start=1.000  step=0.050  min=0.100
$scfdump
$scfintunit
unit=30       size=0        file=twoint
$scfdiis
$scforbitalshift  automatic=.1
$drvopt
   cartesian  on
   basis      off
   global     off
   hessian    on
   dipole     on
   nuclear polarizability
$interconversion  off
   qconv=1.d-7
   maxiter=25
$coordinateupdate
   dqmax=0.3
   interpolate  on
   statistics    5
$forceupdate
   ahlrichs numgeo=0  mingeo=3 maxgeo=4 modus=<g|dq> dynamic fail=0.3
   threig=0.005  reseig=0.005  thrbig=3.0  scale=1.00  damping=0.0
$forceinit on
   diag=default
$energy    file=energy
$grad    file=gradient
$forceapprox    file=forceapprox
$dft
   functional pbe
   gridsize   m4
$scfconv   7
$ricore      500
$rij
$jbas    file=auxbasis
$marij
$rundimensions
   natoms=61
   nbf(CAO)=1562
   nbf(AO)=1384
$denconv 1d-7
$last SCF energy change = 0.97039755E-02
$charge from ridft
          0.000 (not to be modified here)
$dipole from ridft
  x     2.16809595588296    y    -1.19016599840967    z     2.20284027328950    a.u.
   | dipole | =    8.4184416844  debye
$last excitation energy change= -.67274572E-02
$optinfo       file=optinfo
$hessapprox   file=hessapprox
$subenergy  Etot         E1                  Ej                Ex                 Ec                 En
-26949.73103244    -83619.59738273     35160.31172441    -1114.672858209    -30.95467005923     22655.18215414
$restart dscf
$closed shells
a       1-344                                  ( 2 )
$scfinstab rpas
$soes
a            1
$rpacor      274   in_total
#$maxcor    500 MiB  per_core
$last step     define
$end
提交命令是:jobex -ex > job.out

用TDDFT如何进行激发态结构优化以及相关参数设置
wrong.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ohnesie

金虫 (正式写手)

建议用薛定谔软件的 jaguar模块试试,很简单

发自小木虫Android客户端
微信:FreiheitJun
2楼2023-02-11 22:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Han_ping

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by ohnesie at 2023-02-11 22:38:51
建议用薛定谔软件的 jaguar模块试试,很简单

好的,谢谢,我不想更换软件了,我在摸索一下试一试
3楼2023-02-12 22:10:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

季秋十八

新虫 (初入文坛)

你好你现在问题解决了吗,我最近使用turbomole也有些问题,可以讨论下吗

发自小木虫IOS客户端
4楼2023-03-29 15:52:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Han_ping 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 272求调剂 +3 材紫有化 2026-02-28 3/150 2026-02-28 22:52 by ms629
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 290求调剂 +5 材料专硕调剂; 2026-02-28 6/300 2026-02-28 21:40 by gaoxiaoniuma
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见