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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】含时密度泛函计算金团簇的激发态输出出错
51/1返回列表
查看: 2919  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[考研] 284求调剂 +8 天下熯 2026-02-28 8/400 2026-03-02 00:15 by 暮雨星晴
[考研] 292求调剂 +6 yhk_819 2026-02-28 6/300 2026-03-01 23:23 by 向上的胖东
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[考研] 0856求调剂285 +10 吕仔龙 2026-02-28 10/500 2026-03-01 21:37 by 公瑾逍遥
[考研] 306分材料调剂 +4 chuanzhu川烛 2026-03-01 5/250 2026-03-01 19:48 by 无际的草原
[考研] 0856化工专硕求调剂 +12 董boxing 2026-03-01 12/600 2026-03-01 19:45 by 材子momo
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 290求调剂 +9 材料专硕调剂； 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭