Observation of metastable structures of the ethylene glycol-water dimer in helium nanodroplets

Ethylene glycol (EG) is the simplest organic diol. Here we measure infrared spectra of the EG monomer and its dimer with water, the complex, EG(H2O), embedded in superfluid helium nanodroplets. For the monomer, only a single, gauche, conformation is observed. For EG(H2O), no trace of the global energy minimum is seen, a structure that would maximize the hydrogen bonding contacts. Instead, only metastable structures are formed, suggesting that dimerization in a superfluid environment leads to kinetic trapping in local energy minima. In addition, we obtain evidence for a dimer where the conformation of EG switches from gauche to trans on account of dimerization with a water molecule. This observation is assumed to be driven over an energy barrier by utilizing the energy released as hydrogen bonding occurs.

Structural Characterization of Pyruvic Acid Dimers Formed inside Helium Nanodroplets by Infrared Spectroscopy and Ab Initio Study

The structural arrangements of α-keto acid complexes hold significant interest across various fields of chemistry such as enzyme modeling, drug design, or polymer blending. Herein, we report mass-selective infrared (IR) spectra of pyruvic acid monomers and dimers in the range 1720-1820 cm-1 recorded in helium nanodroplets at 0.37 K. The monomer features IR bands at 1807.1 and 1734.5 cm-1, which are assigned to the carboxylic and ketonic C═O stretching vibrations, respectively. Furthermore, the pyruvic acid dimers generated inside the helium nanodroplets are characterized by carboxylic and ketonic C═O stretch vibrations appearing at 1799.2 and 1737.0 cm-1, respectively. This frequency shift of ±7 cm-1 for both C═O stretching bands from the monomer to the dimer demonstrates that the structural motif of the monomer is maintained upon dimer aggregation in helium nanodroplets. The structural assignments were supported by a comparison of the MP2/aug-cc-pVDZ-predicted harmonic vibrational spectra at the C═O stretching region with the experiments. The global minimum monomer structure with an intramolecular hydrogen bond and its dimer stabilized by both inter- and intramolecular hydrogen bonding interactions reproduce the experimental spectra from the monomer and dimer. This assigned dimer structure lies ca.11 kJ/mol above the corresponding global minimum and is favored in helium nanodroplets due to the long-range realignment of molecules via dipole-dipole interaction, followed by short-range stabilization upon intermolecular hydrogen bond formation. The barrier for reconfiguration of the precooled monomer conformer leading to the formation of the most stable dimer structure is around 58 kJ/mol, which is infeasible at 0.37 K.

Cooperativity between Intermolecular Hydrogen and Carbon Bonds in ZY···CH3CN/CH3NC···HX Trimers (ZY = H2O, H2S, HF, HCl, HBr, NH3, and H2CO; HX = HF, HCl, and HBr)

Hydrogen-bonding and carbon-bonding interactions are widespread in nature. We studied the cooperativity between these interactions in 42 trimeric complexes ZY···CH3CN/CH3NC···HX, where ZY molecules are H2O, H2S, HF, HCl, HBr, NH3, and H2CO, and HX molecules are HF, HCl, and HBr. Acetonitrile (CH3CN) and isoacetonitrile (CH3NC) act as hydrogen bond acceptors as well as carbon bond donors in these trimers. Various theoretical methods, such as electronic structure calculations, quantum theory of atoms in molecule (QTAIM), natural bond orbital (NBO), and reduced density gradient analysis, are employed to study these trimers, and the results are compared with the corresponding ZY···CH3CN/CH3NC and CH3CN/CH3NC···HX dimers. Electronic structure calculations are performed at the second-order Mo̷ller-Plesset perturbation theory using the 6-311++G(2d,2p) basis set. We show that both the interactions act synergistically in these trimers leading to an increase in their bond strength as compared to the strength in the individual dimers. The cooperative energies for these trimers are in the range of 0.69 to 3.22 kJ/mol. It is seen that the carbon bonds benefit more from the cooperativity than the hydrogen bonds. The trends of cooperativity and correlations of interaction energies and cooperative energies with relevant QTAIM and NBO parameters are reported.

Secondary ionization of pyrimidine nucleobases and their microhydrated derivatives in helium nanodroplets

Radiation damage in biological systems by ionizing radiation is predominantly caused by secondary processes such as charge and energy transfer leading to the breaking of bonds in DNA. Here, we study the fragmentation of cytosine (Cyt) and thymine (Thy) molecules, clusters and microhydrated derivatives induced by direct and indirect ionization initiated by extreme-ultraviolet (XUV) irradiation. Photofragmentation mass spectra and photoelectron spectra of free Cyt and Thy molecules are compared with mass and electron spectra of Cyt/Thy clusters and microhydrated Cyt/Thy molecules formed by aggregation in superfluid helium (He) nanodroplets. Penning ionization after resonant excitation of the He droplets is generally found to cause less fragmentation compared to direct photoionization and charge-transfer ionization after photoionization of the He droplets. When Cyt/Thy molecules and oligomers are complexed with water molecules, their fragmentation is efficiently suppressed. However, a similar suppression of fragmentation is observed when homogeneous Cyt/Thy clusters are formed in He nanodroplets, indicating a general trend. Penning ionization electron spectra (PIES) of Cyt/Thy are broad and nearly featureless but PIES of their microhydrated derivatives point at a sequential ionization process ending in unfragmented microsolvated Cyt/Thy cations.

Nanostructured supramolecular networks from self-assembled diamondoid molecules under ultracold conditions

Diamondoid molecules and their derivatives have attracted attention as fascinating building blocks for advanced functional materials. Depending on the balance between hydrogen bonds and London dispersion interactions, they can self-organize in different cluster structures with functional groups tailored for various applications. Here, we present a new approach to supramolecular aggregation where self-assembly of diamondoid acids and alcohols in the ultracold environment of superfluid helium nanodroplets (HNDs) was analyzed by a combination of time-of-flight mass spectrometry and computational tools. Experimentally observed magic numbers of the assembled cluster sizes were successfully identified and computed cluster structures gave valuable insights into a different conglomeration mode when compared to previously explored less-polar diamondoid derivatives. We have confirmed that functional groups acting as good hydrogen bond donors completely take over the self-organization process, resulting in fascinating pair-wise or cyclic supramolecular assemblies. Particularly noteworthy is that mono- and bis-substituted diamondoid derivatives of both series engage in completely different modes of action, which is reflected in differing non-covalent cluster geometries. Additionally, formed cyclic clusters with a polar cavity in the center and a non-polar diamondoid outer layer can be of high interest in porous material design and provide insights into the structural requirements needed to produce bulk materials with desired properties.

Dimerization dynamics of carboxylic acids in helium nanodroplets

The dimerization of molecules in helium nanodroplets is known to preferentially yield structures of higher energy than the global energy minimum structure for a number of quite different monomers. Here, we explore dimerization in this environment using an atomistic model within statistically converged molecular dynamics (MD) trajectories, treating the solvent implicitly through the use of a thermostat, or more explicitly by embedding one monomer in a He₁₀₀ cluster. The focus is on the two simplest carboxylic acids, formic and acetic, both of which have been studied experimentally. While the global minimum structure, which comprises two CO⋯HO hydrogen bonds, is predicted to be the most abundant dimer in the absence of the helium solvent, this is no longer the case once helium atoms are included. The simulations confirm the importance of kinetic trapping effects and also shed light on the occurrence of specific dynamical effects, leading to the occasional formation of high-energy structures away from minima, such as saddle configurations. Theoretically predicted infrared spectra, based on the MD statistics, are in good agreement with the experimental spectra.

An infrared spectroscopic study of trifluoromethoxybenzene⋯methanol complexes formed in superfluid helium nanodroplets

We have studied the intermolecular complex formation between trifluoromethoxybenzene and methanol (CD₃OD) in superfluid helium droplets by infrared spectroscopy in the spectral range of 2630-2730 cm^-1, covering the O-D stretches of methanol-d₄ (CD₃OD). The cluster size associated with the observed bands is deduced from the variation of infrared intensity of a particular band with the partial pressures of trifluoromethoxybenzene and methanol. Quantum chemical calculations are performed at the MP2/6-311++G(d,p) level of theory to complement the experimental results. As a result, we have identified six different conformers of the trifluoromethoxybenzene⋯methanol intermolecular complex: three bound via O-H⋯O hydrogen bonds and the other three via O-H⋯π hydrogen bonds. Furthermore, to access the effect of fluorination on the methyl unit of anisole molecules, we compare the IR spectrum of trifluoromethoxybenzene (C₆H₅OCF₃)⋯methanol with our earlier reported spectrum of anisole (C₆H₅OCH₃)⋯methanol.

Rovibronic signatures of molecular aggregation in the gas phase: subtle homochirality trends in the dimer, trimer and tetramer of benzyl alcohol

Molecular aggregation is of paramount importance in many chemical processes, including those in living beings. Thus, characterization of the intermolecular interactions is an important step in its understanding. We describe here the aggregation of benzyl alcohol at the molecular level, a process governed by a delicate equilibrium between OH⋯O and OH⋯π hydrogen bonds and dispersive interactions. Using microwave, FTIR, Raman and mass-resolved double-resonance IR/UV spectroscopic techniques, we explored the cluster growth up to the tetramer and found a complex landscape, partly due to the appearance of multiple stereoisomers of very similar stability. Interestingly, a consistently homochiral synchronization of transiently chiral monomer conformers was observed during cluster growth to converge in the tetramer, where the fully homochiral species dominates the potential energy surface. The data on the aggregation of benzyl alcohol also constitute an excellent playground to fine-tune the parameters of the most advanced functionals.

Vibrational Spectroscopy of Benzonitrile-(Water)1-2 Clusters in Helium Droplets

Polycyclic aromatic hydrocarbons are considered as primary carriers of the unidentified interstellar bands. The recent discovery of the first interstellar aromatic molecule, benzonitrile (C₆H₅CN), suggests a repository of aromatic hydrocarbons in the outer earth environment. Herein, we report an infrared (IR) study of benzonitrile-(D₂O)_n clusters using mass-selective detection in helium nanodroplets. In this work, we use isotopically substituted water, D₂O, instead of H₂O because of our restricted IR frequency range (2565-3100 cm^-1). A comparison of the experimental and predicted spectra computed at the MP2/6-311++G(d,p) level of benzonitrile-(water)_1-2 clusters reveals the formation of a unique local minimum structure, which was not detected in previous gas-phase molecular beam experiments. Here, the solvent water forms a nearly linear hydrogen bond (H-bond) with the nitrile nitrogen of benzonitrile, while the previously reported most stable cyclic H-bonded isomer is not observed. This can be rationalized by the stepwise aggregation process of precooled monomers. The addition of a second water molecule results in the formation of two different isomers. In one of the observed isomers, a H-bonded water chain binds linearly to the nitrile nitrogen similar to the monohydrated benzonitrile-water complex. In the other observed isomer, the water dimer forms a ring-type structure, where a H-bonded water dimer simultaneously interacts with the nitrile nitrogen and the adjacent ortho CH group. Finally, we compare the water-binding motif in the neutral benzonitrile-water complex with the corresponding positively and negatively charged benzonitrile-water monohydrates to comprehend the charge-induced alteration of the solvent binding motif.

Reduced Two-Electron Interactions in Anharmonic Molecular Vibrational Calculations Involving Localized Normal Coordinates

Spatially localized vibrational normal mode coordinates are shown to reduce the importance of calculating the full set of two-electron terms in the molecular electronic Schrödinger equation. Electron correlation and dispersion interactions become less significant in (E,E)-1,3,5,7-octatetraene vibrational self-consistent field calculations when displacing remote atoms along multiple coordinates. Electron correlation interactions between spatially remote modes are also found to be less important compared to their corresponding uncorrelated interaction terms. Attenuation of the Coulomb operator indicates that the two-electron terms between remote electrons become less important for accurately describing the strongly contributing mode-coupling terms between sets of localized vibrational modes.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

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.

IR Spectroscopy of the Cesium Iodide-Water Complex

There has been much interest in I^-(H₂O) as a simple model for a hydrated iodide ion. Here we explore how this fundamental ion-solvent interaction is modified by the presence of a counterion, specifically Cs⁺. This has been achieved by forming the CsI(H₂O) complex in superfluid helium nanodroplets and then probing this system using infrared spectroscopy. The complex retains the ionic hydrogen bond between the I^- and a water OH group seen in I^-(H₂O), but the Cs⁺ ion substantially alters the anion-water interaction through formation of a cyclic Cs⁺-O-H-I^- bonding motif. As with I^-(H₂O), the OH stretching band derived from the hydrogen-bonded OH group shows substructure, splitting into a clear doublet. However, in contrast to I^-(H₂O), where a tunneling splitting arising from hydrogen atom exchange plays a role, the doublet we observe is attributed solely to an anharmonic vibrational coupling effect.

Dimerization of Acetic Acid in the Gas Phase-NMR Experiments and Quantum-Chemical Calculations

Due to the nature of the carboxylic group, acetic acid can serve as both a donor and acceptor of a hydrogen bond. Gaseous acetic acid is known to form cyclic dimers with two strong hydrogen bonds. However, trimeric and various oligomeric structures have also been hypothesized to exist in both the gas and liquid phases of acetic acid. In this work, a combination of gas-phase NMR experiments and advanced computational approaches were employed in order to validate the basic dimerization model of gaseous acetic acid. The gas-phase experiments performed in a glass tube revealed interactions of acetic acid with the glass surface. On the other hand, variable-temperature and variable-pressure NMR parameters obtained for acetic acid in a polymer insert provided thermodynamic parameters that were in excellent agreement with the MP2 (the second order Møller–Plesset perturbation theory) and CCSD(T) (coupled cluster with single, double and perturbative triple excitation) calculations based on the basic dimerization model. A slight disparity between the theoretical dimerization model and the experimental data was revealed only at low temperatures. This observation might indicate the presence of other, entropically disfavored, supramolecular structures at low temperatures.

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.

Shifting formic acid dimers into perspective: vibrational scrutiny in helium nanodroplets

A metastable dimer of formic acid has been prepared inside superfluid helium nanodroplets and examined using IR spectroscopy and quantum chemical calculations.

Concerted Pair Motion Due to Double Hydrogen Bonding: The Formic Acid Dimer Case

Formic acid dimer as the prototypical doubly hydrogen-bonded gas-phase species is discussed from the perspective of the three translational and the three rotational degrees of freedom which are lost when two formic acid molecules form a stable complex. The experimental characterisation of these strongly hindered translations and rotations is reviewed, as are attempts to describe the associated fundamental vibrations, their combinations, and their thermal shifts by different electronic structure calculations and vibrational models. A remarkable match is confirmed for the combination of a CCSD(T)-level harmonic treatment and an MP2-level anharmonic VPT2 correction. Qualitatively correct thermal shifts of the vibrational spectra can be obtained from classical molecular dynamics in CCSD(T)-quality force fields. A detailed analysis suggests that this agreement between experiment and composite theoretical treatment is not strongly affected by fortuitous error cancellation but fully converged variational treatments of the six pair or intermolecular modes and their overtones and combinations in this model system would be welcome.

Direct detection of polar structure formation in helium nanodroplets by beam deflection measurements

Long-range intermolecular forces are able to steer polar molecules submerged in superfluid helium nanodroplets into highly polar metastable configurations. We demonstrate that the presence of such special structures can be identified, in a direct and determinative way, by electrostatic deflection of the doped nanodroplet beam. The measurement also establishes the structures' electric dipole moments. In consequence, the introduced approach is complementary to spectroscopic studies of low-temperature molecular assembly reactions. It is enabled by the fact that within the cold superfluid matrix the molecular dipoles become nearly completely oriented by the applied electric field. As a result, the massive (tens of thousands of helium atoms) nanodroplets undergo significant deflections. The method is illustrated here by an application to dimers and trimers of dimethyl sulfoxide (DMSO) molecules. We interpret the experimental results with ab initio theory, mapping the potential energy surface of DMSO complexes and simulating their low temperature aggregation dynamics.

Molecular aggregation in liquid acetic acid: insight from molecular dynamics/quantum mechanics modelling of structural and NMR properties

The 1H NMR signal of the acidic proton in acetic acid molecules shows a marked upfield shift in the neat liquid as compared to that in low-concentration acetic acid solution in inert solvents where acetic acid cyclic dimers predominate. The underlying reasons for this phenomenon are analyzed in this work by considering classical molecular dynamics simulations and combined quantum mechanics/molecular mechanics calculations of the 1H NMR chemical shift of the acidic proton in the neat liquid and in the cyclic dimer of acetic acid in cyclohexane solution. Recorded trajectories were quantitatively analyzed in terms of different types of molecular aggregates formed in the neat liquid by using a geometrical definition of the hydrogen bond. Both the geometrical analysis and the computational NMR results indicate that the cyclic dimer cannot be the dominating aggregation pattern for acetic acid molecules in the neat liquid. The applied computational approach reproduces the lowering of the 1H NMR chemical shift of the acidic proton in acetic acid when going from cyclohexane solution to the neat liquid very well. The presence of acetic acid aggregates with hydrogen bonding between hydroxyl moieties and of monomeric acetic acid molecules in the neat liquid is found to lead to the observed lowering of the chemical shift, with lesser contribution from the formation of open acetic acid aggregates.

Probing Elusive Cations: Infrared Spectroscopy of Protonated Acetic Acid

Protonated carboxylic acids, (RCOOH)H⁺, are the initial intermediates in acid-catalyzed (Fischer) esterification reactions. However, the identity of the isomeric form has been debated. Surprisingly, no optical spectra have been reported for any isomer of the protonated carboxylic acid monomer, despite it being a fundamental organic cation. Here, we address these issues by using a new approach to prepare cold He-tagged cations of protonated acetic acid (AA), which entails electron ionization of helium nanodroplets containing metastable dimers of AA. The protonated species is subsequently probed using infrared photodissociation spectroscopy, and following a comparison with calculations, we identify the two isomers whose roles in Fischer esterification are debated. These are the carbonyl-protonated E, Z isomer and the metastable hydroxyl-protonated isomer. Our technique provides a novel approach that can be applied to other elusive ionic species.