Microsolvation Clusters Research Articles

Structure and Thermodynamics of Li+Arn Clusters beyond the Second Solvation Shell

Small Li+Arn clusters are employed in this work as model systems to study microsolvation. Although first and second solvation shells are expected to be the most relevant ones for this type of atomic solvents, it is also interesting to explore larger clusters in order to identify the influence of external atoms on structural and thermodynamic properties. In this work, we perform a global geometry optimization for Li+Arn clusters (with n = 41–100) and parallel tempering Monte Carlo (PTMC) simulations for some selected sizes. The results show that global minimum structures of large clusters always have 6 argon atoms in the first solvation shell while maintaining the number of 14 or 16 argon atoms in the second one. By contrast, third and fourth solvation shells vary significantly the number of argon atoms with the cluster size, and other shells can hardly be assigned due to the reduced influence of Li+ on the external argon atoms for large clusters. In turn, PTMC calculations show that the melting of the most external solvation shells of large microsolvation clusters occurs at T∼50K, which is independent of cluster size. Structural transitions can be observed between quasi-degenerated structures at low temperatures. Moreover, the present results highlight the fluxional character of the external solvation shells of these large Li+Arn clusters, which may be seen as typical “snowball” structures.

Conformations of borneol and isoborneol in the gas phase: Their monomers and microsolvation clusters.

Borneol is a natural monoterpene with significant applications in various industries, including medicine and perfumery. It presents several diastereomers with different physical and chemical properties, influenced by their unique structures and interactions with molecular receptors. However, a complete description of its inherent structure and solvent interactions remains elusive. Here, we report a detailed investigation of the gas-phase experimental structures of borneol and isoborneol, along with the description of their microsolvation complexes with the common solvents water and dimethyl sulfoxide. The molecules and complexes were studied using chirped-pulse Fourier transform microwave spectroscopy coupled to a supersonic expansion source. Although three rotamers are potentially populated under the conditions of the supersonic expansion, only one of them was observed for each monomer. The examination of the monohydrated complexes revealed structures stabilized by hydrogen bonds and non-covalent C-H⋯O interactions, with water as the hydrogen bond donor. Interestingly, in the clusters with dimethyl sulfoxide, borneol and isoborneol change their roles acting as donors. We further identified a higher-energy rotamer of the borneol monomer in one of the complexes with dimethyl sulfoxide, while that rotamer was missing in the experiment for the monomer. This observation is not common and highlights a specific position in borneol especially favorable for forming stable complexes, which could have implications in the understanding of the unique physical and chemical properties of the diastereomers.

Effect of charge and solvation shell on non-radiative decay processes in s-block cationic metal ion water clusters.

Intermolecular Coulombic decay or electron transfer-mediated decay are the autoionization processes through which a molecule can relax. This relaxation is only possible if the inner valence's ionization potential (IP) exceeds the system's double ionization potential (DIP). To study the effects of charge and solvation shell, we have calculated the IP, DIP values, and lifetime of Na-2s and Mg-2s temporary bound states in various optimized structures of Na+-(H2O)n and Mg2+-(H2O)n (n = 1-5) micro-solvated clusters, where n water molecules are distributed in a way that some are directly bound to the metal ion and the rest to the water molecules. The first and second solvation shells are the names for the former and the latter water-binding positions, respectively. For a given n, the lifetime of decaying states is longer when water molecules are in the second solvation shell. We found that the Mg-2p state can decay for all n values in Mg2+-(H2O)n clusters, whereas the Na-2p state's decay is possible for n ≥ 2 in Na+-(H2O)n clusters. Our findings highlight the influence of metal ions' charge, different solvation shell structures, and the number of water molecules on the decay rate. These systems are relevant to the human body, which makes this study significant.

Structural Sampling and Solvation Models for the Simulation of Electronic Spectra: Pyrazine as a Case Study

The impact of sampling methods on spectral broadening in the gas phase and on the convergence of spectra in aqueous solution when using microsolvation, continuum solvation, and hybrid models is studied using pyrazine as a test case. For the sake of comparing classical Maxwell-Boltzmann and Wigner samplings in the gas phase, static and time-resolved X-ray absorption spectra after photoexcitation to the lowest 1B2u(ππ*) state, as well as the static UV-vis absorption spectrum, are considered. In addition, the UV-vis absorption spectrum of pyrazine in aqueous solution is also computed in order to systematically investigate its convergence with the number of explicitly included solvent shells with and without taking bulk solvation effects into account with the conductor-like screening model to represent implicit water beyond such explicit solute complexes. Concerning the static and time-resolved X-ray absorption spectra of pyrazine at the carbon K-edge as well as its UV-vis absorption spectrum in the gas phase, we find that these spectra obtained with Wigner and Maxwell-Boltzmann samplings substantially agree. For the UV-vis absorption spectrum in the aqueous solution, only the first two energetically low-lying bands converge quickly with the size of the explicitly included solvation shells, either with or without an additional continuum solvation taken into account. In stark contrast, calculations of the higher-lying excitations relying on finite microsolvated clusters without additional continuum solvation severely suffer from unphysical charge-transfer excitations into Rydberg-like orbitals at the cluster/vacuum interface. This finding indicates that computational UV-vis absorption spectra covering sufficiently high-lying states converge only if continuum solvation of the explicitly microsolvated solutes is included in the models.

Decay Processes in Cationic Alkali Metals in Microsolvated Clusters: A Complex Absorbing Potential Based Equation-of-Motion Coupled Cluster Investigation.

We have employed the highly accurate complex absorbing potential based ionization potential equation-of-motion coupled cluster singles and doubles (CAP-IP-EOM-CCSD) method to study the various intermolecular decay processes in ionized metals (Li+, Na+, K+) microsolvated by water molecules. For the Li atom, the electron is ionized from the 1s subshell. However, for Na and K atoms, the electron is ionized from 2s and both 2s and 2p subshells, respectively. We have investigated decay processes for the Li+-(H2O)n (n = 1-3) systems, as well as Na+-(H2O)n (n = 1, 2), and K+-H2O. The lithium cation in water can decay only via electron transfer mediated decay (ETMD) as there are no valence electrons in lithium. We have investigated how the various decay processes change in the presence of different alkali metal atoms and how the increasing number of water molecules play a significant role in the decay of microsolvated systems. To see the effect of the environment, we have studied Li+-NH3 in comparison to Li+-H2O. In the case of Na+-H2O, we have studied the impact of bond distance on the decay width. The effect of polarization on decay width was checked for the X+-H2O (X = Li, Na) systems. We used the PCM model to study the polarization effect. We have compared our results with existing theoretical and experimental results wherever available in the literature.

UV–Vis, Raman spectroscopic and density functional theoretical studies on microsolvation 1, 2, 4-triazole-3-thione clusters

Photophysical and photochemical reactions of microsolvation clusters are attracting an increasing attention due to their wide applications in materials science and biology. In this paper, 1, 2, 4-triazole-3-thione (3TT) is investigated in solid, protic, and aprotic solvents using FT-Raman, resonance Raman and electronic absorption spectroscopic experiments. The structures of microsolvation clusters in solvents were confirmed by 488 nm Raman spectroscopy combining with density functional theory (DFT) calculation. Steady-state absorption and resonance Raman spectra of 3TT in different environments indicate that the intermolecular hydrogen bonding may reveal important insights in the photophysical and photochemical process. With the aid of DFT and time-dependent density functional theory (TDDFT) calculations, we assigned the observed Raman spectra to the microsolvation clusters in acetonitrile, water and methanol, and carried out preliminary investigations on spectrum shifts of UV and Raman spectra due to the hydrogen bonding with the solvent molecules. The intermolecular >NH···O and >=S···H hydrogen bonding interactions, which are the key constituents of stable thione structure of 3TT, revealed the obvious spectrum shifts of 3TT, including Raman and absorption shifts in CH3CN, CH3OH and H2O. The hydrogen bond sites were further confirmed to be located on the functional group SCNH of 3TT with CH3CN, H2O and CH3OH.

Modeling microsolvation clusters with electronic-structure calculations guided by analytical potentials and predictive machine learning techniques.

We propose a new methodology to study, at the density functional theory (DFT) level, the clusters resulting from the microsolvation of alkali-metal ions with rare-gas atoms. The workflow begins with a global optimization search to generate a pool of low-energy minimum structures for different cluster sizes. This is achieved by employing an analytical potential energy surface (PES) and an evolutionary algorithm (EA). The next main stage of the methodology is devoted to establish an adequate DFT approach to treat the microsolvation system, through a systematic benchmark study involving several combinations of functionals and basis sets, in order to characterize the global minimum structures of the smaller clusters. In the next stage, we apply machine learning (ML) classification algorithms to predict how the low-energy minima of the analytical PES map to the DFT ones. An early and accurate detection of likely DFT local minima is extremely important to guide the choice of the most promising low-energy minima of large clusters to be re-optimized at the DFT level of theory. In this work, the methodology was applied to the Li+Krn (n = 2-14 and 16) microsolvation clusters for which the most competitive DFT approach was found to be the B3LYP-D3/aug-pcseg-1. Additionally, the ML classifier was able to accurately predict most of the solutions to be re-optimized at the DFT level of theory, thereby greatly enhancing the efficiency of the process and allowing its applicability to larger clusters.

Microsolvation of Histidine—A Theoretical Study of Intermolecular Interactions Based on AIM and SAPT Approaches

Histidine is unique among amino acids because of its rich tautomeric properties. It participates in essential enzymatic centers, such as catalytic triads. The main aim of the study is the modeling of the change of molecular properties between the gas phase and solution using microsolvation models. We investigate histidine in its three protonation states, microsolvated with 1:6 water molecules. These clusters are studied computationally, in the gas phase and with water as a solvent (Polarizable Continuum Model, PCM) within the Density Functional Theory (DFT) framework. The structural analysis reveals the presence of intra- and intermolecular hydrogen bonds. The Atoms-in-Molecules (AIM) theory is employed to determine the impact of solvation on the charge flow within the histidine, with emphasis on the similarity of the two imidazole nitrogen atoms—topologically not equivalent, they are revealed as electronically similar due to the heterocyclic ring aromaticity. Finally, the Symmetry-Adapted Perturbation Theory (SAPT) is used to examine the stability of the microsolvation clusters. While electrostatic and exchange terms dominate in magnitude over polarization and dispersion, the sum of electrostatic and exchange term is close to zero. This makes polarization the factor governing the actual interaction energy. The most important finding of this study is that even with microsolvation, the polarization induced by the presence of implicit solvent is still significant. Therefore, we recommend combined approaches, mixing explicit water molecules with implicit models.

Application of the Approximate 3D-Reference Interaction Site Model (RISM) Molecular Solvation Theory to Acetonitrile as Solvent.

The integral equation formalism based reference interaction site model (RISM) molecular solvation theory is applied to pure liquid acetonitrile and water-acetonitrile binary mixtures of different compositions. Solvate formation of d- and f-block ions by ACN is also calculated to check applicability of the RISM theory. The generalized Amber force field (GAFF) and the universal force field (UFF) parameters were found to be suitable for applications with the RISM theory for acetonitrile solvent. The presence of local microsolvated clusters for water-acetonitrile mixtures is validated in the RISM calculations.

Aqueous Micro-hydration of Na+(H2O)n=1-7 Clusters: DFT Study

AbstractSodium ion micro-solvated clusters, [Na(H2O) n]+, n = 1–7, were completed by (DFT) density functional theory at B3LYP/6-311+G(d,p) level in the gaseous phase. At the ambient situation, the four, five and six micro-solvated configurations can convert from each other. The investigation of the sequential water binding energy on Na+ obviously indicates that the influence of Na+ on the neighboring water molecules goes beyond the first solvation layer with the hydration number of 5. The hydration number of Na+ is 5 and the hydration space (rNa-O) is 2.43 Å. The current study displays that all our simulations have an brilliant harmony with the diffraction result from X-ray scattering study. The vibration frequency of H2O solvent was also determined. This work is important for additional identification of the Na+(H2O)n clusters in aqueous medium.

Microsolvation of Li+ in a Mixture of Argon and Krypton: Unveiling the Most Stable Structures of the Clusters.

The microsolvation of Li+ by both argon and krypton atoms has been studied based on a new potential energy surface that includes two- and three-body interactions; the potential terms involving the lithium ion were calibrated with CCSD(T)/aug-cc-pVQZ energies after being corrected for the basis-set superposition error. The structures of the Li+Ar nKr m ([Formula: see text]) clusters arising from global optimization show a first solvation shell preferentially occupied by krypton atoms. These binary-solvent microsolvation clusters are most stable when the total number of krypton (argon) atoms occupy the first (second) solvation shell.

X-Ray absorption spectra of microsolvated metal cations.

Core excited states in clusters or bulk medium are known to undergo a process of internal ionisation, whereby the excited electron delocalises throughout the medium. This delocalisation is visible in the shifting and broadening of lines in X-ray absorption spectra, and it impacts the electronic decay initiated by photoabsorption. In this paper we study the delocalisation of electrons excited from the 1s core orbital of Na(+) and Mg(2+) ions in microsolvated Na(+)(H2O)m and Mg(2+)(H2O)m clusters (m = 1-6) by computing the X-ray absorption spectra and electron distributions in different core excited states. We show that addition of water ligands to the ion leads to more and more pronounced delocalisation of the core-to-valence 1s → 3p and core-to-Rydberg 1s → 4p excitations. Even for the compact 1s → 3p excitation the excited electron is mostly located on the water molecules when the solvation shell is complete. We also found that the degree of delocalisation strongly depends on the cluster geometry and the ionic charge. These results indicate that even in small microsolvated clusters delocalisation of core excited electrons is substantial and will affect the following electronic decay. The accuracy and transferability of our results are corroborated by the good agreement between our XAS spectra of microsolvated Na(+) and experimental X-ray absorption spectra of dilute NaCl solutions.

Computational comparison of oxidation stability: Solvent/salt monomers vs solvent–solvent/salt pairs

A fundamental understanding of the anodic stabilities of electrolytes is important for the development of advanced high-voltage electrolytes. In this study, we calculated and systematically compared the oxidation stabilities of monomeric solvents and anions, and bimolecular solvent–solvent and anion–solvent systems that are considered to be high-voltage electrolyte components, using ab initio calculations. Oxidation stabilities of solvent or anion monomers without considering specific solvation molecules cannot represent experimental oxidation stabilities. The oxidation of electrolytes usually forms neutral or cationic radicals, which immediately undergo further reactions stabilizing the products. Oxidatively driven intermolecular reactions are the main reason for the lower oxidation stabilities of electrolytes compared with those of monomeric compounds. Electrolyte components such as tetramethylene sulfone (TMS), ethyl methyl sulfone (EMS), bis(oxalate)borate (BOB−), and bis(trifluoromethane)sulfonamide (TFSI−) that minimize such intermolecular chemical reactions on oxidation can maintain the oxidation stabilities of monomers. In predictions of the theoretical oxidation stabilities of electrolytes, simple comparisons of highest occupied molecular orbital energies can be misleading, even if microsolvation or bulk clusters are considered. Instead, bimolecular solvent complexes with a salt anion should be at least considered in oxidation calculations. This study provides important information on fundamental and applied aspects of the development of electrolytes.

Microsolvation of the Potassium Ion with Aromatic Rings: Comparison between Hexafluorobenzene and Benzene

We employ a recently developed methodology to study structural and energetic properties of the first solvation shells of the potassium ion in nonpolar environments due to aromatic rings, which is important to understand the selectivity of several biochemical phenomena. Our evolutionary algorithm is used in the global optimization study of clusters formed of K(+) solvated with hexafluorobenzene (HFBz) molecules. The global intermolecular interaction for these clusters has been decomposed in HFBz-HFBz and in K(+)-HFBz contributions, using a potential model based on different decompositions of the molecular polarizability of hexafluorobenzene. Putative global minimum structures of microsolvation clusters up to 21 hexafluorobenzene molecules were obtained and compared with the analogous K(+)-benzene clusters reported in our previous work (J. Phys. Chem. A 2012, 116, 4947-4956). We have found that both K(+)-(Bz)n and K(+)-(HFBz)n clusters show a strong magic number around the closure of the first solvation shell. Nonetheless, all K(+)-benzene clusters have essentially the same first solvation shell geometry with four solvent molecules around the ion, whereas the corresponding one for K(+)-(HFBz)n is completed with nine HFBz species, and its structural motif varies as n increases. This is attributed to the ion-solvent interaction that has a larger magnitude for K(+)-Bz than in the case of K(+)-HFBz. In addition, the ability of having more HFBz than Bz molecules around K(+) in the first solvation shell is intimately related to the inversion in the sign of the quadrupole moment of the two solvent species, which leads to a distinct ion-solvent geometry of approach.

Calculation of the Gibbs free energy of solvation and dissociation of HCl in water via Monte Carlo simulations and continuum solvation models

The Gibbs free energy of solvation and dissociation of hydrogen chloride in water is calculated through a combined molecular simulation/quantum chemical approach at four temperatures between T = 300 and 450 K. The Gibbs free energy is first decomposed into the sum of two components: the Gibbs free energy of transfer of molecular HCl from the vapor to the aqueous liquid phase and the standard-state Gibbs free energy of acid dissociation of HCl in aqueous solution. The former quantity is calculated using Gibbs ensemble Monte Carlo simulations using either Kohn-Sham density functional theory or a molecular mechanics force field to determine the system's potential energy. The latter Gibbs free energy contribution is computed using a continuum solvation model utilizing either experimental reference data or micro-solvated clusters. The predicted combined solvation and dissociation Gibbs free energies agree very well with available experimental data.

Alkali-Ion Microsolvation with Benzene Molecules

The target of this investigation is to characterize by a recently developed methodology, the main features of the first solvation shells of alkaline ions in nonpolar environments due to aromatic rings, which is of crucial relevance to understand the selectivity of several biochemical phenomena. We employ an evolutionary algorithm to obtain putative global minima of clusters formed with alkali-ions (M(+)) solvated with n benzene (Bz) molecules, i.e., M(+)-(Bz)(n). The global intermolecular interaction has been decomposed in Bz-Bz and in M(+)-Bz contributions, using a potential model based on different decompositions of the molecular polarizability of benzene. Specifically, we have studied the microsolvation of Na(+), K(+), and Cs(+) with benzene molecules. Microsolvation clusters up to n = 21 benzene molecules are involved in this work and the achieved global minimum structures are reported and discussed in detail. We observe that the number of benzene molecules allocated in the first solvation shell increases with the size of the cation, showing three molecules for Na(+) and four for both K(+) and Cs(+). The structure of this solvation shell keeps approximately unchanged as more benzene molecules are added to the cluster, which is independent of the ion. Particularly stable structures, so-called "magic numbers", arise for various nuclearities of the three alkali-ions. Strong "magic numbers" appear at n = 2, 3, and 4 for Na(+), K(+), and Cs(+), respectively. In addition, another set of weaker "magic numbers" (three per alkali-ion) are reported for larger nuclearities.

Opposite effects of successive hydration shells on the aqua ion structure of metal cations

The conformation adopted by water molecules around a metal ion in aqueous solution can be approached from a discrete representation of the solvent by studying the structure and stability of microsolvation clusters including a different number of water molecules. In this contribution, we present a quantum mechanics and ab initio molecular dynamics study on how the arrangement of water molecules around and shows preferentially a T h or S 6 symmetry depending on the number of hydration shells explicitly included in the calculation. The behaviour observed for both cases is the following: (1) for a cluster composed of the metal ion plus a first hydration shell, the preferred geometry shows a T h symmetry; (2) when a second hydration shell is added, a minimum can be found with an S 6 disposition of water molecules around the metal centre; (3) however, if a third hydration shell is considered a T h arrangement of water molecules around the central ion is strongly favoured. The structures of the different solvation complexes studied in this work are characterised. The preference for a T h or S 6 distribution of water ligands is rationalised in terms of total energies, interaction energies and the hydrogen bond network.

Electronic decay following ionization of aqueous Li+ microsolvation clusters

Ionization in the energetical range between 35 eV and 75 eV of aqueous Li(+) microsolvation clusters may initialize several different electronic decay processes. Electronic decay following H(2)O 2s ionization in a cationic cluster is reported. Li ionization probes the efficiency of electron transfer mediated decay (ETMD) processes. We report estimated ETMD lifetimes in the range of 20-100 fs for clusters with one to five water monomers. Furthermore, tertiary electron emission may occur via a combined cascade of electron transfer mediated decay and intermolecular Coulombic decay.

Microsolvation of Li+in Water Analyzed by Ionization and Double Ionization

The microsolvation of Li+ in water is investigated. The ionization and double-ionization spectra of the series of (H2O)nLi+ (n = 1−5) clusters are calculated ab initio by Green's function methods and discussed in detail. The impact of the solvation of the lithium cation with an increasing number of water molecules on the spectral characteristics is revealed. In the context of microsolvation, the discussion of the results for the (H2O)5Li+ cluster is particularly important because the second solvation shell is accessed. Ionization- and double-ionization spectra may be considered to be very useful for the study of microsolvation clusters.

Dodecahedral clathrate structures and magic numbers in alkali cation microhydration clusters.

Using global geometry optimization based on our specialized version of Genetic Algorithms, we have examined the global and most important local minimum energy structures of water microsolvation clusters of potassium and cesium cations within the common TIP4P/OPLS model. Together with our earlier results on the corresponding sodium case, this work constitutes a first step towards a theoretical elucidation of "magic numbers" of solvating molecules and proposed special structures occurring in these systems. In particular, the actual role of dodecahedral cage structures is examined. Within the present model, they do not occur in sodium microsolvation, in agreement with the absence of the magic number 20 for this system. For potassium and cesium microsolvation, dodecahedral cages do occur but their actual structures are far from ideal and their importance appears to be overrated. We offer simple explanations for structural features and trends, and for magic numbers smaller than 20.