First-principles and machine-learning approaches for interpreting and predicting the properties of MXenes

First-principles and machine-learning approaches for interpreting and predicting the properties of MXenes

Introduction

Over the past decade, we have been witnessing a transformative wave in the search of novel materials, mainly propelled by an explosion of interest in two-dimensional (2D) materials after the discovery of graphene, their exceptional properties and many potential applications. Researchers worldwide have been investigating the unique properties and applications of these ultrathin structures, notably, of a family of two-dimensional transition metal carbides, nitrides, and carbonitrides, named MXenes, which were discovered in 20111.

The first synthesis of an MXene considered the room temperature exfoliation of a powder of Ti3AlC2 MAX phase with hydrofluoric acid (HF), followed by ultrasonication of the reaction products in methanol, originating a material composed of Ti3C2 layers terminated by fluorine and hydroxyl groups2. This work had a significant impact since many tens of MAX phases were already known, suggesting that several other 2D crystals could be obtained using a similar synthetic approach.

The MAX phases are hexagonal layered ternary transition metal nitrides and/or carbides, with a Mn+1AXn formula, where n = 1–4, where M corresponds to an early transition metal (e.g., from groups 3–6 in the Periodic Table), A stands for an element from a subset of groups 12–16, but mostly from groups 13 and 14, and X corresponds to either carbon (C) and/or nitrogen (N). Traditionally, an etchant is used to selectively remove the A-layer elements, which contributes to the weakening of the interactions between the Mn+1Xn layers, allowing them to be readily separated. However, less hazardous (HF-free) synthetic pathways were proposed for obtaining MXenes, such as the NaOH‐assisted hydrothermal process3, or the electrochemical etching with diluted HCl4 or NH4Cl/TMA-OH (TMA = tetramethylammonium)5. The general formula of MXenes is Mn+1XnTx, where M and X are the same as in the MAX phases, and T stands for the surface termination, usually determined by the synthesis conditions (e.g., -OH, -O, or -F are very common). The subscript x in Tx is used to denote that the surface is not necessarily uniformly or fully functionalized. Heating treatments have been devised to remove the hydroxyl and fluoro surface terminations, followed by an H2 treatment to remove the oxo groups on the surface, yielding bare MXenes (i.e., without surface terminations and formula Mn+1Xn)6,7. Kamysbayev et al.8 successfully defunctionalized the Ti3C2Br2 and Ti2CBr2 to their bare counterparts, Ti3C2 and Ti2C, respectively, using LiH at 300 oC. Please note that while it is possible to attain MXenes experimentally, they cannot exist outside a vacuum environment, which severely limits its application.

In 2017, layered ternary transition metal borides (MAB phases) were used to synthesize a family of 2D materials analogous to that of MXenes, the so-called MBenes9, that rely on layers of boron rather than on layers of carbon/nitrogen.

MXenes have garnered significant attention for their exceptional electronic, mechanical, and thermal properties. They display high electrical and metallic conductivity, high oxidation resistance, and hydrophilicity that can be tuned upon hybridization with other materials which, when combined with concomitant large surface areas and versatile structures, make them very attractive for potential applications in various fields, e.g., energy storage, catalysis, etc. The interest in MXenes, according to Google Trends (accessed November 19, 2024), is shown in Fig. 1. The interest has been growing exponentially by 34% per year, with a correlation coefficient ({R}^{2}=0.97).

Fig. 1: Relative interest in MXenes over time.
First-principles and machine-learning approaches for interpreting and predicting the properties of MXenes

Exponential relative interest in MXenes over time (blue dots, one per month). Data obtained from Google Trends.

Full size image

MXenes have become focal points in both experimental and computational studies. So much is known about MXenes and so much research has been going on that there are already countless reviews on MXenes, focused on very specific areas, such as biomedical applications10, electromagnetic interference shielding11,12, quantum dots13, or polymer nanocomposites14, to name a few very recent ones. However, reviews focused on theoretical advances only are rare. To the extent of our knowledge, three of the most recent such reviews15,16,17 deal solely with electrocatalysis and energy applications. In contrast, there are countless reviews on MXenes as catalysts, with little or no focus on theoretical advances in this area18,19,20,21,22,23,24.

In this review, we aim to provide a more general view of how computational techniques can contribute to predict or explain the general properties and applicability of MXenes, both the ones that have already been obtained experimentally and those that have not yet been synthesized, and to motivate the synthesis of new MXenes. In particular, we focus on the intersection between density functional theory (DFT) and machine learning, and on the potential of these approaches to unravel the intricacies of MXenes. As we delve into the geometric and electronic structure, and functionalization of MXenes, a comprehensive understanding of their behavior at the atomic and molecular levels emerges. These computational approaches are shown to be essential to accelerate the pace of 2D material research, providing crucial insights into MXene properties and in deciphering complex relationships within MXene systems, facilitating the design of materials with tailored properties and guiding experimental efforts.

Methods based on density functional theory

Computational methods based on the DFT circumvent approximately solving the Schrödinger equation of a quantum system, by directly obtaining its electron density and, from it, extracting information on physical quantities of interest for the ground state of the system. These include relative energies between different system configurations, bond lengths and angles, atomic charges, or the band structure. Moreover, the first and second derivatives of the energy with respect to the position of the atoms can be used to find special points (minima or N-order saddle-points) at the potential energy surfaces, which are essential for obtaining a complete picture of the energy landscapes of chemical reactions.

In DFT, the poly-electronic wavefunction composed by 3N spatial and N spin variables is replaced by the electron density consisting of three spatial variables only. This critical step is anchored on two theorems by Hohenberg and Kohn25, one stating that the ground state for any system of interacting particles in an external potential, with a given fixed inter-particle interaction, is a unique functional of the electron density and another stating that the ground state energy of the system is the overall minimum value of this functional, i.e., the ground state energy is obtained through minimization of the energy functional with respect to the electron density, hence yielding the energy and the density of the ground state. However, despite considerable efforts since then, the exact form of the energy functional is not known and approximations had to be introduced. One year later, Kohn and Sham25 proposed to represent the many-body problem of interactions in a poly-electronic system into a problem of non-interacting electrons evolving in an external potential. This strategy enables to flyover the terms that are very difficult to calculate, corresponding to kinetic and electronic repulsion, through one-electron Kohn-Sham equations that yield the energies of the corresponding Kohn–Sham orbitals, which are further employed to obtain the density for an N-particle system. The total energy of the system is obtained upon minimization, with respect to the electron density, of the sum of the Kohn-Sham (non-interacting) kinetic energy, electron–electron interaction energy (or Coulomb energy), and the exchange-correlation energy. The latter combines the effects of electron exchange and correlation, i.e., the residual part of the true kinetic energy (because of the approximation of a non-interacting system) and the non-classical electrostatic contributions. In other words, the exchange-correlation energy includes all that we do not know how to calculate exactly and that must be approximated somehow. It is precisely the complexity of the form that is used to approximate the calculation of the exchange-correlation energy that leads to the classification of methods based on the density functional theory in different classes.

Depending on the complexity (and accuracy) of the way the exchange-correlation term is approximated, Perdew et al. proposed a conceptual hierarchy of density functionals, the so-called “Jacob’s ladder of density functional approximations for the exchange-correlation energy”26, which starts at the Hartree world and ascends in terms of sophistication and accuracy from lower to higher rungs.

  1. i.

    Local Spin Density Approximation (LSDA):

  • The simplest form of exchange-correlation functional.

  • Approximates the exchange-correlation energy density at each point in space based solely on the local electron spin density.

  • Adequate to treat densities that vary slowly over space but underestimates certain properties due to neglecting gradient corrections (suitable for solids but not for molecules).

  1. ii.

    Generalized Gradient Approximation (GGA):

  • An advancement over LSDA that considers not only the electron spin density at each point but also the gradient of the density.

  • Incorporates corrections based on the density gradient, improving accuracy in describing molecular geometries and some properties.

  • A widely used example is the PBE functional27.

  1. iii.

    Meta-Generalized Gradient Approximation (meta-GGA):

  • Goes beyond GGA by including the kinetic energy density in addition to the electron density and its gradient.

  • Provides more accurate descriptions of molecular properties such as dispersion interactions and reaction barriers compared to GGAs.

  • Some examples are TPSS or the M06 family of functionals28,29.

  1. iv.

    Hybrid Functionals:

  1. i.

    Global Hybrid Functionals

  • Combine the GGA or meta-GGA exchange-correlation functionals with a fraction of exact Hartree-Fock exchange.

  • Offer improved descriptions of electronic structure, energetics, and band gaps compared to pure GGA or meta-GGA functionals.

  • Examples include the B3LYP30,31 and PBE032,33 functionals.

  1. ii.

    Range-Separated and Screened Hybrid Functionals:

  • Treat short-range and long-range electron-electron interactions separately.

  • Offer better descriptions of charge-transfer excitations and solid-state properties.

  • Some examples are ωB97X34, CAM-B3LYP35, HSE0636 and LC-PBE37.

  1. v.

    Double-Hybrid Functionals:

  • Include functionals that use not only occupied Kohn-Sham orbitals but also virtual orbitals, via methods such as MP2 or the random phase approximation (RPA).

  • Improve with unprecedented accuracy some chemical or physical properties of large-gap compounds, e.g., reaction energies, reaction barrier heights, and nonbonded interactions.

  • Examples include B2PLYP38 and PWPB9539.

The development of hybrid and double-hybrid functionals is part of an effort to improve the accuracy of calculations for molecular and solid-state systems by addressing the limitations of earlier functionals. Hybrid functionals have adjustable parameters and were introduced to mitigate self-interaction and improve exchange-correlation accuracy, while the double-hybrid functionals were developed to enhance accuracy when applied to complex systems with short- and long-range correlation effects. The choice of density functional often depends on the system being studied and the properties of interest. While higher-rung functionals tend to be more accurate, they also come with increased computational cost, making them less practical for large systems. In fact, most studies published in the literature consider semi-local methods, either functionals from rungs 2 and 3 when combined with plane-waves, or hybrid approaches from rung 4 when combined with atomic orbitals. Moreover, either with plane-waves or atomic orbitals, approaches that are capable of accurately describing dispersion are often combined with the functionals for describing non-covalent interactions. A notable example is the atom–atom pairwise empirical dispersion corrections by Grimme et al.40,41,42 which found widespread application. Also, for overcoming the intrinsic self-interaction errors associated with e.g., the functionals from rung 2 and affecting the correct description of localized electrons in many transition metal oxides as CeO2 or NiO, another empirical approach43 using the so-called Hubbard term (+U) that penalizes double occupancy of d or f orbitals has been widely employed. We performed an analysis of the articles published in 2023, where DFT methods were employed to study MXenes, and we found that nine in ten studies considered the Perdew-Burke-Ernzerhof (PBE) functional (3/5 of them augmented with Grimme et al. dispersion corrections). The remaining works considered functionals from different rungs of the ladder such as PZ44, HSE0636, TB-mBJ45, wB97xD46, or B97M-V47.

In summary, different DFT functionals may be used to perform a computer simulation but they must be chosen wisely depending on the property of interest because of their intrinsic limitations. By being based on first principles, DFT provides results at the electronic level, which is both a strength and a limitation of the method. This is because the finer the level of detail of a computational model, the smaller system sizes can be successfully simulated in a reasonable amount of time. Indeed, DFT provides incredibly reliable results for quantum systems, at the cost of only being able to simulate a few thousand atoms with the algorithms and computational capabilities of today, when compared with the millions of particles that can be studied, for instance, with force field-based approaches48. Still, these approaches have been proven key in studies of the properties of materials, such as localized defects, adsorption and catalysis. This information is not only useful by itself, but it can also be fed to higher-level approaches, so that they do not need to go into the electronic detail and can therefore simulate much larger systems. In the majority of the DFT studies reported in the literature, the DFT calculations are performed at T = 0 K, i.e., the effects of the temperature are not considered. The inclusion of the temperature may be accomplished by combining the classical Newtonian equations to describe the motion of the atomic nuclei in the system with the DFT calculated forces acting on the atomic nuclei under specific initial conditions (positions and velocities of the nuclei) and boundary conditions to study how the dynamics of the molecular system evolves with time. This approach constitutes the ab initio molecular dynamics method, which is very time-consuming. Alternatively, standard statistical thermodynamics combined with data from a second-derivative analysis of the system can be used to estimate energies at specific pressure and temperature conditions, significantly reducing the computational cost.

No matter how much technology advances, computational resources are always finite and the calculations cannot be indefinitely parallelized. Therefore, any computational simulation aims at obtaining results as accurate as possible at the lowest computational cost possible. In the case of quantum mechanical simulations, several parameters can be tuned to achieve a compromise between accuracy and cost. The two most relevant ones are the size of the basis set used to expand the Kohn-Sham orbitals and the density of the grid of k-points used for reciprocal space integration49. The higher these two quantities are, the more accurate and resource-consuming the results become. Therefore, one of the first steps in a first-principles study is to find the minimum values of these parameters that yield converged results for a test calculation, i.e., to find the point beyond which the increased computational cost does not lead to significant changes in the results. Other parameters of interest include energy and force convergence criteria, and performing spin-polarized calculations or not50.

Methods based on machine learning

Machine learning (ML) allows to develop statistical models that enable computers to learn from and analyze existing data, thereby identifying patterns and making predictions51,52. Based on data obtained experimentally or by DFT calculations, ML has proven instrumental in predicting MXene properties and proposing new materials tailored for specific uses, as discussed throughout this review. The integration of ML in materials science research not only facilitate the prediction of material suitability for specific applications, but also provides statistical insight into key material features, elucidating their aptness for particular uses and potentially reducing experimental time and costs, particularly when a large number of predictions have to be made in order to explore a broader portion of the materials’ chemical space.

To support such advances, several databases have been developed to provide comprehensive datasets for MXenes. Some examples are: i) the aNANt MXene database53,54, hosted by the Indian Institute of Science, which contains information on over 23,000 MXene structures, including calculated properties; ii) the MXenes Database by MEM-CES55, which includes 730 structural data files (cif format) for MXenes and aims to provide simulated adsorbent metrics of MXenes for CO2/H2 separations; and iii) the MXene-db56,57, which offers structural and computational band gap data for over 4000 MXene structures. These databases serve as foundational tools for ML-driven exploration and design of MXene materials, facilitating deeper insights and accelerating discovery processes.

The development of an ML model involves the following steps58:

  1. i.

    Data preparation: Machine learning models require well-prepared data for calibration and performance evaluation, often sourced through data mining or computational simulations. These data must be meticulously curated, involving processes such as loading, processing, cleaning, and refining. Key challenges in materials science include addressing missing values, correcting inaccuracies, and standardizing units. Subsequently, conducting an exploratory data analysis is crucial to uncover underlying correlations, identify anomalies, and test hypotheses, typically involving a combination of summary statistics and visual representations. At the end of the data treatment, the obtained dataset should be split into three disjoint groups: training, validation, and testing. The training set is employed to calibrate the internal model parameters for each algorithm by fitting this part of the dataset. The validation set is utilized to choose the best combination of algorithm choice and its hyperparameters (settings that are configured before training the model and are not learned from the data), and also the best features. The test set is used to evaluate how the model performs for new data. A significant portion of machine learning projects is dedicated to data preparation, underscoring its importance for the success of the model. A word of caution is necessary when using data obtained through different DFT methods, as certain properties, such as band gaps, can exhibit significant variability. This variability often stems from the underestimation of experimental values by specific functionals. Therefore, combining data from studies employing different functionals or computational setups must be approached carefully, as machine learning models rely on consistent and homogeneous datasets to ensure reliability. To address this challenge, we propose using the digital object identifier (DOI) of the publication from which the data were retrieved as an additional descriptor. Including this information allows the integration of results from computational studies with slightly different parameters while maintaining traceability and enhancing the interpretability of machine learning predictions.

  2. ii.

    Feature Engineering: Features serve as input variables used to predict output properties. In materials design, they are also referred to as descriptors and can include a large range of material properties, from experimental data to computed values obtained by methods like DFT. In addressing materials-related problems, prioritizing highly predictive features can be as crucial as the choice of algorithm, with expert knowledge playing an important role. Although a vast array of features, ranging from hundreds to thousands, may be initially considered for a given material, typically only a select few prove to be significant for the effectiveness of the final model. Consequently, the process of feature selection becomes indispensable. This process may incorporate various methodologies, such as filter methods that rely on specific statistical metrics, wrapper methods like backward feature elimination, or specialized algorithms including LASSO (Least Absolute Shrinkage and Selection Operator) regression59, which are instrumental in refining the model for optimal performance.

  3. iii.

    Model optimization: Model optimization aims to find the best combination of algorithm and hyperparameters. Most ML tasks can be categorized as unsupervised or supervised learning. Unsupervised learning works with unlabeled raw data to reveal relationships and patterns, while supervised learning uses labeled data to make predictions. Unsupervised learning makes use of clustering and association techniques, while supervised learning can be divided into regression (to predict numerical outputs) and classification (to identify data classes). Some common algorithms for these tasks are linear regression, k-nearest neighbors, decision trees, random forest, gradient boosting, support vector machines, and neural networks. A description of these algorithms can be found in ref.60. XGboost (an implementation of gradient boosting)61 is often the top choice for standard tabular numerical data, while deep neural networks excel in handling images, audio, and text, but can also be a top choice for tabular data62. Nevertheless, usually, the higher the complexity of the algorithm, the lower the explanatory power63. Despite various factors being in play, a recent empirical study suggests that the explanatory power of decision trees and random forest takes over other ML models64. Errors associated with model optimization can usually be categorized in two sources: underfitting and overfitting. Underfitting occurs when a model does not fit the training data well, while overfitting makes a model too specific to the training data, not being able to be generalized for new data. Ideally, k-fold cross-validation should be used to evaluate the model, since it helps estimate how well the model is likely to perform for new data, and can help mitigate issues such as overfitting. Both k = 5 and 10 have been verified to have low bias towards the validation set at a reasonable computational cost, and are commonly employed65. For example, in tenfold cross-validation, the model gets tested ten times on 10% validation subsets after ten training rounds on the remaining 90% of the cross-validation data.

  4. iv.

    Feature importance: When satisfactory performance is obtained from cross-validation during the model optimization phase, feature importance should be investigated in order to interpret the model and understand what the algorithm has found in the data, understand if more features related to those with higher predictive power should be added to the model or, otherwise, to eliminate non-important features that are adding noise to the model. There are several feature importance methods. Mean Decrease in Impurity (MDI) is a specific method for tree-based methods (decision trees, random forest, XGBoost, etc.) and computationally inexpensive to obtain. Feature Permutation Importance (FPI), where each column representing a feature is randomly shuffled one by one, and the model evaluated each time, with the more important features, when shuffled, resulting in a greater degradation of the performance metrics of the model. FPI can be applied to any ML algorithm, but it is also more computationally expensive than MDI, for example66. SHAP (Shapley Additive exPlanations) values is a methodology based on cooperative game theory, where these values are calculated by considering all possible combinations of features and assessing the change in the output of the model when each feature is added to these combinations. It involves retraining the model numerous times, which makes it even more computationally intensive than FPI, but can also be applied to any algorithm67.

  5. v.

    Final test: The final test set is intended to evaluate the real performance of the model, using new data, unseen by the model during the previous feature engineering and model optimization steps. Typically, it comprises 10% to 30% of the whole dataset. If further model refinement is needed, new data should be used for this step. The minimum requirements of a useful ML model can vary greatly depending on the domain of application, the nature of the data, and the problem being solved. In terms of statistical performance metrics, our experience tells us that, as a rule of thumb, an R2 of at least 0.5 should be considered for regression problems, and an accuracy of at least 70% would be expected for classification problems, but even this depends on each problem and other statistical metrics might tell a different story for better or worse. For regression problems, the R2, mean absolute error and root mean squared error should be reported68, while for classification problems, the accuracy, precision, recall, balanced accuracy or F1 score, and kappa value or Matthews correlation coefficient69, should be considered. The robustness to overfitting (drop in the quality of the metrics from cross-validation to final test), training time and computational resource efficiency, scalability with the amount of data, and interpretability, are also factors to consider.

Structure of MXenes

MXenes were first synthesized just over 10 years ago, and that is enough time for plenty of properties of these materials to have been studied extensively. This includes MXenes that have been manufactured, and ones that have not been obtained in the laboratory yet but have been predicted to be stable by computational studies70. Many computational studies have taken advantage of the experimental observation of bare MXenes to simplify the models used by neglecting surface termination. Comprehensive overviews of which MXenes have been experimentally reported, and of the known properties of MXenes in general, can be found in refs.71,72. Here, we present an overview of a selection of recent works in which DFT calculations were employed to study properties of already known MXenes or newly introduced ones.

Effect of surface termination on MXene properties

After the first MXene, Ti3C2, was isolated, computational works were quick to predict that many others, with different M or X elements or with different thickness, should be stable as well, and efforts immediately started to create them in the experimental laboratory. Indeed, for example, around five years after the advent of MXenes, the Ti2N and Ti4N3 MXenes had been synthesized73,74. The Tin+1Nn MXene variant with three Ti layers, Ti3N2, has apparently been elusive so far, but it is still a representative model of the titanium nitride MXenes. Although upon synthesis MXenes are initially covered by O, F, and OH terminal groups75,76, the oxygen group is by far the one that adsorbs the strongest, thus being the most likely to be found77. Furthermore, post-processing methods easily replace the whole surface termination by O76. This motivated Aghaei et al., to analyze the properties of the Ti3N2 MXene versus those of its oxygen-covered counterpart, Ti3N2O2, (see Fig. 2) using DFT78.

Fig. 2: Molecular models for bare and O-covered Ti3N2 surfaces.
figure 2

Side views of models of the Ti3N2 (left-hand side) and Ti3N2O2 (right-hand side) MXene surfaces. The values shown are bond lengths in Å, taken from ref.78.

Full size image

The chosen exchange-correlation functional was the generalized gradient approximation (GGA) variant introduced by Perdew, Burke, and Ernzerhof (PBE)27, supplemented by dispersion corrections by Grimme41, to account for longer-range van der Waals interactions. This combination has become known as PBE-D, which may be followed by a number depending on the version of the Grimme’s correction, and, as we will see, has become one of the standard choices in computational studies of MXenes.

Using PBE-D2, Aghaei et al. found that Ti3N2O2 displays cohesive energy that is 30% lower than that of the clean MXene, signifying greater stability of the functionalized surface. Structurally, the innermost Ti-N ionic bond lengths were found to remain almost unchanged by the addition of O, while the outer Ti-N covalent bonds are stretched by 5%. This is accompanied by an accumulation of electrons both on N and on O atoms, so that all Ti-N bonds become ionic, and the only remaining covalent bonds are Ti-O. Both Ti3N2 and Ti3N2O2 were confirmed to be metallic77, but oxidation was predicted to negatively impact the conductivity of the material. The O surface termination was also predicted to considerably reduce radiation absorption and reflectivity in the visible region, while increasing them at higher frequencies.

Bjork and Rosen studied the competition between twelve surface terminating species (O, OH, N, NH, NH2, S, SH, H, F, Cl, Br, and I) on six experimentally obtained MXenes (Ti2C, Nb2C, V2C, Mo2C, Ti3C2, and Nb4C3) using DFT-PBE79. They confirmed that the chemical elements that adsorb the strongest on both nearly bare and fully covered MXenes are O and F. This result held for both vacuum and a simulated oxygen atmosphere. Infact, it was found that, in an oxygen atmosphere, it is thermodynamically advantageous to have no surface termination than to have a sulfur one, as O2 reacts with the S termination to form SO3 and then O replaces S as the surface termination. Furthermore, O and F were predicted to be the most thermodynamically favorable species to fill any kind of termination vacancy of any of the studied MXenes. In particular, the exposition of an OH-covered MXene to O2 was predicted to remove the hydrogen from all OH groups, a conclusion that was also reached by Gouveia and Gomes80. The calculations by Bjork and Rosen also showed that a termination vacancy of an O-covered MXene is most likely to be filled by F than by any other species, indicating that an oxygen treatment does not replace the F surface termination and that a hypothetical fully F-terminated synthesized MXene should be very resistant against oxidation at ambient temperature. This oxidation resistance was observed both experimentally81 and theoretically82, although it has been demonstrated that high temperatures can remove the fluorine termination but not the oxygen one7.

In a later study83, Hou et al. focused their attention on the oxidation process of MXenes in aqueous systems. They considered the vanadium carbide MXene with oxo terminations as a prototype system and molecular models in which the empty space—often used in DFT studies to separate MXenes from their replicas in the periodically repeated images—was replaced by a varying number of water molecules. They performed DFT-based molecular dynamics (MD) simulations (PBE-D3), within the NVT ensemble (T = 300 K) and using small molecular models with up to 180 atoms (V24C12O24·40 H2O model), which were employed to develop a neural network potential for MD simulations using much larger models, with up to 20250 atoms (V5400C2700O5400 · 2250 H2O). Their simulations are supportive of a spontaneous energy minimization process in which water molecules first interact directly with V atoms before proton release occurs. At the atomic scale, an oxygen atom from a water molecule moves toward the V2CO2 MXene surface, displacing a V atom from its equilibrium position and forming a V–O bond. Supported by two neighboring water molecules, the O–H bonds of the adsorbed water break one after the other as the V–O bond shortens, ultimately leading to the formation of a [VO4]δ− structure. Importantly, this incipient stable form of vanadium oxide together with the formed H3O+ species where found to inhibit subsequent oxidation reactions in nearby regions, hence contributing to the stabilization of the MXene against oxidation and preventing its very fast degradation83.

Dahlqvist and Rosen performed PBE-D3 and rev-vdW-DF2 calculations with the VASP code to analyze how the nature and degree of surface termination of MXenes can affect, and hence used to control, their structural, electronic, and mechanical properties84. They considered halogen- and chalcogen-terminated MXenes, specifically M2CTx (M = Ti, Zr, V, Nb, Ta; T = S, Se, Te, Cl, Br, I; x = 1, 2). The findings show that MXenes with larger terminal atoms (such as Br and Cl on Ti2C, Nb2C, and Ta2C) favor non-ideal termination coverage (x < 2), resulting in mixed termination sites and lower binding energies (i.e., energy difference ongoing from multilayer to single sheet MXene models), with concomitant propensity to delamination into single layers, though very low coverage can lead to structural collapse. Their study also revealed that non-ideal coverage increases the density of electronic states at the Fermi level, potentially boosting conductivity.

Atomic layer stacking and ordering

A peculiar way to tune the properties of MXenes is to change the way their atomic layers are stacked. Sun et al. first introduced this concept for the Mo2C, Tc2C, and Os2C MXenes85. The underlying idea is that, for some MXene compositions, the surface is more stable if all the metallic atomic layers are aligned with each other, in an ABA configuration with D3h symmetry, rather than staggered with respect to each other, as in the traditional MXene ABC configuration with D3d symmetry (see Fig. 3).

Fig. 3: Molecular models for the W2N MXene surface with ABC and ABA stacking.
figure 3

Top and side views of models of the W2N MXene showing the two known stable atomic layer stacking possibilities for MXenes: ABC (left-hand side) and ABA (right-hand side).

Full size image

Gouveia et al. then conducted a systematic DFT study which covered nine M element possibilities, X being C or N, the M2X, M3X2, and M4X3 stoichiometries, both atomic layer stacking possibilities, and bare or O-terminated surfaces, for a total of 216 MXenes, in order to determine what makes different MXenes prefer different stackings86. They used the PBE functional with the Grimme dispersion corrections in their PBE-D3 version40, because this study served as a reference to additional studies to further analyze how the change in atomic layer stacking can affect the adsorptive and catalytic properties of the MXene that are reviewed in separate subsections below. All MXenes made of M in group 6 of the Periodic Table and the MXene nitrides with M in group 5, apart from V2N, were found to thermodynamically prefer the ABA stacking. Regardless of the preferred configuration, the energetic difference between the ABC and ABA configurations was predicted to change linearly with the thickness of the material, indicating that the relative stacking stability is per MXene layer. Furthermore, the surface O coverage bolsters the preference for the ABA stacking, making it even more stable when it is the most favorable one, or less unstable when ABC is preferred. The energy barriers for the transition from ABC to ABA stacking were calculated using the dimer method87, and led to the conclusion that this transition is kinetically feasible, especially on O-terminated MXene surfaces. Nearly all energy barriers were found to be below 1 eV, the lowest ones being 0.09 and 0.10 eV, corresponding to Cr2CO2 and W2NO2, respectively. The shift to the thermodynamic stabilization of ABA stacking was ultimately attributed to a higher electronic density of the material, since stronger preference for the ABA alignment was found for M elements with more d electrons, for MXenes made of N rather than C, and in the presence of O surface groups.

Based on this study, Jurado et al. devised a strategy to experimentally identify the stacking of the Mo2C MXene by using CO, CO2, or H2O as probe molecules88. A difference of almost 1 eV was calculated for the adsorption energy of CO2 on ABC-Mo2C or ABA-Mo2C. Moreover, the frequency of the asymmetric stretching vibrational mode of CO2 and of the stretching mode of CO were found to differ by 150 and 300 cm−1, respectively, which should be more than enough to clearly distinguish the two stacking possibilities experimentally. This provides a spectroscopic method to assess whether an ABC-to-ABA transition (or vice versa) has taken place in practical applications.

Triggered by the discovery of a new sub-family of M5X4Tx-like MXenes, more specifically, of the Ti2.5Ta2.5C4Tx and Ti2.675Nb2.325C4Tx MXenes89, Downes et al. performed non-spin-polarized PBE calculations to predict stable compositions of MAX phases with potential for to add new MXenes to the M5X4 sub-family. More specifically, these authors studied (M′1–yM″y)5AlC4 MAX phases, either with the traditional P63/mmc or the twinned P(bar{6})m2 structures, with M′ and M″ elements randomly arranged, keeping out-of-plane chemical ordering but with in-plane disorder. Interestingly, while in most cases the systems displaying solid-solution disorder were more stable than when having ordered arrangements, a few cases, namely, Mo3Sc2AlC4 (P63/mmc), Mo3Ti2AlC4 (P(bar{6})m2) and Mo3Sc2AlC4 (P(bar{6})m2) have a stable ordered Mo-M″-Mo-M″-Mo (M″ = Sc or Ti) stacking as concluded from the calculated Gibbs free energy of formation at T = 2000 K89.

Several research groups have focused on the distinct structural and compositional characteristics of high-entropy (HE) 2D carbide MXenes, suggested as an exotic family of these promising materials. In particular, Nemani et al. synthesized two high-entropy MXenes, TiVNbMoC3 and TiVCrMoC3, highlighting their stability and complex elemental arrangements, as well complemented with results from DFT calculations with the PBE functional90. The latter allowed to calculate the enthalpy of formation of different MAX phase compositions, from the total energies of these structures and of the constituting M, A, and X elements, per atom, in their standard state (bulk phases), as well as the Bader charges of the M elements. Two important observations in their work were that i) the configurational entropy—the number of discrete representative positions of the alloy constituents, calculated as −0.1773 and −0.2238 eV/f.u. (formula unit) for three-M and four-M systems, respectively—is an important contributor to the stabilization of the phases, and ii) the synthesizability of MAX candidates is higher in those candidates with an electron concentration close to 0.3 e/atom.

Building on the results from the two studies reviewed above—namely, that certain combinations of M’ and M” (e.g., Mo and Ti) lead to ordered MXene phases89, and that HE MXenes can be experimentally synthesized90, Leong et al. conducted a high-throughput first-principles study of the atomic configurations in the TiVNbMoC3 and TiVCrMoC3 MXenes, focusing on possible metal segregation91. To explore the vast compositional space, their approach involved constructing surrogate models for the MXene configurational energies using a generalized Ising model. These models were trained and parameterized with PBE, and subsequently employed to perform accurate calculations of configurational energies for various atomic arrangements in HE MXenes, across different compositions and temperatures, using Monte Carlo simulations. Their calculations predicted a distinct preferential occupancy of metal atoms across the four layers of single-phase HE MXenes, even at temperatures as high as 2900 K. Over a broad compositional range, the outer/inner M layers are predominantly occupied by two of the four metals, with a residual amount of a third element in the M’M”M”’M””C3 MXene, with Cr most preferentially occupying the outer layers, followed by Mo, V, Nb, and Ti. A clear separation of Nb and V atoms was found to occur in the inner layers and, within each M layer, the atoms largely form a solid solution due to the increased configurational entropy.

A novel HE MXene with composition Ti1.1V1.2Cr0.8Nb1.0Mo0.9C4Tx was successfully synthesized by Ma et al.92, who performed in tandem a computational study with the PBE functional and considering different configurations of Ti1V1Cr1Nb1Mo1C4Tx HE MXene model with and without (random positioning) ordered M layers. As reported by other authors89, they found a propensity for structures with M layers arranged in a solid solution to display lower formation energy than structures with ordered M layers. Also, the calculated density of states of the 1Ti: 1V: 1Cr: 1Nb: 1 Mo composition is higher at the Fermi level than those of other Nb2C, V2C, and Ti3C2 MXenes92, which was interpreted as an indication that the former MXene has better conductivity and electron migration than the latter.

Magnetic ordering

MXenes offer significant promise for advancing 2D magnetism research due to the wide range of potential transition-metal compositions. However, experimental investigations into their magnetic properties are still limited. Studies have primarily focused on titanium-based carbide MXenes93,94,95. While titanium itself is non-magnetic, it can contribute to magnetic behavior when alloyed with other elements (e.g., Cr in ref.93) or structurally modified (e.g., by applying various surface terminations through acid treatments as in ref.94), where effects like magnetic transitions and Pauli paramagnetism have been observed.

García-Romeral et al. considered the PBE, PBE0, and HSE06 functionals to carefully analyze the electronic ground state of the Ti₂C MXene96, namely, to investigate the relative stability of antiferromagnetic (AFM), ferromagnetic (FM) and non-magnetic (NM) states. An AFM ground state for Ti₂C was predicted with the three functionals but the PBE0 functional was found to stabilize FM states more than the PBE or HSE06 approaches. The DFT-derived spin exchange parameters (J1, J2, J3) revealed that the intralayer FM coupling (J2) is the dominant interaction, but the other interactions were significant and opposite in sign. The results also suggested a spin model with one unpaired electron per Ti atom, highlighting a strong ionic bonding component, contrasting with prior assumptions of two unpaired electrons per Ti in a triplet state97. Hence, García-Romeral et al.’s analysis establishes an AFM ground state for Ti2C with robust FM intralayer coupling, which can be considered promising for 2D magnetic applications. An important technical aspect arising from their work is that while spin polarization minimally impacts structural parameters, it notably influences total energy stability among different magnetic configurations96. PBE calculations on p(1 × 1) unit cells of Ti2C, Zr2C, and Hf2C confirmed that these MXenes exhibit magnetic ground states, with ordering of stability AFM > FM > NM solutions98.

The topic had already attracted the attention of several computational researchers who considered DFT to predict the magnetic behavior of MXenes incorporating metals like Cr (Cr2C99,100, and Cr2N100, with H, O, OH, F or Cl terminations), or Mn (Mn2C101, with F, Cl, OH, O, or H terminations). In the latter, He et al. performed calculations with the DFT+U approach (U = 3 eV) and with the HSE06 functional, and found that the Mn2C MXene retained a ferromagnetic (FM) ground state even upon symmetrical functionalization with singly negative charged groups (e.g., F, Cl, OH). Furthermore, the same authors found that surface functionalization enables tunable magnetic states, producing FM half metals, antiferromagnetic (AFM) metals, or AFM semiconductors depending on the electronegativity of the functional groups. The change in the properties upon functionalization were previously reported by Si et al. for the Cr2C MXene99, but in the latter case upon surface functionalization with the same −1 charged groups, there was an FM-to-AFM transition only. Hence, the Mn₂CT₂ MXenes show unique electronic and magnetic properties that are promising for spintronic applications, notably, the Mn₂CF₂ MXene which exhibits intrinsic half-metallicity with a Curie temperature above room temperature (TC = 520 K)101.

More recently, García-Romeral et al. conducted further computational studies to investigate how increasing the number of metal layers in titanium carbide MXene—from 2 to 3 to 4—affects its magnetic properties102. They considered the PBE, PBE0, and HSE06 approaches, already used in their previous work96, as well as the PBE+U (U = 3 eV) one. All methods indicated that the Ti₂C, Ti3C₂ and Ti4C3 display magnetic ground states and that the magnetic coupling constants increase monotonically with the thickness of the MXene. These authors proposed a spin model where the surface titanium atoms are paramagnetic (Ti+ with one unpaired electron) while those in the middle layers are diamagnetic (Ti2+ closed-shell). These authors extended their calculations to investigate the electronic ground state of M2C (M = Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, and W) MXenes103, revealing that Ti2C, Cr2C, Zr2C, and Hf2C MXenes exhibit open-shell conducting ground states, with spin polarization localized on the metal atoms. In contrast, V2C, Nb2C, Mo2C, Ta2C, and W2C MXenes display diamagnetic conducting ground states. Further analysis of trends across the Periodic Table groups indicated that the energy difference between FM and NM states decreases when descending along Group IV, whereas it increases when ascending along Groups V and VI.

Vacancy defects

Materials in the laboratory are never as crystalline as models depict them, and MXenes are no exception. In fact, atomic vacancies were quickly identified on the Ti3C2 MXene, the first to be synthesized104,105. The understanding of how likely these defects are to be formed and how they change the properties of the MXene has always been accompanied by theoretical DFT-based studies, such as refs.106,107 for the Ti2X family of MXenes. More recently, higher-throughput studies have been conducted to find the most likely types of vacancies on many MXenes108.

Parey et al. studied three kinds of vacancies on 18 different M2X-type MXenes, with M = Sc, Ti, V, Cr, Mn, Y, Zr, Nb or Mo and X = C or N. The DFT calculations were made using the PBE-D3 approach (the D3 correction was included to later analyze the influence of vacancies on CO2 activation)109. The three types of vacancies considered were metal vacancies (VM), X-element vacancies (VC or VN), and metal-X vacancies (VMC or VMN). Thermodynamically, carbide MXenes were predicted to be more likely to contain vacancies than nitride ones and, within each of these two groups, the most stable vacancy is VC or VN, respectively.

The results of Parey et al. are supported by those of Gouveia et al., who studied the structural and energetic properties of vacancies in 84 MXenes, each composed of one among nine metallic elements, C or N, pristine or with O or F surface terminating groups, and with ABC or ABA stacking108. To calculate the formation energy of the different vacancies, the energy references were as follows: for metal atoms, the energy per atom of one unit cell of hexagonal-close-packed (hcp) Ti, Zr, and Hf metals or body-centered-cubic (bcc) V, Nb, Ta, Cr, Mo, and W metals; for carbon, the energy per atom of one unit cell of diamond; for nitrogen and oxygen, half the energy of gaseous nitrogen (N2) and oxygen (O2), respectively; and for fluorine, the energy difference between gaseous hydrogen fluoride (HF) and half of gaseous hydrogen (H2). Gouveia et al. found that in O-terminated MXenes, C or N vacancies are still the most likely to be found, while M atoms become much harder to remove. For example, on Ti2C, the VTi formation energy was calculated at just over 2 eV, while this value increased to more than 8 eV on Ti2CO2, as per previous studies106. This is a consequence of both the well-known oxophilicity of MXenes and stability of O-terminated MXenes. Indeed, for MXenes with F surface groups, it was predicted that the relatively lower VM formation energies are restored, VTi lying at around 3.5 eV. According to the calculations, the F groups should be much easier to remove than O groups, with F vacancy formation energies being lower than O ones by several eV. Here, the authors raise a crucial point on the calculation of defect formation energies: the resulting value varies depending on the reference used for the energies of single atoms. As an example, F vacancy formation energies calculated using ({rm{HF}}-frac{1}{2}{{rm{H}}}_{2}) or (frac{1}{2}{{rm{F}}}_{2}) as reference for the energy of one F atom yield formation energy differences of circa 3 eV. These references should be considered when comparing the formation energies of atoms of different chemical elements obtained in different studies.

Gouveia et al. also found MXenes that appeared to display unrealistically negative vacancy formation energies of a few tens of eV, along with severe surface deformation, leading to a structure that could hardly be identified as an MXene. This was attributed to the instability of the ABC atomic layer stacking on these materials, motivating the analysis of the ABA-stacked counterparts, which yielded much more moderate vacancy formation energies and maintained the integrity of the MXene lattice108. The unstable MXenes in question include the already synthesized Mo2CF2 and Mo2NF2. Thus, the work of Gouveia et al. showed that, even if these MXenes are produced through exfoliation of an ABC-stacked parent MAX phase, the resulting material will display ABA atomic layer stacking, as observed experimentally81,82.

Li et al. compared the stability of C vacancies in Ti2CO2 with that of clusters of C vacancies, namely C vacancy lines110. Four types of vacancies were studied, as shown in Fig. 4: single C vacancies, nearest-neighbor C vacancies, 3-C vacancies along a symmetry plane of the slab, and a full line of C vacancies spanning the whole material. Of these four, the one predicted to be thermodynamically less favorable is the double C vacancy, with a formation energy per C atom of 7.76 eV. On the opposite end of the stability ranking is the full-line vacancy, with a formation energy of 7.63 eV. Furthermore, the calculated binding energies, i.e., the energy of the systems with respect to all its isolated atoms, indicated that the defective Ti2CO2 slabs were thermodynamically more stable than the pristine MXene, and therefore might be found experimentally. The calculations also hinted that the electronic and magnetic behavior of the Ti2CO2 MXene changes upon introduction of carbon vacancy lines. The pristine material is a semiconductor and displays no magnetism, whereas the introduction of either of the four types of studied defect induces metallic behavior, and magnetism in the case of the one- and three-C vacancies. The absence of a full line of C atoms (Fig. 4) led to the highest conductivity due to the presence of a large amount of dangling bonds around the vacancies. All the defective systems appeared to display higher radiation absorption than pristine Ti2CO2 in the infrared region, while the opposite was found for the visible region, and their potential application as anode materials in aqueous or ionic/organic systems was suggested.

Fig. 4: Molecular model for the Ti2CO2 MXene surface.
figure 4

Top view of the Ti2CO2 MXene, assumed infinite along its plane, with values of formation energies of carbon vacancy clusters, per carbon atom, in eV, calculated in ref.110. The silver-, gray- and red-colored spheres denote atoms of titanium, carbon, and oxygen, respectively. The blue areas denote the studied types of C vacancy clusters: a full line of carbon atoms spanning the whole material (long horizontal rectangle), three carbon atoms contained in a plane of symmetry of the material (slanted blue area), two adjacent carbon atoms (bottommost blue area), or a single carbon vacancy (blue circle).

Full size image

The calculated electronic properties of defective MXenes have been found to depend on the exchange-correlation functional used108,111. For instance, both the PBE and the HSE06 functionals predict the Ti2CO2 MXene to be a semiconductor but, with any kind of single-atom vacancy, the PBE predicts a metallic behavior, while HSE06 yields a small band gap of a few hundredths of an eV. As discussed in ref.111, discrepancies between semi-local and hybrid functionals arise from two main sources. First, lingering strain in pre-relaxed structures can lead to errors in single-point hybrid-DFT calculations, particularly for extended defects, where strain accumulates across bulk-like regions. Second, the underestimation of band gaps by semi-local functionals can cause spurious defect level mixing with host states, leading to artificial structural changes. These factors suggest that fully relaxing defective MXene structures with hybrid-DFT approaches is crucial to accurately capture their electronic and structural properties. For Ti2CO2, the metallic behavior predicted by PBE may reflect such artifacts, while the HSE06 results may provide a more reliable description of defect-induced properties.

The works reviewed above clearly indicate that the band gaps of MXenes can be influenced by the presence of vacancies, which, when introduced in a controlled manner, can lead to tailor-made materials for specific applications. This is exemplified by the work of Tao et al.112, who successfully synthesized a molybdenum carbide MXene with a significant proportion of metal vacancies, designated as Mo1.33C, from a chemically in-plane ordered quaternary MAX phase, specifically (Mo2/3Sc1/3)2AlC. Through detailed experimental and computational analyses (using the PBE approach), a close qualitative and quantitative agreement was observed, allowing the conclusion that the structure of Mo1.33C closely resembles that of Mo₂C but with ordered divacancies at every third Mo site. Compared to Mo₂C, the new Mo1.33C MXene exhibits a reduced mass density, along with significantly higher conductivity and volumetric capacitance112.

Interaction with atoms of transition metal elements

The properties of MXenes can be modified not only by changing their overall structure, like mixing different X or M elements or functionalizing the whole surface, as we have seen, but also by slightly doping the surface with extra atoms (i.e., off-lattice doping). This is a widely recognized method in catalysis, since doping the surface of materials with atoms of another element can induce favorable changes in the electronic structure of the materials that often improve the overall catalytic activity113. More recently, the consideration of single atoms of a metallic element deposited onto a substrate, forming a so-called single-atom catalyst (SAC), to improve reaction efficiency and minimize metallic waste114. With this in mind, theoretical115 and experimental116,117 studies have examined the interaction of MXenes with single transition metal atoms and their effect on catalysis. More recently, systematic first-principles studies, which consider many MXenes and many transition metal elements, have led to the establishment of patterns in the way MXenes interact with atoms of these elements.

In 2021, Oschinski et al. employed a PBE-D3 approach to investigate the interaction of an atom of each of ten different transition metal elements (fourth period of the Periodic Table) with each of nine MXenes with M2C stoichiometry (M = Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W)118. The authors of that work discovered that all MXenes behave qualitatively similarly. For example, on all nine MXenes, as one moves along the fourth period of the Periodic Table, the transition metal atom adsorption strength increases up to vanadium, then decreases abruptly for chromium, gradually increases again up to iron, cobalt or nickel, and dips again to its lowest value for zinc. The adsorption was found to be spontaneous and very exothermic, by 2 to 6 eV, for every MXene/metal combination. The Bader charges on the metal adatoms also follow the same pattern for all MXenes: Sc is the most positively charged atom, while Co or Ni are the most negatively charged ones. Intriguingly, despite the very strong transition metal atom adsorption on the MXenes, the potential energy profile on the surface is almost flat in many cases, and their diffusion energy barriers along the surface are rather low, ranging from 0.02 eV for Sc diffusion on Ta2C, to 0.69 eV for Ti on Cr2C. Apart from Ti diffusion on Cr2C, all diffusion energy barriers are lower than 0.36 eV.

Later, Keyhanian et al. conducted a similar study, using the same approach, on nine oxygen-terminated MXenes, with stoichiometry M2CO2 (M = Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W)119. Again, all surfaces examined presented the same qualitative behavior towards transition metal atom adsorption. The preferred adsorption site is the one that allows the adatom to form the most bonds with the surface: three bonds with O and three bonds with M surface atoms. In the presence of the O layer, the adsorption is strongest for Sc (ranging from approximately 4 eV on Hf2CO2 and 9.5 eV on Cr2CO2), becomes weaker along the first row of transition metals, reaching its lowest for Zn (between approximately 0 eV on Hf2CO2 and 2 eV on Cr2CO2). The magnitude of the adsorption energies calculated is consistent with those found for adsorption of metal atoms on metal oxide surfaces120. All metal atoms become positively charged when adsorbed on M2CO2 MXenes, with the strongest adsorption generally being associated to more charge transferred. The calculated diffusion energy barriers are much higher than the ones obtained by Oschinski et al. for the bare MXenes118, especially for Sc, Ti, and V diffusion. Not surprisingly, the diffusion of the barely adsorbed Zn involves virtually no energy barrier.

The trends in MXene-metal interactions were generalized by Gouveia et al. to transition metal elements of the second and third rows121. Taking advantage of the fact that all MXenes display the same qualitative behavior towards transition metal adsorption, Gouveia et al. considered only the Ti2C and Ti2CO2 MXenes, since these are the most studied ones, and thirty different metal elements. On Ti2C, metals of the second and third rows adsorb on the same surface sites as the corresponding ones of the first row: a bridge site for groups 3–5 of the Periodic Table, a hollow site above a C surface atom for groups 6–10, and a hollow site above a Ti atom for groups 11 and 12. This suggests that the most favorable adsorption site is related to the filling of the d subshell of the transition metal atoms. On Ti2CO2, most atoms adsorb in a position that allows them to form three bonds with Ti and three bonds with O, in agreement with the findings of Keyhanian et al.119.

In terms of adsorption strength on Ti2C, as shown in Fig. 5, the behavior of the atoms of the three transition metal rows predicted by Gouveia et al. was found to be quite different but, on each row, the adsorption is still strongest for the Fe, Co, and Ni group, and weakest for atoms of the Zn group121. The latter is also true for Ti2CO2, but on this MXene the adsorption is strongest for the first metals of each period, and simply becomes weaker along each period. Figure 5 also shows the amount of charge transferred between the MXene surfaces and each transition metal atom. One can observe that, for all three periods, stronger adsorption is roughly accompanied by more charge transferred. This is especially noticeable on Ti2CO2, where all Bader charges are positive, indicating a transfer of electrons from the metal atom to the MXene. This is justified by the presence of oxygen, the second most electronegative element of the Periodic Table.

Fig. 5: Adsorption energies and charges for metals on titanium carbide surfaces.
figure 5

Adsorption energies (Eads, top panels) and Bader charge (ΔQ, bottom panels) of transition metal atoms adsorbed on Ti2C (left-hand side) and Ti2CO2 (right-hand side) MXene surfaces. Values taken from ref.121.

Full size image

The diffusion energy barriers of transition metal atoms on Ti2C were predicted to be lower for atoms at the leftmost and rightmost ends of each period of the Periodic Table and were calculated at 0.25 eV (for Ir) at most121. All studies agree that, on Ti2CO2, these diffusion barriers are much higher, but decrease along each period, thus following the same trend as the adsorption strength. The high barriers predicted for transition metal atom diffusion on Ti2CO2 were attributed to the difficult process of the metal adatom breaking three metal-Ti in order to travel over the O layer and bond with three distinct Ti atoms in a neighboring unit cell (see Fig. 6). This difficulty does not occur in the bare Ti2C, on which diffusion only requires breaking one metal-Ti bond at a time while keeping the other two intact. The stability of the MXene/metal systems demonstrated in all the studies presented in this section paved the way for the potential development of MXene-supported single-atom catalysts with unique properties.

Fig. 6: Diffusion mechanism of metals on titanium carbide surfaces.
figure 6

Illustrations of the diffusion mechanism of a transition metal atom adsorbed on Ti2C (left-hand side) or Ti2CO2 (right-hand side), starting from a hollow site above a carbon atom, and ending in a neighboring equivalent position, following the path represented by the arrows. Titanium atoms are represented in silver color, carbon in gray, oxygen in red, and the metal atom in blue. The semi-transparent blue circles are intermediate metastable states of the diffusion.

Full size image

Prediction of synthesizable MAX phases and MXenes

Although tens of MAX phases and MXenes have been realized, hundreds have been predicted to be synthesizable. In 2021, Khaledialidusti et al. used a DFT-based computational method, with the PBE functional, to estimate the stability of MAX phases, each composed of one of 11 M elements, one of 17 A elements, C or N, and (n=1), 2 or 3, for 1122 MAX phase candidates. Khaledialidusti To this end, the phonon dispersion and formation energy of the MAX structures were calculated. The first of these quantities is related to the second-order derivative of the energy of the system with respect to the atomic positions and was employed to assess the dynamic stability of the materials. In turn, the formation energy is a measure of how stable the configuration of the system is when compared to other configurations. In practice, to be considered locally and thermodynamically stable, a MAX phase should display the absence of modes with imaginary wavenumbers and negative formation energy.

The formation energy criterium was found to be much more restrictive than the phonon dispersion, in the sense that all the M/A/X/n combinations that did not meet the phonon dispersion criterium, did not meet the formation energy one either. This suggests that future studies could in principle restrict themselves to calculating formation energies to arrive at reliable conclusions. In short, MAX phases were predicted to be more stable when i) the M element belongs to group 4 or 5 of the Periodic Table; ii) the A element is not Si, P, S, Cu, or Ir; iii) the X element is C. A total of 466 MAX phases were predicted to be synthesizable, where the number of carbides is more than twice that of nitrides. This is around 3 times as many as the number of MAX phases that have actually been realized.

Encouragingly, from the calculated force constants between the M and A elements and between the M and X elements, Khaledialidusti et al.70 observed that the M-A bonds are much weaker, as expected, since the way MXenes are traditionally produced is precisely by breaking the M-A bonds in MAX phases2. The thermodynamic stability of the resulting MXenes, relative to the precursor MAX phase, was assessed by calculating the exfoliation energy. Weak M-A forces and thermodynamic stability of MXenes with respect to the precursor MAX phases were used as criteria to establish 136 MAX phases as the ones that can be exfoliated to produce 26 MXenes, shown in Table 1. Note that 12 of these 26 had already been manufactured at the time of the prediction, but others came to be synthesized later, such as Ta2C122, Cr2C123, V2N124, and Ti4C3125.

Table 1 The 26 synthesizable MXenes, and corresponding precursor MAX phase A elements, predicted by the DFT calculations in ref. 70
Full size table

DFT calculations have been employed to propose alternative ways of synthesizing MXenes, such as epitaxial growth. By taking advantage of the excellent capabilities of MXenes for dissociating nitrogen126, Gouveia et al. suggested obtaining thicker MXenes experimentally by depositing alternating layers of nitrogen and a metallic element over a thinner precursor MXene98. For this to be realizable, a series of steps must be thermodynamically and kinetically favorable, such as exothermic metallic atom adsorption. The limiting step was determined to be the adsorption and dissociation of N2 on the nearly fully N-covered MXenes, predicted to be achievable for nearly half of the 18 considered MXene seeds. These results are useful as guidelines for experimentalists, as they provide the most likely obtainable MXenes, as well as the most probable paths for obtaining them.

Prediction of new MXenes and their properties

Despite not having been synthesized to date, the Hf2C MXene has attracted the attention of plenty of theoretical studies, possibly because Ti3C2 was the first to be synthesized, Zr3C2 was realized a few years later127, and therefore MXenes based on the next transition metal of this group of the Periodic Table, Hf, are expected to be stable.

Liu et al. compared the electronic properties of Hf2CT2, with six different surface terminations (T = O, S, Se, F, Cl, OH), with those of the bare Hf2C MXene, using the PBE-D3 method128. The geometries of the six terminated materials were predicted to be qualitatively the same, with the functional groups adsorbed on hollow sites over Hf atoms. The lattice constant of Hf2C was predicted to increase with surface termination, from 3.20 Å on pristine Hf2C, up to 3.45 Å on Hf2CSe2. Results for Hf2CSe2 were later confirmed by the calculations of Jasani et al., using the PBE exchange-correlation functional without dispersion corrections, and some elastic properties of this theoretically predicted material were calculated129. By comparing the binding energies, any Hf2CT2 surface was predicted to be more stable than the clean one, with emphasis on Hf2CO2, which displays the lowest binding energy by far. This was an expected result, as MXenes are known to be very oxophilic. The Hf2CT2 MXene was found to be metallic both with and without any surface termination, apart from Hf2CO2, which has a semiconductor indirect band gap of 1.0 eV. Functional groups other than O shift the Fermi level into the conduction band, inducing a semiconductor-to-metal transition. The Hf2CCl2 and Hf2CF2 systems have work functions (the difference between the vacuum energy and the Fermi energy) considerably smaller (3.78 and 3.59 eV) than those of the other studied functionalized surfaces (5.36 to 8.84 eV), and therefore were proposed as applicable in thermionic and electron emission devices. The functionalization of Hf2C was found to eliminate the magnetism inherent to the clean surface and change the type of electrode. In particular, the Cl-, OH-, and F-terminated Hf2C MXenes are potential anode materials, while O-, S-, and Se-terminated ones are potential cathode materials.

The authors of ref.128. additionally considered the effect of mixed O, OH, and F surface terminations, which are known to occur experimentally. Three different ratios were studied: while keeping the OH content at 0.11, the O/F stoichiometric coefficients considered were 0.28/1.61, 0.78/1.11, and 1.61/0.28. Because of the presence of F and OH, all three mixed-termination MXenes displayed metallic behavior, implying that, in reality, Hf2CT2 is most likely a metal, even if most of its surface termination is composed of O adatoms. The lattice constant, density of states and work function of the Ti3C2T2 and V2CT2 MXenes with mixed O, OH, and F terminations were predicted to be very well modeled by a weighted average of the same properties of their uniformly terminated counterparts130.

In a subsequent study, Liu et al. investigated the effect of substituting C atoms by N ones (i.e., in-lattice doping) in the Hf2CO2 MXene131. This was again done using the DFT-D3 method and a 3 × 3 supercell of the material, so that ten different N/C ratios (0/9, 1/8, 2/7, …, 9/0) were considered. In terms of structural stability, the lattice constant was found to decrease as the N doping amount increases, so that when all C’s are replaced by N, the lattice constant is 1.5% smaller than on pristine Hf2CO2. However, in terms of binding energy, the calculations showed that the stability of the system is reduced as N doping increases, ultimately making Hf2NO2 less stable than Hf2CO2 by around 0.2 eV per atom. The semiconductor behavior of Hf2CO2 was found to disappear even with the smallest finite N/C ratio, due to the Fermi level rising above the conduction band minimum. Calculations of the surface storage charge suggested that Hf2CO2 is a potential cathode material, while the slightest N doping turns the MXene into a potential anode material. Interestingly, as the N doping amount increases, the surface tends to become a symmetric electrode, implying that tuning the N concentration should yield any desired type of electrode. Once again, a mixture of surface terminations was considered as well, Hf2CO0.42F1.19(OH)0.39, for comparison with the purely O-terminated MXene. The semiconductor character of the Hf2CO2 MXene was concluded to be exclusive to this pristine material, since the replacement of even a minor amount of C by N, or of the O surface termination by other species or a mixture of O, F, and OH turns the MXene into a metal.

In 2022, 10 years after the isolation of the first MXene, Michałowski et al. synthesized samples of the Ti3AlC2, Mo2TiAlC2, and Cr2TiAlC2 MAX phases and the Ti3C2Tx MXene using traditional methods, and analyzed them using secondary ion mass spectroscopy (SIMS), thus studying their structure at the atomic level132. They discovered that there was oxygen present amid the carbon layers of both classes of materials. The oxygen content was considerable even on the MXene, amounting to 30% replacement of the carbon atoms with oxygen ones, yielding a Ti3C1.4O0.6 stoichiometry. This led to the suggestion that, most likely, nearly all carbide MXenes produced to date are in fact oxycarbides. Taking inspiration from this finding, Gouveia and Gomes studied the structural and energetic stability of substitutional oxygen on the Ti3C2 MXene for several oxygen concentrations, using the PBE-D3 approach133. The calculations showed that, at least up to 75% replacement of carbon content with oxygen one, the substitution is exothermic and the resulting Ti3C2−xOx MXene structure is indistinguishable from that of the pristine Ti3C2 MXene, explaining the difficulty and delay in experimentally identifying the oxycarbide. The Ti3C2-xOx MXenes were also predicted to display metallic behavior independently of the oxygen concentration on the carbon layer. If all carbon is replaced by oxygen, the resulting Ti3O2 2D material is no longer stable with its original ABC atomic layer stacking and should transition to an ABA phase, which is more stable than ABC by 0.37 eV per formula unit. As a consequence, Ti3O2 became the first material derived from a Ti-based MXene to prefer the ABA layer alignment, confirming the postulate by Gouveia et al., that the ABA phase of MXenes is stabilized and driven by a higher electronic density of the material86.

In a more recent study, Michałowski et al. developed a protocol for deconvoluting SIMS data, enabling the quantitative analysis of the composition of individual MXene atomic layers134. This advancement allowed for the distinct detection of not only oxycarbide MXenes but also oxynitride and oxycarbonitride MXenes. Notably, analysis of multiple MXene samples revealed that some could achieve high oxygen concentrations (>65%) in the X layer, corroborating earlier theoretical predictions by Gouveia and Gomes133.

Simultaneously, Gouveia and Gomes135 examined the effects of subsurface substitutional oxygen atoms on the structural, energetic, and electronic properties of MXenes with varying M, X, and T elements. They studied MXenes of the form M2XT2, with M = Ti or V; X = C or N; and T = none, O, or F, which, when compared with their earlier work on the Ti3C2 MXene133, also allowed them to analyze the influence of MXene thickness. Their findings demonstrate that subsurface O is thermodynamically stable in any concentration across all analyzed MXenes. On bare Ti2C, subsurface O atoms migrate to the surface, forming carbon vacancies and localized O terminations. The results in their work135, align with experimental evidence suggesting that all synthesized MXenes may, in fact, be oxycarbides or oxynitrides134. Furthermore, they also suggested that conditions that remove surface terminations are likely to induce significant X-layer vacancies. Notably, O-terminated nitride MXenes (M2NO2, M = Ti or V) exhibit spontaneous distortions, whereas pristine O-terminated carbide MXenes (M2CO2, M = Ti or V) retain hexagonal symmetry unless subsurface O exceeds ~75 %. Fluorine terminations stabilize hexagonal structures, even in nitride MXenes, making F-based synthesis routes particularly suitable for these materials. Electronic analyses revealed that pristine MXenes are metallic, except for the semiconducting Ti2CO2, whose band gap is influenced by subsurface O. While PBE predicted a loss of semiconducting behavior at low O concentrations, HSE06 suggested that the material remains semiconducting, albeit with increased conductivity as the O content rises135.

Even though most MXene research focuses on those made of a single M element, multimetallic MXenes have been synthesized. The possibility of controlling the multimetallic composition of MXenes paves the way for even broader applicability of this family of materials. One such MXene is TiVC, which was first synthesized in 2022136. Sometime later, a theoretical study by Zhang et al. examined the structural properties of TiVCTx considering different surface terminations (F, O, and OH), using the PBE functional137. Different relative positions of the Ti and V atoms were considered, and the most asymmetric one, shown in Fig. 7, was deemed the most stable. As is usual with ABC-stacked MXenes, the F, O, or OH terminal groups prefer being placed on surface hollow sites above a metallic atom. The frequencies of the Raman- and infrared-active vibrational modes of the system were calculated for comparison with the ones obtained experimentally. According to the calculations, the functionalization of the TiVC surface introduces new active modes and considerably changes the vibrational frequencies of the existing ones, as shown in Table 2. With this information and the experimental data from ref.136, it was concluded that the surface termination of the as-synthesized TiVC MXene was composed of F and OH groups, as demonstrated by the very small deviation between the frequencies calculated for TiVCF(OH) and the experimentally measured ones.

Fig. 7: Molecular model for the TiVC bimetallic MXene surface.
figure 7

Top and side views of the optimal configuration of the TiVC bimetallic MXene surface, as obtained in ref.136. The spheres represent atoms of Ti (blue), V (green), or C (gray).

Full size image
Table 2 Frequencies of the local vibrational modes, in cm−1, of the TiVCTx MXene with different surface terminations, calculated in ref.137 and measured in ref.136
Full size table

MXenes as catalysts

Every chemical reaction follows a different mechanism, involving relocation of atomic nuclei and redistribution of electronic charge, to break some chemical bonds and form new ones. More often than not, somewhere along the reaction path, there are intermediate states of the system that are less stable than the reactants in their initial state, which is more likely to happen if the initial bonds are stronger. The least stable of these states is called the reaction transition state. This means that, for the products of the reaction to be formed, the system needs to overcome an energy barrier, called the reaction activation energy, which can be defined as the difference between the energy of the system in the reaction transition state and in its initial state. An analogy with an everyday phenomenon could be the toppling of an object: even though the object is more stable when toppled than when sitting upright, one needs to tilt the object by a certain minimum amount, in order to cause gravity to spontaneously pull it towards the toppled position. Likewise, chemical reactions can be facilitated by the presence of a third party, a catalyst. The catalytic process starts with the catalyst forming chemical bonds with the reactants, weakening the existing bonds within the reactants. Secondly, an external energy source provides enough energy to break the weakened bonds. Then, new bonds are formed, giving rise to the products of the reaction. Finally, the catalyst is regenerated by removing the reaction products, and can be used again. Thus, a suitable catalyst for a given reaction should be able to bind to the reactants, should effectively lower the reaction activation energy, should preferably be selective towards the desired reaction products and not competing side reactions that may occur, and should not be too strongly bound to the products, allowing their separation after the reaction has taken place. Different catalytic categories use different energy sources to overcome the reaction activation energy: thermal catalysis uses temperature, photocatalysis uses electromagnetic radiation, and electrocatalysis uses an external voltage.

Thermal catalysis

Since MXenes are solids, they can serve as catalysts in heterogeneous catalysis, i.e., catalysis in which the phase of the catalyst is different from the phase of the reactants or products. Here, the MXene is the (solid) catalyst and the reaction of interest involves gases. From the theoretical point of view, in particular at the DFT level, assessing the suitability of an MXene for use in catalysis consists of calculating the reactant adsorption energies on the surface, calculating the activation energy for the reaction to occur (which should be lower than in the gas phase), and calculating the adsorption energy of the products.

Mathematically, in a catalyzed reaction, the initial and final states of the system correspond to energy minima with respect to atomic positions, but the transition state is a saddle point of the potential energy surface, and these are much harder to find. Two common methods are the dimer and the nudged elastic band (NEB) methods. Both approaches attempt to find saddle-point system configurations. To use the dimer method, one can take the initial state of the system and manually change the coordinates of the atoms of interest for the reaction, moving them in the direction that one expects them to go for the reaction to take place. For example, to study the catalysis of H2 dissociation on a solid surface, one would first find the ground state adsorption configuration of H2 on the surface, then change the atomic coordinates of the H atoms so that they are slightly farther from each other, and use the resulting configuration as input for the dimer calculation. The NEB method is computationally more expensive and requires the knowledge of the final state of the system, but is more likely to find a saddle point than dimer. Instead of using only one configuration as input, NEB uses several configurations (images), resulting from linear interpolation of the atomic positions between the initial and final states. The method then optimizes all the images simultaneously while keeping equal spacing between consecutive ones. Herein, we present some examples of recent works that reveal the sort of information that DFT calculations can provide in this field.

Dehydrogenation of propane

The dehydrogenation of propane catalyzed by the Ti2CO2-x(OH)x MXene, i.e., the Ti2C MXene with a surface termination consisting of some O groups and some OH groups, was studied by Liu et al. using the van der Waals density functional (vdW-DF), an exchange-correlation functional augmented with van der Waals corrections138,139. The reaction transition states were found by a combination of the NEB and dimer methods: first, the initial and final states of the system (hydrogenated and dehydrogenated propane adsorbed on the MXene) were found, then a NEB calculation with 10 intermediate images provided the dehydrogenation reaction path, and finally the central images were used as dimer inputs for obtaining more accurate transition states. On the purely O-terminated surface, Ti2CO2, the dehydrogenation of the CH2 group of propane was predicted to be easier than that of a CH3 group, since the former requires surpassing an energy barrier of 1.60 eV while the latter has a higher one, 2.01 eV140. The authors then built 24 different Ti2CO2-x(OH)x MXene models, with varying OH/O fraction and distribution, and calculated the propane activation energy on each, shown in Fig. 8. As seen in this Figure, different arrangements of the same fraction of OH and O groups can lead to very different activation energy values, implying that the OH fraction is not a good descriptor. However, the activation energies were found to correlate linearly with the corresponding reaction energies, as in a Brønsted–Evans–Polanyi (BEP, also shown in Fig. 8) relationship, and with a more useful physical quantity, the surface hydrogen affinity. This can be roughly defined as the adsorption energy of a hydrogen atom, and was determined to also strongly depend on the surface termination group distribution, since the highest and lowest calculated hydrogen affinities were found for the same fraction of OH groups (43.75%). Unsurprisingly, the Ti2CO2-x(OH)x MXene models that were predicted to more strongly adsorb hydrogen atoms were the ones with the lowest fraction of OH groups, and yielded the lowest activation energies for propane dehydrogenation on average. The linear relation obtained allows the estimation of dehydrogenation activation energies as a function of the hydrogen affinity, with a mean absolute error lower than 0.2 eV, both at the CH3 and CH2 group. Unlike reaction energies, the hydrogen affinity is an intrinsic property of the material and can therefore be measured experimentally and used to predict the suitability of an MXene for alkane dehydrogenation.

Fig. 8: Activation and reaction energies for propane dehydrogenation.
figure 8

(Top panel) Activation energy for propane dehydrogenation via a CH3 (blue points) or the CH2 (orange points) group, catalyzed by the Ti2CO2-x(OH)x MXene as a function of the percentage of the O surface termination that is replaced by OH groups. (Bottom panels) Brønsted–Evans–Polanyi (BEP) correlations between the reaction energies energies (ΔE) and activation energies (Ea) for propane dehydrogenation via a CH3 (left-hand side) or the CH2 (right-hand side) group. Data from ref.140.

Full size image

Dehydrogenation of methylcyclohexane

Bare MXenes, i.e., MXenes without surface functionalization, have their outer metallic layers exposed and are known to be very reactive, spontaneously dissociating many molecules80,141,142. Other species, e.g., CO2, adsorb on the bare surface without dissociating spontaneously as predicted computationally143, despite the low calculated dissociation energy barriers, and observed experimentally144. This can be useful when the molecule dissociates in the desired way, but it becomes harmful if the molecule is completely disintegrated. The latter was predicted to be an expected outcome if one attempts to dehydrogenate methylcyclohexane to toluene using a bare Ti3C2 MXene surface145. According to the DFT-PBE calculations of Obodo et al., attempting to catalyze this reaction using Ti3C2 would result in a too exothermic process that destroys the molecule and deactivates the surface. However, after depositing layers of Pt over the MXene surface, which has been shown to be experimentally feasible146, the calculations showed that the reactivity of the material becomes much more moderate. By simulating the presence of one, two, or three Pt layers over Ti3C2, Obodo et al. noticed that the reaction becomes endothermic, but the endothermicity is reduced as the number of Pt layers increases. With as few as three Pt layers, the catalytic activity is already improved with respect to the Pt(111) surface.

Dissociation of nitrogen

Armed with the knowledge of the high reactivity of bare MXenes towards most molecules, Gouveia et al. built models of 18 MXenes and performed DFT-based calculations, using the PBE-D3 exchange-correlation functional, to study N2 dissociation on their surfaces126. This is an important reaction, involved in the Haber-Bosch process for industrial production of ammonia. In fact, N2 dissociation is the rate-limiting step of ammonia synthesis, because it involves breaking a triple N-N bond, which in gas phase would require around 10 eV of energy. N2 was found to exothermically adsorb and dissociate on all the studied surfaces. The calculated activation energies are shown in Table 3. Notice that even the highest one represents a decrease of nearly 90 % with respect to gas phase dissociation. Some MXenes were therefore predicted to display N2 dissociation energy barriers lower than those found for some of the established catalysts for this reaction, such as Fe(111) or Ru nanoparticles, with energy barriers of circa 0.80 eV and 0.4–1.9 eV, respectively147,148. The lowest value is achieved by the W2N MXene, which has not been experimentally realized yet. The exothermic reaction energies correlate reasonably well with the activation energies, an occurrence known as Brønsted−Evans−Polanyi relationship149,150,151. In order to form NH3 using MXenes, the successive hydrogenations of the adsorbed N adatoms become the rate-limiting steps, which are endothermic and involve energy barriers between 1.24 and 1.70 eV, suggesting that MXenes would have to be used as co-catalysts with, for example, Ru nanoparticles, in the Haber–Bosch process. An alternative pathway, such as the Mars-van Krevelen mechanism152,153, might be feasible for producing ammonia since the successive dissociation of N2 molecules on the MXene surfaces has been predicted to create a nitrogen adlayer98.

Table 3 Energy barriers for N2 dissociation catalyzed by bare MXenes taken from Gouveia et al.126
Full size table

In their study of N2 dissociation on bare MXenes126, Gouveia et al. observed that MXenes composed of an M element of group 6 of the Periodic Table became distorted in the vicinity of the adsorbed molecule, in agreement with previous predictions by Shao et al.154. This led them to investigate whether this happened due to an attempt by the surface to transition to another more stable phase, which they found to consist of an ABA atomic layer stacking, instead of the usual ABC one found on most MXenes86. This stacking is more stable than the traditional one for 8 out of the 18 MXenes that were analyzed and this was found to have consequences on their catalytic activity. For example, according to the calculations, the ABA-stacked W2N MXene presented an N2 dissociation energy barrier, of 0.18 eV, which is even lower than the 0.28 eV found for ABC stacking. The former signifies a reduction of over 96% with respect to the gaseous phase barrier.

Water-gas shift reaction

The dissociation of water molecules is often the rate-limiting step in processes such as the water-gas shift reaction or the hydrogen evolution reaction155,156,157,158. The former is a very important process involved in the industrial production of hydrogen. The DFT calculations of Gouveia et al. showed that 18 studied MXenes greatly facilitate the water dissociation step by reducing the activation energy of the first O-H bond breaking to almost zero159. In particular, MXenes made of an M element from group 4 of the Periodic Table (Ti, Zr, or Hf) display water dissociation energy barriers below 0.1 eV. The energy barriers and reaction energies follow a BEP relationship, as in the case of N2 dissociation by MXenes126. However, it was predicted that water dissociation would be too exothermic, for example when compared with metallic surfaces160,161,162, which would be a major problem for regenerating the catalyst by desorbing the products of the reaction.

The hope of using MXenes to catalyze the water-gas shift reaction was restored by Deeva et al., who synthesized the molybdenum carbide MXene and used it to catalyze this reaction81. These authors prepared three different samples, namely, Mo2CTx, where F has the highest fraction of terminating groups, and two partially reduced Mo2CTx-500 and Mo2CTx-800 obtained from the former sample, using 20 vol% H2 in N2 at 500 and 800 oC, respectively. Under water gas shift conditions, the oxophilic nature of partially reduced samples leads them to be covered with oxygen, which was suggested to reduce the catalytic activity81. With these results in mind, Gouveia et al. conducted another DFT-based study, this time including Mo2C models that had surface functionalization82. The dissociation of a water molecule was studied on five different models, mimicking five different regions of the MXene: a model with no surface termination, a model with Mo2C fully covered by O (Mo2CO2), a model with full F termination (Mo2CF2), and two models in which the F termination is partially replaced by one or two O adatoms (Mo2CO0.13F1.63 or Mo2CO0.25F1.38, respectively). The models were inspired by the results of Deeva et al., which stated that the surface of the synthesized material was functionalized by mostly (70%) fluorine adatoms (i.e., Mo2CF1.4Tx). The quantitative results are summarized in Table 4. On the bare Mo2C, water was predicted to adsorb exothermically and require a low energy barrier to dissociate, but the products (OH and H) remain too tightly bound to the surface. On Mo2CO2 and Mo2CF2, water weakly physisorbs relatively far from the surface, with adsorption energies of −0.20 and −0.13 eV, respectively. On Mo2CO2, this is not enough to activate an O-H bond, since the energy required to break it on this surface is 2.50 eV. In contrast, the Mo2CF2 actually reduces the activation energy to 0.31 eV, although this is higher than the absolute value of the adsorption energy on this surface, implying that water is more likely to desorb again instead of dissociating. The issue with water dissociation on Mo2CF2 is that, when the H atom binds to a surface F atom, an HF molecule is released, and the surface termination is locally replaced by OH. These qualitative findings for water dissociation on Mo2C, Mo2CO2, and Mo2CF2 motivated the development of the mixed-termination models, since the following rules became apparent: i) the MXene cannot be bare, since that leads to irreversible surface poisoning; ii) the surface termination cannot be comprised of too much oxygen, since regions fully functionalized by oxygen are inert towards adsorption and dissociation; and iii) the surface termination cannot be made solely of F adatoms, since this termination is removed in the form of HF upon water dissociation, damaging the surface. According to the calculations of Gouveia et al., the local replacement of F adatoms by one or two O groups creates a section of the Mo2C MXene that is reactive enough to exothermically adsorb water and dissociate it, with energy barriers of 0.14 or 0.20 eV, respectively, without being damaged in the process or leading to reaction products that are too stable. The products are hydroxyl groups that remain adsorbed on the surface, but not on their thermodynamically most favorable adsorption site, which may increase their reactivity for the following steps of processes that begin with the dissociation of water. This study showed that the synthesis of MXenes by HF etching of the A element of a MAX phase leads to a surface that is mostly covered by F, which is partially replaced by O, and that the resulting surface composition is fortuitously adequate for catalyzing the dissociation of water molecules.

Table 4 Adsorption energies (Eads), dissociation energy barriers (Ebarr), and dissociation reaction energies (Ereac), in eV, of a water molecule on the Mo2C MXene, with different surface terminations, as calculated by Gouveia et al.82
Full size table

The previous paragraphs illustrate just a few examples of how easy it is for bare MXenes to dissociate adsorbed molecules, and therefore how hard it usually is to form new molecules on these surfaces. We now proceed to describe this contrast on MXenes, namely for CO2 dissociation and CO oxidation.

Carbon dioxide reduction and carbon monoxide oxidation

Morales-Salvador et al. used the PBE-D3 exchange-correlation functional to show that 18 MXenes very easily dissociate CO2 molecules into CO+O, with Mo2N, Hf2N, and V2C eliminating the energy barrier entirely143. The CO2 adsorption and dissociation were studied not only on the bare MXenes but also on Ti2CO2 and Mo2NO2 with a small bare region where three O adatoms have been removed. Even in this small area, the CO2 dissociation energy barriers were predicted to be as low as 0.14 and 0.01 eV, and the desorption energy of the resulting CO became 0.29 and 1.75 eV, respectively. The results showed that, after dissociating CO2, the MXenes ultimately become covered by an oxygen terminating layer, which is known to be removable by hydrogenation into H2O, regenerating the catalyst. Parey et al. extended these results to MXenes with vacancy defects, again using a DFT-based PDE-D3 methodology109. Three types of vacancies were considered: M vacancies, C/N vacancies and MC/MN vacancies. The latter was predicted to require the most energy to be formed, in agreement with the calculations of Gouveia et al.108, whereas C/N vacancies are thermodynamically the easiest to form. The CO2 molecule was found to adsorb more strongly on most defective MXenes than on the pure ones (i.e., CO2 is attracted to the vacancies), most notably on the C/N-defective V2C, Nb2C, Ti2N, V2N, Zr2N, and Mo2N MXenes, on which CO2 dissociates spontaneously. Since the C/N vacancies are the most likely to occur, this means that the list of MXenes that spontaneously dissociate CO2 in the laboratory should be even larger than the one predicted by Morales-Salvador et al.143

The reverse process, i.e., CO oxidation, is unfeasible on bare MXenes, given the almost non-existing energy barrier and clear exothermicity of CO2 dissociation on these surfaces. Two recent theoretical works used a DFT-D3 approach to predict the feasibility of CO oxidation on transition metal single-atom catalysts based on the surface-functionalized Ti2C MXene. Zhu et al. analyzed the potential of single Fe atoms deposited on the Ti2CO2 MXene surface (Fe/Ti2CO2) for CO and O2 adsorption, and then for CO oxidation163. The CO molecule was found to exothermically adsorb on the Fe/Ti2CO2, with an adsorption energy of −1.16 eV, clearly differing from the undoped Ti2CO2, which is known to barely adsorb even the smallest molecules80,142,164,165. Curiously, Fe/Ti2CO2 can even adsorb two more CO molecules, with subsequent adsorption energies of −1.33 and −0.82 eV. The adsorption of O2 is stronger than that of CO by 0.42 eV, implying that, statistically, O2 should adsorb first. After one O2 adsorbs on Fe/Ti2CO2, the adsorption energies of CO and O2 become −0.57 and −0.23 eV, respectively, thermodynamically leading preferably to CO+O2 co-adsorption. A second CO adsorption, i.e., 2CO+O2 co-adsorption, was predicted to be possible as well, and to yield the energetically cheapest mechanism for CO2 formation:

$$2{{rm{CO}}}^{* }+{{rm{O}}}_{2}^{* }to {{rm{CO}}}_{2}^{* }+{{rm{CO}}}^{* }+{{rm{O}}}^{* }to 2{rm{C}}{{rm{O}}}_{2}^{* }$$
(1)

The first and second CO2 formations require energy barriers of only 0.20 eV and 0.15 eV, respectively. The CO2 desorption energy, of 0.03 eV, is also very low. Peng et al. considered a SAC based on the same MXene but with a different surface functionalization, Ti2CN2, with a single atom of Zn166. Here, the calculated CO and O2 adsorption energies were −0.89 and −0.88 eV, respectively. The CO oxidation mechanism predicted to be the easiest to realize is shown in Eq. 1, and consisted of gaseous CO interacting directly with the adsorbed O2, forming gaseous CO2 and leaving behind an O adatom, which then bonds with a second CO to form a second CO2. The two reactions of this mechanism have very low energy barriers, of 0.06 and 0.03 eV, respectively, and are exothermic by over −2 eV. The desorption energies of the first and second formed CO2 are 0.23 and 0.44 eV, respectively, indicating that the formed CO2 molecules could be desorbed at room temperature. Given the popularity and vast knowledge about the Ti2C MXene and the non-nobility of Fe and Zn, these results unveil the Fe/Ti2CO2 and Zn/Ti2CN2 single-atom catalysts as promising cheap and very active catalysts for CO oxidation.

Fischer-Tropsch synthesis

Kountoupi et al.167 studied by experimental and computational approaches the suitability of the Mo2CTx MXene of the Fischer–Tropsch (FT) process168,169,170 under steady-state conditions (p = 25 bar and T = 600 K), focusing on how the activity and selectivity are affected by the surface oxygen and carbidic carbon contents. From the experimental side, they concluded that the defunctionalization of the MXene using H2 treatment at 400 oC yields a Mo2C MXene-based catalyst with 55% selectivity towards higher hydrocarbons (alkanes with five or more carbon atoms). Furthermore, when the H2 pre-treatment is made at 500 oC, carbon vacancies are created and the resulting catalyst can only form methane, not larger hydrocarbons. According to their DFT calculations, direct CO dissociation into C* and O* is easy avoiding the formation of other oxygenate species beyond CO2 which becomes the major reaction product. The hydrogenation of CH3* species to methane was predicted to involve a very high energy barrier thus leading to the formation of larger alkanes. Both the experimental and computational results led to the conclusion that the surface will become poisoned by O-adatoms due to CO decomposition reaction, with or without H2 assistance, and the high-oxophilicity of the MXene.

Photocatalysis

In photocatalysis, electromagnetic radiation is used to promote electrons from the valence band to the conduction band, allowing their reaction with oxidants to create reduced products, and generate holes in the valence band, which can react with reductants to create oxidized products. In general, for this to successfully happen, two things are needed: the radiation energy must be greater than or equal to the band gap of the catalyst, and the lifetime of the electron–hole pair (called an exciton) must be long enough to keep the pair from recombining. Typically, infrared radiation (photons with energies between 1.24 meV and 1.7 eV) is used for this process because its energy range coincides with the range of energies separating the quantum levels of molecular vibrations and it contains the band gaps of many semiconductors. Infrared also represents more than half of the energy received from sunlight at the surface of the Earth. Studies on photocatalysis using MXenes are relatively scarce, probably because most MXenes are either metallic or have band gaps that are too wide54,171. According to DFT calculations using the very accurate HSE06 exchange-correlation functional172, and machine-learning predictions54, the MXenes most likely to have suitable band gaps for photocatalysis are made of M = Sc, Ti, or Y. We highlight Sc2C, whose Janus-type surface termination was extensively studied in a PBE study by Zhang et al. in search for a photocatalysis-capable Sc2C-based MXene173. Ten surface termination combinations were chosen from T, T’ = O, F, H, OH, or Cl, forming Sc2CTT’ Janus surfaces. Of these, Sc2CFCl, Sc2COHCl, and Sc2CHCl were deemed the most promising for photocatalysis because of their calculated band gaps, 1.79, 0.99, and 1.84, respectively, and large internal electric field. The Sc2COHCl MXene should be suitable for photocatalysis using infrared radiation.

Ontiveros et al. explored MXenes as photocatalysts for driving the water splitting reaction57. Using the PBE functional for optimization runs and the PBE0 hybrid functional for more reliable band gap estimations, they analyzed the band alignment of M2X and M2XO2 MXenes, where M = Sc, Y, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, or W, and X = C or N, with respect to water splitting potentials. Bare MXenes (i.e., without O-termination) were found to be metallic, whereas oxygen termination induced semiconducting behavior, with group III and IV MXenes emerging as the most promising candidates. Carbide MXenes (X = C) exhibited larger and more frequent band gaps when compared to nitride MXenes (X = N). Notably, for carbide MXenes, the band gap decreased for group III metals and increased for group IV metals when moving down the group, whereas an inverse trend was observed for nitride MXenes. Ontiveros et al.57 also identified MXenes with band gaps exceeding 1.23 eV, the minimum photon energy thermodynamically required to drive water splitting, as determined by the oxidation potential of H2O to O2 (1.23 V vs. NHE). Among the O-terminated MXenes, Zr₂CO₂ emerged as a standout candidate due to its favorable band alignment with water-splitting potentials and significant charge density separation between the valence and conduction bands. Interestingly, the effect of stacking on the calculated band gaps was found to be subtle but sufficient to influence the electronic properties under certain conditions.

Electrocatalysis

Climate change and the increasing global energy demands are among the main problems that humanity faces today. The main electrochemical processes involved in the field of renewable energies are hydrogen evolution (HER), oxygen evolution (OER), oxygen reduction (ORR), and carbon dioxide reduction (CO2RR)174. These processes occur at the solid/liquid interface between an electrode and an electrolyte solution. From the theoretical point of view, the mechanism via which an electrochemical processes occurs over an electrocatalyst can be described by a sequence of reaction intermediates and the Gibbs free energies changes between them ((Delta G))175. This information can be condensed in a free energy diagram, which shows how the free energy changes along the reaction coordinate176.

Regardless of the preferred mechanism of an electrochemical reaction over an electrocatalyst, the overall (Delta G) at zero electrode potential is the same. For example, OER is a four-electron-transfer process,

$$2{{rm{H}}}_{2}{rm{O}}to {{rm{O}}}_{2}+4{{rm{H}}}^{+}+4{{rm{e}}}^{-}$$
(2)

Assuming an equilibrium potential of zero for the reversible hydrogen electrode, the equilibrium potential, ({U}^{0}), of OER reaction is 1.23 V. Since the whole process involves transferring four electrons, this means that, if one sums the Gibbs free energy changes along the reaction path of OER over any electrocatalyst, the value thus found is always (4times 1.23,=,4.92{rm{eV}}). In other words, in the free energy diagram of OER at zero applied potential, the final state lies 4.92 eV above the initial state. The application of an electrode potential can increase or decrease (depending on its polarity) the free energy of the intermediate and final states with respect to the initial state, potentially driving the reaction in the desired direction. Taking again the example of OER, an electrode potential of ({U}^{0}) = 1.23 V brings the reaction to an equilibrium, by definition bringing the free energy of the products, from 4.92 eV above the reactants, to 0 eV, (Gleft({rm{reactants}}right)=Gleft({rm{products}}right)). In this situation, the system is at electrochemical equilibrium, and the net current is zero. If one applies an electrode potential, (U), that is greater than the one needed to achieve equilibrium, then the forward direction will be favored, and O2 will be produced. The difference, (eta =U-{U}^{0}), is the definition of the overpotential. Note that the electrode potential can stabilize not only the reaction products but also all the intermediates whose formation involves at least one electron or proton transfer. Obviously, one can almost always apply a very large overpotential and force all electrochemical steps to become exergonic, but this is energetically unprofitable. Realistically, the applied overpotential should not exceed 0.3 V.

The activity of an electrocatalyst towards a certain electrochemical reaction can be assessed from the knowledge of the sequence of (Delta G) connecting the reactants and products, and by using the Sabatier principle and assuming a Brønsted–Evans–Polanyi relationship177. These two results state that the initial, intermediate and final species should not adsorb too weakly or too strongly on the electrocatalyst, and that the energy barrier values depend linearly upon the reaction energies. The descriptor ({G}_{max }left(eta right)), defined as the largest free-energy span between all reaction intermediates, was introduced by Exner to assess the performance of electrocatalysts178. This descriptor accounts for the applied overpotential and kinetic effects, and includes the concept that in catalysis, there are no rate-determining steps, but instead rate-determining states, as concluded by Kozuch and Shaik179.

Hydrogen evolution

Hydrogen evolution is a reaction that produces H2 from two protons. Although it is suitable for producing gaseous H2, it is very commonly a nuisance competing side reaction in processes like nitrogen fixation or CO2 reduction and, for this reason, knowledge on how likely it is to occur on an electrocatalyst is always useful. In general, two possible mechanisms are considered. The first reaction step is common to both and is called the Volmer step, in which a proton is reduced on the electrocatalytic surface,

$${{rm{H}}}^{+}+{{rm{e}}}^{-}+* to {{rm{H}}}^{* }$$
(3)

In the Volmer-Tafel (VT) mechanism, the second (Tafel) step assumes that two Volmer steps have occurred, so that two ({{rm{H}}}^{* }) exist on the surface, which can then undergo a chemical step, generating molecular H2,

$$2{{rm{H}}}^{* }to {{rm{H}}}_{2}+2*$$
(4)

In the Volmer-Heyrovsky (VH) mechanism, the second proton is reduced directly on top of the first one, directly forming H2,

$${{rm{H}}}^{+}+{{rm{H}}}^{* }+{{rm{e}}}^{-}to {{rm{H}}}_{2}+*$$
(5)

The activity of the Ti3C2Tx MXene as a catalyst for the HER was studied by Meng et al., considering a large number of surface termination possibilities180. Surface termination phase diagrams were built, depicting the most stable surface functionalization for each (U, pH) pair, the so-called Pourbaix diagrams. The equilibrium potential of the HER is, by the usual definition, 0 V. At the equilibrium potential, and assuming the whole surface is functionalized by the same species, the most stable surface termination was found to be -O. In more realistic models, allowing binary and ternary functionalization, competition was found between Tx = O2/3(OH)1/3, F1/3O1/3(OH)1/3 and F3/9O4/9(OH)2/9, with similar stability, with the F-containing terminations occurring only at low pH. Gibbs free energy diagrams were obtained for all the surface termination possibilities. The strong H affinity of the pristine Ti3C2 MXene was found to prevent the formation of H2. Interestingly, the three most stable terminations under working conditions were found to be the most promising for HER performance, since they require almost negligible overpotentials (0.01 V) for all the electrochemical steps to become exergonic (Table 5). For the O2/3(OH)1/3 and F3/9O4/9(OH)2/9 terminations, the preferred mechanism is the Volmer-Heyrovsky one, while on Ti3C2F1/3O1/3(OH)1/3, the Volmer-Tafel one is favored. In comparison, the V2C MXene was predicted to achieve maximum electrocatalytic activity for HER when with a full H surface termination181. These findings show that the most favored surface termination for HER is MXene-dependent and that these surfaces can compete with other electrocatalysts commonly used for the HER, such as Pt, MoS2, or WS2.

Table 5 Calculated data for the hydrogen evolution reaction on the Ti3C2Tx MXene surface with different surface terminations
Full size table

A combination of experimental and computational approaches was used to analyze the role of surfaces terminations (Tx) on the HER activity of five different carbide MXenes (viz. Ti2C, Ti3C2, Mo2C, Mo2TiC2, and Mo2Ti2C3)182. Handoko et al. observed that the HER activities of Ti3C2Tx samples produced by different fluorine-containing etchants were largely dependent on the fluorine functionalization of the basal plane, with the samples containing higher coverages exhibiting lower activities (larger overpotentials). Moreover, the large Tafel slopes suggested that H adsorption (Volmer step) could be rate limiting. DFT calculations with the BEEF-vdW functional were performed to understand this behavior from the changes in the free energy for H adsorption (ΔGH) with varying amounts of fluoro groups on the basal plane, as the ideal HER catalyst must hold ΔGH values close to zero. The models Ti3C2O2, Ti3C2O1.5F0.5, Ti3C2OF, Ti3C2O0.5F1.5, and Ti3C2F2 with F:Ti ratios corresponding to 0.00, 0.17, 0.33, 0.50 and 0.67 per unit cell, respectively, were used to mimic MXenes with increasing degrees of fluoro functionalization. The calculated ΔGH values were found to increase monotonically from F:Ti ratio of 0.00 to 0.50, with a more significant increase being observed on going from F:Ti ratio of 0.50 to 0.67, because the H atom adsorbs preferably on top of the oxo terminations which are not present in the basal plane of Ti3C2F2. These authors extended the calculations to carbide MXenes containing molybdenum (Mo2CTx, Mo2TiC2Tx, and Mo2Ti2C3Tx) and observed that in the Mo2TiC2Tx and Mo2Ti2C3Tx MXenes, the increase of the fluoro functionalization induces a significant increase of the ΔGH from F:Mo ratio of 0 to 0.25 and that the ΔGH stays constant for larger functionalizations. Interestingly, in the case of the Mo2TiC2F2 and Mo2Ti2C3F2 MXenes, the adsorption of H was found to be unstable and led to the formation of HF, as it was found by Gouveia and Gomes82 when exploring the feasibility of the water gas shift reaction on the Mo2CF2 MXene. In the case of the Mo2CTx MXene without F functionalization, the BEEF-vdW ΔGH value is 0.048 eV, which is closer to thermoneutrality than the ΔGH calculated for the 2H-phase of MoS2 (cf. 0.08 eV with the RPBE functional183), hence suggesting that the Mo2CTx without fluoro groups on the basal plane may be an excellent HER catalyst.The HER was also analyzed over 24 ordered double transition-metal carbide MXenes by Jin et al.184. As in the work of Handoko et al.182, the BEEF-vdW ΔGH values were used to identify the promising HER catalysts. Under the standard conditions, the most favorable surface terminations found were either the pure O functionalization, or a mix of O and OH groups, depending on the elemental composition of the outermost metallic layers. In particular, if these layers are made up of Mo, Ti or Nb, the surface tends to prefer the simple O termination, whereas when the outer layer contains Cr or V, the mixed O/OH termination is preferred. This mixed surface termination has also been predicted by López et al. for monometallic MXenes185. The activity and chemistry of the bimetallic MXenes was found to be dominated by the composition of the outer metal layer, and strongly dependent on the strength of the outer M-O bond: the weaker this bond in the MXene, the stronger the bond between the MXene O and the adsorbed H atom. The ordered double transition-metal MXene that showed the lowest overpotential (0.003 V) was Mo2NbC2O2.

MXenes have also been mixed with other materials in order to improve their durability for the HER. For example, Lim and collaborators synthesized a 2H-MoS2/Mo2CTx MXene nanohybrid186. Both their experiments and DFT calculations revealed that the sulfidation of the MXene, with an S/O ratio of 1/3, effectively prevents its oxidation and consequent deactivation, allowing the system to sustain current densities over 450 mA/cm2. This not only experimentally demonstrates the capability of this particular nanohybrid, but also suggests that the interfaces between other materials and chalcogenides can be suitable for electrochemical applications.

Nitrogen and nitrate reduction

The electrochemical nitrogen reduction reaction (NRR) to ammonia,

$${{rm{N}}}_{2}left({rm{g}}right)+6left({{rm{H}}}^{+}+{{rm{e}}}^{-}right)to 2{rm{N}}{{rm{H}}}_{3}left({rm{g}}right)$$
(6)

with an equilibrium potential of 0.057 eV, is a considerably more complex process than hydrogen evolution, involving more steps and more possible reaction mechanisms. For this reason, only bare MXenes had been studied theoretically for this application154,187,188,189, until Johnson et al. investigated the role of the surface termination in controlling the NRR electrocatalytic activity of MXenes190. It was shown that the surface termination is especially important for this reaction because the calculated Pourbaix diagrams indicated that, under the NRR operating conditions, the O and H terminations are most favorable, with free energies of several eV below the bare materials. The associative NRR mechanism comprises six electron exchanges,

$${{rm{N}}}_{2}left({rm{g}}right)+{{rm{H}}}^{+}+{{rm{e}}}^{-}to ,*, {{rm{N}}}_{2}{rm{H}}$$
(7)
$$*, {{rm{N}}}_{2}{rm{H}}+{{rm{H}}}^{+}+{{rm{e}}}^{-}to *, {{rm{N}}}_{2}{{rm{H}}}_{2}$$
(8)
$$*, {{rm{N}}}_{2}{{rm{H}}}_{2}+{{rm{H}}}^{+}+{{rm{e}}}^{-}to *, {rm{N}}+{{rm{NH}}}_{3}left({rm{g}}right)$$
(9)
$$*, {rm{N}}+{{rm{H}}}^{+}+{{rm{e}}}^{-}to *, {rm{NH}}$$
(10)
$$*, {rm{NH}}+{{rm{H}}}^{+}+{{rm{e}}}^{-}to ,*, {rm{N}}{{rm{H}}}_{2}$$
(11)
$$*, {rm{N}}{{rm{H}}}_{2}+{{rm{H}}}^{+}+{{rm{e}}}^{-}to ,{{rm{NH}}}_{3}left({rm{g}}right)$$
(12)

In other words, one N atom of the N2 molecule is three-fold hydrogenated and released as NH3, and then the second atom undergoes the same process. Due to the well-known high reactivity of bare MXenes towards nitrogen adsorption126 and the near-zero equilibrium potential, the free energy profiles of NRR on these surfaces are significantly downhill in the first parts of the reaction and equally uphill in later steps, implying catalyst poisoning by the early NRR intermediates. All bare MXenes were predicted to spontaneously form N2H, and four of them (Mo2C, Mo2N, W2C, and W2N) offered a reasonable theoretical overpotential for NRR, which is not enough to make them suitable for nitrogen fixation, as bare MXenes are not stable under NRR conditions. In contrast, H- and O-terminated MXene surfaces are more stable but not as reactive. Of the 65 MXenes compared by Johnson et al.190, only W2CH2 satisfies the three required criteria: i) low theoretical overpotential; ii) stability under NRR conditions; and iii) selectivity towards NRR instead of the unwanted HER.

Since the theoretical prediction of the stability of ordered bimetallic MXenes191 and of the optimal conditions for their synthesis192, many works have dedicated themselves to comparing the behavior of these materials with that of the traditional monometallic MXenes. This is the case of a study by Zhao et al., in which the activity of 18 double-metal MXenes, with stoichiometries ({{rm{M}}}_{2}^{{prime} }{{rm{M}}}^{{prime} {prime} }{{rm{C}}}_{2}) and ({{rm{M}}}_{2}^{{prime} }{{rm{M}}}_{2}^{{prime} {prime} }{{rm{C}}}_{3}), for electrochemical NRR was analyzed193. One important result was that the calculated (Delta G) of the potential-determining step are all below 0.98 eV, which is the benchmark value on the Ru(0001) surface. For most MXenes, this step was the second hydrogenation of the second N atom, Eq. 11. This includes Mo2Nb2C3, the MXene that was estimated to be the most promising one for several reasons. Firstly, the calculated (Delta G) potential-determining step was the lowest one (0.39 eV), with an overpotential of 0.48 V. Secondly, by performing calculations including implicit solvation, the authors found that the solvent has a positive effect on NRR activity. Lastly, the O-functionalized Mo2Nb2C3O2 MXene surface with an O vacancy, a more realistic model compared to the experimentally prepared material, displayed the same activity as the bare one. The calculations revealed that the Mo atoms that directly bind to the reaction intermediates regulate the electron transfers between these species and are probably responsible for the high activity.

Very recently, the nitrate reduction reaction (NO3RR) catalyzed by MXenes has been suggested as a very promising approach for converting nitrate into e.g., ammonia194,195, which begins with the reduction of nitrate to N adatoms and then followed by reactions described by Eqs. 10–12. On transition metal (TM) clusters with 1 to 3 atoms deposited on O-vacancies of M3C2O2 MXenes (TM = V, Cr, Mn, Fe, Co, Ni, Cu, Zr, and Hf; M = Ti, V, Cr, Nb and Mo), Zhao et al.194 performed DFT calculations and found that the deposition of the TM clusters regulates the adsorption strength of the species involved in the NO3RR. The Nb2@V3C2O2 and Nb3@V3C2O2 showed remarkable thermodynamic stability and the best catalytic activity and selectivity towards ammonia.

Gao et al.195 studied the NO3RR catalyzed by single atoms of a transition metal element from Ti to Au replacing a metal of the Ti3C2O2 with a Ti-O vacancy. The most promising systems were predicted to be the ones containing an Ag (representative precious metal) or a Cu (representative of a non-precious metal) atom. High efficiency and selectivity towards ammonia and low limiting potentials were calculated for these two cases.

Oxygen evolution and reduction

OER and its reverse, ORR are two textbook reactions that are commonly tackled in theoretical studies of electrocatalysis. While bare MXenes are extremely oxophilic and can spontaneously dissociate oxygen molecules upon adsorption, the opposite happens on functionalized MXenes, and this can deter the application of any MXene, with or without surface termination, as an electrocatalyst for oxygen evolution or reduction. Table 6 compares the adsorption energy and charge transfers of oxygen molecules on some bare and surface-terminated MXenes.

Table 6 Adsorption energies (({E}_{{rm{ads}}}) in eV) and charge transfers (Q in e) of O2 on some bare (x = 0) and functionalized (x = 2) MXenes, according to the calculations in ref.197
Full size table

Both bare and terminated MXenes are unsuited for OER/ORR, for opposite reasons. The former interact too strongly with reaction intermediates, keeping the products from desorbing, whereas the latter interact too weakly to effectively bind the species and provide electrons for the reaction to occur. One could argue that Nb2C(OH)2 is an exception to this, given the moderate O2 adsorption energy of −0.92 eV. However, this adsorption was observed to be accompanied by the spontaneous extraction of two H atoms from surface OH groups to generate an H2O2 molecule, demonstrating the instability the OH-functionalized MXene, as had been predicted by Gouveia et al.80 and Gao et al.196.

Kan et al. showed that, by doping the surface of the Nb2CT2 MXene with single atoms of Pt or Pd, the MXene-O2 interaction becomes more manageable and both OER and ORR can be conducted using the same surface as a bifunctional electrocatalyst197. This is illustrated in Table 7. The inclusion of single Pd or Pt atoms leads to O2 adsorption energies and charge transfers that are between those of the bare and fully terminated Nb2C MXene. Furthermore, if the single atoms are adsorbed on T vacancies, effectively partially replacing the T atoms, the O2 adsorption energy nearly doubles and the OER and ORR overpotentials are nearly halved, ultimately causing the Pt SACs to display the lowest overpotentials for both reactions, especially when located at O or F vacancies of Nb2CO2 or Nb2CF2, respectively.

Table 7 Adsorption energies (({E}_{{rm{ads}}})) and charge transfers (Q) of O2 on Pt-or Pd-doped Nb2CT2 MXene surfaces, and respective OER and ORR overpotentials, ({eta }_{{rm{OER}}}) and ({eta }_{{rm{ORR}}}), respectively
Full size table

SACs have also been built by replacing an atom of the M element of a MXene with an atom of another transition metal element198,199. For example, Chen et al. did a systematic DFT study of the effect of substitutional transition metal atoms on Mo2CO2 in the activity of this MXene towards the ORR200. Even though Mo2CO2 does not adsorb O2, replacing one Mo with another transition metal can lead to O2 adsorption energies up to −1.31 eV (for substitutional V). Substitutional Mn, Fe, or Ni atoms were projected to be the best for ORR, with overpotentials between 0.27 and 0.35 V and moderate O2 adsorption energies between −0.32 and −0.50 eV, as required by the Sabatier principle. The employment of SACs here has several advantages, such as providing a middle-ground between the reactivity of bare and terminated MXenes, and allowing more efficient use of noble metal atoms—100 % on Nb2CO2-VO-Pt and Nb2CF2-VF-Pt vs. 33.3% on Pt(111) or IrO2(110), the benchmark ORR catalysts.

Carbon dioxide reduction

Carbon dioxide can be electrochemically reduced through many mechanisms and into different products. Motivated partly by the fact that the bare Nb2N MXene surface had been predicted to bind to transition metal atoms201, CO2 electroreduction process was studied by Lu et al. over SACs (V, Cr, Mn, Fe, Co, or Ni) on the Nb2NO2 MXene surface, using the PBE-D3 approach202. They found that CO2 adsorbs preferably on the transition metal atom with moderate adsorption energies, ranging between −0.77 and −0.30 eV, which are suitable for catalysis, and with moderate charge transfers, ranging from 0.60 to 0.25 electrons, as shown in Table 8. Furthermore, the analysis was simplified by the fact that intermediate species with more than one carbon atom are not stable over SACs, so that only mono-carbon species had to be considered. Notably, the Co and Ni SACs exhibit ({G}_{max }) lower than and similar to, respectively, the limiting potential of 0.74 V of the state-of-the-art Cu(211) catalyst203.

Table 8 Adsorption energies (({E}_{{rm{ads}}})) of CO2 and CH4, charge transfers (Q) upon CO2 adsorption, and maximum Gibbs free energy spans (({G}_{max })) of CO2 reduction to CH4 on the Nb2NO2 MXene surface doped with one atom of V, Cr, Mn, Fe, Co or Ni, according to the results of the calculations taken from ref.202
Full size table

These results can be compared with analogous ones obtained for the CO2 electrochemical reduction on SACs over the Ti2CO2 MXene surface204. Zheng et al. found that the catalytic activity of such SACs can be described by the valence state of the transition metal adatom and observed a volcano-shaped relationship between these two quantities, peaked around 0.8 electrons for Co and Ni. The same transition metal elements were indeed estimated to yield the best catalytic activity on Nb2NO2 as well, albeit with a different peak charge transferred of circa 0.5 electrons202. This reaction has also been studied on 10 SACs deposited on the Ti2CN2 MXene surface205. There, Sc-, Ti-, V-, Mn- and Fe-based SACs were predicted to effectively originate CO with overpotentials below 0.5 V, but with poor CO selectivity on the first three transition metals. However, N-functionalized MXenes may be too reactive towards transition metal atom adsorption98.

The effect of the surface termination in the CO2RR was analyzed very recently by Meng et al.206. These authors considered the Ti3C2 MXene with several combinations of the O, OH and F surface terminations, viz. OH, (OH)2/3O1/3, (OH)1/2O1/2, (OH)1/3O2/3 and F1/3(OH)1/3O1/3, which differs from previous studies that focused on uniformly functionalized MXenes. The OH groups were found to serve mainly as H-donors, and the presence of both OH and O groups benefited the catalysis by reducing the limiting potential, in agreement with experimental results available in the literature.

Hydrazine oxidation

The electrochemical hydrazine oxidation follows the global reaction

$${{rm{N}}}_{2}{{rm{H}}}_{4}to {{rm{N}}}_{2}+4{{rm{H}}}^{+}+4{e}^{-}$$
(13)

with standard potential -0.33 V. Zhou et al. reported a joint experimental and theoretical study of the activity of a Ni SAC based on the Ti3C2Tx MXene towards this reaction207.

Experimentally, it was observed that Ni atoms occupy Ti vacancies on the surface of the MXene and that the activity of the Ni SAC surpasses that of Ni nanoparticles, both deposited on the same Ti3C2Tx MXene.

The DFT-based calculations, using the PBE exchange-correlation functional, determined that an Ni atom deposited on a Ti vacancy of the Ti3C2Tx MXene becomes strongly coupled to the neighboring C atoms. In turn, the presence of Ni at a Ti vacancy of the Ti3C2Tx MXene was found to cause a downshift in the d-band center of the C atoms from 1.74 eV to 1.52 eV, allowing for an optimal activation energy of N2H4. This molecule was adsorbed preferably over the Ni single atom, with a moderate adsorption energy of -0.69 eV, which suggests a stronger adsorption than on Ni nanoparticles (Eads = -0.55 eV). The overall activation energy for hydrazine oxidation was estimated at 0.45 eV, or 0.13 eV lower than on Ni nanoparticles deposited on the same MXene, supporting the experimental findings that a Ni single atom prevails over Ni nanoparticles for hydrazine oxidation.

Predictions by machine learning

The MXene family of materials includes a wide range of elemental compositions combined with different surface terminations, which currently leaves the majority of possible MXenes undiscovered or uncharacterized. Traditional experimental trial and error takes a long time to unveil the full potential of these materials. Therefore, DFT calculations have been employed to obtain descriptors and output properties associated with their production and application and, more recently, machine learning (ML) strategies are being developed to digest high-throughput DFT data, and thus accelerate the discovery of new and synthesizable MXenes for targeted applications (Fig. 9)208. As a result, an account of the current state of the art regarding machine learning studies that target the MXene family of materials is provided herein, with the main details of the ML approaches presented in Table 9.

Fig. 9: Computational modeling strategy.
figure 9

DFT and ML complementarity to understand and design applied MXenes.

Full size image
Table 9 Details of the machine learning (ML) approaches used to model and predict different MXene outputsa
Full size table

Synthesizability and stability

A central issue for the application of MXenes is their synthesizability. To tackle this problem, Frey et al.209 employed a positive and unlabeled (PU) machine learning framework, adapting it to predict the likelihood of successfully obtaining these materials. Besides the MXene family, they also built a model for their layered precursor materials known as MAX phases. Elemental information and data from high-throughput DFT calculations were utilized to characterize each material. After a feature selection step based on principal component analysis, the model hyperparameters were optimized. Two positive and unlabeled learning algorithms were employed, with the most helpful being an implementation of transductive bagging coupled with a decision tree classifier, chosen for its interpretability. The model generated synthesizability scores ranging from 0 to 1 for unlabeled samples, with the analysis revealing that features related to thermodynamic stability, bond strength, and charge distribution had the most significant impact on predictions. Notably, the predictions of the model aligned well with experimental findings, indicating that earlier group transition metal MXene and MAX phases with lower charge densities of the M atoms were more likely to be synthesized. The positive and unlabeled learning model successfully identified 111 MAX phases and 18 MXenes as synthesizable materials, including novel systems, such as Hf4N3, Sc3C2, and W4C3, even though monometallic MXenes based on Hf, Sc, or W had not been synthesized yet. Moreover, by considering both MAX phase synthesizability and the MXene score, the authors identified the top 20 most promising MAX systems that could be synthesized and subsequently etched to produce previously unavailable MXenes. Their results were validated based on 10-fold cross-validation, and, in time, these predictions of promising systems allow for a litmus test on their work.

He and Zhang210 explored the stability of Mn+1Xn MXenes, testing several machine learning algorithms and symbolic regression methods to assess and classify MXene material stability, based on a relatively small dataset of 85 materials. The output stability value was defined as the convex hull distance, according to the formula Estability = EformationHf, where Eformation is the formation energy and Hf is the convex hull energy, following the definition from other authors211. After using Pearson regression to select the most appropriate descriptors, random forest, k-nearest neighbors, logistic regression, support vector machine, and gaussian naive bayes were tested, with support vector machine emerging as the most accurate method for classifying materials based on their stability. Beyond the supervised learning classification approach, one interesting aspect of the study was the use of symbolic regression. Symbolic regression aims to find a symbolic expression or mathematical equation that best fits the data, essentially modeling the relationship between the input features and the outcomes. This symbolic expression can then be used to predict new data points or for further analysis. Using this methodology, the authors identified the key descriptors and generated an indicator (from the combination of other descriptors) that can be associated with MXene material stability. In this study, (χ(M)-RvdW(M))3 was the new and most important descriptor, obtained from symbolic regression, where RvdW(M) is the atomic van der Waals radius of the metal, and χ(M) corresponds to its electronegativity. Note that the greater variability on the considered Mn+1Xn MXenes is associated with the nature of the M element; hence, it is not surprising that the key indicators prioritize data for the M element over the X element.

Electronic properties

Rajan et al.54 delved into the electronic properties of MXenes, especially their band gaps, which are important to identify promising candidates for electronic, optoelectronic, and photocatalytic applications. The authors employed kernel ridge regression, support vector machine, Gaussian process, and bootstrap aggregating, using as input variables readily available elemental properties and DFT-optimized structural and electronic properties. From 23 870 generated structures, 7200 were randomly selected to obtain all the DFT and elemental properties necessary to build an initial metal−semiconductor classification model, which filters out finite gap MXenes with 94 % accuracy using the bootstrap aggregation algorithm. From these, 70 were selected as semiconductors, and their band gaps predicted using a ML regression model with the Gaussian process algorithm, which yielded the best results. The authors were able to address the band gap underestimation issue commonly found when using local and semi-local functionals in DFT calculations, without resorting to more time-intensive approaches, such as GW, despite still requiring the calculation of DFT based structural and electronic descriptors.

Besides the many possible MXene compositions and their different terminations, the functionalization of MXenes can result in tens of thousands of possible materials, with potential applications spanning from electronics to energy-related fields. For these types of applications, MXene properties depend on the position of their band edges. Therefore, their determination can allow the tailoring of their structure for specific applications. However, computing the position of these bands can be time-intensive, especially taking into account all the possible configurations that could be considered. To address this challenge, Mishra et al.212 developed a machine-learning model capable of determining the position of the valence or conduction band edges, using gaussian process regression, with a level of accuracy akin to GW calculations, achieving a minimal root-mean-squared error of 0.12 eV. Their model is based on PBE band edges and vacuum potential together with other DFT (formation energy, ground state phase) and elemental (standard deviations of atomic radius, average of molar volume for all the species, and heat of vaporization) descriptors, correlating with GW-derived band edges. Their work also pointed to the role of the electronegativity difference between the functional groups and the transition metals in determining band edges and fine-tuning the electronic properties of MXenes. This factor was considered as a descriptor in the model and reveals the role of creating surface dipoles due to MXene functionalization on the band positioning and modification of vacuum potential.

An electronic surface property that can influence the application of MXenes in catalysis, electronics, and energy storage is the work function. Using a dataset of 275 MXenes encompassing various compositions and surface terminations, including O, OH, F, or none. Roy et al.213 developed a series of ML models to identify MXenes with desirable work function values. Based on 15 fundamental descriptors for the elements and using artificial neural networks, the authors were able to estimate work function values with an associated mean absolute error value of 0.25 eV for the test set. Feature importance analysis revealed that the work function is mostly influenced by the properties of atoms terminating the MXene surface, namely their electronegativity, mutually validating the conclusions obtained years before by Mishra et al.212. In order to increase the transferability of their model to new materials, the authors developed reduced-order models with ten, eight, and five features, which were able to predict the work function of MXenes with various surface terminations beyond the original training set, such as Br, Cl, S, N, and NH (including extra data related to these new models in comparison with the initial dataset), with only a minor increase in mean absolute error.

Mechanical properties

Besides their unique electronic properties, MXenes are also known to be mechanically stronger than other 2D materials like graphene oxide (GO)214,215, being applied as part of protective coatings, composites and membranes215. Venturi et al.214 developed a ML approach to avoid the computation of descriptors of materials, which might require more expensive DFT calculations. Therefore, they used crystal graph convolutional neural networks on systems with planar periodicity (also containing systems besides MXenes), to calculate their thermodynamic, mechanical, and electronic properties. Their screening process considered nearly 45 000 material structures, focusing on two end uses: mechanical strength and photovoltaics. By analyzing statistical data derived from their screening, they unveiled the main structural and compositional design principles that influence the properties of the materials. Their approach reciprocally validated well-established design principles and revealed novel insights. For instance, they reaffirmed that hybrid organic-inorganic perovskites containing lead or tin are promising candidates for solar cell applications, in line with existing literature. Additionally, they identified titanium-based MXenes as having high stiffness coefficients, aligning with their reputation as mechanically robust materials. Moreover, their investigation uncovered promising untested compositions of all-inorganic perovskites, offering unexplored avenues for research in the field of photovoltaics. Lastly, they made the code available, which can foster further developments in the field. Their work has shown the capabilities of crystal graph convolutional neural networks, and the design principles presented by the authors can guide further experimental and computational research.

Tian et al.216 conducted a machine learning-based study focused specifically on the mechanical properties of MXenes, in particular, on the tensile stiffness and strength of 157 types of MXenes, obtained from DFT calculations. Their results reveal a wide range of tensile stiffness values for MXenes, with 42 structures surpassing the stiffness of materials known for their outstanding mechanical properties, such as graphene and transition metal dichalcogenide monolayers. The mechanical properties of MXenes were found to be strongly influenced by their thickness, bond strength, and surface terminations. Remarkably, similarly to the electronic properties, the surface terminations were identified as a key factor in enhancing the tensile stiffness, which can be twice as high for some surface terminations in comparison with bare MXenes. The authors used a machine learning technique for high-dimensional data analysis and feature selection, named Sure Independence Screening and Sparsifying Operator (SISSO), that allows to obtain interpretable relationships between descriptors and outputs. As a result, based on physically interpretable descriptors, the following analytical formula for predicting tensile stiffness (E2D) was obtained:

$${E}^{2D}={E}_{{rm{MX}}}+{E}_{{rm{coupling}}}+{E}_{rm{T}}$$
(14)

where the three quantities on the righthand side of the equation are also derived formulas themselves. EMX represents the effect of thickness and bond stiffness and depends on the MXene thickness and M–X bond stiffness, Ecoupling represents the weakening effect of surface terminations on internal M–X bonds, therefore it is a function of both M–X and M–T bond stiffness, and ET is the contribution of the stiffness of the M–T bonds. Furthermore, the authors found that the tensile strength (σs) scales linearly with the tensile stiffness.

Magnetic properties

Two-dimensional ferromagnetic materials have been playing an increasingly essential role in low-dimensional spintronics217. Since MXenes are emerging as one of the largest families of two-dimensional materials, Song et al.218 developed an ML model to classify the magnetic states of these materials, starting with a database of 23,825 materials. Their approach combined high-throughput DFT calculations with machine learning techniques, and introduced a new descriptor, derived from the averaging of matrices with information about elemental properties, crystal structure (connectivity), and Coulomb interactions. After developing and evaluating their model on the larger dataset, the authors conducted DFT calculations with the PBE functional, generating structures and properties for 4153 new MXene materials, from which 753 materials were found to possess high saturation magnetization (a saturation magnetization criterion of 15 μBnm−2 was employed to classify magnetic from non-magnetic materials). Then, they employed the previously developed ML classification model (AdaBoost, cf. Table 9) to purpose 1432 promising 2D magnetic MXene materials, that should possess high saturation magnetization. Note that the quality of the ML data is strongly determined by the choice of the DFT approach (please see section “Magnetic ordering” above). The PBE functional may struggle to accurately describe magnetic materials, especially those with strong electron-electron interactions, and hybrid functionals are typically regarded as more effective choices for accurately capturing magnetic properties.

Cytotoxicity

For the industrial application of MXenes, the toxicological profile of these materials must be evaluated. Traditional in vitro toxicological studies are time-consuming and expensive, therefore predictive models, such as those being developed with the aid of ML techniques, are often employed in some step of the development process and life cycle analysis of materials. Marchwiany et al.219 identified key surface-specific features associated with the cytotoxic behavior of 2D materials and leveraged machine learning to predict the toxicity for MXenes not previously tested in vitro. They highlighted two factors contributing to the cytotoxicity of 2D materials, notably the presence of transition metal oxides and lithium atoms on the surface. Furthermore, their research underscores the importance of surface characteristics to the toxicity of MXenes, while other structural properties have minimal impact. The machine-learning model successfully assessed the cytotoxicity of 19 MXenes, identifying two as cytotoxic with almost 90% probability, while the rest are considered non-toxic (cytotoxicity probability lower than 20 %) and suitable for various technological applications. The two MXenes that were predicted to be toxic are based on Ti3C2 and are the only ones with metal (titanium) oxide on their surface. Additionally, the research methodology and machine learning models developed by the authors can be extended to other complex 2D materials, including emerging materials like 2D transition metal borides (MBenes)9 and van der Waals heterostructures220.

Hydrogen evolution reaction

For the high-throughput screening of potential hydrogen evolution reaction (HER) MXene catalysts, Zheng et al.221 employed a combination of ML models and DFT calculations, to predict the Gibbs free energy of hydrogen adsorption (ΔGH*), which, according to the authors, is an accepted proxy of HER performance for various catalysts. First, the authors employed four different ML regression algorithms, including Elman artificial neural networks, kernel ridge regression, support vector machine, and random forest, based on DFT calculations for a total of 132 MXenes. The random forest algorithm gave the best performance and, in comparison with other algorithms, also had the advantage of its performance being the same when using DFT-based descriptors (lattice parameter, bond lengths, Bader charges, etc.) together with elemental properties (tabulated properties for the M, X, and T atoms included in the MXene structure), or using solely elemental descriptors with no DFT descriptors. Therefore, the authors relied on the random forest model with elemental properties alone to make further predictions, since it was faster to evaluate new prospective MXenes catalysts. Their work considered bare and S-terminated MXenes and identified Os2B and the S-terminated Scn+1Nn (n = 1, 2, 3) MXenes as active catalysts, since ΔGH* approaches zero for a broad range of hydrogen coverages (from 1/9 to 4/9 monolayer). Furthermore, the authors highlighted the role of S functional groups in modulating HER performance, causing the weakening of the adsorption of H*, which is a pivotal step in the HER process. Furthermore, ΔGH* values were predicted using the developed random forest model considering only elemental properties for 167 new MXenes. This was in addition to the 132 tested initially that also included DFT descriptors, giving a total of 299 MXenes evaluated in this work, and resulting in the identification of 8 new MXenes as potential HER catalysts. This work included a thorough validation and testing methodology, based on tenfold cross-validation, leaving 25% of the dataset for testing, which was randomized 100 times, while also considering leave-one-cluster-out cross-validation, with clusters depending on the physical similarities of the materials.

Wang et al.222 used high-throughput DFT calculations and ML to identify trends in the activity of oxygen-terminated monometallic and ordered bimetallic MXenes as HER catalysts. The materials considered by the authors were evaluated in terms of catalytic activity, thermal stability, and conductivity computations. Their work identified 110 experimentally unexplored MXenes with promising HER activity and thermal stability that surpasses that of the ideal platinum catalyst. Interestingly, titanium was featured in most of the predicted promising bimetallic MXene catalysts, which is supported by experimental results. Furthermore, their work demonstrated that the development of descriptors based on the AdaBoost ensemble learning model is able to accurately predict HER activity and unveil the fundamental geometric and chemical origins of the process, which are also consistent with electronic insights. The authors relied on tenfold cross-validation to optimize the model, and the materials predicted to surpass the ideal platinum catalyst, when they are obtained experimentally, can be used to further evaluate their predictions in the future.

Liang et al.223 explored the effect on the MXene electrocatalytic performance of surface termination tuning and the incorporation of single transition metal (TM) atoms, using DFT) calculations as the primary tool. Their research identified 21 previously uninvestigated TM@MXene catalysts that exhibit superior electrocatalytic activity compared to platinum, among a pool of 264 promising candidates. Furthermore, 7 of these 21 promising HER catalysts were predicted to have both thermal and thermodynamic stability, underlining their potential for practical application. Interestingly, these 7 MXenes are all based on Ti3C2, with varying surface termination and TM doping. This is of interest for practical applications, since Ti3C2-based MXenes are the ones whose method of synthesis has been known for the longest time. Additionally, the authors verified that ML techniques were able to predict the catalyst activity (ΔGH*) and thermal stability (cohesive energy) using only easily accessible elemental properties of single atoms and surfaces group elements, such as molar volumes, atomic radius, boiling points, and electron affinities, for ΔGH*, and the third ionization, ionic radius, atomic mass, crystal radius, and molar volume for the cohesive energy.

Oxygen reduction and evolution reactions

The oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) are critical for electrochemical conversion processes that play a role in various clean energy conversion and energy storage technologies, which are pivotal to unveil sustainable energy solutions. Despite their importance, ORR and OER suffer from slow kinetics and require catalysts such as platinum, ruthenium dioxide or iridium dioxide, to occur. However, the price of the metals in these catalysts prevents their widespread adoption for large-scale industrial applications. Therefore, Chen et al.224 explored the catalytic potential of M-N4-Gr/MXene heterojunction nanosheets for ORR and OER, through a combination of DFT calculations and ML predictions. The M-N4-Gr layer is based on graphene (Gr), replacing two neighboring carbon atoms by one atom of a metallic (M) element and the four surrounding carbons with nitrogens (N4). A total of 78 heterostructure nanosheets, comprising the Ti2C, Nb2C, and V2C MXenes and M-N4-Gr implanted with each of 26 different transition metal atoms, were examined. The random forest regression algorithm demonstrated the best performance among the tested ML methods, resulting in reduced root mean squared root errors. This study also highlighted the potential of Ni-N4-Gr/Nb2C and Ni-N4-Gr/V2C for ORR, Ru-N4-Gr/Nb2C for OER, whereas Co-N4-Gr/V2C was predicted to be a promising bifunctional catalyst. The authors also identified the d-band center of the metal active site and the amount of charge transfer from the metal in each structure as key descriptors, from the correlation between these descriptors and the overpotential, thus providing valuable insights for future catalyst design.

Ma et al.225 studied the catalytic activity of Pt-doped dual transition metal Janus MXenes, which, in this case, are MXenes that have different types of transition metal elements on each side of the surface. The authors took advantage of this to design bifunctional catalysts for both ORR and OER. Through a combination of first-principles calculations and ML models, their study explores the underlying physical and chemical factors that influence the catalytic overpotential. The dynamic and thermal stability of Janus MXenes was initially evaluated using cohesive energies, phonon dispersion, and ab initio molecular dynamics simulations. Subsequently, the electronic properties of Pt-doped Janus MXenes were examined. The analysis confirmed the stable support of Pt atoms on the oxygen vacancies (Pt-VO) of dual transition metal MXenes, while the density of states and projected density of states calculations indicated their favorable electrical conductivity as single-atom catalysts. Building upon these findings, high-performance single-atom catalysts with minimal overpotentials for ORR and OER were identified, such as Pt-VO-CrTiCO2, Pt-VO-MnTiCO2 and Pt-VO-PdTiCO2. Pt-VO-MnTiCO2, in particular, showed the most significant improvement compared to Pt-VO-Ti2CO2, which was considered as a reference system in this study. This achievement was credited to the adjustment of single-atom catalysts electronic properties through dual transition metal modifications. A key aspect of this study was the consideration of machine-learning models (random forest regression and gradient boosted regression) to pinpoint descriptors crucial to the catalytic activity of MXenes-based single-atom catalysts for ORR/OER, concluding that the adsorption energy of OH was the most important descriptor for both ML models as well as for both reactions. Therefore, as future work, the authors propose to focus on the adsorption of OH on the catalyst surface, increasing the dataset and building a new ML model for this purpose. Moreover, the binding energy of Pt onto the substrate and the d-band center of the Pt atom were also found to be important descriptors to model both reactions. Hence, this study provides valuable guidance for the development of novel and efficient bifunctional Janus-MXenes electrocatalysts for ORR and OER.

Energy storage

The mass adoption of electric cars and the critical role of wireless electronic devices in our lives are putting increasing pressure on lithium prices, since it is ubiquitous in batteries due to its excellent balance of energy density, power density, and life cycle. This has led to the development of new high-capacity battery materials that can reduce the dependence on lithium, and MXenes are also candidates for this application. Li and Barnard226 explored the electrochemical energy storage potential of MXenes, in order to establish a relationship between MXene composition and various electrochemical properties. The authors transformed the MXene chemical formulas into categorical descriptors, enabling the simultaneous prediction of multiple electrochemical properties. The ML workflow developed incorporated multi-target regression techniques, to predict the gravimetric capacity, voltage, and induced charge, emphasizing the physicochemical characteristics crucial for battery design. Furthermore, this study involved an inverse design approach by predicting MXene formulas based on predefined battery performance criteria, using multi-target classification techniques. The MXene surfaces studied were of the form Z2M2CT2, where the MXene transition metal element is M = Sc, Ti, V, Cr, Zr, Nb, Mo, Hf, or Ta; the surface termination is T = H, O, OH, F, or none; and an extra ion intercalation layer was considered, Z = Li, Na, K, or Mg. The resulting inverse model suggests Li2M2C and Mg2M2C (M = Sc, Ti, Cr) as promising candidates for more detailed analysis using DFT, according to the predicted gravimetric capacity, voltage, and induced charge properties. On the one hand, the more traditional forward predictions, based on a regression model, allow to establish guidelines for the discovery of new battery materials, identifying as important factors the intercalating ions, transition metals, and functional groups. On the other hand, the inverse model, based on classification, is a tool to predict new battery materials with targeted specifications, while considering three important variables for final application within the same model (multi-target ML). This study shows the capability of the inverse model to classify multi-class materials accurately, expanding the scope of inverse design and offering specific chemical formulas as potential outputs.

Hydrogen storage

Cheng et al.227 performed a computational study exploring the hydrogen storage capabilities of MXene multilayers, motivated by growing interest in hydrogen as a clean energy source and the potential of MXenes multilayers for hydrogen storage, as highlighted by Liu et al.228. Cheng et al. developed a workflow to screen a large number of MXene compounds (23,857 in total, with 12,647 used to train and validate the model) to investigate the relationship between the activated H2 bond length and the adsorption distance. The methodology involved the use of DFT to generate a dataset concerning the adsorption geometries of hydrogen on MXenes, which allowed to understand how hydrogen molecules interact with the MXene surface. Furthermore, the study employed machine-learning techniques, specifically physics-informed atomistic line graph neural networks (ALIGNNs), to predict adsorption parameters. ALIGNNs are a type of neural network designed to model physical systems by incorporating known laws of physics into their structure, which enhances their predictive accuracy and interpretability. They were interpreted within the framework of Pauling’s resonating valence bond (RVB) theory, which describes the quantum mechanical phenomenon where electrons are not localized on individual atoms in a molecule or a crystal lattice but are spread over several atoms, forming a resonance structure. The main conclusion of the study was the derivation of a general formula that describes the relationship between the activated H2 bond length and the adsorption distance on MXenes. This formula, based on DFT, machine learning, and RVB theory, indicates how the ligancy (i.e., coordination number) and valence of transition metals in MXenes affect the activation of dihydrogen, which is a crucial aspect of hydrogen storage.

CO2 activation and CO sensing

In order to create added value by removing the excessive carbon dioxide concentration from the atmosphere of the Earth, CO2 should be utilized at an industrial level instead of just stored. Abraham et al.229 combined high-throughput DFT and ML methodologies to identify the main descriptors associated with the activation of CO2 over MXene surfaces. They analyzed a dataset comprising 114 pure and defective MXenes, offering insights into their adsorption energy, with random forest regression emerging as the best-performing method. Based on feature importance analysis embedded in random forest, the authors found that surface properties such as d-band center, surface metal electronegativity, and the valence electron number of metallic atoms were the key indicators to predict CO2 activation. This work established a framework for the design of MXene-based catalysts, facilitating the prediction of potential markers for CO2 activation and its subsequent utilization.

Carbon monoxide detection is vital for environmental safety, since it is colorless and odorless, but poses a high health risk. Due to their surface accessibility, electrical conductivity, stability, and surface functionalization versatility, MXenes are ideal candidates for gas sensing applications. O- and S- terminated MXenes exhibit limited interaction with CO, but it can be significantly enhanced through transition metal decoration. Therefore, Boonpalit et al.230 presented an innovative approach involving the high-throughput screening of 450 combinations of transition-metal decorated MXenes (TM@MXene) for CO sensing applications. This screening is facilitated by an integrated active learning and DFT pipeline developed by the authors. The active learning component harnesses a crystal graph convolutional neural network (CGCNN) as a surrogate model, enabling the swift identification of CO sensor candidates while minimizing computational resource demands. The CGCNN approach allows to increase the speed of virtual screening by skipping the iterative feature engineering step, due to its automatic feature extraction ability, whereas the active learning approach allows to reduce the number of required training samples compared to a fully supervised learning methodology, starting with a reduced database of 40 adsorption energies and requesting more DFT calculations as the learning process progresses until 450 TM@MXene were analyzed. Four criteria were imposed for suitability towards CO sensing: CO adsorption energy between −0.8 and −0.6 eV, negative surface formation energy, charge transfer of at least 0.15 electrons upon CO adsorption, and visible surface density of states change near the Fermi level in the presence of adsorbed CO. In the end, the authors identified Sc@Zr3C2O2 and Y@Zr3C2O2, the only systems predicted to satisfy all four criteria, as the most promising TM@MXene candidates for CO sensing.

C−N coupling

Electrochemical C−N coupling is a promising way to efficiently convert excess CO2/N2 or harmful N-containing NO3/NO2 anions into valuable products, such as urea, methylamine, and acetamide, for example. This approach is both environmentally beneficial and sustainable, since it decreases the dependence on fossil fuels, while promoting a more circular economy. It is known that the pivotal step in this process is the initial formation of the C−N bond, which determines the activity and selectivity of the C−N coupling. Therefore, Jiao et al.231 built activity and selectivity pathways for 54 MXene surfaces by employing DFT calculations. Their findings revealed that the activity of the C−N coupling step primarily hinges on the strength of *CO adsorption (Ead‑CO), while selectivity is influenced by the co-adsorption strengths of *N and *CO (Ead‑N and Ead‑CO). Consequently, they proposed that an ideal C−N coupling MXene catalyst should display moderate *CO adsorption and stable *N adsorption. Their study further incorporated ML techniques to derive data-driven correlations, describing the relationship between Ead‑CO and Ead‑N, using atomic physical chemistry features. By applying these formulas, the researchers efficiently screened 162 MXene materials without the need for more time-intensive DFT calculations. Several potential catalysts, including Ta2W2C3, were identified as promising for C−N coupling performance, and subsequently validated through further DFT calculations.

Solar spectral absorption

Ding et al.232 investigated the use of MXenes-based metasurface absorbers (MMA) to enhance solar spectral absorption, a crucial aspect of solar cell technology. Despite significant advancements in this field, conventional materials and design methods have not met practical requirements. MXenes, known for their spectral selectivity, are considered as promising solar-absorbing materials, with their optical properties closely linked to their chemistry. The work from these authors focused on examining the absorption characteristics of various MMA configurations with different terminal groups in the solar spectrum. To address this, the authors employed a random forest regression model, built specifically for these systems. The random forest method facilitated the entire process, from parameter selection to structural optimization, with a remarkable low error rate of only 0.02%, surpassing two other classical ML algorithms, linear regression and k-nearest neighbors, which were used for comparison. The simulation results demonstrated that the MMA with Ti3C2O2 achieved exceptional solar absorption, absorbing 93.2% of the solar spectrum and emitting a mere 1.2% in the mid-infrared band (5000–13,000 nm), making it an almost ideal solar absorber. Furthermore, the proposed MMA exhibited advantageous properties, such as polarization insensitivity and angular independence, all while maintaining a slender profile. Notably, the maximum absorptivity and corresponding structural parameters were predicted directly by the random forest regression model, eliminating the complexities and time-consuming aspects associated with traditional manual design processes. This approach points in the direction regarding the design of solar absorber structures.

Machine-learning features

The majority of the machine-learning works targeting MXenes reviewed in this work (Table 9) focused on developing workflows to predict numerical properties (regression), with random forest being the most common selected algorithm. The median dataset had 287 examples and the average number of features among the works that included feature selection was 12, which, by our experience, constitutes a reasonable proportion of selected features to total number examples (~4%). In the works that calculated the feature importance for the respective models, cf. Table 10 and Fig. 10, it is possible to observe that most of the important properties were elemental properties, such as electronegativity, ionization energy, atomic mass, atomic radius, boiling point, valence electrons, molar volume, etc. Properties related to the elements present on the external layer (surface terminations) of the MXene material were also key, since many predicted processes and properties are promoted by the interface of these materials. The most important DFT property that appeared as input was the formation energy.

Table 10 Number of features considered in studies devoted to the prediction of the properties of MXenes by machine learning and, when available, the ones that were considered more important
Full size table
Fig. 10: Word cloud.
figure 10

Word cloud of the most important feature descriptions in the ML works that had this information available.

Full size image

Future perspectives

As the computational power increases, molecular models used in the simulations are becoming more realistic. The inclusion of the MXene surface termination in the computational models was one such improvement, since MXene functionalization, aside from more accurately matching the experimentally observed material, is desired for catalytic applications because bare MXenes do not satisfy the Sabatier principle for most molecules, as they interact too strongly. Another crucial advance of the computational models of MXenes was the consideration of metallic atoms other than the M element of the MXene, both as double transition metal MXenes and as metal atom dopants on the surface to form single-atom catalysts.

The number of MXenes predicted to be stable is expected to increase exponentially each time a new experimental or theoretical discovery increases the tuneability of MXenes. Examples of this are the synthesis of bimetallic MXenes, the observation of previously unknown possible MXene surface terminations and mixtures thereof, or the identification of oxygen within the carbon layers of MXenes, which suggested that most carbide MXenes are in fact oxycarbides.

We expect to see a rise in the number of studies on MXene-deposited single-atom catalysts, since the deposition of single atoms of a metallic element has been shown in many contexts and for many MXenes to lead to an intermediate MXene-adsorbate interaction strength, between the easily poisoned bare MXenes and the almost inactive fully terminated ones. Thus, MXene-based single-atom catalysts seem very promising to make use of the reactivity of the M layer of MXenes towards adsorption and bond activation, while still allowing the surface termination to limit the adsorption strength of reaction intermediates, following the Sabatier principle. Also, the dynamics of mixtures of surface terminations, of ions and solvent molecules within multiple layers of MXenes, and their role in the properties displayed by MXenes, are difficult to prove by experimental approaches only and, consequently, future computational studies are needed, conceivably taking advantage of on-the-fly machine-learned potentials. Another underexplored area is related with defective MXene surfaces (vacancies, areas devoid of terminations).

In order to enhance our understanding of the properties of MXenes and to improve the discovery of new MXenes, it will be essential to integrate high-throughput DFT, data-driven methods, and ML approaches into a cohesive framework. Within this approach, the following strategies can be employed:

  1. 1.

    High-throughput DFT Methods: i) Develop automated pipelines for performing high-throughput DFT calculations, reducing manual intervention and speeding up the process of preparing calculations; ii) continue leveraging high-performance computing clusters to handle large-scale DFT calculations; iii) take advantage of improved DFT methods, with more accuracy to compute specific interactions and properties, while also allowing faster convergence and reduced computational cost, enabling the evaluation of a larger number of MXene possibilities.

  2. 2.

    Data Management: i) Establish a centralized and standardized database to store and manage the vast amount of MXene data that can be simulated using high-throughput DFT; ii) ensure data quality and consistency by implementing rigorous criteria for appropriate MXene DFT calculations; iii) employ data-mining techniques to extract meaningful MXene properties from the increasing number of publications targeting these materials and build a larger database.

  3. 3.

    Machine learning: i) Mutually validate DFT and experimental data using appropriate statistical criteria and/or machine learning methods. ii) use machine learning to identify the most relevant features that influence MXene properties, in order to target these features to be included in databases, thus contributing to improve future models; iii) implement active learning approaches where machine learning models iteratively suggest the most informative MXene systems to simulate or synthesize next, thus saving time and resources.

  4. 4.

    Integrated DFT, data and ML approaches: Create an automated closed-loop system where DFT simulations feed data into machine learning models, and the insights gained inform further simulations and experimental validations.

  5. 5.

    Collaboration and open-source: Foster the collaboration among scientists working with MXenes using different experimental and computational approaches through the use and production of open-source platforms, tools and shared databases.

  6. 6.

    The next MXene frontier: Explore combinations of elements and characteristics of MXene materials that are absent from current databases, and, thus, might not be recommended by expert knowledge or current ML models.

By integrating high-throughput DFT, data-driven methods, and ML approaches, the process of discovering new MXene-based materials can become more efficient and effective, leading to faster advancements in their realistic application.

Related Articles

Chalcogen and halogen surface termination coverage in MXenes—structure, stability, and properties

MXenes are a diverse family of two-dimensional (2D) transition metal carbides, nitrides, and carbonitrides. They can be synthesized through both top-down approaches, such as selective etching of A-layers from MAX phases using acids or molten salts, and bottom-up approaches, such as direct synthesis using chemical vapor deposition. However, the degree of the surface termination coverage depends on the synthesis route and is one key parameter for controlling its properties. This study focuses on halogen- and chalcogen-terminated MXenes, particularly M2CTx where M = Ti, Zr, V, Nb, Ta, and T = S, Se, Te, Cl, Br, I, and with surface termination coverage ranging from 100% (ideal, x = 2) to 50% (x = 1). The incorporation of oxygen on vacant termination sites was also evaluated. Using density functional theory (DFT) calculations, we investigated the structural, electronic, and mechanical properties of these MXenes. Our findings reveal that non-ideal termination coverage (x < 2) is more favorable for MXenes terminated with a larger size of T, such as Ti2CBrx, Nb2CClx, and Ta2CClx, and leads to mixed termination sites and lower binding energies. A reduced binding energy may facilitate delamination into single sheets, however, too low termination coverage may also cause structural collapse. Electronic properties showed an increased number of states at the Fermi level under non-ideal coverage, potentially enhancing the conductivity. Mechanically, we find the moduli of MXenes to be comparable to other 2D materials, such as transition metal chalcogenides and hexagonal boron nitride, indicating their suitability for applications requiring flexibility and durability. This study underscores the potential of tailoring MXene properties through precise control of termination coverage and composition, paving the way for enhanced application-specific performance.

Flash Joule heating for synthesis, upcycling and remediation

Electric heating methods are being developed and used to electrify industrial applications and lower their carbon emissions. Direct Joule resistive heating is an energy-efficient electric heating technique that has been widely tested at the bench scale and could replace some energy-intensive and carbon-intensive processes. In this Review, we discuss the use of flash Joule heating (FJH) in processes that are traditionally energy-intensive or carbon-intensive. FJH uses pulse current discharge to rapidly heat materials directly to a desired temperature; it has high-temperature capabilities (>3,000 °C), fast heating and cooling rates (>102 °C s−1), short duration (milliseconds to seconds) and high energy efficiency (~100%). Carbon materials and metastable inorganic materials can be synthesized using FJH from virgin materials and waste feedstocks. FJH is also applied in resource recovery (such as from e-waste) and waste upcycling. An emerging application is in environmental remediation, where FJH can be used to rapidly degrade perfluoroalkyl and polyfluoroalkyl substances and to remove or immobilize heavy metals in soil and solid wastes. Life-cycle and technoeconomic analyses suggest that FJH can reduce energy consumption and carbon emissions and be cost-efficient compared with existing methods. Bringing FJH to industrially relevant scales requires further equipment and engineering development.

Printable graphene inks with polypropylene carbonate for low-surface-tension solvents and mild-temperature post-processing

For dispersion stability, printable graphene inks commonly employ solvents with limited surface tensions or incorporate dispersant aids that require high-temperature post-processing, restricting printability and substrate compatibility. Here, printable graphene inks are introduced with low-surface-tension solvents and mild-temperature post-processing using polypropylene carbonate (PPC). Graphene is produced by liquid-phase exfoliation with PPC, and the exfoliated graphene/PPC is used to generate printable inks. As a dispersant aid, PPC improves graphene exfoliation, dispersion stability, and redispersability in solvents with low surface tensions (<30 mJ m–2), facilitating the formulation of desirable inks for efficient aerosol jet printing on diverse substrates. Moreover, the low decomposition temperature of PPC eases its thermal removal from printed graphene, allowing high electrical conductivity with a mild post-processing temperature of 220 °C. Consequently, the graphene inks enable the fabrication of fully-printed graphene micro-supercapacitors on heat-sensitive paper substrates, exhibiting high areal capacitances, cycling stability, and mechanical resilience against bending deformation.

Solution-processable 2D materials for monolithic 3D memory-sensing-computing platforms: opportunities and challenges

Solution-processable 2D materials (2DMs) are gaining attention for applications in logic, memory, and sensing devices. This review surveys recent advancements in memristors, transistors, and sensors using 2DMs, focusing on their charge transport mechanisms and integration into silicon CMOS platforms. We highlight key challenges posed by the material’s nanosheet morphology and defect dynamics and discuss future potential for monolithic 3D integration with CMOS technology.

Antibody-functionalized MXene-based electrochemical biosensor for point-of-care detection of vitamin D deficiency

Clinical studies routinely show that individuals suffer from vitamin D deficiency, which can result in health complications that include cardiovascular disease, autoimmune disorders, neurodegenerative diseases, and different skeletal deformities. Given its integral role in homeostasis and connection to many pathologies, early diagnosis of vitamin D deficiency is crucial. However, monitoring vitamin D levels is challenging, particularly in remote regions, due to the cost, time, and complexity of existing methods. Here, we develop an electrochemical biosensor for vitamin D based on antibody-functionalized MXenes, offering clinically relevant sensitivity, specificity, and amenability for point-of-care testing. Ti3C2Tx MXene nanosheets are amine-functionalized by electrostatically-driven modification with polyethylenimine, whose functionalities are then used for covalent conjugation of anti-vitamin D antibodies via glutaraldehyde chemistry. This platform achieves a detection limit of 1 pg mL−1 with a dynamic range (0.1–500 ng mL−1) that covers clinically relevant deficiency, insufficiency, sufficiency, and toxicity.

Responses

Your email address will not be published. Required fields are marked *