周　帥，方志剛，崔遠(yuǎn)東，張　偉，徐詩(shī)浩，劉　琪，馮　天，李雯博

（遼寧科技大學(xué) 化學(xué)工程學(xué)院，遼寧 鞍山　114051）

原子簇是指由原子（或分子）結(jié)合在一起的團(tuán)體結(jié)構(gòu)，它是介于原子（或分子）與固體粒子之間的團(tuán)粒分子。近幾十年來(lái)由于金屬原子簇化合物的催化功能［1-2］、生物金屬原子簇［3-4］、超導(dǎo)及新型材料［5］等方面的研究需要，促使金屬原子簇化學(xué)快速發(fā)展。金屬原子簇具有較高的耐磨及抗腐蝕性［6］、良好的催化活性［7-8］、優(yōu)異的磁學(xué)性能［9-12］及超導(dǎo)性能［13］等，使其在制備新型材料［14］、催化劑［15］、儲(chǔ)氫材料［16］及超導(dǎo)材料［17-18］等領(lǐng)域得到廣泛的應(yīng)用。過(guò)渡金屬元素與類金屬元素硼組成的團(tuán)簇是目前的研究熱點(diǎn)，對(duì)Mn-B二元體系和含Mn、B、P多元體系研究?jī)H有少量報(bào)道［19-22］，對(duì)Mn-B-P三元體系的研究還未見(jiàn)相關(guān)報(bào)道。本文將對(duì)Mn-B-P三元體系的原子簇進(jìn)行研究，以Mn3BP為模型，對(duì)團(tuán)簇Mn3BP的電子自旋密度進(jìn)行較深層次的研究，從而推斷其相關(guān)性質(zhì)，希望能給相關(guān)的研究做出一些有價(jià)值的參考。

1　計(jì)算方法

首先根據(jù)拓?fù)鋵W(xué)原理和化學(xué)鍵理論［23-24］，設(shè)計(jì)出了團(tuán)簇Mn3BP二、四重態(tài)的所有可能構(gòu)型。之后利用含相關(guān)校正的DFT［25］方法，在B3LYP/Lan12dz水平下對(duì)原子簇Mn3BP的所有可能構(gòu)型進(jìn)行全參數(shù)優(yōu)化和頻率計(jì)算，得到其穩(wěn)定的構(gòu)型。計(jì)算時(shí)對(duì)金屬錳原子采用Hay［26］等人的含相對(duì)論校正的有效核電勢(shì)價(jià)電子從頭計(jì)算基組，即采用18-eECP雙ξ基組（3s，3p，3d/2s，2p，2d），且P原子加加極化函數(shù)ξP.d=0.55［27］；所有計(jì)算均在啟天M7150微機(jī)上采用Gaussian09程序完成。

2　結(jié)果與討論

2.1　穩(wěn)定構(gòu)型

采用上述模型和計(jì)算方法分別對(duì)團(tuán)簇Mn3BP二、四重態(tài)構(gòu)型進(jìn)行優(yōu)化計(jì)算，排除含有虛頻的不穩(wěn)定構(gòu)型和幾何構(gòu)型相同的構(gòu)型，最終得到四重態(tài)穩(wěn)定構(gòu)型5種，如圖1中的1（4）～5（4）（右上角括號(hào)數(shù)字表示重態(tài)，下同），得到二重態(tài)穩(wěn)定構(gòu)型4種，如圖1中1（2）～4（2），二、四重態(tài)構(gòu)型分別按能量由低到高排序。圖1中各構(gòu)型中的數(shù)字表示各鍵鍵長(zhǎng)（單位：nm），為方便研究，以能量最低的構(gòu)型1（4）作為基準(zhǔn)，令其能量為0 kJ/mol，然后依次得到其它構(gòu)型的相對(duì)能量（括號(hào)里的數(shù)值表示相對(duì)能量）。

圖1中可以看出，各穩(wěn)定構(gòu)型的幾何構(gòu)型有戴帽三角錐（1（4）和 4（4）），三角雙錐（5（4）、1（2）和 2（2）），平面四邊形（3（4）和3（2））和平面五邊形（2（4）和4（2））。構(gòu)型1（4）是以Mn1-Mn2-Mn3為基準(zhǔn)面，B為錐頂，P為帽頂?shù)拇髅比清F；構(gòu)型4（4）是以B-Mn1-Mn3為基準(zhǔn)面，P為錐頂，Mn2為帽頂?shù)拇髅比清F。構(gòu)型5（4）和 2（2）均是以 Mn1-Mn2-Mn3 為基準(zhǔn)面，P 為錐頂，B為錐底的三角雙錐，優(yōu)化構(gòu)型相同，多重度不同；構(gòu)型1（2）是以Mn1-Mn3-B為基準(zhǔn)面，P為錐頂，Mn2為錐底的三角雙錐。構(gòu)型 3（4）和 3（2）均是B原子作為中心原子，四邊頂點(diǎn)以Mn1-Mn2-P-Mn3的排列順序構(gòu)成的平面四邊形。構(gòu)型2（4）和4（2）是由Mn1、Mn2、Mn3、B、P構(gòu)成的平面五邊形。

2.2　各構(gòu)型的自旋密度分布

構(gòu)型的穩(wěn)定性有很多影響因素，電子自旋密度分布是其中的重要影響因素，分析各構(gòu)型的自旋密度分布可以幫助分析各構(gòu)型的穩(wěn)定性。表1列出了各優(yōu)化構(gòu)型原子的電子自旋密度值，正值的含義是α成單電子出現(xiàn)的凈概率密度，負(fù)值則代表β成單電子出現(xiàn)的凈概率密度。分析表1中數(shù)據(jù)可知，所有構(gòu)型的B原子的電子自旋密度均為負(fù)值，而P原子的電子自旋密度也均為負(fù)值（除了構(gòu)型1（4）外），由此可知各構(gòu)型中（除了構(gòu)型1（4）外）的非金屬原子B、P上的電荷分布都是自旋向下的β電子。從整體來(lái)看，所有構(gòu)型中（除構(gòu)型5（4）外）有一個(gè)Mn原子的電子自旋密度值為負(fù)數(shù)，其余兩個(gè)Mn原子的電子自旋密度為正數(shù)，也就表示每個(gè)構(gòu)型中（除構(gòu)型5（4）外）有一個(gè)Mn原子的電荷分布是自旋向下的β電子，其余兩個(gè)Mn原子的電荷分布是自旋向上的α電子，構(gòu)型4（4）、5（4）、2（2）、3（2）有兩個(gè)Mn原子的電荷分布是自旋向上的α電子，且電子自旋密度值相同或近似相等。

圖1　團(tuán)簇Mn3BP優(yōu)化構(gòu)型圖Fig.1　Optimized configurations of cluster Mn3BP

表1　團(tuán)簇Mn3BP各構(gòu)型原子的電子自旋密度值Tab.1　Spin density of each configuration atom in cluster Mn3BP

只從原子的電子自旋密度分析，存在很大的局限性，并不能得到很準(zhǔn)確的結(jié)論。為了能更好地研究電子自旋密度對(duì)構(gòu)型穩(wěn)定性的影響，需要分析α電子、β電子及其重疊部分對(duì)其穩(wěn)定性的影響。表2列出了團(tuán)簇Mn3BP的各構(gòu)型原子間的電子自旋密度值，正值代表原子間成鍵時(shí)α電子過(guò)剩，負(fù)值代表原子間成鍵時(shí)β電子過(guò)剩。為更方便地研究自旋密度分布對(duì)構(gòu)型穩(wěn)定性的影響，圖2畫出了團(tuán)簇Mn3BP的各構(gòu)型相對(duì)應(yīng)的自旋密度分布圖，按能量高低排序，圖中淺灰色代表α電子，深灰色代表β電子。

依據(jù)表2的電子自旋密度值來(lái)分析內(nèi)部原子間α、β電子分布情況，及其引起的構(gòu)型各原子間的成鍵強(qiáng)弱和均勻度，從而影響各構(gòu)型的穩(wěn)定性。所有構(gòu)型中非金屬原子B與P之間的電子自旋密度值全都為正值，說(shuō)明非金屬原子間成鍵時(shí)α電子過(guò)剩；又通過(guò)表1知所有構(gòu)型中（除構(gòu)型1（4）外）的非金屬原子B、P的電子分布為β電子；以上兩點(diǎn)說(shuō)明了B、P之間成鍵時(shí)的電子很大程度上受金屬原子Mn的影響。結(jié)合表2數(shù)據(jù)和圖2分析，構(gòu)型1（4）與其他構(gòu)型相比，構(gòu)型原子間電子分布密度大且較為均勻，有效降低該構(gòu)型的能量，使構(gòu)型1（4）更加的穩(wěn)定；構(gòu)型4（2）的原子間電子分布密度小且不均勻，且對(duì)稱性也不好，導(dǎo)致該構(gòu)型的能量高，穩(wěn)定性差。

為更好地尋找電子自旋密度分布對(duì)構(gòu)型穩(wěn)定影響的規(guī)律，將各構(gòu)型的電子自旋密度分布圖進(jìn)行分類討論，大致將其分為了三類。第一類：全部構(gòu)型中最為穩(wěn)定的前兩種構(gòu)型，包括構(gòu)型1（4）和構(gòu)型2（4）。這兩種構(gòu)型的電子自旋密度分布比較均勻，且具有對(duì)稱性較好，構(gòu)型1（4）的B-Mn原子間只有B-Mn1有β原子過(guò)剩，P-Mn原子間只有P-Mn2有β原子過(guò)剩，Mn-Mn原子間只有Mn1-Mn2有β原子過(guò)剩；構(gòu)型2（4）的B-Mn原子間只有B-Mn3有α電子過(guò)剩，P-Mn原子間只有P-Mn2有α電子過(guò)剩，Mn-Mn原子間只有β電子過(guò)剩；構(gòu)型 1（4）的 P-Mn2和P-Mn3之間的電子自旋密度值絕對(duì)值近似相等（P-Mn2：-0.095；P-Mn3：0.092），Mn1-Mn2和Mn2-Mn3之間的電子自旋密度值的絕對(duì)值也近似相等（Mn1-Mn2：-0.125；Mn2-Mn3：0.126）；構(gòu)型 2（4）的Mn1-Mn2和Mn2-Mn3的電子自旋密度值近似相等（Mn1-Mn2：-0.026；Mn2-Mn3：-0.029）；以上幾點(diǎn)說(shuō)明了構(gòu)型1（4）的原子間電子自旋密度分布比構(gòu)型2（4）更均勻，原子間成鍵強(qiáng)度大且均勻，能量較低，比構(gòu)型 2（4）的穩(wěn)定性好。

表2　團(tuán)簇Mn3BP各構(gòu)型原子間的電子自旋密度值Tab.2　Spin density between atom and atom in cluster Mn3BP

圖2　團(tuán)簇Mn3BP各優(yōu)化構(gòu)型的電子自旋密度分布圖Fig.2　Spin density distribution map of each configuration in cluster Mn3BP

第二類：多重度不同幾何構(gòu)型很相似的構(gòu)型3（4）和3（2），還有多重度不同幾何構(gòu)型相同的構(gòu)型5（4）和2（2）。構(gòu)型3（4）的B-Mn、P-Mn、Mn-Mn原子間成鍵時(shí)各有一個(gè)為α電子過(guò)剩，兩個(gè)為β電子過(guò)剩；構(gòu)型3（2）的B-Mn原子間成鍵時(shí)有兩個(gè)為α電子過(guò)剩，一個(gè)為β電子過(guò)剩，P-Mn原子間成鍵時(shí)有一個(gè)為α電子過(guò)剩，兩個(gè)為β電子過(guò)剩，Mn-Mn原子間成鍵時(shí)都為β電子過(guò)剩；以上分析說(shuō)明了構(gòu)型3（4）的電子自旋密度分布比構(gòu)型3（2）的分布均勻。再結(jié)合自旋密度分布圖分析，構(gòu)型3（4）的α電子、β電子及兩者重疊部分分布較構(gòu)型 3（2）的更均勻，也證明了構(gòu)型3（4）的電子自旋密度分布更為均勻，使構(gòu)型的能量降低，更加的穩(wěn)定。構(gòu)型 5（4）和 2（2）構(gòu)型作為幾何構(gòu)型相同的構(gòu)型有許多相似之處，構(gòu)型5（4）的B與Mn1、Mn2原子間電子自旋密度值相等（B-Mn1：-0.045；B-Mn2：-0.045），P與Mn1、Mn2原子間電子自旋密度值相等（P-Mn1：-0.056；P-Mn2：-0.056），Mn3與Mn1、Mn2原子間電子自旋密度值相等（Mn1-Mn3：0.028：；Mn2-Mn3：0.028）；構(gòu)型 2（2）的 B 與Mn2、Mn3原子間電子自旋密度值相等（BMn2：-0.068；B-Mn3：-0.068），P與Mn2、Mn3原子間電子自旋密度值相等（P-Mn2：-0.041；PMn3：-0.041），Mn1與Mn2、Mn3原子間電子自旋密度值相等（Mn1-Mn2：-0.007；Mn1-Mn3：-0.007）；以上描述說(shuō)明了構(gòu)型 5（4）和 2（2）外部電子自旋密度分布比較均勻，但結(jié)合圖2來(lái)看構(gòu)型5（4）的原子間電子自旋密度分布要比構(gòu)型2（2）均勻，且對(duì)稱性更好，所以構(gòu)型5（4）就能量較低，更加穩(wěn)定?？偨Y(jié)以上說(shuō)法可得出結(jié)論，幾何構(gòu)型相同或相似但多重度不同的構(gòu)型四重態(tài)的構(gòu)型更加穩(wěn)定。

第三類：自旋密度分布不是很明顯的構(gòu)型，包括構(gòu)型4（4）、1（2）和 4（2）。構(gòu)型 4（4）的 B與Mn2 和 Mn3原子間電子自旋密度值相等（B-Mn2：0.018；BMn3：0.018），P與Mn2、Mn3原子間電子自旋密度值相等（P-Mn2：-0.039；P-Mn3：-0.039），Mn1與Mn2、Mn3原子間電子自旋密度值相等（Mn1-Mn2：0.005；Mn1-Mn3：0.005），而其余兩種構(gòu)型的原子間自旋密度值沒(méi)有相等的，說(shuō)明構(gòu)型4（4）的原子間電子自旋密度分布比較均勻，再結(jié)合圖2來(lái)看構(gòu)型4（4）的電子自旋密度分布的對(duì)稱性比較好，所以構(gòu)型 4（4）能量低，穩(wěn)定性好。構(gòu)型 1（2）和 4（2）相比，構(gòu)型4（2）的電子自旋密度分布比構(gòu)型 1（2）分布均勻，對(duì)稱性也比構(gòu)型 1（2）好，但是實(shí)際上構(gòu)型 1（2）的能量更低，穩(wěn)定性好，說(shuō)明電子自旋密度只是影響構(gòu)型穩(wěn)定性的重要因素之一，還有其它因素也影響構(gòu)型穩(wěn)定性。

3　結(jié)論

團(tuán)簇Mn3BP的優(yōu)化構(gòu)型有：戴帽三角錐（1（4）和4（4）），三角雙錐（1（2）、2（2）和 5（4）），平面四邊形（3（2）和3（4））和平面五邊形（4（2）和2（4））；各構(gòu)型中（除了構(gòu)型1（4）外）的非金屬原子B、P上的電荷分布都是自旋向下的β電子，B-P原子間電子自旋密度均為正值，α電子過(guò)剩，B-P原子間成鍵時(shí)電子分布受金屬M(fèi)n原子影響；幾何構(gòu)型相同或相似多重度不同的構(gòu)型，四重態(tài)電子自旋密度分布均勻，對(duì)稱性較好，能量較低，穩(wěn)定性好；電子自旋密度是影響構(gòu)型穩(wěn)定性的因素之一，還有其它因素也影響構(gòu)型穩(wěn)定性。

參考文獻(xiàn)：

［1］LANG S M，BERNHARDT T M.Gas phase metal cluster model systems for heterogeneous catalysis.［J］.Physical Chemistry Chemical Physics，2012，14（26）：9255-9269.

［2］KULKARNI A，LOBOLAPIDUS R J，GATES B C.Metal clusters on supports：synthesis，structure，reactivity，and catalytic properties.［J］.Chemical Communications，2010，46（33）：5997-6015.

［3］SHANG L，DONG S，NIENHAUS G U.Ultra-small fluorescent metal nanoclusters：synthesis and biological applications［J］.Nano Today，2011，6（4）：401-418.

［4］MATHEW A，PRADEEP T.Noble metal clusters：applications in energy，environment，and biology［J］.Particle&Particle Systems Characterization，2015，31（10）：1017-1053

［5］BONO D，BAKHAREV O N，SCHNEPF A，et al.Magnetization studies of superconductivity in a molecular metal cluster compound［J］.ZeitschriftFürAnorganische Und AllgemeineChemie，2010，633（13-14）：2173-2177.

［6］SOUZA C A C D，BOLFARINI C，BOTTA W J，et al.Corrosionresistanceandglassformingabilityof Fe47Co7Cr15M9Si5B15Y2（M=Mo，Nb）amorphous alloys［J］.Materials Research，2013，16（6）：1294-1298.

［7］CHENG Z，XU J，ZHONG H，et al.Amorphous Co-B/Al nanocomposites prepared by electroless coating technique and their excellent catalytic activity［J］.Rare Metals，2010，29（1）：32-36.

［8］ZHANG J，YAO H R，LI H B，et al.Amorphous alloy ni-fe-b：preparation，characterization and its catalytic ability for hydrogen generation［J］.Applied Mechanics&Materials，2014，672-674：592-596.

［9］YUAN H K，CHEN H，KUANG A L，et al.Structural and magnetic properties of small 4d transition metal clusters：role of spin-orbit coupling［J］.Journal of Physical Chemistry A，2012，116（47）：11673-11684.

［10］GUO P，ZHENG J，GUO X，et al.Electronic and magnetic properties of transition-metal-doped sodium superatom clusters：TM@Na 8，（TM=3d，4d and 5d transition metal）［J］.Computational Materials Science，2014，95：440-445.

［11］LIU D，CHEN Z，HUANG W，et al.Metal clusters constructed by dihydrazide ligand：Synthesis，structure and magnetic properties［J］.InorganicaChimicaActa，2013，400（5）：179-183.

［12］LIANG W，JIA J，LV J，et al.Electronic structure，stability and magnetic properties of small M 1-2 Cr（M=Fe，Co，and Ni）alloy encapsulated inside a（BN）48，cage［J］.Physics Letters A，2015，379（30-31）：1715-1721.

［13］BONO D，BAKHAREV O N，SCHNEPF A，et al.Magnetization studies of superconductivity in a molecular metal cluster compound［J］.ZeitschriftFürAnorganische Und Allgemeine Chemie，2010，633（13-14）：2173-2177.

［14］CUERVA M，GARCíA-FANDI?O R，VáZQUEZVáZQUEZ C，et al.Self-assembly of silver metal clusters of small atomicity on cyclic peptide nanotubes［J］.ACS nano，2015，9（11）：10834-10843.

［15］LI L，LARSEN A H，ROMERO N A，et al.Investigation of catalytic finite-size-effects of platinum metal clusters［J］.Journal of Physical Chemistry Letters，2013，4（1）：222-226.

［16］GIRI S，CHAKRABORTY A，CHATTARAJ P K.Potential use of some metal clusters as hydrogen storage materials—a conceptual DFT approach［J］.Journal of Molecular Modeling，2011，17（4）：777-784.

［17］RODUNER E，JENSEN C.Magnetic properties and the superatom character of 13-atom platinum nanoclusters［J］.Magnetochemistry，2015，1（1）：28-44.

［18］PAGLIONE J，GREENE R L.High-temperature superconductivity in iron-based materials［J］.Nature physics，2010，6（9）：645-658.

［19］SUN W，DU Y，LIU S，et al.Thermodynamic assessment of the Mn-B system［J］.Journal of Phase Equilibria&Diffusion，2010，31（4）：357-364.

［20］KUPCZYK A，?WIERCZEK J，HASIAK M，et al.Microstructure and some thermomagnetic properties of amorphous Fe-（Co）-Mn-Mo-B alloys［J］.Journal of Alloys and Compounds，2018，735：253-260.

［21］ZHENG Z G，ZHU Z R，YU H Y，et al.Large magnetic entropy change and magnetic phase transitions in rapidly quenched bulk Mn-Fe-P-Si alloys［J］.Journal of Alloys and Compounds，2017，725：1069-1076.

［22］REDDI P D，MUKHOPADHYAY N K，MAJUMDAR B，et al.Synthesis of Fe-Si-B-Mn-based nanocrystalline magnetic alloys with large coercivity by high energy ball milling［J］.Bulletin of Materials Science，2014，37（4）：815-821.

［23］NOURY S，COLONNA F，SAVIN A，et al.Analysis of the delocalization in the topological theory of chemical bond［J］.Journal of Molecular Structure，1998，450（1-3）：59-68.

［24］ZHANG L，YING F，WU W，et al.Topology of electron charge density for chemical bonds from valence bond theory：a probe of bonding types［J］.Chemistry-A European Journal，2009，15（12）：2979-2989.

［25］ORIO M，PANTAZIS D A，NEESE F.Density functional theory［J］.Photosynthesis research，2009，102（2-3）：443-453.

［26］HAY P J，WADT W R.Ab initio effective core potentials for molecular calculations.Potentials for K to Au including the outermost core orbitals［J］.The Journal of Chemical Physics，1985，82（1）：299-310.

［27］FANG Z G，HU H Z，GUO J X.Quantum chemical study on geometry and property of cluster Ni4P［J］.結(jié)構(gòu)化學(xué)，2006，25（1）：7-16.