| 查看: 3050 | 回復(fù): 12 | |||
科學(xué)島鐵蟲 (初入文壇)
|
[求助]
如何得到原子中某個(gè)軌道對(duì)分子總磁矩的貢獻(xiàn)? 已有1人參與
|
||
|
在Gaussian 中用pop=full,即可得到磁矩的數(shù)據(jù)。 Do spin unrestricted (or restricted open-shell) calculations (e.g.UB3LYP/6-311+G*) to get spin moments for each orbital printed. Be careful about spin-multiplicity though. Use pop=full or pop=regular in the route card for population analysis as well. Spin moments will be printed below Gross population analysis section. You can calculate total magnetic moment from that. 具體網(wǎng)址:http://www.ccl.net/chemistry/res ... .002-dir/index.html 但我想知道某個(gè)具體軌道對(duì)磁矩的貢獻(xiàn),知道怎么分析嗎?例如下圖中Au的5d軌道對(duì)總磁矩的貢獻(xiàn)怎么分析?謝謝! QQ截圖20131224102322.jpg |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
體系有n個(gè)自旋平行的未配對(duì)電子,分子電子自旋磁矩就為n*μ_B。而自旋相反的電子的自旋磁矩會(huì)抵消,所以閉殼層體系無電子自旋磁矩。通過布居分析(如Mulliken、NPA等),得到每個(gè)原子,或每個(gè)原子軌道上的自旋電子數(shù),就能討論它們各自對(duì)分子磁矩的貢獻(xiàn)。比如某d原子軌道上自旋電子數(shù)為0.4,就說這個(gè)d原子軌道上電子自旋磁矩是0.4μ_B。也可以根據(jù)電子自旋密度討論,某點(diǎn)自旋密度為x,就可以說這個(gè)點(diǎn)的電子對(duì)分子磁矩的貢獻(xiàn)是x*μ_B。也因此,可以將空間劃分為不同子區(qū)域,根據(jù)自旋密度在相應(yīng)區(qū)域積分值討論相應(yīng)區(qū)域內(nèi)的自旋磁矩。μ_B=e*h_bar/(2m_e)這個(gè)量被定義為玻爾磁子。 高斯布居分析只是給出每個(gè)基函數(shù)的布居數(shù),怎么和原子軌道對(duì)應(yīng)上需要基于你對(duì)所用基組的了解。pop=NPA給出的NPA分析結(jié)果則可以直接給出每個(gè)原子軌道的布居數(shù)。由于你用了贗勢(shì),所以NPA分析輸出的哪個(gè)是5d需要自行判斷。 另外,分子磁矩是電子自旋磁矩、電子軌道運(yùn)動(dòng)產(chǎn)生的磁矩和原子核自旋磁矩共同貢獻(xiàn)的,通過電子布居分析討論的只是電子自旋磁矩部分,一般總磁矩也是由此主導(dǎo)。 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
鐵蟲 (初入文壇)
|
謝謝! 一看就知道您對(duì)這方面很了解。我剛接觸這方面的內(nèi)容,按我對(duì)您回復(fù)的理解,以下幾點(diǎn)是否正確: 1.NPA得到各原子的各電子層上的電子布居數(shù)n,相應(yīng)的磁矩就是n*μ_B了?某原子的磁矩就是各軌道磁矩之和,各原子磁矩之和就是團(tuán)簇的總磁矩? 2.http://www.ccl.net/chemistry/res ... .002-dir/index.html 中說是根據(jù)Gross orbital populations之后的數(shù)據(jù)來計(jì)算磁矩,現(xiàn)在主流的是采用“NPA得到各原子的各電子層上的電子布居數(shù)n,相應(yīng)的磁矩就是n*μ_B”還是"Gross orbital populations之后的數(shù)據(jù)來計(jì)算磁矩"思路。 3.網(wǎng)上查詢了一下,也有用MS計(jì)算磁矩的, 4.您自己有沒有這方面的已經(jīng)發(fā)表的文章,到時(shí)可以引用引用。 5.另外,您有沒有聽說過自旋污染? |
鐵蟲 (初入文壇)
|
1.附件中分別是NBO計(jì)算的輸入和輸出文件,是不是總的“ NATURAL POPULATIONS: Natural atomic orbital occupancies” 減去Alpha spin orbitals 的“ NATURAL POPULATIONS: Natural atomic orbital occupancies”就是相應(yīng)的beta 自旋軌道的,之后Alpha,beta spin orbitals 的電子布局?jǐn)?shù)之差就是相應(yīng)原子軌道的未配對(duì)電子? 2.另外,總磁矩不需要考慮激發(fā)態(tài)吧, 例如下面的就不需要考慮了吧, 3 Au 1 S Ryd( 7S) 0.00182 2.01677 4 Au 1 S Ryd( 8S) 0.00003 8.33001 5 Au 1 S Ryd( 9S) 0.00000 70.15820 2.若不對(duì),能否在輸出文件中指出,萬分謝謝! |
鐵蟲 (初入文壇)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
1 不是這么簡(jiǎn)單相乘。比如某軌道上alpha電子布居數(shù)是na,beta是nb,那么自旋布居就是na-nb,(na-nb)*μ_B才是這個(gè)軌道的電子貢獻(xiàn)的磁矩。后面兩句話你的理解正確。 2 高斯輸出的Gross orbital populations是基于Mulliken方法得到的基函數(shù)的布居數(shù),這種計(jì)算方式很不推薦使用。尤其是帶了彌散函數(shù)的情況更不能用(比如CCL上那個(gè)人拿6-311+G*舉例就很不合適),NPA得到的軌道布居數(shù)有意義得多。Mulliken方法的缺陷在《談?wù)勡壍莱煞莸挠?jì)算方法》(http://hi.baidu.com/sobereva/item/d52876160110ca8d89a956c5)、《原子電荷計(jì)算方法的對(duì)比》(http://www.whxb.pku.edu.cn/EN/10.3866/PKU.WHXB2012281)上都有討論 3 問題不在于程序,而在于原理 4 我的文章里沒有專門討論磁矩的 5 自旋污染幾乎是眾所周知的問題 |
鐵蟲 (初入文壇)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
鐵蟲 (初入文壇)
|
您好!我想追問一個(gè)磁矩的問題;我想得到AuC中性團(tuán)簇中Au的5d,6s,6p軌道對(duì)總磁矩的貢獻(xiàn),按照您的思路,alpha軌道電子布居為0.6481,beta軌道電子布局為0.6759,總自旋布居為0.6481-0.6759=-0.0278.總布局為負(fù)值,您如何理解?附件中為輸出文件; |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +6 | Becho359 2026-03-11 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +10 | ADT 2026-03-11 | 10/500 |
|
|
[考研] 266求調(diào)劑 +3 | 陽陽哇塞 2026-03-14 | 7/350 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|