Search

The Martini model, a coarse-grained forcefield for biomolecular simulations, has experienced a vast increase in popularity in the past decade. Its building-block approach balances computational efficiency with high chemical specificity, enabling the simulation of organic and inorganic molecules. The modeling of coarse-grained beads as Lennard-Jones particles poses challenges for the accurate reproduction of liquid–vapor interfacial properties, which are crucial in various applications, especially in the case of water. The latest version of the forcefield introduces refined interaction parameters for water beads, tackling the well-known artifact of Martini water freezing at room temperature. In addition, multiple sizes of water beads are available for simulating the solvation of small cavities, including the smallest pockets of proteins. This work focuses on studying the interfacial properties of Martini water, including surface tension and surface thickness. Employing the test-area method, we systematically compute the liquid–vapor surface tension across various combinations of water bead sizes and for temperatures from 300 to 350 K. These findings are of interest to the Martini community as they allow users to account for the low interfacial tension of Martini water by properly adjusting observables computed via coarse-grained simulations to allow for accurate matching against all-atom or experimental results. Surface tension data are also interpreted in terms of local enrichment of the various mixture components at the liquid–vapor interface by means of Gibbs' adsorption formalism. Finally, the critical scaling of the Martini surface tension with temperature is reported to be consistent with the critical exponent of the 3D Ising universality class. [ABSTRACT FROM AUTHOR]

The probability of finding a spherical "hole" of a given radius r contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions GV(r), have been well studied for hard-sphere fluids in d-dimensional Euclidean space R d . However, little is known about these functions for hard-sphere crystals for values of r beyond the hard-sphere diameter, as large holes are extremely rare in crystal phases. To overcome these computational challenges, we introduce a biased-sampling scheme that accurately determines hole statistics for equilibrium hard spheres on ranges of r that far extend those that could be previously explored. We discover that GV(r) in crystal and hexatic states exhibits oscillations whose amplitudes increase rapidly with the packing fraction, which stands in contrast to GV(r) that monotonically increases with r for fluid states. The oscillations in GV(r) for 2D crystals are strongly correlated with the local orientational order metric in the vicinity of the holes, and variations in GV(r) for 3D states indicate a transition between tetrahedral and octahedral holes, demonstrating the power of GV(r) as a probe of local coordination geometry. To further study the statistics of interparticle spacing in hard-sphere systems, we compute the local packing fraction distribution f(ϕl) of Delaunay cells and find that, for d ≤ 3, the excess kurtosis of f(ϕl) switches sign at a certain transitional global packing fraction. Our accurate methods to access hole statistics in hard-sphere crystals at the challenging intermediate length scales reported here can be applied to understand the important problem of solvation and hydrophobicity in water at such length scales. [ABSTRACT FROM AUTHOR]

Measurement of time evolution of fluorescence of a probe solute has been a quintessential technique to quantify how dipolar solvent molecules dynamically minimize the free energy of an electronically excited probe. During such solvation dynamics in bulk liquids, a substantial part of relaxation was shown to complete within sub-100 fs from time-gated fluorescence measurements, as also predicted by molecular dynamics simulation studies. However, equivalent quantification of solvation timescales by femtosecond pump–probe and broadband fluorescence measurements revealed an exponential nature of this initial relaxation having quite different timescales. Here, we set out to unveil the reason behind these puzzling contradictions. We introduce a method for estimating probe wavelength-dependent instrument response and demonstrate that the observation of the Gaussian vs exponential nature of initial relaxation is indeed dependent on the method of data analysis. These findings call for further experimental investigation and parallel development of theoretical models to elucidate the molecular-level mechanism accounting for different types of early time solvation. [ABSTRACT FROM AUTHOR]

We investigate the vertical electron attachment energy (VAE) of 1-methyl-4-nitroimidazole, a model radiosensitizer, employing quantum mechanics/molecular mechanics (QM/MM) and QM/polarized continuum (QM/PCM) solvation models. We considered the solvent-excluded surface (QM/PCM-SES) and Van der Waals (QM/PCM-VDW) cavities within the PCM framework, the electrostatic embedding QM/MM (EE-QM/MM) model, and the self-consistent sequential QM/MM polarizable electrostatic embedding (scPEE-S-QM/MM) model. Due to slow VAE convergence concerning the number of QM solvent molecules, full QM calculations prove inefficient. Ensemble averages in these calculations do not align with VAEs computed for the representative solute–solvent configuration. QM/MM and QM/PCM calculations show agreement with each other for sufficiently large QM regions, although the QM/PCM-VDW model exhibits artifacts linked to the cavity. QM/MM models demonstrate good agreement between ensemble averages and VAEs calculated with the representative configuration. Notably, the VAE computed with the scPEE-S-QM/MM model achieves faster convergence concerning the number of QM water molecules compared to the EE-QM/MM model, attributed to enhanced efficiency from MM charge polarization in the scPEE-S-QM/MM approach. This emphasizes the importance of QM/classical models with accurate solute–solvent and solvent–solvent mutual polarization for obtaining converged VAEs at a reasonable computational cost. The full-QM approach is very inefficient, while the microsolvation model is inaccurate. Computational savings in QM/MM models result from electrostatic embedding and the representative configuration, with the scPEE-S-QM/MM approach emerging as an efficient tool for describing bulk-solvated anions within the QM/MM framework. Its potential extends to improving transient anion state descriptions in biomolecules and radiosensitizers, especially given the frequent employment of microsolvation models. [ABSTRACT FROM AUTHOR]

The curvature of soft interfaces plays a crucial role in determining their mechanical and thermodynamic properties, both at macroscopic and microscopic scales. In the case of air/water interfaces, particular attention has recently focused on water microdroplets, due to their distinctive chemical reactivity. However, the specific impact of curvature on the molecular properties of interfacial water and interfacial reactivity has so far remained elusive. Here, we use molecular dynamics simulations to determine the effect of curvature on a broad range of structural, dynamical, and thermodynamical properties of the interface. For a droplet, a flat interface, and a cavity, we successively examine the structure of the hydrogen-bond network and its relation to vibrational spectroscopy, the dynamics of water translation, rotation, and hydrogen-bond exchanges, and the thermodynamics of ion solvation and ion-pair dissociation. Our simulations show that curvature predominantly impacts the hydrogen-bond structure through the fraction of dangling OH groups and the dynamics of interfacial water molecules. In contrast, curvature has a limited effect on solvation and ion-pair dissociation thermodynamics. For water microdroplets, this suggests that the curvature alone cannot fully account for the distinctive reactivity measured in these systems, which are of great importance for catalysis and atmospheric chemistry. [ABSTRACT FROM AUTHOR]

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al3+) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al–O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet. [ABSTRACT FROM AUTHOR]

Redox potentials depend on the nature of the solvent/electrolyte through the solvation energies of the ionic solute species. For concentrated electrolytes, ion solvation may deviate significantly from the Born model predictions due to ion pairing and correlation effects. Recently, Ghorai and Matyushov [J. Phys. Chem. B 124, 3754–3769 (2020)] predicted, on the basis of linear response theory, an anomalous trend in the solvation energies of room temperature ionic liquids, with deviations of hundreds of kJ/mol from the Born model for certain size solutes/ions. In this work, we computationally evaluate ionic solvation energies in the prototypical ionic liquid, 1-butyl-3-methylimidazolium tetrafluoroborate (BMIM/BF4), to further explore this behavior and benchmark several of the approximations utilized in the solvation energy predictions. For comparison, we additionally compute solvation energies within acetonitrile and molten NaCl salt to illustrate the limiting behavior of purely dipolar and ionic solvents. We find that the overscreening effect, which results from the inherent charge oscillations of the ionic liquid, is substantially reduced in magnitude due to screening from the dipoles of the molecular ions. Therefore, for the molten NaCl salt, for which the ions do not have permanent dipoles, modulation of ionic solvation energies from the overscreening effect is most significant. The conclusion is that ionic liquids do indeed exhibit unique solvation behavior due to peak(s) in the electrical susceptibility caused by the ion shell structure; redox potential shifts for BMIM/BF4 are of more modest order ∼0.1 V, but may be larger for other ionic liquids that approach molten salt behavior. [ABSTRACT FROM AUTHOR]

Numerous electrochemistry reactions require the precise calculation of the ion solvation energy. Despite the significant progress in the first-principles calculations for crystals and defect formation energies for solids, the liquid system free energy calculations still face many challenges. Ion solvation free energies can be calculated via different semiempirical ways, e.g., using implicit solvent models or cluster of explicit molecule models; however, systematically improving these models is difficult due to their lack of a solid theoretical base. A theoretically sound approach for calculating the free energy is to use thermodynamic integration. Nevertheless, owing to the difficulties of self-consistent convergence in first-principles calculations for unphysical atomic configurations, the computational alchemy approach has not been widely used for first-principles calculations. This study proposes a general approach to use first-principles computational alchemy for calculating the ion solvation energy. This approach is also applicable for other small molecules. The calculated ion solvation free energies for Li+, Na+, K+, Be2+, Mg2+, and Ca2+ are close to the experimental results, and the standard deviation due to molecular dynamics fluctuations is within 0.06 eV. [ABSTRACT FROM AUTHOR]

The analysis of water density fluctuations in the hydration shell of nonpolar solutes provides insights into water-mediated interactions, especially hydrophobic interactions. These fluctuations are sensitive to small perturbations due to changes in thermodynamic conditions, such as temperature and pressure, but also to the presence of cosolutes, such as salts or small organic molecules. Herein, we investigate the effect of two classes of adsorbing cosolutes, using urea and methanol as representatives, on the fluctuations in energy and solvent density within the solvation shell of a model extended hydrophobic solute. We focus on the interactions of the cosolutes with the hydrophobic hydration shell, rather than with the solute itself, which though important remain largely unexplored. We calculate and analyze the interfacial partial molar energy of the cosolute, using a methodology based on the small system method. This approach provides correlated solvent density and energy fluctuations and allows us to decompose them into contributions due to interactions between the different components present in the solvation shell of the solute. The results show that adsorbed urea molecules interact more favorably with water than nonadsorbed urea molecules, which leads to the attenuation of interfacial density fluctuations and thus to the stabilization of the solvation shell. By contrast, the adsorbed methanol molecules interact preferably with other methanol molecules in the solvation shell, leading to a nano-phase segregated structure, which enhances interfacial fluctuations. [ABSTRACT FROM AUTHOR]

Microsolvation of the carbamate moiety delivers precise information on complexation effects on the N–C=O backbone and is of relevance to the peptide bond functionality. In this context, the mono-, di-, and trihydrated complexes of methyl carbamate have been studied in molecular expansion by high-resolution microwave spectroscopy, using chirped-pulse and Fabry–Perot resonator Fourier transform microwave instruments covering the frequency range from 2 to 18 GHz. From the rotational constants of the parent and the 18Ow substituted monoisotopologues, accurate values have been derived for the geometries of the hydrogen bond interactions. The nuclear quadrupole coupling constant χcc of the nitrogen nucleus provides a direct measure of complexation changes and decreases with the degree of hydration, whereas the hindered internal rotation barrier increases slightly with microsolvation. Both tendencies could have a common origin in the π-cooperative inductive effects as the microsolvation series progresses. All transitions are split by the internal rotation of the methyl top and the nuclear quadrupole coupling, and in the largest cluster, they are additionally split by an inversion motion. [ABSTRACT FROM AUTHOR]

The sinking of alkali cations in superfluid 4He nanodroplets is investigated theoretically using liquid 4He time-dependent density functional theory at zero temperature. The simulations illustrate the dynamics of the buildup of the first solvation shell around the ions. The number of helium atoms in this shell is found to linearly increase with time during the first stages of the dynamics. This points to a Poissonian capture process, as concluded in the work of Albrechtsen et al. on the primary steps of Na+ solvation in helium droplets [Albrechtsen et al., Nature 623, 319 (2023)]. The energy dissipation rate by helium atom ejection is found to be quite similar between all alkalis, the main difference being a larger energy dissipated per atom for the lighter alkalis at the beginning of the dynamics. In addition, the number of helium atoms in the first solvation shell is found to be lower at the end of the dynamics than at equilibrium for both Li+ and Na+, pointing to a kinetic rather than thermodynamical control of the snowball size for small and strongly attractive ions. [ABSTRACT FROM AUTHOR]

Current models to understand the reactivity of metal/aqueous interfaces in electrochemistry, e.g., volcano plots, are based on the adsorption free energies of reactants and products, which are often small hydrophobic molecules (such as in CO2 and N2 reduction). Calculations played a major role in the quantification and comprehension of these free energies in terms of the interactions that the reactive species form with the surface. However, solvation free energies also come into play in two ways: (i) by modulating the adsorption free energy together with solute–surface interactions, as the solute has to penetrate the water adlayer in contact with the surface and get partially desolvated (which costs free energy); (ii) by regulating transport across the interface, i.e., the free energy profile from the bulk to the interface, which is strongly non-monotonic due to the unique nature of metal/aqueous interfaces. Here, we use constant potential molecular dynamics to study the solvation contributions, and we uncover huge effects of the shape and orientation (on top of the already known size effect) of small hydrophobic and amphiphilic solutes on their adsorption free energy. We propose a minimal theoretical model, the S.O.S. model, that accounts for size, orientation, and shape effects. These novel aspects are rationalized by recasting the concepts at the base of the Lum–Chandler–Weeks theory of hydrophobic solvation (for small solutes in the so-called volume-dominated regime) into a layer-by-layer form, where the properties of each interfacial region close to the metal are explicitly taken into account. [ABSTRACT FROM AUTHOR]

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca 2 + , Cd 2 + , Cu 2 + , Fe 2 + , Hg 2 + , Mg 2 + , Ni 2 + , Pb 2 + , and Zn 2 + , we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M( H 2 O)2]2+ aqua complexes to [M( H 2 O) 9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M( H 2 O) n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion–water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data. [ABSTRACT FROM AUTHOR]

The entanglement between system and bath often plays a pivotal role in complex systems spanning multiple orders of magnitude. A system–bath entanglement theorem was previously established for Gaussian environments in J. Chem. Phys. 152, 034102 (2020) regarding linear response functions. This theorem connects the entangled responses to the local system and bare bath properties. In this work, we generalize it to correlation functions. Key steps in derivations involve using the generalized Langevin dynamics for hybridizing bath modes and the Bogoliubov transformation that maps the original finite-temperature reservoir to an effective zero-temperature vacuum by employing an auxiliary bath. The generalized theorem allows us to evaluate the system–bath entangled correlations and the bath mode correlations in the total composite space, as long as we know the bare-bath statistical properties and obtain the reduced system correlations. To demonstrate the cross-scale entanglements, we utilize the generalized theorem to calculate the solvation free energy of an electron transfer system with intramolecular vibrational modes. [ABSTRACT FROM AUTHOR]

The catalytic hydrogenation of biomass-derived chemicals is essential in chemical industry due to the growing demand for sustainable and renewable energy sources. In this study, we present a comprehensive theoretical investigation regarding the hydrogenation of glycolaldehyde to ethylene glycol over a Ru/C catalyst by employing density functional theory and ab initio molecular dynamics simulations. With inclusion of explicit solvation, we have demonstrated that the glycolaldehyde hydrogenation is significantly improved due to the fast proton transfer through the hydrogen bond network. The enhanced activity could be attributed to the participation of the solvent water as the hydrogen source and the highly positively charged state of a Ru cluster in an aqueous phase, which are critical for the activation of aldehyde groups and proton-assisted hydrogenation. Overall, our findings provide valuable insights into glycolaldehyde hydrogenation over Ru/C catalysts in the aqueous phase, highlighting the importance of solvation effects in the biomass conversion. [ABSTRACT FROM AUTHOR]

In this work, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory. [ABSTRACT FROM AUTHOR]

Understanding the influence of a crowded intracellular environment on the structure and solvation of DNA functionalized gold nanoparticles (ss-DNA AuNP) is necessary for designing applications in nanomedicine. In this study, the effect of single (Gly, Ser, Lys) and mixture of amino acids (Gly+Ser, Gly+Lys, Ser+Lys) at crowded concentrations is examined on the structure of the ss-DNA AuNP using molecular dynamics simulations. Using the structural estimators such as pair correlation functions and ligand shell positional fluctuations, the solvation entropy is estimated. Combining the AuNP–solvent interaction energy with the solvation entropy estimates, the free energy of solvation of the AuNP in crowded solutions is computed. The solvation entropy favours the solvation free energy which becomes more favourable for larger effective size of AuNP in crowded solutions relative to that in water. The effective size of AuNP depends on the different propensity of the crowders to adsorb on Au surface, with the smallest crowder (Gly) having the highest propensity inducing the least effective AuNP size as compared to other single crowder solutions. In mixed crowded solutions of amino acids of variable size and chemistry, distinctive local adsorption of the crowders on the gold surface is observed that controls the additive or non-additive crowding effects which govern an increase (in Gly+Ser) or decrease (in Gly+Lys) in nanoparticle effective size respectively. The results shed light into the fundamental understanding of the influence of intracellular crowding on structure of ss-DNA AuNP and plausible employability of crowding as a tool to design programmable self-assembly of functionalized nanoparticles. [ABSTRACT FROM AUTHOR]

Delocalization error constrains the accuracy of density functional theory in describing molecular interactions in ion–water systems. Using Na+ and Cl− in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion–water and water–water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular interactions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with many-body MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl− in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected MB formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy. [ABSTRACT FROM AUTHOR]

The solvent effects on the linear and second-order nonlinear optical properties of an aminonaphtylethenylpyridinium (ANEP) dye are investigated by combining experimental and theoretical chemistry methods. On the one hand, deep near infrared (NIR) hyper-Rayleigh scattering (HRS) measurements (1840–1950 nm) are performed on solutions of di-8-ANEPPS in deuterated chloroform, dimethylformamide, and dimethylsulfoxide to determine their first hyperpolarizablity (βHRS). For the first time, these HRS experiments are carried out in the picosecond regime in the deep NIR with very moderate (≤3 mW) average input power, providing a good signal-to-noise ratio and avoiding solvent thermal effects. Moreover, the frequency dispersion of βHRS is investigated for Disperse Red 1 (DR1), a dye commonly used as HRS external reference. On the other hand, these are compared with computational chemistry results obtained by using a sequential molecular dynamics (MD) then quantum mechanics (QM) approach. The MD method allows accounting for the dynamical nature of the molecular structures. Then, the QM part is based on TDDFT/M06-2X/6-311+G* calculations using solvation models ranging from continuum to discrete ones. Measurements report a decrease of the βHRS of di-8-ANEPPS in more polar solvents and these effects are reproduced by the different solvation models. For di-8-ANEPPS and DR1, comparisons show that the use of a hybrid solvation model, combining the description of the solvent molecules around the probe by point charges with a continuum model, already achieves quasi quantitative agreement with experiment. These results are further improved by using a polarizable embedding that includes the atomic polarizabilities in the solvent description. [ABSTRACT FROM AUTHOR]

Structural studies using x-ray scattering methods for investigating molecules in solution are shifting focus toward describing the role and effects of the surrounding solvent. However, forward models based on molecular dynamics (MD) simulations to simulate structure factors and x-ray scattering from interatomic distributions such as radial distribution functions (RDFs) face limitations imposed by simulations, particularly at low values of the scattering vector q. In this work, we show how the value of the structure factor at q = 0 calculated from RDFs sampled from finite MD simulations is effectively dependent on the size of the simulation cell. To eliminate this error, we derive a new scheme to renormalize the sampled RDFs based on a model of the excluded volume of the particle-pairs they were sampled from, to emulate sampling from an infinite system. We compare this new correction method to two previous RDF-correction methods, developed for Kirkwood–Buff theory applications. We present a quantitative test to assess the reliability of the simulated low-q scattering signal and show that our RDF-correction successfully recovers the correct q = 0 limit for neat water. We investigate the effect of MD-sampling time on the RDF-corrections, before advancing to a molecular example system, comprised of a transition metal complex solvated in a series of water cells with varying densities. We show that our correction recovers the correct q = 0 behavior for all densities. Furthermore, we employ a simple continuum scattering model to dissect the total scattering signal from the solvent–solvent structural correlations in a solute–solvent model system to find two distinct contributions: a non-local density-contribution from the finite, fixed cell size in NVT simulations, and a local contribution from the solvent shell. We show how the second contribution can be approximated without also including the finite-size contribution. Finally, we provide a "best-practices"-checklist for experimentalists planning to incorporate explicit solvation MD simulations in future work, offering guidance for improving the accuracy and reliability of structural studies using x-ray scattering methods in solution. [ABSTRACT FROM AUTHOR]

In a recent study, we developed a kinetic-energy density functional that can be utilized in orbital-free quantum mechanical/molecular mechanical (OF-QM/MM) simulations. The functional includes the nonlocal term constructed from the response function of the reference system of the QM solute. The present work provides a method to combine the OF-QM/MM with a theory of solutions based on the energy representation to compute the solvation free energy of the QM solute in solution. The method is applied to the calculation of the solvation free energy Δμ of a QM water solute in an MM water solvent. It is demonstrated that Δμ is computed as −7.7 kcal/mol, in good agreement with an experimental value of −6.3 kcal/mol. We also develop a theory to map the free energy δμ due to electron density polarization onto the coordinate space of electrons. The free energy density obtained by the free-energy mapping for the QM water clarifies that each hydrogen atom makes a positive contribution (+34.7 kcal/mol) to δμ, and the oxygen atom gives the negative free energy (−71.7 kcal/mol). It is shown that the small polarization free energy −2.4 kcal/mol is generated as a result of the cancellation of these counteracting energies. These analyses are made possible by the OF-QM/MM approach combined with a statistical theory of solutions. [ABSTRACT FROM AUTHOR]

The computational modeling of electrochemical interfaces and their applications in electrocatalysis has attracted great attention in recent years. While tremendous progress has been made in this area, however, the accurate atomistic descriptions at the electrode/electrolyte interfaces remain a great challenge. The Computational Hydrogen Electrode (CHE) method and continuum modeling of the solvent and electrolyte interactions form the basis for most of these methodological developments. Several posterior corrections have been added to the CHE method to improve its accuracy and widen its applications. The most recently developed grand canonical potential approaches with the embedded diffuse layer models have shown considerable improvement in defining interfacial interactions at electrode/electrolyte interfaces over the state-of-the-art computational models for electrocatalysis. In this Review, we present an overview of these different computational models developed over the years to quantitatively probe the thermodynamics and kinetics of electrochemical reactions in the presence of an electrified catalyst surface under various electrochemical environments. We begin our discussion by giving a brief picture of the different continuum solvation approaches, implemented within the ab initio method to effectively model the solvent and electrolyte interactions. Next, we present the thermodynamic and kinetic modeling approaches to determine the activity and stability of the electrocatalysts. A few applications to these approaches are also discussed. We conclude by giving an outlook on the different machine learning models that have been integrated with the thermodynamic approaches to improve their efficiency and widen their applicability. [ABSTRACT FROM AUTHOR]

In liquid water, there is a natural tendency to form aggregates that consist of water molecules linked by hydrogen bonds. Such spontaneously formed aggregates are surrounded by a "sea" of disordered water molecules, with both forms remaining in equilibrium. The process of creating water aggregates also takes place in the solvation water of proteins, but in this case, the interactions of water molecules with the protein surface shift the equilibrium of the process. In this paper, we analyze the structural properties of the solvation water in antifreeze proteins (AFPs). The results of molecular dynamics analysis with the use of various parameters related to the structure of solvation water on the protein surface are presented. We found that in the vicinity of the active region responsible for the binding of AFPs to ice, the equilibrium is clearly shifted toward the formation of "ice-like aggregates," and the solvation water has a more ordered ice-like structure. We have demonstrated that a reduction in the tendency to create "ice-like aggregates" results in a significant reduction in the antifreeze activity of the protein. We conclude that shifting the equilibrium in favor of the formation of "ice-like aggregates" in the solvation water in the active region is a prerequisite for the biological functionality of AFPs, at least for AFPs having a well-defined ice binding area. In addition, our results fully confirm the validity of the "anchored clathrate water" concept, formulated by Garnham et al. [Proc. Natl. Acad. Sci. U. S. A. 108, 7363 (2011)]. [ABSTRACT FROM AUTHOR]

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces' degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. Our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates. [ABSTRACT FROM AUTHOR]

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies. [ABSTRACT FROM AUTHOR]

The dynamics of a liquid and its coupling to a solute are crucial for a better understanding of chemical processes in the liquid phase. In isotropic and homogeneous solutions, the time-correlation function of a solute is expected to vanish over time due to the translational and diffusive motions of the solvent. The three-pulse photon echo peak shift (3PEPS) is a third-order nonlinear spectroscopy technique that records the time-correlation function of a solute molecule in a solution, including an offset (inhomogeneity). In this work, we utilized a diffractive optics-based 3PEPS apparatus to fully resolve the dynamics in liquids from femtoseconds to nanoseconds while varying the temperature in the range of 80–298 K and the probe solute molecules. Our observations reveal dynamics slower than the dielectric relaxation of n-alcohols, even at room temperature, consisting of a ∼0.5 ns time constant that persists below the melting points and a static component (offset) on a nanosecond timescale. Based on the experiments, we suggest that locally formed glass-like clusters in liquids can be responsible for the slow dynamics. Our results may provide new insights into the dynamics of liquids and related phenomena such as liquid-glass and liquid-liquid phase transitions. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In this brief review, we consider various characterizations of "monomeric" reversed phases for elucidating the interactions governing adsorbate retention in liquid chromatography. Conventional methods related to the assessment of retention capacity and hydrophobicity (specifically methylene selectivity) using single mobile phase compositions are discussed with a focus on dispersion interactions, along with their inherent strengths and limitations. An alternative approach involving separation maps through relative retention analysis is proposed. It is noted that, in real reversed-phase adsorbents, the density of the attached alkyl chains is typically one half of that of solid n-alkanes. In this case, adsorbate molecules to penetrate into the attached phase, and the process depends on the molecular shape. Consequently, conventional "monomeric" reversed phases exhibit specific selectivity towards substances with specific structures. The review also notes that current analytical methods often do not pay sufficient attention to the difference between the substance retention mechanisms, absorption and adsorption, because the predominant parameters of these mechanisms are quite different. Moreover, in the two most widely used very interesting and informative methods, linear solvation energy relationships (LSERs) and the hydrophobic-subtraction model, this characteristic has not received due attention. Taking into account that the method does not distinguish adsorbates retained by different mechanisms, absorptive versus adsorptive, to the obtained significant discrepancies between the calculated and experimental data do not seem extraordinary. The interpretation of the results of an LSER analysis is also complicated by uncertainties in the contributions of partial properties of adsorbates in both mobile and stationary phases to the total solvation energy, as only their difference is typically calculated. Nonetheless, a comparison of different columns in identical mobile phases can yield informative insights. A drawback of the second approach is the necessity of using multiple columns with substantial qualitative differences in the adsorbate retention among them. Furthermore, a possibility of the decomposition of all interactions into distinct types seems questionable, because the method does not involve any orthogonal (independent of the applied calculation method) properties. [ABSTRACT FROM AUTHOR]

[Display omitted] Hypothesis : Understanding the rules that control the assembly of nanostructured soft materials at interfaces is central to many applications. We hypothesize that electrolytes can be used to alter the hydration shell of amphiphilic oligomers at the air-aqueous interface of Langmuir films, thereby providing a means to control the formation of emergent nanostructures. Experiments : Three representative salts – (NaF, NaCl, NaSCN) were studied for mediating the self-assembly of oligodimethylsiloxane methylimidazolium (ODMS-MIM+) amphiphiles in Langmuir films. The effects of the different salts on the nanostructure assembly of these films were probed using vibrational sum frequency generation (SFG) spectroscopy and Langmuir trough techniques. Experimental data were supported by atomistic molecular dynamic simulations. Findings : Langmuir trough surface pressure – area isotherms suggested a surprising effect on oligomer assembly, whereby the presence of anions affects the stability of the interfacial layer irrespective of their surface propensities. In contrast, SFG results implied a strong anion effect that parallels the surface activity of anions. These seemingly contradictory trends are explained by anion driven tail dehydration resulting in increasingly heterogeneous systems with entangled ODMS tails and appreciable anion penetration into the complex interfacial layer comprised of headgroups, tails, and interfacial water molecules. These findings provide physical and chemical insight for tuning a wide range of interfacial assemblies. [ABSTRACT FROM AUTHOR]

Compared to conventional gels, eutectogels encompassing deep eutectic solvents (DESs) circumvent the shortboards including poor temperature resistance, high cost, and undesired toxicity. However, existing eutectogels suffer weak intra‐ and interchain interactions, which culminates in vanishing autofluorescence, weak mechanic and feeble adhesion. Although myriad strategies have been proposed to boost the polymer/polymer interactions within gels, the processes are time‐consuming, which are not suitable for mass production or emergent deployment. Herein, to deal with it, a facile solvation chemistry is leveraged to build dynamic rigid network via incorporating active polymer backbones and DESs. Within these eutectogels, DESs play a magnet‐like role, and physically cluster polymer chains together via multiple noncovalent forces, thus intensifying polymer/polymer and DESs/polymer interactions. This solvation network bestows rare autofluorescence and unprecedentedly high strength, toughness, and adhesion among reported analogues, even in harsh conditions. Moreover, the mechanism of solvation chemistry is fully scrutinized for the first time. It is believed that this facile strategy can be easily expanded to the orchestrations of other soft materials. [ABSTRACT FROM AUTHOR]

Aqueous Zn–metal batteries have attracted increasing interest for large-scale energy storage owing to their outstanding merits in terms of safety, cost and production. However, they constantly suffer from inadequate energy density and poor cycling stability due to the presence of zinc ions in the fully hydrated solvation state. Thus, designing the dehydrated solvation structure of zinc ions can effectively address the current drawbacks of aqueous Zn–metal batteries. In this case, considering the lack of studies focused on strategies for the dehydration of zinc ions, herein, we present a systematic and comprehensive review to deepen the understanding of zinc-ion solvation regulation. Two fundamental design principles of component regulation and pre-desolvation are summarized in terms of solvation environment formation and interfacial desolvation behavior. Subsequently, specific strategy based distinct principles are carefully discussed, including preparation methods, working mechanisms, analysis approaches and performance improvements. Finally, we present a general summary of the issues addressed using zinc-ion dehydration strategies, and four critical aspects to promote zinc-ion solvation regulation are presented as an outlook, involving updating (de)solvation theories, revealing interfacial evolution, enhancing analysis techniques and developing functional materials. We believe that this review will not only stimulate more creativity in optimizing aqueous electrolytes but also provide valuable insights into designing other battery systems. [ABSTRACT FROM AUTHOR]

The potential energy surfaces for the reactions of 1,3‐butadiene with 2‐hydroxythioacrolein and 2‐aminoacrolein exhibit ambimodal transition states leading to both dipolar (4 + 3) and Diels–Alder (4 + 2) cycloaddition products, thereby demonstrating a post transition state bifurcation feature. We have investigated the bifurcation dynamics of these reactions using three molecular dynamics (MD) methods: quasi‐classical trajectory, classical MD, and ring‐polymer MD simulations. The trajectory calculations were performed with the semiempirical GFN2‐xTB method with the element‐specific parameters optimized to reproduce the density‐functional theory calculations. The effect of water solvation was examined using an implicit solvation model, revealing significant differences in bifurcation dynamic between gas‐phase and solution‐phase reactions. Nuclear quantum effects were found to play a crucial role in the proton‐transfer process from the (4 + 3) intermediate to the (4 + 3) product in the case of the 2‐aminoacrolein reaction. [ABSTRACT FROM AUTHOR]

The utilization of deep learning methodologies for addressing partial differential equations (PDEs) has garnered significant attention in recent years. This paper introduces an improved network structure tailored for the discontinuity-capturing, enabling the resolution of interface problem through a unified neural network framework. Employing the probability space filling argument, we show that our model can generate convergent sequences, where the convergence rate depends on the number of sampling points. Several numerical experiments with regular and irregular interfaces are conducted to elucidate the convergence characteristics, thereby validating the theoretical assertions. Furthermore, we apply our approach to effectively solve the size-modified Poisson-Boltzmann test model, utilizing it for predicting electrostatics and the solvation free energies for proteins immersed in ionic solvents, thus showcasing practical applications of our method. • Improved multi-scale fusion network uses only a single neural network. • Our network model is proposed to solve the elliptical interface problem. • Theoretically give a least convergence rate of our numerical method. • Numerical experiments agree with our theoretical convergence rate. • We calculate the solvation free energies of immersed biomolecules by our method. [ABSTRACT FROM AUTHOR]

Exploiting the van der Waals model of liquids, it is possible to derive analytical formulas for the thermodynamic functions governing solvation, the transfer of a solute molecule from a fixed position in the ideal gas phase to a fixed position in the liquid phase. The solvation Gibbs free energy change consists of two contributions: (a) the high number density of all liquids and the repulsive interactions due to the basic fact that each molecule has its own body leading to the need to spend free energy to produce an appropriate cavity to contain the solute molecule; (b) the ubiquitous intermolecular attractive interactions lead to a gain in free energy for switching-on attractions between the solute molecule and neighboring liquid molecules. Also the solvation entropy change consists of two contributions: (a) there is an entropy loss in all liquids because the cavity presence limits the space accessible to liquid molecules during their continuous translations; (b) there is an entropy gain in all liquids, at room temperature, due to the liquid structural reorganization as a response to the perturbation represented by solute addition. The latter entropy contribution is balanced by a corresponding enthalpy term. The scenario that emerged from the van der Waals model is in qualitative agreement with experimental results. [ABSTRACT FROM AUTHOR]

Practical high‐voltage lithium metal batteries hold promise for high energy density applications, but face stability challenges in electrolytes for both 4 V‐class cathodes and lithium anode. To address this, we delve into the positive impacts of two crucial moieties in electrolyte chemistry: fluorine atom (‐F) and cyano group (‐CN) on the electrochemical performance of polyether electrolytes and lithium metal batteries. Cyano‐bearing polyether electrolytes possess strong solvation, accelerating Li+ desolvation with minimal SEI impact. Fluorinated polyether electrolytes possess weak solvation, and stabilize the lithium anode via preferential decomposition of F‐segment, exhibiting nearly 6000‐h stable cycling of lithium symmetric cell. Furthermore, the electron‐withdrawing properties of ‐F and ‐CN groups significantly bolster the high‐voltage tolerance of copolymer electrolyte, extending its operational range up to 5 V. This advancement enables the development of 4 V‐class lithium metal batteries compatible with various cathodes, including 4.45 V LiCoO2, 4.5 V LiNi0.8Co0.1Mn0.1O2, and 4.2 V LiNi0.5Co0.2Mn0.3O2. These findings provide insights into design principles centered around polymer components for high‐performance polymer electrolytes. [ABSTRACT FROM AUTHOR]

As an emerging secondary battery system, aqueous zinc-ion batteries (AZIBs) show a broad application prospect in the fields of large-scale energy storage and wearable devices. Manganese-based cathode materials have been widely investigated by many researchers due to their high natural abundance, low toxicity, and multiple variable valence states. However, limited active sites, insufficient solvation, and reactivity kinetics of Mn2+ lead to the attenuation of their electrochemical performance. Herein, we introduce appropriate oxygen vacancies into the δ-MnO2 structure by modulating the annealing temperature. The obtained δ-MnO2-400 electrode provided 503 mAh/g capacity at 0.2 A/g and 99% capacity retention after 3000 times cycling at 1 A/g. [ABSTRACT FROM AUTHOR]

Rational design of advanced electrolytes to improve the high‐voltage capability has been attracting wide attention as one critical solution to enable next‐generation high‐energy‐density batteries. However, the limited understanding of electrolyte antioxidant chemistry as well as the lack of valid quantization approaches have resulted in knowledge gap, which hinders the formulation of new electrolytes. Herein, we construct a standard curve based on representative solvation structures to quantify the oxidation stability of ether‐based electrolytes, which reveals the linear correlation between the oxidation potential and the atomic charge of the least oxidation‐resistant solvent. Dictated by the regularity between solvation composition and oxidation potential, a (Trifluoromethyl)cyclohexane‐based localized high‐concentration electrolyte dominated by anion‐less solvation structures was designed to optimize the cycling performance of 4.5 V 30 μm‐Li||3.8 mAh cm−2‐LiCoO2 batteries, which maintained 80 % capacity retention even after 440 cycles. The consistency of experimental and computational results validates the proposed principles, offering a fundamental guideline to evaluate and design aggressive electrochemical systems. [ABSTRACT FROM AUTHOR]

Constructing powerful anode/cathode interphases by modulate ion solvation structure is the principle of electrolyte design. However, the methodological and theoretical design principles of electrolyte/solvation structure and their effect on electrochemical performance are still vague. Here, we propose a cationic weakly coordinating-intervention strategy for modulating the Na+ solvation sheathes and constructing robust anode/cathode interphases in sodium-metal batteries. Unlike the local highly concentrated electrolytes, 1,2-difluorobenzene can weakly coordinate with Na+ thus transforming the solvation structure into Na+-anion-incorporated structures and strengthening anode/cathode interphases formation by combining with salt decomposition. Furthermore, the correlations between the electrode interface properties and solvation structure are revealed, which can be tuned by the weakly coordination. Ultimately, the modulated electrolyte achieves 97.5% Coulombic efficiency for 600 cycles in Na‖Cu cells at 1 mA cm−2 and a beneficial lifetime (2500 h) in Na‖Na cells. Meanwhile, Na‖PB cells have achieved long-term operation at 4.8 V, along with operation at wide temperatures. Here, authors propose a weakly coordinating-intervention strategy to modulate the Na+ solvation sheath and construct a robust interphase in sodium-metal batteries. Correlation between electrode interface properties and solvation structure is also stated. [ABSTRACT FROM AUTHOR]

Electrolytes with anion‐dominated solvation are promising candidates to achieve dendrite‐free and high‐voltage potassium metal batteries. However, it's challenging to form anion‐reinforced solvates at low salt concentrations. Herein, we construct an anion‐reinforced solvation structure at a moderate concentration of 1.5 M with weakly coordinated cosolvent ethylene glycol dibutyl ether. The unique solvation structure accelerates the desolvation of K+, strengthens the oxidative stability to 4.94 V and facilitates the formation of inorganic‐rich and stable electrode‐electrolyte interface. These enable stable plating/stripping of K metal anode over 2200 h, high capacity retention of 83.0 % after 150 cycles with a high cut‐off voltage of 4.5 V in K0.67MnO2//K cells, and even 91.5 % after 30 cycles under 4.7 V. This work provides insight into weakly coordinated cosolvent and opens new avenues for designing ether‐based high‐voltage electrolytes. [ABSTRACT FROM AUTHOR]

This study investigates the application of quantum mechanical (QM) and multiscale computational methods in understanding the reaction mechanisms and kinetics of SN2 reactions involving methyl iodide with NH2OH and NH2O−, as well as the Claisen rearrangement of 8-(vinyloxy)dec-9-enoate. Our aim is to evaluate the accuracy and effectiveness of these methods in predicting experimental outcomes for these organic reactions. We achieve this by employing QM-only calculations and several hybrids of QM and molecular mechanics (MM) methods, namely QM/MM, QM1/QM2, and QM1/QM2/MM methodologies. For the SN2 reactions, our results demonstrate the importance of explicitly including solvent effects in the calculations to accurately reproduce the transition state geometry and energetics. The multiscale methods, particularly QM/MM and QM1/QM2, show promising performance in predicting activation energies. Moreover, we observe that the size of the MM active region significantly affects the accuracy of calculated activation energies, highlighting the need for careful consideration during the setup of multiscale calculations. In the case of the Claisen rearrangement, both QM-only and multiscale methods successfully reproduce the proposed reaction mechanism. However, the activation free energies calculated using a continuum solvation model, based on single-point calculations of QM-only structures, fail to account for solvent effects. On the other hand, multiscale methods more accurately capture the impact of solvents on activation free energies, with systematic error correction enhancing the accuracy of the results. Furthermore, we introduce a Python code for setting up multiscale calculations with ORCA, which is available on GitHub at https://github.com/iranimehdi/pdbtoORCA. [ABSTRACT FROM AUTHOR]

Although interfacial solvation plays an important role in determining carbon dioxide reduction (CO2R) kinetics, present understanding of the potential dependent properties of the electrochemical double layer under conditions relevant for CO2R remains limited. This article summarizes the development and recent applications of plasmon‐enhanced vibrational sum frequency generation (VSFG) spectroscopy to study the effects of cation hydration and interfacial solvation on CO2R using CO as a vibrational Stark reporter. Results show that electrolyte cations retain their entire solvation shell upon adsorption to inactive sites, while active sites retain only a single water layer between the gold surface and the cation. Measurements also show that the total interfacial electric field can be separated into two contributions: one from the electrochemical double layer (Stern field) and another from the polar solvation environment (Onsager field). Surprisingly, correlating VSFG spectra with reaction kinetics reveals that it is the solvation‐mediated Onsager field that governs the chemical reactivity at the electrode/electrolyte interface. Measuring the interfacial water spectra during electrocatalysis also provides evidence for the proton source during H2 evolution, which competes with CO2R in aqueous electrolyte. These findings highlight the importance of directly probing cation hydration and interfacial solvation, which mediates reaction kinetics at electrochemical interfaces. [ABSTRACT FROM AUTHOR]

Aqueous zinc‐ion batteries have garnered renewed interest owing to their inherent safety and cost‐effectiveness. However, Zn anode suffers from notorious side reactions and dendrite growth, compromising the battery performance. Designing anion‐containing solvation structures to facilitate the in situ formation of a solid electrolyte interphase (SEI) layer proves effective in protecting Zn anode, yet balancing the inherent interactions between the cations‐solvents and cations–anions remains a significant challenge. Herein, latent solvent (hexafluoroisopropanol) with specific anion‐philic and cation‐phobic properties is introduced, which can remarkably increase Zn(TFSI)2 solubility in water and induce anion‐containing solvation structures. The introduction of latent solvent enables anions to preferentially decompose, giving rise to an anion‐derived SEI layer. More interestingly, the direct anion‐cation interaction endows Zn2+ species with stepwise dehydration, minimizing the adsorption of water on Zn anode. The anion‐derived SEI layer and optimized de‐solvation process work in synergy to enable homogeneous Zn deposition and ensure high reversibility. When coupled with a NaV3O8·1.5H2O (NVO) cathode, the resultant full cell delivers a high‐capacity retention of 87.8% after 2000 cycles at 1 A g−1. This work provides valuable insights into the advanced latent solvent electrolyte design and beyond. [ABSTRACT FROM AUTHOR]

Octa-acid 1 is a water-soluble cavitand that has been used to investigate hydrophobic solvation and Hofmeister effects, control photophysical and physicochemical properties, modulate the reactivity of encapsulated guests, and as a tool to engender novel separation protocols. The synthesis of 1 has largely centered around its formation from 2, a host that is itself most readily synthesized on the multi-gram-scale in crude form (>75% purity). In this Methods Article we reveal improvements in the synthesis of 2, as well as a new synthetic strategy that efficiently converts crude 2 into pure 1. This provides access to 1 in only a five-step linear sequence, shortening the total reaction time from previous methods, improving the purity, and increasing the final yield of 1. We therefore anticipate that the described protocols will be of interest to researchers seeking to utilize octa-acid in their studies. [ABSTRACT FROM AUTHOR]

We study, using the Molecular Dynamics method, the mobility of Lithium compensation cations in two clay models that are in contact with a water reservoir. The preferential site of Li+ is the center of the hexagonal ring of the structure. The simulations show that Li+ cations can leave their initial positions and move in four different ways. These movements occur over short time intervals, suggesting a jump diffusion mechanism. Furthermore, our simulations have highlighted that Li+ can coordinate with the framework oxygens and with water molecules, and the solvation will depend on the position of the cation. Likewise, the cations close to the surface are the ones that diffuse the most. [ABSTRACT FROM AUTHOR]

The formation of ABC triblock terpolymers through solution casting is still challenging. In this study, core–shell double gyroid network structures are fabricated via solution casting using poly(2,2,2‐trifluoroethyl methacrylate) (PTFEMA) (F)‐b‐[poly(4‐vinylpyridine) (P4VP) (P)]‐b‐[polystyrene (PS) (S)] (FPS) triblock terpolymers in N,N‐dimethylformamide (DMF). Upon heat treatment, the polymer tends to form a sphere‐in‐lamellar structure at the F/S interface. Given the solubility properties of each component in DMF, it is anticipated that the effective volume fraction of F relative to P would increase in concentrated solutions and the effective volume fraction of S would decrease. The microphase‐separated structure derived from the DMF solution consistently results in the formation of a network structure composed of a core–shell double gyroid, with F as the matrix, P as the shell, and S as the core, and their periodic lengths gradually increase to 110.8, 131.8, and 162.7 nm as increase molecular weights of PS blocks to 13.8, 20.7, and 28.8 kg mol−1. Based on the solubility properties of the polymer components highlighted in this study, the solvent selection strategy is broadly applicable to ABC triblock terpolymers featuring various polymer components, offering a more efficient avenue for fabricating core–shell double gyroid structures. [ABSTRACT FROM AUTHOR]

Current potassium‐ion batteries (PIBs) are limited in safety and lifetime owing to the lack of suitable electrolyte solutions. To address these issues, herein, we report an innovative non‐flammable electrolyte design strategy that leverages an optimal moderate solvation phosphate‐based solvent which strikes a balance between solvation capability and salt dissociation ability, leading to superior electrochemical performance. The formulated electrolyte simultaneously exhibits the advantages of low salt concentration (only 0.6 M), low viscosity, high ionic conductivity, high oxidative stability, and safety. Our electrolyte also promotes the formation of self‐limiting inorganic‐rich interphases at the anode surface, alongside robust cathode‐electrolyte interphase on iron‐based Prussian blue analogues, mitigating electrode/electrolyte side reactions and preventing Fe dissolution. Notably, the PIBs employing our electrolyte exhibit exceptional durability, with 80 % capacity retention after 2,000 cycles at high‐voltage of 4.2 V in a coin cell. Impressively, in a larger scale pouch cell, it maintains over 81 % of its initial capacity after 1,400 cycles at 1 C‐rate with high average Coulombic efficiency of 99.6 %. This work represents a significant advancement toward the realization of safe, sustainable, and high‐performance PIBs. [ABSTRACT FROM AUTHOR]

Structural changes induced by water play a pivotal role in chemistry and biology but remain challenging to predict, measure, and control at molecular level. Here we explore size‐governed gas‐phase water aggregation in the flexible molecule 4‐hydroxy‐2‐butanone, modeling the conformational adaptability of flexible substrates to host water scaffolds and the preference for sequential droplet growth. The experiment was conducted using broadband rotational spectroscopy, rationalized with quantum chemical calculations. Two different isomers were observed experimentally from the di‐ to the pentahydrates (4‐hydroxy‐2‐butanone‐(H2O)n=2–5), including the 18O isotopologues for the di‐ and trihydrates. Interestingly, to accommodate water molecules effectively, the heavy atom skeleton of 4‐hydroxy‐2‐butanone reshapes in every observed isomer and does not correspond to the stable conformer of the free monomer. All solvates initiate from the alcohol group (proton donor) but retain the carbonyl group as secondary binding point. The water scaffolds closely resemble those found in the pure water clusters, balancing between the capability of 4‐hydroxy‐2‐butanone for steering the orientation and position of the water molecules and the ability of water to modulate the monomer's conformation. The present work thus provides an accurate molecular description on how torsionally flexible molecules dynamically adapt to water along progressing solvation. [ABSTRACT FROM AUTHOR]

This paper focuses on individual responses to education in crisis, with the strategies of students contextualised and examined within a wider multi‐crisis reality. Offering a conceptual framework of social solvation, the proposed model explains how failures in the education system at the macro (state) and meso (institutional) levels translate into young people's behaviours at the individual (micro) level. Based on the analysis of 26 in‐depth interviews with young students, we identify and discuss three main micro‐rational strategies that the interviewees relied upon to navigate their learning during the pandemic: (1) casualisation (2) customisation and (3) contempt towards the broader system. We argue that the already weak Polish education system was further impaired by the COVID‐19 crisis. In response to growing social and educational risks, young people forge and employ highly individualised strategies. As micro‐rational social actions essentially bypass the unwieldy education system in a variety of ways, they have vital implications for educational futures. [ABSTRACT FROM AUTHOR]

Layered oxides of sodium‐ion batteries suffer from severe side reactions on the electrode/electrolyte interface, leading to fast capacity degradation. Although surface reconstruction strategies are widely used to solve the above issues, the utilization of the low‐cost wet chemical method is extremely challenging for moisture‐sensitive Na‐based oxide materials. Here, the solvation tuning strategy is proposed to overcome the deterioration of NaNi1/3Mn1/3Fe1/3O2 in water‐based solution and conduct the surface reconstruction. When capturing the water molecules by the solvation structure of cations, here is Li+, the structural collapse and degradation of layered oxides in water‐based solvents are greatly mitigated. Furthermore, Li(H2O)3EA+ promotes the profitable Li+/Na+ exchange to build a robust surface, which hampers the decomposition of electrolytes and the structural evolution upon cycling. Accordingly, the lifespan of Li‐reinforced materials is prolonged to three times that of the pristine one. This work represents a step forward in understanding the surface reconstruction operated in a water‐based solution for high‐performance sodium layered oxide cathodes. [ABSTRACT FROM AUTHOR]

Lithium metal is the most promising high‐energy‐density anode. However, it is incompatible with high‐voltage cathodes in ether solvents due to their narrow electrochemical window. Herein, fluoroethylene carbonate (FEC) co‐solvent is introduced to regulate the Li+ solvation structures in ether solvents, including cyclic ether (1,3‐dioxolane [DOL]) and linear glymes with different chain lengths (1,2‐dimethoxyethane [DME], diglyme [G2] and triglyme [G3]). The apparently different effects of ether solvents on solvation ability and interaction strength with FEC are revealed. FEC plays a diverse role and function in 1 m lithium bis(trifluoromethanesulfonyl)imide (LiTFSI)‐ether/FEC electrolyte, thus relevant batteries perform distinct performances due to various ionic dynamics and solid‐electrolyte interphase. The Li+‐solvation structures are explored by Raman and nuclear magnetic resonance spectroscopies. Specifically, part of FEC molecules are inserted into the first solvation shell in 1 m LiTFSI‐DOL/FEC because of the weak solvation ability of DOL and strong interaction of DOL‐FEC, leading to few coordinated TFSI− and sluggish interfacial kinetics. In sharp contrast, FEC as a weak coordinated solvent almost exclusively occupies the second solvation sheath in 1 m LiTFSI‐glyme/FEC, favoring TFSI− coordination and rapid de‐solvation dynamics. Ultimately, the LiNi0.8Co0.1Mn0.1O2/Li battery in G2/FEC presents the most excellent performance, derived from abundant free‐FEC and rapid ionic kinetics. [ABSTRACT FROM AUTHOR]