On the Nature of the Partial Covalent Bond between Noble Gas Elements and Noble Metal Atoms

This article provides a discussion on the nature of bonding between noble gases (Ng) and noble metals (M) from a quantum chemical perspective by investigating compounds such as NgMY (Y=CN, O, NO3, SO4, CO3), [NgM−(bipy)]+, NgMCCH, and MCCNgH complexes, where M=Cu, Ag, Au and Ng=Kr−Rn, with some complexes containing the lighter noble gas atoms as well. Despite having very low chemical reactivity, noble gases have been observed to form weak bonds with noble metals such as copper, gold, and silver. In this study, we explore the factors that contribute to this unusual bonding behavior, including the electronic structure of the atoms involved and the geometric configuration of the concerned fragments. We also investigate the metastable nature of the resulting complexes by studying the energetics of their possible dissociation and internal isomerization channels. The noble gas-binding ability of the bare metal cyanides are higher than most of their bromide counterparts, with CuCN and AgCN showing higher affinity than their chloride analogues as well. In contrast, the oxides seem to have lower binding power than their corresponding halides. In the oxide and the bipyridyl complexes, the Ng-binding ability follows the order Au > Cu > Ag. The dissociation energies calculated, considering the zero-point energy correction for possible dissociation channels, increase as we move down the noble gas group. The bond between the noble gases and the noble metals in the complexes are found to have comparable weightage of orbital and electrostatic interactions, suggestive of a partial covalent nature. The same is validated from the topological analysis of electron density.


Introduction
Chemical bonding is a fundamental concept in chemistry that encompasses the interactions between atoms and the forces that hold them together. This phenomenon plays a crucial role in the formation of chemical compounds and the behavior of their constituent elements. Understanding the nature of chemical bonds is essential for predicting the stability, reactivity, and properties of chemical compounds, as well as for designing new materials with specific characteristics. We have studied 'covalent' and 'ionic' types of bonding in elementary chemistry, each of which are influenced by the electronic structure of the concerned fragments. But how true are the definitions of covalent and ionic bonds? In theory, there can exist a 100% ionic bond between two atoms, where one atom completely donates its valence electrons to the other atom, resulting in the formation of a completely ionic compound. However, in practice, no such bond can exist because all chemical bonds involve some degree of electron sharing as the orbitals of the cations and the anions may not be just spectators. In most cases, even compounds that are considered to have predominantly ionic bonding, e.g., ionic salts such as NaCl, have some degree of covalent character. Again, most 'covalent' bonds have some degree of ionic character, even if it is very small. This is because the electronegativities of the atoms in a bond are rarely exactly the same, leading to an unequal sharing of electrons. Additionally, other factors Molecules 2023, 28,3253 3 of 17 between two low reactive elements has long captured our attention. Here we provide a brief discussion on the stability and bonding of the type of complexes they can form.

Computational Tools to Analyze the Bonding
The bonding nature of these types of compounds are typically analyzed using tools based on computational techniques, viz., Natural Bond Orbital (NBO), Quantum Theory of Atoms-in-Molecules (QTAIM), Electron Localization Function (ELF), and Energy Decomposition Analysis (EDA), among others.
The Natural Bond Orbital (NBO) [46] is based on the concept of localized chemical bonding, where orbitals are constructed for individual chemical bonds in a molecule. It uses the density matrix of a molecule to obtain information about the probability distribution of electrons in the molecule. The NBO analysis, along with the Wiberg bond index (WBI) [47], provides a detailed picture of the bonding pattern in a molecule, including the degree of covalency, polarization, and charge transfer. WBI is a quantitative measure of the degree of electron sharing between two atoms in a molecule and was developed by Professor Steven R. Wiberg. It is based on the molecular orbital theory, specifically, on the concept of natural atomic orbitals, and it is evaluated from the diagonal elements of the density matrix of the molecular orbitals of a molecule. In other words, it is based on the density of electrons shared between two atoms. It ranges from zero to the sum of the formal bond orders between two atoms. A WBI value of zero indicates that there is no electron sharing between the two atoms, while a value equal to the sum of the formal bond orders indicates complete electron sharing. A higher value of WBI indicates a higher degree of covalent character in the concerned bond.
In QTAIM, a zero-flux surface defines an atom to obtain its atomic energy, even while it is a part of a molecule, by assuming the atoms to be the fundamental units rather than the bonds [48][49][50]. Here, the bonding is described in terms of the topology of the electron density in the system by partitioning it through the steepest ascent trajectory or the gradient path (GP). ELF is also based on the same principle [51,52]. Despite being widely used by scientists globally, the QTAIM approach has been subject to reports questioning its accuracy and foundation, specifically in relation to QTAIM atomic charges [53][54][55][56]. There also exist some debates regarding the accuracy of the atomic charges used in the QTAIM approach and their ability to reproduce the dipole/quadrupole moment or the electrostatic potential of a molecule. Some reports indicate a "size dependency" in the QTAIM charges [57] and a wrong indication of the ionic bond in the CO molecule due to the charge on the C center being unexpectedly high [53]. According to Bader and Matta [58], the consideration of atomic polarizations is essential in charge transfer calculations to avoid inaccurate QTAIM charges, as a point-charge model fails to account for them. Additionally, they note that quantum mechanics cannot provide a unique definition of atomic charge [55] as it is not a physical observable. This observation contradicts the principles of an open system in physics [59]. However, Bader and Matta [58] argue that atomic charges can be described as a quantum expectation value for an open system and can therefore be used to calculate external field-induced polarization and permanent moments in molecules. However, they also caution that atomic charges defined outside the context of an open system have limited usefulness, as they cannot be used in quantum mechanical expressions for expectation values and are not related to measurable properties. Electron density (ED, ρ(r)), Laplacian of ED ( 2 ρ(r)), local energy densities (kinetic G(r), and potential V(r)), and their sum (H(r)) calculated at the bond critical point (BCP) are some of the parameters used to describe the bonding. The conventional criterion for identifying the nature of a bond, which involves checking for the sign of 2 ρ(r) at the BCP is often not enough for Ng compounds. This could be due to the weaker orbital involvement of Ng in bonding compared to other elements in the periodic table, or that the criterion tends to fail for heavier elements beyond the first row of the periodic table. The reason for this is that 2 ρ(r) is calculated from the curvatures λ 1 , λ 2 , and λ 3 , of which the first two are negative while the final term is positive. For heavier elements, the positive contribution from the last term can often dominate, resulting in a positive 2 ρ(r) value. In such cases, the H(r) can be a more useful descriptor whose negative value indicates the presence of covalent (or partial covalent) bonding depending on the magnitude, even when 2 ρ(r) is positive. The full topological analysis of H(r) is reported as an effective mode to investigate the nature of the bonds occurring in Ng compounds in cases where they have examined its plotted shape, topology, and value along the AIM bond path, particularly in the vicinity of the BCP. This allows for a systematic classification of the concerned bonds using a range of numerical indices that are specific to each type of bond [60]. A few other interesting publications have been reported on the full bonding analysis in noble gas compounds [61][62][63].
Energy Decomposition Analysis (EDA) is another powerful computational tool used in the field of quantum chemistry to analyze the interaction between molecules [64,65]. It is used to decompose the total interaction energy between fragments of a molecular system into various energy components, which correspond to different types of interactions between them. These energy components can be used to analyze the nature and strength of the intermolecular interactions, thereby explaining the possible bonding therein. The components of the said energy are of attractive and repulsive types. The former constitutes electrostatic (Coulombic interaction between the EDs of the two fragments; E elstat ), orbital (E orb ), and dispersion (attractive van der Waals forces between the fragments which arise due to the fluctuations in the corresponding EDs; E disp ) interaction energies, while the latter is essentially the Pauli repulsion interaction (E Pauli ). To examine the bonding of the Ng compounds, we took into account both neutral and ionic fragmentation methods during the EDA computations. The neutral and ionic fragmentation schemes represent electron-shared covalent interaction and dominant electrostatic interaction, respectively. To identify the better partitioning scheme, we compared the ∆E orb values of both. The partitioning scheme with a smaller attractive ∆E orb value represents the bonding structure more accurately, as it causes less electronic charge reorganization to yield the electronic structure of the compound [66,67].
In our discussions, we have described the bonding in terms of %covalent and %ionic characters. These percentages are calculated as: It may be noted that these definitions are mostly qualitative in nature. In the complexes of Ng with metal cations, the term ∆E orb obtained by the EDA scheme employed is dominated by the inductive interaction, which is "physical" in nature and does not necessarily demand the sharing of electrons that is typical of covalent bonds. Similarly, in the majority of the Ng complexes, the ∆E elstat arises from the interaction between the frozen densities of the interacting fragments, which is different from what is commonly perceived as an ionic bond or an ionic contribution to the bond. For simple molecules in general, those definitions are expected to suffice. It is to be noted, however, that this simplified analysis is not perfectly applicable for bonds that are charge-shifting or eventually spin-polarized. For example, the F 2 molecule is a well-known example of a charge-shifting bond. There are, of course, other ways to classify a chemical bond with detailed bonding analysis, particularly for important cases such as a charge-shifting bond. The chemical bond is not a physical observable, and there exists a high degree of arbitrariness during its analysis, especially when dealing with bonds with a high dispersion component, as is the case with noble gas bonds. Equations (1) and (2) are suggested in a qualitative sense, as they are capable of providing a rough idea about the extent of covalency in a typical bond.

The General Bonding Picture
The notion of a conventional bond is not always well defined when referring to standard compounds with non-covalent interactions, such as traditional ionic compounds. In cases where two ions are attracted to each other by electrostatic forces, there is usually some amount of orbital overlap. Similarly, some amount of ionic character exists in traditional covalent compounds. We have provided a description of the bonding picture in some commonly known diatomic 'covalent' and 'ionic' compounds and compared them with some Ng compounds containing partially covalent bonds (listed in Table 1) [71] with the help of the computational tools discussed in the previous section. While other techniques exist for investigating the nature of bonding, these commonly used and cited tools help eliminate misconceptions among young researchers and advise against the hasty use of terms such as 'covalent' and 'ionic' to describe a bond. It is important to note that most molecules cannot be exclusively classified as either 'covalent' or 'ionic'. Table 1 displays the numerical values of different topological descriptors at the concerned BCPs. For the halides of alkali and alkaline earth metals, ρ(r c ) and H(r c ) values are low and positive, while the Laplacian is positive, and all these values decrease from F to Br. The ratio of G(r c ) and ρ(r c ) is mostly greater than 1, except for KCl, KBr, CaF 2 , and CaBr 2 , where it is slightly less than 1. These parameters suggest that the interactions are predominantly closed-shell. In contrast, hydrogen halides have higher ρ(r c ), negative 2 ρ(r c ) and H(r c ) values and G(r c )/ρ(r c ) ratios that are less than 1 for HCl and greater than 1 for HF and HBr. The initial two parameters suggest a covalent bonding nature, whereas the latter two parameters imply the presence of partially covalent interactions. AIM results for H 2 , N 2 , O 2 , F 2 , Cl 2 , CO, NH 3 , CH 4 , C 2 H 4 , and C 2 H 2 also imply a partially covalent bond character. For the Ng compounds, the BCPs at the M-N and M-S bonds correspond to slightly larger positive ρ(r c ) and 2 ρ(r c ) but negative H(r c ) values, and a greater G(r c )/ρ(r c ) ratio to those of Ng-M bonds. These observations suggest that both bonds possess a partial covalent character, with the Ng-M bond displaying a greater extent of covalency. The ELF values for the halides of alkali and alkaline earth metal are notably low, with Li and Na compounds having the smallest values, followed by the K, Mg, and Ca compounds. Similar results are seen for the Ng compounds as well. These low values indicate minimal covalent character at the BCPs, suggesting predominantly ionic interactions. Conversely, the hydrogen halides, diatomic molecules, hydrocarbons, H 2 O, and NH 3 have very high ELF values, with most values above 0.8 and some even closer to 1. However, CO is an exception with an ELF value of 0.4. These high ELF values suggest a substantial covalent character at the BCPs.
The values of %ionic and %covalent characters obtained from using Equations (1) and (2) are provided in Table 2. The percentages of covalent and ionic interactions in the alkali and alkaline earth metal halides range from 78.9% to 93.3% and 6.7% to 21.2%, respectively. They seem to follow Fajans' rules, where the covalency increases from F to Br for a given metal. In comparison to the metal halides, there is a significant increase in the orbital contribution for the HX molecules, accompanied by a decrease in the electrostatic contribution (∆E elstat and ∆E orb ranging from 43.4% to 61.5% and 38.5% to 56.5%, respectively). These changes suggest a decrease in ionic character and an increase in covalency, where the covalency even exceeds the ionic character in HCl and HBr. In the common diatomic species such as H 2 , N 2 , O 2 , F 2 , Cl 2 , and CO, the orbital contribution seems to provide the highest stabilization, indicating high covalency (73.5% to 85.1%), with an almost pure covalent bond in the H 2 molecule. Hence, compounds that are commonly referred to as 'pure covalent' still possess a significant electrostatic component. H 2 O, NH 3 , CH 4 , C 2 H 4 , and C 2 H 2 show high covalency, with an order of N-H < O-H < C-H, as is expected based on the electronegativities of the atoms involved. The π-bonded systems with C-C bonds have a significant ionic component, with acetylene having a greater degree of covalency than ethylene. The analysis of Ng compounds containing Ng-M bonds (M=Cu, Ag, Au) using EDA involves evaluating two fragmentation schemes, which examine bonds on each side of the metal atoms. In the case of NgMNO 3 and Ng 2 M 2 SO 4 species, the Ng-M bonds display comparable contributions from orbital and electrostatic terms, with electrostatic interactions being slightly greater. As a result, these bonds exhibit both a covalent and ionic nature, with % covalency ranging between 38.5% and 46.1%, which is slightly lower. Conversely, the M-N/S bonds are dominated by ∆E elstat , which accounts for 77.6-86.7% of the interaction and, hence, display a more pronounced % ionic contribution. Even though the dispersion term is the least significant in the total attractive interaction, it is considerably greater in the Ng-M bonds than in the other bonds discussed earlier. As a result, the metal atom is connected on both sides by partially covalent bonds.

Types of Noble Gas Compounds
Noble gas compounds can be broadly categorized into three groups based on the position of the Ng atom in the compound/complex, viz., non-insertion, insertion, and confined. While the non-insertion compounds contain Ng at one open end of the XY moiety (i.e., NgXY type), the insertion type has the Ng atom inserted within a bond between two elements such that it is flanked by two atoms or groups on either side (i.e., XNgY type of compounds, where the Ng is inserted within the X-Y bond). The stability in the former comes from the polarizing power of the atom adjacent to the Ng (i.e., X) which creates an attractive interaction (donor-acceptor type) between X and Ng. The smaller the size of the X atom, the higher is its polarizing power. There exists an electronegativity difference between the X and Y atoms which produces a dipole and can hence cause the said polarization, the extent of which is higher for ionic systems This XNg moiety is then stabilized by the Y atom, which acts as the counterion. In the case of the insertion complex, however, the stability is kinetics-driven rather than thermodynamics-driven. The reason behind this is the low reactivity of Ng with X and Y compared with the more stable X-Y bond, which has to be destroyed in order to form the Ng complex. Hence, the formation of such XNgY complexes is not thermochemically feasible. However, once it is formed, the high free-energy activation barrier along the reaction path makes the system kinetically stable along that dissociation path.
In the cases of such complexes, it is necessary to consider a significant number of dissociation pathways. Thermochemical analyses reveal that, aside from the two dissociation pathways discussed below, the other dissociation channels are highly or moderately endergonic, indicating that XNgY is stable. The stability of XNgY is determined by two competing dissociation paths: 2B and 3B. While the former involves XNgY splitting into Ng and XY and is typically highly exergonic, the latter results in X, Ng, and Y, and is occasionally exhibits low exergonicity at room temperature. Thus, evaluating the activation energy barrier for these paths is crucial to assessing the stability of XNgY. While single reference-based methods are often reliable for studying the TS of 2B dissociation, computing the barrier of 3B dissociation is more complicated and typically requires multireference treatment. At low temperatures, the most favorable systems are those where the 3B dissociation pathway becomes energetically unfavorable. According to Hu et al. [72], the half-life of Ng compounds with an XNgY structure depends on the energy barrier height. Systems with a minimum energy barrier of 6, 13, and 21 kcal mol −1 would correspond to a half-life of approximately 102 s at 100, 200, and 300 K, respectively.
The third category mainly depends on the effects brought about by the confined environment provided by certain cages, e.g., Ng 2 dimers encapsulated within cages, such as fullerene, B 12 N 12 , B 16 N 16 , C 20 H 20 , B 40 , clathrate hydrates, cucurbit[n]uril, octa acids, BN-doped carbon nanotubes, etc., are some of the systems that our group has worked on over the years. Here, the structure, bonding, and reactivity of the encapsulated systems are compared with those before confinement (i.e., in their free states), e.g., as performed for the bonding and movement of Ng 2 units within fullerene cages [73,74].

Case Studies of Noble Gas−Noble Metal Binding
Among the bonding of noble gas atoms with other elements of the periodic table, the ones with the noble (or precious) metals, i.e., Cu, Ag, and Au, are of particular interest to us due to the challenges that come with attempting bond formation between two less reactive elements. Research in this field essentially came to the forefront with Pyykkö's prediction of complexes containing Ng and Au [25,26], followed by pioneering works by several groups in detecting complexes with noble gases and noble metals interacting with each other [29,[41][42][43][44][45]75].
Over the years, we have performed several theoretical investigations on systems containing possible interactions between noble gases and noble metals. The Ng-binding ability of cyanides, oxides, nitrates, sulphates, and carbonates of the noble metals was investigated by studying the stability and bonding of the non-insertion type complexes NgMY, where M=Cu, Ag, Au, and Y=CN, O, NO 3 , SO 4 , CO 3 [75][76][77]. The complexes are optimized at different levels, viz., at the CCSD(T)/def2-TZVPPD [78][79][80] and CCSD(T)/ccpVTZ/cc-pwCVTZ-PP [78,[81][82][83][84] levels for Y=CN and O, respectively, and at the MPW1B95/def2-TZVP level for Y=NO 3 , SO 4 , CO 3 . The optimized geometries of the cyanides and the oxides are depicted in Figure 1. The dissociation energies (D 0 BSSE ) of the Ng-M bonds were calculated for NgMCN, NgMO, and Ng-bound metal nitrates, sulphates, and carbonates, considering the zero-point energy (ZPE) and basis set superposition error (BSSE) corrections. The D 0 BSSE values range from 2.7 to 14.6 kcal/mol for NgMCN, with an increase in value from Ar to Rn. Those for the oxides vary from 2.0 to 13.2 kcal/mol, with the same order down the Ng group in the periodic table. For the Ng-bound metal nitrates, sulphates, and carbonates, the values range from 5.1 to 13.2 kcal/mol, 2.3 to 10.1 kcal/mol, and 7.3 to 19.9 kcal/mol, respectively, for the metals moving from Cu to Au. The structural characteristics of both NgMCN and NgMO complexes were compared with those of the experimentally reported NgMX (X=F, Cl, Br) compounds. The results show that AuCN has the highest Ng-binding ability among the noble metal halides, followed by CuCN and AgCN. The ability of CuCN to bind with Ng is greater than that of CuCl and CuBr, but not as high as that of CuF. The orders for the same in the Ag and Au compounds are AgCN ≈ AgF > AgCl > AgBr and AuF > AuCl > AuCN > AuBr, respectively. For the NgMO complexes, the binding power of MO is slightly lower than that of MBr, leading to an order of MF > MCl > MBr > MO. The D 0 values also increase gradually from Ar to Rn. With respect to thermochemical stability, we see that at room temperature, dissociation of MCN compounds bound to Kr-Rn (except for KrAgCN) and ArAuCN is found to be endergonic (others may need lower temperatures to remain thermochemically stable). That of NgMO complexes into Ng and MO is also endothermic and becomes increasingly so moving from Ar to Rn. The same is calculated for XeAgO, RnAgO, and Kr-Rn-bound CuO and AuO complexes along with most of the sulphates, nitrates, and carbonates, with some exceptions for the Ar and Kr analogues. Lower temperatures may be needed to keep some systems in the bound form.
The bonding picture Is demonstrated in terms of NBO, WBI, EDA, and AIM analyses. The WBI values vary within 0.2-0.3 in the Xe-and Rn-bound MCN complexes. Most of the Ng-M distances, except for all the Ng-Cu bond distances along with that of the Xe and Rn analogues of Ag and Au complexes in NgMCN, are smaller than their respective covalent bond distances. The ∆E HOMO-LUMO is higher in these Ng-bound complexes than in their bare counterparts, suggesting higher electronic stability in the former. For the Ng-M bonds in Ar-Rn-bound MO, the WBI ranges are 0.16-0.28, 0.08-0.22, 0.15-0.34 for Cu, Ag, and Au analogues, respectively, with a gradual increase along the group. The same can be said for the Ng-bound complexes of the noble metal sulphates, nitrates, and carbonates. The concerned bond between Ng and M in all the studied complexes can be considered to contain partially covalent and electrostatic types of interactions. The movement of the electron cloud from the noble gas to the middle of the noble gas and noble metal centers suggests an increased likelihood of bond formation. The degree of covalency is higher from Ar to Rn, as evidenced by the increase in the orbital energy term. Both electrostatic and orbital energy values are almost equal in the Ng−Cu bonds, while the orbital contribution is slightly lower than the electrostatic one in NgAgO and NgAuO (excluding ArAuO) (Table S1). In the case of the cyanides, the orbital term is greater in magnitude than the electrostatic term in the NgCuCN complexes, while the reverse is true for the Ag and Au analogues (except for ArAgCN and ArAuCN) ( Table S2). All Ng-M bonds in the NgMCN complexes have negative H(r c ) values at their bond critical points (except for Ar-Ag), but G(r c )/ρ(r c ) < 1 in the Xe/Rn-Cu and Rn-Ag bonds (Table S3). A certain degree of covalency in the concerned bonds of oxides, nitrates, sulphates, and carbonates is further indicated by the electron density analysis (Tables S4 and S5). The optimized geometries of Ng 2 Cu 2 SO 4 and Ng 2 Ag 2 SO 4 possess a D 2d symmetry, whereas those of Ng 2 Au 2 SO 4 and Ng 2 M 2 CO 3 have C 2 and C 2v symmetries, respectively. The D 0 values for the Ng-M bond vary in the range 2.5−19.9 kcal/mol in Ng 2 M 2 SO 4 , and 2.3-18.2 kcal/mol in Ng 2 M 2 CO 3 , increasing down the Ng group. The corresponding process for the dissociation of Ng atoms is spontaneous for Ar 2 Cu 2 SO 4 , Ar 2 Ag 2 Y, and Kr 2 Ag 2 Y. Similar to the NgMY systems, the HOMO-LUMO gap also increases in the Ng 2 M 2 Y systems compared with their corresponding bare M 2 Y systems. Except for Ar 2 Ag 2 SO 4 and Ar 2 Ag 2 SO 4 all other Ng 2 M 2 Y systems have negative H(r c ) values as calculated from the AIM analysis, suggesting their partial covalent nature (Table S5).
The concerned bond between Ng and M in all the studied complexes can be considered to contain partially covalent and electrostatic types of interactions. The movement of the electron cloud from the noble gas to the middle of the noble gas and noble metal centers suggests an increased likelihood of bond formation. The degree of covalency is higher from Ar to Rn, as evidenced by the increase in the orbital energy term. Both electrostatic and orbital energy values are almost equal in the Ng−Cu bonds, while the orbital contribution is slightly lower than the electrostatic one in NgAgO and NgAuO (excluding ArAuO) (Table S1). In the case of the cyanides, the orbital term is greater in magnitude than the electrostatic term in the NgCuCN complexes, while the reverse is true for the Ag and Au analogues (except for ArAgCN and ArAuCN) (Table S2). All Ng-M bonds in the NgMCN complexes have negative H(rc) values at their bond critical points (except for Ar-Ag), but G(rc)/ρ(rc) < 1 in the Xe/Rn-Cu and Rn-Ag bonds (Table S3). A certain degree of covalency in the concerned bonds of oxides, nitrates, sulphates, and carbonates is further indicated by the electron density analysis (Tables S4 and S5). The optimized geometries of Ng2Cu2SO4 and Ng2Ag2SO4 possess a D2d symmetry, whereas those of Ng2Au2SO4 and Ng2M2CO3 have C2 and C2v symmetries, respectively. The D0 values for the Ng-M bond vary in the range 2.5−19.9 kcal/mol in Ng2M2SO4, and 2.3-18.2 kcal/mol in Ng2M2CO3, increasing down the Ng group. The corresponding process for the dissociation of Ng atoms is spontaneous for Ar2Cu2SO4, Ar2Ag2Y, and Kr2Ag2Y. Similar to the NgMY systems, the HOMO-LUMO gap also increases in the Ng2M2Y systems compared with their corresponding bare M2Y systems. Except for Ar2Ag2SO4 and Ar2Ag2SO4 all other Ng2M2Y systems have negative H(rc) values as calculated from the AIM analysis, suggesting their partial covalent nature (Table S5). Apart from the aforementioned complexes, we have also examined the structure, stability, and bonding nature of monocationic M-bipyridine complexes bound to Ng, Apart from the aforementioned complexes, we have also examined the structure, stability, and bonding nature of monocationic M-bipyridine complexes bound to Ng, optimized at the MPW1B95/cc-pVTZ level of theory ( Figure 2) [85]. The D 0 BSSE values range from 4.0 to 19.7 kcal/mol, which increase as we move towards heavier Ng atoms. Among the three analogues, the order of Ng binding is Au > Cu > Ag, with the exception of the Ar analogue, where the order is Cu > Au> Ag. We have found that all dissociation processes are increasingly endothermic down the Ng group, and the associated ∆G values suggest the thermochemical stability of the Kr-Rn bound complexes. Moving from Ar to Rn, the extent of electron transfer from Ng to the M center increases along with the corresponding WBI values in the Ng-M bonds. Except in the Ar-Ag bond, the H(r c ) value is negative (Table S6). Additionally, EDA results show that the attraction between Ng and M centers consists of approximately 45-55% electrostatic and 41-45% orbital interactions, indicating the presence of both covalent and ionic characters (Table S7). An in-depth examination of the orbital term reveals that Ng→Au σ-donation and Ng←Au σ-and π-back donations contribute 61-69% and 27-35% to the total orbital term, respectively.
Rn, the extent of electron transfer from Ng to the M center increases along with the corresponding WBI values in the Ng-M bonds. Except in the Ar-Ag bond, the H(rc) value is negative (Table S6). Additionally, EDA results show that the attraction between Ng and M centers consists of approximately 45-55% electrostatic and 41-45% orbital interactions, indicating the presence of both covalent and ionic characters (Table S7). An in-depth examination of the orbital term reveals that NgAu σ-donation and NgAu σ-and π-back donations contribute 61-69% and 27-35% to the total orbital term, respectively. Further, the Ng-inserted complexes of MCCH are investigated by considering two possible structures, one where the Ng is inserted within the M-C bond [86] and the other where it is inserted within the C-H bond [87]. The structures are optimized using MPW1B95, [88] MP2 [89], and CCSD(T) [78] methods with basis sets; cc-pVTZ [81,83,90] for H, C, and Ar atoms; and cc-pVTZ-PP with relativistic effective core potential for Kr, Xe, Rn, Cu, Ag, and Au atoms (Figure 3). While for the former, only Xe and Rn were Further, the Ng-inserted complexes of MCCH are investigated by considering two possible structures, one where the Ng is inserted within the M-C bond [86] and the other where it is inserted within the C-H bond [87]. The structures are optimized using MPW1B95, [88] MP2 [89], and CCSD(T) [78] methods with basis sets; cc-pVTZ [81,83,90] for H, C, and Ar atoms; and cc-pVTZ-PP with relativistic effective core potential for Kr, Xe, Rn, Cu, Ag, and Au atoms (Figure 3). While for the former, only Xe and Rn were inserted, Kr is also considered for the latter case. The non-insertion isomer, NgMCCH (Ng=Ar-Rn), is found to be viable, especially under low-temperature conditions. The isomerization between the insertion and the non-insertion complex, MNgCCH → NgMCCH, is spontaneous but is kinetically protected by a high free energy activation barrier ranging within 14.0-34.8 kcal mol −1 , indicating the viability of the insertion complex at room temperature. It is thermochemically stable with respect to other possible dissociation processes, except for a three-body dissociation of the AgXeCCH complex, where a slight lowering of temperature is required to halt the dissociation. Although the Rn analogue is kinetically more stable than the Xe analogue for a given M, the ∆G } value increases in the order of Ag < Cu < Au for a particular Ng. Additionally, non-inserted NgMCCH compounds that are bound to Cu and Au by Kr-Rn, and to Ag by Xe and Rn, are thermally stable and do not dissociate into Ng and MCCH at 298 K. The M-Ng bonds exhibit high WBI values, i.e., high covalency, whereas the Ng-C bonds are primarily ionic, i.e., best represented as (MNg) + (CCH) − . EDA results suggest that the interaction between M and Ng is not purely covalent in nature; significant electrostatic interactions contribute to their formation (Table S8). The same is also evident from the electron density descriptors presented in Table S9. Electronic structure principles associated with conceptual DFT, such as maximum hardness [91][92][93] and minimum electrophilicity [93][94][95][96][97] principles, are observed to hold true along the isomerization process MNgCCH → NgMCCH. These compounds represent the first examples of viable species involving an M-Ng-C bonding motif, which expands the current understanding of Ng insertion compounds with unique bonding units. The large ∆G } values observed in the exergonic dissociation processes of the MNgCCH complex emphasizes the necessity for thorough scrutiny when dealing with species containing the M-Ng-C bonding pattern.
cesses, except for a three-body dissociation of the AgXeCCH complex, where a slight lowering of temperature is required to halt the dissociation. Although the Rn analogue is kinetically more stable than the Xe analogue for a given M, the ∆G ╪ value increases in the order of Ag<Cu<Au for a particular Ng. Additionally, non-inserted NgMCCH compounds that are bound to Cu and Au by Kr-Rn, and to Ag by Xe and Rn, are thermally stable and do not dissociate into Ng and MCCH at 298 K. The M-Ng bonds exhibit high WBI values, i.e., high covalency, whereas the Ng-C bonds are primarily ionic, i.e., best represented as (MNg) + (CCH) − . EDA results suggest that the interaction between M and Ng is not purely covalent in nature; significant electrostatic interactions contribute to their formation (Table S8). The same is also evident from the electron density descriptors presented in Table S9. Electronic structure principles associated with conceptual DFT, such as maximum hardness [91][92][93] and minimum electrophilicity [93][94][95][96][97] principles, are observed to hold true along the isomerization process MNgCCH  NgMCCH. These compounds represent the first examples of viable species involving an M-Ng-C bonding motif, which expands the current understanding of Ng insertion compounds with unique bonding units. The large ∆G ╪ values observed in the exergonic dissociation processes of the MNgCCH complex emphasizes the necessity for thorough scrutiny when dealing with species containing the M-Ng-C bonding pattern. In the cases of the MCCNgH complexes (Figure 3), while the Cu and Ag analogues only prefer to bind with Xe and Rn noble gases, the Au analogue can easily bind with Kr as well. The AuCCH molecule is found to be the most effective in facilitating Ng insertion. In the cases of the MCCNgH complexes (Figure 3), while the Cu and Ag analogues only prefer to bind with Xe and Rn noble gases, the Au analogue can easily bind with Kr as well. The AuCCH molecule is found to be the most effective in facilitating Ng insertion. While the thermochemically most feasible process (MCCNgH breaking off to produce MCCH and Ng) is significantly hindered by large ∆G } values at room temperature, the stability of the complex relies upon two other radical dissociation channels (the 2B splitting of MCCNgH to MCC and NgH, and the 3B splitting of MCCNgH to MCC, Ng, and H), which are exergonic at 298 K for Ng=Kr, Xe and M=Cu, Ag, AuCCKrH. The bonding pattern demonstrates that Ng participates in significant covalent interactions with the neighboring atoms (Tables S10 and S11). Additionally, it is feasible to insert two Xe atoms into the Au−C and C−H bonds of AuCCH, producing a kinetically stable AuXeCCXeH molecule. Another isomer of this compound, the XeAuCCXeH, is also stable against splitting into Xe and AuCCXeH. The insertion of three Xe's, however, is unable to produce a stable complex. These findings thus expand our understanding of the unique bonding present in the Ng insertion compounds.

Conclusions
Contrary to the popularly used terms, viz., 'covalent' and 'ionic' bonds, there truly only exist partially covalent bonds comprising both covalent and ionic characters in various compositions, depending upon the electronic structure of the atoms/groups involved. This review article is aimed at studying the nature of the bonds in different compounds, especially those between two less reactive element groups, i.e., noble gases (Ng) and noble metals (M). We have made use of common computational tools such as Natural Bond Orbital (NBO), Quantum Theory of Atoms-in-Molecules (QTAIM), Electron Localization Function (ELF), and Energy Decomposition Analysis (EDA). Investigating possible dissociation channels of the Ng compounds reveals that most of them are endergonic, whereas the exergonic ones have a very high free energy barrier, indicating the stability of the compounds at room temperature. The Ng-binding ability of the compounds containing the noble metals generally follow the order of Au > Cu > Ag. The dissociation channels are increasingly endergonic as the Ng becomes heavier. The WBI values of the bond between Ng and M are mostly lower than 0.5, with a gradual increase along the group. Those in the C-Ng and Ng-H bonds had a range of 0.33-0.44 and 0.60-0.70, respectively. Detailed bonding analysis shows that the contributions from the orbital and electrostatic energy terms are similar in most of the cases (with a range of 45-55%) in the Ng-M bonds of the complexes. This partial covalent nature in these bonds is further indicated by the electron-density analysis. This unusual bonding between noble gases and noble metals is not a simple van der Waals interaction but is instead influenced by a range of factors, including electrostatic, orbital, charge transfer, and polarization effects, which contribute to stabilizing the insertion and non-insertion type of noble gas compounds.