Evolution of Crystal Structures and Electronic Properties for Ir2P under High Pressure
-
摘要: 在压强为0~100 GPa范围内,运用CALYPSO结构搜索技术,结合基于密度泛函理论中的第一性原理方法,对Ir2P晶体进行结构预测,并对预测出的晶体结构和物理性质进行细致的研究。在常压下,预测得出α-Ir2P相具有立方结构,其空间群为Fm3m,与实验所得结构一致;压强为86.4 GPa时,发生结构相变,由α-Ir2P相转变为β-Ir2P相,为四方结构,其空间群为I4/mmm。在相变过程中,晶体体积发生坍塌,并且出现不连续变化的一级相变。电子性质计算表明,86.4 GPa时,预测的β-Ir2P相中导带和价带在费米面附近发生交叠,表明其结构具有金属性质;电子局域函数计算表明,β-Ir2P相具有丰富的化学键,包括极性共价键、金属键和离子键;Bader电荷转移计算得出,由于Ir原子具有较强的电负性,β-Ir2P相中每个P原子向每个Ir原子电荷转移0.19e。Abstract: The crystals of Ir2P were predicted under the pressure ranging from 0 to 100 GPa using the CALYPSO structure exploration technique with the first-principles method based on the density functional theory. The predicted physical properties and crystal structures were examined in detail. At ambient pressure, the predicted α-Ir2P phase was found to have a cubic structure with Fm3m space group, which is consistent with the experimental structure. The pressure-induced structural transformations were unraveled, from the α-Ir2P phase to the β-Ir2P phase at 86.4 GPa. The predicted β-Ir2P phase has I4/mmm space group. In the process of phase transition, the volume of the crystal collapses and a discontinuous first order phase transition occurred. The calculation of the electronic properties showed that the predicted conduction bands and the valence bands of the β-Ir2P phase overlap near the Fermi surface at 86.4 GPa, indicating that the structure of the β-Ir2P phase has metallic properties. The electron localization function revealed that the β-Ir2P phase has a polar covalent bond, a metallic bond and an ionic bond. The Bader charge transfer calculations showed that each P atom transfers 0.19e to Ir atom, mainly due to the strong electronegativity of the Ir atoms.
-
Key words:
- high pressure /
- first-principles /
- crystal structure prediction /
- Ir2P
-
功能梯度材料(Functionally Graded Materials, FGM)的概念是1984年在航空飞机计划中首次提出的[1],FGM的特性在于它的组成和结构随着体积的变化而变化,从而导致材料相应性质发生改变。因其材料特性呈幂律分布[2-3],FGM被广泛应用于工程领域,如航空航天、机械工程、生物医学等。圆柱壳在联合荷载作用下的屈曲分析备受学术界关注[4-5]。目前,对FGM板壳的研究较为深入[6-8]。Beni等[9]利用改进的偶应力理论,对FGM圆柱壳在不同边界条件下的动力屈曲进行了分析;Kargarnovin等[10]研究了轴向荷载作用下FGM圆柱壳的动力屈曲;Sofiyev等[11]研究了横向压力下功能梯度正交各向异性圆柱壳的动力屈曲,推导出基于一阶剪切变形理论的功能梯度正交各向异性圆柱壳的稳定性和相容性方程;Khazaeinejad等[12]研究了弹性模量在厚度方向上连续变化的FGM圆柱壳在复合外压和轴向压缩载荷作用下的动力屈曲;Khalili等[13]研究了横向冲击载荷作用下FGM圆柱壳的动力屈曲;Alashti等[14]对变厚度FGM圆柱壳外压和轴向压缩的动力屈曲问题进行了分析。
基于以上研究,本研究讨论了FGM圆柱壳在轴向荷载作用下的动力屈曲。根据Donnell壳体理论和经典板壳理论,利用Hamilton变分原理得到FGM圆柱壳的动力屈曲控制方程;采用分离变量法求得动力屈曲临界荷载表达式;通过MATLAB软件计算动力屈曲临界荷载,讨论由不同材料(陶瓷和钛、陶瓷和铁、陶瓷和铜)组成的FGM圆柱壳的径厚比(R/h)、梯度指数(k)、环向模态数(m)、轴向模态数(n)等对临界荷载的影响。
1. FGM圆柱壳的屈曲控制方程
如图 1所示,FGM圆柱壳长度为l,半径为R,总厚度为h,选取柱坐标系(x, θ, z),其相应位移为(u, v, w)。FGM的材料属性(弹性模量E、密度ρ、泊松比μ等)呈幂律分布[2-3],表示为
P(z)=(P1−P2)(2z+h2h)k+P2 (1) 式中:P为物性参数,下标“1”和“2”分别代表组分1和组分2;k为梯度指数,k∈(0, ∞)。圆柱壳内任意点的物性参数为
{E(z)=(E1−E2)(2z+h2h)k+E2ρ(z)=(ρ1−ρ2)(2z+h2h)k+ρ2μ(z)=(μ1−μ2)(2z+h2h)k+μ2 (2) 根据Donnell壳体理论,圆柱壳的小挠度几何方程为
{εx=ε0x+zKxεθ=ε0θ+zKθγxθ=γ0xθ+zKxθ,{u=u0−z∂w∂xv=v0−z∂wR∂θw=w0,{ε0x=∂u0∂xε0θ=∂v0R∂θ−w0Rγ0xθ=∂u0R∂θ−∂v0∂x,{Kx=−∂2w0∂x2Kθ=−∂2w0R2∂θ2Kxθ=−2∂2w0R∂x∂θ (3) 式中:ε为正应变,γ为切应变,上、下标“0”表示壳体中面,K为壳体曲率。
根据经典板壳理论,FGM圆柱壳的内力N与内力矩M可表示为
(NxNθNxθMxMθMxθ)=(A11A12A16B11B12B16A21A22A26B21B22B26A16A26A66B16B26B66B11B12B16D11D12D16B21B22B26D21D22D26B16B26B66D16D26D66)(ε0xε0θγ0xθKxKθKxθ) (4) 式中:Aij、Bij、Dij(i, j=1, 2, 6)分别为FGM圆柱壳的拉伸刚度、耦合刚度和弯曲刚度系数矩阵分量。A11=A22=∫h/2−h/2E(z)1−μ2(z)dz, A12=A21=∫h/2−h/2μ(z)E(z)1−μ2(z)dz, A66=∫h/2−h/2E(z)2[1+μ(z)]dz, B11=B22=∫h/2−h/2E(z)1−μ2(z)zdz, B12=B21=∫h/2−h/2μ(z)E(z)1−μ2(z)zdz, B66=∫h/2−h/2E(z)2[1+μ(z)]zdz, D11=D22=∫h/2−h/2E(z)1−μ2(z)z2dz, D12=D21=∫h/2−h/2μ(z)E(z)1−μ2(z)z2dz, D66=∫h/2−h/2E(z)2[1+μ(z)]z2dz。FGM圆柱壳的力学性能为各向同性[15],那么:A16=A26=B16=B26=D16=D26=0。
对于圆柱壳,系统的应变能(不考虑剪力)为
U=12∫2π0∫l0(Nxε0x+Nθε0θ+Nxθγ0xθ+MxKx+MθKθ+MxθKxθ)Rdxdθ (5) 动能为
T=12∫h/2−h/2∫2π0∫l0ρ(z)[(∂u∂t)2+(∂v∂t)2+(∂w∂t)2]Rdxdθdz (6) 外力功为
W=12∫2π0∫l0N(t)(∂w0∂x)2Rdxdθ (7) Hamilton变分原理为
δ∫t1t2(T−U+W)dt=0 (8) 将(3)式~(7)式代入(8)式中,由Donnell壳体理论可知,圆柱壳内力沿环向均匀分布,忽略中面位移[16],由u0、v0、w0的变分系数为零,整理得到FGM圆柱壳的动力屈曲控制方程为
4I0∂2w0∂t−I2(∂4w0R2∂θ2∂t2+∂4w0∂x2∂t2)=−4A22w0R2−4B12R2∂2w0∂x2−4B22R∂2w0R2∂θ2−D11∂4w0∂x4−2(D12+2D66)∂4w0R2∂x2∂θ2−D22∂4w0R4∂θ4+N(t)∂2w0∂x2 (9) 2. 控制方程的求解
设径向位移表示为[17]
w=Y(x)T(t)eimθ (10) 将(10)式代入(9)式中,分离变量得
{Y(4)=α2Y″+β2Y=0¨T−λT=0 (11) 其中
α2=1D11[4B12R2+2R2(D12+2D66)m2−N(t)−I2λ] (12) β2=1D11[4A22R2+4B22R3+D22R4m4+4I0λ−I2R2m2λ] (13) 当α4>4β2>0且λ>0时,圆柱壳屈曲[18-20],其动力屈曲解为
Y(x)=C1sin(k1x)+C2cos(k1x)+C3sin(k2x)+C4cos(k2x) (14) 式中:C1~C4为系数,k1=√α2−√α4−4β22,k2=√α2+√α4−4β22。
(14)式满足下列两种边界条件:
(1) 对于一端夹支另一端固支的圆柱壳,其边界条件为
{Y(0)=Y′(0)=0Y(l)=Y′(l)=0 (15) (2) 对于一端简支另一端固支的圆柱壳,其边界条件为
{Y(0)=Y″(0)=0Y(l)=Y′(l)=0 (16) 将(14)式代入(15)式中,整理得到如下齐次线性方程组
(0101k10k20sin(k1l)cos(k1l)sin(k2l)cos(k2l)k1cos(k1l)−k1sin(k1l)k2cos(k2l)−k2sin(k2l))(C1C2C3C4)=0 (17) 若要(17)式有非平凡解,其系数行列式必为零,于是
2k1k2−2k1k2cos(k1l)cos(k2l)−(k21+k22)sin(k1l)sin(k2l)=0 (18) 由k1和k2可得
{k21+k22=(n21+n22)π2/l2=α2k21k22=n21n22π4/l4=β2 (19) 将(19)式代入(12)式和(13)式中,得一端夹支另一端固支时FGM圆柱壳动力屈曲临界荷载Ncr,即
Ncr=D11π2(n21+n22)l2−B12R2+2(D12+2D66)m2R2+h2(D11n21n22π4R4−A22l4R2+4B22m2R2l4−D22l4m4)l4R2(h2m2−12R2) (20) 式中:n1=n=1, 2, 3, …;m=1, 2, 3, …;n2=n+2。
同理可得一端简支另一端固支时FGM圆柱壳动力屈曲临界荷载
Ncr=D11π2(n21+n22)l2−B12R2+2(D12+2D66)m2R2+h2(D11n21n22π4R4−A22l4R2+4B22m2R2l4−D22l4m4)l4R2(h2m2−12R2) (21) 此时,n1=n=1, 2, 3, …;m=1, 2, 3, …;n2=n+1。
将FGM退化成金属材料,得到金属材料圆柱壳动力屈曲临界荷载
Ncr=D11π2(n21+n22)l2+2D12m2R2+h2(D11n21n22π4R4−A22l4R2−D22l4m4)l4R2(h2m2−12R2) (22) (22)式与文献[16]中的表达式相同。
根据(10)式,取一端夹支另一端固支时圆柱壳动力屈曲解的表达式[16]
w=T(t)[sin(n1πxl)−n1n2sin(n2πxl)]sin(mθ) (23) 将(23)式代入控制方程(9)式中,计算并化简整理得到临界荷载表达式
Ncr=(n21+n22)D11π2l2+2m2(D12+2D66)R2 (24) (24)式与不考虑转动惯量时用分离变量得到的结果相同,此时n1=n=1, 2, 3, …; m=1, 2, 3, …; n2=n+2。同理可得当边界条件为一端简支另一端固支时的临界荷载表达式,与(24)式相同,此时n1=n=1, 2, 3, …; m=1, 2, 3, …; n2=n+1。
3. 算例分析
采用MATLAB软件编程,对FGM圆柱壳动力屈曲临界荷载进行计算。讨论由不同材料(陶瓷-钛、陶瓷-铁、陶瓷-铜)组成的FGM圆柱壳(见图 2)的径厚比(R/h)、梯度指数(k)、环向模态数(m)、轴向模态数(n)对临界荷载Ncr的影响。基本材料参数如表 1所示。
表 1 材料参数Table 1. Material parametersMaterial E/GPa ρ/(g·cm-3) μ Ceramic 385 3.96 0.230 Ti 109 4.54 0.410 Fe 155 7.86 0.291 Cu 119 8.96 0.326 图 3表示n=1、m=2、k=1、R/h=20时,不同材料组成下Ncr与临界长度l(本研究中临界长度即为圆柱壳长度)的关系曲线。从图 3可以看出:Ncr随l的增加而减小;当l<0.5 m时,Ncr随l的增加而迅速减小;当l>0.5 m时,Ncr随l的增大缓慢减小,且逐渐趋于常数;同一l下,陶瓷-铜的Ncr最大,陶瓷-铁次之,陶瓷-钛的Ncr最小。以下均以陶瓷-钛为例进行讨论。
图 4和图 5分别表示冲击端为夹支和简支时n=1、m=2、k=1时不同R/h下Ncr与l的关系曲线。可以看出:当l增加时,Ncr减小,且逐渐趋于常数;在同一l下圆柱壳的Ncr随着R/h的增大而减小;当R/h和l一定时,冲击端为夹支时的Ncr明显比冲击端为简支时的Ncr大,说明约束条件对Ncr有较大影响。
图 6和图 7分别表示冲击端为夹支和简支条件下n=1、m=2、R/h=20时不同k下Ncr与l的关系曲线。可见:FGM圆柱壳的Ncr随着l的增加而减小;在同一l下,FGM圆柱壳的Ncr随着k的增加而增加;当l<0.5 m时,Ncr随l的增加迅速减小,当l>0.5 m时,Ncr随l的增加缓慢减小并逐渐趋于常数;当k=1且l一定时,冲击端为夹支时的Ncr明显比冲击端为简支时的Ncr大,再次说明约束条件对Ncr的影响较大。
图 8和图 9分别表示冲击端为夹支和简支,n=1、k=1、R/h=20时不同m下Ncr与l的关系曲线。可以看出:当l在一定范围内时,Ncr随l的增加而迅速减小,超出这一范围后Ncr随l的增加而缓慢减小且逐渐趋于常数;同一l下,随着m的增大,FGM圆柱壳的Ncr增大,表明Ncr越大,高阶模态越易被激发。当m=6且l一定时,冲击端为夹支时的Ncr明显比冲击端为简支时的Ncr大,表明约束条件对Ncr有较大影响。
图 10和图 11分别表示冲击端为夹支和简支,m=1、k=1、R/h=20时不同n下Ncr与l的关系曲线。图 10和图 11显示:Ncr随着l的增加而减小;不同n条件下,Ncr随l的增加逐渐趋于同一值;当l<1 m时,在同一l下Ncr随n的增加而增加,说明Ncr越大,高阶模态越容易被激发。
图 12为不同环向模态数m下FGM圆柱壳的动力屈曲模态。可以看出:随着m的增大,圆柱壳的模态变得越来越复杂,俯视图由单一形变为多瓣形;当m=6时,俯视图为12瓣形。
图 13为不同轴向模态数n下FGM圆柱壳的动力屈曲模态图。由图 13可知:随着n的增加,模态图变得越来越复杂。由FGM圆柱壳的俯视图可知,各阶模态数下动力屈曲模态图为轴对称。
4. 结论
(1) 根据Donnell壳体理论和经典板壳理论,由Hamilton变分原理得到轴向荷载作用下FGM圆柱壳的动力屈曲控制方程。
(2) 由圆柱壳周向连续性设出径向位移的周向形式,并用分离变量法得到不同约束条件下FGM圆柱壳动力屈曲临界荷载的表达式和屈曲解式。
(3) 利用MATLAB对临界荷载进行计算,得到:在轴向模态数(n)、环向模态数(m)、梯度指数(k)、径厚比(R/h)一定的情况下,同种材料组成的圆柱壳的临界荷载随着临界长度的增加而减小;在n、m、k、l一定的情况下,临界荷载随着径厚比的增大而减小;在n、m、R/h、l一定的情况下,临界荷载随着梯度指数k的增加而增加;不同约束条件下,冲击端为夹支的临界荷载大于冲击端为简支的临界荷载,表明约束条件对临界荷载有较大影响;圆柱壳的临界荷载随模态数的增加而增大,表明临界荷载越大,越容易激发高阶模态;圆柱壳的动力屈曲模态随模态数的增加变得更为复杂。
-
表 1 α-Ir2P相和β-Ir2P相的平衡态晶格常数和原子位置
Table 1. Lattice parameters and atomic coordinate of α-Ir2P and β-Ir2P
Phase Pressure/GPa Space group Lattice parameters Wyckoff position Atoms Site α-Ir2P 0 Fm3m a=5.622 Å(5.535 Å*), b=c=5.622 Å Ir1 8c(0.250, 0.250, 0.250) α=β=γ=90.0° P1 4a(0, 0, 0) β-Ir2P 86.4 I4/mmm a=b=2.694 Å, c=9.461 Å Ir1 4e(0.500, 0.500, 0.146) α=β=γ=90.0° P1 2a(0.500, 0.500, 0.500) Note: The asterisk represents the experimental data from Ref. [13]. 表 2 86.4 GPa下β-Ir2P相的Bader电荷转移
Table 2. Calculated Bader charges of β-Ir2P phase at 86.4 GPa
Space group Pressure/GPa Atom Number Charge value/e Charge transfer/e I4/mmm 86.4 GPa Ir 2 9.19 –0.19 P 1 4.62 0.38 -
[1] HENKES A E, VASQUEZ Y, SCHAAK R E. Converting metals into phosphides: a general strategy for the synthesis of metal phosphide nanocrystals [J]. Journal of the American Chemical Society, 2007, 129(7): 1896–1897. doi: 10.1021/ja068502l [2] MAUVERNAY B, DOUBLET M L, MONCONDUIT L. Redox mechanism in the binary transition metal phosphide Cu3P [J]. Journal of Physics and Chemistry of Solids, 2006, 67(5/6): 1252–1257. [3] BROCK S L, SENEVIRATHNE K. Recent developments in synthetic approaches to transition metal phosphide nanoparticles for magnetic and catalytic applications [J]. Journal of Solid State Chemistry, 2008, 181(7): 1552–1559. doi: 10.1016/j.jssc.2008.03.012 [4] OYAMA S T, GOTT T, ZHAO H, et al. Transition metal phosphide hydroprocessing catalysts: a review [J]. Catalysis Today, 2009, 143(1/2): 94–107. [5] HALL J W, MEMBRENO N, WU J, et al. Low-temperature synthesis of amorphous FeP2 and its use as anodes for Li ion batteries [J]. Journal of the American Chemical Society, 2012, 134(12): 5532–5535. doi: 10.1021/ja301173q [6] BILTZ W, WEIBKE F, MAY E, et al. Alloyability of platinum and phosphorus [J]. Zeitschrift fur Anorganische und Allgemeine Chemie, 1935, 223(2): 129–143. doi: 10.1002/zaac.v223:2 [7] ZHANG X, QIN J, SUN X, et al. First-principles structural design of superhard material of ZrB4 [J]. Physical Chemistry Chemical Physics, 2013, 15(48): 20894–20899. doi: 10.1039/c3cp53893a [8] KANER R B, GILMAN J J, TOLBERT S H. Designing superhard materials [J]. Science, 2005, 308(5726): 1268–1269. doi: 10.1126/science.1109830 [9] SHI Y, ZHANG B. Correction: Recent advances in transition metal phosphide nanomaterials: synthesis and applications in hydrogen evolution reaction [J]. Chemical Society Reviews, 2016, 45(6): 1781–1781. doi: 10.1039/C6CS90013E [10] CARENCO S, PORTEHAULT D, BOISSIERE C, et al. Nanoscaled metal borides and phosphides: recent developments and perspectives [J]. Chemical Reviews, 2013, 113(10): 7981–8065. doi: 10.1021/cr400020d [11] LI W, DHANDAPANI B, OYAMA S T. Molybdenum phosphide: a novel catalyst for hydrodenitrogenation [J]. Chemistry Letters, 1998, 27(3): 207–208. doi: 10.1246/cl.1998.207 [12] OYAMA S T, GOTT T, ZHAO H, et al. Transition metal phosphide hydroprocessing catalysts: a review [J]. Catalysis Today, 2009, 143(1/2): 94–107. [13] ZUMBUSCH M. Über die strukturen des uransubsulfids und der subphosphide des iridiums und rhodiums [J]. Zeitschrift für Anorganische und Allgemeine Chemie, 1940, 243(4): 322–329. doi: 10.1002/zaac.19402430403 [14] RUNDQVIST S. Phosphides of the platinum metals [J]. Nature, 1960, 185(4705): 31. doi: 10.1038/185031a0 [15] RAUB C J, ZACHARIASEN W H, GEBALLE T H, et al. Superconductivity of some new Pt-metal compounds [J]. Journal of Physics and Chemistry of Solids, 1963, 24(9): 1093–1100. doi: 10.1016/0022-3697(63)90022-2 [16] WANG P, WANG Y, WANG L, et al. Elastic, magnetic and electronic properties of iridium phosphide Ir2P [J]. Scientific Reports, 2016, 6(9): 21787. [17] SUN X W, BIOUD N, FU Z J, et al. High-pressure elastic properties of cubic Ir2P from ab initio calculations [J]. Physics Letters A, 2016, 380(43): 3672–3677. doi: 10.1016/j.physleta.2016.08.048 [18] LIU Z J, SONG T, SUN X W, et al. Thermal expansion, heat capacity and Grüneisen parameter of iridium phosphide Ir2P from quasi-harmonic Debye model [J]. Solid State Communications, 2017, 253: 19–23. doi: 10.1016/j.ssc.2017.01.028 [19] WANG Y, LV J, ZHU L, et al. Crystal structure prediction via particle-swarm optimization [J]. Physical Review B, 2010, 82(9): 094116. doi: 10.1103/PhysRevB.82.094116 [20] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188. doi: 10.1103/PhysRevB.13.5188 [21] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865 [22] BORN M, HUANG K. Dynamical theory of crystal lattices [J]. American Journal of Physics, 1954, 39(2): 113–127. [23] SAVIN A, NESPER R, WENGERT S, et al. ELF: the electron localization function [J]. Angewandte Chemie International Edition in English, 1997, 36(17): 1808–1832. doi: 10.1002/(ISSN)1521-3773 -