CH5+: the infrared spectrum observed

Ion molecule reactions in the HBr+ + CH4 system: a combined experimental and theoretical study

Reactions in the system HBr+ + CH4 have been investigated inside a guided ion-beam apparatus under single-collision conditions. The HBr+ is vibrational and rotational state selected in the electronic X2Π1/2 state created by (2+1)-REMPI. Due to the exitation scheme employed different rotational states of the HBr+ are accessible. Four reaction channels have been observed. The cross section, σ, for the exothermic proton transfer channel (PT) decreases with increasing collision energy, steeper than predicted by the Langevin model. The cross section also decreases with increasing rotational energy in the HBr+, with the effect of the rotational energy being stronger than that of translational energy. The cross section for the endothermic charge transfer (CT) increased with increasing collision energy. The energy dependence is well reproduced by a simple line of center (loc) model. Although the bromine transfer (BT) is exothermic the observed cross section increased with increasing collision energy due to an activation barrier on the potential energy surface (PES). Analysis by a modified loc model suggest the relevance of an angle dependence of σ. The cross section for the endothermic hydrogen atom abstraction (HA) exhibits a maximum at 2 eV Ecm. The measured cross sections are rationalized by means of reaction dynamics simulations which show good agreement with the experimental cross sections. The dynamics simulations are carried out with a machine learning potential that is developed and benchmarked with ab initio molecular dynamics simulation. The absolute cross sections predicted by reaction dynamics simulations are well within the same order of magnitude while reproducing the trends over three different collision energies for all four reaction channels. Furthermore, the simulations demonstrate various reaction mechanisms for these reaction channels, including a very interesting HBr+ orientation selectivity for the BT reaction channel.

Quantum study of the CH3+ photodissociation in full-dimensional neural network potential energy surfaces

C H 3 + , a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by using the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction ab initio calculations are done, and full-dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full-dimensional quantum wave packet method in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids, allowing us to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion bar, result in a lesser destruction of CH3+ compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry.

Strong carbon - noble gas covalent bond and fluxionality in hypercoordinate compounds

Thermodynamic, kinetic, and chemical bonding analysis at the coupled cluster level has been carried out for a series of hypercoordinated carbon compounds with formula CH4Ng2+ (Ng = He-Rn). Results show that these compounds could be stable at room temperature and Born-Oppenheimer molecular dynamics simulations (BOMD) in conjunction with activation energies indicate high kinetic stability. In addition, all chemical bonding descriptors agree with a strong C-Ng covalent bond and a bonding pattern similar to that of CH5+. Finally, BOMD simulations showed that these compounds are fluxional, with a continuous formation/breaking of H-H and C-H bonds. To the best of the authors' knowledge, these results represent the first series of fluxional compounds possessing a covalent bond between a main group element and a noble gas atom.

Quantum Nuclear Delocalization and its Rovibrational Fingerprints

Quantum mechanics dictates that nuclei must undergo some delocalization. In this work, emergence of quantum nuclear delocalization and its rovibrational fingerprints are discussed for the case of the van der Waals complexHHe 3 + ${{\rm{HHe}}_3^ + }$ . The equilibrium structure ofHHe 3 + ${{\rm{HHe}}_3^ + }$ is planar and T-shaped, one He atom solvating the quasi-linear He-H+ -He core. The dynamical structure ofHHe 3 + ${{\rm{HHe}}_3^ + }$ , in all of its bound states, is fundamentally different. As revealed by spatial distribution functions and nuclear densities, during the vibrations of the molecule the solvating He is not restricted to be in the plane defined by the instantaneously bentHHe 2 + ${{\rm{HHe}}_2^ + }$ chomophore, but freely orbits the central proton, forming a three-dimensional torus around theHHe 2 + ${{\rm{HHe}}_2^ + }$ chromophore. This quantum delocalization is observed for all vibrational states, the type of vibrational excitation being reflected in the topology of the nodal surfaces in the nuclear densities, showing, for example, that intramolecular bending involves excitation along the circumference of the torus.

Isotope Effects on Ground and Excited States of Ethyl Cation, H+(C2H4)

The structure and spectra of ethyl cation, H⁺(C₂H₄), and its deuterated analogues are investigated using diffusion Monte Carlo (DMC). These calculations all show that the ground state wave function for H⁺(C₂H₄) is localized near the minimum energy configuration in which the excess proton is in a bridging configuration, although the amplitude of the vibrational motions of the bridging proton is large. Deuteration of the bridging proton reduces the amplitude of this motion, while deuteration of only the ethylenic hydrogen atoms in H⁺(C₂D₄) has little effect on the amplitude of the motion of the bridging proton. Excited states that are accessed by spectroscopically observed transitions in H⁺(C₂H₄) are calculated using fixed-node DMC. The calculated and measured frequencies for the states with one quantum of excitation in the ethylenic CH stretching vibrations show good agreement. We also explore the excited state with one quantum of excitation in the proton transfer vibration of the bridging proton and obtain a frequency of 616 cm^-1 for H⁺(C₂H₄). This frequency increases to 629 cm^-1 in H⁺(C₂D₄). Deuteration decreases this frequency to 491 and 495 cm^-1 in D⁺(C₂H₄) and D⁺(C₂D₄), respectively. The effects of partial deuteration on the frequencies of the CH stretching vibrations, and the corresponding probability amplitudes are also explored. Finally, we report the vibrationally averaged rotational constants for the four isotopologues of ethyl cation considered in this study.

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.

Quantum-Chemical and Quantum-Graph Models of the Dynamical Structure of CH5. J Chem Theory Comput 2023;19:42-50. [PMID: 36534596 DOI: 10.1021/acs.jctc.2c00991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Experimental and computational results about the structure, dynamics, and rovibrational spectra of protonated methane have challenged a considerable number of traditional chemical concepts. Hereby theoretical and computational results are provided about the dynamical structure of CH₅⁺. It is shown that the ground vibrational state investigated thus far by computations, forbidden by nuclear-spin statistics, has a structure similar to the first allowed vibrational state and, in fact, the structures of all vibrational states significantly below 200 cm^-1 are highly similar. Spatial delocalization of the nuclei, determined by nuclear densities computed from accurate variational vibrational wave functions, turns out to be limited when viewed in the body-fixed frame, confirming that the effective structure of CH₅⁺ is well described as a CH₃⁺ tripod with a H₂ unit on top of it. The interesting and unusual qualitative aspects of the sophisticated state-dependent variational results receive full explanation via simple quantum-graph models.

Pentacoordinate Carbon Atoms in a Ferrocene Dication Derivative—[Fe(Si2-η5-C5H2)2]2+

Pentacoordinate carbon atoms are theoretically predicted here in a ferrocene dication derivative in the eclipsed-(1; C2v), gauche-(2; C2) and staggered-[Fe(Si2-η5-C5H2)2]2+(3; C2h) forms for the first time. Energetically, the relative energy gaps for 2 and 3 range from −3.06 to 16.74 and −2.78 to 40.34 kJ mol−1, respectively, when compared to the singlet electronic state of 1 at different levels. The planar tetracoordinate carbon (ptC) atom in the ligand Si2C5H2 becomes a pentacoordinate carbon upon complexation. The ligand with a ptC atom was predicted to be both a thermodynamically and kinetically stable molecule by some of us in our earlier theoretical works. Natural bond orbital and adaptive natural density partitioning analyses confirm the pentacoordinate nature of carbon in these three complexes (1–3). Although they are hypothetical at the moment, they support the idea of “hypercoordinate metallocenes” within organometallic chemistry. Moreover, ab initio molecular dynamics simulations carried out at 298 K temperature for 2000 fs suggest that these molecules are kinetically stable.

Organic Chemistry and Synthesis Rely More and More upon Catalysts

A few months before the COVID-19 pandemic, Pierre Vogel and Kendall N. Houk published with a new textbook Wiley-VCH, “Organic Chemistry: Theory, Reactivity, and Mechanisms in Modern Synthesis”, with a foreword from the late Roberts H. Grubbs. The book demonstrates how catalytic processes dominate all fields of modern organic chemistry and synthesis, and how invention combines thermodynamics, kinetics, spectroscopy, quantum mechanics, and thermochemical data libraries. Here, the authors present a few case studies that should be of interest to teachers, practitioners of organic and organometallic chemistry, and the engineers of molecules. The Vogel–Houk book is both textbook and reference manual; it provides a modern way to think about chemical reactivity and a powerful toolbox to inventors of new reactions and new procedures.

Fast quasi-centroid molecular dynamics

We describe a fast implementation of the quasi-centroid molecular dynamics (QCMD) method in which the quasi-centroid potential of mean force is approximated as a separable correction to the classical interaction potential. This correction is obtained by first calculating quasi-centroid radial and angular distribution functions in a short path integral molecular dynamics simulation and then using iterative Boltzmann inversion to obtain an effective classical potential that reproduces these distribution functions in a classical NVT simulation. We illustrate this approach with example applications to the vibrational spectra of gas phase molecules, obtaining excellent agreement with QCMD reference calculations for water and ammonia and good agreement with the quantum mechanical vibrational spectrum of methane.

Infrared spectra of carbocations and CH4+ in helium

Infrared (IR) spectra of several hydrocarbon cations are reported, namely CH₃⁺, CH₄⁺, CH₅⁺, CH₅⁺(CH₄) and C₂H₅⁺. The spectra were generated from weakly-bound helium-cation complexes formed by electron ionization of helium nanodroplets doped with a neutral hydrocarbon precursor. Spectroscopic transitions were registered by photoexcitation of the complexes coupled with mass spectrometric detection of the bare ions. For CH₃⁺, we provide evidence showing that the helium-bound complexes contain 10-20 helium atoms (on average) and have a rotational temperature of ∼5 K. We show that this technique is well-suited to the study of highly symmetric or fluxional ionic species, as these intrinsic properties are preserved in the helium environment. This is in contrast to conventional tagging methods that use a single atom or molecule, which can change the point group or rigidity of the core ion and therefore the spectral profile. We demonstrate this for the highly fluxional molecular ion CH₅⁺, whose spectrum in the current study matches that of the gas phase ion, whereas the fluxionality is lost when a methane tag is added. Finally, we present the first IR spectrum of methane cation, CH₄⁺. The spectrum of this fundamental organic ion shows CH stretching bands consistent with a non-tetrahedral structure, a consequence of Jahn-Teller distortion.

Rotational action spectroscopy of trapped molecular ions

Rotational action spectroscopy is an experimental method in which rotational spectra of molecules, typically in the microwave to sub-mm-wave domain of the electromagnetic spectrum (∼1-1000 GHz), are recorded by action spectroscopy. Action spectroscopy means that the spectrum is recorded not by detecting the absorption of light by the molecules, but by the action of the light on the molecules, e.g., photon-induced dissociation of a chemical bond, a photon-triggered reaction, or photodetachment of an electron. Typically, such experiments are performed on molecular ions, which can be well controlled and mass-selected by guiding and storage techniques. Though coming with many advantages, the application of action schemes to rotational spectroscopy was hampered for a long time by the small energy content of a corresponding photon. Therefore, the first rotational action spectroscopic methods emerged only about one decade ago. Today, there exists a toolbox full of different rotational action spectroscopic schemes which are summarized in this review.

Infrared spectroscopy of carbocations upon electron ionization of ethylene in helium nanodroplets

The electron impact ionization of helium droplets doped with ethylene molecules and clusters yields diverse C_XH_Y ⁺ cations embedded in the droplets. The ionization primarily produces C₂H₂ ⁺, C₂H₃ ⁺, C₂H₄ ⁺, and CH₂ ⁺, whereas larger carbocations are produced upon the reactions of the primary ions with ethylene molecules. The vibrational excitation of the cations leads to the release of bare cations and cations with a few helium atoms attached. The laser excitation spectra of the embedded cations show well resolved vibrational bands with a few wavenumber widths-an order of magnitude less than those previously obtained in solid matrices or molecular beams by tagging techniques. Comparison with the previous studies of free and tagged CH₂ ⁺, CH₃ ⁺, C₂H₂ ⁺, C₂H₃ ⁺, and C₂H₄ ⁺ cations shows that the helium matrix typically introduces a shift in the vibrational frequencies of less than about 20 cm^-1, enabling direct comparisons with the results of quantum chemical calculations for structure determination. This work demonstrates a facile technique for the production and spectroscopic study of diverse carbocations, which act as important intermediates in gas and condensed phases.

Exactly solvable 1D model explains the low-energy vibrational level structure of protonated methane

A new one-dimensional model is proposed for the low-energy vibrational quantum dynamics of CH5+ based on the motion of an effective particle confined to a 60-vertex graph Γ60 with a single edge length parameter. Within this model, the quantum states of CH5+ are obtained in analytic form and are related to combinatorial properties of Γ60. The bipartite structure of Γ60 gives a simple explanation for curious symmetries observed in numerically exact variational calculations on CH5+.

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.

Probing the Fluxional Bonding Nature of Rapid Cope rearrangements in Bullvalene C10H10 and Its Analogs C8H8, C9H10, and C8BH9

Bullvalene C₁₀H₁₀ and its analogs semibullvalene C₈H₈, barbaralane C₉H₁₀, and 9-Borabarbaralane C₈BH₉ are prototypical fluxional molecules with rapid Cope rearrangements at finite temperatures. Detailed bonding analyses performed in this work reveal the existence of two fluxional π-bonds (2 2c-2e π → 2 3c-2e π → 2 2c-2e π) and one fluxional σ-bond (1 2c-2e σ → 1 4c-2e σ → 1 2c-2e σ) in their ground states and transition states, unveiling the universal π + σ double fluxional bonding nature of these fluctuating cage-like species. The highest occupied natural bond orbitals (HONBOs) turn out to be typical fluxional bonds dominating the dynamics of the systems. The ¹³C-NMR and ¹H-NMR shielding tensors and chemical shifts of the model compound C₈BH₉ are computationally predicted to facilitate future experiments.

Introductory lecture: advances in ion spectroscopy: from astrophysics to biology

This introduction provides a historical context for the development of ion spectroscopy over the past half century by following the evolution of experimental methods to the present state-of-the-art. Rather than attempt a comprehensive review, we focus on how early work on small ions, carried out with fluorescence, direct absorption, and photoelectron spectroscopy, evolved into powerful technologies that can now address complex chemical problems ranging from catalysis to biophysics. One of these developments is the incorporation of cooling and temperature control to enable the general application of "messenger tagging" vibrational spectroscopy, first carried out using ionized supersonic jets and then with buffer gas cooling in radiofrequency ion traps. Some key advances in the application of time-resolved pump-probe techniques to follow ultrafast dynamics are also discussed, as are significant benchmarks in the refinement of ion mobility to allow spectroscopic investigation of large biopolymers with well-defined shapes. We close with a few remarks on challenges and opportunities to explore molecular level mechanics that drive macroscopic behavior.

Tagging effects on the mid-infrared spectrum of microsolvated protonated methane

Although bare protonated methane is by now essentially understood at the level of intramolecular large-amplitude motion, scrambling dynamics and broadband vibrational spectra, the microsolvated species still offer plenty of challenges. One aspect is the effect of the attached solvent molecules on the infrared absorption spectra of microsolvated CH₅⁺ complexes compared to the bare parent molecule. In this study we analyze, based on ab initio molecular dynamics simulations, protonated methane molecules that have been microsolvated with up to three hydrogen molecules, i.e. CH₅⁺·(H₂) _n . In particular, upon introducing a novel multi-channel maximum entropy methodology described herein, we are able to decompose the infrared spectra of these weakly-bound complexes in the frequency window from 1000 to 4500 cm^-1 into additive single mode contributions. Detailed comparisons to the bare CH₅⁺ parent reveal that these perturbed modes encode distinct features that depend on the exact microsolvation pattern. Beyond the specific case, such understanding is relevant to assess tagging artifacts in vibrational spectra of parent molecules based on messenger predissociation action spectroscopy.

Vibrational quantum graphs and their application to the quantum dynamics of CH5+

The first application of quantum-graph theory to molecular vibrations helps understand the low-energy vibrational quantum dynamics of CH₅⁺.

Collective Molecular Superrotation: A Model for Extremely Flexible Molecules Applied to Protonated Methane

The concept of molecular structure is traditionally considered to be virtually fundamental. However, it breaks down in extremely flexible molecules, whose dynamics are governed by large-amplitude motions. For the dynamics of molecules with free internal rotations, we therefore propose a five-dimensional rigid rotor model, with states characterized by two generalized angular momentum quantum numbers and a rotational constant B. The quantum numbers characterize a 5D angular-momentum vector, super-j, which describes collective rotations that involve both the internal and the overall rotation. This model predicts the lowest energy states of the prototypical, extremely flexible molecule, CH_{5}^{+}. Both energies and symmetries compare very favorably to recent experimental results. The respective assignment to the new quantum numbers indicates the validity of our concept of collective rotations in extremely flexible molecules.

Relevance and Significance of Extraterrestrial Abiological Hydrocarbon Chemistry

Astrophysical observations show similarity of observed abiological "organics"-i.e., hydrocarbons, their derivatives, and ions (carbocations and carbanions)-with studied terrestrial chemistry. Their formation pathways, their related extraterrestrial hydrocarbon chemistry originating from carbon and other elements after the Big Bang, their parent hydrocarbon and derivative (methane and methanol, respectively), and transportation of derived building blocks of life by meteorites or comets to planet Earth are discussed in this Perspective. Their subsequent evolution on Earth under favorable "Goldilocks" conditions led to more complex molecules and biological systems, and eventually to humans. The relevance and significance of extraterrestrial hydrocarbon chemistry to the limits of science in relation to the physical aspects of evolution on our planet Earth are also discussed.

CH+5: Symmetry and the Entangled Rovibrational Quantum States of a Fluxional Molecule

Protonated methane, CH5+, is the prototypical example of a fluxional molecular system. The almost unconstrained angular motion of its five hydrogen atoms results in dynamical phenomena not found in rigid or semirigid molecules. Here it is shown that standard concepts to describe rotational quantum states of molecules can not be applied to CH5+ or any other fluxional system of the type ABn or Bn with n > 4 due to fundamental symmetry reasons. Instead, the ro-vibrational states of CH5+ display a unique level scheme, which results from a complex entanglement of rotational and tunneling motions. A detailed analysis of the ro-vibrational quantum states of CH5+ based on full-dimensional quantum dynamics simulations is presented, and the effects of the Pauli principle are considered. The consequences for the interpretation of recent experimental results are highlighted

Dissociative recombination of CH4(+)

CH4(+) is an important molecular ion in the astrochemistry of diffuse clouds, dense clouds, cometary comae, and planetary ionospheres. However, the rate of one of the common destruction mechanisms for molecular ions in these regions, dissociative recombination (DR), is somewhat uncertain. Here, we present absolute measurements for the DR of CH4(+) made using the heavy ion storage ring CRYRING in Stockholm, Sweden. From our collision-energy dependent cross-sections, we infer a thermal rate constant of k(Te) = 1.71(±0.02) × 10(–6)(Te/300)(−0.66(±0.02)) cm3 s(–1) over the region of electron temperatures 10 ≤ Te ≤ 1000 K. At low collision energies, we have measured the branching fractions of the DR products to be CH4 (0.00 ± 0.00); CH3 + H (0.18 ± 0.03); CH2 + 2H (0.51 ± 0.03); CH2 + H2 (0.06 ± 0.01); CH + H2 + H (0.23 ± 0.01); and CH + 2H2 (0.02 ± 0.01), indicating that two or more C–H bonds are broken in 80% of all collisions.