Non-covalent interactions between metal cations and molecular hydrogen: spectroscopic studies of M+–H2complexes

Solvation of cationic copper clusters in molecular hydrogen

Multiply charged superfluid helium nanodroplets are utilized to facilitate the growth of cationic copper clusters (Cun+, where n = 1-8) that are subsequently solvated with up to 50 H2 molecules. Production of both pristine and protonated cationic Cu clusters are detected mass spectrometrically. A joint effort between experiment and theory allows us to understand the nature of the interactions determining the bonding between pristine and protonated Cu+ and Cu2+ cations and molecular hydrogen. The analysis reveals that in all investigated cationic clusters, the primary solvation shell predominantly exhibits a covalent bonding character, which gradually decreases in strength, while for the subsequent shells an exclusive non-covalent behaviour is found. Interestingly, the calculated evaporation energies associated with the first solvation shell markedly surpass thermal values, positioning them within the desirable range for hydrogen storage applications. This comprehensive study not only provides insights into the solvation of pristine and protonated cationic Cu clusters but also sheds light on their unique bonding properties.

Persistence of a Delocalized Radical in Larger Clusters of Hydrated Copper(II) Hydroxide, CuOH+(H2O)3-7

The structures, vibrational spectra, and electronic properties of copper hydroxide hydrates CuOH+(H2O)3-7 were investigated with quantum chemistry computations. As a follow-up to a previous analysis of CuOH+(H2O)0-2, this investigation examined the progression as the square-planar metal coordination environment was filled and as solvation shells expanded. Four-, five-, and six-coordinate structures were found to be low-energy isomers. The delocalized radical character, which was discovered in the small clusters, was found to persist upon continued hydration, although the hydrogen-bonded water network in the larger clusters was found to play a more significant role in accommodating this spin. Partial charges indicated that the electronic structure includes more Cu2+···OH- character than was observed in smaller clusters, but this structure remains decidedly mixed with Cu+···OH· configurations and yields roughly half-oxidation of the water network in the absence of any electrochemical potential. Computed vibrational spectra for n = 3 showed congruence with spectra from recent predissociation spectroscopy experiments, provided that the role of the D2 tag was taken into account. Spectra for n = 4-7 were predicted to exhibit features that are reflective of both the mixed electronic character and proton-/hydrogen-shuttling motifs within the hydrogen-bonded water network.

Probing the binding and activation of small molecules by gas-phase transition metal clusters via IR spectroscopy

Isolated transition metal clusters have been established as useful models for extended metal surfaces or deposited metal particles, to improve the understanding of their surface chemistry and of catalytic reactions. For this objective, an important milestone has been the development of experimental methods for the size-specific structural characterization of clusters and cluster complexes in the gas phase. This review focusses on the characterization of molecular ligands, their binding and activation by small transition metal clusters, using cluster-size specific infrared action spectroscopy. A comprehensive overview and a critical discussion of the experimental data available to date is provided, reaching from the initial results obtained using line-tuneable CO₂ lasers to present-day studies applying infrared free electron lasers as well as other intense and broadly tuneable IR laser sources.

Vibrational spectroscopy of Cu+(H2)4: about anharmonicity and fluxionality

The vibrational spectra of the copper(I) cation-dihydrogen complexes Cu⁺(H₂)₄, Cu⁺(D₂)₄ and Cu⁺(D₂)₃H₂ are studied using cryogenic ion trap vibrational spectroscopy in combination with quantum chemical calculations. The infrared photodissociation (IRPD) spectra (2500-7300 cm^-1) are assigned based on a comparison to IR spectra calculated using vibrational second-order perturbation theory (VPT2). The IRPD spectra exhibit ≈60 cm^-1 broad bands that lack rotational resolution, indicative of rather floppy complexes even at an ion trap temperature of 10 K. The observed vibrational features are assigned to the excitations of dihydrogen stretching fundamentals, combination bands of these fundamentals with low energy excitations as well as overtone excitations of a minimum-energy structure with C_s symmetry. The three distinct dihydrogen positions present in the structure can interconvert via pseudorotations with energy barriers less than 10 cm^-1, far below the zero-point vibrational energy. Ab initio Born-Oppenheimer molecular dynamics (BOMD) simulations confirm the fluxional behavior of these complexes and yield an upper limit for the timeframe of the pseudorotation on the order of 10 ps. For Cu⁺(D₂)₃H₂, the H₂ and D₂ loss channels yield different IRPD spectra indicating non-ergodic behavior.

Electronic Structure and Vibrational Signatures of the Delocalized Radical in Hydrated Clusters of Copper("II") Hydroxide CuOH+(H2O)0-2

The copper hydroxide ion, CuOH⁺, serves as the catalytic core in several recently developed water-splitting catalysts, and an understanding of its chemistry is critical to determining viable catalytic mechanisms. In spite of its importance, the electronic structure of this open-shell ion has remained ambiguous in the literature. In particular, computed values for both the thermodynamics of hydration and the vibrational signatures of the mono- and dihydrates have shown prohibitively large errors compared to values from recent experimental measurements. In this work, the source of this discrepancy is demonstrated to be the propensity of this ion to exist between traditional Cu(I) and Cu(II) oxidation-state limits. The spin density of the radical is accordingly shown to delocalize between the metal center and surrounding ligands, and increasing the hydration serves to exacerbate this behavior. Equation-of-motion coupled-cluster methods demonstrated the requisite accuracy to resolve the thermodynamic discrepancies. Such methods were also needed for spectral simulations, although the latter also required a direct simulation of the role of the deuterium "tag" molecules that are used in modern predissociation spectroscopy experiments. This nominally benign tag molecule underwent direct complexation with the open-valence metal ion, thereby forming a species akin to known metal-H₂ complexes and strongly impacting the resulting spectrum. Thermal populations of this configuration and other more traditional noncovalently bound isomers led to a considerable broadening of the spectral lineshapes. Therefore, at least for the CuOH⁺(H₂O)_0-2 hydrates, these benchmark ions should be considered to be delocalized radical systems with some degree of multireference character at equilibrium. They also serve as a cautionary tale for the spectroscopy community, wherein the role of the D₂ tag is far from benign.

Activation of D2 by Neodymium Cation (Nd+): Bond Energy of NdH+ and Mechanistic Insights through Experimental and Theoretical Studies

The kinetic-energy-dependent cross section for the reaction of Nd⁺ with D₂ was studied by using a guided ion beam tandem mass spectrometer. The formation of NdD⁺ is endothermic, and analysis of the reaction cross section gave an NdH⁺ 0 K bond dissociation energy (BDE) of 1.99 ± 0.06 eV. Theoretical calculations for the NdH⁺ BDE were performed for comparison with the experimental thermochemistry and generally gave accurate results. Additionally, relaxed potential energy surfaces for NdH₂⁺ were performed, and no strongly bound dihydride intermediate was located. The Nd⁺ + D₂ reactivity and BDE of NdH⁺ are compared with analogous results for the lanthanide cations La⁺, Ce⁺, Pr⁺, Sm⁺, Gd⁺, and Lu⁺ to establish periodic trends and insight into the role of the electronic configurations on this reactivity and the lanthanide hydride cation bond strengths.

Cerium Cation (Ce+) Reactions with H2, D2, and HD: CeH+ Bond Energy and Mechanistic Insights from Guided Ion Beam and Theoretical Studies

Reactions of the atomic lanthanide cerium cation (Ce⁺) with H₂, D₂, and HD were studied by using guided ion beam tandem mass spectrometry. Analysis of the kinetic-energy-dependent endothermic reactions to form CeH⁺ (CeD⁺) led to a 0 K bond dissociation energy (BDE) for CeH⁺ of 2.19 ± 0.09 eV. Theoretical calculations for CeH⁺ were performed at the B3LYP, BHLYP, and PBE0 levels of theory and overestimate the experimental BDE. In contrast, extrapolation to the complete basis set limit using coupled-cluster with single, double, and perturbative triple excitations, CCSD(T), gave a value (2.33 eV) in reasonable agreement with the experimental BDE. The branching ratio of the CeH⁺ and CeD⁺ products in the HD reaction suggests that the reaction occurs via a statistical mechanism involving a long-lived intermediate. Relaxed potential energy surfaces for CeH₂⁺ were computed and are consistent with the availability of such an intermediate, but the crossing point between quartet and doublet surfaces helps explain the inefficiency of the association reaction observed in the literature. The reactivity and CeH⁺ BDE are compared with previous results for group 4 transition metal cations (Ti⁺, Zr⁺, and Hf⁺), other lanthanides (La⁺, Sm⁺, Gd⁺, and Lu⁺), and the isovalent actinide Th⁺. Periodic trends and insight into the role of the electronic configuration on metal-hydride bond strength are discussed.

Complexes of Alkali Metal Cations and Molecular Hydrogen: Potential Energy Surfaces and Bound States

Complexes between metal cations and molecular hydrogen are systems quite amenable for precise spectroscopic and theoretical studies, and at the same time, they are relevant for applications in hydrogen storage and astrochemistry. In this work, we report new intermolecular potential energy surfaces and rovibrational states calculations for complexes involving molecular hydrogen and alkaline metal cations, M⁺-H₂ (M⁺ = Na⁺, K⁺, Rb⁺, Cs⁺). The intermolecular potentials, formulated in an internally consistent way to emphasize differences in the properties of the systems, are represented by simple analytical expressions whose parameters have been optimized from comparison with accurate ab initio calculations. Properties of the low-lying bound states-binding energies, frequencies, and rotational constants-are compared with previous measurements or computations and an overall good agreement is achieved, supporting the reliability of the present formulation. Variations of these properties as a function of the cation size and isotopic substitution, with a proper sequence of ortho and para rotational levels, are also discussed.

Snowball formation for Cs+ solvation in molecular hydrogen and deuterium

Interactions of atomic cations with molecular hydrogen are of interest for a wide range of applications in hydrogen technologies. These interactions are fairly strong despite being non-covalent, hence one can ask whether hydrogen molecules would form dense, solid-like, solvation shells around the ion (snowballs) or rather a more weakly bound compound. In this work, the interactions between Cs+ and H2 are studied both experimentally and computationally. Isotopic substitution of H2 by D2 is also investigated. On the one hand, helium nanodroplets doped with cesium and hydrogen or deuterium are ionized by electron impact and the (H2/D2)nCs+ (up to n = 30) clusters formed are identified via mass spectrometry. On the other hand, a new analytical potential energy surface, based on ab initio calculations, is developed and used to study cluster energies and structures by means of classical and quantum-mechanical Monte Carlo methods. The most salient features of the measured ion abundances are remarkably mimicked by the computed evaporation energies, particularly for the clusters composed of deuterium. This result supports the reliability of the present potential energy surface and allows us to recommend its use in related systems. Clusters with either twelve H2 or D2 molecules stand out for their stability and quasi-rigid icosahedral structures. However, the first solvation shell involves thirteen or fourteen molecules for hydrogenated or deuterated clusters, respectively. This shell retains its internal structure when extra molecules are added to the second shell and is nearly solid-like, especially for the deuterated clusters. The role played by three-body induction interactions as well as the rotational degrees of freedom is analyzed and they are found to be significant (up to 15% and 18%, respectively) for the molecules belonging to the first solvation shell.

Samarium cation (Sm+) reactions with H2, D2, and HD: SmH+ bond energy and mechanistic insights from guided ion beam and theoretical studies

Guided ion beam tandem mass spectrometry is used to study the reaction of the lanthanide samarium cation (Sm⁺) with H₂ and its isotopologues (HD and D₂) as a function of collision energy. Modeling the resulting energy dependent product ion cross sections from these endothermic reactions yields 2.03 ± 0.06 eV (two standard deviations) for the 0 K bond dissociation energy of SmH⁺. Quantum chemical calculations are performed to determine stabilities of the ground and low-energy states of SmH⁺ for comparison with the experimentally measured thermochemistry. The calculations generally overestimate the SmH⁺ bond energy, but a better agreement between experiment and theory is achieved after correcting for spin-orbit energy contributions, with coupled-cluster with single, double and perturbative triple excitations/complete basis set [CCSD(T)/CBS] results reproducing the experiment well. In the HD reaction, the SmH⁺ product is observed to be favored over the SmD⁺ by about a factor of three, indicating that the reaction proceeds via a direct mechanism with short-lived intermediates. This is consistent with quantum chemical calculations of relaxed potential energy surface scans of SmH₂ ⁺, which show that there is no strongly bound dihydride intermediate. The reactivity and hydride bond energy of Sm⁺, which has a valence electron configuration typical of most lanthanides, are compared with previous results for the lanthanide cations La⁺, Gd⁺, and Lu⁺, which exhibit configurations more closely related to the group 3 metal cations, Sc⁺ and Y⁺. Periodic trends across the lanthanide series and insights into the role of the electronic configurations on hydride bond strength and reactivity with H₂ are discussed.

Global diabatic potential energy surfaces for the BeH2 + system and dynamics studies on the Be+(2P) + H2(X1Σg +) → BeH+(X1Σ+) + H(2S) reaction

The Be+(2P) + H2(X1Σg +) → BeH+(X1Σ+) + H(2S) reaction has great significance for studying diabatic processes and ultracold chemistry. The first global diabatic potential energy surfaces (PESs) which are correlated with the lowest two adiabatic states 12A' and 22A' of the BeH2 + system are constructed by using the neural network method. Ab initio energy points are calculated using the multi-reference configuration interaction method with the Davidson correction and AVQZ basis set. The diabatic energies are obtained from the transformation of ab initio data based on the dipole moment operators. The topographical characteristics of the diabatic PESs are described in detail, and the positions of crossing between the V d 11 and V d 22 are pinpointed. On new diabatic PESs, the time-dependent quantum wave packet method is carried out to study the mechanism of the title reaction. The results of dynamics calculations indicate the reaction has no threshold and the product BeH+ is excited to high vibrational states easily. In addition, the product BeH+ tends to backward scattering at most collision energies.

Activation of H2 by Gadolinium Cation (Gd+): Bond Energy of GdH+ and Mechanistic Insights from Guided Ion Beam and Theoretical Studies

The energy-dependent reactions of the lanthanide gadolinium cation (Gd⁺) with H₂, D₂, and HD are investigated using guided ion beam tandem mass spectrometry. From analysis of the resulting endothermic product ion cross sections, the 0 K bond dissociation energy for GdH⁺ is measured to be 2.18 ± 0.07 eV. Theoretical calculations on GdH⁺ are performed for comparison with the experimental thermochemistry and generally appear to overestimate the experimental GdH⁺ bond dissociation energy. The branching ratio of the products in the HD reaction suggests that Gd⁺ reacts primarily via a statistical insertion mechanism to form the hydride product ion with contributions from direct mechanisms. Relaxed potential energy surfaces for GdH₂⁺ are computed and are consistent with the availability of both statistical and direct reaction pathways. The reactivity and hydride bond energy for Gd⁺ is compared with previous results for the group three metal cations, Sc⁺ and Y⁺, and the lanthanides, La⁺ and Lu⁺, and periodic trends are discussed.

Sumanene: an efficient π-bowl for dihydrogen storage

The outstanding ability of sumanene, its anionic and dianionic forms and sumanene–M⁺ion-pair complexes (M = Li, Na, K) to bind dihydrogen has been revealed using density functional theory calculations pointing out that these systems could be employed for developing new H₂storage systems.

Ab Initio Characterization of the Electrostatic Complexes Formed by H2 Molecule and Cr(+), Mn(+), Cu(+), and Zn(+) Cations

Equilibrium structures, dissociation energies, and rovibrational energy levels of the electrostatic complexes formed by molecular hydrogen and first-row S-state transition metal cations Cr(+), Mn(+), Cu(+), and Zn(+) are investigated ab initio. Extensive testing of the CCSD(T)-based approaches for equilibrium structures provides an optimal scheme for the potential energy surface calculations. These surfaces are calculated in two dimensions by keeping the H-H internuclear distance fixed at its equilibrium value in the complex. Subsequent variational calculations of the rovibrational energy levels permits direct comparison with data obtained from equilibrium thermochemical and spectroscopic measurements. Overall accuracy within 2-3% is achieved. Theoretical results are used to examine trends in hydrogen activation, vibrational anharmonicity, and rotational structure along the sequence of four electrostatic complexes covering the range from a relatively floppy van der Waals system (Mn(+)···H2) to an almost a rigid molecular ion (Cu(+)···H2).

The remarkable ability of anions to bind dihydrogen

Anions show a noteworthy ability to bind with a large number of hydrogen molecules which can be utilized for the development of novel salt systems for hydrogen storage.

Nuclear motion in the σ-bound regime of metal-H₂ complexes: [Mg(H₂)(n=1-6)]²⁺

The dynamic, quantum structure of [Mg(H2)(n=1-6)](2+)complexes is investigated via ab initio path integral molecular dynamics simulations. These complexes represent the strong, σ-complex regime of metal-H2 interactions and are representative of bonding motifs found in metal-organic frameworks. Significant nuclear motion within the coordination sphere is observed, even though the ligands remain largely intact. Quantum effects are found to be important in the H-H and metal-H2 stretch coordinates, but the remaining motion in the molecule is well represented by classical simulations. Nearly free rotation of the dihydrogen moiety is observed in all complexes. Statistical averages and distributions of structural parameters are found to deviate nontrivially from the same parameters in static, equilibrium structures.