| 查看: 4643 | 回复: 8 | |||
shenzhendsh铁虫 (小有名气)
|
[求助]
Gaussian TDDFT 激发态 频率计算 输入文件
|
|
我要计算单重激发态的频率，输入文件如下： %chk=chk.chk %nprocshared=8 %mem=3200MB # freq td=(nstates=30,root=1,singlet) b3lyp tzvp title Cartesians coordination 请问为何一直算不完，而且有很多SCF Done能量？正确的输入命令是什么呢？ 对于激发态过渡态的优化，输入命令又是什么呢？ 刚接触激发态计算，还望大家不吝赐教！非常感谢！ |
回复此楼» 猜你喜欢
莆田学院环境与生物工程学院资源与环境专硕（环境工程方向） 开放调剂招生
已经有12人回复
-
浙江海洋大学石油化工与环境学院招收工科调剂生（学硕、专硕均可，线上面试）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有245人回复
-
浙江海洋大学石油化工学院接收工科调剂（学硕、专硕均可，近期线上面试）
已经有8人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
浙江海洋大学石油与天然气工程专业接收工科调剂（13日早上截止，线上面试）
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
考研调剂～湖南工程学院材料与化工学院
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高精度量子化学计算软件Molpro 入门及理论精讲课程
已经有3人回复
- Gaussian 03优化结构
已经有12人回复
- gaussian保持对称性优化
已经有8人回复
- 求助有关高斯计算激发态能量准确度问题
已经有6人回复
- gaussian 高性能计算问题
已经有10人回复
- 耦合簇CCSD(T)能否用于计算激发态
已经有4人回复
- 请问在Gaussian中如何固定结构然后再算单点能
已经有6人回复
- Gaussian 计算拉曼光谱最后时刻出错
已经有14人回复
- 求助：gauss09优化H2出错
已经有5人回复
- 【讨论】molpro, molcas和gaussian的优缺点
已经有3人回复
- 【coolrainbow个人文集】从手算理解量子化学：H2体系的量子化学计算（pdf+程序）
已经有65人回复
- 【求助】在gaussian过渡态优化计算出现的问题~请教高手~
已经有3人回复
- 【求助】Gaussian SCRF 优化求助
已经有7人回复
- 【求助】用TDDFT优化出来激发态之后的布居分析
已经有20人回复
- 【求助】Gaussian09作业中断，可否接着断点继续计算？
已经有11人回复
- 【求助】G09使用TDDFT优化激发态结构 输入怎么写
已经有10人回复
- 【求助】咨询单点能计算采用优化计算时得出的结构和基组不能重复结果的问题
已经有3人回复
- 【求助】TD计算激发态时Davidson循环不过
已经有9人回复
- 【求助】单重激发态的高斯算法
已经有29人回复
- 【求助】Gaussian09 优化激发态出错
已经有8人回复
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
2楼2013-11-12 16:21:39
shenzhendsh
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3.3
- 散金: 28
- 帖子: 61
- 在线: 68.3小时
- 虫号: 1843099
- 注册: 2012-06-01
- 性别: GG
- 专业: 高分子物理与高分子物理化
3楼2013-11-12 16:27:56
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
4楼2013-11-12 16:48:48
pwzhou
铁杆木虫 (正式写手)
- QC强帖: 2
- 应助: 51 (初中生)
- 金币: 7551.4
- 散金: 13
- 红花: 14
- 帖子: 412
- 在线: 316.9小时
- 虫号: 56225
- 注册: 2005-01-07
- 性别: GG
- 专业: 化学动力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-11-14 14:27:24
shenzhendsh: 金币+5 2013-11-17 16:27:42
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-11-14 14:27:24
shenzhendsh: 金币+5 2013-11-17 16:27:42
| Gaussian的TDDFT算激发态频率是数值的，所以需要算6N+1次（N是体系的原子数目）TD，所以你会找到6N+1个"SCF Done"，什么时候"SCF Done"的数目到达这个了，什么时候你的激发态频率计算就结束了。2天时间太少，TDDFT算个激发态频率20天到一个月都是正常的，慢慢等吧。没事数数"SCF Done"的数目，大概估算一下需要算多久吧。 |
5楼2013-11-14 13:53:10
shenzhendsh
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3.3
- 散金: 28
- 帖子: 61
- 在线: 68.3小时
- 虫号: 1843099
- 注册: 2012-06-01
- 性别: GG
- 专业: 高分子物理与高分子物理化
6楼2013-11-17 16:27:31
blueybz
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3023.5
- 散金: 359
- 红花: 1
- 帖子: 469
- 在线: 247.7小时
- 虫号: 497342
- 注册: 2008-02-02
- 性别: GG
- 专业: 理论和计算化学
7楼2013-11-18 11:19:46
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
|
大神你好，我在这问您个问题，就是我有个分子，大约含有110额原子，之前优化激发态一直出错，使用的方法是b3lyp/6-31g*。假如说，我计算激发态的频率，能使用b3lyp/6-31g*吗？sob给的解说是这样的：Gaussian没有TDDFT解析二阶导数，这是几乎做不动的，耗时我估计是600~1200次激发能计算的耗时。必须要计算频率的话，只能是在CIS下优化并计算频率，因为CIS下有解析二阶导数，速度和TDDFT计算频率耗时完全不在一个数量级。 我想，基态是b3lyp/6-31g*,激发态用CIS能行吗？b3lyp真的做不了我的这个计算？我曾经计算过一个70个分子，很长很长时间，我该怎么办？谢谢pwzhou的指点。谢谢您 |
8楼2015-01-14 23:41:04
|
你好，请问激发态的频率计算怎么做？为什么我计算老出错 分子的激发态优化后做频率计算 输入文件是 %chk=b.chk %rwf=b.rwf %mem=200MW %nprocshare=12 # freq td=(singlets,nstates=6,root=1) b3lyp/6-31g(d) 计算到下边就停止了 SavETr: write IOETrn= 770 NScale= 10 NData= 16 NLR=1 LETran= 118. IDoAtm=1111111111111111111111111111111111111111 Differentiating once with respect to electric field. with respect to dipole field. Differentiating once with respect to nuclear coordinates. There are 1 degrees of freedom in the 1st order CPHF. IDoFFX=0. LinEq1: Iter= 0 NonCon= 1 RMS=2.02D-03 Max=2.21D-01 AX will form 1 AO Fock derivatives at one time. LinEq1: Iter= 1 NonCon= 1 RMS=4.31D-04 Max=2.25D-02 LinEq1: Iter= 2 NonCon= 1 RMS=2.91D-04 Max=1.36D-02 LinEq1: Iter= 3 NonCon= 1 RMS=2.27D-04 Max=1.60D-02 LinEq1: Iter= 4 NonCon= 1 RMS=1.99D-04 Max=1.26D-02 LinEq1: Iter= 5 NonCon= 1 RMS=1.04D-04 Max=6.16D-03 LinEq1: Iter= 6 NonCon= 1 RMS=5.06D-05 Max=1.92D-03 LinEq1: Iter= 7 NonCon= 1 RMS=2.28D-05 Max=1.17D-03 LinEq1: Iter= 8 NonCon= 1 RMS=1.49D-05 Max=1.22D-03 LinEq1: Iter= 9 NonCon= 1 RMS=1.48D-05 Max=1.31D-03 LinEq1: Iter= 10 NonCon= 1 RMS=7.12D-06 Max=2.68D-04 LinEq1: Iter= 11 NonCon= 1 RMS=3.05D-06 Max=1.36D-04 LinEq1: Iter= 12 NonCon= 1 RMS=2.41D-06 Max=1.46D-04 LinEq1: Iter= 13 NonCon= 1 RMS=1.28D-06 Max=5.58D-05 |
9楼2016-01-19 22:48:51
5