24小时热门版块排行榜    

查看: 585  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工299分求调剂 一志愿985落榜 +4 嘻嘻(*^ω^*) 2026-03-01 4/200 2026-03-01 13:15 by wang_dand
[考研] 0856材料专业298分有科研经历 硕士研究生调剂自荐信 +6 zyf上岸 2026-03-01 6/300 2026-03-01 12:43 by liqiongjy
[考研] 272求调剂 +5 材紫有化 2026-02-28 5/250 2026-03-01 11:51 by gaoxiaoniuma
[考研] 材料学调剂 +7 提神豆沙包 2026-02-28 9/450 2026-03-01 11:50 by 王伟要上岸啊
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +3 虫虫虫虫虫7 2026-03-01 6/300 2026-03-01 11:50 by gaoxiaoniuma
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 材料类求调剂 +8 wana_kiko 2026-02-28 8/400 2026-03-01 11:44 by 王伟要上岸啊
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
信息提示
请填处理意见