mchen10

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 13 (小學(xué)生)
金幣: 1508.2
帖子: 425
在線(xiàn): 159.5小時(shí)
蟲(chóng)號(hào): 695973

[交流] 【求助】MOLPRO SO-ECP 是否需要加權(quán)

一個(gè)關(guān)于用MOLPRO計(jì)算 spin-obit coupling的問(wèn)題，請(qǐng)高手指點(diǎn)，在下感激不盡

molpro manual (40.7) 上面有一個(gè) 碘原子 SO計(jì)算的例子

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 項(xiàng)都進(jìn)行了加權(quán)

然而Xe SO ecp的作者之一Hermann Stoll在[molpro-user]list里面給出的input中scalar部分卻沒(méi)有加權(quán)

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

測(cè)試了一下，如果：
-只對(duì)SO加權(quán)，對(duì)hf 和 mcscf 沒(méi)有影響
-對(duì)SO和scalar加權(quán)，hf 和 mcscf有1~幾個(gè)a.u.的變化，spin-orbit correction 有幾千cm-1的改變

請(qǐng)問(wèn)我該不該對(duì)scalar ecp進(jìn)行加權(quán)？

回復(fù)此樓

» 猜你喜歡

大灣區(qū)大學(xué)劉天輝課題組招收2026級(jí)博士生（與中山大學(xué)物理學(xué)院聯(lián)合培養(yǎng)）已經(jīng)有22人回復(fù)
PbS量子點(diǎn)紫外吸收已經(jīng)有0人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有241人回復(fù)
比利時(shí)魯汶大學(xué)與國(guó)家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship 已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)四氫呋喃溶解的聚合物用甲醇沉淀時(shí)，如何使沉淀過(guò)程加速？已經(jīng)有2人回復(fù)
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生已經(jīng)有11人回復(fù)
Chemical Bonding at Surfaces and Interfaces,最經(jīng)典的一本表面上化學(xué)相互作用教材已經(jīng)有0人回復(fù)
中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復(fù)
誠(chéng)邀有志之士加入江蘇大學(xué)環(huán)境學(xué)院合成生物學(xué)課題組！(長(zhǎng)期接收學(xué)碩專(zhuān)碩、碩博連讀) 已經(jīng)有16人回復(fù)
一志愿南大理論與計(jì)算化學(xué)方向（070300）342分，求調(diào)劑已經(jīng)有1人回復(fù)
歡迎化學(xué)、化工、物理、材料、或計(jì)算機(jī)背景的考研咨詢(xún) 已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

是否硅烷偶聯(lián)劑作為天然高分子的封端劑？已經(jīng)有11人回復(fù)
液相有關(guān)物質(zhì)分析方法驗(yàn)證中用到的雜質(zhì)需要標(biāo)定嗎？已經(jīng)有14人回復(fù)
[VASP編譯]一什么是.so文件、.a文件？已經(jīng)有11人回復(fù)
So tired these days 已經(jīng)有18人回復(fù)
請(qǐng)問(wèn)大家下面這種情況需要出版社的permission嗎？已經(jīng)有7人回復(fù)
怎么判斷提取出的蛋白質(zhì)是否含有需要的蛋白已經(jīng)有5人回復(fù)
訪(fǎng)問(wèn)學(xué)者孩子在上幼兒園是否需要提供國(guó)內(nèi)配偶信息？已經(jīng)有3人回復(fù)
誰(shuí)懂迭代加權(quán)最小二乘法，能否給講下原理已經(jīng)有6人回復(fù)
這種情況，是否需要回信？已經(jīng)有6人回復(fù)
PNIPAm制備過(guò)程中是否需要氮?dú)獗Ｗo(hù)？已經(jīng)有6人回復(fù)
SOFC型家用燃料電池銷(xiāo)量2012年3月底有望達(dá)到近900塊已經(jīng)有5人回復(fù)
請(qǐng)教用molpro打印氫分子離子的重疊積矩陣、單電子積分矩陣、偶極矩陣已經(jīng)有3人回復(fù)
關(guān)于SOFC電解質(zhì)材料的一個(gè)問(wèn)題已經(jīng)有3人回復(fù)
【求助】關(guān)于除SOCl2 已經(jīng)有19人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-04-03 13:04:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

beefly

專(zhuān)家顧問(wèn) (職業(yè)作家)

★
mchen10(金幣+1):謝謝參與
mchen10(金幣+5): 非常感謝 2011-04-04 10:32:39
mchen10(金幣+5): molpro的標(biāo)量部分和矢量部分都已做過(guò)換算，再次感謝 2011-05-26 05:39:08

除了ci模塊中的旋軌耦合計(jì)算外，molpro所有的計(jì)算都是標(biāo)量的，因此修改so部分的權(quán)重不影響標(biāo)量計(jì)算的結(jié)果。只有做旋軌耦合計(jì)算才能看出影響。

權(quán)重因子沒(méi)有必要直接寫(xiě)出來(lái)，所以第二個(gè)例子的標(biāo)量ECP部分直接輸入了換算后的值。SO-ECP部分的表達(dá)式在文獻(xiàn)中有不同的定義，因此有可能需要乘上換算因子（見(jiàn)最后兩段：http://www.nwchem-sw.org/index.php/ECP）。molpro自帶的so-ecp已經(jīng)作了換算，不用再乘上因子了。

贊一下(1人)

回復(fù)此樓

2樓2011-04-03 23:36:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mchen10

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 13 (小學(xué)生)
金幣: 1508.2
帖子: 425
在線(xiàn): 159.5小時(shí)
蟲(chóng)號(hào): 695973

引用回帖:

Originally posted by beefly at 2011-04-03 09:36:30:
除了ci模塊中的旋軌耦合計(jì)算外，molpro所有的計(jì)算都是標(biāo)量的，因此修改so部分的權(quán)重不影響標(biāo)量計(jì)算的結(jié)果。只有做旋軌耦合計(jì)算才能看出影響。

權(quán)重因子沒(méi)有必要直接寫(xiě)出來(lái)，所以第二個(gè)例子的標(biāo)量ECP部分直接輸 ...

謝謝

關(guān)于第二個(gè)例子還有一點(diǎn)疑惑：例子中的標(biāo)量部分的系數(shù)是原paper中的系數(shù)（SO部分也一樣）。原文中，scalar-ecp公式中沒(méi)有加權(quán)系數(shù)，而SO-ecp部分系數(shù)B前面有加權(quán)系數(shù)2/（2l+1），是不是說(shuō)明scalar-ecp的系數(shù)B已經(jīng)吸收了其加權(quán)系數(shù)？

贊一下

回復(fù)此樓

3樓2011-04-04 10:40:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 333求調(diào)劑 +5	87639 2026-03-21	7/350	2026-03-21 19:31 by ColorlessPI
[考研] 求調(diào)劑 +5	十三加油 2026-03-21	5/250	2026-03-21 18:48 by 學(xué)員8dgXkO
[考研] 【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +11	吃吃吃才有意義 2026-03-19	11/550	2026-03-21 18:23 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 材料與化工專(zhuān)碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu