高嶺土結(jié)構(gòu)的多磁場27Al NMR 波譜和ab initio量子理論計(jì)算研究

蘇紹明, 鄒 珍

(中國冶金地質(zhì)總局 中南地質(zhì)勘查院，湖北 武漢 430081)

高嶺土結(jié)構(gòu)的多磁場27Al NMR 波譜和ab initio量子理論計(jì)算研究

蘇紹明, 鄒 珍

(中國冶金地質(zhì)總局 中南地質(zhì)勘查院，湖北 武漢 430081)

自上世紀(jì)80年代以來,試圖使用魔角旋轉(zhuǎn)固體核磁共振(MAS NMR)波譜來分辨高嶺土結(jié)構(gòu)兩個(gè)不同27Al結(jié)晶位置的多次嘗試均未獲成功,所以本次研究試圖通過結(jié)合量子理論計(jì)算和多場MAS NMR和MQMAS波譜再次嘗試分辨出高嶺土結(jié)構(gòu)中的不同27Al位置。超高場如21 T所引起的高各向異性化學(xué)位移作用(CSA)可能會(huì)超過四極核作用(QI),這對(duì)于具較小Cq的27Al位置(如高嶺土結(jié)構(gòu))而言尤為如此。而低場如4.68 T中,QI則是主宰的核磁共振作用,但是,由于高嶺土結(jié)構(gòu)中兩個(gè)27Al位置的環(huán)境很相似,所以,通常27Al MAS NMR也無法分辨出這兩個(gè)27Al 位置。結(jié)合完全勢能線型增廣平面波(FP LAPW)+局域軌道法(lo)、分子軌道(MO)的abinitio量子理論計(jì)算,從多量子魔角旋轉(zhuǎn)波譜、衛(wèi)星躍遷波譜中可得到高嶺土結(jié)構(gòu)中兩個(gè)27Al位置的準(zhǔn)確核磁共振參數(shù)如Cq、η、CSA 和δiso等。

高嶺土;27Al MAS NMR;MQMAS;SATRAS;FP LAPW理論計(jì)算模型;CASTEP;分子軌道計(jì)算模型

1 概述

固體魔角旋轉(zhuǎn)核磁共振 (MAS NMR)波譜能提供各向同性化學(xué)位移 (δiso)、各向異性 (CSA),偶極相互作用(DD)和四極核作用 (QI)等信息。這些核磁共振參數(shù)能對(duì)應(yīng)著結(jié)構(gòu)特征如鍵長、鍵角、配位類型、聚合度和原子間的替換[1-3],所以MAS NMR可以提供礦物在原子、分子尺度上的精細(xì)局域結(jié)構(gòu)信息。例如,鋁硅酸鹽的27Al NMR參數(shù)與配位數(shù)、多面角畸變和聚合度間的相關(guān)性已經(jīng)被成功建立[4]。但是,礦物中的多結(jié)構(gòu)位置、原子間的替換所造成的無序性常造成NMR譜峰的疊合,從而降低波譜分辨率、無法得出準(zhǔn)確的各種NMR參數(shù)。

由于非球形分布的核電荷,自旋數(shù)I>1/2的四極核擁有四極核動(dòng)量(eQ)[5]。eQ與原子核所在位置的電場梯度(EFG)間的靜電作用導(dǎo)致四極核作用,可用兩個(gè)參數(shù)來描述這樣的作用,即四極耦合常數(shù)(Cq)和非對(duì)稱參數(shù)(η),這兩個(gè)參數(shù)可以通過原子核所在位置的EFG張量來計(jì)算:

Cq=eVzzQ/h

(1)

(2)

式中:Vzz、Vxx、Vyy是主軸坐標(biāo)系統(tǒng)(PAS)的EFG 分量,e是電子電荷單位;h普朗克常數(shù)除以2π;Q是核四極動(dòng)量 (27Al 四極動(dòng)量是14.6 fm2)。 QI作用常導(dǎo)致很寬的線型和很低的分辨率,對(duì)多位置疊合峰尤為如此。

高嶺土結(jié)構(gòu)已被詳細(xì)研究[6-10],但是MAS NMR波譜一直沒能分辨出其結(jié)構(gòu)中兩個(gè)不同的27Al結(jié)晶位置[11-14]。 本文試圖使用多場NMR波譜并結(jié)合量子理論計(jì)算,來分辨出這兩個(gè)不同Al結(jié)晶位置,如仍不能分辨,則將試圖解釋NMR波譜缺乏分辨率的原因。這些NMR實(shí)驗(yàn)包括MAS NMR、多量子魔角旋轉(zhuǎn)(MQMAS)NMR、衛(wèi)星躍遷(SATRAS)NMR,希望這些實(shí)驗(yàn)?zāi)軌驕p弱二階QI和其它引起線寬的因素,從而提高NMR波譜分辨率、有助于提取準(zhǔn)確的Cq,η和δiso等參數(shù)[15-20]。完全勢能線性增廣平面波法(FP LAPW)和分子軌道模型(MO)將用于計(jì)算和預(yù)測結(jié)構(gòu)的NMR參數(shù),并能用來優(yōu)化結(jié)構(gòu)、歸屬NMR譜峰等。

1.1 理論計(jì)算

FP LAPW計(jì)算模型基于周期性晶體結(jié)構(gòu)的對(duì)稱和晶胞內(nèi)原子分布,通過密度函數(shù)(DFT)來處理晶體內(nèi)復(fù)雜多體電子體系的交換和相關(guān)作用,而晶胞內(nèi)勢能通過 Muffin-Tin 模型來模擬：即晶胞中勢能分布可劃分成兩種區(qū)域,一種是圍繞每個(gè)原子的球形勢能區(qū),而球間區(qū)域則是類似平面的勢能面。 隨后,由于晶體的周期性結(jié)構(gòu)特征,使用Bloch原理和公式、選擇完全和有效的基函數(shù)(平面波函數(shù)或球面函數(shù))張成基底,并通過Kohn-Sham (KS)公式將晶體的多體問題簡化成單體問題,經(jīng)過自洽方法最終解出薛定諤方程,得到晶體波函數(shù)[21],從而可以計(jì)算出晶體的各種物理化學(xué)性質(zhì)(如電子結(jié)構(gòu)Cq和η),可通過軟件WIEN2k來實(shí)施這樣的計(jì)算[22-26]。

FP LAPW計(jì)算模型對(duì)晶胞內(nèi)原子位置十分敏感,所以,通過調(diào)整晶體XRD結(jié)構(gòu)的原子坐標(biāo)直至晶體系統(tǒng)能量和或原子所受作用力達(dá)最低,此時(shí)晶體結(jié)構(gòu)即達(dá)到平衡態(tài),如此即可優(yōu)化、精修晶體XRD結(jié)構(gòu)[27-29]。NMR實(shí)驗(yàn)結(jié)果如QI參數(shù)可證實(shí)優(yōu)化的精修結(jié)構(gòu)。

CASTEP也是一種abinitio量子計(jì)算軟件,可以計(jì)算和預(yù)測周期性結(jié)構(gòu)的NMR參數(shù)[30-33]。

對(duì)于有著無序特征的晶體結(jié)構(gòu),可使用分子軌道(MO)模型來計(jì)算結(jié)構(gòu)的物理、化學(xué)特性 。對(duì)晶體MO計(jì)算模型而言,首先選擇目標(biāo)原子如Al,在目標(biāo)原子一定距離范圍之內(nèi)的原子組成分子族,設(shè)定收斂條件后,使用MO計(jì)算該分子族的波函數(shù),從而可計(jì)算出相應(yīng)的物理、化學(xué)性質(zhì)[34]。同樣地,可以使用多種方法來近似原子族內(nèi)電子間的交換、相關(guān)作用如Hartree-Fock (HF)和密度函數(shù)[35]。

1.2 高嶺土晶體結(jié)構(gòu)

通過XRD、中子散射確定和精修了高嶺土結(jié)構(gòu),并且最低能量模型也給出了結(jié)構(gòu)中H原子的位置[36-37]。高嶺土晶體結(jié)構(gòu)空間群是C1,晶胞參數(shù)分別為a=5.153 5(3)?,b=8.941 9(5)?,c=7.390 6(4)?,α=91.926(6)°,β=105.046(2)°,γ=89.797(2)°,晶胞體積V=328.70(5)?3。高嶺土具有序的中性1∶1 片層結(jié)構(gòu),片層由Al(O,OH)6八面體和SiO4四面體所組成,片層通過氫鍵連接在一起。

2 實(shí)驗(yàn)方法

2.1 試驗(yàn)樣品

兩個(gè)高嶺土樣品均為白色粉末狀,粉晶XRD 實(shí)驗(yàn)證實(shí)了兩個(gè)高嶺土樣品的純度。室溫下高嶺土樣品粉末XRD數(shù)據(jù)由計(jì)算機(jī)控制的ADP-2衍射儀采集,使用 CuKα 作為射線源,2的步長為0.02°、變化范圍是10°～70°,每步計(jì)數(shù)時(shí)間為5 s。

2.227Al MAS NMR

所有27Al波譜均以 [Al(H2O)6]3+定標(biāo)。最低場27Al MAS NMR 實(shí)驗(yàn)使用4.68 T、7 mm探頭的Tegmac-200 譜儀,Si3N4轉(zhuǎn)子、轉(zhuǎn)速為8 kHz,激發(fā)脈沖寬度為 1 μs (18°)、掃描次數(shù)為 8、d1時(shí)間為0.5 s、系統(tǒng)死時(shí)間是 8 μs、1H去偶功率為 50 kHz。最高磁場NMR實(shí)驗(yàn)使用21.06 T的Bruker Avance II 900 NMR譜儀,樣品轉(zhuǎn)速15 kHz。在21.06 T譜儀上使用標(biāo)準(zhǔn)3QMAS脈沖序列來采集MQMAS波譜,這樣的脈沖序列由兩個(gè)硬脈沖及其后的z-filter脈沖所組成,3Q激發(fā)脈沖寬度是4.5 μs、3Q轉(zhuǎn)換1Q的脈沖寬度1.5 μs、中心躍遷選擇脈沖寬度是12.5 μs。21.06 T場1H去偶功率為100 kHz。

在7.1 T、9.4 T 和14.04 T場的Varian INOVA譜儀上運(yùn)行了1H-去偶27Al MAS NMR實(shí)驗(yàn),使用Si3N4轉(zhuǎn)子、12 kHz轉(zhuǎn)速、1 kW的去偶功率。在Varian INOVA-300譜儀上使用z-filter脈沖序列也采集了3QMAS波譜[38-43]。3Q激發(fā)和3Q-1Q轉(zhuǎn)換脈沖寬度分別是4.0和1.0s,B1場強(qiáng)是60 kHz,90°選擇脈沖寬度是20s、d1為 0.5 s。共采集128次t1增量、480 次掃描。

2.3 波譜擬合程序

DMFIT用于擬合 21T場的MAS NMR波譜和MQMAS切面波譜。δiso、Cq和η參數(shù)由手動(dòng)輸入直到擬合波譜與實(shí)驗(yàn)波譜相吻合。

2.4 FP LAPW計(jì)算

高嶺土結(jié)構(gòu)的晶胞參數(shù)、晶胞內(nèi)原子坐標(biāo)、空間群是FP LAPW WIEN2k計(jì)算的輸入?yún)?shù)。計(jì)算使用cutoff值6.5、200 k點(diǎn),在沒有電荷泄露的前提下能高效運(yùn)行FP LAPW計(jì)算的Al、O、Si、H原子球半徑Rmt分別是1.8、1.6、1.4和 0.7 a.u.,原子球內(nèi)波函數(shù)包含了高達(dá)l=12的角動(dòng)量組分。采用了并行式超算模式,FP LAPW計(jì)算自洽收斂條件為 0.0001 Ry,且為非自旋極化計(jì)算模式。

2.5 CASTEP計(jì)算

使用Bish(1993)的高嶺土結(jié)構(gòu),基于密度函數(shù)理論,在Materials Studio 4.3環(huán)境下使用CASTEP NMR模塊計(jì)算NMR參數(shù),并使用超軟膺勢、平面波基函數(shù)(cut-off值500 eV)、PBE泛函近似電子間的交換—相關(guān)作用,計(jì)算的化學(xué)屏蔽值為絕對(duì)值。

3 結(jié)果

粉晶XRD揭示高嶺土樣品中只有微量雜質(zhì)相(<0.2 wt %),這些雜質(zhì)不會(huì)干擾樣品的NMR波譜。

在最低場4.68 T中,高嶺土中心躍遷(CT)的27Al MAS NMR 實(shí)驗(yàn)波譜為半高寬約75×10-6、重心(CG)約-12.5×10-6、高場位置有肩部的非對(duì)稱譜峰線型(圖1-a)。在7.1 T中,CT 線型(CG 約-4.0×10-6)半高寬減小到30×10-6,但仍保持一些非對(duì)稱特征,仍有高場方向的肩部 (圖1-b)。11.7 T場CT線型的CG已至0.5×10-6,且線型更加平滑。14.04 T場時(shí),CT線型(CG約4.8×10-6)就變成沒有特征的Gaussian線型(圖1-c)。21.06 T時(shí),CT譜峰為很窄(半高寬5×10-6)的Gaussian線型,且CG也移至6.2×10-6(圖1-d)。

圖1 高嶺土多場27Al MAS NMR波譜Fig.1 27Al MAS NMR spectra of kaolinite(a).4.68 T(CG=-12.5×10-6);(b).7.10 T (CG=-4.0×10-6);(c).14.04 T(CG=4.8×10-6);(d).21.06 T(CG=6.2×10-6)。

圖2 高嶺土27Al 衛(wèi)星躍遷MAS NMR實(shí)驗(yàn)波譜Fig.2 27Al SATRAS spectra acquired by MAS NMR(a).21T場、15 kHz 轉(zhuǎn)速;(b).11.7T場、12 kHz轉(zhuǎn)速。 插圖為(*±1/2～±3/2、**±3/2～±5/2)的旋轉(zhuǎn)邊帶放大圖。

11.7 T場中(圖2)使用DMFIT擬合的高嶺土1H去偶27Al MQMAS NMR結(jié)果列入表1,其中27Al(1)位置的δiso、Cq分別是(7.55±0.5)×10-6、3.6 MHz;27Al(2)位置則為9.0×10-6、(4.0±0.1)MHz。27Al(1)位置的η是0.57,但27Al(2)位置的η因不同擬合而給出分別為0.42和1.00不同的結(jié)果。而21 T場MQMAS波譜切面擬合的Cq和η平均值分別約1.7 MHz和0.6,這樣的擬合結(jié)果與同場 MAS波譜的擬合結(jié)果相差較大。

表1 11.7 T場高嶺土1H去偶27Al MQMAS 波譜切面使用DMFIT擬合結(jié)果

4.68、7.1、9.4、11.04和 21 T場衛(wèi)星躍遷的MAS波譜每條旋轉(zhuǎn)邊帶可分辨出一對(duì)譜峰(圖 2-a,圖2-b),但7.1 T 和 21 T 場的27Al MQMAS (3Q)NMR 仍沒能分辨出不同Al位置 (圖3)。

圖3 21T高嶺土樣品1H去偶 27Al 3QMAS NMR波譜(轉(zhuǎn)速15 kHz)Fig.3 1H-decoupled 27Al 3QMAS NMR spectrum at 21T with spinning speed of 15 kHz

使用WIEN2k 的FP LAPW 計(jì)算給出27Al(1)、27Al(2)位置Cq和η分別為-3.15 MHz、0.55和2.14 MHz、0.98(表2)。高嶺土中兩種Al位置的各向同性化學(xué)位移(δiso)、各向異性化學(xué)位移(δσ)及其對(duì)稱參數(shù)(ησ)的CASTEP計(jì)算結(jié)果也列入表2。

表2 高嶺土中兩個(gè)Al位置使用CASTEP計(jì)算的化學(xué)位移(δiso)、各向異性化學(xué)位移參數(shù)(δσ)和對(duì)稱參數(shù)(ησ),以及FP LAPW的QI參數(shù)計(jì)算值(Cq、ηQ)

4 結(jié)論

雖然XRD和中子散射識(shí)別出高嶺土中兩個(gè)不同Al結(jié)晶位置,但27Al NMR波譜仍沒能令人信服地分辨出這兩個(gè)位置[44-45]。

FP LAPW 和CASTEP 計(jì)算使用Bloch定理來考慮和處理晶體周期性結(jié)構(gòu)特征和邊界,使用DFT來近似體系內(nèi)電子間復(fù)雜的相互作用,并通過Kohn-Sham公式將多體問題轉(zhuǎn)化為單體問題,這樣,FP LAPW 和 CASTEP 模型能給出具周期性結(jié)構(gòu)的固體最準(zhǔn)確的EFG、NMR等物理化學(xué)特性的理論計(jì)算結(jié)果[46]。

高嶺土的27Al 每條旋轉(zhuǎn)邊帶的兩個(gè)峰(圖2-a和圖2-b)可能源自內(nèi)、外衛(wèi)星躍遷,即±5/2?±3/2 和±3/2?±1/2,也可能來自高嶺土結(jié)構(gòu)中的兩個(gè)不同Al位置。由于高嶺土中兩個(gè)Al位置的原子數(shù)之比為1∶1,且理論計(jì)算的Cq和η結(jié)果也十分相似,所以如果每條邊帶中的兩個(gè)譜峰源自這兩個(gè)不同Al位置,那么這兩個(gè)峰的寬度、強(qiáng)度應(yīng)該類似,但是每條邊帶中譜峰對(duì)寬度、強(qiáng)度的巨大差別排除了這種可能性,所以作者支持Rocha & Pedrosa de Jesus(1994)的分析,每條邊帶的譜峰對(duì)應(yīng)該源自衛(wèi)星躍遷,而非由不同Al位置所致。對(duì)I=5/2四極核而言,外部(±5/2?±3/2)和內(nèi)部(±3/2?±1/2)衛(wèi)星躍遷的寬度和四極化學(xué)位移分別為中心躍遷的1/3和1/9[47],這就解釋了每條邊帶的譜峰對(duì)的寬度、強(qiáng)度、位置對(duì)比。

如果引起線寬的因素如CSA、DD和QI 被MAS完全去除而僅留下δiso,那么這兩個(gè)Al位置就可能被NMR波譜分辨出來。MQMAS NMR實(shí)驗(yàn)就是通過操縱多量子相互干涉、多量子轉(zhuǎn)移途徑,來試圖去除二階QI,并抑制 CSA和 DD 作用,從而在F1維上只留下δiso信息,F2維度仍是通常的 MAS NMR波譜。但是,高嶺土中很強(qiáng)的同核和異核DD 作用在MQMAS實(shí)驗(yàn)過程中將被乘以量子干涉階數(shù)[48],從而會(huì)降低27Al MQMAS 波譜的分辨率。高嶺土中兩個(gè)Al位置δiso的相似性、高場中極強(qiáng)的 CSA作用也會(huì)降低MQMAS波譜的分辨率。所以,高嶺土的27Al MQMAS 波譜均未能分辨出不同的Al位置(圖3)。

Bish (1993)指出,高嶺土中可能存在H 原子站位的無序性,如此,類似于zunyite 的F/OH無序性,這也將是造成27Al MAS NMR CT譜峰疊合,引起NMR譜寬和缺乏典型QI線型[49]。

總之,由于高嶺土中兩個(gè)Al位置的相似性,很強(qiáng)的剩余CSA和DD作用,以及可能的 H 無序性,二階QI線型特征被覆蓋,導(dǎo)致高場中波譜無特征的Gaussian線型,從而造成需要使用很大的線寬因子(LB)來擬合27Al CT 波譜線型,這就導(dǎo)致了擬合過程對(duì)η和Cq變化的敏感度降低。所以,高嶺土二階CT MAS波譜和MQMAS切面波譜的擬合Cq和η值之間,以及與理論計(jì)算值存在著不一致性(表1、表2)。

然而,根據(jù)下列關(guān)系式,多場27Al MAS NMR能夠?yàn)樘崛≥^為準(zhǔn)確的NMR參數(shù)提供限制條件[50-51]：

(3)

(4)

(5)

(6)

因此,通過測量不同場中27Al MAS NMR譜峰中心 (δmcg),然后對(duì)場強(qiáng)因子Cm/VL2投點(diǎn),使用公式(5)就能夠獲得參數(shù)SOQE和δiso。由于高嶺土結(jié)構(gòu)中兩個(gè)Al位置的相似性,使用4.68、7.1、9.4、11.7、14.04和21.06 T場中CT平均δ1/2cg值與C1/2/VL2投點(diǎn)(圖4),投點(diǎn)圖的最小二乘法分別給出斜率和截距值-8.84和6.25,最小二乘法的相關(guān)系數(shù)高達(dá)0.998(圖4)。由于斜率和截距分別對(duì)應(yīng)著-SOQE2和δiso值(公式(5)),所以這兩個(gè)Al位置的平均Cq值2.89 MHz,與FC LAPW理論計(jì)算值2.65 MHz十分吻合(表2),而截距6.25×10-6(δiso)也與實(shí)驗(yàn)化學(xué)位移和理論計(jì)算結(jié)果很接近(表1、表2),并與21 T的CT譜峰CG值幾乎完全一致。這進(jìn)一步證實(shí)在超高場中,高嶺土QI作用降至次要地位、四極化學(xué)位移對(duì)各向同性化學(xué)位移的移動(dòng)幾乎可以忽略。

圖4 4.68、7.1、9.4、14.04和21.06場CT NMR波譜的CG (10-6)與C1/2/VL2投點(diǎn)圖Fig.4 Centers of gravity (CG in ×10-6)for the NMR central transitions at 4.68,7.1,9.4,14.04 and 21.06 T magnetic fields plotted against C1/2/VL2

SATRAS旋轉(zhuǎn)邊帶總體外形(圖2)也可為提取NMR參數(shù)提供進(jìn)一步的限制,從而得出更為準(zhǔn)確的結(jié)果,這是因?yàn)樾l(wèi)星躍遷主要受QI影響,Cq僅與SATRAS邊帶覆蓋范圍同步變化,因此可以通過SATRAS 所有邊帶的覆蓋頻率范圍來估算Cq(≈ 3Cq/2I)。這種估算方法對(duì)較小Cq如高嶺土中Al位置更為重要,因?yàn)檩^小Cq往往不能給出可供擬合的典型QI線型。 通過類比已經(jīng)準(zhǔn)確測定27AlCq的Al-F礦物如cryolithionite、cryolite的SATRAS邊帶所覆蓋的頻率范圍,高嶺土中27AlCq被限制在2.83～3.30 MHz間;另一方面,所有邊帶總體外形對(duì)η很敏感,據(jù)此高嶺土中27Alη被限制在 0.7～0.9間,均與理論計(jì)算十分吻合。

因此,雖然NMR波譜沒能直接分辨出高嶺土中兩個(gè)不同的Al晶格位置,但通過理論計(jì)算和多種NMR實(shí)驗(yàn)方法,高嶺土中 Al 位置的NMR參數(shù)被準(zhǔn)確測算,這種方法可應(yīng)用推廣至其他固體精細(xì)局域結(jié)構(gòu)的研究之中。值得一提的是這項(xiàng)研究表明,對(duì)固體、尤其是礦物NMR而言,并非場強(qiáng)越高,精度和分辨率也隨之越高,必須在分析固體中各種NMR作用如QI、CSA,根據(jù)這些作用對(duì)場強(qiáng)的響應(yīng)特性,選取最優(yōu)場強(qiáng)以獲得最高分辨率和精度的NMR波譜。

[1] Smith M.E..Application of27Al NMR techniques to structure determination in solids [J].Appl.Magn.Reson.,1993,4,1-64.

[2] Ashenhurst J.,Gang Wu,G.,Wang S..Syntheses,Structures,Solution,and Solid-State27Al NMR Studies of Blue Luminescent Mononuclear Aluminum Complexes:Al(7-azain)2(7-azain-H)(CH3),Al(7-azain)3(7-azain-H),and Al(7-azain)(7-azain-H)(OCH(CF3)2)2(7-azain-H = 7-azaindole) [J].J.Am.Chem.Soc.,2000,122:2541-2547,

[3] Liu Y.,Tossell J.,Nekvasil H..A theoretical study on structural factors correlated with23Na NMR parameters [J].Amer.Mineral.,2004,89:1314-1322.

[4] Kirkpatrick R.J.,Phillips B.L..27Al NMR spectroscopy of minerals and related materials [J].Appl.Magn.Reson.,1993,4:213-236.

[5] Slichter C.P..Principles of Magnetic Resonance [M].Belin:Springer-Verlag,1992.

[6] Rocha J.,Klinowski J..The rehydration of metakaolinite to kaolinite [J]. J.Chem.Soc.Chem.Commun.,1991,3:582-584.

[7] MacKenzie K.J.D.,Meinhold R.H.,White G.V.,Sheppard C.M.,Sherriff B.L..Carbothermal formation of β-sialon from kaolinite and halloysite studied by29Si and27Al solid state NMR [J].J.Mater.Sci,1994,29:2611-2619.

[8] Rocha J..Single and triple quantum27Al MAS NMR study of the thermal transformation of kaolinite [J].J.Phys.Chem.B.,1999,103:9801-9804.

[9] Ashbrook S.E.,McManus J.,Mackenzie K.J.D.,Wimperis,S..Multiple-quantum and cross-polarized27Al MAS NMR of mechanically treated mixtures of kaolinite and gibbsite [J].J.Phys.Chem,2000,104:6408-6414.

[10] Massiot D.S.,Dion P.,Alcover J.F.,Bergaya F..27Al and29Si MAS NMR study of kaolinite thermal decomposition by controlled rate thermal analysis [J].J.Amer.Ceramic Soc.,1995,78:2940-2944.

[11] Meinhold R.H.,Mackenzie K.J.D.,Brown I.W.M..Thermal reactions of kaolinite studied by solid state 27-Al and 29-Si NMR [J]. J.Mater.Sci,1985,4:163-166.

[12] Watanabe T.,Shimizu H.,Nagasawa K.,Masuda A.,Saito H..29Si-and27Al-MAS/NMR study of the thermal transformations of kaolinite [J].Clay Minerals,1987,22:37-48.

[13] Sanz J.,Madani A.,Serratosa J.M.,Moya J.S.,Aza S..Aluminum-27 and Silicon-29 Magic-Angle Spinning Nuclear Magnetic Resonance study of the kaolinite-Mullite transformation [J].J.Amer.Ceramic Soc.,1988,71:C418-C421.

[14] Lambert J.F.,Millman W.S.,Fripiat J.J..Revisiting kaolinite dehydroxylation:A29Si and27Al MAS NMR study [J].J.Am.Chem.Soc.,1989,111:3517-3522.

[15] Wu Y.,Chmelka,B.F.,Pines A.,Davis M.E.,Grobet P.J.,Jacobs P.A..High-resolution27Al NMR spectroscopy of the aluminophosphate molecular sieve VPI-5 [J].Nature.,1990a,346:550-552.

[16] Wu Y.,Sun B.Q.,Pines A.,Samoson A.,Lippmaa,E..NMR experiments with a new double rotor [J].J.Magn.Reson.,1990b,4:203-211.

[17] Müeller K.T.,Sun B.Q.,Chingas C.G.,Zwanziger J.W.,Terao T.,Pines A..Dynamic-angle spinning of quadrupolar nuclei [J].J.Magn.Reson.,1990,84:470-487.

[18] Müeller K.T.,Wooten E.W.,Pines A..Pure-absorption-phase dynamic-angle spinning [J].J.Magn.Reson.,1991,92:620-627.

[19] Xu Z.,Sherriff B.L..27Al double-rotation NMR study of Al2SiO5polymorph minerals [J].Appl.Magn.Reson.,1993,4:203-211.

[20] Xu Z.,Sherriff B.L..23Na27Al9Be29Si solid state NMR study of tugtupite [J].Can.Mineral,1994,32:935-943.

[21] Kohn W.,Sham L.J..Self-consistent equation including exchange and correlation effects [J].Phys.Rev.,1965,140(4A):1133-1138.

[22] Blaha P.,Schwarz K.,Herzig P..First-principles calculation of the electric filed gradient of Li3N [J].Phys.Rev.Lett.,1985,54:1192-1195.

[23] Blaha P.,Schwarz K.,Dederichs P.H..First-principles calculation of the electric-field gradient in hcp metals [J].Phys.Rev.,1988,B37:2792-2796.

[24] Blaha P.,Schwarz K.,Sorantin P..Full potential linearized augmented plane wave programs for crystalline systems [J].Comp.Phys.Comm.,1990,59:39-415.

[25] Blaha P.,Singh D.J.,Sorantin P.I.,Schwarz K..Electric-field-gradient calculations for systems with large extended-core-state contributions [J].Phys.Rev.,1992,B46:1321-1325.

[26] Blaha P.,Schwarz K.,Madsen G.K.,Kvasnicka D.,Luitz J..WIEN2k,an augmented plane wave+local orbitals program for calculating crystal properties[M].Wien:Karlheinz Schwarz,Techn.Universit?t,2001.

[27] Bryant P.L.,Harwell C.R.,Wu K.,Fronczek F.R..Single-Crystal27Al NMR of Andalusite and Calculated Electric Field Gradients:the First Complete NMR Assignment for a 5-Coordinate Aluminum Site [J].J.Phys.Chem.,1999,103:5246-5252.

[28] Zhou B.,Giavani T.,Bildsoe H.,Skibsted J.,Jakobsen H..Structural refinement of CsNO3(II) by coupling of14N MAS NMR experiments with WIEN2k DFT calculations [J].Chem.Phys.Lett.,2005,402:133-137.

[29] Zhou B.,Sherriff B.,Hartman J.S.,Wu W..27Al and23Na NMR spectroscopy and Structural modeling of aluminofluoride minerals [J].Am.Mineral,2007,92:34-43.

[30] Profeta M.,Mauri F.,Pickard C.J..Accurate First Principles Prediction of17O NMR Parameters in SiO2:Assignment of the Zeolite Ferrierite Spectrum [J].J.Am.Chem.Soc.,2003,125:541.

[31] Clark S.J.,Segall M.D.,Pickard C.J.,Hasnip P.J.,Probert M.J.,Refson K.,Payne M.C..First principles methods using CASTEP [J].Zeitschrift fuer Kristallographie,2005,220:567-570.

[32] Pickard C.J.,Mauri F..All-electron magnetic response with pseudopotentials:NMR chemical shifts [J].Phys.Rev.,2001,245(B63):101.

[33] Yates J.R.,Pickard C.J.,Mauri F..Calculation of NMR Chemical Shifts for extended systems using Ultrasoft Pseudopotentials [J].Phys.Rev.,2007,24(B76):401.

[34] Foresman J.B.,Frisch A..Exploring chemistry with electronic structure methods,2ndEdition [M].Pittsburgh:Gaussian,Inc.,1996.

[35] Richards W.G.,Horsley J.A..AB Initio Molecular Orbital Calculations for Chemists [M].Oxford:Clarendon Press,1970.

[36] Collins D.R.,Catlow C.R..Energy-minimized hydrogen-atom positions of kaolinte [J].Acta Cryst.,1991,B47:678-682.

[37] Bish D.L..Rietvveld refinement of the kaolinite structure at 1.5 K [J].Clays and Clay Minerals,1993,41:738-744.

[38] Wu G.,Rovnyak D.,Sun B.,Griffin G..High-resolution multiple quantum MAS NMR spectroscopy of half-integer quadrupolar nuclei [J].Chem.Phys.Lett.,1995,249:210-217.

[39] Wu G.,Rovnyak,Griffin G..Quantitative multiple-quantum magic-angle-spinning NMR spectroscopy of quadrupolar nuclei in solids [J].J.Am.Soc.,1996,118:9326-9332.

[40] Wu G.,Rovnyak D.,Huang P.C.,Griffin G..High-resolution oxygen-17 spectroscopy multiple-quantum magic-angle-spinning [J].Chem.Phys.Lett.,1997,277:79-83.

[41] Amoureux J.P.,Fernandez C.,Steuernagel S..Z Filtering in MQMAS NMR [J].J.Magn.Reson.,1996,A123:116-118.

[42] Amoureux J.P,Pruski M.,Lang D.P.,Fernandez C..The effect of RF power and spinning speed on MQMAS [J].J.Magn.Reson.,1998,131:170-175.

[43] Amoureux J.P,Huguenard C.,Engelke F.,Taulelle F..:Unified representation of MQMAS and STMAS NMR of half-integer quadrupolar nuclei [J].Chem.Phys.Lett.,2002,356:497-504.

[44] Rocha J.,Pedrosa de Jesus J.D..27Al satellite transition MAS-NMR spectroscopy of kaolinite [J].Clay Minerals,1994,29:287-291.

[45] Newman R.H.,Childs C.W.,Churchman G.J..Aluminium coordination and structural disorder in halloysite and kaolinite by27Al NMR spectroscopy [J].Clay Minerals,1994,29:305-312.

[46] Zhou B.,Sherriff B.L..Quantum calculations of the electronic structure and NMR quadrupolar interaction parameters for tugtupite [J]. Am.Mineral,2004,89:377-381.

[47] J?ger C..How to get more from27Al MAS NMR by high-speed satellite-transition spectroscopy [J].J.Magn.Reson.,1992,99:353-362.

[48] Lacassagne V.,Florian P.,Montouillout V.,Gervais C.,Babonneau F.,Massiot D..Resolution enhancement in solid-state MQ-MAS experiments achieved by composite decoupling [J].Magn.Reson.Chem.,1998,36:956-959.

[49] Zhou B.,Sherriff B.L.,Taulelle F.,Wu G..Nuclear Magnetic Resonance Study of Al:Si and F:OH order in Zunyite [J].Can.Mineral,2003,41:891-903.

[50] Samoson A..Satellite transition high-resolution NMR of quadrupolar nuclei in powders [J].Chem.Phys.Lett.,1985,119:29-32.

[51] Skibsted J.,Nielsen N.,Bilds?e H.,Jakobsen H.,Satellite transitions in MAS spectra of quadrupolar nuclei [J].J.Magn.Reson.,1991,95:117-132.

(責(zé)任編輯：陳文寶)

Calculation study of27Al NMR Spectrum in Muti-magietic Field with KaoliniteStructure and ab initio Quantum Theory

SU Shaoming, ZOU zhen

(CentralSouthInstituteofGeologicalExploration,ChinaGeneralBureauofMetallurgyGeology,Wuhan,Hubei430081)

Many unsuccessful attempts have been made to resolve the two different Al sites in kaolinite by magic angle spinning nuclear magnetic resonance (MAS NMR) spectroscopy since the early 1980’s. This study attempted to resolve these sites and explain the reasons for the lack of separation by combining NMR spectroscopy at various fields with ab initio quantum modeling. Ata high field such as 21.06 T, large chemical shift anisotropy (CSA) and dipolar effects (DD) can overwhelm the quadrupolar (QI) effects for sites with relatively smallCqvalues such as the two Al sites in kaolinite, so the spectra display featureless line-shapes lacking the typical QI singularities and shoulders. At lower fields such as 4.68 T, the QI effects become dominant but due to the similarity of the two Al environments in kaolinite, the27Al MAS NMR peaks could not be separated either. Theoretical ab initio quantum modeling using Full Potential Linearized Augmented Plane Wave (FP LAPW) and Molecular Orbital (MO) allowed experimental27Al MAS NMR, multiple quantum (MQ) MAS and satellite transition (SATRAS) spectra to be interpreted.Cq,η, CSA andδisofor the two Al sites were obtained with high accuracy by simulating these spectra, while constraining the QI parameters.

Kaolinite;27Al MASNMR; MQMAS; SATRAS; Calculation model of FP LAPW; CASTEP; Calculation model of molecular orbital

2015-05-07;改回日期:2016-03-21

本文得到湖北省國土資源廳(編號(hào)為ETZ2015A05)科研項(xiàng)目的資助。

蘇紹明(1965-),男,高級(jí)工程師,礦產(chǎn)普查與勘探專業(yè),從事地質(zhì)找礦工作。E-mail:78737567@qq.com

O641

A

1671-1211(2016)02-0211-06

10.16536/j.cnki.issn.1671-1211.2016.02.017

數(shù)字出版網(wǎng)址:http://www.cnki.net/kcms/detail/42.1736.X.20160303.1642.038.html 數(shù)字出版日期:2016-03-03 16:42