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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】含时密度泛函计算金团簇的激发态输出出错
51/1返回列表
查看: 2998  |  回复: 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

铁杆木虫 (正式写手)
我是看文献上采用的这个组合的：
“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的回帖
查看全部 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的回帖

qzhaosdu

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lanl2dz 不能直接写，要gen，pesudo=read...
赞一下(1人) 回复此楼
Anewday,anewhour,anewminute,anewpeople.
4楼2010-11-28 18:15:27
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工学硕294分，求导师收留 +14 yzyzx 2026-04-12 14/700 2026-04-13 08:46 by lhj2009
[考研] 考研求调剂 +6 子木呐 2026-04-12 6/300 2026-04-13 07:54 by wp06
[考研] 求调剂，262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 316求调剂 +5 想读研究生( ?∵ 2026-04-07 5/250 2026-04-12 00:43 by 蓝云思雨
[考研] 电子信息279求调剂，有书读就行 +8 wwwooden 2026-04-08 11/550 2026-04-11 20:22 by cq2548
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 283求调剂，工科！ +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 291 求调剂 +29 化工2026届毕业� 2026-04-09 29/1450 2026-04-10 22:55 by dick_runner
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 生物学308分求调剂（一志愿华东师大）做过分子实验 +8 相信必会光芒万� 2026-04-07 9/450 2026-04-10 21:03 by zhouxiaoyu
[考研] 287求调剂 +15 Fnhc 2026-04-07 21/1050 2026-04-10 19:09 by chemisry
[考研] 生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 材料科学与工程320求调剂，080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭