24小时热门版块排行榜    

查看: 1071  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ideallzc

铜虫 (小有名气)

[求助] Molpro 计算CO+ 不同π轨道 transition dipole时候出现跳变?

输入文件如下:
***COplus_pi_trans_dipole

memory,100,m;
symmetry,x,y
geometry={c;o,c,r(i)}

!basis=aug-cc-pV5Z
basis=AVQZ
!basis=cc-pVQZ

distances=[1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.65,2.7,2.75,2.8,2.9,3.0,3.2,3.4,3.6,3.8,3.9,4.0,4.2,4.6,5.0,5.5,6.0]

do i=1,#distances
r(i)=distances(i)

{hf;wf,13,1,1;}


{casscf; ! do casscf calculaion
orbital,ignore_error;
OCC,6,2,2,0; ! specify occupied space
FROZEN,2,0,0,0; ! specify Frozen core orbitals
WF,13,2,1;state,4;
} ! define wavefunction symmetry

{mrci;
orbital,ignore_error;
maxiter,100,100;
wf,13,2,1;state,4;}
pE1(i)=energy(1);
PE2(i)=energy(2);
PE3(i)=energy(3);
PE4(i)=energy(4);

enddo

根据*.out文件画出的跃迁矩阵随原子核间距的变化在附件中(以及1π,2π,3π,4π的势能曲线),我发现在某些原子核间距R的时候跃迁矩阵DMZ会有一些突变产生,例如2Pi~3Pi的DMZ在R=2.7a.u.附近会从正的变为负的,还有3Pi~4Pi的DMZ在R=2.2a.u.附近也会从负的变为正的,物理直觉上DMZ应该连续变化才对呀,感觉应该是计算上的问题,不知道怎么解决,求各位大神帮助。

Molpro 计算CO+ 不同π轨道 transition dipole时候出现跳变?
CO+_pppp_energy_surface.jpg


Molpro 计算CO+ 不同π轨道 transition dipole时候出现跳变?-1
CO+_pppp_trans_dipole.jpg
回复此楼
ideal-neil
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ideallzc

铜虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by virtualzx at 2016-01-14 02:04:51
透热态(diabatic state)的dipole才会平滑。绝热态的特性会不断变化,尤其在交叉点附近是跳跃的。跃迁偶极距符号本身也需要调整,尤其在交叉附近受到berry phase的影响,符号无法唯一确定。
但解离时即使看起来完全 ...

那这样的话,我能否用molpro求解透热态的势能曲线,以及transition dipole呢?还有一个问题,如果我考虑解离的话,分子是沿绝热势能曲线,还是沿透热曲线 解离?有什么依据么?  多谢!!!
ideal-neil
3楼2016-01-20 22:17:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

virtualzx

木虫 (著名写手)

透热态(diabatic state)的dipole才会平滑。绝热态的特性会不断变化,尤其在交叉点附近是跳跃的。跃迁偶极距符号本身也需要调整,尤其在交叉附近受到berry phase的影响,符号无法唯一确定。

但解离时即使看起来完全平滑的势能面,也只是假象,它的电子结构也是剧烈变化的,从一个闭壳层结构变成了两个单电子,所以偶极也会大幅变化。的到平滑偶极或跃迁

发自小木虫IOS客户端

» 本帖已获得的红花(最新10朵)

2楼2016-01-14 02:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见