Complete Assignment of the Infrared Spectrum of the Gas-Phase Protonated Ammonia Dimer

Photochemical Evolution of Alanine in Association with the Martian Soil Analog Montmorillonite: Insights Derived from Experiments Conducted on the International Space Station

The Photochemistry on the Space Station (PSS) experiment was part of the European Space Agency's EXPOSE-R2 mission and was conducted on the International Space Station from 2014 to 2016. The PSS experiment investigated the properties of montmorillonite clay as a protective shield against degradation of organic compounds that were exposed to elevated levels of ultraviolet (UV) radiation in space. Additionally, we examined the potential for montmorillonite to catalyze UV-induced breakdown of the amino acid alanine and its potential to trap the resulting photochemical byproducts within its interlayers. We tested pure alanine thin films, alanine thin films protected from direct UV exposure by a thin cover layer of montmorillonite, and an intimate combination of the two substances forming an organoclay. The samples were exposed to space conditions for 15.5 months and then returned to Earth for detailed analysis. Concurrent ground-control experiments subjected identical samples to simulated solar light irradiation. Fourier-transform infrared (FTIR) spectroscopy quantified molecular changes by comparing spectra obtained before and after exposure for both the space and ground-control samples. To more deeply understand the photochemical processes influencing the stability of irradiated alanine molecules, we performed an additional experiment using time-resolved FTIR spectroscopy for a second set of ground samples exposed to simulated solar light. Our collective experiments reveal that montmorillonite clay exhibits a dual, configuration-dependent effect on the stability of alanine: while a thin cover layer of the clay provides UV shielding that slows degradation, an intimate mixture of clay and amino acid hastens the photochemical decomposition of alanine by promoting certain chemical reactions. This observation is important to understand the preservation of amino acids in specific extraterrestrial environments, such as Mars: cover mineral layer depths of several millimeters are required to effectively shield organics from the harmful effects of UV radiation. We also explored the role of carbon dioxide (CO2), a byproduct of alanine photolysis, as a tracer of the amino acid. CO2 can be trapped within clay interlayers, particularly in clays with small interlayer ions such as sodium. Our studies emphasize the multifaceted interactions between montmorillonite clay and alanine under nonterrestrial conditions; thus, they contribute valuable insights to broader astrobiological research questions.

Isomer-Specific, Cryogenic Ion Vibrational Spectroscopy Investigation of D2- and N2-Tagged, Protonated Formic Acid Complexes Using Two-Color, IR-IR Photobleaching

Here we analyze cryogenic ion vibrational spectra of tagged protonated formic acid (PFA) with electronic structure and anharmonic vibrational calculations to establish the isomers generated by electrospray ionization (ESI) followed by buffer gas cooling to ∼25 K. Two isomers are identified (the trans form (E,Z) and the cis form (E,E)) and generated in comparable abundance despite the fact that the calculated E,E structure lies 6.40 kJ mol-1 above the E,Z form. A large (∼60 kJ mol-1) barrier separates them such that the E,E form can be kinetically trapped upon cooling in the ion trap. The anticooperativity between the H-bonds of the OH groups is explored by measuring the shift in the D2-bound OH fundamental when a second D2 is attached. Both isomers are observed in the N2-tagged counterparts, displaying the expected red-shifted OH bands. These results indicate that ESI generates both isomers and both must be considered when analyzing cluster spectra based on the PFA core ion.

Anharmonic IR spectra of solvated ammonium and aminium ions: resemblance between water and bisulfate solvations

In this work, we analyze the vibrational spectra of ammonium, methylammonium, and dimethylammonium ions solvated by either water molecules or bisulfate anions using anharmonic vibrational algorithms. Rich and complicated spectral features in the 2700-3200 cm^-1 region of the experimental spectra of these clusters are attributed to originate from strong Fermi resonance between hydrogen-bonded NH stretching fundamentals and NH bending overtones. Additional weaker bands around 2500-2600 cm^-1 in solvated aminium ions are assigned to the combination tones involving the CH-NH (methyl-amino) rocking modes. Furthermore, the qualitative resemblance in band positions and spectral patterns between two-water-solvated and two-bisulfate-solvated cations suggest a common vibrational coupling scheme beneath the two seemingly different micro-solvation environments.

Molecular vibrational spectral simulation connects theoretical cluster structure identification and vibrational spectral evidence

Structure identification of molecular clusters has long been a fundamental and challenging issue for cluster science. The traditional theoretical optimization on the potential energy surface heavily depends on the levels of theory and sometimes diverse identifications were reported. A solution to these disputations is to reinspect the theoretical results with the experimental data such as vibrational predissociation spectra with high sensitivity to the molecular cluster structures. Herein, the combination of global low-lying structure search and vibrational predissociation spectral simulation is proposed as an accurate and reliable approach for cluster structure identification, by which the assignments can be validated using experimental measurements. The qualitative agreement between simulated and measured vibrational spectra lends solid experimental evidence to the assignment of the cluster structures. Taking NH₄⁺(H₂O)_n (n = 2-4) as an example, we have unambiguously identified their structures and directly demonstrated the coexistence of two NH₄⁺(H₂O)₄ isomers (with 3 and 4 water molecules directly linked to NH₄⁺, respectively), which were debatable in previous studies. The developed methods would pave the way to the structure determination of the molecular clusters.

Molecular Dynamics with Constrained Nuclear Electronic Orbital Density Functional Theory: Accurate Vibrational Spectra from Efficient Incorporation of Nuclear Quantum Effects

Nuclear quantum effects play a crucial role in many chemical and biological systems involving hydrogen atoms yet are difficult to include in practical molecular simulations. In this paper, we combine our recently developed methods of constrained nuclear-electronic orbital density functional theory (cNEO-DFT) and constrained minimized energy surface molecular dynamics (CMES-MD) to create a new method for accurately and efficiently describing nuclear quantum effects in molecular simulations. By use of this new method, dubbed cNEO-MD, the vibrational spectra of a set of small molecules are calculated and compared with those from conventional ab initio molecular dynamics (AIMD) as well as from experiments. With the same formal scaling, cNEO-MD greatly outperforms AIMD in describing the vibrational modes with significant hydrogen motion characters, demonstrating the promise of cNEO-MD for simulating chemical and biological systems with significant nuclear quantum effects.

Molecular vibrational frequencies from analytic Hessian of constrained nuclear-electronic orbital density functional theory

Nuclear quantum effects are important in a variety of chemical and biological processes. The constrained nuclear-electronic orbital density functional theory (cNEO-DFT) has been developed to include nuclear quantum effects in energy surfaces. Herein, we develop the analytic Hessian for cNEO-DFT energy with respect to the change in nuclear (expectation) positions, which can be used to characterize stationary points on energy surfaces and compute molecular vibrational frequencies. This is achieved by constructing and solving the multicomponent cNEO coupled-perturbed Kohn-Sham (cNEO-CPKS) equations, which describe the response of electronic and nuclear orbitals to the displacement of nuclear (expectation) positions. With the analytic Hessian, the vibrational frequencies of a series of small molecules are calculated and compared to those from conventional DFT Hessian calculations as well as those from the vibrational second-order perturbation theory (VPT2). It is found that even with a harmonic treatment, cNEO-DFT significantly outperforms DFT and is comparable to DFT-VPT2 in the description of vibrational frequencies in regular polyatomic molecules. Furthermore, cNEO-DFT can reasonably describe the proton transfer modes in systems with a shared proton, whereas DFT-VPT2 often faces great challenges. Our results suggest the importance of nuclear quantum effects in molecular vibrations, and cNEO-DFT is an accurate and inexpensive method to describe molecular vibrations.

Understanding the Hydrogen-Bonded Clusters of Ammonia (NH3) n (n = 3-6): Insights from the Electronic Structure Theory

Although it is well known that hydrogen bonds commonly exist in ammonia clusters and play an important role, there are still many challenges in understanding the electronic structure properties of hydrogen bonds. In this paper, the geometric and electronic structure properties of cyclic ammonia clusters are investigated by using first-principles density functional theory (DFT) and the Møller-Plesset perturbation theory (MP2). The calculation results show that the pentamer and hexamer have deviated from the perfect plane, while the trimer and tetramer present planarization that has been confirmed by infrared (IR) spectra. The electronic structure analysis further shows that the covalent properties play a non-negligible role in hydrogen bonding. The results also indicate that the electronic structure facilitates structure planarization. Our work not only provides insight into the role and nature of hydrogen bonds in ammonia clusters but also provides a theoretical basis for frontier science in fields such as atmospheric haze and biomolecular functions.

Quantum chemical accuracy from density functional approximations via machine learning

Kohn-Sham density functional theory (DFT) is a standard tool in most branches of chemistry, but accuracies for many molecules are limited to 2-3 kcal ⋅ mol-1 with presently-available functionals. Ab initio methods, such as coupled-cluster, routinely produce much higher accuracy, but computational costs limit their application to small molecules. In this paper, we leverage machine learning to calculate coupled-cluster energies from DFT densities, reaching quantum chemical accuracy (errors below 1 kcal ⋅ mol-1) on test data. Moreover, density-based Δ-learning (learning only the correction to a standard DFT calculation, termed Δ-DFT ) significantly reduces the amount of training data required, particularly when molecular symmetries are included. The robustness of Δ-DFT  is highlighted by correcting "on the fly" DFT-based molecular dynamics (MD) simulations of resorcinol (C6H4(OH)2) to obtain MD trajectories with coupled-cluster accuracy. We conclude, therefore, that Δ-DFT  facilitates running gas-phase MD simulations with quantum chemical accuracy, even for strained geometries and conformer changes where standard DFT fails.

Perturbation of Pyridinium CIVP Spectra by N2 and H2 Tags: An Experimental and BOMD Study

In cryogenic ion vibrational predissociation (CIVP) spectroscopy, the influence of the tag on the spectrum is an important consideration. Whereas for small ions several studies have shown that the tag effects can be significant, these effects are less understood for large ions or for large numbers of tags. Nevertheless, it is commonly assumed that if the investigated molecular ion is large enough, the perturbations arising from the tag are small and can therefore be neglected in the interpretation. In addition, it is generally assumed that the more weakly bound the tag is, the less it perturbs the CIVP spectrum. Under these assumptions, CIVP spectra are claimed to be effectively IR absorption spectra of the free molecular ion. Having observed unexpected splittings in otherwise unproblematic CIVP spectra of some tagged ions, we report Born-Oppenheimer molecular dynamics (BOMD) simulations that strongly indicate that mobility among the more weakly bound tags leads to the surprising splittings. We compared the behavior of two tags commonly used in CIVP spectroscopy (H₂ and N₂) with a large pyridinium cation. Our experimental results surprisingly show that under the appropriate circumstances, the more weakly bound tag can perturb the CIVP spectra more than the more strongly bound tag by not just shifting but also splitting the observed bands. The more weakly bound tag had significant residence times at several spectroscopically distinct sites on the molecular ion. This indicates that the weakly bound tag is likely to sample several binding sites in the experiment, some of which involve interaction with the reporter chromophore.

Anharmonic coupling behind vibrational spectra of solvated ammonium: lighting up overtone states by Fermi resonance through tuning solvation environments

Fascinating Fermi resonance bands emerge from anharmonic couplings between NH stretching fundamentals and bending overtones in ammonium-centered clusters.

Temperature Dependence of Intramolecular Vibrational Bands in Small Water Clusters

Cyclic water clusters are pivotal for understanding atmospheric reactions as well as liquid water, yet the temperature (T) dependence of their dynamics and spectroscopy is poorly studied. The development of highly accurate water potentials, such as MB-pol, partly rectifies this. It remains to account for the quantum nuclear effects (NQE), because quantum nuclear dynamics become increasingly inaccurate at low temperatures. From a practical point of view, we find that NQE can be accounted for simply by subtracting a constant from the frequencies obtained from the velocity autocorrelation functions (VACF) of classical nuclear dynamics, resulting in unprecedented agreement with experiment, mostly within 5 cm^-1. We have performed classical simulations of (H₂O)_n clusters (n = 2-5) from 20 K and up to their melting temperature, calculating both all-atom and partial VACF, thus generating the temperature dependence of the vibrational frequencies (IR and Raman bands). Focusing on the hydrogen-bonded (HBed) OH stretch and HOH bend, we find opposing T dependencies. The HBed OH modes blue shift linearly with T, attributed to ring expansion rather than any specific conformational change. The lowest-frequency Raman concerted mode is predicted to show the largest such shift. In contrast, the HOH bend undergoes a red-shift, with the highest frequency concerted band undergoing the largest red-shift. These results can be explained by a coupled-oscillator model for n hydrogen atoms on a ring, constrained to move either tangentially (stretch) or perpendicularly (bend) to the ring. With increasing temperature and weakening of HBs, the intrinsic force constant increases (stretch) or remains constant (bend), while the nearest-neighbor coupling constant decreases, and this results in the interesting behavior revealed herein. T-dependent Raman studies are required for testing some of these predictions.

Spectroscopy of Proton Coordination with Ethylenediamine

Protonated ethylenediamine monomer, dimer, and trimer were produced in the gas phase by an electrical discharge/supersonic expansion of argon seeded with ethylenediamine (C₂H₈N₂, en) vapor. Infrared spectra of these ions were measured in the region from 1000 to 4000 cm^-1 using laser photodissociation and argon tagging. Computations at the CBS-QB3 level were performed to explore possible isomers and understand the infrared spectra. The protonated monomer exhibits a gauche conformation and an intramolecular hydrogen bond. Its parallel shared proton vibration occurs as a broad band around 2785 cm^-1, despite the formally equivalent proton affinities of the two amino groups involved, which usually leads to low frequency bands. The barrier to intramolecular proton transfer is 2.2 kcal mol^-1 and does not vanish upon addition of the zero-point energy, unlike the related protonated ammonia dimer. The structure of the dimer is formed by chelation of the monomer's NH₃⁺ group, thereby localizing the excess proton and increasing the frequency of the intramolecular shared proton vibration to 3157 cm^-1. Other highly fluxional dimer structures with facile intermolecular proton transfer and concomitant structural reorganization were computed to lie within 2 kcal mol^-1 of the experimentally observed structure. The spectrum of the trimer is rather diffuse, and a clear assignment is not possible. However, an isomer with an intramolecular proton transfer like that of the monomer is most consistent with the experimental spectrum.

Intra-cavity proton bonding and anharmonicity in the anionophore cyclen

Proton bonding drives the supramolecular chemistry of a broad range of materials with polar moieties. Proton delocalization and electronic charge redistribution have a profound impact on the structure of proton-bound molecular frameworks, and pose fundamental challenges to quantum chemical modelling. This study provides insights into the structural and spectral signatures of the intramolecular proton bond formed in a benchmark polyazamacrocycle anionophore (cyclen, 1,4,7,10-tetraazacyclododecane). Infrared action spectroscopy is employed to characterize the macrocycle, isolated in protonated form. In its most stable configuration, protonated cyclen adopts an open arrangement of Cs symmetry with a particularly strong NHδ+N bond across the cavity. The quantum chemical analysis of the infrared spectrum reveals intrinsic difficulties for the accurate description of the vibrational modes of the system. The reconciliation of the computational predictions with experiment demands a careful anharmonic treatment of the proton motion, which exposes the limitations of current methods. Best results are obtained with the incorporation of anharmonicity only to the fundamental modes directly related to motions of the proton. However, the full anharmonic treatment of the system fails to describe correctly the vibrations related to the macrocycle backbone. The results should serve as motivation for new developments in the modelling of proton bonded systems.

Spin-state dependence of the structural and vibrational properties of solvated iron(ii) polypyridyl complexes from AIMD simulations: aqueous [Fe(bpy)3]Cl2, a case study

LS and HS IR spectra of aqueous [Fe(bpy)₃]²⁺ and corresponding HS–LS difference IR spectrum as obtained from state-of-the-art ab initio molecular dynamics simulations applied to the determination of the structural and vibrational properties of the solvated complex.

Hydrogen bond induced enhancement of Fermi resonances in N–H⋯N hydrogen bonded complexes of anilines

Enhancement of Fermi resonance intensities due to the formation of N–H⋯N hydrogen bonding of anilines with alkyl amines is analyzed using a two-state deperturbation model.

Vibrational Spectroscopy and Proton Transfer Dynamics in Protonated Oxalate

The dynamics and infrared spectroscopic signatures of proton transfer in protonated oxalate (p-Oxa) are studied using classical and quantum dynamics. The intermolecular interactions are described by a force field suitable to follow proton transfer. This allows to carry out multiple extended classical molecular dynamics (MD) and ring polymer MD simulations from which the infrared spectrum is determined. Simulations at 600 K sample the quantum mechanical ground state probability distribution and best reproduce the experimentally observed maximum absorption wavelength and part of the line shape. Comparison with the experimentally measured spectrum provides an estimate for the barrier height for proton transfer which can not be determined directly from experiment. A barrier of 4.2 kcal/mol is found to best reproduce the position and width of the infrared absorption of the transferring proton in p-Oxa and also leads to an infrared (IR) spectrum in good agreement with experiment for the deuterated species d-Oxa. A novel means to capture the two resonance forms of oxalate depending on the localization of the excess proton on either CO moiety is found to yield improved results for the spectroscopy in the framework region between 1000 and 2000 cm^-1.

Protons and Hydroxide Ions in Aqueous Systems

Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.