|
各位老师好,学生想用molpro计算PN分子的跃迁电偶极矩。学生参考文章中(参考文献已上传)的设置,计算中根据别人的工作将PN分子的9 个外层分子轨道(5a1, 2b1和2b2)设为活性空间,即5-9σ, 2π和3π分子轨道, PN分子的10 个价电子放入活化空间中; 剩下的 12 个内层电子放入6个封闭的分子轨道(4a1, 1b1和1b2),即1-4σ和1π分子轨道。先用casscf方法生成初始波函数,再用mrci方法计算PN从A1π态到X1Σ+态的跃迁电偶极矩(第一激发态到基态),但是我计算出来的结果和文献中的数据正负号完全相反,请问老师这是哪里出现了问题,非常感谢!(casscf方法和mrci方法初始输入和输出文件已上传,由于mrci输出文件过大无法上传,所以out里删除了前几次循环计算的内容以便上传)
这里是文献中的数据
A1π-X1Σ+
R(bohr) D_TRDMX R/Å TDM/a.u.
2.1 -0.20763137 1.05 0.2468816
2.2 -0.18588492 1.1 0.2166693
2.3 -0.16880757 1.15 0.1927694
2.4 -0.15513501 1.2 0.1735582
2.5 -0.14397937 1.25 0.1579853
2.6 -0.13469129 1.3 0.1454142
2.7 -0.12678991 1.32 0.1411181
2.8 -0.11991061 1.34 0.1372056
2.9 -0.11375339 1.36 0.1336237
3 -0.10809092 1.38 0.1303247
3.1 -0.10273228 1.4 0.1272568
3.2 -0.09746393 1.42 0.1242431
3.3 -0.07763021 1.44 0.1214592
3.4 -0.07241967 1.46 0.1188112
3.6 -0.06221643 1.48 0.1162642
3.8 -0.05243071 1.5 0.1138211
4 -0.04370465 1.52 0.1114757
4.2 -0.03694786 1.54 0.1092294
1.56 0.1071143
1.58 0.1054368
1.6 0.1043805
1.62 0.1022112
1.64 0.1000466
1.66 0.0979169
1.68 0.0958326
1.7 0.0937863
1.72 0.091766
1.74 0.0897859
1.76 0.0758602
1.78 0.0741823
1.8 0.0725279
1.82 0.0708977
1.84 0.0692932
1.86 0.0677155
1.88 0.0661676
1.9 0.0646519
1.95 0.0610076
2 0.0576458
|
|