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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】含时密度泛函计算金团簇的激发态输出出错
51/1返回列表
查看: 2999  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程085601，270求调剂 +36 @ASDF1234 2026-04-08 39/1950 2026-04-13 08:37 by Delta2012
[考研] 复试调剂 +14 积极向上； 2026-04-10 16/800 2026-04-12 20:02 by gruyclewee
[考研] 2本，初试303，0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 086000调剂 +6 十七sa 2026-04-07 6/300 2026-04-12 11:05 by 大力水手力大无�
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 求调剂 +13 雪逢冬 2026-04-10 13/650 2026-04-11 09:58 by 猪会飞
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 一志愿华南理工大学331分材料求调剂 +9 天下ww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 生物学求调剂 一志愿沪9，326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
信息提示
关闭
请填处理意见
关闭