应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 【求助】MOLPRO SO-ECP 是否需要 加权
31/1返回列表
查看: 556  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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进行加权？
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-04-03 13:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)
★
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已经作了换算，不用再乘上因子了。
赞一下(1人) 回复此楼
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 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭