Synthesis of 6H-Type Hexagonal Perovskite Phase of BaGeO3 at High Temperature and High Pressure
-
摘要: 利用金刚石对顶砧高压装置,结合显微激光双面加热技术,对BaGeO3开展了高温高压实验研究。常温常压下赝硅灰石相的BaGeO3于12 GPa左右开始非晶化。进一步加压到22 GPa并对已完全非晶化的BaGeO3样品进行(1800 ± 200) K的高温处理,拉曼光谱显示其转变成一种未见报道的高压新相。在0~17.4 GPa压力范围对BaGeO3高压新相开展同步辐射X射线衍射测试,其衍射谱可以用6H型六方钙钛矿相进行指标化,并且卸压到常压时仍保持稳定。以6H型钙钛矿相为结构模型,分别对17.4 GPa和常压下的X射线衍射谱进行Rietveld结构精修,获得其结构参数。应用二阶Birch-Murnaghan状态方程拟合实验体积-压力数据,得到其体弹模量K0 = 150(2) GPa和零压晶胞体积V0 = 373.0(3) Å3。在实验研究的基础上,对6H型钙钛矿相BaGeO3进行第一性原理理论计算,所得不同压力下的晶格常数和体积数据与实验结果符合得很好,状态方程参数K0 = 153(1) GPa,V0 = 374.2(1) Å3。20.0 GPa时计算的拉曼光谱也很好地描述了拉曼实验测量结果。研究结果补充了赝硅灰石相BaGeO3在更高温压条件下的结构相转变。6H型钙钛矿相BaGeO3的获得为进一步表征该相的物理化学性质奠定了基础,为开发高性能钙钛矿结构锗酸盐材料提供了可能性,同时对于理解硅酸盐钙钛矿结构的相变规律及稳定性、地球下地幔物理化学性质及其变化等具有重要的指示意义。
-
关键词:
- BaGeO3 /
- 6H型六方钙钛矿相 /
- Rietveld结构精修 /
- 拉曼光谱 /
- 高温高压
Abstract: Barium germanate (BaGeO3) was studied using double-sided laser-heating diamond anvil cell (LHDAC). At ambient conditions, BaGeO3 has a pseudowollastonite structure. At about 12 GPa, BaGeO3 crystal begin to translate into an amorphous phase. The amorphous BaGeO3 was further pressurized to about 22 GPa and then heated at (1800 ± 200) K conditions. Raman spectra shows the amorphous BaGeO3 transforms into a new high pressure phase, which has not been reported so far. The new high pressure phase of BaGeO3 was further measured with the synchrotron radiation X-ray diffraction in the pressure ranges of 0−17.4 GPa. The diffraction patterns can be indexed with a 6H-type hexagonal perovskite structure and this structure remained stable as the pressure unloading to ambient pressure. In order to obtain the structural parameters of the new high pressure phase of BaGeO3, the X-ray diffraction patterns of 17.4 GPa and ambient pressure were refined with a model structure of 6H-type perovskite using the Rietveld method. The experimental pressure-volume data was fitted with the second-order Birch-Murnaghan equation of state, and obtained the volume bulk modulus and zero-pressure unit-cell volume are K0 = 150(2) GPa and V0 = 373.0(3) A3 respectively. On the basis of the experimental results in this study, we also carried out the first-principle theoretical calculation on the 6H-type perovskite BaGeO3. The calculated lattice constants and volume with the corresponding pressures are good agreement with the experimental results. Furthermore, the calculated volume bulk modulus and zero-pressure unit-cell volume are K0 = 153(1) GPa, V0 = 374.2(1) A3 respectively. The calculated Raman spectra at 20.0 GPa is also well consistent with the experimental results. This study not only complements the structural phase transition of pseudowallastonite BaGeO3 at high temperature and high pressure, but also builds a solid foundation for further characterizing the physical and chemical properties of pseudowallastonite BaGeO3, and gives a chance to develop the perovskite structured germanate functional materials. In addition, this study has an important indicative significance for us to understand the phase transition rule and stability of silicate perovskite, the physical and chemical properties and changes of Earth's lower mantle. -
半穿甲导弹是在穿甲弹基础上发展起来的用于攻击舰船目标的武器,兼具穿甲和爆破毁伤能力,是舰船防护领域的重点设防对象[1]。前舱物作为半穿甲导弹的重要组成部分,在导弹的制导和控制方面起着至关重要的作用。半穿甲导弹侵彻靶板的过程中,前舱物先于战斗部接触靶板,对靶板产生影响。因此,研究前舱物对战斗部侵彻靶板的影响具有实际意义。迄今为止,对于半穿甲导弹侵彻舰船结构等问题,研究人员通常忽略前舱物的影响,即只考虑战斗部的动能穿甲作用,将其抽象为低速大质量弹体的穿甲问题。陈斌等[2]利用忽略前舱物的半穿甲模拟弹,开展了侵彻陶瓷/钢复合装甲试验,明确了弹着角对半穿甲弹侵彻效能有明显影响;朱锡等[3]采用与实际半穿甲导弹战斗部几何形状相似的弹丸,开展了侵彻多种装甲结构试验,得到了不同装甲结构形式抗半穿甲导弹动能穿甲的特性。目前,针对前舱物对半穿甲导弹侵彻靶板影响的研究相对较少,这是由于前舱物的组成部件复杂,无法完全模拟所有器件,因此,通常对前舱物进行力学等效。
对前舱物的等效主要有两种方法。第一种方法为采用质量块对前舱物进行等效。楼建锋等[4]通过有限元方法,分别使用金属材料和脆性材料质量块等效模拟前舱物,得到了高速(675 m/s)和低速(250 m/s)两种速度条件下前舱物对战斗部侵彻过程的影响。陈刚[5]利用质量块模拟前舱物,对比分析了前舱物对战斗部高速撞击靶板的影响,得到了前舱物改进战斗部弹体在斜穿甲过程中的荷载环境条件。Chen等[6]考虑前舱物作用导致的靶板结构预响应,将前舱物剩余部分等效为附着在战斗部前部的质量块,并对剪切冲塞模型进行修正,得到了与实验结果相吻合的终点弹道公式。该等效方法的关键是质量块材料的选取,只有选择与前舱物名义力学特性相近的材料,才能正确表征前舱物对导弹弹体动能穿甲的影响。第二种方法为采用薄壁圆筒结构对前舱物进行等效。徐钰巍等[7]采用硬铝圆筒结构对前舱物进行等效,建立了变截面刚度/含缺陷的自由梁模型,并结合模型试验得到了前舱物在弹体侵彻靶板初期对弹体姿态变化的影响规律。该等效方法忽略了前舱物内部仪器设备强度的影响,仅考虑前舱段外壳在侵彻过程中所起的作用,在研究半穿甲导弹正侵彻时,与实际情况的差异较大。对于正侵彻而言,将前舱物等效为质量块的方法是可行的,但现有的等效方法为了体现前舱物稀松多孔、中空易碎的结构特点,所选取的等效材料模型和前舱物在撞击过程中所表现出的力学特性存在不相符的情况。
本研究采用轻质泡沫铝材料模拟半穿甲导弹前舱物,建立变截面刚度弹体侵彻金属薄板的有限元模型,分析带前舱物平头弹侵彻金属薄板的动态力学过程,研究不同靶板厚度、弹速条件下前舱物对平头战斗部侵彻金属薄板的影响。
1. 计算模型与材料参数
1.1 前舱物特性
前舱物主要由雷达导引头、导航系统、电源系统、高度表等仪器设备组成,其特点是质量较大(约为战斗部总质量的1/4~1/3)且分布均匀,整体结构中空、易碎、延展性差,结构强度远低于战斗部。前舱物由高分子聚合物、铝合金、单晶硅等多种材料构成,平均密度低(小于1 g/cm3)。前舱物在撞击靶板时,会经历“接触—屈服—压实—破碎飞散”的破坏过程,其中屈服阶段持续时间长,名义应力基本保持不变,即存在较长的平台应力,随后前舱物被压实,名义应力大幅增加,前舱段整体长度大幅减小,这与泡沫铝材料在冲击荷载作用下的位移变化曲线类似[8-10],如图1所示。因此,本研究中前舱物材料模型采用低强度、低密度、有长应力平台的泡沫铝材料进行等效。
1.2 带前舱物战斗部等效计算模型
半穿甲反舰导弹通常由前舱物、战斗部和推进系统3部分组成。战斗部是弹体中强度最高、刚度最大的部分,有较大的药室,装填炸药量较多,大多为钝头[1]。本研究建立的前舱物战斗部等效计算模型如图2所示。战斗部壳体模型外径为428.40 mm,内径为384.79 mm,内置装药模型直径为384.79 mm,前舱物等效为与战斗部截面尺寸相等的圆柱体。
采用有限元软件LS-DYNA模拟不同工况下带前舱物战斗部侵彻靶板的过程。计算中,结构采用拉格朗日实体单元进行网格划分,时间积分采用显式中心差分法。靶板平面尺寸为3 000 mm × 3 000 mm,边界条件设置为固支约束。计算模型设置为沿中轴面对称,建立1/2模型,战斗部与前舱物采用六面体网格划分。靶板自中心采用过渡网格,中心区域网格尺寸为2.5 mm × 2.5 mm,距离撞击区域越远,网格尺寸越大。
1.3 本构方程与材料参数
数值仿真中主要涉及两种材料模型。战斗部壳体和金属靶板采用金属材料常用的Johnson-Cook材料模型,战斗部内置装药和前舱物等效材料采用弹塑性材料常用的弹塑性硬化模型。
Johnson-Cook材料模型考虑了塑性强化、应变率强化、温升软化效应,能够反映材料在高温、高应变率下力学特性的变化,其具体形式为
s=(A+Bˉepn)(1+Cln˙ε˙ε0)[1−(T−T0Tm−T0)m] (1) 式中:s为应力;A为静态屈服强度;B、n通过拟合等效应力-应变曲线数据获得;
${\bar {e}_{\rm{p}}}$ 为有效塑性应变;C为应变率系数;$\dot \varepsilon$ 为有效塑性应变率;${\dot \varepsilon _0}$ 为参考塑性应变率,一般取${\dot \varepsilon _0}$ = 1 s–1;T0为参考温度(取室温);Tm为材料熔点;m为热软化指数。战斗部壳体和金属靶板的材料参数根据文献[11-12]并结合工程实践调整得到,如表1所示,其中cp为定压比热容,ρ为密度,μ为泊松比。
表 1 弹体和靶板的J-C模型参表Table 1. J-C model parameters of projectile body and targetMaterial ρ/(g·cm–3) E/GPa μ cp/(J·kg–1·K–1) T0/K Tm/K Warhead shell 7.80 205 0.28 400.90 300 1 765 Target plate 7.80 200 0.30 452.00 298 1 881 Material $\dot \varepsilon $/s–1 A/MPa B/MPa n C m Warhead shell 1 760 500 0.53 0.014 1.13 Target plate 1 355 450 0.36 0.022 1.00 弹塑性硬化模型考虑了材料应变的等向或随动塑性强化。其应变率效应由Cowper-Symonds模型描述
σy=[1+(˙εc)1/p](σ0+βEEtE−Etεeff) 式中:σy为屈服强度,
$\dot \varepsilon $ 为等效塑性应变率;c、p、β为材料常数;σ0为静态屈服强度,E为杨氏模量,Et为硬化模量,εeff为有效塑性应变。内置装药和前舱物等效泡沫材料的模型参数如表2所示,其中fs为失效应变。考虑到不同导弹的前舱物内包含仪器设备的数量和种类不同,前舱物等效材料的屈服强度无法准确给出,本研究只给出取值范围。表 2 弹塑性硬化模型材料参数Table 2. Material parameters of the Elastoplastic hardening modelMaterial ρ/(kg·m–3) E/GPa μ σy/MPa Et/MPa β c/s–1 p fs Charge 1 750 5 0.3 20 300 0 0 0 5 Nose cabin 700 0.5 0.3 0.5–10.0 0 0 0 0 0.4 2. 结果分析
2.1 带前舱物战斗部侵彻过程分析
图3为带前舱物战斗部侵彻金属薄板的过程。图4为侵彻过程中典型时刻的靶板内应力分布。典型工况下,带前舱物战斗部与不计前舱物战斗部侵彻金属薄板的速度时程曲线如图5所示。根据金属薄板受力和变形的特点,可以将带前舱物战斗部侵彻金属薄板的过程分为3个阶段。
(1)金属薄板结构预变形阶段。从前舱物接触靶板开始,接触区域内最大压强远大于前舱物等效材料的极限强度和金属薄板的屈服强度。因此,前舱物发生变形破坏,金属薄板产生大范围碟形塑性变形,并获得一定的横向速度,战斗部与靶板之间的速度差减小,弹靶接触区边缘的应力集中程度减轻。由图4可知,在金属薄板结构预变形阶段,前舱物与靶板接触区域内正应力较大,且随着与中心距离的增大而降低;剪应力很小,且只存在于接触区域及其附近区域。在此阶段,战斗部的速度变化可以分为两个时期,如图5所示。AB时期,由于前舱物与金属薄板撞击产生的压缩应力波尚未传递到战斗部,战斗部以初速度向前运动;BC时期,随着前舱物与金属薄板撞击产生的压缩应力作用于战斗部,战斗部的速度开始下降,但由于作用于战斗部的压缩应力强度较小,战斗部速度下降较慢。
(2)冲塞破坏阶段。从战斗部接触金属薄板开始,接触区域由于速度差产生剪切变形,在拉应力与剪应力的联合作用下,接触区域周围产生裂纹,裂纹直径略大于战斗部直径。如图4所示,此阶段战斗部与靶板接触区域内正应力最大,且由于裂纹拓展释放压缩波,接触区域内出现拉伸应力与压缩应力混合的情况;接触区域内剪应力强度远大于结构预变形阶段,剪应力在接触边缘附近达到峰值,并随着与中心距离的增大迅速衰减。如图5所示,在此阶段,战斗部直接侵彻金属薄板,弹体受到的应力强度大,速度下降明显,可近似视为匀减速运动。
(3)弹体贯穿阶段。随着弹体的进一步运动,裂纹贯穿整个金属薄板厚度方向,形成与战斗部直径相似的帽形失效块,同时,金属薄板由于惯性作用继续发生碟形变形。如图4所示,该阶段由于战斗部与靶板接触的区域形成失效块与靶板脱离,接触区域内正应力与剪应力均为零。靶板动能转化为变形能,挠度增加,靶板仍有较大正应力,由于接触区域附近靶体材料存在速度差,因此靶板接触区域附近存在一定的剪应力。如图5所示,在该阶段,战斗部的速度变化呈现两个时期:DE阶段,由于冲塞块附着在战斗部头部,战斗部推动冲塞块前进,战斗部的一部分动能转化为冲塞块的动能,导致战斗部速度略有下降;EF阶段,冲塞块与战斗部速度达到一致,战斗部向前匀速运动。
由图5中带前舱物战斗部与不计前舱物战斗部的速度时程曲线对比可知,不计前舱物战斗部侵彻靶板的贯穿时间为90 μs,带前舱物战斗部的侵彻贯穿时间为2 180 μs,侵彻作用时间大幅增加,但增加的侵彻时间主要是前舱物被压溃破坏的时间。在金属薄板预变形阶段,金属薄板与前舱物接触区域获得的速度较小,与战斗部之间的速度差较大,因此,战斗部侵彻贯穿金属薄板所需要的时间并未明显增加。此外,相比于不计前舱物战斗部速度变化情况,前舱物的存在使得战斗部在侵彻过程中的受力更加复杂,战斗部速度变化的阶段性和多样性更加明显。
从塑性变形范围来看,通过有限元仿真软件的计算结果可知,不计前舱物战斗部侵彻靶板过程[13] 中金属薄板塑性变形的直径为883.8 mm,带前舱物战斗部侵彻靶板过程中,由于存在金属薄板预变形阶段,金属薄板迎弹面塑性变形范围大幅增加,直径达到2 629.7 mm。其原因在于:带前舱物战斗部侵彻金属薄板的过靶时间长,弹靶接触区域产生的塑性波传播范围远大于不计前舱物战斗部。因此,前舱物的“缓冲垫”作用明显,在大幅延长战斗部过靶时间的同时大幅扩大了金属薄板塑性变形范围。
2.2 前舱物等效模型屈服强度的影响
由于不同导弹的前舱物内包含仪器设备的种类和数量不同,因此,无法定量给出前舱物等效材料的屈服强度。本研究给定前舱物等效材料的屈服强度范围为0.5~10.0 MPa,在此屈服强度范围内研究前舱物对平头战斗部侵彻金属薄板的影响。
图6给出了初速度为340 m/s的不同前舱物等效屈服强度战斗部侵彻20 mm厚金属薄板的应变云图。由图6可知,前舱物等效屈服强度在0.5~10.0 MPa范围内时,金属薄板的塑性应变分布相似,战斗部贯穿金属薄板形成的帽形失效块直径相等。表3列出了不同前舱物等效屈服强度战斗部侵彻靶板时金属薄板的塑性变形范围和战斗部剩余速度。由表3可知,前舱段等效屈服强度增加时,金属薄板塑性变形范围并未出现明显变化,而是在一定范围内波动,最大波动幅度为0.994%。其原因在于,模型计算结果的输出时间间隔为20 μs,侵彻作用时间选取最接近战斗部贯穿靶板的时间,时间选取存在不准确的情况。随着前舱物等效屈服强度的增加,战斗部剩余速度增大,但增大幅度很小,仅为0.147%。因此,可以认为前舱物等效屈服强度对带前舱物平头弹侵彻金属薄板的影响很小,可以忽略。这是由于前舱物存在破坏强度低、延展性差的特点,变形破坏所耗散的变形能远小于带前舱物战斗部所具有的动能,因此当前舱物等效模型的屈服强度远小于靶板时,前舱物的惯性力起主要作用,结构强度的影响可以忽略不计。
表 3 不同等效屈服强度下侵彻结果的比较Table 3. Comparison of penetration results with different yield stressEquivalent yield strength/MPa Distortion range/mm Residual velocity/(m·s–1) 0.5 2 629.70 325.88 1.0 2 624.41 326.30 2.0 2 629.31 326.37 5.0 2 625.23 326.75 10.0 2 659.20 326.86 2.3 靶板厚度和侵彻速度的影响
为了分析不同靶板厚度时前舱物对战斗部侵彻的影响,分别模拟了靶板厚度在10~40 mm范围内的7种工况。不同靶板厚度条件下,战斗部剩余速度如图7所示。由图7可知:带前舱物战斗部剩余速度与不计前舱物战斗部的剩余速度均随靶板厚度的增大而减小,且变化规律一致,均近似呈线性变化;靶板厚度相同时,两者的剩余速度相差很小,最大偏差为0.67%。因此,靶板厚度在10~40 mm范围内时,前舱物对战斗部侵彻的影响很小。
常见反舰半穿甲导弹的飞行马赫数Ma 在 0.8~1.0之间[14],为了分析不同侵彻速度下前舱物对平头战斗部侵彻金属薄板的影响,分别模拟了侵彻速度在238~340 m/s范围内的7种情况。不同侵彻速度条件下,战斗部的剩余速度如图8所示。由图8可知,随着侵彻速度的增大,不计前舱物平头战斗部和带前舱物平头战斗部的剩余速度均近似呈线性增长,且带前舱物平头战斗部的剩余速度略大于不计前舱物平头战斗部(最大偏差为1.15%)。Chen等[6]提出的带前舱物钝头弹侵彻金属靶的计算公式可得到与本模型相似的结论:侵彻速度相同时,带前舱物平头战斗部的剩余速度略大于不计前舱物平头战斗部。
通过对不同靶板厚度和侵彻速度条件下战斗部的剩余速度进行分析可知,带前舱物战斗部的剩余速度均大于不计前舱物战斗部,原因在于:与不计前舱物战斗部相比,带前舱物战斗部的质量更大,侵彻速度相同时动能更大。前舱物变形破坏和金属靶板结构预变形消耗部分前舱物具有的动能,使带前舱物战斗部以更大动能侵彻靶板。带前舱物战斗部对靶板的侵彻破坏模式为冲塞破坏,与不带前舱物战斗部侵彻靶板的破坏模式相同,在侵彻过程中消耗的能量相差不大。因此,带前舱物战斗部的剩余速度略大于不计前舱物战斗部。
3. 结 论
本研究建立了变截面刚度弹体模型,对带前舱物战斗部侵彻金属薄板进行仿真分析。通过对比不同前舱物等效屈服强度战斗部的侵彻性能,以及不同靶板厚度与初速度条件下带前舱物战斗部与不计前舱物战斗部的剩余速度,得出如下主要结论。
(1)根据靶板的受力变形特征,带前舱物平头战斗部侵彻金属薄板可以分为3个阶段,即靶板预变形阶段、冲塞破坏阶段和弹体贯穿阶段。虽然带前舱物平头战斗部与不计前舱物平头战斗部侵彻靶板的过程存在明显差异,但靶板破坏模式均为剪切冲塞破坏。
(2)当前舱物等效材料的强度远小于靶板的屈服强度时,前舱物等效泡沫铝材料的屈服强度对战斗部侵彻性能的影响可以忽略不计。
(3)在导弹侵彻目标的过程中,一方面,前舱物的存在增大了弹体的总动能,且先于战斗部撞击目标,造成目标结构变形甚至破坏,对战斗部侵彻靶板起到了积极作用;另一方面,由于前舱物的强度远小于战斗部,前舱物的存在延长了弹体的过靶时间,降低了靶板过载,起到了“缓冲垫”的作用,不利于平头战斗部对目标的侵彻。在这两方面综合影响下,当飞行马赫数Ma 在 0.8~1.0之间,侵彻厚度为10~40 mm的靶板时,带前舱物平头战斗部的剩余速度比不计前舱物平头战斗部均有所增加,但增加幅度小于2%。
-
表 1 6H型六方钙钛矿相BaGeO3在常压和17.4 GPa的结构参数
Table 1. Structural parameters of hexagonal perovskite phase BaGeO3 (6H-type) at ambient pressure and 17.4 GPa
Structural parameters at 17.4 GPa Structural parameters at 0.1 MPa Atom Site x y z Atom Site x y z Ba1 2b 0 0 0.2500 Ba1 2b 0 0 0.2500 Ba2 4f 0.3333 0.6667 0.0896(5) Ba2 4f 0.3333 0.6667 0.0930(5) Ge1 2a 0 0 0 Ge1 2a 0 0 0 Ge2 4f 0.3333 0.6667 0.8392(5) Ge2 4f 0.3333 0.6667 0.8427(5) O1 6h 0.5178(27) 0.0360(5) 0.2500 O1 6h 0.4976(23) −0.0050(5) 0.2500 O2 12k 0.8342(27) 0.6680(5) 0.0730(5) O2 12k 0.8140(23) 0.6280 0.0765 Bond lengths at 17.4 GPa/Å Bond lengths at 0.1 MPa/Å Ge1―O2 Ge2―O2 Ge2―O1 Ge1―O2 Ge2―O2 Ge2―O1 1.837(19) × 6 1.961(16) × 3 1.835(17) × 3 2.088(17) × 6 1.808(14) × 3 2.075(15) × 3 Note: Numbers in parentheses indicate standard deviation. 表 2 计算得到的6H型六方钙钛矿相BaGeO3在20.0 GPa下的拉曼振动模
Table 2. Calculated Raman vibrational modes of hexagonal perovskite phase BaGeO3 (6H-type) at 20.0 GPa cm−1
E2g E1g A1g E2g E1g E2g E2g E1g A1g E1g 88 131 142 158 202 213 242 264 381 395 E2g E2g E1g A1g E2g E2g E1g A1g A1g 399 429 431 444 529 576 578 698 840 -
[1] MIZOGUCHI H, KAMIYA T, MATSUISHI S, et al. A germanate transparent conductive oxide [J]. Nature Communications, 2011, 2: 470. doi: 10.1038/NCOMMS1484 [2] HORIUCHI H, ITO E, WEIDNER D J. Perovskite-type MgSiO3: single-crystal X-ray diffraction study [J]. American Mineralogist, 1987, 72(3/4): 357–360. [3] MAO H K, CHEN L C, HEMLEY R J, et al. Stability and equation of state of CaSiO3-perovskite to 134 GPa [J]. Journal of Geophysical Research: Solid Earth, 1989, 94(B12): 17889–17894. doi: 10.1029/JB094iB12p17889 [4] XIAO W S, TAN D Y, ZHOU W, et al. A new cubic perovskite in PbGeO3 at high pressures [J]. American Mineralogist, 2012, 97(7): 1193–1198. doi: 10.2138/am.2012.4021 [5] XIAO W S, TAN D Y, ZHOU W, et al. Cubic perovskite polymorph of strontium metasilicate at high pressures [J]. American Mineralogist, 2013, 98(11/12): 2096–2104. doi: 10.2138/am.2013.4470 [6] GIBBS G V, BOISEN M B, HILL F C, et al. SiO and GeO bonded interactions as inferred from the bond critical point properties of electron density distributions [J]. Physics and Chemistry of Minerals, 1998, 25(8): 574–584. doi: 10.1007/s002690050150 [7] AKAOGI M, KOJITANI H, YUSA H, et al. High-pressure transitions and thermochemistry of MGeO3 (M = Mg, Zn and Sr) and Sr-silicates: systematics in enthalpies of formation of A2+B4+O3 perovskites [J]. Physics and Chemistry of Minerals, 2005, 32(8/9): 603–613. doi: 10.1007/s00269-005-0034-1 [8] NAKATSUKA A, ARIMA H, OHTAKA O, et al. Crystal structure of SrGeO3 in the high-pressure perovskite-type phase [J]. Acta Crystallographica Section E: Crystallographic Communications, 2015, 71(5): 502–504. doi: 10.1107/S2056989015007264 [9] ROSS N L, ANGEL R J. Compression of CaTiO3 and CaGeO3 perovskites [J]. American Mineralogist, 1999, 84(3): 277–281. doi: 10.2138/am-1999-0309 [10] RUNGE C E, KUBO A, KIEFER B, et al. Equation of state of MgGeO3 perovskite to 65 GPa: comparison with the post-perovskite phase [J]. Physics and Chemistry of Minerals, 2006, 33(10): 699–709. doi: 10.1007/s00269-006-0116-8 [11] YUSA H, AKAOGI M, SATA N, et al. Letter: unquenchable hexagonal perovskite in high-pressure polymorphs of strontium silicates [J]. American Mineralogist, 2005, 90(5/6): 1017–1020. doi: 10.2138/am.2005.1835 [12] YUSA H, SATA N, OHISHI Y. Rhombohedral (9R) and hexagonal (6H) perovskites in barium silicates under high pressure [J]. American Mineralogist, 2007, 92(4): 648–654. doi: 10.2138/am.2007.2314 [13] HIRAMATSU H, YUSA H, IGARASHI R, et al. An exceptionally narrow band-gap (~4 eV) silicate predicted in the cubic perovskite structure: BaSiO3 [J]. Inorganic Chemistry, 2017, 56(17): 10535–10542. doi: 10.1021/acs.inorgchem.7b01510 [14] YANG H X, PREWITT C T. Crystal structure and compressibility of a two-layer polytype of pseudowollastonite (CaSiO3) [J]. American Mineralogist, 1999, 84(11/12): 1902–1905. doi: 10.2138/am-1999-11-1217 [15] NISHI F. Strontium metagermanate, SrGeO3 [J]. Acta Crystallographica Section C: Crystal Structure Communications, 1997, 53(4): 399–401. doi: 10.1107/S0108270196013960 [16] WAN S M, ZENG Y, YAO Y N, et al. BaGeO3: a mid-IR transparent crystal with superstrong raman response [J]. Inorganic Chemistry, 2020, 59(6): 3542–3545. doi: 10.1021/acs.inorgchem.0c00155 [17] GSPAN C, KAHLENBERG V, KOTHLEITNER G, et al. Atomic and domain structure of the low-temperature phase of barium metagermanate (BaGeO3) [J]. Acta Crystallographica Section A: Foundations of Crystallography, 2006, 62(6): 1002–1009. doi: 10.1107/S0108768106039140 [18] SHIMIZU Y, SYONO Y, AKIMOTO S. High-pressure transformations in SrGeO3, SrSiO3, BaGeO3, and BaSiO3 [J]. High Temperatures-High Pressures, 1970, 2(1): 113–120. [19] OZIMA M, SUSAKI J I, AKIMOTO S I, et al. The system BaO-GeO2 at high pressures and temperatures, with special reference to high-pressure transformations in BaGeO3, BaGe2O5, and Ba2Ge5O12 [J]. Journal of Solid State Chemistry, 1982, 44(3): 307–317. doi: 10.1016/0022-4596(82)90378-4 [20] GASPARIK T, WOLF K, SMITH C M. Experimental determination of phase relations in the CaSiO3 system from 8 to 15 GPa [J]. American Mineralogist, 1994, 79(11/12): 1219–1222. [21] AKAOGI M, YANO M, TEJIMA Y, et al. High-pressure transitions of diopside and wollastonite: phase equilibria and thermochemistry of CaMgSi2O6, CaSiO3 and CaSi2O5-CaTiSiO5 system [J]. Physics of the Earth and Planetary Interiors, 2004, 143/144: 145–156. doi: 10.1016/j.pepi.2003.08.008 [22] KATZ L, WARD R. Structure relations in mixed metal oxides [J]. Inorganic Chemistry, 1964, 3(2): 205–211. doi: 10.1021/ic50012a013 [23] CHENG J G, ALONSO J A, SUARD E, et al. A new perovskite polytype in the high-pressure sequence of BaIrO3 [J]. Journal of the American Chemical Society, 2009, 131(21): 7461–7469. doi: 10.1021/ja901829e [24] SASAKI S, PREWITT C T, LIEBERMANN R C. The crystal structure of CaGeO3 perovskite and the crystal chemistry of the GdFeO3-type perovskites [J]. American Mineralogist, 1983, 68(11/12): 1189–1198. [25] SHANNON R D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides [J]. Acta Crystallographica Section A: Foundations and Advances, 1976, 32(5): 751–767. doi: 10.1107/S0567739476001551 [26] MAO H K, XU J, BELL P M. Calibration of the ruby pressure gauge to 800 kbar under quasi-hydrostatic conditions [J]. Journal of Geophysical Research: Solid Earth, 1986, 91(B5): 4673–4676. doi: 10.1029/JB091iB05p04673 [27] HAMMERSLEY A P, SVENSSON S O, HANFLAND M, et al. Two-dimensional detector software: from real detector to idealised image or two-theta scan [J]. High Pressure Research, 1996, 14(4/6): 235–248. doi: 10.1080/08957959608201408 [28] HOLLAND T J B, REDFERN S A T. Unit cell refinement from powder diffraction data: the use of regression diagnostics [J]. Mineralogical Magazine, 1997, 61(404): 65–77. doi: 10.1180/MINMAG.1997.061.404.07 [29] TOBY B H, VON DREELE R B. GSAS-Ⅱ: the genesis of a modern open-source all purpose crystallography software package [J]. Journal of Applied Crystallography, 2013, 46(2): 544–549. doi: 10.1107/S0021889813003531 [30] GONZALEZ-PLATAS J, ALVARO M, NESTOLA F, et al. EosFit7-GUI: a new graphical user interface for equation of state calculations, analyses and teaching [J]. Journal of Applied Crystallography, 2016, 49(4): 1377–1382. doi: 10.1107/S1600576716008050 [31] 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 [32] LIN C C, SHEN P Y. Pressure-induced metastable phase transformations of calcium metasilicate (CaSiO3): a Raman spectroscopic study [J]. Materials Chemistry and Physics, 2016, 182: 508–519. doi: 10.1016/j.matchemphys.2016.07.065 [33] KRONBO C H, MENESCARDI F, CERESOLI D, et al. High pressure structure studies of three SrGeO3 polymorphs: amorphization under pressure [J]. Journal of Alloys and Compounds, 2021, 855: 157419. doi: 10.1016/j.jallcom.2020.157419 -