Tina Chen1,2
Gopalakrishnan Sai Gautam
Pieremanuele Canepa- 1Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, CA, United States
- 2Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
- 3Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ, United States
- 4Department of Materials Science and Engineering, National University of Singapore, Singapore, Singapore
Mg batteries utilizing a Mg metal anode with a high-voltage intercalation cathode define a potential pathway toward energy storage with high energy density. However, the making of Mg batteries is plagued by the instability of existing electrolytes against the Mg-metal anode and high-voltage cathode materials. One viable solution to this problem is the identification of protective coating materials that could effectively separate the distinct chemistries of the metal-anode and the cathode materials from the electrolyte. Using first-principles calculations we mapped the electrochemical stability windows for non-redox-active Mg binary and ternary compounds in order to identify potential coating materials for Mg batteries. Our results identify Mg-halides and Mg(BH4)2 as promising anode coating materials based on their significant reductive stability. On the cathode side, we single out MgF2, Mg(PO3)2, and MgP4O11 as effective passivating agents.
Introduction
Multivalent batteries, such as those based on Mg, present a potential alternative to Li-ion batteries, particularly in terms of increased energy density (Canepa et al., 2017a). Mg batteries are able to use Mg metal as an anode at reasonable current densities (<0.5 mA/cm2) (Yoo et al., 2013), which in combination with the higher oxidation state of Mg (+2 rather than Li's +1) can provide a significant increase in the energy density of Mg batteries compared to Li-ion batteries. So far, prototypes of Mg batteries have utilized electrolytes, such as MgCl2 with AlCl3, Mg(ClO4)2, Mg(NO3)2, Mg(TFSI)2, and more complex molecules dissolved in acetonitrile, THF, or glymes-based solvents, in combination with Mg metal as the anode and a low voltage sulfide cathode (MgxMo6S8 and MgxTiS2) (Aurbach et al., 2000; Cohen et al., 2000; Pour et al., 2011; Mohtadi et al., 2012; Muldoon et al., 2012, 2014; Yoo et al., 2013; Carter et al., 2014; Doe et al., 2014; Canepa et al., 2015a; Tutusaus et al., 2015, 2016; Sun et al., 2016; Hahn et al., 2018).
Typical Mg electrolytes have significantly narrower electrochemical stability windows (~1.5–3.0 V vs. Mg) (Lipson et al., 2016) compared to what is available in the Li-ion battery space (~1.5–5 V vs. Li) (Marom et al., 2011). Indeed, most electrolytes, including the solvents used in commercial Li electrolytes (e.g., PC and DMC) (Goodenough and Kim, 2009), have poor reductive stability (i.e., cathodic stability) and tend to decompose at the Mg metal anode (Lu et al., 1999; Muldoon et al., 2012). In addition, the utilization of high-voltage cathodes (e.g., oxides) is greatly impeded by the limited oxidative stability (i.e., anodic stability) of Mg electrolytes (Rosenberg and Nicolau, 1964; Cohen et al., 2000; Pour et al., 2011; Mohtadi et al., 2012; Muldoon et al., 2012, 2014; Yoo et al., 2013; Carter et al., 2014; Doe et al., 2014; Canepa et al., 2015a; Liu et al., 2015; Tutusaus et al., 2015, 2016; Chen et al., 2017; Hahn et al., 2018). Thus, the reactivity of the electrolyte against both Mg-anode and a high-voltage cathode results in electrolyte decomposition, often producing a passivating layer primarily containing a binary Mg-salt, such as MgO (and Mg(OH)2 if moisture is present) (Gofer et al., 2003; Ling et al., 2015; Ling and Zhang, 2017; Hannah et al., 2018). The presence of MgO greatly inhibits Mg2+ transport (Canepa et al., 2017b) and eventually the ability of the battery to store energy reversibly (Levi et al., 2009). Further work is still being done to develop Mg electrolytes that can reversibly strip and deposit Mg at the anode and cathode (Muldoon et al., 2012, 2014; Canepa et al., 2015b). For example, a class of carboranes has recently been proposed as promising electrolytes, stable against Mg metal and high voltage cathodes (up to 4.6 V vs. Mg) (Hahn et al., 2018). However, more work is required to elucidate the mechanisms of reversible Mg transfer at the cathode and develop strategies to mitigate electrolyte decomposition (Shao et al., 2013; Keyzer et al., 2016).
In analogous Li-systems, several approaches have been utilized to address the safety and electrochemical stability limitations of typical Li electrolytes (Aurbach et al., 2004; Guerfi et al., 2010). For example, solid electrolytes have been shown to be safer compared to typical solvent-based electrolytes, which may experience thermal runaway issues (Kamaya et al., 2011; Masquelier, 2011; Bachman et al., 2015; Kato et al., 2016). Another ongoing field of research is the application of protective coating layers to shield one or both electrodes from an incompatible electrolyte, while providing sufficient ionic mobility and preferably low electronic conductivity. Indeed, the solid electrolyte interphase (SEI) that forms at the graphitic anode-electrolyte interface is a good example of a protective layer with sufficient Li mobility that enables the reversible operation of Li-ion batteries (Verma et al., 2010). Therefore, similar solutions can be envisioned for Mg-batteries as well. To accomplish this goal, we searched for materials that can act as either protective coatings or even solid electrolytes by analyzing the electrochemical stability of various Mg-containing compounds.
Using a combination of density functional theory (DFT) calculations and thermodynamics, we assessed the electrochemical stability of various Mg-binary and ternary compounds, which may form as a result of electrolyte decomposition at either the Mg-metal anode or a high-voltage cathode. Specifically, we considered all Mg binaries and ternaries that do not contain redox-active metal ions (except Ti4+) and that are known to be electronic insulators. The choice of Mg compounds is also motivated by the highly reducing conditions that appear when in contact with Mg metal. For example, Li binaries and ternaries, such as Li3N, Li3P, LiH, Li2S, Li2O, and LiCl, tend to form (and be stable) at the Li electrolyte-anode interface in Li-ion batteries (Richards et al., 2015).
By calculating the electrochemical stability windows of candidate compounds, we identified their oxidative and reductive voltages. Our findings provide general guidelines for developing, via either in situ or ex situ deposition techniques, protective coating materials that are compatible with the anode or the cathode or both. Provided good bulk Mg2+ mobility exists (Sai Gautam et al., 2017), some of these materials may be investigated as protective coating materials or even solid electrolytes.
Methodology
The set of elements from which we evaluated Mg binaries and ternaries is shown in Figure 1, with Mg colored in red and the other elements colored based on their respective group numbers (a complete list of all Mg-binaries and ternaries investigated is provided in Tables S1, S2). In addition to the highlighted elements, we considered borohydrides, niobates, titanates, titanium phosphates, and zirconium phosphates which have been reported to be promising coating materials in Li-ion batteries (Richards et al., 2015). Also, we included Mg-(Sc/In)-(S/Se) compounds since they have been explored as potential Mg solid-electrolyte materials in prior studies (Canepa et al., 2017b,c), apart from Mg-(Al/Ga/In)-(O/S/Se).
Figure 1. Periodic table highlighting the non-transition-metal elements that form binary (and ternary) compounds with Mg (red), including triels (Group IIIA, green), tetrels (Group IVA, light blue), pnictogens (Group VA, yellow), chalcogens (Group VIA, gray), halogens (Group VIIA, orange), and other elements (Hydrogen, purple). We considered all Mg-X binaries and stable Mg-X-Y ternaries, where X and Y are highlighted elements, with the exception of the Mg-X-H chemical space where only Mg-B-H compounds were considered. In addition, we evaluated some compounds containing a non-Mg metal, such as Sc, Ti, Nb, Zr, Al, Ga, and In, because either they are commonly used as coating materials in Li-ion batteries or have been considered as Mg ionic conductors in prior studies.
The electrochemical stability windows of each compound are calculated using the approach developed by (Richards et al., 2015) by constructing the corresponding grand potential (ϕ) phase diagram by means of the pymatgen library (Jain et al., 2011; Ong et al., 2013), where ϕ is defined as:
For all μMg, we constructed the convex hull in the grand potential composition-space and identified compounds that are stable at each μMg. The Mg chemical potential μMg relates directly to the voltage vs. Mg/Mg2+ via (Equation 2):
where F is the Faraday constant, z is the number of electrons transferred (z = 2 for Mg) and μMg is referenced to the energy of Mg metal. The internal energy of each compound [E in Equation (1)], in the relevant chemical space, was either obtained from the Materials Project (Jain et al., 2011, 2013) database or calculated directly using DFT (Kohn and Sham, 1965; Hohenberg and Kohn, 1973) (see Input parameters for DFT calculations in SI for more details on the calculation parameters used). For each compound, we utilized the atomic coordinates reported in the Inorganic Crystal Structure Database (ICSD) (Bergerhoff and Brown, 1987) as initial guesses during our DFT structure relaxation. For Mg0.5Zr2(PO4)3 and Mg0.5Ti2(PO4)3, which are disordered structures in the ICSD database, we enumerated possible configurations within the respective unit cell (Hart and Forcade, 2008, 2009; Hart et al., 2012; Ong et al., 2013) and included the lowest energy configuration.
Results and Discussion
Electrochemical Stability Windows of Mg-Binaries
Figure 2 shows the voltage windows of all Mg-X binaries considered, where the compounds are grouped by the anion column (Figure 1) and sorted within each group by increasing electronegativity.
Figure 2. Electrochemical stability windows of non-metallic Mg-binaries, indicating the voltages (vs. Mg metal) at which the compound is stable against decomposition. Compounds that are not stable at any voltage, such as the Mg-carbides, are not shown. Compounds are grouped by the anion column and compounds within each group are ordered by increasing electronegativity. For systems with multiple compositions, compounds are ordered by decreasing ratio of Mg to anion. The number at the end of each bar indicates the width of the voltage window. The MgxCrO2 spinel is shown above the plot at its calculated average voltage (~3.6 V) for reference.
To form a binary system with multiple stable compounds (e.g., Mg-B), we ordered the compounds according to a decreasing ratio of Mg to anion (Mg:B). Only binaries that were thermodynamically stable (i.e., with negative formation energy at 0 K) are shown. Unstable compounds have been removed from Figure 2 because they will not be stable at any μMg. For example, MgC2 has a formation energy of 173 meV/atom at 0 K. The left and right ends of the bar for each compound indicate the lower and upper voltage limits, corresponding to the reductive (cathodic) and oxidative (anodic) stabilities, respectively. Lower reductive stabilities and higher oxidative stabilities imply better resistance against reduction and oxidation, respectively. Thus, the width of the bar (text annotation to the right of each bar in Figure 2) for a given compound signifies its electrochemical stability window. The zero on the voltage axis is referenced to bulk Mg metal (i.e., V vs. Mg/Mg2+). Higher voltage values mimic the open circuit voltages of cathode materials, such as Chevrel-Mo6S8 (~1.1 V) (Aurbach et al., 2000), layered-V2O5 (~3.3 V) (Sai Gautam et al., 2015), or MgxCrO2 (~3.6 V) (Chen et al., 2017).
Significantly, all of the Mg-halides, Mg-chalcogenides, and Mg-pnictides (except MgP4) are stable at 0 V vs. Mg/Mg2+ and thus stable against Mg metal. Among the Mg-triels and Mg-tetrels, only MgB2, Mg2Ge, and Mg2Si are stable vs. Mg metal. However, the widths of the stability windows of MgB2, Mg2Ge, Mg2Si are small (<0.1 V), and thus Mg-triels and Mg-tetrels do not appear to be viable coating materials against typical electrolytes. The poor stability windows of MgB2, and Mg2Ge, Mg2Si may be attributed to the weak electronegativity of the anions (i.e., B, Ge, and Si) and a consequent low resistance to oxidation. Additionally, B forms three thermodynamically stable compounds at various oxidation states with Mg, namely MgB2 (oxidation state of B is −1), MgB4 (B−0.5), and MgB7 (B−0.28). While MgB2 is stable against Mg metal (highest reducing conditions), at increasing voltages (~0.05 V vs. Mg/Mg2+), compounds with higher B oxidation states become stable, limiting the oxidative stability of MgB2. On the other hand, Cl and Mg only form MgCl2 as a stable binary, which oxidizes directly to Cl2 gas at ~3.39 V vs. Mg/Mg2+. Notably, MgCl2 is used as a precursor for Mg-Al-Cl-based electrolytes and its limited solubility in an ether-based solvent (typically used in Mg batteries) is well-documented (Doe et al., 2014; Canepa et al., 2015a). Therefore, MgCl2 may already be present in existing electrolytes, given its stability against Mg-metal (Figure 2), and may inherently protect the anode against further reactions with the electrolyte. In light of this, the role of MgCl2 as a potential protective coating layer on the Mg metal electrode needs to be further investigated.
Within each group of compounds of Figure 2 (i.e., each column of Figure 1), there is a strong correlation between the electronegativity of the anion and the oxidative stability. For example, within halogen compounds (orange bars), the oxidative stability rigorously follows the order MgF2 > MgCl2 > MgBr2 > MgI2, which correlates with the relative order of electronegativity of F > Cl > Br > I. Analogous trends can be observed among chalcogens (gray bars) and pnictogens (yellow). From this analysis we concluded that the electronegativity of the anion can be used as a proxy for the oxidative potential of Mg binary compounds since it describes the ability of the anion to limit an oxidation reaction.
Electrochemical Stability Windows of Mg-Ternaries
Figure 3 shows the voltage windows of Mg ternary and quaternary oxides, while Figure 4 shows the voltage windows of Mg ternary non-oxides (i.e., sulfides, selenides, tellurides, and a hydride).
Figure 3. Electrochemical stability windows of Mg-ternary and quaternary oxides, indicating the voltages (vs. Mg metal) at which the compound is stable against decomposition. Compounds that are not stable at any voltage are not shown. Ternaries are grouped by the periodic table column of the non-Mg, non-anion elements and ordered within each group by increasing electronegativity of the non-Mg cation. For systems with multiple compositions, compounds are ordered by increasing reductive stability. The text next to each bar indicates the width of the voltage window and the decomposition products at the reductive and oxidative limits. Compounds sharing common decomposition products, such as MgO or O2 are grouped together by brackets.
Figure 4. Electrochemical stability windows of Mg-ternary non-oxides, indicating the voltages (vs. Mg metal) at which the compound is stable against decomposition. Compounds that are not stable at any voltage are not shown. Ternaries are grouped by anion in order of increasing electronegativity and ordered within each group by increasing electronegativity of the non-Mg cation (e.g., P, S). The text next to each bar indicates the width of the voltage window and the decomposition products at the reductive and oxidative limits. Compounds sharing common decomposition products such as MgS or S are grouped together by brackets.
The widths of the voltage windows are written next to the respective bars on either the left or the right side. Decomposition products at the reductive (oxidative) stability limits are written to the left (right) of the bars. For compounds within a group that share a common decomposition product (such as, MgO, O2 in Figure 3, and MgS, S or MgSe, Se in Figure 4), the common compounds are factored out and indicated in brackets. The compounds shown are based on the elements highlighted in Figure 1 and a set of stable Mg-niobates, Mg-titanates, Mg-titanium-phosphates and Mg-zirconium-phosphates are plotted vs. Mg/Mg2+ as the reference. Compounds that are not thermodynamically stable (i.e., with a non-zero decomposition energy or energy above the convex hull) are not plotted. For example, Mg14Si5O24 is calculated to decompose into Mg2SiO4 and MgO and hence not included in Figure 3. Of note, Mg does not form ternary halides where the halogen is the anion, according to the structures available in the ICSD. Thus, no compounds in the ternary phase spaces of Mg-P-Cl, Mg-N-F, Mg-B-F are known to exist. Instead, we found that the stable Mg-ternaries are ternary chalcogenides, where the anion is oxygen, sulfur, selenium, or tellurium (except for the Mg-borohydride).
Based on Figures 3, 4, we observed that Mg ternaries do not show reductive stability against Mg metal, as indicated by the lack of reductive stability down to 0 V for any compound considered. The ternary with the best reductive stability is Mg(BH4)2 (purple bar in Figure 4), which is stable up to 0.01 V against Mg metal. Further, none of the ternary compounds exceed the anodic stability of MgF2 (~5.8 V, Figure 2). Among the ternaries, the Mg-B-O-based compounds, Mg2SiO4 and Mg3(PO4)2 have the widest stability windows, with voltage window widths >2.5 V. Additionally, there exist a few ternary oxides, such as MgP4O11 (~4.55 V), MgS2O7 (~4.45 V), and a quaternary Mg0.5Ti2(PO4)3 (~3.82 V) which have significantly high oxidative stability and may represent potential protective coatings for high-voltage oxide cathodes (Chen et al., 2017).
In general, trends in oxidative and reductive stability from Figures 3, 4 can be explained by analyzing the species most prone to oxidation and reduction, respectively. In most of the ternary compounds considered, the species most prone to oxidation is the anion since the other elements are already at their highest oxidation states (e.g., P5+ in Mg3(PO4)2). Thus, the susceptibility of the anion to be oxidized dictates the oxidative stability of the compound. For example, among the phosphates (yellow bars in Figure 3), thio-phosphates (orange bar in Figure 4), and seleno-phosphates (gray bar in Figure 4), phosphates exhibit the highest oxidative stabilities compared to MgPS3 and MgPSe3 because O2− is more difficult to oxidize than S2− or Se2−. Given that the electronegativity of the anion directly corresponds to the tendency of the anion to attract electrons and its resistance to oxidation, there is a high degree of correlation between increasing anion electronegativity (e.g., O > S > Se) (Pauling, 1932) and higher oxidative stabilities of binary (Figure 2) and ternary (Figures 3, 4) oxides compared to sulfides and selenides, respectively. Additionally, the hybridization of the anion (e.g., between O2− and P5+ in moieties) tends to stabilize it by lowering the energy of its electronic states, making the anion more difficult to oxidize. For example, binary MgO, where O2− hardly hybridizes with Mg2+, oxidizes at ~3.10 V vs. Mg. On the other hand, most Mg-ternary oxides (including the phosphates) oxidize at higher voltages (i.e., exhibit superior oxidative stability) due to the hybridization of the O2− by the non-Mg cation, such as P5+, S6+, etc.
The reductive stability of ternary compounds depends primarily on two key metrics: (i) the electronegativity of the species that undergoes reduction, which is the non-Mg cation in ternary compounds, and (ii) the electronegativity of the anion that does not undergo reduction but regulates the thermodynamic stability of the ternary compound vs. the corresponding binary compounds. Notably, reductive stability correlates inversely with the electronegativity of the non-Mg cation species, since larger electronegativities reflect higher attraction toward electrons and a higher propensity for reduction. For example, the reductive stability of ternary compounds (Figure 3) follows the order Mg-Cl-O (~2.78 V vs. Mg) < Mg-S-O (~1.65 V) < Mg-P-O (~0.9 V) < Mg-Si-O (~0.47 V) < Mg-Al-O (~0.14 V), which is the inverse of the electronegativity trends, namely Cl (3.16) > S (2.58) > P (2.19) > Si (1.90) > Al (1.61) (Pauling, 1932). In the case of quaternary systems, such as Mg-Ti-P-O and Mg-Zr-P-O, we predicted that P5+ reduces in preference to Ti4+ and Zr4+ (brown bars in Figure 3), which is consistent with the larger electronegativity of P (2.19) vs. Ti (1.54) and Zr (1.33).
Importantly, higher electronegativity of the anion results in poorer reductive stability of the ternary compound. For example, the reductive stability among Mg-Ge-, Mg-Sc-, Mg-In-ternary oxides follows Mg-Ge-O (~1.32 V) < Mg-Ge-S (~1.13 V) < Mg-Ge-Se (1.08 V), Mg-Sc-S (~0.36 V) < Mg-Sc-Se (0.18 V), and Mg-In-S (~1.04 V) < Mg-In-Se (~0.93 V) < Mg-In-Te (~0.67 V), respectively, consistent with the anion electronegativity trend (O > S > Se > Te). Note that higher anion electronegativity leads to more stable Mg-binary compounds, i.e., Mg-binaries with larger stability windows (Figure 2), which are common decomposition products under reducing conditions. A more stable Mg-binary reflects a larger thermodynamic driving force for reduction, as quantified by the corresponding formation energy (MgO ~ −3.06 eV/atom, MgS ~ −1.76 eV/atom, MgSe ~ −1.25 eV/atom, and MgTe ~ −0.87 eV/atom) (Jain et al., 2013), resulting in a lower reductive stability. Interestingly, the compound with the highest reductive stability, Mg(BH4)2, is composed of a low electronegative anion and a non-Mg cation, H (2.20) and B (2.04), respectively. Thus, minimizing the electronegativities of both the non-Mg-cations and the anions could be the key to discovering ternary compounds that are stable against Mg-metal.
Notable exceptions to the aforementioned trends in reductive stability vs. (non-Mg cation/anion) electronegativity can be observed across different chemistries in Figures 3, 4. For example, electronegativity of B (2.04) > Ga (1.81) > Al (1.61), but the reductive stability of Mg-Al-O (~0.14 V) > Mg-B-O (~0.58 V) > Mg-Ga-O (~1.13 V). Similar trends can be observed among Mg-IVA-O, and Mg-VA-O compounds (Figure 3). Such anomalies can be attributed to two factors that override non-Mg-cation electronegativity trends: (i) stability of Mg-(IIIA/IVA/VA) binaries (signifying the thermodynamic driving force to form decomposition products), and (ii) the relative position of the empty electronic states of IIIA/IVA/VA elements, as influenced by the extend of hybridization with oxygen (difficulty in reducing the ternary compound). For example, the highest oxidative stability of binary Mg-Al alloys [~0.06 V (Jain et al., 2013), not shown in Figure 2)] is lower than both Mg-B compounds (~0.53 V, Figure 2) and Mg-Ga alloys (~0.19 V, not shown). On the other hand, the significant hybridization of the electronic states of P with O likely pushes the empty (anti-bonding) P states to higher energy levels, making P difficult to reduce in ternary Mg-P-O, compared to As in Mg-As-O and N in Mg-N-O.
In the case of reductive stability vs. anion electronegativities, the stability of Mg-Al-O (~0.14 V) > Mg-Al-S (~0.52 V), and Mg-P-O (~0.9 V) > Mg-P-S (~1.20 V), despite the electronegativity of O > S is another notable exception. Here, the discrepancy can be attributed to the stability of Al-O and P-O bonds in comparison to Al-S and P-S bonds, as quantified by the formation energies (Al2O3 ~ −3.44 eV/atom, Al2S3 ~ −1.46 eV/atom and P2O5 ~ −2.46 eV/atom and P2S5 ~ −0.64 eV/atom) (Jain et al., 2013). The higher stability of Al-O and P-O bonds is possibly due to better hybridization of Al and P among the oxides vs. sulfides, respectively. Thus, despite MgO creating a larger thermodynamic driving force for reduction than MgS (as indicated by the stability windows in Figure 2), the lack of affinity for S from Al and P in Mg-Al-, and Mg-P-ternaries facilitates the reduction of Al3+ and P4+/5+, respectively, in the ternary sulfides compared to the oxides.
Potential Candidate Materials
Based on the voltage windows of the Mg binaries, ternaries and quaternaries in Figures 2–4, we suggest potential coatings on both the Mg metal//Mg electrolyte and the Mg electrolyte//cathode interfaces. At the cathode interface, the oxidative stability should be high for candidate compounds. Among the binaries, only MgF2 has an oxidation limit above 4.0 V, whereas among the ternaries, including Mg(PO3)2, MgP4O11, Mg(NO3)2, and MgS2O7 show oxidation limits above 4.0 V. Note that among the candidate materials, those with the widest voltage windows should be given preference, which may enable compatibility with liquid electrolytes that are stable against Mg metal. Therefore, among the high-oxidation-limit compounds, MgF2, Mg(PO3)2, MgP4O11, and Mg0.5Ti2(PO4)3, which have the widest voltage windows (all > 2.0 V), should be considered the most promising candidate materials.
For the Mg metal//Mg electrolyte interface, the reductive stability of a candidate compound should ideally be ~0 V vs. Mg metal. In this context, Mg(BH4)2, with a reductive stability of ~0.01 V vs. Mg is a promising candidate for a protective anode coating. Previous experiments utilizing Mg(BH4)2-containing electrolytes have reported the formation of a Mg-conducting interphase layer against Mg-metal with an oxidative stability of 1.7 V vs. Mg, which is generally in accordance with our computational results (1.25 V vs. Mg) (Mohtadi et al., 2012; Arthur et al., 2017). The higher oxidative stability of Mg(BH4)2 observed in experiments (1.7 V vs. Mg) compared to in theory (1.25 V vs. Mg) could be due to kinetic stability, which is not accounted for in our calculations. Thus, Mg(BH4)2 should be further investigated as a protective coating on the Mg-metal anode. Additionally, in scenarios where the reductive stability is < ~0.5 V, such as MgAl2O4, Mg2SiO4 (Figure 3), Mg(ScS2)2, and Mg(ScSe2)2 (Figure 4), the compounds may exist in a metastable manner and may still be valid candidates. For example, in Li-ion batteries the solid electrolyte, garnet-Li7La3Zr2O12, has an estimated reductive stability of ~0.1 V vs. Li but has been shown to be metastable against Li metal (Richards et al., 2015; Ma et al., 2016). However, recent theoretical and experimental studies have shown that Mg(ScS2)2 and Mg(ScSe2)2 tend to decompose to binary MgS/MgSe and ScS/ScSe against Mg metal, ruling out any metastable existence (Canepa et al., 2017b,c). Another case to consider is when the Mg metal anode is replaced by Bi (or Sb or their alloys) as the reductive potential of the anode is shifted by up to ~+0.32 V vs. Mg metal (Arthur et al., 2012). In case these alternative anodes are used, several coating materials, such as MgAl2O4 or Mg(ScSe2)2, could be envisioned as potential coating materials. Nevertheless, changing the anode chemistry can not only change the overall energy density of the cell but also introduce additional over-potentials for Mg alloying at the anode. Notably, all binaries considered should be stable vs. Mg metal, except for MgP4, MgB4, and MgB7 (Figure 2), and are candidates for protective coatings at the anode//electrolyte interface. Specifically, Mg-halides, including MgF2, MgCl2, MgBr2, which have voltage windows wider than 2.0 V, should be considered as the most promising candidates.
A number of studies have suggested that the Cl− in magnesium-aluminum-chloride-based electrolytes can protect the Mg-metal anode during Mg deposition via adsorption on the Mg-metal surface (Aurbach et al., 2002; Doe et al., 2014; Canepa et al., 2015a,b; See et al., 2015, 2017; Salama et al., 2017). Our results suggest that MgCl2 is stable against the highly reductive environment of Mg-metal, showing a wide stability window ~3.39 V. We speculate that a layer of MgCl2 may form in situ as a protective coating, which is further justified by the sparing solubility of this salt in ether-based solvents (Doe et al., 2014; Canepa et al., 2015a; Salama et al., 2017). Therefore, a careful experimental characterization of the Mg//electrolyte interface will shed light on the role of the speciation of Cl in the form of MgCl2 or as a free ion.
For all of the suggested anode or cathode coating materials, a thorough evaluation of Mg2+ mobility is required to verify their viability as actual coating materials. Mobility evaluations are especially necessary to demonstrate proof-of-concept oxidative coatings that can enable high voltage cathodes [such as MgxCr2O4 (Chen et al., 2017), MgxMn2O4, (Sai Gautam et al., 2017), and MgxV2O5 (Sai Gautam et al., 2015)] in conjunction with current liquid electrolytes and Mg-metal. Note that the Mg2+ migration barrier has been calculated for a number of Mg-binaries in a prior study (Canepa et al., 2017b), including MgO (~1,800 meV), MgS (~900 meV), and MgSe (~700 meV) of Figure 5, and a few ternaries, such as Mg(ScSe2)2 (~375 meV), Mg(InS2)2 (~488 meV), and Mg(ScS2)2 (~415 meV), while more studies are in progress for other candidates listed in this work. The poor bulk Mg mobility causes MgO and MgS to be inactive passivating materials that limit any Mg transference, despite their wide stability ranges (0–3.1 V for MgO and 0–1.6 V for MgS). Similarly, poor Mg mobility in bulk Mg0.5Ti2(PO4)3 [>1 eV (Canepa et al., 2017a)] will hinder its use as a protective oxidative coating. Nevertheless, our study identifies a tractable list of possible coating and electrolyte candidates in which Mg2+ mobility must be estimated, based on their calculated electrochemical stabilities.
Figure 5. Plot of migration barriers of MgO (red), MgS (blue), and MgSe (green) as calculated in Canepa et al. (2017b). The high migration barriers of MgO, which is predicted to be stable vs. Mg metal and to have a reasonable oxidation limit (3.08 V vs. Mg metal), demonstrate the necessity of Mg2+ diffusivity data in determining the viability of potential coating and electrolyte materials.
Conclusion
In this study, we evaluated, using density functional theory calculations, the electrochemical stability windows for non-redox-active Mg binary, ternary, and selected quaternary compounds in order to identify potential coating materials for Mg batteries. From the Mg binaries considered, we identified Mg-halides, specifically MgCl2 and MgBr2, as potential anode coating materials based on their reductive stability (at 0 V vs. Mg/Mg2+). We also suggested Mg(BH4)2, MgAl2O4, and Mg2SiO4, as possible ternary anode coating materials, given their reductive stability below 0.5 V, with MgAl2O4 and Mg2SiO4 exhibiting a voltage window that is >2.0 V wide. Additionally, we expect MgF2, Mg(PO3)2, and MgP4O11 to be promising candidates for protecting high-voltage cathodes against typical Mg electrolytes. While careful evaluation of Mg mobility in candidate materials is essential, this work identifies specific chemistries as well as general guidelines on compound stabilities that will be useful to design practical coating materials in Mg batteries.
Author Contributions
PC and GS conceived the project. TC, GS, and PC performed the simulation and analyzed the data. TC, GS, and PC wrote the manuscript and discussed it with GC.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
PC acknowledges support from the Singapore Ministry of Education Academic Fund Tier 1 (R-284-000-186-133). This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE 1106400. The work is also supported by the Joint Center for Energy Storage Research (JCESR), an Energy Innovation Hub funded by the U.S. Department of Energy, Office of Science and Basic Energy Sciences, through Subcontract 3F-31144. Any opinion, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00024/full#supplementary-material
References
Arthur, T. S., Glans, P.-A., Singh, N., Tutusaus, O., Nie, K., Liu, Y.-S., et al. (2017). Interfacial insight from operando XAS/TEM for magnesium metal deposition with borohydride electrolytes. Chem. Mater. 29, 7183–7188. doi: 10.1021/acs.chemmater.7b01189
Arthur, T. S., Singh, N., and Matsui, M. (2012). Electrodeposited Bi, Sb and Bi1-xSbx alloys as anodes for Mg-ion batteries. Electrochem. Commun. 16, 103–106. doi: 10.1016/j.elecom.2011.12.010
Aurbach, D., Gizbar, H., Schechter, A., Chusid, O., Gottlieb, H. E., Gofer, Y., et al. (2002). Electrolyte solutions for rechargeable magnesium batteries based on organomagnesium chloroaluminate complexes. J. Electrochem. Soc. 149, A115–A121. doi: 10.1149/1.1429925
Aurbach, D., Lu, Z., Schechter, A., Gofer, Y., Gizbar, H., Turgeman, R., et al. (2000). Prototype systems for rechargeable magnesium batteries. Nature 407, 724–727. doi: 10.1038/35037553
Aurbach, D., Talyosef, Y., Markovsky, B., Markevich, E., Zinigrad, E., Asraf, L., et al. (2004). Design of electrolyte solutions for Li and Li-ion batteries: a review. Electrochim. Acta 50, 247–254. doi: 10.1016/j.electacta.2004.01.090
Bachman, J. C., Muy, S., Grimaud, A., Chang, H.-H., Pour, N., Lux, S. F., et al. (2015). Inorganic solid-state electrolytes for lithium batteries: mechanisms and properties governing ion conduction. Chem. Rev. 116, 140–162. doi: 10.1021/acs.chemrev.5b00563
Bergerhoff, G., and Brown, I. (1987). “Crystallographic databases,” in Inorganic Crystal Structure Database (Bonn; Cambridge; Chester: International Union of Crystallography), 77–95.
Canepa, P., Bo, S.-H., Gautam, G. S., Key, B., Richards, W. D., Shi, T., et al. (2017b). High magnesium mobility in ternary spinel chalcogenides. Nat. Commun. 8:1759. doi: 10.1038/s41467-017-01772-1
Canepa, P., Gautam, G. S., Malik, R., Jayaraman, S., Rong, Z., Zavadil, K. R., et al. (2015b). Understanding the initial stages of reversible Mg deposition and stripping in inorganic nonaqueous electrolytes. Chem. Mater. 27, 3317–3325. doi: 10.1021/acs.chemmater.5b00389
Canepa, P., Jayaraman, S., Cheng, L., Rajput, N. N., Richards, W. D., Gautam, G. S., et al. (2015a). Elucidating the structure of the magnesium aluminum chloride complex electrolyte for magnesium-ion batteries. Energy Environ. Sci. 8, 3718–3730. doi: 10.1039/C5EE02340H
Canepa, P., Sai Gautam, G., Broberg, D., Bo, S.-H., and Ceder, G. (2017c). Role of point defects in spinel Mg chalcogenide conductors. Chem. Mater. 29, 9657–9667. doi: 10.1021/acs.chemmater.7b02909
Canepa, P., Sai Gautam, G., Hannah, D. C., Malik, R., Liu, M., Gallagher, K. G., et al. (2017a). Odyssey of multivalent cathode materials: open questions and future challenges. Chem. Rev. 117, 4287–4341. doi: 10.1021/acs.chemrev.6b00614
Carter, T. J., Mohtadi, R., Arthur, T. S., Mizuno, F., Zhang, R., Shirai, S., et al. (2014). Boron clusters as highly stable magnesium-battery electrolytes. Angewandte Chem. Int. Ed. 53, 3173–3177. doi: 10.1002/anie.201310317
Chen, T., Sai Gautam, G., Huang, W., and Ceder, G. (2017). First-principles study of the voltage profile and mobility of Mg intercalation in a chromium oxide spinel. Chem. Mater. 30, 153–162. doi: 10.1021/acs.chemmater.7b04038
Cohen, Y. S., Cohen, Y., and Aurbach, D. (2000). Micromorphological studies of lithium electrodes in alkyl carbonate solutions using in situ atomic force microscopy. J. Phys. Chem. B 104, 12282–12291. doi: 10.1021/jp002526b
Doe, R. E., Han, R., Hwang, J., Gmitter, A. J., Shterenberg, I., Yoo, H. D., et al. (2014). Novel, electrolyte solutions comprising fully inorganic salts with high anodic stability for rechargeable magnesium batteries. Chem. Commun. 50, 243–245. doi: 10.1039/C3CC47896C
Gofer, Y., Turgeman, R., Cohen, H., and Aurbach, D. (2003). XPS investigation of surface chemistry of magnesium electrodes in contact with organic solutions of organochloroaluminate complex salts. Langmuir 19, 2344–2348. doi: 10.1021/la026642c
Goodenough, J. B., and Kim, Y. (2009). Challenges for rechargeable Li batteries. Chem. Mater. 22, 587–603. doi: 10.1021/cm901452z
Guerfi, A., Dontigny, M., Charest, P., Petitclerc, M., Lagacé, M., Vijh, A., et al. (2010). Improved electrolytes for Li-ion batteries: Mixtures of ionic liquid and organic electrolyte with enhanced safety and electrochemical performance. J. Power Sources 195, 845–852. doi: 10.1016/j.jpowsour.2009.08.056
Hahn, N. T., Seguin, T. J., Lau, K.-C., Liao, C., Ingram, B. J., Persson, K. A., et al. (2018). Enhanced stability of the carba-closo-dodecaborate anion for high-voltage battery electrolytes through rational design. J. Am. Chem. Soc. 140, 11076–11084. doi: 10.1021/jacs.8b05967
Hannah, D. C., Sai Gautam, G., Canepa, P., and Ceder, G. (2018). On the balance of intercalation and conversion reactions in battery cathodes. Adv. Energy Mater. 8:1800379. doi: 10.1002/aenm.201800379
Hart, G. L., and Forcade, R. W. (2008). Algorithm for generating derivative structures. Phys. Rev. B 77:224115. doi: 10.1103/PhysRevB.77.224115
Hart, G. L., and Forcade, R. W. (2009). Generating derivative structures from multilattices: algorithm and application to hcp alloys. Phys. Rev. B 80:014120. doi: 10.1103/PhysRevB.80.014120
Hart, G. L., Nelson, L. J., and Forcade, R. W. (2012). Generating derivative structures at a fixed concentration. Comput. Mater. Sci. 59, 101–107. doi: 10.1016/j.commatsci.2012.02.015
Jain, A., Hautier, G., Moore, C. J., Ong, S. P., Fischer, C. C., Mueller, T. A., et al. (2011). High-throughput infrastructure for density functional theory calculations. Comput. Mater. Sci. 50, 2295–2310. doi: 10.1016/j.commatsci.2011.02.023
Jain, A., Ong, S. P., Hautier, G., Chen, W., Richards, W. D., Dacek, S., et al. (2013). Commentary: the materials project: a materials genome approach to accelerating materials innovation. Apl. Mater. 1:011002. doi: 10.1063/1.4812323
Kamaya, N., Homma, K., Yamakawa, Y., Hirayama, M., Kanno, R., Yonemura, M., et al. (2011). A lithium superionic conductor. Nat. Mater. 10, 682–686. doi: 10.1038/nmat3066
Kato, Y., Hori, S., Saito, T., Suzuki, K., Hirayama, M., Mitsui, A., et al. (2016). High-power all-solid-state batteries using sulfide superionic conductors. Nat. Energy 1:16030. doi: 10.1038/nenergy.2016.30
Keyzer, E. N., Glass, H. F., Liu, Z., Bayley, P. M., Dutton, S. N. E., et al. (2016). Mg (PF6) 2-based electrolyte systems: understanding electrolyte–electrode interactions for the development of mg-ion batteries. J. Am. Chem. Soc. 138, 8682–8685. doi: 10.1021/jacs.6b04319
Kohn, W., and Sham, L. J. (1965). Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, 4A113. doi: 10.1103/PhysRev.140.A1133
Levi, E., Gofer, Y., and Aurbach, D. (2009). On the way to rechargeable Mg batteries: the challenge of new cathode materials. Chem. Mater. 22, 860–868. doi: 10.1021/cm9016497
Ling, C., and Zhang, R. (2017). Manganese dioxide as rechargeable magnesium battery cathode. Front. Energy Res. 5:30. doi: 10.3389/fenrg.2017.00030
Ling, C., Zhang, R., Arthur, T. S., and Mizuno, F. (2015). How general is the conversion reaction in Mg battery cathode: a case study of the magnesiation of α-MnO2. Chem. Mater. 27, 5799–5807. doi: 10.1021/acs.chemmater.5b02488
Lipson, A. L., Han, S.-D., Pan, B., See, K. A., Gewirth, A. A., Liao, C., et al. (2016). Practical stability limits of magnesium electrolytes. J. Electrochem. Soc. 163, A2253–A2257. doi: 10.1149/2.0451610jes
Liu, M., Rong, Z., Malik, R., Canepa, P., Jain, A., Ceder, G., et al. (2015). Spinel compounds as multivalent battery cathodes: a systematic evaluation based on ab initio calculations. Energy Environ. Sci. 8, 964–974. doi: 10.1039/C4EE03389B
Lu, Z., Schechter, A., Moshkovich, M., and Aurbach, D. (1999). On the electrochemical behavior of magnesium electrodes in polar aprotic electrolyte solutions. J. Electroanalyt. Chem. 466, 203–217. doi: 10.1016/S0022-0728(99)00146-1
Ma, C., Cheng, Y., Yin, K., Luo, J., Sharafi, A., Sakamoto, J., et al. (2016). Interfacial stability of Li metal–solid electrolyte elucidated via in situ electron microscopy. Nano Lett. 16, 7030–7036. doi: 10.1021/acs.nanolett.6b03223
Marom, R., Amalraj, S. F., Leifer, N., Jacob, D., and Aurbach, D. (2011). A review of advanced and practical lithium battery materials. J. Mater. Chem. 21, 9938–9954. doi: 10.1039/c0jm04225k
Masquelier, C. (2011). Solid electrolytes: lithium ions on the fast track. Nat. Mater. 10, 649–650. doi: 10.1038/nmat3105
Mohtadi, R., Matsui, M., Arthur, T. S., and Hwang, S. J. (2012). Magnesium borohydride: from hydrogen storage to magnesium battery. Angewandte Chem. Int. Ed. 51, 9780–9783. doi: 10.1002/anie.201204913
Muldoon, J., Bucur, C. B., and Gregory, T. (2014). Quest for nonaqueous multivalent secondary batteries: magnesium and beyond. Chem. Rev. 114, 11683–11720. doi: 10.1021/cr500049y
Muldoon, J., Bucur, C. B., Oliver, A. G., Sugimoto, T., Matsui, M., Kim, H. S., et al. (2012). Electrolyte roadblocks to a magnesium rechargeable battery. Energy Environ. Sci. 5, 5941–5950. doi: 10.1039/c2ee03029b
Ong, S. P., Richards, W. D., Jain, A., Hautier, G., Kocher, M., Cholia, S., et al. (2013). Python materials genomics (pymatgen): a robust, open-source python library for materials analysis. Comput. Mater. Sci. 68, 314–319. doi: 10.1016/j.commatsci.2012.10.028
Pauling, L. (1932). The nature of the chemical bond. Iv. The energy of single bonds and the relative electronegativity of atoms. J. Am. Chem. Soc. 54, 3570–3582.
Pour, N., Gofer, Y., Major, D. T., and Aurbach, D. (2011). Structural analysis of electrolyte solutions for rechargeable Mg batteries by stereoscopic means and DFT calculations. J. Am. Chem. Soc. 133, 6270–6278. doi: 10.1021/ja1098512
Richards, W. D., Miara, L. J., Wang, Y., Kim, J. C., and Ceder, G. (2015). Interface stability in solid-state batteries. Chem. Mater. 28, 266–273. doi: 10.1021/acs.chemmater.5b04082
Rosenberg, M., and Nicolau, P. (1964). Electrical properties and cation migration in MgMn2O4. Phys. Status Solidi B 101, 101–110. doi: 10.1002/pssb.19640060107
Sai Gautam, G., Canepa, P., Abdellahi, A., Urban, A., Malik, R., and Ceder, G. (2015). The intercalation phase diagram of Mg in V2O5 from first-principles. Chem. Mater. 27, 3733–3742. doi: 10.1021/acs.chemmater.5b00957
Sai Gautam, G., Canepa, P., Urban, A., Bo, S.-H., and Ceder, G. (2017). Influence of inversion on Mg mobility and electrochemistry in spinels. Chem. Mater. 29, 7918–7930. doi: 10.1021/acs.chemmater.7b02820
Salama, M., Shterenberg, I., Shimon, L. J. W., Keinan-Adamsky, K., Afri, M., et al. (2017). Structural analysis of magnesium chloride complexes in dimethoxyethane solutions in the context of Mg batteries research. J. Phys. Chem. C 121, 24909–24918. doi: 10.1021/acs.jpcc.7b05452
See, K. A., Chapman, K. W., Zhu, L., Wiaderek, K. M., Borkiewicz, O. J., Barile, C. J., et al. (2015). The interplay of Al and Mg speciation in advanced Mg battery electrolyte solutions. J. Am. Chem. Soc. 138, 328–337. doi: 10.1021/jacs.5b10987
See, K. A., Liu, Y.-M., Ha, Y., Barile, C. J., and Gewirth, A. A. (2017). Effect of concentration on the electrochemistry and speciation of the magnesium aluminum chloride complex electrolyte solution. ACS Appl. Mater. Interfaces 9, 35729–35739. doi: 10.1021/acsami.7b08088
Shao, Y., Liu, T., Li, G., Gu, M., Nie, Z., Engelhard, M., et al. (2013). Coordination chemistry in magnesium battery electrolytes: how ligands affect their performance. Sci. Rep. 3:3130. doi: 10.1038/srep03130
Sun, X., Bonnick, P., Duffort, V., Liu, M., Rong, Z., Persson, K. A., et al. (2016). A high capacity thiospinel cathode for Mg batteries. Energy Environ. Sci. 9, 2273–2277. doi: 10.1039/C6EE00724D
Tutusaus, O., Mohtadi, R., Arthur, T. S., Mizuno, F., Nelson, E. G., and Sevryugina, Y. V. (2015). An efficient halogen-free electrolyte for use in rechargeable magnesium batteries. Angewandte Chem. 127, 8011–8015. doi: 10.1002/ange.201412202
Tutusaus, O., Mohtadi, R., Singh, N., Arthur, T. S., and Mizuno, F. (2016). Study of electrochemical phenomena observed at the Mg metal/electrolyte interface. ACS Energy Lett. 2, 224–229. doi: 10.1021/acsenergylett.6b00549
Verma, P., Maire, P., and Novák, P. (2010). A review of the features and analyses of the solid electrolyte interphase in Li-ion batteries. Electrochim. Acta 55, 6332–6341. doi: 10.1016/j.electacta.2010.05.072
Keywords: Mg batteries, first-principles calculation, density functional theory, coating materials, intercalation batteries, solid electrolytes, multivalent ion batteries
Citation: Chen T, Ceder G, Sai Gautam G and Canepa P (2019) Evaluation of Mg Compounds as Coating Materials in Mg Batteries. Front. Chem. 7:24. doi: 10.3389/fchem.2019.00024
Received: 12 November 2018; Accepted: 10 January 2019;
Published: 30 January 2019.
Edited by:
Federico Cesano, University of Turin, ItalyReviewed by:
Timothy Sean Arthur, Toyota, United StatesFride Vullum-Bruer, Norwegian University of Science and Technology, Norway
Michael Salama, Bar-Ilan University, Israel
Copyright © 2019 Chen, Ceder, Sai Gautam and Canepa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gopalakrishnan Sai Gautam, Z2F1dGFtOTFAcHJpbmNldG9uLmVkdQ==
Pieremanuele Canepa, cGNhbmVwYUBudXMuZWR1LnNn