Next Article in Journal
Innovative Purification Method of Ovatodiolide from Anisomeles indica to Induce Apoptosis in Human Gastric Cancer Cells
Next Article in Special Issue
Linker Functionalization Strategy for Water Adsorption in Metal–Organic Frameworks
Previous Article in Journal
Microplastics (Polystyrene) Exposure Induces Metabolic Changes in the Liver of Rare Minnow (Gobiocypris rarus)
Previous Article in Special Issue
Cinnamal Sensing and Luminescence Color Tuning in a Series of Rare-Earth Metal−Organic Frameworks with Trans-1,4-cyclohexanedicarboxylate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Redox Hyperactive MOF for Li+, Na+ and Mg2+ Storage

1
Faculty of Chemistry and Pharmacy, University of Sofia, 1164 Sofia, Bulgaria
2
Institute of General and Inorganic Chemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
3
Faculty of Chemistry, Wroclaw University of Science and Technology, 50-373 Wroclaw, Poland
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(3), 586; https://doi.org/10.3390/molecules27030586
Submission received: 26 December 2021 / Revised: 8 January 2022 / Accepted: 12 January 2022 / Published: 18 January 2022

Abstract

:
To create both greener and high-power metal-ion batteries, it is of prime importance to invent an unprecedented electrode material that will be able to store a colossal amount of charge carriers by a redox mechanism. Employing periodic DFT calculations, we modeled a new metal-organic framework, which displays energy density exceeding that of conventional inorganic and organic electrodes, such as Li- and Na-rich oxides and anthraquinones. The designed MOF has a rhombohedral unit cell in which an Ni(II) node is coordinated by 2,5-dicyano-p-benzoquinone linkers in such a way that all components participate in the redox reaction upon lithiation, sodiation and magnesiation. The spatial and electronic changes occurring in the MOF after the interaction with Li, Na and Mg are discussed on the basis of calculated electrode potentials versus Li0/Li+, Na0/Na+ and Mg0/Mg2+, respectively. In addition, the specific capacities and energy densities are calculated and used as a measure for the electrode applicability of the designed material. Although the highest capacity and energy density are predicted for Li storage, the greater structural robustness toward Na and Mg uptake suggests a higher cycling stability in addition to lower cost. The theoretical results indicate that the MOF is a promising choice for a green electrode material (with <10% heavy metal content) and is well worth experimental testing.

Graphical Abstract

1. Introduction

Matching the high-power and large-scale energy storage with environmental constraints entails to consider replacement of the current lithium-ion batteries (LIBs) with technologies, in which more abundant, cheaper and safer elements are used [1,2,3,4]. In this respect, sodium- and magnesium-ion batteries (SIBs and MIBs) are the most prospective alternatives [5,6]. Both SIBs and MIBs operate by the same mechanism as LIBs, comprising the reversible exchange of Na+ and Mg2+ ions between two electrodes via electrolyte [7,8]. Given that sodium and magnesium have higher redox potentials versus SHE than lithium (−2.71 V for Na0/Na+ and −2.37 V for Mg0/Mg2+ compared to −3.04 V for Li0/Li+ [9]), the enhancement of the energy density of SIBs and MIBs can only be achieved through the advancement of a new class of electrode materials that have the unique property of colossal specific capacity (i.e., able to exchange a huge amount of Na+ or Mg2+ ions)—a phenomenon not feasible for conventional electrodes based on transition metal oxides and phosphates. This particularity follows from the ability of the electrode material to provide simultaneously flexible structural sites for storage and diffusion of large amounts of Na+ or Mg2+ ions, as well as to compensate through high redox activity the charges of all transportable metal ions.
Metal-organic frameworks (MOFs) are a specific class of materials that could meet the above requirements due to their structural and redox versatility [10,11,12,13,14]. In contrast to purely organic battery materials [15,16,17], through the bonding of metal-containing nodes with polyfunctional organic linkers, MOFs are prone to form highly ordered structures displaying an enormous variety of arrangements, providing a diversity of diffusion paths for both alkali and alkaline earth ions [18,19,20,21]. The redox activity of MOFs can be differentiated depending on whether the metal centers or the organic ligands or both the metal and the organic entities are active [22,23,24,25]. For MOFs with redox-active metal centers (such as Fe2+/Fe3+, Cu+/Cu2+, V4+/V5+, Co2+/Co3+, etc.), the number of inserted ions (e.g., alkali or alkaline earth ions) is limited up to 1 mole per metal center [26,27]. This type of redox reaction mimics the operation of the conventional inorganic electrode materials. When the organic ligands take part in the redox reaction, more transportable ions can be inserted, and higher potentials can be achieved [23]. Further advantages of the redox activity of the organic ligands are associated with restricted bond rearrangement during metal ion insertion resulting in small structural changes in comparison with inorganic oxide and phosphate electrodes. Moreover, the participation of the organic ligands in the redox reaction is less dependent on the radius and charge of alkali or alkaline earth ions than in the case of inorganic electrodes.
The most suitable organic ligands included in redox reactions are those containing one or more electron-withdrawing groups like carbonyl and nitrile ones [28]. Of particular interest are the quinones—a class of carbonyl compounds derived from aromatic hydrocarbons. Both theoretical and experimental investigations have shown that quinone-based structures have excellent qualities for application in batteries with high stability/energy density/capacity [29,30,31,32,33,34,35,36,37]. Thus, a viable strategy to build a MOF for energy storage application is to use quinone derivatives functionalized with coordinating groups as organic linkers in combination with suitable metal ions. In order to increase the energy density of the material, the metal ions, if appropriately selected, can also participate in the redox reactions. This has been shown to work in a practical cell comprising Cu+/Cu2+ redox couple and two different organic linkers—anthraquinone dicarboxylate [23] and tricarboxyphenylamine [38], the latter achieving the very high potential of 4.3 V. Furthermore, Gu et al. [39] applied the previously synthesized conductive MOF—Cu3(2,3,6,7,10,11-hexahydroxytriphenylene)2 [40] as cathode material. The combination of Li-ion insertion capability with electronic conductivity resulted in an outstanding rate performance. A Cu(II)-based 2D MOF was shown to be able to accommodate lithium as well, while demonstrating high conductivity in both pristine and lithiated forms [41]. Other examples of a redox-active metal-ion couples are V4+/V5+ in a layered oxalatophosphate [42] and Fe2+/Fe3+ in a 1,4-benzenedicarboxylate-based MOF [26] for Li-ion batteries; Fe2+/Fe3+ [27,43] and Co2+/Co3+ [27] in Prussian-blue analogues for Na-ion batteries. Moreover, new MOF-based electrode materials for anodes in LIB (with low potential vs. Li0/Li+) have been synthesized and tested. Li et al. [44] reported a Zn/1,3,5-benzenetribenzoate-based MOF but electrochemical testing of the material resulted in a large irreversible capacity loss after the first cycle. The 2,6-naphthalene(COOLi)2 intercalation MOF has proven to have better cycling performance (96% capacity retention after 100 cycles) and low potential of 0.8 V vs. Li-metal [45].
Although nickel metal is widely used in battery technologies, its integration into redox-active MOFs is rare. Until now, there have been few reports, including the paper of Hameed et.al. [46], in which a Ni-based 1D-MOF/reduced graphene oxide composite was tested as electrode material in LIB. Nickel(II)-terephthalate coordination complexes have been shown to participate in a multi-electron redox process due to the change in the oxidation state of Ni and interaction of additional Li with the conjugated aromatic ligands [47].
Given the enormous variety of conceivable MOF structures, the exploration of their application as components of metal-ion batteries is still in its infancy. Here, we present a theoretical model of a MOF that is able to store both alkali and alkaline earth metals. Looking for small and light building elements, we have suggested an unreported (to our knowledge) MOF. The building elements are 2,5-dicyano-p-benzoquinone (2CN-BQ) used as a redox-active linker and Ni(II) serving as a metal node. The advantage of the proposed MOF over other MOF structures has three aspects. First, we have selected p-benzoquinone instead of o-benzoquinone in order to achieve consistent output voltages [31]. Then, the p-benzoquinone functionalization with CN-groups increases the overall electrochemical potential for metal ion storage together with the opportunity for monodentate coordination with Ni(II), allowing a 3D architecture of the MOF. Finally, the designed MOF is tested for its ability to store Li+, Na+ and Mg2+—a combination that has not been reported so far. The structural characterization of the MOF and its affinity towards Li, Na and Mg are quantified by DFT periodic calculations. The electrode potentials versus Li0/Li+, Na0/Na+ and Mg0/Mg2+, as well as the specific capacities and energy densities, are evaluated and used as a measure for the electrode applicability of the designed material. It has been shown that the suggested MOF can accommodate a huge number of Li and Na atoms (and a more modest number of Mg) with a remarkable charge transfer, exhibiting very promising capacities and energy densities with a great potential for electrode application.

2. Results and Discussion

2.1. Structure of Pristine MOF-S21

The unit cell of the pristine 3D-MOF model was constructed as a cubic cell with an octahedral Ni(II) at the vertices and 2,5-dicyano-p-benzoquinone (2CN-BQ) molecules as edges along the three lattice axes (Figure 1). The structure was neutralized by two Cl as counterions. Periodic boundary conditions were applied. After geometry optimization, a rhombohedral structure ( R 3 ¯ ) was obtained with unit cell parameters a = 11.67 Å and α = 56.16°, and volume 1024.8 Å3; such distortion has been witnessed experimentally in other coordination polymers of Ni(II) with linear bifunctional ligands (Au(CN)2) [48]. The nickel ion retained an octahedral coordination with the nitrogens from the cyano-groups (Figure 1). Because of the trigonal distortion, the initial pore volume underwent a strong shrinkage (i.e., of around ¼). This is an indication for unfeasibility of interpenetration pathways for the designed MOF (Figure 1). Although the pore sizes are not suitable for interpenetration, there are still a plethora of sites for metal ion insertion. Further on, the 3D-structure will be denoted as MOF-S21.
The constructed MOF was computationally tested as an electrode via successive introduction of three different metal atoms—lithium, sodium and magnesium. Important characteristics that can be estimated by means of molecular modelling are the capacity of the material toward the tested metals, the electrode potential and the charge distribution in the system, depending on the degree of metal uptake.

2.2. Metal Loading into MOF-S21

The computational procedure to test the prospect of MOF-S21 to serve as electrode material has the following main stages: (i) finding the most favorable structure for every degree of metal insertion; (ii) calculating the internal energy change between two consecutive degrees of metal insertion according to Equation (1); (iii) estimating the electrode potential against the respective metal according to Equation (2).

2.2.1. Lithiation of MOF-S21

The main redox-active parts of the pristine MOF are the cyano-quinones; they are expected to coordinate Li+ via their oxygens, so that in the LinMOF structures, n = 1 through 6, Li-atoms were placed next to the oxygens. For each next degree of metal insertion, several initial structures were optimized; only the results for the lowest-energy ones were considered for further discussion. The optimized lithiated structures are visualized on Figure 2 and Figure S1 (Supplementary Materials). All electro-negative elements (O, N, Cl) participate in the coordination of Li+.
In the first step of lithiation (two Li atoms), lithium is coordinated only to oxygens from the C=O groups. In the next step (four Li atoms), the chloride counterions are included in the coordination as well. It is noteworthy that in this step, two Li are coordinated by four oxygens—two individual and two shared for each, while the other two Li share the two chloride ions and coordinate to one individual oxygen (Figure 2). The strive of Li+ to be tetracoordinated with four covalently bound oxygens or to be tricoordinated if a halide ion occupies one of the coordination sites has been registered already in an earlier study [49]. Only after the insertion of 6 or more Li atoms all functional groups and counterions are involved in the coordination. The MOF was able to accommodate 24 Li atoms before the energy change became unfavorable.
The inserted Li atoms reduce initially only the organic molecules but upon loading of 6 Li the Ni charge gradually starts to drop (Table 1) and decreases from 1.25 in the empty MOF to 0.24 at 24 Li. For 26 Li atoms the Ni charge becomes close to zero, meaning that it is decoordinated, which is one of the reasons further lithiation is energetically unfavorable.
Li charge varies in a narrow interval of values throughout the insertion of the first 24 Li atoms and sharply decreases thereafter (Table 1, Figure 3). Electron density is transferred from Li mainly to the electronegative atoms—oxygens and nitrogens (Figure 3), in accordance with the pattern of Li+ coordination. Figure 3 shows that O and N bear a similar amount of electron density up to 16 Li and, thereafter, more negative charge is allocated on the nitrogens. The active involvement of nitrogen in charge exchange with Li is the reason for its gradual partial decoordination from the Ni ions—starting at 8 Li, the Ni coordination number drops from 6 to 4 and decreases further down to 2 (Table 1). Parallel to this, close contacts Li-Li (<2.67 Å, the bond length in Li2) are registered, growing in number with lithium loading. These non-ionic clusters bear no charge; this is one of the reasons for the gradual decrease of the average Li charge (Table 1, Figure 3 and Figure 4—bottom). This is an interesting result if we refer to the experimentally observed Li clusters that are formed inside the pores of carbon-based anodes (such as hard carbons, microporous carbon with open porosity) upon overlithiation [50,51]. The energy storage by the formation of Li clusters is of practical importance in order to achieve electrode materials with colossal capacity. In this respect, further studies will be needed to understand whether the Li clustering is a specific feature of MOF-S21.
The effect of Li insertion on the periodic box proportions is presented in Table 2. Metal insertion causes lowering of the lattice symmetry due to variation of all lattice parameters. In more detail, the volume of the box shrinks during the insertion of the first 6 Li atoms, more prominently at the expense of angle reduction, then expands again due to angle growth until 16 Li, and finally contracts again due to both edge and angle decrease. Owing to this ‘breathing’ behavior, the volume oscillates (Table 1) but stays lower than the initial one at any degree of Li loading. At the maximum Li uptake, the volume is the same as for 2Li—~20% less than that of the pristine MOF. Irrespective of this dimension variation, the designed MOF is sufficiently stable to accommodate a huge amount of Li.
The changes in the edge lengths were due to the changes in the geometry of the ligands resulting from the charge redistribution in the system. A characteristic feature of the initial structure of the MOF edges is the bond length alternation (BLA), which is defined in principle as the average difference between the longer and the shorter bonds in a π-conjugated molecule. We have considered two values of BLA: (i) the BLA of the C-C bonds in the cyclic portion of the quinones (Q-M) and (ii) the alternation of the lengths of the triple C≡N and the adjacent simple C-C exocyclic bonds (CN-M). The results are shown in Figure 4 and reveal that, upon Li loading, the quinoid ring gradually becomes aromatic, as the BLA reaches zero at 10 Li, and, further on, becomes negative, i.e., the former shorter bonds become longer and vice versa. On the other hand, the BLA of the exocyclic fragment drops drastically, indicating the tendency towards equalization of the two bonds. As a result, the quinoid character of the cycle is restored to a certain degree but in another direction and the overall effect is shortening of the edges. Alongside the change in geometry, the organic ligand changes slightly its orientation due to partial decoordination from Ni and coordination to Li+. This results in reduction of the Ni-Ni distances and in overall volume decrease.
Based on the above geometric and electronic analysis, it is now possible to understand the reduction properties of the lithiated MOF-S21. Figure 5 shows the calculated electrochemical potential of the material versus Li0/Li+. The MOF reduction by Li proceeds at three distinct potential steps: at ~3.5 V interact the first 2 Li atoms, followed by a potential at 2.8 V (4 Li to 8 Li), and finally a stable potential at 0.4 V is developed (beyond 8 to 24 Li). This ’terrace-wise’ change in the electrochemical potential upon lithiation mirrors the Li charge profile (Figure 5-bottom), which, on its turn, reflects the pattern of Li+ coordination in the MOF. It is interesting to note that a two-step reduction has been experimentally observed upon the lithiation of 2,3-dicyano-o-benzoquinone (o-2CN-BQ): 3.37 and 2.46 V, respectively [52]. The good agreement between the electrochemical potentials calculated for the lithiated MOF and the measured potential of o-2CN-BQ gives further evidence that the MOF reduction is dominated by the 2CN-BQ structural element during the initial degrees of lithiation. In addition, the MOF accommodates extra Li atoms in the third potential step (i.e., at 0.4 V). This is achieved through the activation of the redox properties of the Ni ion in addition to the CN-nitrogens. At 8 Li the formation of Li clusters starts (Figure 4). It is clear that a combination between the redox properties of Ni and 2CN-BQ, as well as the MOF flexible structure, bear promise for achieving a very high lithium capacity.

2.2.2. Sodiation of MOF-S21

The initial placement of the Na atoms followed the same logic as Li—next to oxygens. However, the optimized geometries did not coincide ideally—as early as the first pair of Na atoms, the Cl counterion was already involved in the coordination of Na+ and at 4 Na the CN-nitrogens joined in as well (Figure 6, Figure S2 and S3). Sodium has a flexible coordination number but is definitely higher than 4, which explains the coordination to more electron-acceptors.
The effect of Na insertion on Ni is both similar to and different from Li insertion. The charge of Ni decreases faster until 8Na but keeps almost constant thereafter, even when insertion becomes unfavorable (Table 3). The Na charge maintains a high essentially constant value upon loading of the first 14 Na and only decreases in the last steps (Table 3, Figure 7 and Figure 8—bottom). The electron density is again transferred from Na to oxygen and nitrogen but more uniformly than in the case of Li. Only in the last steps do the nitrogens accept more charge than the oxygens. As a result, the Ni coordination number drops from 6 to 4 and does not decrease further (Table 3). This is the reason why the effects of Ni decoordination on the geometry of the periodic box are different (Table 2) and tend to rectify the structure back to a cubic box. The volume again shrinks during the insertion of the first 6 Na by 17% but then expands and maintains an almost constant value at about 95% of the empty box (Table 3). The BLA values for both Q-Na and CN-Na (Figure 4) have exactly the same conduct as the respective Q-Li and CN-Li—aromatization of the quinoid ring upon insertion of 8–12 Na and inversion of the direction of quinoidization thereafter. However, a major difference from the Li uptake behavior is the absence of clusters—only one close Na-Na contact is found in the last three favorable insertion steps (Table 3).
The electrochemical potential versus Na0/Na+ electrode of each step of sodiation is calculated and shown in Figure 8. As in the case of Li, two features can be outlined: the appearance of two high-voltage steps at 2.9 and 2.4 V at the beginning of the MOF sodiation, followed by a drastic decrease in the potential (by more than 1.0 V) for the last sodiation degrees. The high-voltage steps of MOF-sodiation are shifted in comparison with those corresponding to the MOF lithiation, which is in agreement with the standard redox potentials of sodium and lithium. This allows assigning the high-voltage steps to the preferential reduction of 2CN-BQ ligands by Na. However, the total amount of Na needed to reduce 2CN-BQ is lower than that for Li (i.e., 6 versus 8, Figure 5 and Figure 8). For the sake of comparison, the o-2CN-BQ grafted CNTs have been shown to exhibit lower capacity in a sodium half-cell than that in a lithium one [52]. The observed difference in the reduction of 2CN-BQ ligands by Na and Li is, most probably, associated with their coordination dissimilarity.
However, the MOF is still able to uptake extra Na due to the complementary redox properties of the Ni ion and CN nitrogens. The high degrees of sodiation take place at potentials lower than 1.2 V through several well-discernable steps. The total number of sodium uptake is 16 Na, which is lower than that of Li (i.e., 24). This difference arises from the specific feature of Li to form metal clusters in the reduced MOF. In comparison with Li, it appears that the formation of metal clusters is restricted for Na. This can be explained with the structural changes of more pronounced volume compactization.

2.2.3. Magneziation of MOF-S21

As in the case of Li, Mg displays a selective coordination towards the functional groups and counterions in MOF-S21. The first two magnesium atoms are coordinated solely to C=O oxygens. Upon the introduction of the next two Mg atoms, first the chlorides and then the CN-nitrogens are engaged in the coordination (Figure 9). The maximum number of 6 Mg could be added energetically favorably.
The insertion of the first 4 Mg invoked practically no effect on the charges of Mg and Ni, nor on the coordination number of the latter (Table 4). The next two Mg atoms caused gradual drop in the charges of both metals and the coordination number of Ni, resulting in inefficient further magneziation (Table 4, Figure 10 and Figure 11—bottom). The BLA values reveal that the quinoid ring was aromatized as early as the insertion of 2 Mg and stayed aromatic, while BLA of the exocyclic fragment started decreasing at 4 Mg (Figure 4). Upon Mg insertion the volume of the unit cell underwent volume contraction, followed by a rapid expansion thereafter, most probably due to the partial Ni decoordination (Table 4). The expansion did not go beyond 95% of the initial unit cell volume. Moreover, the averaged periodic box parameters are the closest to those of the empty box, i.e., the unit cell withstands the least deformation upon Mg loading. No aggregation of Mg atoms is registered (Table 4).
As usual, the electron density from Mg went to oxygen and nitrogen; however, unlike Li and Na, the shares of O and N were equal (Figure 10) and kept essentially constant after 2 Mg. Therefore, at 6Mg and beyond, the Mg charge started decreasing and reduction of Ni was witnessed (Table 4, Figure 10).
The Mg charge variation profile is synchronized with the electrochemical potential profile, except upon the addition of 3 Mg (Figure 11). In comparison with Li and Na, the electrochemical curve of Mg consists of two high-potential steps at 2.4 and 2.0 V, followed by a large jump in the potential of over 1.5 V after 3 Mg. The high potentials corresponding to the uptake of 2 and 3 Mg are lower than those calculated for low degrees of MOF lithiation and sodiation, as could be expected from the magnitudes of the standard redox potentials of magnesium, sodium and lithium. This means that the high-voltage steps occurring during the MOF magnesiation are also a consequence of the preferential reduction of the 2CN-BQ ligands. Contrary to Li and Na, the MOF accommodates a limited number of extra Mg (i.e., only 3 Mg) at potentials lower than 1 V.

2.3. Generalized Electrochemical Results

Based on the results for the electrochemical potential, the theoretical specific capacity and energy density are calculated (Table 5). The highest capacity is delivered upon lithiation, while the lowest—upon magnesiation. The same order is obeyed for the energy density. To appraise these results, we referred to data for metal-organic batteries employing Li and Mg metal anodes, and anthraquinone as a cathode: the theoretical energy density is 530 W h/kg for the Li-cell and 345 W h/kg for the Mg-cell [53]. In addition, it has recently been found that a ferrocene-based MOF with composition iron (III) 1,1′-ferrocenedicarboxylate showed a high energy density of 549 W h/kg as a cathode in a lithium-ion battery [54]. The comparison manifests that the designed MOF exhibits twice as high energy density in respect to both Li+ and Mg2+. Moreover, the energy density of MOF-S21 towards Li+ exceeds that of Li-rich oxides, which are now considered as next-generation electrode materials for LIBs (i.e., energy density is higher than 400 W h/kg) [55].
In comparison with Li+ and Mg2+, the capacity and energy density upon MOF-S21 sodiation are intermediate. However, the magnitudes of both capacity and energy density outperform those predicted for the highest-performing layered oxides Na2RhO3 and Na2PdO3 [56]: energy density of 747 and 734 W h/kg, respectively.
The suitability of the designed MOF as an electrode material in Li, Na and Mg-ion batteries can also be evaluated by the chemical and geometrical stability upon metal insertion. MOF-S21 appears to be chemically stable since no covalent bond cleavage is observed at any level of metal uptake. The structural changes invoked by the maximum amount of favorably added Li, Na and Mg into MOF-S21 are compared on Figure 12. Li24-MOF exhibits the largest structural transformation: the unit cell deviates most significantly from the initial rhombohedral cell with a volume change of around 20%. In contrast, the unit cells of both Na16-MOF and Mg6-MOF are closer to the initial one with only 5 and 8% volume decrease, respectively. These are exceptionally small volume changes if we consider the ionic radius of Na+ versus that of Li+, as well as the charge of Mg2+ versus that of Li+. The limited volume changes will generate little stress and strain in the electrode during cycling, thus preventing its performance failure. Therefore, one can expect a better cycling stability of the MOF-S21 during sodiation and magnesiation in comparison with that of lithiation.

3. Models and Methods

All calculations were performed in periodic boundary conditions (PBC) with the Vienna Ab initio Simulation Package (VASP 5.4.4) [57,58,59]. Spin polarized density functional theory (DFT) was employed, specifically the PBE generalized gradient approximation exchange and correlation functionals [60]. The core electrons were described with pseudopotentials using the projector-augmented wave (PAW) method [61,62]; for the outer electrons, a plane-wave basis set with an energy cut-off of 600 eV was utilized. The used cores for the elements were as follows: [He] for Li, C, N and O; [Be] for Na and Mg; [Ne] for Cl; and [Ar] for Ni. The Γ-point was used to sample the Brillouin zone. The electronic partial occupancies were calculated according to the Gaussian smearing scheme with a smearing parameter of 0.05 eV. Geometry optimization, including the unit cell shape and volume, was carried out until all forces acting on the atoms became less than 0.01 eV/Å. The charge distribution was calculated according to Bader’s quantum theory of atoms in molecules (QTAIM) [63] with the Bader program [64]. The VESTA program [65] was used for structure visualization. The models with a small to medium amount of metal atoms were built by testing different chemically sensible interaction sites and proceeding with the most stable structure after optimization. Upon a certain degree of metalation, all intuitively chosen sites became occupied. In order to generate further initial structures, short molecular dynamics simulations were carried out. Computational details, unique to the MD simulations, included: Canonical ensemble; Nose-Hoover thermostat [66,67,68]; integration of the equations of motion according to the Verlet algorithm [69]; time step of 0.5 fs; temperature = 400 K; energy cut-off of 400 eV; 4000 or 6000 MD steps. The energies of bulk metals were calculated using the experimental cell parameters and an 11 × 11 × 11 k-point mesh.
The step-wise energy differences used for the estimation of the potential were calculated according to the equation:
Δ E DFT   =   E ( M y MOF ) [ E ( M x MOF ) + ( y x ) E ( M ) ]
where x and y are the number of metal atoms in two consecutive steps of metal insertion and E(M) is the energy of one atom in the bulk metal.
It has been shown that, for solid systems, the difference between the internal energy and the free energy is negligibly small and the electrode potential can be estimated using the equation [70,71]:
U   [ V ] = Δ E DFT [ e V ] / z
The capacity and the energy density of the electrode material are calculated according to the formulas:
Capacity = z F / ( 3.6 m M O F )
Energy   density = F / ( 3 . 6 m M O F ) U d z
where z is the number of exchanged electrons, F is the Faraday constant, mMOF is the molar mass of the electrode material in g.mol−1.

4. Conclusions

Coordination of Ni(II) with CN-groups from 2,5-dicyano-p-quinone ligands leads to the formation of a rhombohedral MOF, the charge neutralization being achieved by Cl¯ counterions. The designed MOF is able to uptake an extremely high amount of alkali and alkaline earth metals, i.e., 24 Li, 16 Na and 6 Mg. The inserted metal atoms further lower the cell symmetry and convert it into a triclinic one, but in a different way—Li invokes maximum flattening and the smallest volume (about 20% volume drop), whereas Mg, and especially Na, have the opposite effect towards rectifying the angles (and no more than 8 and 5% volume decrease). This could be considered as a guarantee for a good cycling stability during the sodiation and magnesiation of the designed material. In all cases, the addition of metal atoms inflicts compactization of the MOF owed to structural changes in its edges—the quinone molecules. The latter gradually lose their quinoid character caused by carboxyl oxygen reduction and become aromatic. This is accompanied by strengthened delocalization in the exocyclic fragment due to CN reduction, leading to a subsequent requinoidization of the organic molecule in the N-N direction. The latter process results in partial (Na, Mg) or substantial (Li) reduction and decoordination of Ni(II). The charge redistribution and the electrochemical potential profiles confirm this description of the processes occurring upon introduction of Li, Na, and Mg. Each stage of involvement of MOF fragments in the reduction are reflected in the charge and potential profile of the respective added atom. The first quotas of added metals (i.e., 8 Li, 6 Na and 3 Mg) are stored in the MOF at two high-potential steps (i.e., 3.55/2.8 V vs. Li0/Li+, 2.9/2.4 V vs. Na0/Na+ and 2.4/2.0 V vs. Mg0/Mg2+) due to reduction of the 2CN-BQ ligands, while below 1.0 V, the MOF accommodated an extra number of metals through the activation of the redox properties of the Ni ion in addition to the CN-nitrogens. An additional feature of MOF-S21 is the overlithiation due to Li clusters formation.
The designed MOF stores Li+, Na+ and Mg2+ by delivering high specific capacity and energy density that significantly exceeds those typical for Li- and Na-rich layered oxides, as well as the organic anthraquinone-based electrodes. Moreover, MOF-S21 is definitely more environmentally friendly, as the contents of Ni is 9.7 wt%, while, in the conventional layered oxides, the transition metals (Ni, Co) constitute more than 50 wt%. Moreover, as reported recently [72], quinones could be extracted from biomass, which would make the production of quinone-based MOFs a green chemistry project. The predicted excellent electrochemical properties of MOF-S21 regarding Li+, Na+ and Mg2+ deserve to be tested experimentally.

Supplementary Materials

The following are available online, Figure S1: MOF-S21 at different degrees of lithiation; Figure S2: MOF-S21 at different degrees of sodiation. Figure S3: Moderately lithiated and sodiated MOF-S21.

Author Contributions

Conceptualization, A.T. and R.S.; methodology, H.R.; software, H.R. and A.S.; investigation H.R. and A.S.; data curation, H.R. and A.S.; writing—original draft preparation, H.R. and A.S.; writing—review and editing, A.T., H.R. and R.S.; visualization, H.R. and A.S.; project administration, R.S. and A.T.; funding acquisition, R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bulgarian National Science Fund, project CARiM-VIHREN, grant number KP06-DV-6/16.12.2019.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The input and output files from the calculations (including geometries and energies) are available from the corresponding author upon reasonable request.

Acknowledgments

This work was supported by the European Regional Development Fund within the Operational Programme “Science and Education for Smart Growth 2014–2020” under the Project CoE “National center of mechatronics and clean technologies” BG05M2OP001-1.001-0008, and by the CARiM-VIHREN project, KP06-DV-6/16.12.2019. A.S. acknowledges the Erasmus+ support. The authors are indebted to Ilia Kichev and Lyuben Borislavov for providing tools for automated data processing.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Arroyo-de Dompablo, M.; Ponrouch, A.; Johansson, P.; Palacín, M. Achievements, Challenges, and Prospects of Calcium Batteries. Chem. Rev. 2020, 120, 6331–6357. [Google Scholar] [CrossRef]
  2. Elia, G.A.; Marquardt, K.; Hoeppner, K.; Fantini, S.; Lin, R.; Knipping, E.; Peters, W.; Drillet, J.-F.; Passerini, S.; Hahn, R. An Overview and Future Perspectives of Aluminum Batteries. Adv. Mater. 2016, 28, 7564–7579. [Google Scholar] [CrossRef] [PubMed]
  3. Qin, K.; Huang, J.; Holguin, K.; Luo, C. Recent advances in developing organic electrode materials for multivalent rechargeable batteries. Energy Environ. Sci. 2020, 13, 3950–3992. [Google Scholar] [CrossRef]
  4. Tian, Y.; Zeng, G.; Rutt, A.; Shi, T.; Kim, H.; Wang, J.; Koettgen, J.; Sun, Y.; Ouyang, B.; Chen, T.; et al. Promises and Challenges of Next-Generation “Beyond Li-ion” Batteries for Electric Vehicles and Grid Decarbonization. Chem. Rev. 2021, 121, 1623–1669. [Google Scholar] [CrossRef] [PubMed]
  5. Hwang, J.-Y.; Myung, S.-T.; Sun, Y.-K. Sodium-ion batteries: Present and future. Chem. Soc. Rev. 2017, 46, 3529–3614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Muldoon, J.; Bucur, C.B.; Gregory, T. Quest for nonaqueous multivalent secondary batteries: Magnesium and beyond. Chem. Rev. 2014, 114, 11683–11720. [Google Scholar] [CrossRef] [PubMed]
  7. Yabuuchi, N.; Kubota, K.; Dahbi, M.; Komaba, S. Research Development on Sodium-Ion Batteries. Chem. Rev. 2014, 114, 11636–11682. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, F.; Fan, X.; Gao, T.; Sun, W.; Ma, Z.; Yang, C.; Han, F.; Xu, K.; Wang, C. High-Voltage Aqueous Magnesium Ion Batteries. ACS Cent. Sci. 2017, 3, 1121–1128. [Google Scholar] [CrossRef] [Green Version]
  9. Vanysek, P. Electrochemical series. In CRC Handbook of Chemistry and Physics; CRC Press: Boca Raton, FL, USA, 2000. [Google Scholar]
  10. Furukawa, H.; Cordova, K.E.; O’Keeffe, M.; Yaghi, O.M. The Chemistry and Applications of Metal-Organic Frameworks. Science 2013, 341, 1230444. [Google Scholar] [CrossRef] [Green Version]
  11. Zhou, H.-C.; Kitagawa, S. Metal–Organic Frameworks (MOFs). Chem. Soc. Rev. 2014, 43, 5415–5418. [Google Scholar] [CrossRef] [Green Version]
  12. Yaghi, O.M.; O’Keeffe, M.; Ockwig, N.W.; Chae, H.K.; Eddaoudi, M.; Kim, J. Reticular synthesis and the design of new materials. Nature 2003, 423, 705–714. [Google Scholar] [CrossRef]
  13. Lu, W.; Wei, Z.; Gu, Z.-Y.; Liu, T.-F.; Park, J.; Park, J.; Tian, J.; Zhang, M.; Zhang, Q.; Iii, T.G.; et al. Tuning the structure and function of metal–organic frameworks via linker design. Chem. Soc. Rev. 2014, 43, 5561–5593. [Google Scholar] [CrossRef]
  14. Odoh, S.O.; Cramer, C.; Truhlar, D.; Gagliardi, L. Quantum-Chemical Characterization of the Properties and Reactivities of Metal–Organic Frameworks. Chem. Rev. 2015, 115, 6051–6111. [Google Scholar] [CrossRef] [Green Version]
  15. Lu, Y.; Chen, J. Prospects of organic electrode materials for practical lithium batteries. Nat. Rev. Chem. 2020, 4, 127–142. [Google Scholar] [CrossRef]
  16. Lu, Y.; Zhang, Q.; Li, L.; Niu, Z.; Chen, J. Design Strategies toward Enhancing the Performance of Organic Electrode Materials in Metal-Ion Batteries. Chem 2018, 4, 2786–2813. [Google Scholar] [CrossRef] [Green Version]
  17. Esser, B.; Dolhem, F.; Becuwe, M.; Poizot, P.; Vlad, A.; Brandell, D. A perspective on organic electrode materials and technologies for next generation batteries. J. Power Sources 2021, 482, 228814. [Google Scholar] [CrossRef]
  18. Cui, X.; Dong, H.; Chen, S.; Wu, M.; Wang, Y. Progress and Perspective of Metal- and Covalent-Organic Frameworks and their Derivatives for Lithium-Ion Batteries. Batter. Supercaps 2020, 4, 72–97. [Google Scholar] [CrossRef]
  19. Wang, Z.; Tan, R.; Wang, H.; Yang, L.; Hu, J.; Chen, H.; Pan, F. A Metal-Organic-Framework-Based Electrolyte with Nanowetted Interfaces for High-Energy-Density Solid-State Lithium Battery. Adv. Mater. 2018, 30, 1704436. [Google Scholar] [CrossRef]
  20. Shrivastav, V.; Sundriyal, S.; Goel, P.; Kaur, H.; Tuteja, S.K.; Vikrant, K.; Kim, K.-H.; Tiwari, U.K.; Deep, A. Metal-organic frameworks (MOFs) and their composites as electrodes for lithium battery applications: Novel means for alternative energy storage. Coord. Chem. Rev. 2019, 393, 48–78. [Google Scholar] [CrossRef]
  21. Jaschin, P.W.; Gao, Y.; Li, Y.; Bo, S.-H. A materials perspective on magnesium-ion-based solid-state electrolytes. J. Mater. Chem. A 2020, 8, 2875–2897. [Google Scholar] [CrossRef]
  22. Calbo, J.; Golomb, M.J.; Walsh, A. Redox-active metal–organic frameworks for energy conversion and storage. J. Mater. Chem. A 2019, 7, 16571–16597. [Google Scholar] [CrossRef]
  23. Zhang, Z.; Yoshikawa, H.; Awaga, K. Monitoring the Solid-State Electrochemistry of Cu(2,7-AQDC) (AQDC = Anthraquinone Dicarboxylate) in a Lithium Battery: Coexistence of Metal and Ligand Redox Activities in a Metal–Organic Framework. J. Am. Chem. Soc. 2014, 136, 16112–16115. [Google Scholar] [CrossRef]
  24. Baumann, A.E.; Burns, D.A.; Liu, B.; Thoi, V.S. Metal-organic framework functionalization and design strategies for advanced electrochemical energy storage devices. Commun. Chem. 2019, 2, 86. [Google Scholar] [CrossRef] [Green Version]
  25. Li, C.; Wang, K.; Li, J.Z.; Zhang, Q. Nanostructured potassium–organic framework as an effective anode for potassium-ion batteries with a long cycle life. Nanoscale 2020, 12, 7870–7874. [Google Scholar] [CrossRef] [PubMed]
  26. Férey, G.; Millange, F.; Morcrette, M.; Serre, C.; Doublet, M.-L.; Grenèche, J.-M.; Tarascon, J.-M. Mixed-Valence Li/Fe-Based Metal–Organic Frameworks with Both Reversible Redox and Sorption Properties. Angew. Chem. Int. Ed. 2007, 46, 3259–3263. [Google Scholar] [CrossRef]
  27. Luo, D.; Lei, P.; Tian, G.; Huang, Y.; Ren, X.; Xiang, X. Insight into Electrochemical Properties and Reaction Mechanism of a Cobalt-Rich Prussian Blue Analogue Cathode in a NaSO3CF3 Electrolyte for Aqueous Sodium-Ion Batteries. J. Phys. Chem. C 2020, 124, 5958–5965. [Google Scholar] [CrossRef]
  28. Schon, T.B.; McAllister, B.T.; Li, P.-F.; Seferos, D.S. The rise of organic electrode materials for energy storage. Chem. Soc. Rev. 2016, 45, 6345–6404. [Google Scholar] [CrossRef] [Green Version]
  29. Huynh, M.; Anson, C.; Cavell, A.C.; Stahl, S.S.; Hammes-Schiffer, S. Quinone 1 e and 2 e/2 H+ Reduction Potentials: Identification and Analysis of Deviations from Systematic Scaling Relationships. J. Am. Chem. Soc. 2016, 138, 15903–15910. [Google Scholar] [CrossRef]
  30. Araujo, R.B.; Banerjee, A.; Panigrahi, P.; Yang, L.; Strømme, M.; Sjödin, M.; Araujo, C.M.; Ahuja, R. Designing strategies to tune reduction potential of organic molecules for sustainable high capacity battery application. J. Mater. Chem. A 2017, 5, 4430–4454. [Google Scholar] [CrossRef]
  31. Emanuelsson, R.; Sterby, M.; Strømme, M.; Sjödin, M. An All-Organic Proton Battery. J. Am. Chem. Soc. 2017, 139, 4828–4834. [Google Scholar] [CrossRef] [PubMed]
  32. Miao, L.; Liu, L.; Shang, Z.; Li, Y.; Lu, Y.; Cheng, F.; Chen, J. The structure–electrochemical property relationship of quinone electrodes for lithium-ion batteries. Phys. Chem. Chem. Phys. 2018, 20, 13478–13484. [Google Scholar] [CrossRef] [PubMed]
  33. Oka, K.; Strietzel, C.; Emanuelsson, R.; Nishide, H.; Oyaizu, K.; Strømme, M.; Sjödin, M. Characterization of PEDOT-Quinone conducting redox polymers in water-in-salt electrolytes for safe and high-energy Li-ion batteries. Electrochem. Commun. 2019, 105, 106489. [Google Scholar] [CrossRef]
  34. Kim, K.C.; Liu, T.; Jung, K.H.; Lee, S.W.; Jang, S.S. Unveiled correlations between electron affinity and solvation in redox potential of quinone-based sodium-ion batteries. Energy Storage Mater. 2019, 19, 242–250. [Google Scholar] [CrossRef]
  35. Rasheev, H.G.; De Araújo, R.B.N.; Tadjer, A.; Johansson, P. Fundamental promise of anthraquinone functionalized graphene based next generation battery electrodes: A DFT study. J. Mater. Chem. A 2020, 8, 14152–14161. [Google Scholar] [CrossRef]
  36. Jung, K.H.; Jeong, G.; Go, C.Y.; Kim, K.C. Conjugacy of organic cathode materials for high-potential lithium-ion batteries: Carbonitriles versus quinones. Energy Storage Mater. 2019, 24, 237–246. [Google Scholar] [CrossRef]
  37. An, S.Y.; Schon, T.B.; Seferos, D.S. Stable, Dual Redox Unit Organic Electrodes. ACS Omega 2020, 5, 1134–1141. [Google Scholar] [CrossRef] [Green Version]
  38. Peng, Z.; Yi, X.; Liu, Z.X.; Shang, J.; Wang, D. Triphenylamine-Based Metal–Organic Frameworks as Cathode Materials in Lithium-Ion Batteries with Coexistence of Redox Active Sites, High Working Voltage, and High Rate Stability. ACS Appl. Mater. Interfaces 2016, 8, 14578–14585. [Google Scholar] [CrossRef]
  39. Gu, S.; Bai, Z.; Majumder, S.; Huang, B.; Chen, G. Conductive metal–organic framework with redox metal center as cathode for high rate performance lithium ion battery. J. Power Sources 2019, 429, 22–29. [Google Scholar] [CrossRef]
  40. Hmadeh, M.; Liu, Z.; Gándara, F.; Furukawa, H.; Wan, S.; Augustyn, V.; Chang, R.; Liao, L.; Zhou, F.; Perre, E.; et al. New Porous Crystals of Extended Metal-Catecholates. Chem. Mater. 2012, 24, 3511–3513. [Google Scholar] [CrossRef]
  41. Wu, Z.; Adekoya, D.; Huang, X.; Kiefel, M.J.; Xie, J.; Xu, W.; Zhang, Q.; Zhu, D.; Zhang, S. Highly Conductive Two-Dimensional Metal–Organic Frameworks for Resilient Lithium Storage with Superb Rate Capability. ACS Nano 2020, 14, 12016–12026. [Google Scholar] [CrossRef]
  42. Hameed, A.S.; Nagarathinam, M.; Schreyer, M.; Reddy, M.V.; Chowdari, B.V.R.; Vittal, J.J. A layered oxalatophosphate framework as a cathode material for Li-ion batteries. J. Mater. Chem. A 2013, 1, 5721–5726. [Google Scholar] [CrossRef]
  43. Wessells, C.D.; Peddada, S.V.; Huggins, R.A.; Cui, Y. Nickel Hexacyanoferrate Nanoparticle Electrodes for Aqueous Sodium and Potassium Ion Batteries. Nano Lett. 2011, 11, 5421–5425. [Google Scholar] [CrossRef]
  44. Li, X.; Cheng, F.; Zhang, S.; Chen, J. Shape-controlled synthesis and lithium-storage study of metal-organic frameworks Zn4O(1,3,5-benzenetribenzoate)2. J. Power Sources 2006, 160, 542–547. [Google Scholar] [CrossRef]
  45. Ogihara, N.; Yasuda, T.; Kishida, Y.; Ohsuna, T.; Miyamoto, K.; Ohba, N. Organic Dicarboxylate Negative Electrode Materials with Remarkably Small Strain for High-Voltage Bipolar Batteries. Angew. Chem. Int. Ed. 2014, 53, 11467–11472. [Google Scholar] [CrossRef] [PubMed]
  46. Hameed, S.; Petnikota, S.; Hassan, N.; Al-Qaradawi, S.; Karim, Z.; Reddy, M. Synthesis of Nickel Fumarate and Its Electrochemical Properties for Li-Ion Batteries. Electrochem 2021, 2, 439–451. [Google Scholar] [CrossRef]
  47. Lee, H.H.; Park, Y.; Kim, S.H.; Yeon, S.-H.; Kwak, S.K.; Lee, K.T.; Hong, S.Y. Mechanistic Studies of Transition Metal-Terephthalate Coordination Complexes upon Electrochemical Lithiation and Delithiation. Adv. Funct. Mater. 2015, 25, 4859–4866. [Google Scholar] [CrossRef]
  48. Lefebvre, J.; Chartrand, D.; Leznoff, D.B. Synthesis, structure and magnetic properties of 2-D and 3-D [cation]{M[Au(CN)2]3} (M=Ni,Co) coordination polymers. Polyhedron 2007, 26, 2189–2199. [Google Scholar] [CrossRef]
  49. Rasheev, H.; Stoyanova, R.; Tadjer, A. Rivalry at the Interface: Ion Desolvation and Electrolyte Degradation in Model Ethylene Carbonate Complexes of Li+, Na+, and Mg2+ with PF6 on the Li4Ti5O12 (111) Surface. ACS Omega 2021, 6, 29735–29745. [Google Scholar] [CrossRef]
  50. Tian, R.; Duan, H.; Guo, Y.; Li, H.; Liu, H. High-Coulombic-Efficiency Carbon/Li Clusters Composite Anode without Precycling or Prelithiation. Small 2018, 14, e1802226. [Google Scholar] [CrossRef] [PubMed]
  51. Bloi, L.M.; Hippauf, F.; Boenke, T.; Rauche, M.; Paasch, S.; Schutjajew, K.; Pampel, J.; Schwotzer, F.; Dörfler, S.; Althues, H.; et al. Mechanistic insights into the reversible lithium storage in an open porous carbon via metal cluster formation in all solid-state batteries. Carbon 2022, 188, 325–335. [Google Scholar] [CrossRef]
  52. Lyu, H.; Jafta, C.J.; Popovs, I.; Meyer, H.M.; Hachtel, J.A.; Huang, J.; Sumpter, B.G.; Dai, S.; Sun, X.-G. A dicyanobenzoquinone based cathode material for rechargeable lithium and sodium ion batteries. J. Mater. Chem. A 2019, 7, 17888–17895. [Google Scholar] [CrossRef]
  53. Bitenc, J.; Pavčnik, T.; Košir, U.; Pirnat, K. Quinone Based Materials as Renewable High Energy Density Cathode Materials for Rechargeable Magnesium Batteries. Materials 2020, 13, 506. [Google Scholar] [CrossRef] [Green Version]
  54. Li, C.; Zhang, C.; Xie, J.; Wang, K.; Li, J.; Zhang, Q. Ferrocene-based metal-organic framework as a promising cathode in lithium-ion battery. Chem. Eng. J. 2020, 404, 126463. [Google Scholar] [CrossRef]
  55. Cui, S.; Gao, M.; Li, G.; Gao, X. Insights into Li-Rich Mn-Based Cathode Materials with High Capacity: From Dimension to Lattice to Atom. Adv. Energy Mater. 2021, 2003885. [Google Scholar] [CrossRef]
  56. Do, J.; Kim, I.; Kim, H.; Jung, Y. Towards stable Na-rich layered transition metal oxides for high energy density sodium-ion batteries. Energy Storage Mater. 2020, 25, 62–69. [Google Scholar] [CrossRef]
  57. Kresse, G.; Hafner, J. Ab initiomolecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558–561. [Google Scholar] [CrossRef] [PubMed]
  58. Kresse, G.; Furthmüller, J. Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  59. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  60. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  61. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [Green Version]
  62. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758. [Google Scholar] [CrossRef]
  63. Bader, R.F.W. Atoms in molecules. Acc. Chem. Res. 1985, 18, 9–15. [Google Scholar] [CrossRef]
  64. Tang, W.; Sanville, E.; Henkelman, G. A grid-based Bader analysis algorithm without lattice bias. J. Phys. Condens. Matter 2009, 21, 084204. [Google Scholar] [CrossRef] [PubMed]
  65. Momma, K.; Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  66. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef] [Green Version]
  67. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Nosé, S. Constant Temperature Molecular Dynamics Methods. Prog. Theor. Phys. Suppl. 1991, 103, 1–46. [Google Scholar] [CrossRef] [Green Version]
  69. Verlet, L. Computer “Experiments” on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules. Phys. Rev. (Ser. I) 1967, 159, 98–103. [Google Scholar] [CrossRef]
  70. Aydinol, M.; Kohan, A.F.; Ceder, G.; Cho, K.; Joannopoulos, J. Ab initiostudy of lithium intercalation in metal oxides and metal dichalcogenides. Phys. Rev. B 1997, 56, 1354–1365. [Google Scholar] [CrossRef] [Green Version]
  71. Araujo, R.B.; Banerjee, A.; Ahuja, R. Divulging the Hidden Capacity and Sodiation Kinetics of NaxC6Cl4O2: A High Voltage Organic Cathode for Sodium Rechargeable Batteries. J. Phys. Chem. C 2017, 121, 14027–14036. [Google Scholar] [CrossRef]
  72. Lambert, F.; Danten, Y.; Gatti, C.; Frayret, C. A tool for deciphering the redox potential ranking of organic compounds: A case study of biomass-extracted quinones for sustainable energy. Phys. Chem. Chem. Phys. 2020, 22, 20212–20226. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Construction of the unit cell of MOF-S21 and 3D (3 × 3 × 3) multiple cell representation of the optimized MOF. Color code: brown—C; ecru—H; red—O; blue—N; green—Cl; grey—Ni.
Figure 1. Construction of the unit cell of MOF-S21 and 3D (3 × 3 × 3) multiple cell representation of the optimized MOF. Color code: brown—C; ecru—H; red—O; blue—N; green—Cl; grey—Ni.
Molecules 27 00586 g001
Figure 2. The MOF-S21 unit cell early stages of lithiation. Color code pink—Li, the remaining colors as in Figure 1.
Figure 2. The MOF-S21 unit cell early stages of lithiation. Color code pink—Li, the remaining colors as in Figure 1.
Molecules 27 00586 g002
Figure 3. Average AIM charges upon lithiation of MOF-S21.
Figure 3. Average AIM charges upon lithiation of MOF-S21.
Molecules 27 00586 g003
Figure 4. Bond length alternation (BLA) of the C-C bonds in the benzene ring (Q-M) and of the exocyclic C-C≡N bonds (CN-M) upon metal insertion.
Figure 4. Bond length alternation (BLA) of the C-C bonds in the benzene ring (Q-M) and of the exocyclic C-C≡N bonds (CN-M) upon metal insertion.
Molecules 27 00586 g004
Figure 5. Calculated electrode potential profile (U vs. Li0/Li+) compared to Li charge variation as a function of lithiation.
Figure 5. Calculated electrode potential profile (U vs. Li0/Li+) compared to Li charge variation as a function of lithiation.
Molecules 27 00586 g005
Figure 6. The MOF-S21 unit cell early stages of sodiation. Color code: violet—Na, the remaining colors as in Figure 1.
Figure 6. The MOF-S21 unit cell early stages of sodiation. Color code: violet—Na, the remaining colors as in Figure 1.
Molecules 27 00586 g006
Figure 7. Average AIM charges upon sodiation of MOF-S21.
Figure 7. Average AIM charges upon sodiation of MOF-S21.
Molecules 27 00586 g007
Figure 8. Calculated electrode potential profile (U vs. Na0/Na+) compared to Na charge variation as a function of sodiation.
Figure 8. Calculated electrode potential profile (U vs. Na0/Na+) compared to Na charge variation as a function of sodiation.
Molecules 27 00586 g008
Figure 9. The MOF-S21 unit cell early stages of magnesiation. Color code: yellow—Mg, the remaining colors as in Figure 1.
Figure 9. The MOF-S21 unit cell early stages of magnesiation. Color code: yellow—Mg, the remaining colors as in Figure 1.
Molecules 27 00586 g009
Figure 10. Average AIM charges upon magnesiation of MOF-S21.
Figure 10. Average AIM charges upon magnesiation of MOF-S21.
Molecules 27 00586 g010
Figure 11. Calculated electrode potential profile (vs. Mg0/Mg2+) compared to Mg charge variation as a function of magnesiation.
Figure 11. Calculated electrode potential profile (vs. Mg0/Mg2+) compared to Mg charge variation as a function of magnesiation.
Molecules 27 00586 g011
Figure 12. Multiple cell representation of the maximum exergonic uptake of MOF-S21 with the respective metal atoms.
Figure 12. Multiple cell representation of the maximum exergonic uptake of MOF-S21 with the respective metal atoms.
Molecules 27 00586 g012
Table 1. AIM metal charges and some structural descriptors of the lithiated MOF-S21.
Table 1. AIM metal charges and some structural descriptors of the lithiated MOF-S21.
n(Li)qave(Li)q(Ni)Volume of Unit Cell, ųCoord. # of Ni R(CN…Ni) < 2.3 Å# Li-Li Bonds R(Li-Li) < 2.67 Å
01.251024.86
20.881.21795.360
40.871.21753.960
80.870.75763.945
180.860.32770.737
240.840.24794.1210
260.750.05872.9213
Table 2. Structural parameters of the periodic cell upon maximum exergonic metal insertion.
Table 2. Structural parameters of the periodic cell upon maximum exergonic metal insertion.
Metala, Åb, Åc, Åα, °β, °γ, °
011.6711.6711.6756.1656.1656.16
24 Li10.7911.7712.0850.2743.8554.06
16 Na9.1911.1510.6491.7768.65108.54
6 Mg10.0810.8311.1561.9261.1565.21
Table 3. AIM metal charges and some structural descriptors of the sodiated MOF-S21.
Table 3. AIM metal charges and some structural descriptors of the sodiated MOF-S21.
n(Na)qavr(Na)Volume of Unit Cell, ųq(Ni)Coord. # of Ni
R(CN…Ni) < 2.3 Å
# Na-Na Bonds
R(Na-Na) < 3.08 Å
01024.81.256
40.88869.01.1860
60.88851.30.9440
80.88946.30.6140
100.88963.30.5140
120.86985.40.4641
140.86993.30.4441
160.76959.40.4341
180.69974.00.4041
Table 4. AIM metal charges and some structural descriptors of the magnesiated MOF-S21.
Table 4. AIM metal charges and some structural descriptors of the magnesiated MOF-S21.
n(Mg)qave(Mg)Volume of Unit Cell, ųq(Ni)Coord. # of Ni
R(CN…Ni) < 2.3 Å
# Mg-Mg
R(Mg-Mg) < 2.85 Å
0-1024.81.256
21.77944.41.2460
31.71802.81.2460
41.71770.60.8250
51.69750.10.7940
61.47910.90.7440
81.30950.60.1420
Table 5. Electrochemical properties of MOF-S21 as an electrode material for different charge carriers. Capacity in mA h g−1 and specific energy in W h/kg.
Table 5. Electrochemical properties of MOF-S21 as an electrode material for different charge carriers. Capacity in mA h g−1 and specific energy in W h/kg.
Charge CarrierMaximum CapacityGravimetric Energy Density
24 Li10651343
16 Na7101025
6 Mg533753
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rasheev, H.; Seremak, A.; Stoyanova, R.; Tadjer, A. Redox Hyperactive MOF for Li+, Na+ and Mg2+ Storage. Molecules 2022, 27, 586. https://doi.org/10.3390/molecules27030586

AMA Style

Rasheev H, Seremak A, Stoyanova R, Tadjer A. Redox Hyperactive MOF for Li+, Na+ and Mg2+ Storage. Molecules. 2022; 27(3):586. https://doi.org/10.3390/molecules27030586

Chicago/Turabian Style

Rasheev, Hristo, Agnieszka Seremak, Radostina Stoyanova, and Alia Tadjer. 2022. "Redox Hyperactive MOF for Li+, Na+ and Mg2+ Storage" Molecules 27, no. 3: 586. https://doi.org/10.3390/molecules27030586

APA Style

Rasheev, H., Seremak, A., Stoyanova, R., & Tadjer, A. (2022). Redox Hyperactive MOF for Li+, Na+ and Mg2+ Storage. Molecules, 27(3), 586. https://doi.org/10.3390/molecules27030586

Article Metrics

Back to TopTop