| 查看: 3048 | 回復: 12 | ||
[求助]
如何得到原子中某個軌道對分子總磁矩的貢獻? 已有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 但我想知道某個具體軌道對磁矩的貢獻,知道怎么分析嗎?例如下圖中Au的5d軌道對總磁矩的貢獻怎么分析?謝謝! QQ截圖20131224102322.jpg |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
體系有n個自旋平行的未配對電子,分子電子自旋磁矩就為n*μ_B。而自旋相反的電子的自旋磁矩會抵消,所以閉殼層體系無電子自旋磁矩。通過布居分析(如Mulliken、NPA等),得到每個原子,或每個原子軌道上的自旋電子數(shù),就能討論它們各自對分子磁矩的貢獻。比如某d原子軌道上自旋電子數(shù)為0.4,就說這個d原子軌道上電子自旋磁矩是0.4μ_B。也可以根據(jù)電子自旋密度討論,某點自旋密度為x,就可以說這個點的電子對分子磁矩的貢獻是x*μ_B。也因此,可以將空間劃分為不同子區(qū)域,根據(jù)自旋密度在相應區(qū)域積分值討論相應區(qū)域內(nèi)的自旋磁矩。μ_B=e*h_bar/(2m_e)這個量被定義為玻爾磁子。 高斯布居分析只是給出每個基函數(shù)的布居數(shù),怎么和原子軌道對應上需要基于你對所用基組的了解。pop=NPA給出的NPA分析結(jié)果則可以直接給出每個原子軌道的布居數(shù)。由于你用了贗勢,所以NPA分析輸出的哪個是5d需要自行判斷。 另外,分子磁矩是電子自旋磁矩、電子軌道運動產(chǎn)生的磁矩和原子核自旋磁矩共同貢獻的,通過電子布居分析討論的只是電子自旋磁矩部分,一般總磁矩也是由此主導。 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
謝謝! 一看就知道您對這方面很了解。我剛接觸這方面的內(nèi)容,按我對您回復的理解,以下幾點是否正確: 1.NPA得到各原子的各電子層上的電子布居數(shù)n,相應的磁矩就是n*μ_B了?某原子的磁矩就是各軌道磁矩之和,各原子磁矩之和就是團簇的總磁矩? 2.http://www.ccl.net/chemistry/res ... .002-dir/index.html 中說是根據(jù)Gross orbital populations之后的數(shù)據(jù)來計算磁矩,現(xiàn)在主流的是采用“NPA得到各原子的各電子層上的電子布居數(shù)n,相應的磁矩就是n*μ_B”還是"Gross orbital populations之后的數(shù)據(jù)來計算磁矩"思路。 3.網(wǎng)上查詢了一下,也有用MS計算磁矩的, 4.您自己有沒有這方面的已經(jīng)發(fā)表的文章,到時可以引用引用。 5.另外,您有沒有聽說過自旋污染? |
|
1.附件中分別是NBO計算的輸入和輸出文件,是不是總的“ NATURAL POPULATIONS: Natural atomic orbital occupancies” 減去Alpha spin orbitals 的“ NATURAL POPULATIONS: Natural atomic orbital occupancies”就是相應的beta 自旋軌道的,之后Alpha,beta spin orbitals 的電子布局數(shù)之差就是相應原子軌道的未配對電子? 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.若不對,能否在輸出文件中指出,萬分謝謝! |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
1 不是這么簡單相乘。比如某軌道上alpha電子布居數(shù)是na,beta是nb,那么自旋布居就是na-nb,(na-nb)*μ_B才是這個軌道的電子貢獻的磁矩。后面兩句話你的理解正確。 2 高斯輸出的Gross orbital populations是基于Mulliken方法得到的基函數(shù)的布居數(shù),這種計算方式很不推薦使用。尤其是帶了彌散函數(shù)的情況更不能用(比如CCL上那個人拿6-311+G*舉例就很不合適),NPA得到的軌道布居數(shù)有意義得多。Mulliken方法的缺陷在《談談軌道成份的計算方法》(http://hi.baidu.com/sobereva/item/d52876160110ca8d89a956c5)、《原子電荷計算方法的對比》(http://www.whxb.pku.edu.cn/EN/10.3866/PKU.WHXB2012281)上都有討論 3 問題不在于程序,而在于原理 4 我的文章里沒有專門討論磁矩的 5 自旋污染幾乎是眾所周知的問題 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
您好!我想追問一個磁矩的問題;我想得到AuC中性團簇中Au的5d,6s,6p軌道對總磁矩的貢獻,按照您的思路,alpha軌道電子布居為0.6481,beta軌道電子布局為0.6759,總自旋布居為0.6481-0.6759=-0.0278.總布局為負值,您如何理解?附件中為輸出文件; |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +4 | 手機用戶 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考博] 環(huán)境領域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 3/150 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風 2026-03-17 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 材料與化工一志愿南昌大學327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[考研] 0703 物理化學調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|