|
|
各位老师好，学生想用mrci方法计算PN分子激发态的势能曲线。用casscf方法生成初始波函数，计算中根据别人的工作将PN分子的9 个外层分子轨道（5a1, 2b1和2b2）设为活性空间，即5-9σ, 2π和3π分子轨道, PN分子的10 个价电子放入活化空间中; 剩下的 12 个内层电子放入6个封闭的分子轨道(4a1, 1b1和1b2)，即1-4σ和1π分子轨道。现在想用mrci计算波函数B2对称性下的激发态，我在该对称性下计算了4个态，并把最终结果用Table指令统计了出来，发现B2对称性下第四个态的势能曲线的趋势有问题，体系能量先降低后升高，然后体系能量在5.8 bohr核间距的位置又降低了，请问各位老师这是哪里出了问题？非常感谢!(输入和输出文件已上传,势能曲线的out文件过大无法上传，所以删除了部分循环计算的内容以便上传)
R E_TOL1 E_TOL2 E_TOL3 E_TOL4
2.0 -394.7960693 -394.6202303 -394.5614218 -394.5413943
2.1 -394.9758221 -394.8063949 -394.7647298 -394.7456134
2.2 -395.1095300 -394.9451382 -394.9161222 -394.8985210
2.3 -395.2075182 -395.0486008 -395.0265980 -395.0117334
2.4 -395.2779509 -395.1253264 -395.1057249 -395.0940079
2.5 -395.3274754 -395.1807096 -395.1624178 -395.1517617
2.6 -395.3606295 -395.2199678 -395.2018742 -395.1910812
2.7 -395.3816114 -395.2460983 -395.2305942 -395.2130661
2.8 -395.3933570 -395.2632614 -395.2484316 -395.2291928
2.9 -395.3984466 -395.2741439 -395.2602142 -395.2353865
3.0 -395.3983175 -395.2810775 -395.2658198 -395.2380777
3.1 -395.3948624 -395.2855828 -395.2639059 -395.2386576
3.2 -395.3887073 -395.2881241 -395.2623046 -395.2450716
3.3 -395.3808484 -395.2887381 -395.2593242 -395.2486056
3.4 -395.3718856 -395.2877354 -395.2561053 -395.2499641
3.5 -395.3622708 -395.2854326 -395.2532928 -395.2496996
3.6 -395.3523425 -395.2821443 -395.2508060 -395.2482598
3.8 -395.3324580 -395.2735749 -395.2430852 -395.2451734
4.0 -395.3135101 -395.2637214 -395.2429916 -395.2366985
4.2 -395.2958954 -395.2541979 -395.2381284 -395.2312316
4.4 -395.2800679 -395.2454079 -395.2315513 -395.2259069
4.6 -395.2660081 -395.2380071 -395.2244469 -395.2217313
4.8 -395.2537853 -395.2319202 -395.2175049 -395.2182999
5.0 -395.2434907 -395.2269204 -395.2165183 -395.2100083
5.2 -395.2351804 -395.2229249 -395.2149985 -395.2033530
5.4 -395.2287342 -395.2199219 -395.2139671 -395.1974318
5.6 -395.2239465 -395.2177905 -395.2133475 -395.1929824
5.8 -395.2214066 -395.2170292 -395.2132797 -395.2102262
6.0 -395.2188173 -395.2163994 -395.2131985 -395.2111741
6.2 -395.2170299 -395.2158801 -395.2132330 -395.2118790
6.4 -395.2158626 -395.2153716 -395.2133148 -395.2124141
6.6 -395.2151212 -395.2148750 -395.2133948 -395.2128230
6.8 -395.2146376 -395.2144435 -395.2134413 -395.2131315
7.0 -395.2143074 -395.2141346 -395.2134393 -395.2133374
7.2 -395.2140734 -395.2139584 -395.2134060 -395.2134320
7.4 -395.2138979 -395.2138679 -395.2133776 -395.2134381
7.6 -395.2137633 -395.2138150 -395.2133655 -395.2133941
7.8 -395.2136568 -395.2137765 -395.2133640 -395.2133296
8.0 -395.2135707 -395.2137457 -395.2133677 -395.2132631
|
|
发表于 Post on 2023-4-5 20:51:56
收藏 Add to favorites