24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 598  |  回复: 2

mchen10

金虫 (正式写手)


[交流] 【求助】MOLPRO SO-ECP 是否需要 加权

一个关于用MOLPRO计算 spin-obit coupling的问题,请高手指点,在下感激不尽

molpro manual (40.7) 上面有一个 碘原子 SO计算的例子
CODE:
ecp,I,46,4,3;
1; 2,   1.00000000,       0.00000000;                                       ! lokal term = 0
2; 2,   3.50642001,      83.09814545; 2,   1.74736492,       5.06370919;    ! s-terme
4; 2,   2.99860773, 1/3* 81.88444526; 2,   3.01690894, 2/3* 83.41280402;    ! p-terms with weights
   2,   1.59415934, 1/3*  2.32392477; 2,   1.19802939, 2/3*  2.72079843;
4; 2,   1.03813792, 2/5*  6.40131754; 2,   1.01158599, 3/5*  6.21328827;    ! d-terms with weights
   2,   2.04193864, 2/5* 19.11604172; 2,   1.99631017, 3/5* 19.08465909;
4; 2,   2.64971585,-3/7* 24.79106489; 2,   2.75335574,-4/7* 24.98147319;    ! f-terms with weights
   2,   0.49970082,-3/7*  0.27936581; 2,   0.79638982,-4/7*  0.70184261;
4; 2,   2.99860773,-2/3* 81.88444526; 2,   3.01690894, 2/3* 83.41280402;    ! ECP-SO for p-terms
   2,   1.59415934,-2/3*  2.32392477; 2,   1.19802939, 2/3*  2.72079843;
4; 2,   1.03813792,-2/5*  6.40131754; 2,   1.01158599, 2/5*  6.21328827;    ! ECP-SO for d-terms
   2,   2.04193864,-2/5* 19.11604172; 2,   1.99631017, 2/5* 19.08465909;
4; 2,   2.64971585, 2/7* 24.79106489; 2,   2.75335574,-2/7* 24.98147319;    ! ECP-SO for f-terms
   2,   0.49970082, 2/7*  0.27936581; 2,   0.79638982,-2/7*  0.70184261;

其中scalar 和 SO 项都进行了加权

然而Xe SO ecp的作者之一Hermann Stoll在[molpro-user]list里面给出的input中scalar部分却没有加权
CODE:
ECP,XE,46,5,3;
1; ! LOCAL POTENTIAL
2, 1.00000000,.000000000;
2; ! S POTENTIAL
2, 3.94026300,122.763829;
2, 2.27726400,8.30885115;
2; ! P POTENTIAL
2, 3.02837300,68.8230044;
2, 1.39431900,3.64674223;
2; ! D POTENTIAL
2, 2.12260500,23.6520785;
2, 0.79866900,3.25844113;
2; ! F POTENTIAL
2, 6.16436000,-47.703199;
2, 1.54237400,-6.5411399;
1; ! G POTENTIAL
2, 1.84789200,-7.1058506;
2; ! P SO POTENTIAL
2, 3.02837300,-1.7460911*2/3;
2, 1.39431900,2.19487257*2/3;
2; ! D SO POTENTIAL
2, 2.12260500,-1.0456759*2/5;
2, 0.79866900,0.27199291*2/5;
1; ! F SO POTENTIAL
2, 6.16436000,0.20772207*2/7;

Hermann Stoll
University of Stuttgart
Germany

测试了一下,如果:
-只对SO加权,对hf 和 mcscf 没有影响
-对SO和scalar加权,hf 和 mcscf有1~几个a.u.的变化,spin-orbit correction 有几千cm-1的改变

请问我该不该对scalar ecp进行加权?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

mchen10(金币+1):谢谢参与
mchen10(金币+5): 非常感谢 2011-04-04 10:32:39
mchen10(金币+5): molpro的标量部分和矢量部分都已做过换算,再次感谢 2011-05-26 05:39:08
除了ci模块中的旋轨耦合计算外,molpro所有的计算都是标量的,因此修改so部分的权重不影响标量计算的结果。只有做旋轨耦合计算才能看出影响。

权重因子没有必要直接写出来,所以第二个例子的标量ECP部分直接输入了换算后的值。SO-ECP部分的表达式在文献中有不同的定义,因此有可能需要乘上换算因子(见最后两段:http://www.nwchem-sw.org/index.php/ECP)。molpro自带的so-ecp已经作了换算,不用再乘上因子了。
2楼2011-04-03 23:36:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)


引用回帖:
Originally posted by beefly at 2011-04-03 09:36:30:
除了ci模块中的旋轨耦合计算外,molpro所有的计算都是标量的,因此修改so部分的权重不影响标量计算的结果。只有做旋轨耦合计算才能看出影响。

权重因子没有必要直接写出来,所以第二个例子的标量ECP部分直接输 ...

谢谢

关于第二个例子还有一点疑惑:例子中的标量部分的系数是原paper中的系数(SO部分也一样)。原文中,scalar-ecp公式中没有加权系数,而SO-ecp部分系数B前面有加权系数2/(2l+1),是不是说明scalar-ecp的系数B已经吸收了其加权系数?
3楼2011-04-04 10:40:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mchen10 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 352分 化工与材料 +4 海纳百川Ly 2026-03-27 4/200 2026-03-28 01:08 by 飞行日记西
[考研] 一志愿211院校 344分 东北农业大学生物学学硕,求调剂 +5 丶风雪夜归人丶 2026-03-26 8/400 2026-03-27 19:22 by 丶风雪夜归人丶
[考研] 307求调剂 +8 超级伊昂大王 2026-03-24 9/450 2026-03-27 15:34 by 超级伊昂大王
[考研] 274求调剂 +17 顾九笙要谦虚 2026-03-24 23/1150 2026-03-27 15:16 by caszguilin
[考研] 279 分 求调剂 +4 睡个好觉_16 2026-03-24 4/200 2026-03-27 15:05 by 醉在风里
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[考研] 317求调剂 +7 蛋黄咸肉粽 2026-03-26 7/350 2026-03-27 02:29 by fmesaito
[考研] 349求调剂 +5 杰斯塔里斯 2026-03-21 5/250 2026-03-27 00:31 by wxiongid
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 324求调剂 +4 wysyjs25 2026-03-21 4/200 2026-03-26 20:38 by fmesaito
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 生物学学硕,一志愿湖南大学,初试成绩338 +4 YYYYYNNNNN 2026-03-26 4/200 2026-03-26 19:00 by macy2011
[考研] 085601求调剂总分293英一数二 +4 钢铁大炮 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
[考研] 081700 调剂 267分 +11 迷人的哈哈 2026-03-23 11/550 2026-03-26 15:41 by zzll406
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4 小鱼爱有机 2026-03-25 4/200 2026-03-26 10:19 by aa331100
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 289材料与化工(085600)B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
信息提示
请填处理意见