应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】含时密度泛函计算金团簇的激发态输出出错
汕头大学海洋科学接受调剂
51/1返回列表
查看: 3013  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bjwang

铁杆木虫 (正式写手)

[交流] 【求助】含时密度泛函计算金团簇的激发态输出出错 已有5人参与

我计算的是两个金原子的二聚体，计划利用含时密度泛函算它激发态的能量和振荡强度（重复别人的工作练习gaussian），运行后，输出文件有错误信息，没有正常结束！运行在服务器上的，g09-linux版。
         哪位同学帮忙看看，可能比较低级，还望指点一二呀！谢谢！
  
   附 输入文件：

%chk=au2.chk
%rwf=au2.rwf
%mem=2GB
#p pw91/LANL2DZ  td(Triplets,nstates=5)

The TD calc of au2

0  1
Au 0.00 0.00 0.00
Au 2.55 0.00 0.00



附 输出文件如下：

Gaussian 09:  EM64L-G09RevA.01  8-May-2009
                28-Nov-2010
******************************************
%chk=au2.chk
%rwf=au2.rwf
%mem=2GB
--------------------------------------
#p pw91/LANL2DZ td(Triplets,nstates=5)
--------------------------------------
QPErr --- An ambiguous keyword was detected.
#p pw91/LANL2DZ td(Triplets,nstates=5)
    '
Last state="DFT"
TCursr=19940 LCursr=    3
Error termination via Lnk1e in /home/bjwang/g09/l1.exe at Sun Nov 28 15:55:44 2010.
Job cpu time:  0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1
回复此楼

» 收录本帖的淘帖专辑推荐

不懂计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
1楼 2010-11-28 16:43:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)
楼主说的有些麻烦呀，我参考例子test405.com（#p pw91PW91 6-31g(d,p) opt=tight freq test）和例子test588.com（#p pbepbe/lanl2dz opt=(modred,expert) test
），把上面的改为：
#p pw91pw91/LANL2DZ  td(Singlets,nstates=25)
可以正常运行了，但结果和文献的好像不太一致。

我把一部分结果贴出来，还望大家帮忙指点一二,谢谢！

附计算结果：
Excitation energies and oscillator strengths:

Excited State   1:      Singlet-SGG    2.4542 eV  505.20 nm  f=0.0000  =0.000
      19 -> 20         0.70804
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -270.944678611
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited state symmetry could not be determined.
Excited State   2:      Singlet-?Sym    2.4692 eV  502.12 nm  f=0.0068  =0.000
      16 -> 20         0.54652
      17 -> 20         0.44957

Excited state symmetry could not be determined.
Excited State   3:      Singlet-?Sym    2.4692 eV  502.12 nm  f=0.0068  =0.000
      16 -> 20        -0.44957
      17 -> 20         0.54652

Excited state symmetry could not be determined.
Excited State   4:      Singlet-?Sym    2.5208 eV  491.85 nm  f=0.0000  =0.000
      15 -> 20         0.70723

Excited state symmetry could not be determined.
Excited State   5:      Singlet-?Sym    2.5208 eV  491.84 nm  f=0.0000  =0.000
      14 -> 20         0.70723

Excited State   6:      Singlet-SGU    2.8115 eV  440.99 nm  f=0.1547  =0.000
       9 -> 20        -0.13683
      18 -> 20         0.69666
Excited state symmetry could not be determined.
Excited State   7:      Singlet-?Sym    2.8126 eV  440.81 nm  f=0.0000  =0.000
      13 -> 20         0.70718

Excited state symmetry could not be determined.
Excited State   8:      Singlet-?Sym    2.8127 eV  440.80 nm  f=0.0000  =0.000
      12 -> 20         0.70718

Excited state symmetry could not be determined.
Excited State   9:      Singlet-?Sym    3.8587 eV  321.31 nm  f=0.0000  =0.000
      11 -> 20         0.70619

Excited state symmetry could not be determined.
Excited State  10:      Singlet-?Sym    3.8587 eV  321.31 nm  f=0.0000  =0.000
      10 -> 20         0.70619

Excited State  11:      Singlet-SGU    5.2842 eV  234.63 nm  f=0.0181  =0.000
       9 -> 20         0.47641
      16 -> 21        -0.36704
      17 -> 22        -0.36704

Excited state symmetry could not be determined.
Excited State  12:      Singlet-?Sym    5.4335 eV  228.19 nm  f=0.0000  =0.000
      19 -> 22         0.68975
symmetry could not be determined.
Excited State  13:      Singlet-?Sym    5.4335 eV  228.19 nm  f=0.0000  =0.000
      19 -> 21         0.68975

Excited state symmetry could not be determined.
Excited State  14:      Singlet-?Sym    5.4625 eV  226.98 nm  f=0.0000  =0.000
      16 -> 21         0.49989
      17 -> 22        -0.49989

Excited state symmetry could not be determined.
Excited State  15:      Singlet-?Sym    5.4625 eV  226.97 nm  f=0.0000  =0.000
      16 -> 22         0.49989
      17 -> 21         0.49989

Excited state symmetry could not be determined.
Excited State  16:      Singlet-?Sym    5.4661 eV  226.82 nm  f=0.0000  =0.000
      16 -> 22        -0.50000
      17 -> 21         0.50000

Excited state symmetry could not be determined.
Excited State  17:      Singlet-?Sym    5.6570 eV  219.17 nm  f=0.0000  =0.000
      14 -> 22        -0.49989
      15 -> 21         0.50005

Excited state symmetry could not be determined.
Excited State  18:      Singlet-?Sym    5.6570 eV  219.17 nm  f=0.0000  =0.000
      14 -> 21         0.49989
      15 -> 22         0.50005

Excited state symmetry could not be determined.
Excited State  19:      Singlet-?Sym    5.7361 eV  216.15 nm  f=0.1459  =0.000
      12 -> 22        -0.26016
      13 -> 21        -0.26020
      18 -> 21         0.53933
      18 -> 22        -0.20694

Excited state symmetry could not be determined.
Excited State  20:      Singlet-?Sym    5.7361 eV  216.15 nm  f=0.1459  =0.000
      12 -> 21        -0.26016

.................................................................

文献上给的五个最低的the energies of excited states 分别为：2.47,2.82，3.87，
5.28,5.74。振子强度分别为：0.01,0.15,0.00,0.02,0.15
      
       看结果里面第2激发态（激发能2.469，振子强度0.0068）；
                                         6（2.8115,0.1547）；
                                        9（3.8587,0.000）；
                                         11（5.2842,0.118）；
                                          19（5.7361,0.1459）。
这五个无论是能量还是强度和文献上的基本完全符合，但是文献上说是五个最低的the energies of excited states，这是里面跳出来的，还有更低的，算法有问题还是结果分析有问题呢？

    现在我还有点疑惑就是文献上的the energies of excited states和计算输出中的Excitation energies 是一样的吗？激发态的能量和激发能量？感觉激发态的能量应该等于基态的能量加上激发能量，还望各位兄弟给予指导呀！！

[ Last edited by bjwang on 2010-11-29 at 11:25 ]
赞一下 回复此楼
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
5楼2010-11-28 20:59:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
QPErr --- An ambiguous keyword was detected.
#p pw91/LANL2DZ td(Triplets,nstates=5)
错误在此，楼主查查手册有这个泛函组合吗？
赞一下(1人) 回复此楼
前途光明，出路难觅！
2楼2010-11-28 17:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)
我是看文献上采用的这个组合的：
“To assess the quality of the method, we calculate the five
lowest excitation energies of gold dimer at the level of PW91
with the LANL2DZ basis set”

     我在手册上确实没见到过这个组合，但我在网上搜了一下，也有好多人用呀
http://www.google.com.hk/search? ... p;q=pw91%2FLANL2DZ+

     pw91是DFT 方法的关键字，Gaussian内置的Lanl2dz基组是使用最广泛的ECP基组之一，广泛用于包含过渡金属体系的计算。

[ Last edited by bjwang on 2010-11-28 at 18:01 ]
赞一下 回复此楼
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
3楼2010-11-28 17:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lanl2dz 不能直接写，要gen，pesudo=read...
赞一下(1人) 回复此楼
Anewday,anewhour,anewminute,anewpeople.
4楼2010-11-28 18:15:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297工科调剂? +13 河南农业大学-能 2026-04-13 13/650 2026-04-14 16:46 by Art1977
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +25 慕绝cc 2026-04-09 29/1450 2026-04-14 16:45 by zhouxiaoyu
[考研] 一志愿085502，267分求调剂 +19 再忙也要吃饭啊 2026-04-08 20/1000 2026-04-14 16:03 by zs92450
[考研] 各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 求调剂 +3 我爱高数高数爱� 2026-04-12 3/150 2026-04-14 01:00 by 王珺璞
[考研] 求调剂，985材料与化工348分 +9 涵竹刘 2026-04-11 14/700 2026-04-13 22:26 by 涵竹刘
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +8 慕绝cc 2026-04-09 8/400 2026-04-13 14:08 by 张zhihao
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +5 弱水听文 2026-04-09 5/250 2026-04-12 13:16 by shengxi123
[考研] 0854调剂 +12 长弓傲 2026-04-09 13/650 2026-04-12 09:56 by 逆水乘风
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH（等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
信息提示
关闭
请填处理意见
关闭