Nickel, the 28th chemical element of Mendeleev’s periodic table of chemical elements, is relatively widespread; its prevalence in the Earth’s crust is 0.009 wt% ^[1]. The study of the catalytic activity of nickel began in 1879 by the French chemists Paul Sabatier and Jean-Baptiste Senderens, and the reaction that first showed its activity, namely the direct one-stage hydrogenation of ethylene to ethane, is known as the Sabatier-Senderens reaction. Sabatier was the first to widely demonstrate the catalytic activity of nickel, which he described in detail in his famous book “Catalysis in Organic Chemistry” ^[2]. It seems that many researchers who have studied the catalytic reactions with this particular metal have heard Sabatier’s famous description of a nickel catalyst: “It can be compared to a spirited horse, delicate, difficult to control, and incapable of sustainable work”. However, few of them specify that later in his book, Sabatier determined that the pretreatment of nickel greatly affects its catalytic activity. Sabatier clearly formulated the conditions and features of the catalyst’s preparation, and the reactions of hydrogenation and dehydrogenation of organic compounds served as the starting point for the development of nickel catalysts.

Nowadays, nickel-based catalysts are widely used both at laboratory and industrial scales. Ethylene carbonylation in the presence of nickel carbonyl makes it possible to synthesize propionic acid (Reppe reaction) ^[3], which, in turn, is widely used in the production of drugs (ibuprofen, etc.), plastics (for example, polyvinyl propionate), surfactants, and also as a preservative E280. The presence of the Ni[(ArO)₃P]₄ catalyst ensures the hydrocyanation of 1,3–butanediene followed by the isomerization of the reaction products to adiponitrile ^[4]. The Raney Ni catalyst is used in many hydrogenation reactions ^[5], for instance, in the conversion of benzene to cyclohexane, in the production of sorbitol from glucose, and in the production of amines from nitro-containing compounds. In the examples given, nickel, with the exception of Raney Ni, acts as a homogeneous catalyst, which is usually inferior to heterogeneous catalysts in industrial processes due to the difficulties associated with separating the catalyst and the reaction mixture. However, homogeneous catalysts have an undeniable advantage since each Ni-containing molecule is a catalytically active center, which allows them to be used in a small concentration, in contrast to heterogeneous catalysts. To offset this shortcoming, the current development of heterogeneous catalysis is focused on single-atom catalysts (SACs).

The development of SACs is a trend of the last decade in modern catalysis, that is just beginning to gain momentum. Many scientists around the world are reporting that the use of catalytically active single atoms improves the efficiency of some reactions by increasing conversion, increasing selectivity to target products, and prolonging the operation of catalysts due to their exceptional stability. Carrying out a catalytic reaction on a specific single atom stabilized by a support makes it possible to reduce the loading of the metal as a whole. Such an approach allows one to reduce overall costs and save non-renewable resources (for example, the content of the most commonly used platinum group metals in the Earth’s crust is only about 10⁻⁶–10⁻⁸ wt% ^[1]). The change in catalytic activity is associated with the nature of active sites.

Figure 1 shows how metal-based catalysts can be classified by size, from bulk metal, which does not require the presence of a support, to single atom, which does. In each case presented, the surface atoms can catalyze the reaction; however, the particle surface structure, as well as the total number of metal atoms on the surface, will determine the catalytic activity. SincThe this review research is devoted to nickel, weresearchers should consider its crystal—a face-centered cubic lattice ^[6]. In this structure, each bulk Ni atom has a coordination number of 12 (saturated coordination), while the surface atoms located in the center of the cube faces are bonded to eight atoms, the atom in the middle of the edge has a coordination number of five, and the nickel atom located at the corner of the lattice is connected to three atoms. Thus, the highest degree of unsaturation is observed for the corner surface atoms, which makes them energetically and often catalytically more attractive. However, this reasoning is applicable only to particles with more than 13 atoms, having at least one internal bond with the maximum saturation, and having the Ni crystal lattice. In the case of amorphous particles or somehow ordered clusters consisting of less than 12 atoms, bond unsaturation for each atom, especially for surface ones, will increase, reaching a maximum for a single atom. Having the highest surface energy, a single Ni atom, like any other metal, needs to be strongly stabilized on the support to prevent its migration for the purpose of binding (for example, agglomeration) in order to reduce the free surface energy. Thus, with the transition from bulk metal to a single atom, there is a loss of stability, often accompanied by an increase in the activity of the catalytic center due to an increase in the unsaturation degree.

The study of metals in the atomic state is associated with difficulties caused by the resolution of characterization methods. Many reviews are devoted to the characterization of metal in a single-atom state [7,8,9,10,11,12,13,14]^{[7][8][9][10][11][12][13][14]}. In the simplest case, indirect confirmation can be used based on a set of methods that confirm the presence of metal in the sample, but its absence in the form of bulk or nanoparticles. However, researchers most often use the three methods depicted: X-ray photoelectron spectroscopy (XPS), electron microscopy, and X-ray absorption fine structure (XAFS) spectroscopy. XPS can be used to determine the electronic state of nickel. Using XAFS spectroscopy, it is possible to determine the nickel oxidation state (X-ray absorption near edge structure, XANES region) and atom environment (Fourier transform of extended X-ray absorption fine structure, EXAFS region). To determine the surrounding atoms, bond lengths, and coordination numbers, it is necessary to carry out additional fitting of the spectra. To directly show the monatomic distribution of the metal, the authors usually use electron microscopy. Thanks to the resolution achieved by modern microscopes, the single metal atoms and the lattice in which they are located can be observed. Thus, by combining the three characterization methods presented, the resulting metal single atoms can be fully investigated, which is necessary to further establish the ‘material-properties-application’ relationship.

Over the past few years, the number of publications devoted to Ni SACs has increased. More than half of them are related to the study of electrocatalytic reactions and corresponding density functional theory (DFT) calculations. This tendency is expected due to humanity’s attempts to avoid the production of energy partially or completely by releasing it from non-renewable, non-ecofriendly oil and gas, while many advanced clean energy technologies require highly active catalysts to lower the energy barrier and increase the reaction rate with an efficient and stable route.

2. Application of Supported Ni Single-Atom Catalysts

In the Introduction section, wresearche rs briefly discussed the historical and modern applications of Ni-based catalysts, which are mainly related to the catalysis of hydrogenation and dehydrogenation processes. When discussing catalysts for such reactions, one cannot fail to mention the competition between palladium and nickel, which was well covered in a recent review by Chernyshov and Ananikov [106]^[15]. Indeed, being one under the other in the same group of the periodic table, these metals are analogues; however, the difference in the periods causes differences in redox potentials and the availability of different oxidation states, which provides outstanding opportunities for Ni catalysis for complex multistep transformations and the discovery of new catalytic reactions. Thus, wresearchers will consider both the theoretical and practical results of using Ni SACs in various hydrogenation and dehydrogenation reactions. It should be noted that the development of Ni SACs’ application field is not much different from the historical ones for Ni bulk and Ni nanoparticles: this is the breaking of X–H bonds, where X is various atoms, for example, C or O, and H is a hydrogen atom. In this regard, for convenience, weresearchers will further consider separately the activity of heterogeneous catalysts based on Ni single atoms in breaking various X–H bonds, and we researchers will start with the C–H bond.

2.1. C–H Activation

Currently, the most common catalysts for the activation of C–H bonds are Ni-based catalysts. The activity of nickel is associated with the presence of an unfilled d orbital, which can accept the σ electron of the C–H bond, thus weakening or breaking it [107]^[16]. Such an approach to the activation of hydrocarbons can facilitate cross-coupling reactions, leveling the disadvantages known for this reaction: the use of an expensive catalyst (Pd-based), the presence of several stages of pre-functionalization, which entails a multi-stage reaction, and the presence of by-products. However, the activation of the C–H bond also has its drawbacks, in particular, high reaction temperature and high metal loading are needed. The best-known nickel-catalyzed reaction involving C–H activation is the dry reforming of methane (DRM).

Dry Reforming of Methane

DRM is a process for the simultaneous synthesis of H₂ and CO (syngas) via the conversion of carbon dioxide and methane. Even though this reaction may not actually involve hydrogenation and dehydrogenation, its mechanism includes the C–H bond cleavage as an important step of the reaction. The general reaction (1) describing this process is endothermic with a standard enthalpy of 247 kJ mol^–1 and, therefore, requires high temperatures (>500 °C) [59,60,108]^[17][18][19]. The mechanism of this reaction has been thoroughly studied [109,110]^[20][21]. The reaction begins with the adsorption of a methane molecule on the catalytically active center (2), followed by its activation (3), which is the rate-determining stage of the reaction [111,112]^[22][23]: where ‘(g)’ stands for the gas phase and ‘*’ here and below in the text means the adsorbed species. The methane molecule is a highly symmetrical nonpolar molecule, which should be considered when choosing a catalyst. As shown in the literature [113]^[24], the adsorption of methane on a metal single-atom site occurs due to the hybridization of bonding orbitals of the methane with antibonding orbitals of the metal and vice versa. It should be noted that a support on which metal single atoms are located may also play a role in methane adsorption. The stability of a metal atom is ensured by the strength of its bonding with a support, which is related to overlapping their orbitals. The lengthening of the C–H bond, leading to its further cleavage, depends on the adsorption energy of methane and on the charge given back from metal. This bond elongation is most noticeable for the Ni atom, which is the rationale for the effectiveness of Ni SACs in the DRM reaction. Indeed, according to the reports of the past few years [114^[25][26],115], Ni single-atom sites exhibit superior activity, in particular, as compared to Ni nanoparticles. However, their use often leads to catalyst deactivation, which is associated with the formation of carbon as a result of the Boudouard reaction (4) and the complete cracking of methane (5), catalyzed by nickel [110]^[21]. It is interesting that the equilibrium in the Boudouard reaction can be shifted towards the formation of CO if the chosen support can donate electrons to the nickel atom or Lewis basic sites, which are attractive for CO₂ adsorption, are introduced into the catalyst. The former factor was widely studied in [107]^[16], where Ni nanoparticles supported by CeO₂ were used as a catalyst after a reduction in hydrogen. As a result of the reduction, oxygen vacancies were formed, which easily released free electrons. These free electrons increased the d-electron density on the nickel atoms. Such electron-saturated atoms prevented the migration of σ–electrons of the CH₄ molecule to the d orbital of the Ni atom, thereby reducing the adsorption of CH₄ and, thus, the degree of carbon deposition, and also increasing the adsorption of CO₂ due to including its unoccupied π–orbital. Quite recently, these assumptions were partly confirmed for a similar Ni/M–CeO₂ system, where M = Mg, Co, or Zn [114]^[25]. In this case, CO₂ molecules were adsorbed on oxygen vacancies generated due to the presence of metal cations, whose ionic radius is smaller than that of Ce⁴⁺. The carbon dioxide molecule dissociated into carbon monoxide and surface active oxygen, which easily removed the deposited carbon, thereby maintaining high catalyst activity for more than 150 h at 800 °C. A similar effect was achieved by Ni and Ce co-doping of a hydroxyapatite support [60]^[18] and co-anchoring of Rh and Ni atoms on a CeO₂ support [59]^[17]. In addition, ZrO₂ possesses the same properties of producing oxygen vacancies [116]^[27]. Thus, the first way to modify Ni SACs is based on the selection of supports containing a large number of oxygen vacancies, which are considered to be closely linked to the content of surface adsorbed oxygen species. The second way to modify the catalyst by introducing additional basic centers into it is to add a metal and/or its oxide to the support. The creation of such centers will enhance the chemisorption of CO₂ on them [109]^[20]. As for the Zr and Ce atoms mentioned above, their addition to oxygen-containing supports also increased the mobility of oxygen ions on the surface, making Ni catalysts more stable [117]^[28]. As for other metals, there is a report of better catalytic activity and stability of Ni catalysts for 40 h at 700 °C, in which the support was partially coated with MgO [118]^[29]. Interestingly, such doping initiated better Ni dispersion on the support surface. However, this method did not work when using MgO as a support for Ni single-atoms [119]^[30].  The use of SACs, including Ni catalysts, is not sufficiently developed since this is a rather young area of catalysis. In this regard, there is no confirmation/denial in the literature of the promotion of the DRM reaction previously obtained for Ni nanoparticles in the case of single atoms (adding K [120]^[31], Ca [121]^[32], Fe [122]^[33], Sr [123]^[34] or using CaO, BaO [124]^[35], La₂O₃ [125]^[36], etc. as supports [126,127]^[37][38]). However, significant progress is expected in this field.

2.2. H–H Activation

Hydrogen is the simplest molecule, and its properties are the most studied. Since this pure resource is available in abundance at a very low cost, catalytic hydrogenation is a mainstream technology in both research and industry. Molecular hydrogen is not very active under ambient conditions, but many positively charged transition metal atoms are capable of bonding and activating H₂. As is well known, bulk Ni and Pd have a very strong hydrogen adsorption affinity, can easily form metal hydrides, and are widely considered good catalysts due to their excellent hydrogen solubility, corrosion resistance, and diffusivity. As indicated in the literature, the adsorption energy of the hydrogen molecule on the Ni single atom [128,129]^[39][40] is higher than that on bulk Ni [130]^[41]. The initiation of the hydrogenation reaction most often occurs with the dissociation of the hydrogen molecule. Homolytic decomposition yields adsorbed H* atoms, while heterolytic decomposition produces partially charged H^δ+ and H^δ– species. The type of decomposition, the energy barrier, and the overall energy profile of the reaction strongly depend on the structure of the catalytically active Ni site. For example, according to DFT calculations of a catalytically active Ni atom coordinated to three and four carbon atoms in single and double graphene vacancies, H₂ dissociation is an endothermic process with energy barriers of 0.69 and 0.33 eV, respectively [131]^[42]. As can be seen, the dissociation of the adsorbed H₂ is an energy-consuming process [128,131,132]^[39][42][43]. However, with an increase in the number of Ni atoms to three in a catalytically active site, the dissociation is almost barrier-free [130]^[41]. Interestingly, a similar effect was observed when considering the reverse reaction, where the recombination of hydrogen atoms into a molecule took place [133]^[44]. After decomposition, hydrogen atoms can spontaneously migrate to support atoms. This effect is called the hydrogen ‘spillover’ process, and it may take place in the metal-carbon support system [134]^[45]. Such behavior, in particular, was observed for some nitrogen-containing systems [128]^[39] and also for MoS₂ [132]^[43]. For instance, the addition of Ni single-atoms on MoS₂ support causes an excess of electron density on the nearest sulfur atoms, thus enhancing their activity towards hydrogen adsorption [32,49]^[46][47]. Although the activation of a hydrogen molecule can also occur due to an increase in the H–H bond length after adsorption onto the Ni atom. Furthermore, during the co-adsorption of the reagents on the metal atom, the breaking of this bond for the further course of the reaction will be easier.2.3. O–H Activation

2.3. O–H Activation

The use of hydrogen as a hydrogenating agent has some disadvantages associated with its gaseous state under ambient conditions. In the reactions discussed above, the typical hydrogen pressure was in the range of 1–3 MPa, which introduces difficulties in carrying out reactions both on laboratory and industrial scales. To overcome this, other molecules are used from which the hydrogen atom can be obtained, such as alcohols. Wang et al. studied the effect of the composition and structure of the Ni–N_xC_y site in graphene on the adsorption of isopropyl alcohol and the subsequent detachment of hydrogen from it [149]^[48]. The Ni–N₄ and Ni–N₃ sites turned out to be the most energetically favorable. It was shown that, despite the energetic preference for its formation, the flat structure of the Ni–N₄ site sterically prevents the dissociation of the O–H bond and makes the process highly endothermic. On the contrary, the Ni–N₃ site, in which the Ni atom protrudes above the graphene surface, not only easily adsorbs the isopropyl alcohol molecule due to the overlapping of the p-orbitals of the O atom and the d-orbitals of the Ni atom but also easily breaks the O–H bond.

2.4. N–H Activation

In addition to the reactions of the formation of hydrogen atoms or molecules by activation of C–H, H–H, and O–H bonds, which are widely covered in the literature, some authors considered the possibility of using single-atom Ni sites for activation of other bonds, for example, N–H. Feng et al. studied the mechanism of reduction of Br–C₆H₅–NO₂ to Br–C₆H₅–NH₂ at Ni–N₄ and Ni–N₃ sites using N₂H₄·H₂O as a hydrogen donor [155]^[49]. Since the former catalytically active site exhibited significantly higher activity than the second one, the energy values given below correspond to the Ni–N₃ site. Hydrazine was adsorbed with the formation of the Ni–N bond, and the N–H bond breaking energy was equal to 0.93 eV, with both N₂H₃* and H* species remaining coordinated to the nickel atom. The subsequent adsorption of the R–NO₂ molecule was exothermic. The limiting stage of the reaction was the chemo-selective reduction process of the −NOH group through the recombination of the adsorbed hydrogen atom and the hydrogen atom of the −OH group—with an activation energy of 0.96 eV. The presence of the Br substituent did not adversely affect the reduction of the nitro group; moreover, the authors calculated that 1.18 eV would be required to break the C–Br bond, which makes this process kinetically unfavorable. In practice, the catalyst achieved 100% conversion and selectivity at 60 °C within 3 h of operation, and the calculated TOF was higher than that for the catalyst containing Ni nanoparticles. It is interesting that theoretical calculations of the decomposition of N₂H₄ at the Ni–C₃ site did not show such good results as in the work [156]^[50]. The adsorption of the molecule on the catalytically active Ni atom still occurred due to the overlapping of the 3d and 2p orbitals of Ni and N, respectively. However, the authors calculated not only the dissociation of the N–H bond but also the dissociation of the N–N bond and showed that the second way is energetically preferable (0.86 and 0.51 eV, respectively). Moreover, the final state for the system after the N–H bond cleavage was endothermic, while after the N–N bond cleavage, it was exothermic, which indicates the potential unsuitability of the Ni–C₃ site as catalytically active in the production of H* from hydrazine. Although this site still showed potential activity and high selectivity in the formation of ammonia.