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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 【求助】MOLPRO SO-ECP 是否需要 加权
51/1返回列表
查看: 590  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

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的回帖
查看全部 3 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +3 pin8023 2026-02-28 5/250 2026-03-02 00:24 by 花YOU重开日
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 材料类求调剂 +10 wana_kiko 2026-02-28 12/600 2026-03-01 22:10 by 海嵙Y
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 6/300 2026-03-01 21:24 by cgyzqwn
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 化工专硕348，一志愿985求调剂 +5 弗格个 2026-02-28 8/400 2026-03-01 17:25 by sunny81
[考研] 290求调剂 +9 材料专硕调剂； 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 化工专硕342，一志愿大连理工大学，求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭