Atomic Clusters: Structure, Reactivity, Bonding, and Dynamics

Atomic clusters lie somewhere in between isolated atoms and extended solids with distinctly different reactivity patterns. They are known to be useful as catalysts facilitating several reactions of industrial importance. Various machine learning based techniques have been adopted in generating their global minimum energy structures. Bond-stretch isomerism, aromatic stabilization, Rener-Teller effect, improved superhalogen/superalkali properties, and electride characteristics are some of the hallmarks of these clusters. Different all-metal and nonmetal clusters exhibit a variety of aromatic characteristics. Some of these clusters are dynamically stable as exemplified through their fluxional behavior. Several of these cluster cavitands are found to be agents for effective confinement. The confined media cause drastic changes in bonding, reactivity, and other properties, for example, bonding between two noble gas atoms, and remarkable acceleration in the rate of a chemical reaction under confinement. They have potential to be good hydrogen storage materials and also to activate small molecules for various purposes. Many atomic clusters show exceptional opto-electronic, magnetic, and nonlinear optical properties. In this Review article, we intend to highlight all these aspects.


INTRODUCTION
A cluster is defined as a finite aggregation of atoms, starting with as low as two atoms and extending up to an upper bound of some hundred thousand atoms (Herr, 1993;Corrigan and Dehnen, 2017). A theoretical study elucidates that cluster properties strongly depend on the geometry of the isolated cluster and the topology of the cluster sample. F. A. Cotton first used the term cluster for compounds with metal-metal bonds. Nonmetallic atomic clusters came into the limelight later, from both theoretical and experimental studies. The unique size-dependent properties of clusters are distinct from the molecules and bulk solids (Sergeeva et al., 2014). Experimental and theoretical methods on cluster research have seen substantial amount of improvement over the years in discovering a diversity of size-specific phenomena and physicochemical cluster properties (Jena and Castleman, 2010). Various physicochemical properties of metal clusters such as optical, magnetic, thermal, chemical properties differ remarkably from their bulk counterparts (Tsukuda and Hakkinen, 2015).
In the last few decades, several studies have been reported on the ligand protected metal clusters, viz. phosphine-protected small Au cluster, thiolate (RS)-protected Au nanoparticles (Brust et al., 1994), etc. For small metal clusters (<∼100 atoms), the electronic structures are not continuous as in the bulk metals, but rather discretized, which is the primary reason for different physicochemical properties and functionalities in small clusters and bulk metal (Sharma et al., 2017). Studies have shown that clusters made of atoms with appropriate size and composition could potentially mimic the chemistry of elemental atoms in periodic table, and hence are known as superatoms . Various experimental techniques such as laser ablation coupled with mass spectrometry, photoelectron spectroscopy have been employed to get insights into the atomic clusters. Along with experimental studies, theoretical investigations are required to get a better understanding of their geometric arrangement and corresponding properties (Srivastava, 2021).
A large number of theoretical studies have been reported dealing with finding the minimum energy structures of pure elemental clusters by using different optimization algorithms. The potential energy surface (PES) of various atomic clusters is explored to minimize their energy functional with the final objective of locating their global minimum (GM). Different optimization methods such as genetic algorithm (Holland, 1992) (GA), basin hopping (BH) algorithm (Wales and Doye, 1997), particle swarm optimization (PSO) algorithm (Bai, 2010), adaptive particle swarm optimization (APSO) (Zhan and Zhang, 2008), simulated annealing (SA) (Woodley et al., 1999), artificial bee colony optimization (ABC) (Karaboga and Basturk, 2007), honey bee mating optimization (HBMO) (Pham et al., 2005), ant colony optimization (ACO) (Colorni et al., 1991), heuristic algorithm combined with the surface and interior operators (HA-SIO), fast annealing evolutionary algorithm (FAEA), firefly algorithm (FA) (Yang, 2010). etc. are increasingly being used to solve the optimization problem in a time-and costefficient manner. Various models/empirical potentials (EPs) such as Lennard-Jones (LJ), Born-Mayer, Sutton-Chen, Gupta and Murrell-Mottram potentials can effectively explain the bonding within various clusters. A number of studies performed by Chattaraj et al. reveal that PSO is more efficient than commonly used techniques such as GA, SA, and BH for finding the GM of small clusters (Mitikiri et al., 2018;Jana et al., 2019). Further developments over PSO algorithm have been accomplished. Global optimization of boron clusters (B 5 and B 6 ) has been studied using an advanced PSO approach by Mitikiri et al. (2018). Jana et al. (2019) performed a similar study on carbon clusters, C n (n 3-6,10) by using a modified PSO algorithm. Mitra et al. (2020) reported the global optimization of Al 4 2clusters by using firefly algorithm along with DFT. The concept of trapping atoms and small molecules into the hollow cavity of clusters has shown several applications in biology (Cagle et al., 1996;Thrash et al., 1999;Wilson et al., 1999) and electrical engineering (Cioslowski and Nanayakkara, 1992). Pan et al. (2018) reported the encapsulation of noble gas (Ng) atoms into the B 40 host moiety, which is shown to have a fluxional character . The dynamical study of the aforementioned system showed that the fluxional character persists even after the encapsulation. Smaller cages such as C 20 H 20 (Cross et al., 1999;Jiménez-Vázquez et al., 2001), C 10 H 16 (Haaland et al., 2004), BN cages (B 12 N 12 and B 16 N 16 ) (Khatua et al., 2014a), Pb 12 2-, and Sn 12 2- (Sekhar et al., 2017) can act as host molecules to endohedrally trapped noble gas atoms. Cucurbit[n]urils, abbreviated as CB[n]s, n being the number of gycoluril units, can also act as host for different guest molecules including metal cations, organic dyes, drugs, halide ions, etc. (Lagona et al., 2005;Pan et al., 2013a). CB[7] was also reported to bind different guest compounds such as organic dyes (e.g., Stilbenes, naphthalene), viologens, and metal complexes (e.g., Oxaliplatin) applicable in cancer treatment (Wagner et al., 2003). In 2017, Pan et al. (2017) have reported the adsorption of 14 molecules, viz., CH 4 , C 2 H 2 , C 2 H 4 , C 2 H 6 , F 2 , Cl 2 , NO 2 , NO, CO, CO 2 , SO 2 , H 2 S, N 2 , H 2 endohedrally within the hydrophobic inner cavity of CB [7]. CB[6] is also known to encapsulate noble gas atoms . Chakraborty et al. (Chakraborty and Chattaraj, 2015) reported the accommodation of noble gas atoms within the BN-doped (3, 3) single-walled carbon nanotubes. The Ng binding ability of BeX (X SO 4 , CO 3 , O) has been reported by Saha et al. (2015). Pan et al. (2014) have explored the stability of Ng-bound SiH 3 + cluster ions. An emerging host molecule is the basket-shaped octa acid (OA) cavitand that can encapsulate different gas molecules (Chakraborty et al., 2016). Various steroids (Liu and Gibb, 2008), hydrophobic moieties such as ethane, ethylene, acetylene (Stang and Diederich, 2008;Hu et al., 2009;Florea and Nau, 2011;Zhang et al., 2011;Yang et al., 2015;Chong et al., 2016) have been confined inside OA. The encapsulation of gas molecules, especially the hazardous ones, by the cluster cavity has great applications in environmental chemistry. Encapsulation of greenhouse gases (CO 2 ) (Zhang and Chen, 2009;Jin et al., 2010;Kim et al., 2010;Jin et al., 2011;Mastalerz et al., 2011;Lüet al., 2014), air pollutants (NO 2 ) , and poisonous gases (CO, NO) using molecular cages has applications in reducing their negative impact on the atmosphere. N 2 encapsulation (Msayib et al., 2009;Akhtar et al., 2012;Schneider et al., 2012) is yet another important research topic in environmental chemistry. Chakraborty et al. (2016) reported a set of small gaseous molecules (C 2 H 2 , C 2 H 4 , C 2 H 6 , CO, CO 2 , NO 2 , NO, N 2 , H 2 − ), and rare gas atoms as guest molecules for the OA host system. In fact, OA is a very efficient reaction vessel for accommodating various different guest molecules. Li + , Na + , K + , Be 2+ , Mg 2+ , Ca 2+ , Li 3 O + , Na 3 O + , K 3 O + and various nucleobases can occupy the basket-shaped octa acid cavitand as reported by Chakraborty and Chattaraj (2018).
In recent days, the conservation of the atmosphere and the desire to save up fuel for the upcoming generations has been a major concern in the scientific community that led to the search for alternatives of fossil fuel. Hydrogen, being renewable, recyclable, environment friendly, and abundantly available in nature, is now a globally acceptable fuel source with the potential to replace fossil fuels in the near future. The challenge, however, is designing compatible storage and transport materials. To that end, hydrogen-storing capacity of metal-organic-frameworks (MOF) (Rosi et al., 2003;Rowsell and Yaghi, 2005), covalentorganic-frameworks (COF) (Kuc et al., 2007;Cabria et al., 2008), clathrate hydrates (Lee et al., 2005;Chattaraj et al., 2011), polymers (McKeown et al., 2006, carbon nanotubes, BN cages, fullerene, grapheme-like materials (Froudakis, 2001;Deng et al., 2004;Heine et al., 2004;Sun et al., 2005;Wu et al., 2008), metal hydrides (Lee et al., 2005) have been explored. Pan et al. (2012a) performed a theoretical study on the H 2 -storing capability of some Li-doped clusters and super-alkalis. Zhu et al. (2010) have shown cucurbiturils acting as a promising candidate for hydrogen storage. Pan et al. (2013a) have discussed the hydrogen storage capability of the CB[7] system. A different class of compound, alkali-doped carbon materials (graphene sheet and single-walled carbon nanotubes), have been designed for reversible hydrogen storage for transportation purposes by Wei-Qiao Deng et al. (2004). On the other hand, very explosive acetylene can be stored within the porous MOF-505 analogue as reported by Yunxia Hu et al. (2009).
The confinement effect on atoms and molecules has intrigued both theoreticians and experimentalists alike. It brings out interesting changes in the energy levels of the confined systems, their bonding, reactivity, and properties (Grochala et al., 2007;Schettino and Bini, 2007;Sabin and Brandas, 2009;Gubbins et al., 2011;Chakraborty and Chattaraj, 2019;Pal and Chattaraj, 2021). Khatua et al. (2014b) have performed a theoretical investigation on the entrapment of (HF) 2 in C n (n 60, 70, 80, 90) cages. Although CO and N 2 are isoelectric species, the latter is known to be pretty inert owing to its high ionization potential, low electron affinity, and high frontier orbital energy gap (ΔE HOMO-LUMO ). Thus, N 2 capture in various transition metal complexes has proven to induce bond activation that has various industrial applications (Chatt and Leigh, 1972;Yandulov and Schrock, 2003;Latysheva et al., 2012;Bergman et al., 2013;Hoffman et al., 2013). In this regard, Saha et al. (2017) reported the CO and N 2 bound metal supported boron clusters (MB 12 -, M Co, Rh, Ir) which form a spinning umbrella-like structure and activate the bound molecules. Boron clusters have found profound applications in material science owing to their ability to act as nanomaterial building blocks. Their property to act as such is due to a bowl-like structure with an outer rim (B 9 ) and an inner well (B 3 ).
Along with the coordination and inorganic cages, various organic cavitands can also be considered host system for encapsulation. One such class of compounds is the cucurbiturils. Chemical reactions catalyzed by host-guest interactions are comparable to those catalyzed by enzymes (Dong et al., 2012). Hennig et al. (2007) successfully induced protease inhibition using the host-guest interaction with CB [7]. In addition to CB[7], ExBox +4 (Barnes et al., 2013) can also act as an organic host molecule that can encapsulate a wide array of guest moieties.  performed a theoretical study on [4+2] cycloaddition reaction confined within CB[7] and ExBox +4 host systems.
In recent times, low dimensional materials are being given more and more attention to be used as host moiety. Graphene has provided us with a plethora of highly efficient devices such as gas sensors, spintronic devices, nanoelectronics, and optoelectronic devices ). An inorganic counterpart of graphene is the boron nitride doped system that can be functionalized with OLi 4 , CLi 6 , NLi 5 , BLi 7 , Al 12 Be to achieve some interesting properties . The M 3 O + (M Li, Na, K) functionalized graphene nanoflakes (Chakraborty and Chattaraj, 2016a) are known to sequestrate various polar molecules such as CO, NO, and CH 3 OH. Sequestration of gas molecules such as H 2 , O 2 , O 3 , CO, NO, and H 2 O through bare boron nitride flakes (BNF) and metal oxide, MO (M Cu, Ag, Au) functionalized BNF are also reported (Chakraborty and Chattaraj, 2016b). In this review we report some optimization techniques for the generation of minimum energy structures of some selected clusters and also the bonding, reactivity, and different properties of some selected confined systems. Aromatic behavior and electride properties of some clusters are also investigated.

THEORETICAL BACKGROUND AND COMPUTATIONAL DETAILS
Before optimizing the geometry of any system, we carefully ponder over the requirement of the study and select the level of theory maintaining a parity between the level of accuracy required and the computational cost to be incurred. Most often, the easiest way is to take into consideration the experimental data (if available) and select accordingly. The systems discussed in this article are optimized using the computational chemistry software package, Gaussian 09 (Frisch et al., 2009). We have used B3LYP (Lee et al., 1988;Becke, 1992), BP86 (Perdew, 1986a;Perdew, 1986b;Becke, 1988), ωb97X-D (Chai and Head-Gordon, 2008), PBE Perdew et al., 1997), TPSSTPSS Tao et al., 2003), M06, M06-2X (Zhao and Truhlar, 2008), and M05-2X (Zhao et al., 2006) functionals for carrying out DFT calculations of various systems. The exact level of theory (method and basis set) chosen for the individual case studies is mentioned in the Results and Discussion section. Relativistic effects for heavier atoms are taken care of by effective core potentials (ECPs). The stationary states are better understood from the harmonic vibrational frequencies. The minimum energy structure and the transition state (TS) are identified with the presence of zero and one imaginary frequency, respectively.
The atomic charges, nature of interactions present within the systems, and the possible bond formation are analyzed with the help of natural population analysis (NPA) (Reed et al., 1985), Wiberg bond indices (WBI) (Wiberg, 1968) in the NBO scheme (Reed et al., 1988). The electron density topology is mapped using Bader's quantum theory of atoms-in-molecules (QTAIM) (Bader, 1985) in Multiwfn (Lu and Chen, 2012). Parameters such as electron density [ρ(r c )], total electron energy density [H(r c )], local kinetic energy density [G(r c )], local potential energy density [V(r c )], and Laplacian of electron density [∇ 2 ρ(r)] are computed at the bond critical points (BCPs) and they help analyze the extent of covalent or ionic character present along that bond path. The NCI index reveals the localized binding interaction in a system, and the plot can be visualized as red, blue, or green regions in the NCIPLOT program (Contreras-García et al., 2011) depending on whether the interaction is repulsive, H-bond, or van der Waals, respectively. The nonlinear optical (NLO) properties are evaluated in terms of average linear polarizability (α), first (β) and second (c ) hyperpolarizabilities.
ADF 2013.01 software (Baerends et al., 2013) is utilized to perform the energy decomposition analysis (EDA) (Morokuma, 1971) with the natural orbitals for chemical valence (NOCV) (Mitoraj et al., 2009). The interaction between two selected fragments of the studied system is represented in terms of three attractive and one repulsive energy terms. The attractive term includes electrostatic energy (ΔE elstat ), orbital interaction energy (ΔE orb ), dispersion interaction energy (ΔE disp ), while the repulsive term is known as Pauli repulsion energy (ΔE Pauli ).
In NOCV, the orbital term is represented as the sum of ΔE k orb (pairwise orbital energies) which is related to Δρ k (r) (pairwise charge contributions).
Atom-centered density matrix propagation (ADMP) Schlegel et al., 2001;Schlegel et al., 2002) in Gaussian 09, and Born-Oppenheimer molecular dynamics (BOMD) in deMon2K software (Koster et al., 2011) are used to perform the dynamic study of the systems under discussion in this article.

Global Optimization Using Machine Learning Techniques
Minimization of a system's energy functional is the most fundamental step in the determination of its ground state. Reaching the global minimum (GM) geometry, however, poses a number of challenges, the most important being the high probability of getting stuck in local minima in the PES. Swarm intelligence (SI)-based algorithms have turned out to be very effective in searching for optimal solutions in a given search space. Here we discuss three different techniques, PSO combined with DFT, PSO with CNN, and DFT-integrated FA. They do not need to implement any symmetry constraint, or consider bond characterization. DFT-PSO adjusts each particle's trajectory at every time stamp while following the convergence criteria. We have successfully implemented these techniques to find the GM configurations for small-sized nonmetallic clusters such as Boron (B 5 and B 6 ) (Yuan et al., 2014), Carbon (C 5 ) (Jana et al., 2019), and polynitrogen clusters (N 4 2and N 6 4-) (Mitra et al., 2021), and metallic clusters such as Al 4 2- (Mitra et al., 2020), Au n (n 2-8), and Au n Ag m (2 ≤ n+m ≤ 8) (Mitra et al., 2021).

PSO Combined With DFT (DFT-PSO)
We start off with 14 and 15 random structures for B 5 and B 6 ( Figure 1), respectively, with velocity set at zero and maximum number of iterations set at 1,000. No significant change in bond length is detected after reaching the global best configuration at the end of the PSO run. Following this, the optimization is performed at the B3LYP/6-311+G(d,p) level as the postprocessing step that helps align the symmetry of the PSOobtained final structure and obtain the corresponding exact energy. The post-processing in this case takes only 20 s to complete, and the final geometry obtained is energetically very close (0.0015 eV difference) to that obtained at the end of the PSO run. The zero-point energy (ZPE) corrected energy, free energy, and enthalpy for B 5 (C 2v ) are −123.9873, −124.0135, and −123.9821 a.u., respectively, while those for B 6 (C 2h ) are −148.8100, −148.8381, and −148.8038 a.u., respectively. A comparison drawn between our method and other popular algorithms such as DFT-SA and DFT-BH reveals that while these two require a CPU time of 369.64 and 455.43 min to locate the minimum energy structure of B 5 , respectively, our method takes only 80.50 min. It is also observed that while the BH and the SA require 600 (unconverged) and 324 number of iterations, respectively, our modified PSO converges after only 138 number of iterations.
For the carbon clusters, 10 random configurations are chosen with initial velocity zero and 1,000 number of iterations. For the C n (n 3-6) clusters, linear geometries are obtained with D ∞h point group as the GM. For n 4-6, a cyclic isomer for each of them is also obtained with point groups D 2h , C 2v , and D 3h for C 4 , C 5 , and C 6 , respectively, whereas for the relatively larger C 10 cluster, the GM geometry is a D 10h ring structure. It is to be noted that the geometries and corresponding energies reported here match with those obtained from the experimental reports (Raghavachari and Binkley, 1987;Watts et al., 1992;Hutter and Lüthi, 1994;Pless et al., 1994;Martin and Taylor, 1996;Van Orden and Saykally, 1998). Again, comparing with DFT-SA and DFT-BH, we get encouraging results for our modified PSO approach. The total execution time for our technique is 143.30 min versus 215.98 and 5085.67 min for DFT-SA and DFT-BH, respectively.
For N 4 2and N 6 4clusters, convergence takes place after 483 and 627 numbers of iterations with 228 and 323 min of execution time, respectively. The geometries obtained before and after the post-processing step are energetically very close ( Figure 2). The N 6 4cluster, however, shows a higher order saddle point at the post-processing symmetry constrained optimization with point group D 6h . In case of the binary gold-silver clusters (Au n Ag m ), AuAg 2 has a ring (C s ) doublet GM, those with n+m 5 have trapezoidal GM, n+m 6, 7 have triangular 3D geometry (C 1 ), and n+m 8 failed to converge within the initial coordinates range of [−4, 4].

CNN With PSO (CNN-PSO)
Supervised learning (CNN) is performed on an initial guess set generated with the help of (ADMP) simulation. Their corresponding single point energies (SPEs) are calculated and stored to be read by the PSO to search for the GM geometry. Statistically relevant analysis is derived by making the method 8fold (each with 30 files containing number of iterations and the SPEs). Remarkably high success rate (∼77-90%) is observed for this combined ADMP, CNN, PSO technique, indicating its efficiency in finding GMs. This technique is tested with the C 5 cluster and we have obtained the previously reported (Van Orden and Saykally, 1998) linear geometry as the GM with a higher convergence rate. However, two local minima are also detected due to a premature convergence (Figure 3).  Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 5

FA With DFT (DFT-FA)
A comparative study of the DFT-integrated FA algorithm with PSO is performed on Al 4 2cluster considering planar and nonplanar structures, and it turns out that the former performs better than the PSO. The mean convergence times for the PSO and FA are 68.22 and 61.25 min for the planar, and 85.13 and 74.40 min for the nonplanar approach, respectively. The corresponding success rates are also higher for the FA. Since we know that the GM of Al 4 2is planar, we have also investigated the search space of only the planar geometry to get a faster convergence since the number of variables decreases in this problem. Again, the modified FA performs better. A relation is also drawn between the stabilization energy of the system and the change in its aromaticity and number of iteration steps it takes to converge to the GM. The energy functional optimization and the NICS (0) value of Al 4 2− are scanned and depicted in Figure 4. It is observed that the aromaticity increases with the decrease in the energy, i.e., with the increase in the stability of the system.

Aromaticity of Clusters From a CDFT Perspective
It is well known from Hückel's (4n+2) π electron theory (Hückel, 1931) and Pauling's quantum mechanical description (Pauling and Sherman, 1933) that aromaticity of conjugated systems with cyclic and planar geometry is associated with their increased stability. Various structural, energetic, electronic, and magnetic behavior-based parameters are considered to analyze the aromaticity of the systems, among which the Nucleus-Independent Chemical Shift (NICS) (Schleyer et al., 1996) is, perhaps, the most widely used criterion for aromaticity. CDFT also plays an important role in aromaticity determination (Chattaraj et al., 2005;Chattaraj et al., 2006;Chattaraj et al., 2007), and it does so with the help of reactivity descriptors and associated electronic structure principles (Chattaraj et al., 1993;Chattaraj and Maiti, 2001;Pan et al., 2013b;Chakraborty and Chattaraj, 2021). The relative aromaticity index ΔX, where X could be energy (E), polarizability (α), electrophilicity (ω), or hardness (η), is defined as the difference between the respective indices in the cyclic and open (or localized) systems, i.e., ΔX X CYCLIC -X OPEN(/ LOCALIZED) . They are observed to show a similar performance as those of the NICS and MCI values. They also provide valuable insights into the stability, reactivity, and electronic properties of the associated cluster. All-metal aromatic Al 4 2− (Li et al., 2001) and antiaromatic Al 4 4− (Kuznetsov et al., 2003) clusters are investigated at the B3LYP/6-311G(d, p) level of theory to analyze their aromaticity from a CDFT perspective (Chattaraj et al., 2005;Chattaraj et al., 2006;Chattaraj et al., 2007). While Δη > 0 represents aromaticity and Δη < 0 represents antiaromaticity, the reverse is true in the cases of E, α, and ω indices. Their values for the aluminum clusters are compared with those of benzene (C 6 H 6 ) and cyclobutadiene (C 4 H 4 ). While Al 4 2− and C 6 H 6 show positive Δη and negative ΔE, Δα, and Δω values indicating an aromatic behavior, C 4 H 4 shows an exact opposite trend to account for its antiaromaticity. For the Al 4 4− cluster, however, we have obtained somewhat contradictory results. ΔE and Δα values indicate the cluster's antiaromatic nature, whereas Δη, Δω, and NICS values reflect its aromatic nature. It is in conformity with the current knowledge that this cluster exhibits conflicting aromaticity. Experimental synthesis and theoretical studies reported along with ELF analysis (Li et al., 2001;Santos et al., 2005) suggest the cluster to be antiaromatic. Its σ-aromaticity directs the overall aromaticity by dominating over its π-counterpart as studied through NICS (Chen et al., 2003) and magnetic field induced current density (Havenith et al., 2004) analyses. Such conflicting aromaticity, along with other varieties of multiple aromaticity and antiaromaticity, δand Φ-aromaticity, bond stretch isomerism, etc. are exhibited by several other allmetal clusters (Zubarev et al., 2009;Zubarev and Boldyrev, 2011). Be 3 2− , Ca 3 2− , and Mg 3 2− clusters are classified as aromatic in terms of the ΔX indices Giri et al., 2010). Other applications of aromatic clusters are studied by our group through molecular electronic transport (Khatua et al., 2008), hydrogen storage (Havenith et al., 2005;Giri et al., 2011a;Giri et al., 2011b;Das and Chattaraj, 2012;Srinivasu et al., 2012;Pan et al., 2012b), and Zn-Zn and Be-Be bond stabilization Roy and Chattaraj, 2008) in the domain of CDFT. The ΔX and NICS parameters show their versatility in quantifying the aromaticity of not just planar and cyclic systems, but also any other nonplanar closed structure. They are both easily computable and ΔX also has a conceptual lucidity since it originates from the electronic structure principles of CDFT.

Structure, Bonding, and Reactivity of Various Molecular Electrides
Certain chemical entities contain loosely bound electrons not directly connected to any atom(s) within the cluster, but trapped Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 6 within a hollow space (cavity of cage compounds or packing void in crystals) that act as anions. Such entities are known as electrides and they are known for their nonlinear optical (NLO) properties. Besides showing NLO properties, which is considered to be an identifiable character of an electride, it is also widely applicable in electron emission, catalysis, reversible hydrogen storage, super conductivity, etc. (Toda et al., 2007;Xu et al., 2007;Kitano et al., 2012). Organic electrides of crown ethers and cryptands (Ellaboudy et al., 1983;Ward et al., 1988;Dawes et al., 1991;Xie et al., 2000), and inorganic electrides such as [Ca 24 Al 28 O 64 ] 4+ (4e − ) (Matsuishi et al., 2003), Y 5 Si 3 (Lu et al., 2016), [Ba 2 N 2 ](e − ), [Li 2 (Lu et al., 2018) are well reported in the literature. A different class of electrides, known as molecular electrides, are basically guest@host complexes containing a significant amount of localized electron cloud within the void of the host. Cavity-containing molecular structures such as decaborane (Muhammad et al., 2009;Muhammad et al., 2011), pyrrole (Chen et al., 2005), tetracyanoquinodimethane (TCNQ) (Li et al., 2009), fullerene cages , C 20 F 20 (Wang et al., 2012), C 60 F 60 , extended (3.1.3.1) porphyrin (EP) (Saha and Chattaraj, 2018), and many more are utilized for this purpose. The guest atoms are usually alkali and alkaline earth metals. For a guest@host complex to be characterized as a molecular or cluster electride, certain criteria need to be fulfilled such as the presence of a non-nuclear attractor/maximum (NNA/NNM) (a non-nuclear critical point with a local maximum of electron density), a negative Laplacian of electron density [∇ 2 ρ(r c )], presence of an ELF basin near the Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 NNM, high NLO properties, and a green region in the NCI plot showing accumulation of electron density. The Mg 2 EP complex studied at the M06-2X-D3/6-311G(d,p) level of theory (Saha and Chattaraj, 2018) shows the presence of NNA in between the two Mg atoms where the value of ∇ 2 ρ(r c ) < 0, with an ELF basin nearby. The electron population at said NNA is 1.02 e with 46% localization. The NLO properties in terms of α, β, and c are calculated and compared with other electride systems that show that α is higher while β and c are lower in the studied systems. The donation of the loosely trapped electron to antibonding MOs of certain bonds in small molecules (H-H in H 2 , C-O in CO 2 , N-O in N 2 O, and C-H in CH 4 and C 6 H 6 ) results in the activation followed by the dissociation of the respective bonds (Saha and Chattaraj, 2018). Another study (Saha et al., 2019) performed by our group utilizes the modified form of β-diketiminate ligand ( Dipp Nacnac) to hold four Mg atoms with two equivalent Mg(I)-Mg(I) bonds ([Mg 4 ( Dipp H L) 2 ] 2-) in its lower energy singlet state. Two NNAs are found to be present at the center of each Mg(I)-Mg(I) bond with an electron population of 1.18 |e| and 52% localization. The calculated values of α, β, and c are 891.7, 0.0, and 8.9 × 10 5 a.u., respectively, which clearly indicate the system to be classified as an electride. This system is further stabilized by sandwiching it between two K@crown-6ether + [K@CE] + counter cations. Again, the C 60 cage trapping a magnesium dimer (Mg 2 @C 60 ) ) and a lithium trimer (Li 3 @C 60 ) (Das and Chattaraj, 2021a) act as electrides as indicated by the presence of NNA at the center of the Mg-Mg bond path and Li 3 cluster in the respective complexes. The Li 3  Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 9 Chattaraj, 2020) contain NNAs at the center of the M-M bonds with population varying between 0.95 and 1.39 and percentage localization ranging within 43-62%. Discernible substituent effects on electride characterizers have also been reported (Das and Chattaraj, 2021b).

Noble Gas Encapsulated B 40 Cage
The encapsulation of noble gas (Ng) atoms in B 40 cavitand along with their structure and interactions of Ng with the Ng and B atoms in the host-guest complex is discussed with the help of DFT-based computations (Pan et al., 2018). Dissociation energy (ΔE diss ) and Gibbs free energy change (ΔG diss ) are calculated to study the stability of the encapsulated complexes. NBO, EDA, and NOCV calculations have been done for studying the nature of bonding. The optimized structure of the B 40 cage along with the Ng@B 40 and Ng 2 @B 40 systems (at ωB97X-D/def2-TZVP level) are depicted in Figure 5. The size of the B 40 cage is suitable to accommodate He and Ne atoms at its center, whereas for heavier atom (Ar-Rn) encapsulation, the cavity diameter expands with the help of a certain amount of energy (preparation energy, ΔE prep ). Although these complexes are thermochemically unstable, they remain in their encapsulated form on account of their high kinetic barrier. Despite that, the lighter Ng encapsulated complexes have very low ΔE diss (−1.8 kcal mol −1 for He and −7.1 kcal mol −1 for Ne) as compared with the experimentally identified He@C 20 H 20 complex (−33.8 kcal mol −1 ). For the heavier Ng atom-encapsulated Ng@ B 40 systems (Ng Kr-Rn), the ΔE diss increases with the size of Ng. The possibility of releasing Ng atoms in the dissociation process is either through B 7 or B 6 holes for the lighter He-Ar atoms, whereas the heavier ones can only escape through the B 7 holes due to their larger size. Decapsulation through the B 7 hole has ΔG ≠ values ranging within 84.7-206.3 kcal mol −1 . The rate constant (k) calculated at 298 K for the dissociation through either the B 7 or the B 6 hole comes out to be pretty low suggesting that all the Ng@B 40 systems are kinetically stable.
In the case of two Ng atoms encapsulation, the inter atomic distance decreases from that of their free state. For Ar 2 @B 40 and Kr 2 @B 40 , larger repulsion results in the exergonic dissociation of Ng 2 @B 40 into Ng and Ng@B 40 . Along with this, large ΔE prep indicates the nonviability of Ar 2 B 40 and Kr 2 B 40 . For even heavier Ng atoms (Xe and Rn), no minimum energy structures are obtained for the corresponding dimer encapsulation. The corresponding transition states of the Ng release process of Ng 2 @B 40 (Ng He-Ar) suggest that one Ng atom approaches and leaves through the B 7 rings, while the other remains near the center of the cavity. The associated ΔG ≠ values for the He and Ne dimer encapsulated complexes are high enough for them to be kinetically viable.
The B 40 cage exhibits a fluxional behavior due to the continuous interconversion between the B 6 and B 7 rings caused by the transfer of one B center from B 7 to B 6 . The encapsulated system, Ng@B 40 , also shows similar dynamic behavior as that of the bare B 40 . The Ng atom inside the cage does not have any substantial influence on the fluxionality of the cage which is reflected in the free energy barrier values (16.4 kcal mol −1 for the bare cage, and a range of 16-18.

Small Gas Molecule Encapsulation Within Octa Acid Cavitand
Small gas molecules such as C 2 H 2 , C 2 H 4 , C 2 H 6 , CO, CO 2 , NO 2 , NO, N 2 , H 2 , and Ng atoms (He n -Xe n , n 1,2) are selected as guest molecules encapsulated in OA cavitand ( Figure 6) and analyzed via DFT approach (Chakraborty et al., 2016). The systems under study are optimized at the ωB97X-D/6-311G(d,p) level of theory (LanL2DZ basis set with ECP is used for Xe). There are two possible cavities for the accommodation of the guest atoms, the inner cavity of OA is more suitable as it increases the host-guest interaction. In the case of Ng@OA, due to encapsulation of Ng atoms, no notable distortion is observed in the OA. For Ng 2 @OA systems, one guest atom can occupy the center of the host cavity, whereas the other one remains in the outer cavity. The Ng-Ng bond distances inside OA are 3.6, 3.7, 3.8, 3.9, and 4.1 Ǻ, respectively, for He 2 , Ne 2 , Ar 2 , Kr 2 , and Xe 2 . In the cases of CO@OA and NO@OA, the guest molecules prefer to stay in the inner cavity of OA. Their orientation with respect to the two nearest benzene-like fragments is almost perpendicular, whereas it is parallel to the rest of the benzene fragments of OA. The presence of these π electron clouds close to the encapsulated guests is expected to have a significant impact on the stability of the complexes, the nature of interaction, and dynamical behavior as well. For N 2 @OA and H 2 @OA, N 2 and H 2 remain well inside the inner cavity. Thermochemical study reveals that except for He@OA, in all cases, D 0 value is positive indicating the stability of the host-guest complexes concerning their dissociation into the corresponding individual components. Going from lighter to heavier Ng atoms (also for Ng 2 @OA), D 0 value increases, i.e., the host-guest interaction increases. This could most likely be due to the increasing polarizability of the Ng atoms down the group. The dissociation channels for all the encapsulated complexes have positive D 0 values. Most of them, however, dissociate spontaneously at room temperature except in the cases of Kr, Kr 2 , Xe, C 2 H 2 , C 2 H 4 , C 2 H 6 , and N 2 guest molecules. Thus, the encapsulation of the mentioned guest molecules is favorable at 298 K. In the cases of NO/NO 2 @OA and CO/CO 2 @OA, an increase in D 0 value is observed from NO to NO 2 and from CO to CO 2 . It can thus be deduced that the encapsulation of NO 2 and CO 2 inside OA forms more stable complexes compared with NO and CO, respectively. Having said that, it is to be noted that all the four complexes have favorable dissociation channels at room temperature. The hydrocarbons have better interaction with the OA and hence are not prone to dissociation at ambient temperatures. The interaction between the host and guest moieties is analyzed with the help of NBO, NCI, and EDA. All the Ng atoms acquire some positive charges, i.e., transfer of electron density occurs from the Ng atoms to the host OA surface. The donation primarily occurs from the lone pair (LP) of Ng to the C-H antibonding orbital of OA for all the Ng atoms [except Ne where it takes place from LP to antibonding Rydberg state (Ry*)] and for CO, NO, N 2 , and CO 2 molecules. WBI values for Ng-Ng and Ng-OA interactions indicate a purely noncovalent character therein. From this discussion it is clear that the confinement brings about an increase in the reactivity of all the guest atoms/ molecules within the OA. NCI isosurfaces show the presence of green surface around the guest molecules indicating van der Waals interaction that stabilizes the host-guest complexes. EDA results show that there exists closed shell type of interaction between guest Ng/H 2 and OA. The ΔE disp and ΔE orb are the largest and the smallest contributors toward the total attractive interaction, respectively. The former increases in magnitude with increasing size of the encapsulated Ng atoms. For the guest hydrocarbons, the contribution from ΔE disp increases and ΔE elstat decreases with the increasing number of H atoms. This is because the molecules such as C 2 H 2 and C 2 H 4 containing labile electron cloud can accumulate enough positive charge to favorably interact with the electron-rich fragments of the OA. This makes the contribution from ΔE elstat very important in stabilizing these complexes. The contribution from charge transfer and polarization are very less toward ΔE tot as indicated by the very low values of ΔE orb . A similar type of situation is observed for CO/CO 2 @OA complexes. For nitrogencontaining guests, the ΔE orb outweighs the ΔE elstat contribution toward ΔE tot . Since the main stabilizing factor for all the complexes is ΔE disp , the nitrogen-containing guest molecules are prone to be affected by the polarization or charge transfer by OA, in comparison with the other guest molecules. ADMP simulation performed at 298 K shows that all the Ng atoms (except Ne) are prone to leaving the cavity, whereas at 50 K they remain within the OA. Polar molecules such as CO, CO 2 , NO, and NO 2 have a higher tendency to stay within the OA than the nonpolar H 2 . Again, C 2 H 2 , C 2 H 4 , and N 2 containing π electron cloud also prefer to stay inside OA. From the above thermochemical, kinetic, and dynamical analyses, it can be said that OA makes a reasonably good choice for accommodating gas molecules.

Cucurbit[n]urils as a Host Moiety
Cucurbiturils are methylene-linked macrocyclic molecules having glycoluril unit [ C 4 H 2 N 4 O 2 ] as a building block. This repeating glycoluril unit can bind with hydrogen with sufficient amount of binding energy. Thus CB[n] can be designed as an effective hydrogen storage compound. The nitrogen and the oxygen centers are found to be the most active centers to bind with hydrogen with positive binding energy. It was found that (CH 3 ) 2 C 4 H 2 N 4 O 2 (CH 3 ) 2 unit can interact with total 13H 2 atoms (Pan et al., 2013a). Since hydrogen has an electric quadrupole moment, a charge-quadrupole interaction plays a pivotal role in binding the hydrogen with the host. NPA charge analysis reveals that the charge transfer occurs from the N and O centers to the σ* orbital of H 2 molecule.
Among the CB[n] family, CB[7] can accommodate five H 2 molecules endohedrally. The O centers can adsorb a total of 28 H 2 molecules (two per O atom) and 19 H 2 molecules get adsorbed at the N centers exohedrally, making it a total of 52 hydrogen molecules. The binding energy and adsorption enthalpy are positive and negative, respectively, indicating CB[7] to be a potentially promising H 2 storage material. The gravimetric wt % of hydrogen for the CB[7] adsorbing 52 H 2 molecules comes out to be 8.3 with an average binding energy of 7.8 kJ mol −1 . These values are very encouraging when compared with various other potential hydrogen storage materials such as α-cyclodextrin (Zhu et al., 2010), COFs (Klontzas et al., 2008;Li et al., 2010), MOFs (Bhatia and Myers, 2006;Vitillo et al., 2008), Li-doped nanotubes , and polyacetylenes (Li and Jena, 2008), with average binding energies ranging within 5-8 kJ mol −1 .
The endohedral adsorption of gas molecules such as C 2 H 2 , C 2 H 4 , C 2 H 6 , CH 4 , CO, CO 2 , NO 2 , NO, N 2 , H 2 , F 2 , Cl 2 , H 2 S, and SO 2 into the CB[7] cavitand is depicted in Figure 7. Geometry optimizations are performed at the ωB97X-D/6-31G(d,p) and ωB97XD/6-311+G(d,p) levels . Both the CB[7] cage and the encapsulated gas molecules remain unaffected by the encapsulation. For SO 2 encapsulation, the binding enthalpy shows the highest value (14.3 kcal mol −1 ) followed by Cl 2 and C 2 H 2 (11.6 and 10.4 kcal mol −1 , respectively). CB[7] also encapsulates C 2 H 4 and C 2 H 6 more favorably than CO 2 , NO 2 , and H 2 S. The binding enthalpies for NO/F 2 /N 2 /CO/CH 4 @CB[7] systems vary from 4.7 to 5.8 kcal mol −1 , the highest value corresponding to CH 4 . From the enthalpy values it is clear that CB[7] can selectively adsorb SO 2 among various gas molecules, and hence can be applicable in the SO 2 separation process from gas mixtures. ΔG value suggests that C 2 H 6 is less prone to be encapsulated than CO 2 inside CB[7]. SO 2 , Cl 2 , and C 2 H 2 adsorb with negative ΔG values, whereas those of C 2 H 4 and CO 2 are slightly endergonic (0.6-0.7 kcal mol −1 ). The corresponding ΔG values for the adsorption of C 2 H 6 , N 2 , F 2 , NO 2 , NO, and H 2 S vary Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 within 1.3-2.8 kcal mol −1 at 298 K temperature. EDA results show that for all the discussed complexes, ΔE disp contributes more toward the stabilization of the host-guest systems. ΔE elstat term also plays an important role here. In the hydrocarbons, as the number of H atoms increases, ΔE elstat decreases gradually due to the reduction in the acidic character of the H atoms. Higher ΔE elstat and lower ΔE pauli values in the SO 2 encapsulated complex make the interaction between the host and the guest stronger than for C 2 H 4 and C 2 H 6 analogues. This again validates the higher SO 2 selectivity of CB [7]. Cucurbit[6]uril, among the CB[n] family, is found to be a compatible host for the encapsulation of Ngs within its cavity . The optimized geometries at the ωB97X-D/6-311G(2d,p) level of theory of the Ng n @CB[6] complexes are provided in Figure 8. CB[6] effectively accommodates three Ne atoms, but can only trap two of the large Ar and Kr atoms. No significant distortion in the cage is observed for trapping all the three Ne atoms or for the first atom of Ar and Kr, whereas inserting a second atom deforms the shape of the host. The Ng dissociation process becomes more endothermic as we move from Ne to Kr. At 298 K, the dissociations of all Ng n @CB[6] are exergonic except Kr@ CB [6]. At 77 K, apart from the second Ng (Ar and Kr) atom dissociation from Ng 2 encapsulated CB[6], all dissociations become endergonic. Kr encapsulation at 298 K and 1 atm pressure is thermochemically favorable, whereas for Ne and Ar encapsulation, high pressure and moderately low temperature are preferred.
NPA charge analysis reveals N and O to be negatively charged in Ng n @CB[6], while C and H have positive charges. Slight charge transfer (∼0.01 e -) occurs from Ng→CB[6] moiety. Small ρ(r c ) value and positive ∇ 2 ρ(r c ) and H(r c ) values from topological analysis suggest the interaction to be of closed shell type. ELF analysis shows an absence of electron localization between the Ng-Ng and Ng-cage atoms, corroborating the result obtained from AIM. EDA analysis reveals the contribution from ΔE disp to be the largest, followed by ΔE elstat , and the smallest contribution is from ΔE orb , all of which gradually increases going from Ne to Kr. The green surfaces observed between the Ng and CB[6] units in the NCI isosurface ( Figure 9) are an indication of a small van der Waals interaction, which increases with the size of the Ng atoms. The dynamical study (ab initio MD) for 1 ps and at 298 K reveals that Ne and Ar remain inside the cavity, whereas Kr and all the Ng 2 in Ng 2 @CB[6] move toward the open end but do not leave the cage. At 77 K, all guests stay inside the host.

Small Molecules Bound Metal Coordinated Boron Cluster
The activation of small molecules by a metal-supported boron cluster is studied through DFT calculations (Saha et al., 2017). They are known to be applicable in nanomaterial building blocks, automobiles (Norbye, 1971;Jiménez-Halla et al., 2010) etc. The global minimum energy structures calculated at PBE/def2-TZVPPD Weigend and Ahlrichs, 2005) level of MB 12 -, CO@MB 12 -, N 2 @MB 12 clusters and their corresponding TSs for the internal rotation of the B 3 ring are provided in Figure 10. The coordination of the small molecules with the MB 12 cluster forms an umbrella-shaped structure in which the M-L bonds act like the stick of the umbrella. The coordination of CO with the metal center can take place through both the C and O ends, the former producing a more stable isomer. For both OCMB 12 _ and NNMB 12 _ , the Ir-L bond has the highest strength, followed by Co and Rh, while for a particular M center, CO forms a stronger bond than N 2 . ΔG values of these complexes are highly positive which suggest that the corresponding complexes are thermodynamically stable concerning the dissociation process. The O-side bound isomers, however, have low positive ΔG for Co and Ir, and become slightly negative for Rh. They can be made viable by lowering the temperature. The N−N and C-O bonds get lengthened due to complexation in the order COMB 12 -< N 2 MB 12 -< OCMB 12 causing a red shift in their bond stretching frequencies which is the highest for the Ir analogues.
NBO analysis reveals that upon complexation, the metal centers get more negatively charged, apart from COIrB 12 _ and NNIrB 12 _ , where Ir still contains positive charge (less than that in IrB 12 − ). L→M and M→L back transfers take place, and in certain complexes the latter completely compensates (or overcompensates) the former which is indicated by the zero (or negative) charge on the ligand. The Wiberg bond indices suggest that the covalent character follows the order M-C in OCMB 12 -> M-N in N 2 MB 12 -> M-C in COMB 12 -. From EDA-NOCV, it is seen that the bonding between the metal and the ligand is predominantly orbital and electrostatic interactions , however, the contribution from ΔE elstat is higher than ΔE orb . Figure 11 shows the deformation densities [Δρ(r)] for the pairwise orbital interactions for the LMB 12 _ complexes, where a shift in the electron density occurs from the red to the blue region. The Δρ(σ 1 ) plot reveals that the shift of electron density occurs through the L→M→B scheme. Δρ(π 1 ) and Δρ(π 2 ) explain the π electron density shift from the d L→M . The extent of L←M π-back-donation is greater than the L→M σ-donation, causing a red-shift in its stretching frequency. The σ-donation occurs from the HOMO (CO) to the LUMO (LMB − 12 ) fragment and π-back-donations occur from the degenerate HOMO (MB − 12 ) to the degenerate π* LUMO (CO) . An internal rotation of the inner B 3 ring with respect to the outer B 9 ring occurs within the MB 12 cluster. The energy barrier associated with this rotation is reported Liu et al., 2016) to follow the order Co > Rh > Ir. BOMD simulation at 800 K shows the L-M bonds to be intact during the rotation. This makes the complex seem like a spinning umbrella with the L-M bond as the stick.

Hydrogen Storage in Clathrate Hydrates, Li-Doped Clusters, and Super Alkalis
Clathrate hydrates, a class of inclusion compounds, are known to encapsulate guest compounds within its hydrogen bonded polyhedral cage (Mao et al., 2002;   Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 and Mao, 2004;Lee et al., 2010). They constitute a very effective host for hydrogen storage. Four types of clathrate hydrates and their maximum possible hydrogenencapsulated complexes studied at B3LYP/6-31G(d) level (Chattaraj et al., 2011) are depicted in Figure 12. 5 12 represents the cavity having 12 pentagonal faces, whereas 5 12 6 k (k 2,4,8) represents 12 pentagonal faces along with k hexagonal faces. Here we discuss the structure, bonding, and stability of the bare and hydrogen-encapsulated complexes from a density functional theory perspective. For nH 2 @5 12 complexes it is seen that for the first H 2 encapsulation, the process is energetically favorable although the overall nH 2 encapsulation is method dependent. Owing to the small size of 5 12 cavity, it can accommodate a maximum of five H 2 molecules, after which a deformation in the cavity is observed. The GM for H 2 confinement in the 5 12 cavity occurs endohedrally. The H 2 molecules favor the inside of 5 12 more than the outside. In the case of 5 12 6 2 cage, it can also take up a maximum of five H 2 molecules. Slight distortion is observed in the system that becomes more noticeable during the third H 2 encapsulation which slowly decreases for the fourth and the fifth hydrogen molecule encapsulation. This is reflected in the slightly conflicting trend in the corresponding interaction energies. Now in the case of 5 12 6 4 clathrate, obtaining the minimum energy structure was difficult. It is fascinating to note that the encapsulation of one H 2 into the cage stabilizes the structure although it could not provide with the minimum energy structure. Further incorporation of guest molecules deforms the structure of the system. For 5 12 6 8 , the interaction energy for all the six H 2 encapsulation is negative making the process favorable. The large size of the host cavity makes it feasible to accommodate all the six guest molecules efficiently. Positive ΔG value suggests that the complexes are kinetically stable. Finally, it can be concluded that the 5 12 and 5 12 6 2 clathrates can encapsulate up to two hydrogen molecules without undergoing any structural distortions, whereas the 5 12 6 8 clathrate may entrap up to six H 2 molecules depending upon the level of theory used. Calculation of CDFT-based reactivity descriptors of the complexes with and without H 2 encapsulation suggests that for most of the systems, stability increases with the increase in number of trapped hydrogen molecules. This is concluded from the increasing hardness and decreasing electrophilicity values. Li ion is popularly known to bind well with hydrogen molecule owing to its positive charge (Pan et al., 2012a). Inspired by this, a number of efforts have been made to effectively polarize the Li center of various clusters to increase its hydrogen adsorbing ability. Here we study the H 2 storage potential of the Li-doped clusters, M 5 Li 7 + (M C, Si, Ge), M 4 Li 4 (M Si, Ge) at the M06/6-311+G(d,p) level, and some super-alkali ions at the M052X/6-311+G(d) level (Figures 13, 14). The Li centers attain a net positive charge due to the high polarizability of the clusters, facilitating electrostatic interactions to bind with the H 2 molecules. The negative values of interaction energies and enthalpies

(HF) 2 Confinement in Fullerene Cages
The influence of encapsulation on the hydrogen bond strength in (HF) 2 within the fullerene cages is studied using DFT and ab initio MD (Khatua et al., 2014b). The optimized geometries of (HF) 2 @C n (n 60, 70, 80, 90) complexes wB97X-D/6-31G are depicted in Figure 15. The dissociation energy, enthalpy, and change in free energy are negative for the (HF) 2 @C 60 system which indicates that the encapsulation process is thermodynamically unfavorable, whereas positive values for the rest of the HF encapsulated C n cages imply them to be favorable (highest being for the C 80 cage). Owing to the smaller size of the C 60 , the HF units orient themselves antiparallelly to reduce repulsion at the cost of hydrogen bond strength. Thus, the energy associated with the HF-HF interaction is observed to be highest in the C 60 cage (positive, and hence repulsive in nature). For all the studied cases, upon encapsulation, the hydrogen bond distance reduces from that in the free state, the least being inside the C 70 cage. The EDA study reveals that the contribution from ΔE pauli increases and the ΔE int value decreases with decreasing the C n cage cavity except for C 80 cage. For the C 60 cage, a very large value of ΔE pauli makes the overall ΔE int value positive. On account of the smaller H-bond distance within the C 70 and C 90 cages compared to the same within the C 80 cavity, Frontiers in Chemistry | www.frontiersin.org August 2021 | Volume 9 | Article 730548 both the ΔE elstat and ΔE orb contribute more to the attractive interaction than those in C 80 . AIM analysis reveals that for all these confined systems, ∇ 2 ρ(r c ) > 0 and H(r c ) < 0 implying the partial covalent nature of the hydrogen bonds. The hydrogen bond is mostly covalent in case of (HF) 2 @C 70 ELF analysis that also supports this observation.

CONCLUDING REMARKS
There exists an appreciable amount of interest in the field of cluster chemistry, especially in the gas-phase and surfaceadsorbed studies, and for good reasons. Common curiosities in this area include the difference between the properties exhibited by the bulk and individual clusters, how the cluster size affects the overall behavior of the bulk, etc. Other important branches of this cluster chemistry include solving their global optimization problem in a fast and cost-effective way, and investigating the effect of confinement on the cluster-encapsulated systems. The global optimizers discussed in this review are shown to locate the global minimum configurations for small metallic and nonmetallic clusters with less execution time and higher success rate than commonly used optimization algorithms, without having the need to impose any symmetry constraint or any other external restrictions. The only requirement is to adjust the local and global best parameters at each iteration. Comparisons made between our modified PSO with other DFT-integrated BH and SA reveal the superiority of the former with respect to the total execution time and number of iterations the program takes to converge. Again, the DFTintegrated FA turns out to be more efficient than the modified PSO. Furthermore, the ADMP-CNN-PSO technique is well suited for locating the global solution from a huge dataset of initial configurations.
The effect of adsorption and confinement of hydrogen, noble gas atoms, and various other small molecules on their stability, reactivity, nature of interactions, and dynamics are studied from a DFT perspective. The concept of aromaticity is analyzed in terms of CDFT-based descriptors such as E, α, ω, and η, where a lower value of the first three parameters and a higher value of hardness in comparison with that of a reference system characterize an aromatic molecule. The reverse is true for antiaromatic compounds. Certain guest@host complexes containing loosely bound electrons acting as anions and showing high NLO properties, known as molecular electrides, are capable of bond activation in small molecules. Other host-guest complexes exhibit fluxionality. One such example is the B 40 cage whose fluxional property remains unaltered even after Ng atoms encapsulation. The complexation ability of the B 40 cage is also studied in some sandwich complexes and it is seen than the presence of Xe within the cage enhances its complexation ability. The gas molecules accommodated within the Octa acid cavitand become slightly more reactive compared to their free state. Most of the OA-guest complexes are stable with respect to dissociation. OA can thus be designated as a reasonably good storage material for a variety of small gas molecules. Cucurbiturils form another class of compounds which is well known for its hosting capabilities. CB [6] can act as an efficient noble gas carrier and CB [7] can bind up to 52 hydrogen molecules (8.3 wt%). CB[7] is also found to be highly selective toward the adsorption of SO 2 and hence can be used in separating SO 2 from a gas mixture. It is also known to accelerate the otherwise slow [4+2] cycloaddition reaction. The binding ability of the transition metal boron cluster (MB 12 -) with isoelectronic species, CO and N 2 , is studied along with its fluxionality. Bond activation in both CO and N 2 is observed, and the rotation of the ligand-bound complex makes it look like a spinning umbrella. Hydrogen storage capabilities of clathrate hydrates, Li-doped clusters, and super alkali are investigated and it is found that the former can accommodate 2-6 hydrogen molecules, whereas the Li systems show a gravimetric wt% range of 7.1-28.0% for the star-like clusters and 13.2-40.9% for the super-alkali systems. The (HF) 2 encapsulation by the fullerene cages describes the confinement effect on the H-bond therein. Apart from C 60 , all the cages form the complexes in a thermodynamically favorable process. Also, a partial covalent character is observed in the H-bonds upon confinement.

AUTHOR CONTRIBUTIONS
PKC came up with the concept and design of the review, wrote the abstract, reviewed the final manuscript. RP and AP contributed towards the literature survey, writing the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.