| 查看: 1173 | 回复: 4 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
javy金虫 (正式写手)
|
[交流]
【求助】激发态振子强度的问题
|
|||
|
在激发态计算中 （分子中含Ｉ Ｐ Ｎ Ｈ Ｃ Ｃ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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 光谱 |
» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有99人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有2人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
招收博士一名-宁波大学
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
javy
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1442.5
- 红花: 1
- 帖子: 340
- 在线: 123.2小时
- 虫号: 419009
- 注册: 2007-07-04
- 性别: MM
- 专业: 理论和计算化学
5楼2009-02-18 15:26:01
2楼2009-02-17 20:45:06
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13319.5
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★
jghe(金币+2,VIP+0):感谢你参与讨论！ 2-18 20:32
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 ] |
3楼2009-02-18 00:29:49
javy
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1442.5
- 红花: 1
- 帖子: 340
- 在线: 123.2小时
- 虫号: 419009
- 注册: 2007-07-04
- 性别: MM
- 专业: 理论和计算化学
4楼2009-02-18 10:56:57
