InOOH壓致氫鍵對稱化及其彈性性質?

康端 巫翔

1)(北京大學地球與空間科學學院,造山帶與地殼演化教育部重點實驗室,北京 100871)

2)(中國地質大學(武漢),地質過程與礦產資源國家重點實驗室,武漢 430074)

InOOH壓致氫鍵對稱化及其彈性性質?

康端1)巫翔2)?

1)(北京大學地球與空間科學學院,造山帶與地殼演化教育部重點實驗室,北京 100871)

2)(中國地質大學(武漢),地質過程與礦產資源國家重點實驗室,武漢 430074)

InOOH,氫鍵對稱化,彈性性質

1 引 言

氫鍵廣泛存在于DNA、水、納米材料、地球深部含水物質中,對生命、科技、地球的演化等都具有至關重要的作用,在生物、物理、材料、地學等領域得到了廣泛的研究.氫鍵在高壓或高溫下易發生畸變,且高壓對氫鍵的影響十分顯著.當壓力施加于含有氫鍵的物質上時,氫鍵的變化會使物質表現出一些新的性質,甚至引起相變的發生.高壓下氫鍵對稱化轉變是一種典型的行為,是指與H原子最鄰近的兩個O原子之間的勢能面從雙勢阱勢能面轉變為單勢阱勢能面的過程,即H原子從偏離兩個O原子中間的位置變為兩個O原子正中間的位置.自從Holzapfel[1]在理論上預測出氫鍵對稱化之后,壓致氫鍵對稱化行為就在冰、甲酸、MOOH(M=In,Al,Ga,Fe,Cr)等[2?7]物質中得到了廣泛的實驗和理論研究.例如,實驗研究結果表明冰VII相會在約62 GPa壓力下由于氫鍵的對稱化變為冰X相[8].

InOOH,δ-AlOOH,β-GaOOH,ε-FeOOH,β-CrOOH羥氧化物具有相似構型[2,4,9?11],本文統一稱之為MOOH(M=In,Al,Ga,Fe,Cr).MOOH中O原子的排列具有畸變的金紅石型結構,O原子近乎完美地平行于(101)面呈緊密排列,共邊的MO6八面體沿(001)面呈鏈狀分布,這些鏈之間通過O原子共角頂相連,如此便形成了平行于(101)面的八面體層.大量的實驗和理論研究結果表明MOOH會在高壓下發生氫鍵對稱化,且氫鍵的對稱化對MOOH的性質有重要影響.實驗結果表明β-GaOOH,InOOH,β-CrOOD及δ-AlOOH的軸比率對壓力的斜率正負性在具有非對稱氫鍵的相和具有對稱氫鍵的相中不同,表明物質在不同方向上的可壓縮性相對大小在氫鍵非對稱相和氫鍵對稱相中有差異[2,12,13].實驗和第一性原理計算結果表明氫鍵對稱化使δ-AlOOH的體積模量、體波速以及縱波波速等參數異常增大[14,15].另外,實驗和理論研究表明ε-FeOOH中的氫鍵對稱化還對Fe的自旋轉變有促進作用[4].

InOOH為帶隙為3.6 eV的n型半導體材料,是一種能在紫外光下有效分解有機污染物的光催化劑[16?18],吸引了大量學者對InOOH合成方法進行實驗研究[16,19,20].文獻[16]研究結果表明可利用溶劑熱合成法加工In(NO3)3和乙二胺合成InOOH.有研究人員曾通過實驗和理論計算研究了畸變金紅石型InOOH相變為黃鐵礦型InOOH的溫壓條件,第一性原理的熱力學計算結果表明畸變金紅石型InOOH在0 K,15 GPa時會相變為氫鍵非對稱的黃鐵礦型InOOH[3],但X射線衍射(XRD)測試結果表明常溫下畸變金紅石型InOOH加壓至35 GPa時都未相變為黃鐵礦型InOOH[2],當加溫至1300 K時需在14 GPa的壓力下才能轉變為黃鐵礦型InOOH[21].

研究InOOH彈性性質在高壓下的變化及由高壓導致的氫鍵對稱化對InOOH彈性性質的影響,對探索InOOH在材料領域的應用及預測氫鍵對稱化對其他畸變金紅石型MOOH性質的影響具有重要的作用,但目前關于InOOH在高壓下的彈性性質,及氫鍵對稱化對InOOH彈性性質影響的詳細研究鮮有報道.第一性原理計算是研究晶體結構、力學性質、熱學性質的有效方法,應用越來越廣泛[22?28].本文基于第一性原理方法開展了對高壓下InOOH的彈性性質、氫鍵對稱化及氫鍵對稱化對彈性性質影響的研究,計算結果表明氫鍵對稱化的壓力約為18 GPa,得到了InOOH在0—40 GPa范圍內的彈性性質,以及氫鍵對稱化對彈性性質的影響.

2 計算模型和方法

常溫常壓下InOOH為正交結構,空間群為P21nm.晶胞參數為a=5.26 ? (1 ? =0.1 nm),b=4.56 ?,c=3.27 ?,O—H鍵長為1.08 ?,O···H鍵長為1.46 ?,氫鍵H—O···H呈非對稱構型(圖1(a))[29].參照δ-AlOOH,ε-FeOOH等相似物的高壓行為,推測InOOH中的氫鍵在高壓下呈對稱構型,對應的空間群為Pnnm(圖1(b)).以InOOH單胞為研究對象,每個單胞中含有2個In原子,4個O原子和2個H原子.為了表述方便,將具有非對稱氫鍵的MOOH晶體標記為A-MOOH,將具有對稱氫鍵的MOOH晶體標記為S-MOOH.

第一性原理計算采用VASP(Viennaab-initioSimulation Package)軟件包完成[30,31].贗勢采用投影綴加平面波(projected augmented wavefunction,PAW)描述.電子和電子之間的交換關聯采用廣義梯度近似(generalized gradient approximation,GGA)下的Perdew-Burke-Ernzerhof(PBE)泛函進行處理[32].平面波的截斷能量為850 eV.布里淵區求和采用6×6×10的Monkhorst-Pack型K點網格.更大的平面波截斷能和K空間采樣密度對每個原子總能的影響不超過0.001 eV.電子自洽的收斂條件為系統總能量變化小于1×10?7eV,結構弛豫的收斂標準為每個離子上的最大作用力小于0.001 eV/?.利用應力-應變法,通過對晶格施加6次應變且允許離子弛豫計算得到彈性常數[33,34].

圖1 (網刊彩色)(a)具有非對稱氫鍵的A-InOOH晶體結構,空間群為P21nm;(b)具有對稱氫鍵的S-InOOH晶體結構,空間群為PnnmFig.1. (color online)(a)Crystal structure of AInOOH with asymmetric hydrogen bond,space group P21nm;(b)crystal structure of S-InOOH with symmetric hydrogen bond,space group Pnnm.

3 計算結果

基于GGA-PBE交換關聯泛函計算得到InOOH,在0 GPa時的晶胞參數為a=5.385 ?,b=4.662 ?,c=3.352 ?, 與文獻[3]采用GGAPW91交換關聯泛函得到的晶胞參數(a=5.381 ?,b=4.649 ?,c=3.349 ?)基本一致,但稍大于實驗得到的晶胞參數(a=5.26 ?,b=4.56 ?,c=3.27 ?).這是GGA泛函會低估晶體的結合能,因而高估晶體的晶胞參數導致的.

0 GPa時InOOH的O···H鍵 長 為1.519 ?,O—H鍵長為1.046 ?.隨著壓力的增加,O···H鍵縮短,O—H鍵變長(圖2(a)).18 GPa時O···H鍵長為1.203 ?,O—H鍵長為1.186 ?,兩者相差0.017 ?,氫鍵O—H···O基本對稱. 文獻[3]第一性原理計算結果表明InOOH中的氫鍵約在25 GPa時對稱,本文計算得到的氫鍵對稱化壓力與文獻[3]結果基本一致.a/c對壓力的斜率在0—40 GPa內為正且變化不大(圖2(b)).b/c對壓力的斜率在18 GPa之前為負值,在18 GPa之后為正值(圖2(b)),這表明在18 GPa時InOOHb軸和c軸壓縮性的相對大小發生了變化.實驗結果表明InOOH中的a/c對壓力斜率的正負性在0—35 GPa內沒有發生變化,b/c對壓力的斜率約在15 GPa時由負值變為正值.本文計算結果與實驗結果基本一致,但b/c對壓力斜率正負性發生變化的壓力稍小于XRD實驗得到的壓力(15 GPa)[2],這是由于未考慮熱效應和核量子效應導致的.

圖2 (網刊彩色)(a)InOOH中O—H鍵和O···H鍵鍵長隨壓力的變化;(b)InOOH軸比率a/c和b/c隨壓力的變化Fig.2.(color online)(a)Pressure dependence of O—H and O···H bond lengths in InOOH;(b)pressure dependence of axial ratios a/c and b/c in InOOH.

A-InOOH和S-InOOH都屬于正交晶系,有9個獨立的彈性常數:壓縮彈性常數C11,C22,C33,剪切彈性常數C44,C55,C66和非對角彈性常數C12,C13,C23.由計算結果可知,InOOH在0—40 GPa內的彈性常數均滿足Born穩定性判據[35],說明0 K時InOOH在0—40 GPa內力學性能穩定.

InOOH的壓縮彈性常數和非對角彈性常數在0—40 GPa內都隨著壓力的增加而增大,在15 GPa之前對壓力的斜率與在20 GPa之后對壓力的斜率較為接近,但在15—20 GPa之間對壓力的斜率顯著大于15 GPa之前和20 GPa之后的斜率(圖3(a)).這表明壓縮彈性常數和非對角彈性常數在15—20 GPa內發生了突變.C66隨著壓力的增加而增加,其對壓力的斜率在0—40 GPa內變化不大.C55在15 GPa之前隨著壓力的增加稍有增大,在20 GPa之后隨著壓力的增加稍有減小,但在0—40 GPa內的變化幅度不大.C44在0—40 GPa內基本都隨著壓力的增加而減小,且在20 GPa之后的減小速率比在15 GPa之前大.C44代表晶體抵抗剪切變形的能力,C44隨壓力變小表明壓力增加時InOOH抵抗剪切變形的能力下降.

根據Voigt-Reuss-Hill近似方法[36]可以計算晶體的體積模量B和剪切模量G,進而可以算出楊氏模量E和泊松比ν[37].壓力小于15 GPa或大于20 GPa時,InOOH的體積模量B隨壓力的增加幾乎都呈線性增大(圖3(b)).這是由于隨著壓力的增加,晶胞體積縮小,原子間的相互排斥力增大,InOOH更加難以被壓縮.壓縮彈性常數和非對角彈性常數在15—20 GPa內的異常增加,引起了體積模量在15—20 GPa內的突變.15 GPa前InOOH的剪切模量隨著壓力的增加而增加,但20 GPa后剪切模量隨著壓力的增加而減小.這表明15 GPa前InOOH隨著壓力的增加抵抗剪切變形的能力增強,20 GPa后InOOH隨著壓力的增加抵抗剪切變形的能力變弱.

根據體積模量B和剪切模量G,可通過公式計算得到縱波波速VP和橫波波速VS,其中ρ為密度.InOOH的縱波波速VP隨著壓力的增加而增加,且在15—20 GPa之間異常增加;橫波波速VS在15 GPa之前隨著壓力的增加緩慢減小,在20 GPa之后隨著壓力的增加而減小,且減小的速率顯著大于15 GPa之前的速率(圖3(c)).雖然InOOH的剪切模量在15 GPa前隨著壓力的增加而增加,但由于剪切模量隨壓力增加的速率很小,小于密度隨壓力的增加速率,導致G/ρ隨壓力的增加而減小,因此剪切波速在15 GPa之前隨著壓力的增加而減小.

圖3 (網刊彩色)InOOH在不同壓力下的(a)彈性常數,(b)體積模量B和剪切模量G,(c)縱波波速VP和橫波波速VSFig.3.(color online)Pressure dependence of(a)elastic constants,(b)bulk modulus B and shear modulus G,(c)longitudinal wave velocity VPand shear wave velocity VSin InOOH.

圖4 InOOH在不同壓力下的B/G和泊松比νFig.4.B/G and Poisson ratio ν of InOOH under different pressures.

對于固體物質,如果其泊松比ν大于0.26,則表現為韌性;反之則表現為脆性.根據Pugh準則[38],B/G＞1.75時物質表現為韌性;B/G＜1.75時物質表現為脆性.由圖4可見,InOOH在0—40 GPa內的泊松比ν大于0.26,B/G值大于1.75,表明InOOH在0—40 GPa內都表現為韌性.隨著壓力的增加,泊松比ν和B/G值增大,InOOH的韌性增強,且在15—20 GPa內發生了韌性的異常增加.

楊氏模量E在15 GPa前隨著壓力的增加而增加,在20 GPa后隨著壓力的增加而減小,但在0—40 GPa內變化幅度不大.楊氏模量E在各個方向上的三維表示可以直觀地用于分析晶體的彈性各向異性.對于具有正交結構的晶體,楊氏模量在任意方向的值可根據如下公式計算[39]:

式中Sij為彈性柔順常數,l1,l2,l3為方向余弦.圖5為不同壓力下InOOH在各個方向上的楊氏模量.對于各向同性晶體,楊氏模量表面為完美的球形.在0 GPa時,InOOH的楊氏模量表面不為球形,有沿著 [001],[00ˉ1],[110],[ˉ110],[ˉ1ˉ10]和 [1ˉ10]6 個方向的突起,這表明InOOH呈現各向異性.[100],[010]和[001]方向的楊氏模量在0—15 GPa內隨著壓力的增加而增大,但在20—40 GPa內隨著壓力的增加而減小,其中[001]方向的楊氏模量隨壓力的增加變化最小.[110],[ˉ110],[ˉ1ˉ10]及 [1ˉ10]方向的楊氏模量在0—40 GPa內都隨著壓力的增加而增大.由圖5可見楊氏模量表面偏離球形的程度越來越大,表明InOOH的各向異性程度越來越深.

圖5 (網刊彩色)InOOH在不同壓力下的楊氏模量三維立體圖Fig.5.(color online)The variation of 3D Young’s modulus surface with pressure.

4 討 論

以上計算結果表明InOOH的氫鍵約在18 GPa時發生了對稱化,且軸比率b/c對壓力的導數、彈性性質、波速等在15—20 GPa間發生了突變,兩種壓力呈現出很好的一致性.表明InOOH中氫鍵的對稱化導致了15—20 GPa間部分性質的突變.

本文和前人的研究結果表明,具有相似晶體結構的InOOH,δ-AlOOH和β-GaOOH中的氫鍵對稱化均導致b/c對壓力的導數由負變為正[2,13],這說明氫鍵對稱化對MOOH的b/c對壓力的斜率有重要影響.MOOH中的氫鍵O—H···O 位于ab平面內,且在b軸上的投影最大,因此氫鍵強度的變化對b軸的影響最為顯著.對稱氫鍵的強度比非對稱氫鍵大,在氫鍵由非對稱構型變為對稱構型時b軸的可壓縮性降低,從而導致A-MOOH晶體中b軸比c軸易壓縮,但在S-MOOH晶體中b軸比c軸難壓縮的現象.在A-InOOH中a軸本身就難以壓縮,表現為a/c隨著壓力的增加而增大(圖2(b)),且氫鍵對稱化對a軸可壓縮性的影響較小,因此氫鍵對稱化發生時a/c對壓力的斜率變化不大.

本文InOOH和文獻[3]中δ-AlOOH的第一性原理計算結果表明氫鍵的對稱化會導致InOOH和δ-AlOOH壓縮和非對角彈性常數的異常增加,但不會引起剪切彈性常數的異常增加[15].據此可推測MOOH中氫鍵的對稱化對彈性常數的影響都與此類似.由于MOOH中的氫鍵對稱化對b,a,c軸的影響依次減小,因此引起的C22,C11,C33增加幅度依次減小,這在InOOH及δ-AlOOH[15]的計算結果中得到了證實.Phase D(MgSi2O6H2)是一種穩定壓力最高的致密含水含鎂硅酸鹽[40],能穩定至40 GPa.由于Phase D中的氫鍵平行于c軸,因此當Phase D中的氫鍵發生對稱化時,C33的增加幅度最大[41].由此可知,在含有氫鍵的物質中,氫鍵對稱化對各個方向上壓縮彈性常數的影響取決于氫鍵在各軸上的投影,投影越大,影響越大.

對15 GPa前A-InOOH的體積模量BA與壓力P采用線性擬合得到的關系為BA(GPa)≈4.126×P(GPa)+126(GPa);20 GPa后S-InOOH的BS與P的關系為BS(GPa)≈3.101×P(GPa)+177(GPa).對比可知InOOH的BS隨壓力的變化速率比BA的小.這與δ-AlOOH的計算結果一致[5],可推測MOOH的體積模量BS隨壓力的變化速率都小于BA的變化速率.A-InOOH與SInOOH體積模量之間的差值為BS?BA(GPa)≈?1.025×P(GPa)+51(GPa),20 GPa之后InOOH只能以S-InOOH的形式存在,因此本文僅討論20 GPa之前BS與BA的差異.0 GPa時BS約比BA大51 GPa,約為BA的40%;BS與BA之間的差異隨著壓力的增加而減小,15 GPa時BS比BA大36 GPa,約為BA的19%.以上數據表明氫鍵對稱化是圖3(b)中InOOH的體積模量在15—20 GPa內發生異常增加的原因.XRD實驗結果表明當體積模量對壓力的一階導數B0′固定為4時A-InOOH的零壓體積模量為159 GPa[2],比S-InOOH的零壓體積模量(192 GPa)約小33 GPa,約為A-InOOH零壓體積模量的21%.因此由本文的第一性原理計算及XRD實驗結果可知,InOOH中氫鍵的對稱化會導致體積模量20%—40%的異常增加.

第一性原理計算結果表明S-δ-AlOOH的零壓體積模量為207 GPa,比A-δ-AlOOH的零壓體積模量(153 GPa)大54 GPa,增加量約為35%[15].實驗結果表明氫鍵對稱化導致δ-AlOOH的零壓體積模量由183 GPa變為226 GPa,約增加23%[42].第一性原理計算結果顯示氫鍵的對稱化導致Phase D的零壓體積模量由159 GPa變為191 GPa,約增加20%[41].InOOH,δ-AlOOH和Phase D的研究結果均表明氫鍵對稱化會引起物質體積模量20%—40%的增加,這是由氫鍵對稱化時壓縮彈性常數的異常增加導致的.氫鍵的對稱化不會導致剪切模量的異常增加.

實驗和第一性原理計算結果均表明InOOH的體積模量小于δ-AlOOH的體積模量.In3+的半徑(0.81 ?)顯著大于Al3+的半徑(0.50 ?),In3+與O2?之間的相互作用弱于Al3+與O2?之間的相互作用.In—O鍵的強度比Al—O鍵弱,In—O鍵的鍵能為(346±30)kJ/mol,Al—O鍵的鍵能為(502±11)kJ/mol[43],這是InOOH的體積模量小于δ-AlOOH的體積模量的原因之一.δ-AlOOH中Al—O鍵的平均鍵長為1.954 ?,AlO6八面體體積為9.850 ?3[15];InOOH中In—O鍵的平均鍵長為2.230 ?,InO6八面體體積為14.622 ?3,比δ-AlOOH中AlO6的體積大4.772 ?3.每個MOOH單胞中含有2個MO6八面體,因此InOOH和δ-AlOOH晶體中金屬八面體的體積差異為9.544 ?3.δ-AlOOH的晶胞參數為a=4.788 ?,b=4.275 ?,c=2.877 ?,體積為58.888 ?3[15];InOOH的晶胞參數為a=5.385 ?,b=4.662 ?,c=3.352 ?,體積為84.151 ?3,比δ-AlOOH的體積大25.263 ?3,遠大于MO6八面體的體積差異(9.544 ?3).這表明與δ-AlOOH相比,InOOH晶體中未被MO6八面體占據的空隙更大,使得InOOH更容易被壓縮,這也是導致InOOH體積模量比δ-AlOOH小的因素之一.

5 結 論

基于第一性原理研究了畸變金紅石型InOOH的氫鍵在高壓下的對稱化行為,以及壓力和氫鍵對稱化對InOOH彈性性質的影響.計算結果表明InOOH的氫鍵約在18 GPa處對稱,氫鍵的對稱化使b/c對壓力的斜率由負值變為正值.壓縮彈性常數和非對角彈性常數、體積模量以及縱波波速隨著壓力的增加而增加,且伴隨著氫鍵的對稱化異常增加,例如伴隨著氫鍵對稱化的發生體積模量發生了20%—40%的異常增加.剪切模量和楊氏模量在A-InOOH中隨壓力的增加而增加,在S-InOOH中隨壓力的增加而減小.彈性常數C44和橫波波速VS在氫鍵對稱化前后的InOOH中都隨著壓力的增加而減小,但氫鍵對稱化使其減小速率增大.InOOH在常壓時呈現韌性,隨著壓力的增加韌性增強,且伴隨著氫鍵對稱化韌性異常增加.氫鍵的對稱化使[100],[010]和[001]方向的楊氏模量隨壓力的增加而增大的變化趨勢變為隨壓力的增加而減小.結合其他畸變金紅石型MOOH的研究結果分析可知,MOOH化合物的軸比率、彈性性質、波速等性質隨壓力變化的規律與InOOH相似.

[1]Holzapfel W B 1972J.Chem.Phys.56 712

[2]Sano-Furukawa A,Yagi T,Okada T,Gotou H,Kikegawa T 2012Phys.Chem.Miner.39 375

[3]Tsuchiya J,Tsuchiya T,Sano A,Ohtani E 2008J.Miner.Petrol.Sci.103 116

[4]Gleason A E,Quiroga C E,Suzuki A,Pentcheva R,Mao W L 2013Earth Planet.Sci.Lett.379 49

[5]Goncharov A F,Manaa M R,Zaug J M,Gee R H,Fried L E,Montgomery W B 2005Phys.Rev.Lett.94 065505

[6]Lu X Z,Zhang Y,Zhao P,Fang S J 2011J.Phys.Chem.B115 71

[7]Benoit M,Marx D,Parrinello M 1998Nature392 258

[8]Aoki K,Yamawaki H,Sakashita M,Fujihisa H 1996Phys.Rev.B54 15673

[9]Suzuki A,Ohtani E,Kamada T 2000Phys.Chem.Miner.27 689

[10]Bolotina N B,Molchanov V N,Dyuzheva T I,Lityagina L M,Bendeliani N A 2008Crystallogr.Rep.53 960

[11]Christensen A N,Hansen P,Lehmann M S 1976J.Solid State Chem.19 299

[12]Kuribayashi T,Sano-Furukawa A,Nagase T 2013Phys.Chem.Miner.41 303

[13]Sano-Furukawa A,Kagi H,Nagai T,Nakano S,Fukura S,Ushijima D,Iizuka R,Ohtani E,Yagi T 2009Am.Mineral.94 1255

[14]Mashino I,Murakami M,Ohtani E 2016J.Geophys.Res.Sol.Ea.121 595

[15]Tsuchiya J,Tsuchiya T 2009Phys.Earth.Planet.In.174 122

[16]Li Z,Xie Z,Zhang Y,Wu L,Wang X,Fu X 2007J.Phys.Chem.C111 18348

[17]Zhao H,Yin W,Zhao M,Song Y,Yang H 2013Appl.Catal.B130 178

[18]Chen X,Lin P Y,Shi X C,Zhang Z L,Xie Z P,Li Z H 2009Chin.J.Inorg.Chem.25 1917

[19]Tao X J,Zhao Y B,Sun L,Zhou S M 2015Mater.Chem.Phys.149 275

[20]Muruganandham M,Lee G J,Wu J J,Levchuk I,Sillanpaa M 2013Mater.Lett.98 86

[21]Sano A,Yagi T,Okada T,Gotou H,Ohtani E,Tsuchiya J,Kikegawa T 2008J.Miner.Petrol.Sci.103 152

[22]Wang X F,Ma J J,Jiao Z Y,Zhang X Z 2016Acta Phys.Sin.65 206201(in Chinese)[王雪飛,馬靜婕,焦照勇,張現周2016物理學報65 206201]

[23]Wang H Y,Hu Q K,Yang W P,Li X S 2016Acta Phys.Sin.65 077101(in Chinese)[王海燕,胡前庫,楊文朋,李旭升2016物理學報65 077101]

[24]Hao J,Zhou G G,Ma Y,Huang W Q,Zhang P,Lu G W 2016Acta Phys.Sin.65 113101(in Chinese)[郝娟,周廣剛,馬躍,黃文奇,張鵬,盧貴武 2016物理學報 65 113101]

[25]Pan X D,Wei Y,Cai H Z,Qi X H,Zheng X,Hu C Y,Zhang X X 2016Acta Phys.Sin.65 156201(in Chinese)[潘新東,魏燕,蔡宏中,祁小紅,鄭旭,胡昌義,張詡翔 2016物理學報65 156201]

[26]Yao B D,Hu G Q,Yu Z S,Zhang H F,Shi L Q,Shen H,Wang Y X 2016Acta Phys.Sin.65 026202(in Chinese)[姚寶殿,胡桂青,于治水,張慧芬,施立群,沈皓,王月霞2016物理學報65 026202]

[27]Wu R X,Liu D J,Yu Y,Yang T 2016Acta Phys.Sin.65 027101(in Chinese)[吳若熙,劉代俊,于洋,楊濤2016物理學報65 027101]

[28]Huang A,Lu Z P,Zhou M,Zhou X Y,Tao Y Q,Sun P,Zhang J T,Zhang T B 2017Acta Phys.Sin.66 016103(in Chinese)[黃鰲,盧志鵬,周夢,周曉云,陶應奇,孫鵬,張俊濤,張廷波2017物理學報66 016103]

[29]Lehmann M S,Larsen F K,Poulsen F R,Christensen A N,Rasmussen S E 1970Acta.Chem.Scand.24 1662

[30]Kresse G,Furthmüller J 1996Phys.Rev.B54 11169

[31]Kresse G,Furthmüller J 1996Comp.Mater.Sci.6 15

[32]Perdew J P,Burke K,Ernzerhof M 1996Phys.Rev.Lett.77 3865

[33]Page Y L,Saxe P 2002Phys.Rev.B64 104104

[34]Liu L L,Wu X Z,Wang R,Nie X F,He Y L,Zou X 2017Crystals7 111

[35]Wu Z,Zhao E J,Xiang H P,Hao X F,Liu X J,Meng J 2007Phys.Rev.B76 054115

[36]Hill R 1952Proc.Phys.Soc.Sect.A65 349

[37]Wen Y F,Wang L,Liu H L,Song L 2017Crystals7 39[38]Pugh S F 1954The London,Edingburgh,and Dublin Philos.Mag.J.Sci.45 823

[39]Nye J F 1985Clarencon45 391

[40]Wu X,Wu Y,Lin J F,Liu J,Mao Z,Guo X,Yoshino T,McCammon C,Prakapenka V B,Xiao Y 2016J.Geophys.Res.Sol.Ea.121 6411

[41]Tsuchiya J,Tsuchiya T 2008Phys.Earth.Planet.In.170 215

[42]Tsuchiya J,Tsuchiya T,Tsuneyuki S,Yamanaka T 2002Geophys.Res.Lett.29 1909

[43]Luo Y R 2007Comprehensive Handbook of Chemical Bond Energies(Boca Raton:CRC Press)pp1057,1080

Pressure-induced hydrogen bond symmetrization of InOOH and its elastic properties?

Kang Duan1)Wu Xiang2)?

1)(Key Laboratory of Orogenic Belts and Crustal Evolution,Ministry of Education,School of Earth and Space Sciences,Peking University,Beijing 100871,China)
2)(State Key Laboratory of Geological Processes and Mineral Resources,China University of Geosciences,Wuhan 430074,China)

21 July 2017;revised manuscript

14 August 2017)

Pressure-induced hydrogen bond symmetrization of InOOH as well as its effects on the elastic properties is investigated by first-principles simulation.The results indicate that the hydrogen bond in InOOH symmetrized at about 18 GPa,resulting in the pressure derivative of theb/caxial ratio changing from negative to positive.While thea/caxial ratio increases with the increasing pressure over a range of 0—40 GPa,its pressure derivative does not change signi ficantly across the hydrogen bond symmetrization.

In the text,‘A-InOOH’denotes the asymmetric hydrogen bond phase and ‘S-InOOH’refers to the symmetric hydrogen bond phase.The compressional and off-diagonal elastic constants,bulk modulusB,Poisson’s ratioν,B/G(Grepresents shear modulus)and longitudinal wave velocityVPincrease with the increasing pressure in both A-InOOH and S-InOOH.These properties of A-InOOH are signi ficantly smaller than those of S-InOOH,and therefore they increase abnormally during the hydrogen bond symmetrization,such as a 20%—40%increase of the bulk modulus.Shear modulusGand Young’s modulusEincrease with the increasing pressure in A-InOOH,but decrease with the increasing pressure in S-InOOH,implying that hydrogen bond symmetrization would change their pressure evolution trends obviously.Shear elastic constantC44and shear wave velocityVSdecrease with the increasing pressure in both A-InOOH and S-InOOH,and more quickly in the latter,indicating that the structure change of hydrogen bond would change their pressure evolution rates.The Young’s moduli along the[100],[010]and[001]directions increase with the increasing pressure in A-InOOH,while decrease with the increasing pressure in S-InOOH,and those along the[110],[ˉ110],[ˉ1ˉ10]and[1ˉ10]directions always increase with the increasing pressure over a range of 0—40 GPa.The anisotropy and toughness of InOOH increase with the increasing pressure in both A-InOOH and S-InOOH,and the hydrogen bond symmetrization results in abnormal increase.In the materials containing hydrogen bonds,the effects of hydrogen bond symmetrization on different compressional elastic constants depend on the hydrogen bond projection on corresponding axes:the bigger the projection,the more significant the effect is.

InOOH has an obviously smaller bulk modulus than δ-AlOOH.The dominant reason is that the In3+radius(0.81 ?,1 ? =0.1 nm)is larger than Al3+radius(0.50 ?),resulting in the weaker interaction between In3+and O2?than that between Al3+and O2?.In addition,InOOH has more vacancies than δ-AlOOH.Combining with previous investigations on other rutile-distortedMOOH(M=Al,Ga,Fe,Cr),we can infer that the axial ratios,elastic properties and wave velocities of allMOOH materials have similar pressure evolutions to those of InOOH,and the hydrogen bond symmetrization has similar effects on the properties ofMOOH.

InOOH,hydrogen bond symmetrization,elastic properties

PACS:62.20.—x,82.30.Rs,62.20.de,62.20.dqDOI:10.7498/aps.66.236201

*Project supported by the National Natural Science Foundation of China(Grant No.41473056).

?Corresponding author.E-mail:wuxiang@cug.edu.cn

(2017年7月21日收到;2017年8月14日收到修改稿)

利用第一性原理研究了InOOH在高壓下的氫鍵對稱化行為及其對InOOH彈性等性質的影響.結果表明約在18 GPa時InOOH中的氫鍵發生了對稱化轉變,導致軸比率b/c對壓力的斜率由負值變為正值;壓縮彈性常數、非對角彈性常數、體積模量和縱波波速出現異常增加,如體積模量增加了20%—40%.高壓下InOOH彈性性質呈現出更加明顯的各向異性.常壓下InOOH呈現韌性,且伴隨著氫鍵對稱化韌性異常增加.對畸變金紅石型MOOH(M=Al,In,Ga,Fe,Cr)化合物在高壓下的彈性性質轉變與氫鍵性質轉變的耦合規律進行了初探.

10.7498/aps.66.236201

?國家自然科學基金(批準號:41473056)資助的課題.

?通信作者.E-mail:wuxiang@cug.edu.cn