24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1008  |  回复: 4
当前主题已经存档。

javy

金虫 (正式写手)

[交流] 【求助】激发态振子强度的问题

在激发态计算中 (分子中含I P N H C Cu等原子)非金属用全电子基组,金属用赝势。计算的结果中,出现振子强度为负值,如下:

Excited State   1:    654.19 nm  f=-4.9722  
Excited State   2:    625.65 nm  f=1.5388   
Excited State   3:    605.57 nm  f=-5.5431  
Excited State   4:    577.40 nm  f=0.0000   
Excited State   5:    535.99 nm  f=0.0000   
Excited State   6:    519.18 nm  f=0.0000   
Excited State   7:    426.77 nm  f=0.0000   
Excited State   8:    396.77 nm  f=-0.0000  
Excited State   9:    395.85 nm  f=0.0000   
Excited State  10:    386.23 nm  f=-4.7455  
Excited State  11:    382.79 nm  f=1.4145   
Excited State  12:    380.43 nm  f=-5.4315  
Excited State  13:    371.83 nm  f=-0.0000  
Excited State  14:    364.33 nm  f=0.0000   
Excited State  15:    363.70 nm  f=-0.0000  
Excited State  16:    362.45 nm  f=0.0000   
Excited State  17:    350.43 nm  f=0.0000   
Excited State  18:    348.03 nm  f=-0.0000  
Excited State  19:    344.12 nm  f=-6.0471  
Excited State  20:    333.77 nm  f=-0.6528  
Excited State  21:    332.42 nm  f=-6.0805  
Excited State  22:    327.30 nm  f=-0.0000  
Excited State  23:    325.74 nm  f=-0.0000  
Excited State  24:    324.54 nm  f=-0.0000  
Excited State  25:    319.39 nm  f=0.0000   
Excited State  26:    318.68 nm  f=-0.0000  
Excited State  27:    316.25 nm  f=-0.0000  
Excited State  28:    308.97 nm  f=-5.8379  
Excited State  29:    305.80 nm  f=-0.5624  
Excited State  30:    305.51 nm  f=-6.0389  
Excited State  31:    304.97 nm  f=-0.0000  
Excited State  32:    304.44 nm  f=-0.0000  
Excited State  33:    300.88 nm  f=-0.0000  
Excited State  34:    298.04 nm  f=0.0000   
Excited State  35:    297.40 nm  f=0.0000   
Excited State  36:    295.14 nm  f=-0.0000  
Excited State  37:    293.29 nm  f=-5.1809  
Excited State  38:    289.74 nm  f=1.8682   
Excited State  39:    287.79 nm  f=-4.0837  
Excited State  40:    286.83 nm  f=0.0000   
Excited State  41:    286.39 nm  f=0.0000   
Excited State  42:    281.96 nm  f=-0.0000  
Excited State  43:    281.20 nm  f=-0.0000  
Excited State  44:    278.86 nm  f=-0.0000  
Excited State  45:    275.54 nm  f=0.0000   
Excited State  46:    272.79 nm  f=-5.2225  
Excited State  47:    272.13 nm  f=-1.1349  
Excited State  48:    270.90 nm  f=-4.1924  
Excited State  49:    267.11 nm  f=0.0000   
Excited State  50:    265.32 nm  f=-0.0000  
  
以前没碰到过这种情况。 想问下为什么这里的 f 有的会是负值

谢谢大家

[ Last edited by erylingjet on 2009-2-17 at 21:57 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)


xuefei06(金币+1,VIP+0):thanks! 2-17 22:58
怎么会有这么多负值阿,你们的程序好像是吴东升改过的,你用原版的程序重新算算
2楼2009-02-17 20:45:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
jghe(金币+2,VIP+0):感谢你参与讨论! 2-18 20:32
振子强度(f)的计算公式是:(2/3)*跃迁偶极矩(TDM)的平方*垂直激发能(dE)

全部是au单位

如果激发能为负值(也即,参考态不是基态),f就是负的。但你列出的dE都是正值。

还一个可能是TDM为虚数,导致负的f,但这种情况只有在奇数电子体系的旋轨耦合计算中出现(这时候上面公式中的TDM**2要改为|TDM|**2,所以f的符号只由dE决定),绝大多数程序的TDDFT不能考虑旋轨耦合,因此不应该出现虚数的TDM,所以肯定是程序源代码或编译器的问题。

建议先检查一下程序打印的TDM,然后根据上面的公式,重算一下f。如果担心结果不可靠,还可以用官方版本的gaussian或其他程序(如gamess)验证一下。

[ Last edited by beefly on 2009-2-18 at 00:42 ]
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2009-02-18 00:29:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

javy

金虫 (正式写手)

谢谢楼上的两位
看了一下TDM,都是正常的。
只能重新算一下看看
4楼2009-02-18 10:56:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

javy

金虫 (正式写手)

好象说这个问题是某些版本gaussian的bug。
问题解决了,谢谢
5楼2009-02-18 15:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 javy 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
[考研] 310求调剂 +17 争取九点睡 2026-03-30 17/850 2026-04-02 16:40 by guanxin1001
[考研] 学硕化学工程与技术,一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 301求调剂 +14 骆驼男人 2026-04-02 14/700 2026-04-02 14:08 by baoball
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 266分求材料化工冶金矿业等专业的调剂 +7 哇呼哼呼哼 2026-03-26 7/350 2026-04-01 20:47 by cqupH
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 324求调剂 +10 hanamiko 2026-03-26 12/600 2026-04-01 16:59 by oooqiao
[考研] 086502化学工程342求调剂 +7 阿姨复古不过 2026-03-27 7/350 2026-04-01 16:14 by yanflower7133
[考研] 070300求调剂306分 +5 26要上岸 2026-03-27 5/250 2026-04-01 11:09 by oooqiao
[考研] 288资源与环境专硕求调剂,不限专业,有学上就行 +25 lllllos 2026-03-30 26/1300 2026-04-01 09:52 by 一只好果子?
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 289求调剂 +3 Acesczlo 2026-03-29 4/200 2026-03-31 14:48 by 热情沙漠
[考研] 本科双非材料,跨考一志愿华电085801电气,283求调剂,任何专业都可以 +6 芝士雪baoo 2026-03-28 8/400 2026-03-29 08:16 by 松花缸1201
[考研] 数一英一271专硕(085401)求调剂,可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
信息提示
请填处理意见