Compressibility and equation of state of beryl (Be3Al2Si6O18) by using a diamond anvil cell and in situ synchrotron X-ray diffraction

High-pressure single-crystal synchrotron X-ray diffraction was carried out on a single crystal of natural beryl compressed in a diamond anvil cell. The pressure-volume (P-V) data from room pressure to 9.51 GPa were fitted by a third-order Birch-Murnaghan equation of state (BM-EoS) and resulted in unit-cell volume V-0 = 675.5 +/- 0.1 angstrom(3), isothermal bulk modulus K-0 = 180 +/- 2 GPa, and its pressure derivative K-0' = 4.2 +/- 0.5. We also calculated V-0 = 675.5 +/- 0.1 angstrom(3) and K-0 = 181 +/- 1GPa with fixed K-0' at 4.0 and then obtained the axial moduli for a (K-a0)-axis and c (K-c0)-axis of 209 +/- 1 and 141 +/- 2 GPa by "linearized" BM-EoS approach. The axial compressibilities of a-axis and c-axis are beta(a) = 1.59 x 10(-3) GPa(-1) and beta(c) = 2.36 x 10(-3) GPa(-1) with an anisotropic ratio of beta(a):beta(c) = 0.67:1.00. On the other hand, the pressure-volume-temperature (P-V-T) EoS of the natural beryl has also been measured at temperatures up to 750 K and at pressures up to 16.81 GPa, using diamond anvil cell in conjunction with in situ synchrotron angle-dispersive powder X-ray diffraction. The P-V data at room temperature and at a pressure range of 0.0001-15.84 GPa were then analyzed by third-order BMEoS and yielded V-0 = 675.3 +/- 0.1 angstrom(3), K-0 = 180 +/- 2 GPa, K-0' = 4.2 +/- 0.3. With K-0' fixed to 4.0, we also obtained V-0 = 675.2 +/- 0.1 angstrom(3) and K-0 = 182 +/- 1 GPa. Consequently, we fitted the P-V-T data with high-temperature BM-EoS approach using the resultant K-0' (4.2) from room-temperature BM-EoS and then obtained the thermoelastic parameters of V-0 = 675.3 +/- 0.2 angstrom(3), K-0 = 180 +/- 1 GPa, temperature derivative of the bulk modulus (partial derivative K/partial derivative T)(P) = -0.01 7 +/- 0.004 GPa K-1, and thermal expansion coefficient at ambient conditions alpha(0) = (2.82 +/- 0.74) x 10(-6) K-1. Present results were also compared with previous studies for beryl. From the comparison of these fittings, we propose to constrain K-0 = 180 GPa and K-0' = 4.2 for beryl. And we also observed that beryl exhibits anisotropic thermal expansion at relatively low temperatures, which is very consistent with previous studies. Furthermore, no phase transition was observed in the entire pressure and temperature range (up to 16.84 GPa and 750 K) of this study for the natural beryl.