Miaoting
Xu
a,
Cuicui
Wang
a,
Benjamin J.
Morgan
b and
Lee A.
Burton
*a
aInternational Centre for Quantum and Molecular Structures, Department of Physics, Shanghai University, Shanghai 200444, China. E-mail: leeburton@shu.edu.cn
bDepartment of Chemistry, University of Bath, Claverton Down, BA2 7AY, UK
First published on 31st March 2022
The electride Sr3CrN3 has a one-dimensional channel of electron density, which is a rare feature that offers great potential for fast ion conduction. Using density functional theory, we find that Sr3CrN3 is an excellent hydride conductor within this channel, with a diffusion barrier as low as 0.30 eV and an estimated diffusion coefficient of 5.37 × 10−8 cm2 s−1. This diffusion barrier is lower than those reported for the best hydride conductors to date. We also show the most-stable amount of hydride in the host material under standard conditions and the corresponding change in electronic structure from metal to wide-gap insulator. Our results highlight the potential offered by 1D electride materials for ion-transport applications such as energy storage or gas separation.
Solid-state ion-conducting materials have been long-studied for their potential use in energy storage systems.4 Recent years have seen particular advances in the development of highly-conducting lithium-ion solid electrolytes with potential application in all-solid-state lithium-ion batteries.5–8 The limited global availability of lithium, however, means that lithium-ion batteries are expected to meet only part of the projected future energy-storage needs.9 This motivates the continuing search for alternate materials, including solid-state electrolytes, that could be used in non-lithium energy-storage devices.
Hydride ions have small ionic radii, large electronic polarizability and a high standard redox potential for H2/H− (−2.3V).10 These characteristics make H− ions promising for applications in next-generation electrochemical energy storage with high voltage and high energy density. Furthermore, hydrogen is highly abundant, globally available and low cost. To date, materials with high ion conduction for H− have been reported but these were only achieved at relatively high temperatures.11 Thus, there is still ample scope for improvement in this area.
Electrides are class of rare ionic compounds that possess free electrons localized within cavities in the host structure, with these electrons acting as anions. To date, only a few inorganic electrides have been experimentally identified,12 but they are typically classified according to dimensionality of their free electron density. For example, if the excess electron density occupies cavities or pores the electrides are classified as 0D (e.g., Ca12Al14O3213), if the density is in a one-dimensional channel they are 1D (e.g. Y5Si314) or if the density is continuous in a plane they are 2D (e.g. Ca2N15). The anionic electrons of electrides have been shown to effectively interact with external hydrogen, leading to excellent hydride absorption and desorption properties,16,17 even in 1D electrides.14
In this work we consider Sr3CrN3, a recently identified 1D electride, as a hydride conductor.18 Sr3CrN3 can be described with the formal oxidation states Sr32+Cr4+N33−:e− and the excess electrons aggregating in one-dimensional channels (see Fig. 1). The assignment of Cr(IV) was demonstrated with previous DFT calculations and subsequently confirmed by experimental analysis.18,19 Here, we explore the optimal diffusion path of hydride ions in the channel using the climbing-image nudged-elastic-band (CI-NEB) method. We also investigate the hydrogen capacity by calculating insertion energies of hydride ions in the structure. Finally, we report the electronic nature of the host material before and after hydride intercalation, showing a shift from metallic to insulating behaviour. Overall, we find that Sr3CrN3 has a low capacity for hydrogen under standard conditions but the diffusion barrier of 0.3 eV within the 1D channels is significantly lower than barriers reported for best-in-class hydride ion conductors previously e.g. 0.52 eV for barium hydride.11
For the Sr3CrN3Hx (0 < x < 3) calculations, an energy cut-off of 520 eV was employed. Any H atoms were added to the system as neutral species so no charge image or charge balancing considerations are required and all computations were spin polarized. Magnetic ions were initialised in a high-spin ferromagnetic configuration and then allowed to relax during each calculation. To sample k-space we used a 6 × 6 × 9 Monkhorst–Pack set of k-points with the tetrahedron method.26 Electronic convergence criterion was set to 1 × 10−6 eV and ionic convergence criterion was set to 1 × 10−5 eV Å−1 in all cases. The H2 gas was calculated in a unit cell of 10 Å3 with a gamma centred single point k-grid.
To investigate potential diffusion pathways for mobile H-we performed a series of Climbing Image Nudged Elastic Band (CI-NEB) calculations with fixed mid-points, as implemented in VASP combined with the VTST-Tools27 by Henkelman et al.28
Finally the screened–exchange hybrid density functional of Heyd–Scuseria–Ernzerhof (HSE06)29 was used for calculations of electronic properties. This method is known to correct the underestimation of band gap and over delocalisation of conventional DFT.30 For plotting the band structures we use Pymatgen's “electronic_structure” module BSPlotter.31
Lattice constant | a (Å) | b (Å) | c (Å) | α (°) | β (°) | γ (°) |
---|---|---|---|---|---|---|
DFT relaxed | 7.84 | 7.84 | 5.24 | 90.0 | 90.0 | 120.0 |
Experiment | 7.72 | 7.72 | 5.25 | 90.0 | 90.0 | 120.0 |
In Sr3CrN3, Sr and N ions occupy Wyckoff 6h sites and Cr occupies 2c sites. The remaining high-symmetry Wyckoff sites are vacant in the stoichiometric material and are therefore available as potential sites to accommodate anionic hydrogen.
To analyse the ability of Sr3CrN3 to accommodate intercalated hydrogen, we consider all possible H− configurations for x(H) = 1 to 6, where x(H) = 6 corresponds to all available Wyckoff sites being occupied by H. The intercalation energy as a function of x(H) is calculated as
(1) |
The energy change (ΔE) for all 25 total combinations of hydride position are listed in Table S1 (ESI†) and plotted in Fig. 2, although symmetry equivalence means fewer than 25 data points are visible in the plot. These results show that under standard conditions one intercalated hydride per formula unit has the lowest energy overall compared to the pure electride. There are two formula units of Sr3CrN3 in a unit cell, each one providing an excess electron that makes the material an electride. Thus, the addition of two hydrogen atoms to the unit cell creates the material Sr3CrN3H, which agrees with experimental reports on the related material Ba3CrN3H.33 The positions of the hydride ions corresponding to the lowest energy are on the 2b Wyckoff positions with the fractional unit cell coordinates (0, 0, 0) and (0, 0, 0.5). The continued addition of hydrides is still favourable up to the addition of the 3rd hydride ion per formula unit, at which point the change in energy becomes positive. This indicates that increasing chemical potential of H would be needed for additional uptake of H to proceed.
Fig. 2 The energy change as a function of hydrogen atoms added to Sr3CrN3 relative to the phase pure material. The lowest energies are connected to form a convex hull. The full data for this figure, corresponding to 25 total calculations, can be found in the ESI.† |
McColm et al.34 previously reported promise for electrides in the area of hydrogen storage. However, we find that under standard conditions at most 2.5 hydride ions per formula unit can be placed in the Sr3CrN3. Given that the volume of the unit cell is 282.68 Å3 we calculate the corresponding volumetric H2 density (V(H2)) and relative gravimetric H2 density (ω(H2)) as 29.61 kg m−3 and 0.70% respectively. These are calculated using eqn (2) and (3), below. Compared with other hydrogen storage materials, the hydrogen storage density of Sr3CrN3is low—even lower than pressurised H2 gas—owing to the elemental composition and structure of this material.35 Firstly, there are two rather heavy metal elements present, which makes the relative molecular mass large. Secondly, Sr3CrN3 has a 1D cavity channel, so the available volume is relatively low compared to other open-pore materials.
(2) |
(3) |
Fig. 3 The lowest energy hydride ion pathway in the Sr3CrN3 structure. The green spheres are Sr, the blue spheres are Cr, the light grey spheres are N, and the red spheres are H. |
The CI-NEB method generates a direct linear path connecting the start and end points with a very low energy barrier of 0.30 eV as shown in Fig. 4. Even if the atoms of the host material remain fixed at the phase-pure positions, the energy is still found to be remarkably low at 0.35 eV. Our calculated energy barrier is smaller than those reported in the literature for the fastest reported hydride ion conducting materials: 0.52 eV for barium hydride,11 or 1.2 eV for oxygen substituted lanthanum hydride.10 Furthermore, proton migration is usually associated with an activation energy of higher than 0.5 eV in oxide materials,37 despite protons having a smaller ionic radii than hydrides. As a result, these are among the most favourable properties for ionic diffusion reported previously.
Fig. 4 The energy barriers from the CI-NEB method for the in-channel migration corresponding between (0, 0, 0) and (0, 0, 0.5) internal coordinates as shown in Fig. 3. |
We sample increasing displacement from the linear path within the channel to confirm the route calculated here is the ground state. We also analysed alternate migration routes outside of the one-dimensional channel. When we select (0.10, 0.25, 0.25) as midpoint, the migration barrier is 2.40 eV; for (0.25, 0, 0.25) the barrier is 2.64 eV and for (0.15, 0.10, 0.23), the barrier is 3.94 eV. Finally, we also sampled the paths along the other 2 unit-cell axes with a direction (0, 0, 0 → 0, 0, 0.5) giving a barrier of 5.04 eV and b direction (0, 0, 0 → 0, 0.5, 0) giving a barrier of 4.92 eV. Thus, we conclude that migration within the electride channel is highly favourable relative to alternate routes through the material. The large difference between energy barriers means that this material will be highly directional in its ionic migration, suggesting that so-called superionic fast conduction may be possible.38
Given that the NEB barrier of 0.3 eV is still much greater than KBT at 298 K the H diffusion process falls into a so-called “independent hopping” regime, meaning that transition-state theory can be applied to calculate diffusivity. This allows us to effectively account for temperature.39 The full equation for calculating diffusivity (D) is shown below in eqn (4):
(4) |
Furthermore, ionic conductivity can be estimated using the Nernst–Einstein approximation from the diffusion coefficient, as shown in eqn (5).
(5) |
Using the relaxed volume of the unit cell we get a V of 322.1 cubic angstroms. For T we can use 300 K, N as 2 ions per unit cell and Z of 1. Finally, for D0 we use = 6.3 × 10−3 cm2 s−1 as before. This gives an estimated conductivity of 242.4 Siemens per centimetre. Finally, the ionic transference number will be exactly 1, because the only mobile species in this material is the H− ion.42
Fig. 5 The band structure for (top) Sr3CrN3 showing metallic behaviour and (bottom) Sr3CrN3H showing a relatively wide band gap. The solid lines are spin up and dashed lines are spin down. |
We also perform HSE06 calculations on supercell structures with partial H ion occupancy on the 2b Wyckoff sites i.e. the location of the anionic electrons. The data for these calculations is shown in full in Table S2 (ESI†). When not all anionic electrons are used in the reduction of the H atoms there is still an opening of a band gap, although smaller in magnitude than for the full case. We believe that this is because the delocalised anionic electron density in the channel is disrupted by the presence of the hydride, even in relatively dilute amounts, leading to semiconductor rather than metallic behaviour.
Druffel et al.44 have suggested that electrides might find use in high-capacity electrodes that realise electron–anion reversible exchange at room temperature. In the case of Sr3CrN3 we predict no significant structural reorganisation and very small changes in unit cell volume upon hydrogen intercalation. Our analysis of the electronic nature of the material precludes such an application as the intercalated Sr3CrN3H is strongly insulating, whereas an electrode is typically sought to be metallic in nature. What's more the flat bands mean a high charge effective mass and low conductivity, which is again unfavourable for application as an electrode. The hydride capacity is also relative low for such an application. However, the wide-band gap in the presence of hydride ions suggests that the material could be employed as an electrolyte in a fuel cell configuration, wherein a constant flow of hydride ions could maintain electrically insulating behaviour, which is what forces the electrons around the external circuit to produce electricity available for work. Finally, the band structure shows that the spin-up and spin-down states are virtually degenerate, which indicates no intrinsic ferromagnetism for either the hydride or phase pure material.
On the other hand, we find Sr3CrN3H to be an excellent conductor of hydride ions. Within the 1D channel hydrogen ions show a migration barrier of 0.30 eV, which is significantly lower than for even some of the best performing hydride-ion conductors reported to date, e.g., 0.52 eV for barium hydride,11 or 1.2–1.3 eV for oxygen substituted lanthanum hydride.10 While it is tempting to attribute the low diffusion barrier to the anionic electrons in a one-dimensional channel in the material facilitating conduction, the electronic structure indicates that the excess electrons are used in reducing the H to H−; becoming a wide-band gap ionic material in the process. Thus, we assume that the partially-coordinated Sr2+ cations surrounding the 1D channel provide a uniquely favourable environment to unimpeded anion conduction. We foresee ample scope for further study of ionic conduction in this material and other electrides including a broader range of intercalated ions; molecular dynamics or phonon calculations to more extensively explore temperature effects and surface or kinetic considerations to better understand absorption/desorption processes. Overall, we find that Sr3CrN3H exhibits exceptional ionic transport properties and could be particularly useful for energy storage devices, catalysis or gas separation applications.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1tc05850a |
This journal is © The Royal Society of Chemistry 2022 |