Soft-sphere continuum solvation models for nonaqueous solvents

Solvation effects profoundly influence the characteristics and behavior of chemical systems in liquid solutions. The interaction between solute and solvent molecules intricately impacts solubility, reactivity, stability, and various chemical processes. Continuum solvation models gained prominence in quantum chemistry by implicitly capturing these interactions and enabling efficient investigations of diverse chemical systems in solution. In comparison, continuum solvation models in condensed matter simulation are very recent. Among these, the self-consistent continuum solvation (SCCS) and the soft-sphere continuum solvation models (SSCS) have been among the first to be successfully parameterized and extended to model periodic systems in aqueous solutions and electrolytes. As most continuum approaches, these models depend on a number of parameters that are linked to experimental or theoretical properties of the solvent, or that can be tuned based on reference data. Here, we present a systematic parameterization of the SSCS model for over 100 nonaqueous solvents. We validate the model's efficacy across diverse solvent environments by leveraging experimental solvation-free energies and partition coefficients from comprehensive databases. The average root means square error over all the solvents was calculated as 0.85 kcal/mol which is below the chemical accuracy (1 kcal/mol). Similarly to what has been reported by Hille et al. (J. Chem. Phys. 2019, 150, 041710.) for the SCCS model, a single-parameter model accurately reproduces experimental solvation energies, showcasing the transferability and predictive power of these continuum approaches. Our findings underscore the potential for a unified approach to predict solvation properties, paving the way for enhanced computational studies across various chemical environments.

Water Arrangements upon Interaction with a Rigid Solute: Multiconfigurational Fenchone-(H2O)4-7 Hydrates

Insight into the arrangements of water molecules around solutes is important to understand how solvation proceeds and to build reliable models to describe water-solute interactions. We report the stepwise solvation of fenchone, a biogenic ketone, with 4-7 water molecules. Multiple hydrates were observed using broadband rotational spectroscopy, and the configurations of four fenchone-(H2O)4, three fenchone-(H2O)5, two fenchone-(H2O)6, and one fenchone-(H2O)7 complexes were characterized from the analysis of their rotational spectra in combination with quantum-chemical calculations. Interactions with fenchone deeply perturb water configurations compared with the pure water tetramer and pentamer. In two fenchone-(H2O)4 complexes, the water tetramer adopts completely new arrangements, and in fenchone-(H2O)5, the water pentamer is no longer close to being planar. The water hexamer interacts with fenchone as the least abundant book isomer, while the water heptamer adopts a distorted prism structure, which forms a water cube when including the fenchone oxygen in the hydrogen bonding network. Differences in hydrogen bonding networks compared with those of pure water clusters show the influence of fenchone's topology. Specifically, all observed hydrates except one show two water molecules binding to fenchone through each oxygen lone pair. The observation of several water arrangements for fenchone-(H2O)4-7 complexes highlights water adaptability and provides insight into the solvation process.

Electronic and vibrational spectroscopies of aromatic clusters with He in a supersonic jet: The case of neutral and cationic phenol-Hen (n = 1 and 2)

Van der Waals clusters composed of He and aromatic molecules provide fundamental information about intermolecular interactions in weakly bound systems. In this study, phenol-helium clusters (PhOH-Hen with n ≤ 2) are characterized for the first time by UV and IR spectroscopies. The S1 ← S0 origin and ionization energy both show small but additive shifts, suggesting π-bound structures of these clusters, a conclusion supported by rotational contour analyses of the S1 origin bands. The OH stretching vibrations of the PhOH moiety in the clusters match with those of bare PhOH in both the S0 and D0 states, illustrating the negligible perturbation of the He atoms on the molecular vibration. Matrix shifts induced by He attachment are discussed based on the observed band positions with the help of complementary quantum chemical calculations. For comparison, the UV and ionization spectra of PhOH-Ne are reported as well.

A five-carbon unsaturated Criegee intermediate: synthesis, spectroscopic identification, and theoretical study of 3-penten-2-one oxide

Biogenic alkenes, such as isoprene and α-pinene, are the predominant source of volatile organic compounds (VOCs) emitted into the atmosphere. Atmospheric processing of alkenes via reaction with ozone leads to formation of zwitterionic reactive intermediates with a carbonyl oxide functional group, known as Criegee intermediates (CIs). CIs are known to exhibit a strong absorption (π* ← π) in the near ultraviolet and visible (UV-vis) region due to the carbonyl oxide moiety. This study focuses on the laboratory identification of a five-carbon CI with an unsaturated substituent, 3-penten-2-one oxide, which can be produced upon atmospheric ozonolysis of substituted isoprenes. 3-Penten-2-one oxide is generated in the laboratory by photolysis of a newly synthesized precursor, (Z)-2,4-diiodopent-2-ene, in the presence of oxygen. The electronic spectrum of 3-penten-2-one oxide was recorded by UV-vis induced depletion of the VUV photoionization signal on the parent m/z 100 mass channel using a time-of-flight mass spectrometer. The resultant electronic spectrum is broad and unstructured with peak absorption at ca. 375 nm. To complement the experimental findings, electronic structure calculations are performed at the CASPT2(12,10)/aug-cc-pVDZ level of theory. The experimental spectrum shows good agreement with the calculated electronic spectrum and vertical excitation energy obtained for the lowest energy conformer of 3-penten-2-one oxide. In addition, OH radical products resulting from unimolecular decay of energized 3-penten-2-oxide CIs are detected by UV laser-induced fluorescence. Finally, the experimental electronic spectrum is compared with that of a four-carbon, isoprene-derived CI, methyl vinyl ketone oxide, to understand the effects of an additional methyl group on the associated electronic properties.

An investigation into the adsorption mechanism of n-butanol by ZIF-8: a combined experimental and ab initio molecular dynamics approach

The zeolitic imidazolate framework, ZIF-8, has been shown by experimental methods to have a maximum saturation adsorption capacity of 0.36 g g^-1 for n-butanol from aqueous solution, equivalent to a loading of 14 butanol molecules per unit cell or 7 molecules per sodalite β-cage. Diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) shows the presence of hydrogen bonding between adsorbed butanol molecules within the cage; the presence of three different O-H stretching modes indicates the formation of butanol clusters of varying size. Ab initio molecular dynamics simulations show the formation of intermolecular hydrogen bonding between the butanol molecules, with an average hydrogen-bond coordination number of 0.9 after 15 ps simulation time. The simulations also uniquely demonstrate the presence of weaker interactions between the alcohol O-H group and the π-orbital of the imidazole ring on the internal surface of the cage during early stages of adsorption. The calculated adsorption energy per butanol molecule is -33.7 kJ mol^-1, confirming that the butanol is only weakly bound, driven primarily by the hydrogen bonding. Solid-state MAS NMR spectra suggest that the adsorbed butanol molecules possess a reasonable degree of mobility in their adsorbed state, rather than being rigidly held in specific sites. 2D ¹³C-¹H heteronuclear correlation (HETCOR) experiments show interactions between the butanol aliphatic chain and the ZIF-8 framework experimentally, suggesting that O-H interactions with the π-orbital are only short lived. The insight gained from these results will allow the design of more efficient ways of recovering and isolating n-butanol, an important biofuel, from low-concentration solutions.

Structures of (Pyrazine)2 and (Pyrazine)(Benzene) Dimers Investigated with Infrared-Vacuum Ultraviolet Spectroscopy and Quantum-Chemical Calculations: Competition among π-π, CH···π, and CH···N Interactions

The structures of a pyrazine dimer (pyrazine)₂ and (pyrazine)(benzene) hetero-dimer cooled in a supersonic beam were investigated by the measurement of the infrared spectra in the C-H stretching region with infrared-vacuum ultraviolet (IR-VUV) spectroscopy and quantum-chemical calculations. The stabilization energy calculation at the CCSD(T)/aug-cc-pVTZ level of theory predicted three isomers for (pyrazine)₂ and three for (pyrazine)(benzene) with energy within 6 kJ/mol. Among them, the cross-displaced π-π stacked structure is the most stable in both dimers. In the observed IR spectra, both dimers exhibited two intense bands near 3065 cm^-1, with intervals of 8 cm^-1 in (pyrazine)₂ and 11 cm^-1 in (pyrazine)(benzene), while only one band appeared in the monomer. For (pyrazine)(benzene), we also measured the IR spectrum of (pyrazine)(benzene-d₆), where the interval of the two bands was unchanged. The analysis of the observed IR spectra with anharmonic calculations suggested the coexistence of three isomers of (pyrazine)₂ and (pyrazine)(benzene) in a supersonic jet. For (pyrazine)₂, the two isomers which were previously assigned to the H-bonded planar and the π-π stacked structures respectively were reassigned to the cross-displaced π-π stacked and T-shaped structures, respectively. In addition, the quantum chemical calculation and IR-VUV spectral measurement suggested the coexistence of the H-bonded planar isomer in the jet. For (pyrazine)(benzene), the IR spectrum of the (pyrazine) site showed a similar spectral pattern to that of (pyrazine)₂, especially the split at ∼3065 cm^-1. However, the anharmonic analysis suggested that they are assigned to the different vibrational motions of (pyrazine). The anharmonic vibrational analysis is essential to associate the observed IR spectra with the correct structures of the dimer.

Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State

The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.

Microhydration of the Pyrrole Cation (Py+) Revealed by IR Spectroscopy: Ionization-Induced Rearrangement of the Hydrogen-Bonded Network of Py+(H2O)2

Microhydration of heterocyclic aromatic molecules can be an appropriate fundamental model to shed light on intermolecular interactions and functions of macromolecules and biomolecules. We characterize herein the microhydration process of the pyrrole cation (Py⁺) by infrared photodissociation (IRPD) spectroscopy and dispersion-corrected density functional theory calculations (B3LYP-D3/aug-cc-pVTZ). Analysis of IRPD spectra of mass-selected Py⁺(H₂O)₂ and its cold Ar-tagged cluster in the NH and OH stretch range combined with geometric parameters of intermolecular structures, binding energies, and natural atomic charge distribution provides a clear picture of the growth of the hydration shell and cooperativity effects. Py⁺(H₂O)₂ is formed by stepwise hydration of the acidic NH group of Py⁺ by a hydrogen-bonded (H₂O)₂ chain with NH···OH···OH configuration. In this linear H-bonded hydration chain, strong cooperativity, mainly arising from the positive charge, strengthens both the NH···O and OH···O H-bonds with respect to those of Py⁺H₂O and (H₂O)₂, respectively. The linear chain structure of the Py⁺(H₂O)₂ cation is discussed in terms of the ionization-induced rearrangement of the hydration shell of the neutral Py(H₂O)₂ global minimum characterized by the so-called "σ-π bridge structure" featuring a cyclic NH···OH···OH···π H-bonded network. Emission of the π electron from Py by ionization generates a repulsive interaction between the positive π site of Py⁺ and the π-bonded OH hydrogen of (H₂O)₂, thereby breaking this OH···π hydrogen bond and driving the hydration structure toward the linear chain motif of the global minimum on the cation potential.

Light, Water, and Melatonin: The Synergistic Regulation of Phase Separation in Dementia

The swift rise in acceptance of molecular principles defining phase separation by a broad array of scientific disciplines is shadowed by increasing discoveries linking phase separation to pathological aggregations associated with numerous neurodegenerative disorders, including Alzheimer’s disease, that contribute to dementia. Phase separation is powered by multivalent macromolecular interactions. Importantly, the release of water molecules from protein hydration shells into bulk creates entropic gains that promote phase separation and the subsequent generation of insoluble cytotoxic aggregates that drive healthy brain cells into diseased states. Higher viscosity in interfacial waters and limited hydration in interiors of biomolecular condensates facilitate phase separation. Light, water, and melatonin constitute an ancient synergy that ensures adequate protein hydration to prevent aberrant phase separation. The 670 nm visible red wavelength found in sunlight and employed in photobiomodulation reduces interfacial and mitochondrial matrix viscosity to enhance ATP production via increasing ATP synthase motor efficiency. Melatonin is a potent antioxidant that lowers viscosity to increase ATP by scavenging excess reactive oxygen species and free radicals. Reduced viscosity by light and melatonin elevates the availability of free water molecules that allow melatonin to adopt favorable conformations that enhance intrinsic features, including binding interactions with adenosine that reinforces the adenosine moiety effect of ATP responsible for preventing water removal that causes hydrophobic collapse and aggregation in phase separation. Precise recalibration of interspecies melatonin dosages that account for differences in metabolic rates and bioavailability will ensure the efficacious reinstatement of the once-powerful ancient synergy between light, water, and melatonin in a modern world.

Fragment imaging in the infrared photodissociation of the Ar-tagged protonated water clusters H3O+-Ar and H+(H2O)2-Ar

Infrared photodissociation of protonated water clusters with an Ar atom, namely H₃O⁺-Ar and H⁺(H₂O)₂-Ar, was investigated by an imaging technique for mass-selected ions, to reveal the intra- and intermolecular vibrational dynamics. The presented system has the advantage of achieving fragment ion images with the cluster size- and mode-selective photoexcitation of each OH stretching vibration. Translational energy distributions of photofragments were obtained from the images upon the excitation of the bound (ν_b) and free (ν_f) OH stretching vibrations. The energy fractions in the translational motion were compared between ν_b^I and ν_f^I in H₃O⁺-Ar or between ν_b^II and ν_f^II in H⁺(H₂O)₂-Ar, where the labels "I" and "II" represent H₃O⁺-Ar and H⁺(H₂O)₂-Ar, respectively. In H₃O⁺-Ar, the ν_f^I excitation exhibited a smaller translational energy than ν_b^I. This result can be explained by the higher vibrational energy of ν_f^I, which enabled it to produce bending (ν₄) excited H₃O⁺ fragments that should be favored according to the energy-gap model. In contrast to H₃O⁺-Ar, the ν_b^II excitation of an Ar-tagged H₂O subunit and the ν_f^II excitation of an untagged H₂O subunit resulted in very similar translational energy distributions in H⁺(H₂O)₂-Ar. The similar energy fractions independent of the excited H₂O subunits suggested that the ν_b^II and ν_f^II excited states relaxed into a common intermediate state, in which the vibrational energy was delocalized within the H₂O-H⁺-H₂O moiety. However, the translational energy distributions for H⁺(H₂O)₂-Ar did not agree with a statistical dissociation model, which implied another aspect of the process, that is, Ar dissociation via incomplete energy randomization in the whole H⁺(H₂O)₂-Ar cluster.

Molecular Structure, Electronic Properties, Reactivity (ELF, LOL, and Fukui), and NCI-RDG Studies of the Binary Mixture of Water and Essential Oil of Phlomis bruguieri

Essential oils are volatile oil-like liquids with a characteristic strong smell and taste. They are formed in plants and are then extracted. Essential oils have extremely strong physiological and pharmacological properties, which are used in the medicine, cosmetics, and food industries. In this study, the molecules caryophyllene oxide, β-pinene, 1,8-cineol, α-cubebene, and β-caryophyllene, which are the molecules with the highest contents in the essential oil of the plant mentioned in the title, were selected and theoretical calculations describing their interactions with water were performed. Because oil–water mixtures are very important in biology and industry and are ubiquitous in nature, quantum chemical calculations for binary mixtures of water with caryophyllene oxide, β-pinene, 1,8-cineol, α-cubebene, and β-caryophyllene were performed using the density functional theory (DFT)/B3LYP method with a basis of 6–31 G (d, p). Molecular structures, HOMO–LUMO energies, electronic properties, reactivity (ELF, LOL, and Fukui), and NCI-RDG and molecular electrostatic potential (MEP) on surfaces of the main components of Phlomis bruguieri Desf. essential oil were calculated and described.

Laser Spectroscopy of Helium Solvated Clusters of Methanol and Methanol-Water in the Symmetric Methyl Stretching Band

Mid-infrared spectra of methanol and methanol-water clusters have been investigated in the symmetric CD₃ stretching band of CD₃OH and CD₃OD. We find that the position of this band provides a useful signature of the general type of hydrogen-bonded cluster it is associated with. Our results are consistent with those previously reported in the OH stretching region (Sulaiman, M. I.; Yang, S.; Ellis, A. M. J. Phys. Chem. A 2017, 121, 771-776) in that methanol clusters from the trimer to the pentamer are cyclic and that mixed clusters with one water molecule (and at least two methanol molecules) are also cyclic. We additionally provide evidence that the methanol trimer adopts a chair-like structure (as opposed to bowl-like), that mixed clusters with a larger number of water molecules are also cyclic, and that branched methanol clusters contribute to the depletion signal in larger methanol clusters. We performed double-hybrid DFT calculations which support these interpretations.

Non-covalent interactions in the monohydrated complexes of 1,2,3,4-tetrahydroisoquinoline

The eleven monohydrates of 1, 2, 3, 4-tetrahydroisoquinoline (THIQ) are analyzed through natural bond orbital (NBO) analysis and QTAIM methods employing M06-2X functional in DFT and MP2 methods. Here, the role of OH bonds as an acceptor and donor is critically analyzed. The role of lone pairs of O is critically monitored in two of the complexes, where N-H···O hydrogen bonds are present. The relative contributions of rehybridisation and hyperconjugation are compared in detail. Popelier criteria are satisfied in all the complexes barring a few exceptions involving weak hydrogen bonds. At the bond critical points (BCP), four monohydrates show higher values of electron density (ρ_C) and negative values of total electron energy density (H_C), while Laplacian [Formula: see text] remains positive. These complexes satisfy the criteria of partial covalency. All these are O-H⋅⋅⋅N-type bonds. Remaining h-bonds are weaker in nature. These are also confirmed by the smaller values of ρ_C at the respective BCP. The variation of potential energy density (V_C) among the complexes seems to be the most important factor in determining the nature of non-covalent interactions.

The solvation of SOH group in hydrated HSO4−(H2O)n clusters

The S=O stretching and SOH bending peaks in the vibrational spectra of HSO4−(H2O) n, with n up to 6, are analyzed by both harmonic analysis and ab initio molecular dynamics simulation. The SOH bending mode is found to be much more sensitive to the extent of hydration and to the fluctuation of hydrogen bonds than the S=O stretching mode. The SOH donor hydrogen bond is gradually stabilized by n = 4, and further shortened up to n = 6, which is the key factor to understand the trend of evolution observed in the infrared multiple photon dissociation spectra.

Microhydration of Phenyl Formate: Gas-Phase Laser Spectroscopy, Microwave Spectroscopy, and Quantum Chemistry Calculations

Herein, we have investigated the structure of phenyl formate⋅⋅⋅water (PhOF⋅⋅⋅H₂ O) dimer and various non-covalent interactions present there using gas-phase laser spectroscopy and microwave spectroscopy combined with quantum chemistry calculations. Two conformers of PhOF⋅⋅⋅H₂ O (C1 and T1), built on the two cis/trans conformers of the bare molecule, have been observed in the experiment. In cis-PhOF, there is an n_CO → π A r * ${{{\rm \pi }}_{{\rm A}{\rm r}}^{{\rm {^\ast}}}}$ interaction between the lone-pair orbital of the carbonyl oxygen atom and the π* orbital of the phenyl ring, which persists in the monohydrated C1 conformer of PhOF⋅⋅⋅H₂ O according to the NBO and NCI analyses. On the other hand, this interaction is absent in the trans-PhOF conformer as the C=O group is away from the phenyl ring. The C1 conformer is primarily stabilized by an interplay between O-H⋅⋅⋅O=C hydrogen bond and O-H⋅⋅⋅π interactions, while the stability of the T1 conformer is primarily governed by the O-H⋅⋅⋅O=C hydrogen bond. The most important finding of the present work is that the conformational preference of the PhOF monomer is retained in its monohydrated complex.

Ultrafast light-induced dynamics in the microsolvated biomolecular indole chromophore with water

Interactions between proteins and their solvent environment can be studied in a bottom-up approach using hydrogen-bonded chromophore-solvent clusters. The ultrafast dynamics following UV-light-induced electronic excitation of the chromophores, potential radiation damage, and their dependence on solvation are important open questions. The microsolvation effect is challenging to study due to the inherent mix of the produced gas-phase aggregates. We use the electrostatic deflector to spatially separate different molecular species in combination with pump-probe velocity-map-imaging experiments. We demonstrate that this powerful experimental approach reveals intimate details of the UV-induced dynamics in the near-UV-absorbing prototypical biomolecular indole-water system. We determine the time-dependent appearance of the different reaction products and disentangle the occurring ultrafast processes. This approach ensures that the reactants are well-known and that detailed characteristics of the specific reaction products are accessible - paving the way for the complete chemical-reactivity experiment.

Competition between In-Plane vs Above-Plane Configurations of Water with Aromatic Molecules: Non-Covalent Interactions in 1,4-Naphthoquinone-(H2O)1-3 Complexes

Non-covalent interactions between aromatic molecules and water are fundamental in many chemical and biological processes, and their accurate description is essential to understand molecular relative configurations. Here we present the rotational spectroscopy study of the water complexes of the polycyclic aromatic hydrocarbon 1,4-naphthoquinone (1,4-NQ). In 1,4-NQ-(H₂O)_1,2, water molecules bind through O-H···O and C-H···O hydrogen bonds and are located on the plane of 1,4-NQ. For 1,4-NQ-(H₂O)₃, in-plane and above-plane water configurations are observed exhibiting O-H···O, C-H···O, and lone pair···π-hole interactions. The observation of different water arrangements for 1,4-NQ-(H₂O)₃ allows benchmarking theoretical methods and shows that they have great difficulty in predicting energy orderings due to the strong competition of C-H···O binding with π and π-hole interactions. This study provides important insight into water interactions with aromatic systems and the challenges in their modeling.

Infrared Spectroscopy of Stepwise Hydration Motifs of Sulfur Dioxide

Experimental characterization of microscopic events and behaviors of SO₂-H₂O interactions is crucial to understanding SO₂ atmospheric chemistry but has been proven to be very challenging due to the difficulty in size selection. Here, size-dependent development of SO₂ hydrate structure and cluster growth in the SO₂(H₂O)_n (n = 1-16) complexes was probed by infrared spectroscopy based on threshold photoionization using a tunable vacuum ultraviolet free electron laser. Spectral changes with cluster size demonstrate that the sandwich structure initially formed at n = 1 develops into cycle structures with the sulfur and oxygen atoms in a two-dimensional plane (n = 2 and 3) and then into three-dimensional cage structures (n ≥ 4). SO₂ is favorably bound to the surface of larger water clusters. These stepwise features of SO₂ hydration on various sized water clusters contribute to understanding the reactive sites and electrophilicity of SO₂ on cloud droplets, which may have important atmospheric implications for studying the SO₂-containing aerosol systems.

Conformer-Specific Spectroscopy and IR-Induced Isomerization of a Model γ-Peptide: Ac-γ4-Phe-NHMe

Single-conformation IR and UV spectroscopy of the prototypical capped γ-peptide Ac-γ⁴-Phe-NHMe (γ⁴F) was carried out under jet-cooled conditions in the gas phase in order to understand its innate conformational preferences in the absence of a solvent. We obtained conformer-specific IR and UV spectra and compared the results with calculations to make assignments and explore the differences between the γ²- and γ⁴-substituted molecules. We found four conformers of γ⁴F in our experiment. Three conformers form nine-membered hydrogen-bonded rings (C9) enclosed by an NH···O═C H-bond but differing in their phenyl ring positions (a, g+, and g-). The fourth conformer forms a strained seven-membered hydrogen-bonded ring in which the amide groups lie in a nominally anti-parallel arrangement stacked on top of one another (labeled S7). This conformer is a close analogue of the amide-stacked conformer (S) found previously in γ²F, in which the Phe side chain is substituted at the γ² position, Ac-γ2-Phe-NHMe (J. Am. Chem. Soc. 2009, 131, 14243-14245). IR population transfer spectroscopy was used to determine the fractional abundances of the γ⁴F conformers in the expansion. A combination of force field and density functional theory calculations is used to map out the conformational potential energy surfaces for γ⁴F and compare it with its γ²F counterpart. Based on this analysis, the phenyl ring prefers to take up structures that facilitate NH···π interactions in γ⁴F or avoid phenyl interactions with the C═O group in γ²F. The disconnectivity graph for γ⁴F reveals separate basins associated with the C9 and amide-stacked conformational families, which are separated by a barrier of about 42 kJ/mol. The overall shape of the potential energy surface bears a resemblance to peptides and proteins that have a misfolding pathway that competes with the formation of the native structure.

Optimal control of photodissociation of phenol using genetic algorithm

Photodissociation dynamics of the OH bond of phenol is studied with an optimally shaped laser pulse. The theoretical model consists of three electronic states (the ground electronic state, ππ* state, and πσ* state) in two nuclear coordinates (the OH stretching coordinate as a reaction coordinate, r, and the CCOH dihedral angle as a coupling coordinate, θ). The optimal UV laser pulse is designed using the genetic algorithm, which optimizes the total dissociative flux of the wave packet. The latter is calculated in the adiabatic asymptotes of the S0 and S1 electronic states of phenol. The initial state corresponds to the vibrational levels of the electronic ground state and is defined as | n r, n θ⟩, where n r and n θ represent the number of nodes along r and θ, respectively. The optimal UV field excites the system to the optically dark πσ* state predominantly over the optically bright ππ* state with the intensity borrowing effect for the |0, 0⟩ and |0, 1⟩ initial states. For the |0, 0⟩ initial condition, the photodissociation to the S1 asymptotic channel is favored slightly over the S0 asymptotic channel. Addition of one quantum of energy along the coupling coordinate increases the dissociation probability in the S1 channel. This is because the wave packet spreads along the coupling coordinate on the πσ* state and follows the adiabatic path. Hence, the S1 asymptotic channel gets more ([Formula: see text]11%) dissociative flux as compared to the S0 asymptotic channel for the |0, 1⟩ initial condition. The |1, 0⟩ and |1, 1⟩ states are initially excited to both the ππ* and πσ* states in the presence of the optimal UV pulse. For these initial conditions, the S1 channel gets more dissociative flux as compared to the S0 channel. This is because the high energy components of the wave packet readily reach the S1 channel. The central frequency of the optimal UV pulse for the |0, 0⟩ and |0, 1⟩ initial states has a higher value as compared to the |1, 0⟩ and |1, 1⟩ initial states. This is explained with the help of an excitation mechanism of a given initial state in relation to its energy.

Setting up the HyDRA blind challenge for the microhydration of organic molecules

The procedure leading to the first HyDRA blind challenge for the prediction of water donor stretching vibrations in monohydrates of organic molecules is described. A training set of 10 monohydrates...

Isomer-Selective Spectroscopy and Dynamics of Phenol-Arn (n ≤ 5) Clusters

Structures and ionization-induced solvation dynamics of phenol-(argon)_n clusters, PhOH-Ar_n (n ≤ 5), were studied by using a variety of isomer-selective photoionization and vibrational spectroscopic techniques. Several higher-energy isomers were found and assigned for the first time by systematically controlling the experimental conditions of the supersonic expansion. This behavior is also confirmed for the PhOH-Kr₂ cluster. Solvation structures are elucidated by evaluating systematic shifts in the S₁ ← S₀ origin and ionization energies obtained by resonance-enhanced photoionization, in addition to the OH stretching frequency obtained by IR photodissociation. Isomer-selective picosecond time-resolved IR spectroscopy for the n = 2 clusters revealed that the dynamics for the ionization-induced intermolecular π → H site-switching reaction strongly depends on the initial isomeric structure. In particular, the reaction time for the (1|1) isomer is 7 ps, while that for (2|0) is <3 ps. This difference shows that the switching time is determined by the distance of the reaction coordinate between the initial π-site and the final OH-site.

Stacked but not Stuck: Unveiling the Role of π→π* Interactions with the Help of the Benzofuran-Formaldehyde Complex

The 1:1 benzofuran–formaldehyde complex has been chosen as model system for analyzing π→π* interactions in supramolecular organizations involving heteroaromatic rings and carbonyl groups. A joint “rotational spectroscopy–quantum chemistry” strategy unveiled the dominant role of π→π* interactions in tuning the intermolecular interactions of such adduct. The exploration of the intermolecular potential energy surface led to the identification of 14 low‐energy minima, with 4 stacked isomers being more stable than those linked by hydrogen bond or lone‐pair→π interactions. All energy minima are separated by loose transition states, thus suggesting an effective relaxation to the global minimum under the experimental conditions. This expectation has been confirmed by the experimental detection of only one species, which was unambiguously assigned owing to the computation of accurate spectroscopic parameters and the characterization of 11 isotopologues. The large number of isotopic species opened the way to the determination of the first semi‐experimental equilibrium structure for a molecular complex of such a dimension.

Active learning of potential-energy surfaces of weakly bound complexes with regression-tree ensembles

Several pool-based active learning (AL) algorithms were employed to model potential-energy surfaces (PESs) with a minimum number of electronic structure calculations. Theoretical and empirical results suggest that superior strategies can be obtained by sampling molecular structures corresponding to large uncertainties in their predictions while at the same time not deviating much from the true distribution of the data. To model PESs in an AL framework, we propose to use a regression version of stochastic query by forest, a hybrid method that samples points corresponding to large uncertainties while avoiding collecting too many points from sparse regions of space. The algorithm is implemented with decision trees that come with relatively small computational costs. We empirically show that this algorithm requires around half the data to converge to the same accuracy in comparison to the uncertainty-based query-by-committee algorithm. Moreover, the algorithm is fully automatic and does not require any prior knowledge of the PES. Simulations on a 6D PES of pyrrole(H₂O) show that <15 000 configurations are enough to build a PES with a generalization error of 16 cm^-1, whereas the final model with around 50 000 configurations has a generalization error of 11 cm^-1.

Dispersion-Bound Isolated Dimers in the Gas Phase: Observation of the Shortest Intermolecular CH⋅⋅⋅H-C Distance via Stimulated Raman Spectroscopy

The triphenylmethane and all‐meta tert‐butyl triphenylmethane dimers, (TPM)₂ and (T ^tBuPM)₂, respectively, were studied with ionization loss stimulated Raman spectroscopy in molecular beam experiments to resolve structure sensitive vibrations. This answers the question whether the recently reported linear head‐to‐head arrangement in (T ^tBuPM)₂ results from crystal packing or prevails also in the gas phase, and therefore must result from extraordinarily strong London dispersion (LD) interactions. Our study clearly demonstrates that the head‐to‐head arrangement is maintained even under isolated molecular beam conditions in the absence of crystal packing effects. The central Raman‐active aliphatic C−D vibration of appropriately deuterated (T ^tBuPM)₂ associated with an unusually short C−D⋅⋅⋅D−C distance exhibits a strong blue‐shift compared to the undisturbed case. As the LD stabilizing tert‐butyl groups are absent in (TPM)₂, it displays an approximately S₆‐symmetric tail‐to‐tail arrangement.

Revisiting the Spectroscopy of Water Dimer in Jets

Laser spectroscopy in jets is one of the main sources of structural data from molecular aggregates. Consequently, numerous and sophisticated experimental systems have been developed to extract precise information, which is usually interpreted in the light of quantum mechanical calculations. However, even with the most sophisticated experiments, it is sometimes difficult to interpret the experimental results. We present here the example of water dimer and how after almost 70 years, the assignment of its mass-resolved IR spectrum still generates controversy that extends toward the mechanism of ionization of water aggregates.

Pickup and reactions of molecules on clusters relevant for atmospheric and interstellar processes

In this perspective, we review experiments with molecules picked up on large clusters in molecular beams with the focus on the processes in atmospheric and interstellar chemistry. First, we concentrate on the pickup itself, and we discuss the pickup cross sections. We measure the uptake of different atmospheric molecules on mixed nitric acid-water clusters and determine the accommodation coefficients relevant for aerosol formation in the Earth's atmosphere. Then the coagulation of the adsorbed molecules on the clusters is investigated. In the second part of this perspective, we review examples of different processes triggered by UV-photons or electrons in the clusters with embedded molecules. We start with the photodissociation of hydrogen halides and Freon CF₂Cl₂ on ice nanoparticles in connection with the polar stratospheric ozone depletion. Next, we mention reactions following the excitation and ionization of the molecules adsorbed on clusters. The first ionization-triggered reaction observed between two different molecules picked up on the cluster was the proton transfer between methanol and formic acid deposited on large argon clusters. Finally, negative ion reactions after slow electron attachment are illustrated by two examples: mixed nitric acid-water clusters, and hydrogen peroxide deposited on large Ar_N and (H₂O)_N clusters. The selected examples are discussed from the perspective of the atmospheric and interstellar chemistry, and several future directions are proposed.

Unravelling the microhydration frameworks of prototype PAH by infrared spectroscopy: naphthalene–(water)1–3

Microhydration structures of the prototypical PAH, naphthalene, are probed by IR spectroscopy in helium droplets. The sequential water addition produces an extended hydrogen-bonded hydration network bound via π hydrogen bond to the aromatic ring.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Rotational Spectroscopy Meets Quantum Chemistry for Analyzing Substituent Effects on Non-Covalent Interactions: The Case of the Trifluoroacetophenone-Water Complex

The most stable isomer of the 1:1 complex formed by 2,2,2-trifluoroacetophenone and water has been characterized by combining rotational spectroscopy in supersonic expansion and state-of-the-art quantum-chemical computations. In the observed isomer, water plays the double role of proton donor and acceptor, thus forming a seven-membered ring with 2,2,2-trifluoroacetophenone. Accurate intermolecular parameters featuring one classical O-H···O hydrogen bond and one weak C-H···O hydrogen bond have been determined by means of a semi-experimental approach for equilibrium structure. Furthermore, insights on the nature of the established non-covalent interactions have been unveiled by means of different bond analyses. The comparison with the analogous complex formed by acetophenone with water points out the remarkable role played by fluorine atoms in tuning non-covalent interactions.

IR-VUV spectroscopy of pyridine dimers, trimers and pyridine-ammonia complexes in a supersonic jet

The infrared spectra of the C-H stretching vibrations of (pyridine)m, m = 1-3, and the N-H stretching vibrations of (pyridine)m-(NH3)n, m = 1, 2; n = 1-4, complexes were investigated by infrared (IR)-vacuum ultraviolet (VUV) spectroscopy under jet-cooled conditions. The ionization potential (IP0) of the pyridine monomer was determined to be 74 546 cm-1 (9.242 eV), while its complexes showed only smooth curves of the ionization thresholds at ∼9 eV, indicating large structural changes in the ionic form. The pyridine monomer exhibits five main features with several satellite bands in the C-H stretching region at 3000-3200 cm-1. Anharmonic calculations including Fermi-resonance were carried out to analyze the candidates of the overtone and combination bands which can couple to the C-H stretching fundamentals. For (pyridine)2 and (pyridine)3, most C-H bands are blue-shifted by 3-5 cm-1 from those of the monomer. The structures revealed by random searching algorithms with density functional methods indicate that the π-stacked structure is most stable for (pyridine)2, while (pyridine)3 prefers the structures stabilized by dipole-dipole and C-Hπ interactions. For the (pyridine)m-(NH3)n complexes, the mass spectrum exhibited a wide range distribution of the complexes. The observed IR spectra in the N-H stretching vibrations of the complexes showed four main bands in the 3200-3450 cm-1 region. These features are very similar to those of (NH3)n complexes, and the bands are assigned to the anti-symmetric N-H stretching band (ν3), the symmetric N-H stretching (ν1) band, and the first overtone bands of the N-H bending vibrations (2ν4). The anharmonic calculations including the Fermi-resonance between ν1 and 2ν4 well reproduced the observed spectra.

Coexistence of Left- and Right-Handed 12/10-Mixed Helices in Cyclically Constrained β-Peptides and Directed Formation of Single-Handed Helices upon Site-Specific Methylation

The inherent conformational preferences of the neutral β-peptide foldamer series, Ac-(ACHC)_n-NHBn, n = 2-4, are studied in the gas phase using conformation-specific IR-UV double resonance methods. The cyclically constrained chiral β-amino acid cis-2-aminocyclohexane carboxylic acid (ACHC) is designed to bring both right- and left-handed helices into close energetic proximity. Comparison of the infrared spectra in the NH stretch and amide I/II regions with the predictions of DFT calculations lead to the unambiguous assignment of four out of the six observed conformations of the molecules in this series, while corroborating computational and spectral evidence, affords tentative assignments of the remaining two conformers for which IR data were not recorded. The observed structures fall into one of two conformational families: a right-handed 12/10-mixed helix or its "cap-disrupted" left-handed helical analogue, which coexist with significant populations. Site-specific and stereospecific methylation on the cyclohexane backbone at the dipeptide (n = 2) level is also tested as a means to sterically lock in a predetermined cyclohexane chair conformation. These substitutions are proven to be a means of selectively driving formation of one helical screw sense or the other. Calculated relative energies and free energies of all possible structures for the molecules provide strong supporting evidence that the rigid nature of the ACHC residue confers unusual stability to the 12/10-mixed helix conformation, regardless of local environment, temperature, or C-terminal capping unit. The simultaneous presence of both handed helices offers unique opportunities for future studies of their interconversion.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

Migrations and Catalytic Action of Water Molecules in the Ionized Formamide-(H2O)2 Cluster

Isomerization dynamics involving the migrations, proton transfer reaction, and catalytic actions of water molecules upon vertical ionization of the formamide (FA)-(H₂O)₂ cluster is investigated by the infrared spectroscopy and theoretical reaction path search calculation. The infrared spectroscopic result indicates the [FA-(H₂O)₂]⁺ cation has the hydrogen-bonded structure of the enol isomer cation of formamide and the water dimer. This structure is formed by proton transfer from the CH bond to the carbonyl group through the catalytic action of the water molecules. The isomerization paths involving this enolization in ionized FA-(H₂O)₂ are explored by using the anharmonic downward distortion following method. We found multiple enolization paths which accompany proton exchanges among the formamide moiety and water molecules through the catalytic actions of the water molecules.

A close competition between O–H⋯O and O–H⋯π hydrogen bonding: IR spectroscopy of anisole–methanol complex in helium nanodroplets

Anisole forms O–H⋯O as well O–H⋯π bound complexes with methanol.

Structures and internal dynamics of diphenylether and its aggregates with water

We report on a detailed multi-spectroscopic analysis of the structures and internal dynamics of diphenylether and its aggregates with up to three water molecules by employing molecular beam experiments.

Long-distance proton transfer induced by a single ammonia molecule: ion mobility mass spectrometry of protonated benzocaine reacted with NH3

A long-distance proton transfer via the vehicle mechanism in the absence of a hydrogen-bonded solvent-bridge in molecules.

Conformer Specific Ultraviolet and Infrared Detection of Nicotine in the Vapor Phase

A gas-phase electronic spectrum of nicotine in a supersonic expansion has been recorded using two-color resonant two-photon ionization spectroscopy. Efficient photoionization was achievable only via the pyridine chromophore owing to poor Franck-Condon overlap in the N-methylpyrrolidine moiety. Two conformers of nicotine have been characterized and assigned by infrared-ultraviolet (IR-UV) ion depletion and IR-UV hole-burning spectroscopy, in combination with quantum chemical techniques. Trans-A with nitrogen atoms further apart is more stable by 2 kJ mol^-1 and the most populated conformer in the supersonic jet, owing this stability to a stronger inter-ring CH···N hydrogen bond than the trans-B counterpart.