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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】含时密度泛函计算金团簇的激发态输出出错
51/1返回列表
查看: 3520  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +12 koalala 2026-08-24 14/700 2026-08-24 12:41 by S09300097
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭