On the Proton-Bound Noble Gas Dimers (Ng-H-Ng)+ and (Ng-H-Ng’)+ (Ng, Ng’ = He-Xe): Relationships between Structure, Stability, and Bonding Character

The structure, stability, and bonding character of fifteen (Ng-H-Ng)+ and (Ng-H-Ng’)+ (Ng, Ng’ = He-Xe) compounds were explored by theoretical calculations performed at the coupled cluster level of theory. The nature of the stabilizing interactions was, in particular, assayed using a method recently proposed by the authors to classify the chemical bonds involving the noble-gas atoms. The bond distances and dissociation energies of the investigated ions fall in rather large intervals, and follow regular periodic trends, clearly referable to the difference between the proton affinity (PA) of the various Ng and Ng’. These variations are nicely correlated with the bonding situation of the (Ng-H-Ng)+ and (Ng-H-Ng’)+. The Ng-H and Ng’-H contacts range, in fact, between strong covalent bonds to weak, non-covalent interactions, and their regular variability clearly illustrates the peculiar capability of the noble gases to undergo interactions covering the entire spectrum of the chemical bond.


Method of Bonding Analysis
The method that we recently proposed to analyze the bonding character of noble-gas compounds is extensively discussed in key references [75][76][77] so we briefly recall here only the most salient features.
Our taken approach relies on the examination of the plotted shape of the local electron energy density H(r) [75,78,79], and on the values that this function takes over the volume Ω s enclosed by the s(r) = 0.4 reduced density gradient (RDG) isosurface [80,81] associated with the bond critical point (BCP) that is located for a given Ng-X bond (X = binding partner) from the topological analysis of the electron density ρ(r) [82]. Ancillary indices include the size of Ω s , the total electronic charge enclosed by Ω s , N(Ω s ), the average electron density over Ω s , ρ s (ave) = N(Ω s )/Ω s , and the average, maximum, and minimum value of H(r) over Ω s , H s (ave, max, min). As discussed previously [75], the H(r) partitions the atomic space in two well-recognizable regions, namely an inner one of negative values, indicated as H − (r), and an outer one of positive values, indicated as H + (r). The boundary of these regions falls at a distance R − , that is typical of each atom; at this distance, H(r = R − ) = 0. Interestingly, when two atoms form a chemical bond, their H − (r) and H + (r) regions combine in modes that signal the nature of the interaction. Particularly for the Ng-X bonds, it is possible to distinguish three types of interactions, indicated as A, B, or C. In interactions of type A, the atoms overlap all the contour lines of their H + (r) regions, and part of the contour lines of their inner H − (r) regions, the bond appearing as a continuous region of negative values of H(r), plunged in a zone of positive values. The interaction is topologically-signed by a (3 + 1) critical point of the H(r) (denoted as the HCP) falling on the bond axis. In interactions of type B, the H − (r) region of Ng is, again, overlapped with the H − (r) region of the binding partner, but (i) no HCP does exist on the bond axis, and (ii) the Ng-X internuclear region does include a (more or less wide) region of positive H(r). In interactions of type C, the Ng and the binding partner overlap only part of their H + (r) regions, their H − (r) regions remaining, instead, perfectly closed, and separated by a (more or less wide) region of positive H(r). The bond thus appears as two clearly distinguishable H − (r) regions, separated by a region of positive values of H(r).
Any Ng-X is assigned as covalent (Cov) if (i) it is of type A, and (ii) the electron density at the BCP, ρ(BCP) is at least 0.08 ea 0 −3 . The strength of a Cov bond is also quantified in terms of the bond degree (BD). Borrowing a concept introduced so far by Espinosa et al. [83], this index is defined by the equation: where H(HCP) and ρ(HCP) are, respectively, the H(r) and the ρ(r) at the HCP of the Ng-X bond. Any Ng-X not fulfilling the criteria of covalency [i.e., it is of type B or C, or, if of type A, its ρ(BCP) is lower than 0.08 e a 0 −3 ] is further assayed by integrating the H(r) over Ω s . If the function is, invariably, positive over the entire volume, the interaction is assigned as non-covalent (nCov). If the function is partially or fully negative, the Ng-X is assigned as partially-covalent (pCov), and distinguished as H +/− , H −/+ , and H − , the superscript indicating that, over Ω s , the H(r) is ranging from negative to positive, but, on the average, it is positive (H +/− ) or negative (H −/+ ), or that it is, invariably, negative (H − ). Interactions of type C are also assayed in terms of the degree of polarization of Ng, DoP(Ng), an index that measures, in essence, the deformation of the H − (r) region of Ng arising from the interaction with X [76]. It is, in particular, defined by the equation: is the radius of the H − (r) region of Ng along the axis formed by Ng and the Ng-X BCP, and R − Ng is the radius of the H − (r) region of the free atom. Numerous illustrative examples of bonds classification are given in [77].

Computational Details
The calculations were performed at the coupled cluster level of theory, with inclusion of single and double substitutions, and an estimate of connected triples, CCSD(T) [84], using the Dunning's correlation consistent aug-cc-pVnZbasis sets (n = T, Q, 5) [85] (henceforth denoted as aVnZ). The Xe atom was treated by the Stuttgart/Cologne small-core (28 electrons), scalar-relativistic effective core potential (ECP-28) [86], and the jointly-designed aVnZ-PP basis sets. The calculations were performed by explicitly correlating the 2s 2 2p 6 (Ne), 3s 2 3p 6 (Ar), 4s 2 4p 6 (Kr), and 5s 2 5p 6 (Xe) outer electrons. The CCSD(T) correlation energies were extrapolated to the complete basis set (CBS) limit using the cubic extrapolation formula [87]: where E corr 4 and E corr 5 are, respectively, the CCSD(T)/aVQZ and CCSD(T)/aV5Z correlation energies. All the CCSD(T) calculations were performed with the CFOUR program (V2.1) [88]. The bonding analysis was accomplished at the CCSD(T)/aVTZ level of theory, the main analyzed functions being the ρ(r) [82], the H(r) [75,78,79], and the RDG, and itsrelated non-covalent interactions (NCI) indices [80,81]. The ρ(r) is defined by the equation [82]: where η i is the occupation number of the natural orbital ϕ i , in turn expanded as a linear combination of the basis functions.
The H(r) is the sum of the kinetic energy density G(r) and the potential energy density V(r): The presently-employed definition [82,89] of the G(r) is given by the equation: where the sum runs over all the occupied natural orbitals ϕ i of occupation numbers η i . The potential energy density V(r) is evaluated [82] from the local form of the virial theorem: The RDG is defined by the equation [80,81]: Low-value s(r) isosurfaces (typically 0.3-0.6) appear among atoms undergoing any type of interaction, the NCIs emerging, in particular, by considering the spatial regions of low ρ(r). The low-s(r)/low-ρ(r) isosurfaces are, in turn, mapped in terms of the sign (λ 2 ) × ρ(r), λ 2 being the second eigenvalue (λ 1 < λ 2 < λ 3 ) of the Hessian matrix of ρ(r). In essence, the sign of λ 2 is used to distinguish between attractive (λ 2 < 0) and repulsive (λ 2 > 0) interactions, and the value of ρ(r) is exploited to rank the corresponding strength. In the present study, we calculated also the integral of a given property P [particularly the ρ(r) and the H(r)] over the volume Ω s enclosed by the s(r) = 0.4 isosurface (s) at around the BCP located on any Ng-H bond path, P(Ω s ). This integration was accomplished by producing an orthogonal grid of points that encloses the isosurface and applying the formula: where P(r i ) is the value of P at the grid point r i , and d x , d y , and d z are the grid step sizes in the x, y, and z directions, respectively (d x = d y = d z = 0.025 a 0 ). The summation is carried out on all grid points r i having RDG <s. The ρ(r), the H(r), and the s(r) were analyzed with the Multiwfn program (version 3.7.dev) [90] using the wave functions stored in the molden files generated with CFOUR, and properly formatted with the Molden2AIM utility program [91]. The two-(2D) plots of the H(r) were as well produced with Multiwfn, and include the standard contour lines belonging to the patterns ±k × 10 n (k = 1, 2, 4, 8; n = −5 ÷ 6), together with the contour lines corresponding to the critical points specifically located from the topological analysis of the H(r).

Results and Discussion
The presently-investigated species include the five homonuclear (Ng-H-Ng) + and the ten heteronuclear (Ng-H-Ng') + . All these ions possess linear structures, and their CCSD(T)/aVTZ bond distances and harmonic vibrational frequencies are shown in Figure 1 and quoted in Table 1. Table 2 shows the values of ∆R(Ng-H), namely the elongation with respect to the diatomic NgH + of the Ng-H bond of the various (Ng-H-Ng') + .

Results and Discussion
The presently-investigated species include the five homonuc the ten heteronuclear (Ng-H-Ng') + . All these ions possess linear CCSD(T)/aVTZ bond distances and harmonic vibrational frequenc ure 1 and quoted in Table 1. Table 2 shows the values of ΔR(Ng-H tion with respect to the diatomic NgH + of the Ng-H bond of the vari    The thermochemical stabilities of the investigated ions were assessed by computing the CCSD(T)/CBS energy change of the two-body (2B) dissociations:

Ng') + R(Ng-H) R(Ng'-H)
(Ng-H-Ng') + → NgH + /Ng'H + + Ng'/Ng (10) distinguished here using the notations 2B(Ng) and 2B(Ng'), and the three-body (3B) dissociation: The obtained values are quoted in Table 1. The results of the CCSD(T)/aVTZ bonding analysis are given in Figure 2 and in Table 3. Table 1. CCSD(T)/aVTZ bond distances (Å), harmonic vibrational frequencies (cm −1 ), and CCSD(T)/CBS dissociation energies (kcal mol −1 ) of the (Ng-H-Ng') + and NgH + .  The diatomic NgH + are typically covalent species. As shown in Figure 2 and Table 3, their plotted H(r) is of type A, and their ρ(BCP) is definitely higher than 0.08 e a 0 −3 . The values of BD, ranging between 1.217 hartree e −1 (HeH + ) and 0.824 hartree e −1 (XeH + ), are also relatively high, and positively-correlated with rather high stretching frequencies ranging between 3206 cm −1 (HeH + ) and 2297 cm −1 (XeH + ). Compared with the NgH + , the homonuclear (Ng-H-Ng) + are structurally less compact and thermochemically less stable. Thus, the Ng-H distance is invariably longer, ∆R(Ng-H) being computed as ca. 0.15 Å for Ng = He and Ne, and ca. 0.22-0.26 Å for Ng = Ar, Kr, and Xe, and the ν 3 stretching frequency reduces in the range between 1616 cm −1 (Ne-H-Ne) + and 846 cm −1 (Xe-H-Xe) + . In addition, the half of the energy change of the 3B dissociation [Equation (11)], representing, in practice, the average proton affinity (aPA) of two Ng atoms sharing the (formal) H + , is, invariably, lower than the PA of Ng, namely the energy change of the reaction:   , and BD (hartreee −1 ) are, respectively, the electron density at the BCP, and the electron density, the energy density, and the bond degree at the HCP. For pCov or nCov bonds,the quotedΩs (a0 3 ),N(Ωs) (me), ρs(ave) (e a0 −3 ), and Hs(ave/max/min) (hartree a0 −3 ) are, respectively, the volume enclosed by the s(r) = 0.4 isosurface at around the BCP, and the total electronic charge, the average electron density, and the average, maximum and minimum value of H(r) over Ωs.  Based on the data quoted in Table 1, the difference between PA(Ng) and aPA(Ng), 17.0 kcal mol −1 for Ng = He, 20.5 kcal mol −1 for Ng = Ne, 45.8 kcal mol −1 for Ng = Ar, 51.9 kcal mol −1 for Ng = Kr, and 59.8 kcal mol −1 for Ng = Xe, unravel, in particular, a thermochemical weakening effect that progressively increases when going from He to Xe. In any case, the bonds of any (Ng-H-Ng) + are still assigned as Cov. As shown in Figure 2 and Table 3, the plotted H(r) is, invariably, of type A, and the ρ(BCP) is, invariably, higher than 0.08 ea 0 −3 . The BD is, however, lower, than that of the corresponding NgH + , with percentage decreases of ca. 23% for Ng = Ne, and ca. 33-40% for Ng = He, Ar, Kr, and Xe. In this regard, it is of interest to note that, when going from H 2 (R = 0.743 Å) to the linear centrosymmetric (H-H-H) − (R = 1.059 Å), isoelectronic with (He-H-He) + , the BD of the H-H bond decreases as well by ca. 37% (1.150 vs. 0.722 hartree e −1 ). However, at variance with (He-H-He) + , the covalent H 3 − is a first-order saddle point on the potential energy surface, connecting two equivalent linear van der Waals energy minima H − (H 2 ), and less stable than H − + H 2 by ca. 10 kcal mol −1 [92][93][94]. Major differences referable to the positive vs. negative charge of the two ions.

ρ(BCP) ρ(HCP) H(HCP) BD
To summarize, when going from NgH + to (Ng-H-Ng) + , the strength of the Ng-H interaction decreases, but its character does not change. In the heteronuclear (Ng-H-Ng') + , the mutual effects of Ng/Ng' on the adjacent Ng'-H/Ng-H bond are, instead, more complex and variegated. This is best discussed in the subsequent paragraph.   The structure, stability, and bonding character of the heteronuclear (Ng-H-Ng') + depend, essentially, on the difference between the PA of Ng and Ng'. The PAs of He and Ne are comparable (47.1 vs. 52.8 kcal mol −1 , see Table 1) and, in fact, in the (He-H-Ne) + the (formal) H + is, essentially, shared between the two atoms. The He-H and Ne-H distances are comparable (0.9589 vs. 1.1082 Å), the 2B(Ne) channel is only slightly more endothermic than the 2B(He) (17.6 vs. 11.9 kcal mol −1 ), and the endothermicity of the 3B channel, 64.7 kcal mol −1 , is definitely lower than the sum of the PA of He and Ne, 99.9 kcal mol −1 , thus confirming a mutual weakening effect. The values of ∆R(He-H) and ∆R(Ne-H) are also comparable, and predicted as ca. 0.18 and ca. 0.12 Å, respectively. Consistently, the bonding situation of the (He-H-Ne) + is similar to that occurring in the homonuclear (Ng-H-Ng) + , and accounted by the two nearly-equivalent resonance structures (He-H + )(Ne) and (He)(H-Ne + ). Both the He-H and Ne-H bonds are, in fact, assigned as Cov, with a BD of 0.650 and 0.898, respectively. When going from He and Ne to Ar, the PA appreciably increases up to 93.6 kcal mol −1 . Consistently, (He-H-Ar) + and (Ne-H-Ar) + are best described, respectively, by the resonance structures (He)(H-Ar + ) and (Ne)(H-Ar + ), accounting for weakly-bound complexes of He and Ne with the covalent Ar-H + . The He-H/Ne-H distances are, in fact, as long as 1.5157/1.5790 Å, the 2B(He)/2B(Ne) channel is endothermic by only 2.08/3.83 kcal mol −1 , and the ν 3 stretching frequencies of 2569 and 2399 cm −1 , respectively, are lower than, but comparable with the stretching absorption of ArH + , 2729 cm −1 . These structural assignments anticipate the results of the bonding analysis. Thus, as shown in Figure 2, the plotted H(r) of both (He-H-Ar) + and (Ne-H-Ar) + clearly shows the H − (r) regions of Ng and ArH + that are just touching (Ng = Ne) or even separate(Ng = He). The adjacent He or Ne only slightly perturb the Ar-H bond, actually assigned as Cov, with a BD of 0.964 hartree e −1 (Ng = He) and 0.944 hartree e −1 (Ng = Ne) that is only slightly lower than the BD of ArH + , 0.976 hartree e −1 (see Table 3). The He-H and Ne-H bonds are, instead, definitely weaker than those occurring in the HeH + and NeH + , even though they still maintain a contribution of covalency. They are, in fact, assigned as pCov(C/H − ) and pCov(A/H − ), respectively, with corresponding ρ s (ave) of 0.0249 and 0.0383 e a 0 −3 , and H s (ave) of −0.0024 and −0.0062 hartree a 0 −3 , respectively. The PA of Kr (105.0 kcal mol −1 ) and Xe (120.2 kcal mol −1 ) is definitely higher than that of He and Ne, and in fact, in the four complexes (Ng-H-Kr) + and (Ng-H-Xe) + (Ng = He, Ne), the (formal) H + is definitely closer to Kr or Xe, with Kr-H or Xe-H distances that are only slightly elongated with respect to the diatomic KrH + or XeH + [the ∆R(Kr-H) and ∆R(Xe-H) are, invariably, less than 0.01 Å, see Table 2]. Consistently, as shown in Table 1, the 2B(He) and 2B(Ne) channels are endothermic by only 0.6-2.4 kcal mol −1 , and the endothermicity of the 3B channels are, invariably, nearly coincident with the PA of Kr or Xe. The ν 3 stretching frequencies of (He-H-Kr) + , 2482 cm −1 , (Ne-H-Kr) + , 2408 cm −1 , (He-H-Xe) + , 2291 cm −1 , and (Ne-H-Xe) + , 2271 cm −1 , are also close to the stretching absorptions of KrH + , 2528 cm −1 , or XeH + , 2297 cm −1 . In essence, these systems should be considered as a diatomic ion separated from a nearly isolated Ng atom. Consistently, the plots of Figure 2 and the data quoted in Table 3 unravel a clearly distinguishable covalent KrH + or XeH + moiety, whose BD is coincident or quite close to that of the diatomic ions, and whose H − (r) region is well separated from that of He or Ne. The three (He-HKr) + , (He-HXe) + , and (Ne-HXe) + bonds are, indeed, all assigned as nCov, and their predicted N(Ω s ) (less than 1 me) and ρ s (ave) (less than 0.015 e a 0 −3 ) are typical of non-covalent contacts [77]. The computed DoP(He) and DoP(Ne) of 7.92, 4.49, and 3.22, respectively, are, however, rather large [76], and clearly signal He or Ne atoms that are appreciably polarized toward the adjacent cation. As a matter of fact, in the (Ne-H-Kr) + , the polarization of Ne (DoP = 4.53) is already sufficient to make the (Ne-HKr) + bond partially covalent [pCov(C/H +/− )], with a ρ s (ave) of 0.0228 e a 0 −3 . The ν 3 absorptions of (Ar-H-Kr) + , 1390 cm −1 , (Ar-H-Xe) + , 1849 cm −1 , and (Kr-H-Xe) + , 1410 cm −1 , are significantly lower than the stretching frequencies of ArH + , KrH + , and XeH + . This resembles the situation occurring in (He-H-Ne) + , and in the five homonuclear (Ng-H-Ng) + (Ng = He-Xe) (vide supra), and could suggest relatively similar bonding situations. As a matter of fact, the plotted H(r) of the (Ar-H-Kr) + , (Ar-H-Xe) + , and (Kr-H-Xe) + (see Figure 3) invariably shows the continuing overlapping of the H − (r) regions that is typical of covalent or partially-covalent bonds. The quantitative analysis, however, unravels mutual weakening effects of Ng and Ng' that are peculiar of the three ions. Thus, in the (Ar-H-Kr) + , the ∆R(Kr-H) amounts to only 0.1102 Å, and the Kr-H bond is assigned as Cov, even though its BD of 0.782 hartree e −1 is lower than the BD of KrH + , 0.937 hartree e −1 . The ∆R(Ar-H), 0.4129 Å, is, instead, large enough to produce a Ar-H bond assigned as pCov(A,H − ), even though still featuring a definitely negative H s (ave) of −0.0333 hartree a 0 −3 , and a relatively high ρ s (ave) of 0.0715 e a 0 −3 . In essence, the difference of 11.4 kcal mol −1 between the PA of Ar and Kr produces a resonance structure (Ar)(H-Kr + ) that is more weighting than, but not by far prevailing on the (Ar-H + )(Kr). When going from (Ar-H-Kr) + to (Ar-H-Xe) + , the ∆R(Ar-H) further increases up to 0.6701 Å, and the endothermicity of the 2B(Ar) channel decreases from 10.4 to 5.70 kcal mol −1 (see Table 1). Consistently, the Ar-H bond, still assigned as pCov(A/H − ), has a lower degree of covalency, signaled by an increased (less negative) H s (ave) of -0.0086 hartree a 0 −3 , and a dereased ρ s (ave) of 0.0387 e a 0 −3 . The resonance structure (Ar)(H-Xe + ) is, thus, by far prevailing on the (Ar-H + )(Xe). The difference between the PA of Kr and Xe, 15.2 kcal mol −1 , is only slightly larger than that between the PA of Ar and Kr, and the bonding situation of (Kr-H-Xe) + is, indeed, similar to that occurring in ( Overall, the structure, stability, and bonding character of the fifteen proton-bound (Ng-H-Ng) + and (Ng-H-Ng') + are strictly related, their qualitative and also quantitative trends mirroring the periodically-variable mutual effects of Ng and Ng'. A summarizing overview of these relationships is given in the subsequent paragraph.

4.
3. An Overview of the (Ng-H-Ng') + (Ng = Ng' or Ng = Ng') As already noted previously [63], an effective mode to overview the properties of the fifteen (Ng-H-Ng') + (Ng = Ng' or Ng = Ng') is to analyze the effects produced on the diatomic NgH + by the ligation with Ng'. The sharing of the (formal) H + between Ng and Ng' generally produces an elongation of the Ng-H bond, the effect depending on the resistance of Ng to donate the proton, and on the ability of Ng' to accept it. These competing tendencies both increase by increasing the PA of Ng and Ng', and, therefore, increase on going from He to Xe.Consistently, the values of ∆R(Ng-H) quoted in Table 2 clearly unravel i) the elongation of any Ng-H bond by any Ng' ii) for any Ng', a decreased elongation when going from He-H to Xe-H, and iii) for any Ng-H, an increased elongation when going from Ng' = He to Ng' = Xe. The quantitative effects of the nature of Ng and Ng', and, in particular, of the values of their PA are best caught by inspecting the graphs given in Figure 3, showing the trends of the ∆R(Ng-H) as a function of Ng'.
Thus, the bond distances of HeH + and NeH + , only slightly and comparably elongated by ligation with He or Ne (by ca. 0.12 and ca. 0.18 Å), suddenly increase upon ligation with Ar (by ca. 0.6-0.7 Å), and further increase nearly linearly by ligation with Kr (by ca. 0.8-1.0 Å), and Xe (by ca. 1.0 and 1.3 Å). These structural effects clearly mirror the changes occurring in the thermochemical stability and bonding character of the (He-H-Ng') + and (Ne-H-Ng') + . Thus, the endothermicty of the 2B(He) and 2B(Ne) channels, ranging between ca. 12 and ca. 18 kcal mol −1 for Ng' = He or Ne, drastically reduces to ca. 2-4 kcal mol −1 for Ng' = Ar, and arrives down to ca. 1 kcal mol −1 or even less for Ng' = Kr or Xe. Meanwhile, as schematically shown in Figure 4, when going from Ng' = He or Ne to Ng' = Ar, the character of the (He-HNg') + and (Ne-HNg') + bonds features a major shift from the Cov to the pCov domain, and arrives up to the nCov for (He-HKr') + , (He-HXe') + , and (Ne-HXe') + .
As shown in Figure 3, due to their lowest PA, both He and Ne are, essentially, unable to elongate the bond distances of ArH + , KrH + , and XeH + . Appreciable effects on these ions are, instead, exerted by Ar, and, especially, Kr and Xe. One also notes that, on going from Ng' = Ne to Ng' = Kr, the dependencies of ∆R(Ar-H), ∆R(Kr-H), and ∆R(Xe-H) are, essentially, linear but progressively less pending, and become slightly more pending on going to Ng' = Xe. These trends mirror, in essence, the progressively-increased reluctance of ArH + , KrH + , and XeH + to share the (formal) H + , and the highest ability of Xe to accept it. As a matter of fact, as shown in Figure 4, on going from the (Ne-H-Ng') + to the (Xe-H-Ng') + , the pCov domain of the Ng-H bond becomes progressively less populated, and all the five (Xe-HNg') + bonds are, actually, assigned as Cov. In any case, as also eye-caught by inspecting the right part of Figure 4, the intrinsic strength of the various covalent Ng-H bonds progressively decreases on going from the He-H to the Xe-H; the complementary aspect of the progressively increased tendency of the Ng atoms to share their valence electrons when going from He to Xe.  As shown in Figure 3, due to their lowest PA, both He and Ne are, essentially, unable to elongate the bond distances of ArH + , KrH + , and XeH + . Appreciable effects on these ions are, instead, exerted by Ar, and, especially, Kr and Xe. One also notes that, on going from Ng' = Ne to Ng' = Kr, the dependencies of ΔR(Ar-H), ΔR(Kr-H), and ΔR(Xe-H) are, essentially, linear but progressively less pending, and become slightly more pending on going to Ng' = Xe. These trends mirror, in essence, the progressively-increased reluctance of ArH + , KrH + , and XeH + to share the (formal) H + , and the highest ability of Xe to accept it. As a matter of fact, as shown in Figure 4, on going from the (Ne-H-Ng') + to the (Xe-H-Ng') + , the pCov domain of the Ng-H bond becomes progressively less populated, and all the five (Xe-HNg') + bonds are, actually, assigned as Cov. In any case, as also eye-caught by inspecting the right part of Figure 4, the intrinsic strength of the various covalent Ng-H bonds progressively decreases on going from the He-H to the Xe-H; the complementary aspect of the progressively increased tendency of the Ng atoms to share their valence electrons when going from He to Xe.

Conclusions
One of the most fascinating aspects of the chemistry of the noble gases is their ability to form chemical bonds of quite diverse type, ranging from weak non-covalent contacts to strong covalent bonds. Exemplary in this regard are the presently-investigated proton-

Conclusions
One of the most fascinating aspects of the chemistry of the noble gases is their ability to form chemical bonds of quite diverse type, ranging from weak non-covalent contacts to strong covalent bonds. Exemplary in this regard are the presently-investigated protonbound dimers (Ng-H-Ng) + and (Ng-H-Ng') + . The progressively-variable nature of the Ng-H and Ng'-H contacts parallels the concomitant variations of the bond distances and thermochemical stabilities, featuring as well regular trends as a function of Ng and Ng'. The structure, stability, and bonding character of these ionic complexes are, indeed, strictly related, and this suggests similar relationships for other groups of cationic noble-gas hydrides. Their further investigation is also stimulated by the current active interest for the conceivable role of various Ng m H n + (m, n ≥ 1) in processes of astrochemical interest.