Acid solvation versus dissociation at "stardust conditions": Reaction sequence matters

Mechanism of ionic dissociation of HCl in the smallest water clusters

The dissociation of strong acids into water is a fundamental process in chemistry and biology. Determining the minimum number of water molecules that can result in an ionic dissociation of hydrochloric acid (HCl → H+ + Cl-) remains a challenging subject. In this study, the reactions of H2O with HCl(H2O)n-1 (HCl-H2O cluster), i.e., HCl(H2O)n-1 + H2O (n = 3-7), were investigated by using the direct ab initio molecular dynamics (AIMD) method. Direct AIMD calculations were performed to set the collision energy of H2O to zero for all trajectories. For n = 3, no reaction occurred. In contrast, HCl dissociated to H+ + Cl- at n = 4, forming a contact ion pair (cIP) and solvent-separated ion pair (ssIP) as products. The reactions were expressed as HCl(H2O)3 + H2O → H3O+(H2O)2Cl- (ssIP), and HCl(H2O)3 + H2O → H3O+(Cl-)(H2O)2 (cIP). The ion pair (IP) products were dependent on the collision site of H2O relative to HCl(H2O)3. For n = 5-7, both IPs were formed through the reaction between H2O and HCl(H2O)n-1 (n = 5-7). The reaction between HCl and (H2O)4 (HCl + (H2O)4 → HCl(H2O)4) was non-reactive in IP formation. The reaction mechanism was discussed based on the theoretical results.

Cryovacuum setup for optical studies of astrophysical ice

This paper presents a cryovacuum setup for the study of substances under near-space conditions. The setup makes it possible to study the infrared spectra, refractive index, and density of substances that are condensed from the vapor phase onto a cooled substrate in the temperature range from 11 to 300 K. At the same time, it is possible to obtain the ultimate pressure of 1 × 10-10 Torr in the vacuum chamber. The presented setup is based on FTIR spectroscopy (the spectral measurement range is 400-7800 cm-1) and laser interference, through which the important physical and optical parameters are determined. A number of experiments allow us to point out that the data obtained using this setup correlate well with the experiments of other authors. Due to the non-directional deposition of substances from the vapor phase, the ice formed resembles the one formed under cosmic conditions as closely as possible, which makes the presented setup particularly valuable. The presented cryovacuum setup can be used for the interpretation of data obtained during astrophysical observations, providing a means to determine the properties of cosmic objects.

Evolution of Solute-Water Interactions in the Benzaldehyde-(H2O)1-6 Clusters by Rotational Spectroscopy

The investigation on the preferred arrangement and intermolecular interactions of gas phase solute-water clusters gives insights into the intermolecular potentials that govern the structure and dynamics of the aqueous solutions. Here, we report the investigation of hydrated coordination networks of benzaldehyde-(water)n (n = 1-6) clusters in a pulsed supersonic expansion using broadband rotational spectroscopy. Benzaldehyde (PhCHO) is the simplest aromatic aldehyde that involves both hydrophilic (CHO) and hydrophobic (phenyl ring) functional groups, which can mimic molecules of biological significance. For the n = 1-3 clusters, the water molecules are connected around the hydrophilic CHO moiety of benzaldehyde through a strong CO···HO hydrogen bond and weak CH···OH hydrogen bond(s). For the larger clusters, the spectra are consistent with the structures in which the water clusters are coordinated on the surface of PhCHO with both the hydrophilic CHO and hydrophobic phenyl ring groups being involved in the bonding interactions. The presence of benzaldehyde does not strongly interfere with the cyclic water tetramer and pentamer, which retain the same structure as in the pure water cluster. The book isomer instead of cage or prism isomers of the water hexamer is incorporated into the microsolvated cluster. The PhCHO molecule deviates from the planar structure upon sequential addition of water molecules. The PhCHO-(H2O)1-6 clusters may serve as a simple model system in understanding the solute-water interactions of biologically relevant molecules in an aqueous environment.

Acid Enriched with Methanol: Formation of a Prebiotic Cluster in the Interstellar Medium

Organic molecules are a potential source of prebiotic chemistry in the interstellar medium (ISM). Methanol (MetOH) is a very important source of more complex molecules. H 3 O + (aq) and Cl - (aq) are fundamental to living organisms and can be generated in the ISM from the dissociation of HCl with just four water molecules, yielding the (H 3 O) + (H 2 O) 3 Cl - ion-pair. Here, a detailed mechanism, based on density functional theory (DFT) and ab-initio (2 nd order Mfller-Plesset perturbation theory, MP2) calculations, is suggested for the substitution reactions of these water molecules by MetOH. The time required for formation of an appreciable amount of the product ((H 3 O) + (MetOH) 3 Cl - ) can be only few years. Such reaction can take place in Sagittarius B2, where HCl, H 2 O and MetOH have already been identified and it can be an important source for the formation of more complex prebiotic structures.

Internal Electric Field-Induced Formation of Exotic Linear Acetonitrile Chains

The application of external electric and magnetic fields is a powerful tool for aligning molecules in a controlled way, if the thermal fluctuations are small. Here we demonstrate that the same holds for internal electric fields in a molecular cluster. The electric field of a single molecular dipole, HCl, is used to manipulate the aggregation mechanism of subsequently added acetonitrile molecules. As a result, we could form exotic linear acetonitrile (CH₃CN) chains at 0.37 K, as confirmed by infrared spectroscopy in superfluid helium nanodroplets. These linear chains are not observed in the absence of HCl and can be observed only when the internal electric field created by an HCl molecule is present. The accompanying simulations provide mechanistic insights into steric control, explain the selectivity of the process, and show that non-additive electronic polarization effects systematically enhance the dipole moment of these linear chains. Thus, adding more CH₃CN monomers even supports further quasi-linear chain growth.

Zwitter Ionization of Glycine at Outer Space Conditions due to Microhydration by Six Water Molecules

We investigate glycine microsolvation with water molecules, mimicking astrophysical conditions, in our laboratory by embedding these clusters in helium nanodroplets at 0.37 K. We recorded mass selective infrared spectra in the frequency range 1500-1800  cm^{-1} where two bands centered at 1630 and 1724  cm^{-1} were observed. By comparison with the extensive accompanying calculations, the band at 1630  cm^{-1} was assigned to the COO^{-} asymmetric stretching mode of the zwitter ion and the band at 1724  cm^{-1} was assigned to redshifted C=O stretch within neutral clusters. We show that zwitter ion formation of amino acids readily occurs with only few water molecules available even under extreme conditions.

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

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.

Effective Two-Body Interactions

Cooperative or nonadditive effects contribute to the pairwise noncovalent interaction of two molecules in a cluster or the condensed phase in ways that depend on the specific arrangements and interactions of the other surrounding molecules that constitute their environment. General expressions for an effective two-body interaction are presented, which are correct to increasing orders in the many-body expansion. The simplest result, correct through third order, requires only seven individual calculations, in contrast to a linear number of three-body contributions. Two applications are presented. First, an error analysis is performed on a model (H₂O)₈ cluster which completes the first solvation shell of a central water-water hydrogen bond. Energy decomposition analysis is performed to show that the largest effects of cooperativity on the central hydrogen bond arise from electrical polarization. Second, the nature of cooperative effects on proton transfer in an HCl + (H₂O)₄ cluster is characterized.

Unveiling Zwitterionization of Glycine in the Microhydration Limit

Charge separation under solvation stress conditions is a fundamental process that comes in many forms in doped water clusters. Yet, the mechanism of intramolecular charge separation, where constraints due to the molecular structure might be intricately tied to restricted solvation structures, remains largely unexplored. Microhydrated amino acids are such paradigmatic molecules. Ab initio simulations are carried out at 300 K in the frameworks of metadynamics sampling and thermodynamic integration to map the thermal mechanisms of zwitterionization using Gly(H₂O) _n with n = 4 and 10. In both cases, a similar water-mediated proton transfer chain mechanism is observed; yet, detailed analyses of thermodynamics and kinetics demonstrate that the charge-separated zwitterion is the preferred species only for n = 10 mainly due to kinetic stabilization. Structural analyses disclose that bifurcated H-bonded water bridges, connecting the cationic and anionic sites in the fluctuating microhydration network at room temperature, are enhanced in the transition-state ensemble exclusively for n = 10 and become overwhelmingly abundant in the stable zwitterion. The findings offer potential insights into charge separation under solvation stress conditions beyond the present example.

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid-base chemistry depending on temperature in small water assemblies.

Slowing Down of the Molecular Reorientation of Water in Concentrated Alkaline Solutions

It is generally accepted that the hydroxide ion (OH^-) is a strong hydrogen bond acceptor and that its anomalously high diffusion constant in water results from a Grotthuss-like structural diffusion mechanism. However, the spatial extent over which OH^- ions influence the dynamics of the hydrogen-bond network of water remained largely unclear. Here, we measure the ultrafast dynamics of OH groups of HDO molecules interacting with the deuterated hydroxide ion OD^-. For solutions with OD^- concentrations up to 4 M, we find that HDO molecules that are not directly interacting with the ions have a reorientation time constant of ∼2.7 ps, similar to that of pure liquid water. When the concentration of OD^- ions is increased, the reorientation time constant increases, indicating a strong slowing down of the structural dynamics of the solution.

Ultrafast photoinduced dynamics of single atoms solvated inside helium nanodroplets

Helium nanodroplets can serve as reaction containers for photoinduced time-resolved studies of cold, isolated molecular systems that are otherwise inaccessible. Recently, three different dynamical processes, triggered by photoexcitation of a single atom inside a droplet, were observed in their natural time scale: Expansion of the He solvation shell (He bubble) within 600 fs initiates a collective bubble oscillation with a ∼30 ps oscillation period, followed by dopant ejection after ∼60 ps. Here, we present a systematic investigation of these processes by combining time-resolved photoelectron and photoion spectroscopy with time-dependent He density functional theory simulations. By variation of the photoexcitation energy, we find that the full excess excitation energy, represented by the blue-shifted in-droplet excitation band, is completely transferred to the He environment during the bubble expansion. Surprisingly, we find that variation of the droplet size has only a minor influence on the ejection time, providing insight into the spatial distribution of the ground-state atoms before photoexcitation. Simulated particle trajectories after photoexcitation are in agreement with experimental observations and suggest that the majority of ground-state atoms are located at around 16 Å below the droplet surface. Bubble expansion and oscillation are purely local effects, depending only on the ultimate dopant environment. These solvation-induced dynamics will be superimposed on intramolecular dynamics of molecular systems, and a mechanistic description is fundamental for the interpretation of future experiments.

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.

Ion Pairing in HCl-Water Clusters: From Electronic Structure Investigations to Multiconfigurational Force-Field Development

In the bulk, condensed-phase HCl exists as a dissociated Cl^- ion and a proton that is delocalized over solvating water molecules. However, in the gas phase, HCl is covalent, and even on the introduction of hydrating water molecules, the HCl covalent state dominates small clusters and is relevant at larger clusters including 21 water molecules. Electronic structure calculations (at the MP2 level) and ab initio metadynamics simulations (at the DFT level) have been carried out on HCl-(H₂O)_n clusters with n = 2-22 to investigate distinct solvation environments in clusters from covalent HCl structure, to contact ion pairs and solvent-separated ion pairs. The data were further used to train and validate a multiconfigurational force-field for HCl-water clusters that incorporates covalent HCl states into the MS-EVB3.2 formalism. Additionally, the many-body interaction of the Cl^- ion with water and the excess proton was modeled by the introduction of two geometric three-body terms that incorporates the dominant many-body interaction in an efficient noniterative manner.

Quantum nature of the hydrogen bond from ambient conditions down to ultra-low temperatures

Quantum simulations reveal strong temperature effects for weak hydrogen bonds and differences in quantum delocalization between various hydrogen-bonded systems.