24小时热门版块排行榜

返回列表

当前主题已经存档。

javy

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 1442.5
红花: 1
帖子: 340
在线: 123.2小时
虫号: 419009
注册: 2007-07-04
性别: MM
专业: 理论和计算化学

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

在激发态计算中　（分子中含Ｉ　Ｐ　Ｎ　Ｈ　Ｃ　Ｃu等原子）非金属用全电子基组，金属用赝势。计算的结果中，出现振子强度为负值，如下：

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 ]

回复此楼

» 收录本帖的淘帖专辑推荐

光谱

» 猜你喜欢

Ni元素XPS分析已经有1人回复
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有98人回复
求一份origin2019以上版本的origin软件压缩包已经有1人回复
东方理-中科大联合博士生招聘已经有0人回复
电化学基础知识与资源网站已经有3人回复
中国科大电池方向任晓迪课题组招收2026级博士生-电解液/电池安全性/人工智能方向已经有19人回复
26年秋季博士申请已经有0人回复
推荐给英语教学者的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有33人回复
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有24人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2009-02-17 16:53:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jghe

金虫 (著名写手)

应助: 10 (幼儿园)
贵宾: 0.225
金币: 112.2
散金: 5
红花: 3
帖子: 1481
在线: 123小时
虫号: 322338
注册: 2007-03-11
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 量子化学

★
xuefei06(金币+1,VIP+0):thanks！ 2-17 22:58

怎么会有这么多负值阿，你们的程序好像是吴东升改过的，你用原版的程序重新算算

赞一下(10人)

回复此楼

2楼2009-02-17 20:45:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +457
QC强帖: 12
应助: 86 (初中生)
金币: 12778.4
散金: 27307
红花: 300
帖子: 4664
在线: 925.5小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

★ ★
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 ]

赞一下(11人)

回复此楼

beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地