An official website of the United States government
We explored the B–C–O system at pressures in the range 0–50 GPa by ab initio variable-composition evolutionary simulations in the hope of discovering new stable superhard materials. A new tetragonal thermodynamically stable phase B 4 CO 4 , space group I 4̄, and two low-enthalpy metastable compounds (B 6 C 2 O 5 , B 2 CO 2 ) have been discovered. Computed phonons and elastic constants show that these structures are dynamically and mechanically stable both at high pressure and zero pressure. B 4 CO 4 is thermodynamically stable at pressures above 23 GPa, but should remain metastable under ambient conditions. Its computed hardness is about 38–41 GPa, which suggests that B 4 CO 4 is potentially superhard.
more »
« less
Novel magnesium borides and their superconductivity
With the motivation of searching for new superconductors in the Mg–B system, we performed ab initio evolutionary searches for all the stable compounds in this binary system in the pressure range of 0–200 GPa. We found previously unknown, yet thermodynamically stable, compositions MgB 3 and Mg 3 B 10 . Experimentally known MgB 2 is stable in the entire pressure range 0–200 GPa, while MgB 7 and MgB 12 are stable at pressures below 90 GPa and 35 GPa, respectively. We predict a reentrant behavior for MgB 4 , which becomes unstable against decomposition into MgB 2 and MgB 7 at 4 GPa and then becomes stable above 61 GPa. We find ubiquity of phases with boron sandwich structures analogous to the AlB 2 -type structure. However, with the exception of MgB 2 , all other magnesium borides have low electron–phonon coupling constants λ of 0.32–0.39 and are predicted to have T c below 3 K.
more »
« less
- Award ID(s):
- 1723160
- PAR ID:
- 10094674
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 19
- Issue:
- 22
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 14486 to 14494
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The high-pressure structure and stability of the calcic amphibole tremolite (Ca2Mg5Si8O22(OH)2) was investigated to ~40 GPa at 300 K by single-crystal X-ray diffraction using synchrotron radiation. C2/m symmetry tremolite displays a broader metastability range than previously studied clinoamphiboles, exhibiting no first-order phase transition up to 40 GPa. Axial parameter ratios a/b and a/c, in conjunction with finite strain versus normalized pressure trends, indicate that changes in compressional behavior occur at pressures of ~5 and ~20 GPa. An analysis of the finite strain trends, using third-order Birch-Murnaghan equations of state, resulted in bulk moduli (𝐾) of 72(7), 77(2), and 61(1) GPa for the compressional regimes from 0-5 GPa (regime I), 5-20 GPa (II), and above 20 GPa (III), respectively, and accompanying pressure-derivatives of the bulk moduli (𝐾′) of 8.6(42), 6.0(3), and 10.0(2). The results are consistent with first-principle theoretical calculations of tremolite elasticity. The axial compressibility ratios of tremolite, determined as 𝛽a : 𝛽b : 𝛽c = 2.22:1.0:0.78 (regime I), 2.12:1.0:0.96 (II), and 1.03:1.0:0.75 (III), demonstrate a substantial reduction of the compressional anisotropy of tremolite at high pressures, which is a notable contrast with the increasingly anisotropic compressibility observed in the high-pressure polymorphs of the clinoamphibole grunerite. The shift in compression-regime at 5 GPa (I-II) transition is ascribed to stiffening along the crystallographic a-axis corresponding to closure of the vacant A-site in the structure, and a shift in the topology of the a-oriented surfaces of the structural I-beam from concave to convex. The II-III regime shift at 20 GPa corresponds to an increasing rate of compaction of the Ca-polyhedra and increased distortion of the Mg-octahedral sites, processes which dictate compaction in both high-pressure compression-regimes. Bond-valence analyses of the tremolite structure under pressure show dramatic overbonding of the Ca-cations (75% at 30 GPa), with significant Mg-cation overbonding as well (40%). These imply that tremolite’s notable metastability range hinges on the calcium cation’s bonding environment. The 8-fold coordinated Ca-polyhedron accommodates significant compaction under pressure, while the geometry of the Ca-O polyhedron becomes increasingly regular and inhibits the reorientation of the tetrahedral chains that generate phase transitions observed in other clinoamphiboles. Peak/background ratio of diffraction data collected above 40 GPa and our equation of state determination of bulk moduli and compressibilities of tremolite in regime III, in concert with the results of our previous Raman study, suggest that C2/m tremolite may be approaching the limit of its metastability above 40 GPa. Our results have relevance for both the metastable compaction of tremolite during impact events, and for possible metastable persistence of tremolite within cold subduction zones within the Earth.more » « less
-
Abstract Iron nitrides are possible constituents of the cores of Earth and other terrestrial planets. Pressure‐induced magnetic changes in iron nitrides and effects on compressibility remain poorly understood. Here we report synchrotron X‐ray emission spectroscopy (XES) and X‐ray diffraction (XRD) results for ε‐Fe7N3and γ′‐Fe4N up to 60 GPa at 300 K. The XES spectra reveal completion of high‐ to low‐spin transition in ε‐Fe7N3and γ′‐Fe4N at 43 and 34 GPa, respectively. The completion of the spin transition induces stiffening in bulk modulus of ε‐Fe7N3by 22% at ~40 GPa, but has no resolvable effect on the compression behavior of γ′‐Fe4N. Fitting pressure‐volume data to the Birch‐Murnaghan equation of state yieldsV0 = 83.29 ± 0.03 (Å3),K0 = 232 ± 9 GPa,K0′ = 4.1 ± 0.5 for nonmagnetic ε‐Fe7N3above the spin transition completion pressure, andV0 = 54.82 ± 0.02 (Å3),K0 = 152 ± 2 GPa,K0′ = 4.0 ± 0.1 for γ′‐Fe4N over the studied pressure range. By reexamining evidence for spin transition and effects on compressibility of other candidate components of terrestrial planet cores, Fe3S, Fe3P, Fe7C3, and Fe3C based on previous XES and XRD measurements, we located the completion of high‐ to low‐spin transition at ~67, 38, 50, and 30 GPa at 300 K, respectively. The completion of spin transitions of Fe3S, Fe3P, and Fe3C induces elastic stiffening, whereas that of Fe7C3induces elastic softening. Changes in compressibility at completion of spin transitions in iron‐light element alloys may influence the properties of Earth's and planetary cores.more » « less
-
Abstract The response of forsterite, Mg2SiO4, under dynamic compression is of fundamental importance for understanding its phase transformations and high‐pressure behavior. Here, we have carried out an in situ X‐ray diffraction study of laser‐shocked polycrystalline and single‐crystal forsterite (a‐,b‐, andc‐orientations) from 19 to 122 GPa using the Matter in Extreme Conditions end‐station of the Linac Coherent Light Source. Under laser‐based shock loading, forsterite does not transform to the high‐pressure equilibrium assemblage of MgSiO3bridgmanite and MgO periclase, as has been suggested previously. Instead, we observe forsterite and forsterite III, a metastable polymorph of Mg2SiO4, coexisting in a mixed‐phase region from 33 to 75 GPa for both polycrystalline and single‐crystal samples. Densities inferred from X‐ray diffraction data are consistent with earlier gas‐gun shock data. At higher stress, the response is sample‐dependent. Polycrystalline samples undergo amorphization above 79 GPa. For [010]‐ and [001]‐oriented crystals, a mixture of crystalline and amorphous material is observed to 108 GPa, whereas the [100]‐oriented forsterite adopts an unknown phase at 122 GPa. The first two sharp diffraction peaks of amorphous Mg2SiO4show a similar trend with compression as those observed for MgSiO3in both recent static‐ and laser‐driven shock experiments. Upon release to ambient pressure, all samples retain or revert to forsterite with evidence for amorphous material also present in some cases. This study demonstrates the utility of femtosecond free‐electron laser X‐ray sources for probing the temporal evolution of high‐pressure silicate structures through the nanosecond‐scale events of shock compression and release.more » « less
-
null (Ed.)Diffusional isotope fractionation occurs in geochemical processes (such as magma mixing, bubble growth, and crystal growth), even at magmatic temperatures. Isotopic mass dependence of diffusion is commonly expressed as Di Dj ¼ mj mi b , where Di and Dj are diffusion coefficients of two isotopes whose masses are mi and mj. How the dimensionless empirical parameter b depends on temperature, pressure, and composition remains poorly constrained. Here, we conducted a series of first-principles molecular dynamics simulations to evaluate the b factor of Mg isotopes in MgSiO3 and Mg2SiO4 melts using pseudo-isotope method. In particular, we considered interactions between Mg isotopes by simultaneously putting pseudo-mass and normalmass Mg atoms in a simulation supercell. The calculated b for Mg isotopes decreases linearly with decreasing temperature at zero pressure, from 0:158 0:004 at 4000 K to 0:121 0:017 at 2200 K for MgSiO3 melt and from 0:150 0:004 at 4000 K to 0:101 0:012 at 2200 K for Mg2SiO4 melt. Moreover, our simulations of compressed Mg2SiO4 melt along the 3000 K isotherm show that the b value decreases linearly from 0:130 0:006 at 0 GPa to 0:060 0:011 at 17 GPa. Based on our diffusivity results, the empirically established positive correlation between b and solvent-normalized diffusivity (Di/DSi) seems to be applicable only at constant temperatures or in narrow temperature ranges. Analysis of atomistic mechanisms suggests that the calculated b values are inversely correlated with force constants of Mg at a given temperature or pressure. Good agreement between our first principles results with available experimental data suggests that interactions between isotopes of major elements must be considered in calculating b for major elements in silicate melts. Also, we discuss diffusion-controlled crystal growth by considering our calculated b values.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/c7cp00840f
