Theoretical investigations of proton-bound cluster ions

Gas-Phase Infrared Action Spectroscopy of CH2Cl+ and CH3ClH+: Likely Protagonists in Chlorine Astrochemistry

Two fundamental halocarbon ions, CH2Cl+ and CH3ClH+, were studied in the gas phase using the FELion 22-pole ion trap apparatus and the Free Electron Laser for Infrared eXperiments (FELIX) at Radboud University, Nijmegen (the Netherlands). The vibrational bands of a total of four isotopologs, CH235,37Cl+ and CH335,37ClH+, were observed in selected wavenumber regions between 500 and 2900 cm-1 and then spectroscopically assigned based on the results of anharmonic force field calculations performed at the CCSD(T) level of theory. As the infrared photodissociation spectroscopy scheme employed probes singly Ne-tagged weakly bound complexes, complementary quantum-chemical calculations of selected species were also performed. The impact of tagging on the vibrational spectra of CH2Cl+ and CH3ClH+ is found to be virtually negligible for most bands; for CH3ClH+-Ne, the observations suggest a proton-bound structural arrangement. The experimental band positions as well as the best estimate rotational molecular parameters given in this work provide a solid basis for future spectroscopic studies at high spectral resolutions.

Structural and vibrational characterization of HCO+ and Rg-HCO+, Rg = {He, Ne, Ar, Kr, and Xe}

The structures of the formyl ion (HCO⁺) and its rare gas tagged counterparts (Rg-HCO⁺, Rg = He, Ne, Ar, Kr, and Xe) were studied at the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)]/aug-cc-pVTZ level of theory and basis set. A linear structure for these tagged complexes was predicted. The Rg binding energies for Rg-HCO⁺ are also examined at the CCSD(T) level. It was found that the binding interaction increases from He-HCO⁺ to Xe-HCO⁺. A multilevel potential energy surface built at the CCSD(T) and second-order Møller-Plesset perturbation levels of theory were used to study these species' vibrational spectra. By changing the Rg in the first-solvation shell for HCO⁺, the Fermi resonance interaction between the first H⁺ bend overtone and the asymmetric and symmetric H-C-O stretches can be modulated. This Fermi resonance modulation is demonstrated by examining a series of rare gas solvated HCO⁺.

Analysis of the Proton Transfer Bands in the Infrared Spectra of Linear N2H+···OC and N2D+···OC Complexes Using Electric Field-Driven Classical Trajectories

In this work, we describe ab initio calculations and assignment of infrared (IR) spectra of hydrogen-bonded ion-molecular complexes that involve a fluxional proton: the linear N₂H⁺···OC and N₂D⁺···OC complexes. Given the challenges of describing fluxional proton dynamics and especially its IR activity, we use electric field-driven classical trajectories, i.e., the driven molecular dynamics (DMD) method that was developed by us in recent years and for similar applications, in conjunction with high-level electronic structure theory. Namely, we present a modified and a numerically efficient implementation of DMD specifically for direct (or "on the fly") calculations, which we carry out at the MP2-F12/AVDZ level of theory for the potential energy surface (PES) and MP2/AVDZ for the dipole moment surfaces (DMSs). Detailed analysis of the PES, DMS, and the time-dependence of the first derivative of the DMS, referred to as the driving force, for the highly fluxional vibrations involving H⁺/D⁺ revealed that the strongly non-harmonic PES and non-linear DMS yield remarkably complex vibrational spectra. Interestingly, the classical trajectories reveal a doublet in the proton transfer part of the spectrum with the two peaks at 1800 and 1980 cm^-1. We find that their shared intensity is due to a Fermi-like resonance interaction, within the classical limit, of the H⁺ parallel stretch fundamental and an H⁺ perpendicular bending overtone. This doublet is also observed in the deuterated species at 1360 and 1460 cm^-1.

Structure and Vibrational Spectra of ArnH+ (n = 2-3)

The structure and vibrational spectra of protonated Ar clusters Ar_nH⁺ (n = 2-3) are studied using potential energy surfaces at the CCSD(T)/aug-cc-pVTZ level and basis set. Ar binding energies, as well as position isomerism in Ar₃H⁺, were investigated. In our previous work, the spectra of Ar₂H⁺ reveal a strong progression of combination bands, which involves the asymmetric Ar-H⁺ stretch with multiple quanta of the symmetric Ar-H⁺ stretch. In this work, insights on the origin of such progression were examined using an adiabatic model. In addition, contributions from mechanical and electrical anharmonicity on the progressions' intensities were also examined. Comparison of the calculated spectrum for the bare and Ar-tagged ions reveals that the reduction of the symmetry group, from D_∞h to either C_∞v or C_2v, results in a richer vibrational structure in the 500-1700 cm^-1 region. When compared with previously reported action spectra (D. C. McDonald III, D. T. Mauney, D. Leicht, J. H. Marks, J. A. Tan, J.-L. Kuo, and M. A. Duncan, J. Chem. Phys., 2016, 145, 231,101), it appears that the position isomers, because of the binding of the weakly bound Ar messenger, are needed to account for the additional bands in the infrared photodissociation spectrum for Ar₃H⁺. These findings demonstrate the active role of the messenger atom in relaxing some of the selection rules for the bare ion's vibrational transitions - resulting in an augmentation of the bands in the action spectrum.

Assignment of Infrared-Active Combination Bands in the Vibrational Spectra of Protonated Molecular Clusters Using Driven Classical Trajectories: Application to N4H+ and N4D. J Phys Chem A 2019;123:5613-5620. [PMID: 31244119 DOI: 10.1021/acs.jpca.9b04480]

We investigate the utility of the driven molecular dynamics (DMD) approach to complex molecular vibrations by applying it to linear clusters with several degenerate vibrational modes and infrared (IR) intense combination bands. Here, the prominent features in N₄H⁺ and N₄D⁺ IR spectra, reported and described by others previously, have been characterized for the first time by DMD using recently published high-level potential and dipole moment surfaces. Namely, the calculations closely correlate the parallel proton stretch vibration in N₄H⁺, at 750 cm^-1, with the one observed experimentally at 743 cm^-1. Second, the intense IR-active combination bands found in experimental spectra within 900-1100 cm^-1 have been properly recovered by DMD at 950 cm^-1 as strongly IR-active and confirmed as consisting of H⁺ asymmetric stretch and N₂···N₂ intermolecular symmetric stretch modes. Furthermore, we show that certain combination bands involving overtone transitions may be recovered by DMD using a hard-driving regime, such as the 1409 cm^-1 band measured in N₄H⁺, revealed by DMD at 1375 cm^-1, and assigned to a progressive combination of the parallel H⁺ stretch and two quanta of N₂···N₂ stretch, in agreement with quantum mechanical studies reported previously by others.

IR Spectroscopy of Protonated Acetylacetone and Its Water Clusters: Enol-Keto Tautomers and Ion→Solvent Proton Transfer

Protonated ions of acetylacetone, H⁺(Hacac), and their argon-tagged analogues are produced via a pulsed discharge and cooled in a supersonic expansion. These ions are mass analyzed, selected in a time-of-flight spectrometer, and studied with infrared laser photodissociation spectroscopy using the method of rare-gas atom tagging. Computational studies at the DFT/B3LYP level are employed to elucidate the structures and spectra of these ions, which are expected to exist as either enol- or keto-based tautomers. The protonated acetylacetone ion is found to form a single enol-based isomer. Adding one or two water molecules to this ion, for example, H⁺(Hacac)(H₂O)_1,2, produces primarily enol-based structures, although a small concentration of keto structures also contribute to the spectra. The vibrational patterns resulting from hydrogen bonding in these systems are not well-described by theory. Addition of a third water molecule to form the H⁺(Hacac)(H₂O)₃ ion causes a significant change in the spectroscopy, attributed to proton transfer from the H⁺(Hacac) ion into the water solvent.

Communication: Trapping a proton in argon: Spectroscopy and theory of the proton-bound argon dimer and its solvation

Ion-molecule complexes of the form H⁺Ar_n are produced in pulsed-discharge supersonic expansions containing hydrogen and argon. These ions are analyzed and mass-selected in a reflectron spectrometer and studied with infrared laser photodissociation spectroscopy. Infrared spectra for the n = 3-7 complexes are characterized by a series of strong bands in the 900-2200 cm^-1 region. Computational studies at the MP2/aug-cc-pVTZ level examine the structures, binding energies, and infrared spectra for these systems. The core ion responsible for the infrared bands is the proton-bound argon dimer, Ar-H⁺-Ar, which is progressively solvated by the excess argon. Anharmonic vibrational theory is able to reproduce the vibrational structure, identifying it as arising from the asymmetric proton stretch in combination with multiple quanta of the symmetric argon stretch. Successive addition of argon shifts the proton vibration to lower frequencies, as the charge is delocalized over more ligands. The Ar-H⁺-Ar core ion has a first solvation sphere of five argons.

Infrared spectra and anharmonic coupling of proton-bound nitrogen dimers N2–H+–N2, N2–D+–N2, and 15N2–H+–15N2 in solid para-hydrogen

The proton-bound nitrogen dimer, N₂–H⁺–N₂, and its isotopologues were investigated by means of vibrational spectroscopy.

A theoretical investigation of the energetics and spectroscopic properties of the gas-phase linear proton-bound cation-molecule complexes, XCH(+)-N2 (X = O, S)

The structural features, spectroscopic properties, and interaction energies of the linear proton-bound complexes of OCH(+) and its sulfur analog SCH(+) with N2 were investigated using the high-level ab initio methods MP2 and CCSD(T) as well as density functional theory with the aug-cc-pVXZ (X = D, T) basis sets. The rotational constants along with the vibrational frequencies of the cation-molecule complexes are reported here. A comparison of the interaction energies of the OCH(+)-N2 and SCH(+)-N2 complexes with those of the OCH(+)-CO and OCH(+)-OC complexes was also performed. The energies of all the complexes were determined at the complete basis set (CBS) limit. CS shows higher proton affinity at the C site than CO does, so the complex OCH(+)-N2 is relatively strongly bound and has a higher interaction energy than the SCH(+)-N2 complex. Symmetry-adapted perturbation theory (SAPT) was used to decompose the total interaction energies of the complexes into the attractive electrostatic interaction energy (E elst), induction energy (E ind), dispersion energy (E disp), and repulsive exchange energy (E exch). We found that the ratio of E ind to E disp is large for these linear proton-bound complexes, meaning that inductive effects are favored in these complexes. The bonding characteristics of the linear complexes were elucidated using natural bond orbital (NBO) theory. NBO analysis showed that the attractive interaction is caused by NBO charge transfer from the lone pair on N to the σ*(C-H) antibonding orbital in XCH(+)-N2 (X = O, S). The quantum theory of atoms in molecules (QTAIM) was used to analyze the strengths of the various bonds within and between the cation and molecule in each of these proton-bound complexes in terms of the electron density at bond critical points (BCP). Graphical Abstract Linear proton-bound complexes of OCH(+)-N2 and SCH(+)-N2. In these complexes, inductive effect is favored over dispersive effect. The attractive interaction is the NBO charge transfer from N-lone pair of N2 to CH σ* antibonding orbital of XCH(+) (X = O, S).

Infrared Spectra and ab initio Calculations for Fluoride-acetylene Clusters: F - -(HCCH)n, n=3 - 6

Infrared spectra and ab initio calculations are presented for gas phase clusters consisting of a fluoride anion attached to acetylene ligands. Spectra obtained in the C-H stretch region contain a single strong band, consistent with cluster structures in which roughly equivalent acetylene ligands are hydrogen bonded to a central fluoride anion core. Minimum energy structures predicted from ab initio calculations at the MP2 level of theory are highly symmetric with acetylene ligands equally spaced about a central fluoride anion core. The predicted H-bonded C-H stretching frequencies, after scaling to correct for anharmonicity, agree well with the experimental band positions.

Infrared spectra of HC2−∙(C2H2)n and O2−∙(C2H2)n clusters (n=2–5)

We study the solvation of HC2- and O2- with acetylene ligands by means of midinfrared photodissociation spectroscopy in the CH stretching region, monitoring C2H2 evaporation upon infrared photon absorption by the parent cluster ions. Our findings are interpreted with the help of density functional theory. The infrared spectra indicate that while the binding generally occurs through ionic H bonds, there are two different classes of ligands which differ in their binding strength. This holds true for both core ions, even though their electronic structures and charge distributions are very different.

Structures and infrared spectra of fluoride–hydrogen sulfide clusters from ab initio calculations: F–-(H2S)n, n= 1–5

Clusters formed between a fluoride anion and several hydrogen sulfide molecules have been investigated via ab initio calculations at the MP2 level of theory, using Dunning's augmented correlation consistent basis sets. Optimised geometries, vibrational frequencies, and enthalpy changes for the ligand association reactions are presented for clusters with up to five H2S ligands interacting with a F- anion. The minimum energy structure for the 1:1 F(-)-H2S complex features proton transfer from the H2S to the F- anion, forming a planar C(s) symmetry FH...SH- structure. For the F(-)-(H2S)2 cluster, the FH...SH- core remains and is solvated by a perturbed H2S ligand. For the larger F(-)-(H2S)(3-5) clusters, in addition to the FH...SH(-)-(H2S)n cluster forms, other minima featuring a 'solvated F-' anion are predicted. Calculated infrared spectra for the minima of each cluster size are presented to aid in assigning spectra from future experimental studies.

Intermolecular interaction in an open-shell π-bound cationic complex: IR spectrum and coupled cluster calculations for C2H2+-Ar

The intermolecular potential energy surface (PES) of Ar interacting with the acetylene cation in its (2)Pi(u) ground electronic state is characterized by infrared photodissociation (IRPD) spectroscopy and quantum chemical calculations. In agreement with the theoretical predictions, the rovibrational analysis of the IRPD spectrum of C(2)H(2) (+)-Ar recorded in the vicinity of the antisymmetric CH stretching fundamental (nu(3)) is consistent with a vibrationally averaged T-shaped structure and a ground-state center-of-mass separation of R(c.m.) = 2.86 +/- 0.09 A. The nu(3) band experiences a blueshift of 16.7 cm(-1) upon complexation, indicating that vibrational excitation slightly reduces the interaction strength. The two-dimensional intermolecular PES of C(2)H(2) (+)-Ar, obtained from coupled cluster calculations with a large basis set, features strong angular-radial coupling and supports in addition to a global pi-bound minimum also two shallow side wells with linear H-bound geometries. Bound state rovibrational energy level calculations are carried out for rotational angular momentum J = 0-10 (both parities) employing a discrete variable representation-distributed Gaussian basis method. Effective spectroscopic constants are determined for the vibrational ground state by fitting the calculated rotational energies to the standard Watson A-type Hamiltonian for a slightly asymmetric prolate top.