Many industrial processes, several natural processes involving non-living matter, and all the processes occurring within living organisms take place in solution. This means that the molecules playing active roles in the processes are present within another medium, called solvent. The solute molecules are surrounded by solvent molecules and interact with them. Understanding the nature and strength of these interactions, and the way in which they modify the properties of the solute molecules, is important for a better understanding of the chemical processes occurring in solution, including possible roles of the solvent in those processes. Computational studies can provide a wealth of information on solute-solvent interactions and their effects. Two major models have been developed to this purpose: a model viewing the solvent as a polarizable continuum surrounding the solute molecule, and a model considering a certain number of explicit solvent molecules around a solute molecule. Each of them has its advantages and challenges, and one selects the model that is more suitable for the type of information desired for the specific system under consideration. These studies are important in many areas of chemistry research, from the investigation of the processes occurring within a living organism to drug design and to the design of environmentally benign solvents meant to replace less benign ones in the chemical industry, as envisaged by the green chemistry principles. The paper presents a quick overview of the modelling approaches and an overview of concrete studies, with reference to selected crucial investigation themes.

Meta-aminobenzoic acid, an important model system in the study of polymorphism and crystallization of active pharmaceutical ingredients, exist in water in both the nonionic (mABA) and zwitterionic (mABA^±) forms. However, the constituent molecules of the polymorph that crystallizes from aqueous solutions are zwitterionic. This study reports atomistic simulations of the events surrounding the early stage of crystal nucleation of meta-aminobenzoic acid from aqueous solutions. Ab initio molecular dynamics was used to simulate the hydration of mABA± and mABA, and to quantify the interaction of these molecules with the surrounding water molecules. Density functional theory calculations were conducted to determine the low-lying energy conformers of meta-aminobenzoic acid dimers and compute the Gibbs free energies in water of nonionic, (mABA)₂, zwitterionic, (mABA^±)₂ and nonionic-zwitterionic, (mABA)(mABA^±), species. Classical molecular dynamics simulations of mixed mABA–mABA^± aqueous solutions were carried out to examine the aggregation of meta-aminobenzoic acid. According to these simulations the selective crystallization of the polymorph which constituent molecules are zwitterionic is driven by the formation of zwitterionic dimers in solution, which are thermodynamically more stable than (mABA)₂ and (mABA)(mABA^±) pairs. This work represents a paradigm of the role of molecular processes during the early stages of crystal nucleation in affecting polymorph selection during crystallization from solution.

Asymmetric lithium salts, such as Lithium (difluoromethanesulfonyl)(trifluoromethanesulfonyl)imide (LiDFTFSI), have been demonstrated to surpass traditional symmetric lithium salts with improved Li+ conductivity and the capacity to generate stable solid electrolyte interphase (SEI) while maintaining compatibility with aluminum (Al0) current collector. However, the intrinsic reductive mechanism through which LiDFTFSI influences battery performance remains unclear and under debate. Herein, the detailed SEI reactions of LiDFTFSI-based electrolytes were investigated by combining density functional theory and molecular dynamics, aiming to clarify the formation process and atomic structure of SEI. Our results show that asymmetric DFTFSI– weakens the interaction between carbonate solvents and Li+, and substantially alter the solvation structure, exhibiting a well-balanced coordination capacity to bis(trifluoromethanesulfonyl)imide (TFSI–). Nanoseconds hybrid molecular dynamics simulation further reveals that the preferential decomposition of LiDFTFSI produces sufficient LiF and Li2O to facilitate a robust SEI. Moreover, the abundant F– generated from LiDFTFSI decomposition accumulates on the Al surface and subsequently combines with Al3+ from the current collector to form AlF3, potentially inhibiting corrosion of the current collector. Overall, these findings elucidate how LiDFTFSI regulates the solvation sheath and SEI structure, advancing the development of high-performance electrolytes compatible with current collector.

For condensed systems, the incorporation of quantum chemical solvent effects into molecular dynamics simulations has been a major concern. To this end, quantum mechanical/molecular mechanical (QM/MM) techniques are popular and powerful options to treat gigantic systems. However, they cannot be directly applied because of temporal and spatial discontinuity problems. To overcome these problems, in a previous study, we proposed a corrective QM/MM method, size-consistent multipartitioning (SCMP) QM/MM, and successfully demonstrated that, using SCMP, it is possible to perform stable molecular dynamics simulations by effectively taking into account solvent quantum chemical effects. The SCMP method is characterized by two original features: size-consistency of a QM region among all QM/MM partitioning and partitioning update. However, in our previous study, the performance was not fully elicited compared to the theoretical upper bound, and the optimal partitioning update protocol and parameters were not fully verified. To elicit the potential performance, in the present study, we simplified the theoretical framework and modified the partitioning protocol.

The Lewis basicity of selected organic bases, modeled by the enthalpies of adduct formation between gaseous BF3 and the bases in dichloromethane (DCM) solution, is critically examined. Although experimental enthalpies for a large number of molecules have been reported in the literature, it may be desirable to estimate missing or uncertain data for important Lewis bases. We have decided to use high-level ab initio procedures, combined with a polarized continuum solvation model, in which the solvated species are the clusters formed by specific hydrogen bonding of DCM with the Lewis base and the Lewis base/BF3 adduct. This mode of interaction with DCM corresponds to a specific solvation model (SSM). The results actually show that the enthalpy of BF3 adduct formation in DCM solution is clearly influenced by specific interactions, DCM acting as hydrogen-bonding donor (HBD) molecule in two ways: base/DCM and adduct/DCM, confirming that specific solvation is an important contribution to experimentally determined Lewis basicity scales. This analysis allows us to conclude that there are reasons to suspect some gas-phase values to be in error by more than the stated experimental uncertainty. Some experimental values in DCM solution that were uncertain because of identified reasons can be complemented by the computed values.

Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, Most of important methodological advancements in more than half century of molecule modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science based on force fields parameterization by coarse graining, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but no transferability is available. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science and a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.

Reactive oxygen species (ROS) play a pivotal role in cellular function, and their overproduction contributes to oxidative stress, fatigue, aging, and diseases. To counteract these effects, incorpo-rating antioxidant-rich nutrients is recommended. Functional waters, including electrolyzed hydrogen water (EHW), are gaining attention for potential antioxidant effects through daily in-take. This study employs aquaphotomics and near-infrared spectroscopy (NIRS) to elucidate the molecular mechanisms behind the antioxidant activity of EHW. Previous investigation reported a 100% loss in antioxidant functionality through autoclaving, and a 40% decrease with hydrogen gas removal, prompting the current investigation into the mechanisms of antioxidant effects and especially the role of water molecular structure in it. Aquaphotomics analysis of near infrared (NIR) spectral data revealed that the primary characteristic of EHW associated with anti-radical function is presence of water solvation shells. Autoclaving and degassing led to loss of these water molecular species and a significant increase in free water molecules. The study suggests that water shells (solvation shells) in EHW play a crucial role in neutralizing ROS, although the specific ions or nanoparticles responsible for creating these shells remain unidentified. This concept aligns well with existing proposals of the crucial role of water solvation shells in ROS neutralization.

The calculation of the standard enthalpies of vaporization, sublimation and solvation of organic molecules is presented using a common computer algorithm on the basis of a group-additivity method. The same algorithm is also shown to enable the calculation of their entropy of fusion as well as the total phase-change entropy of liquid crystals. The present method is based on the complete break-down of the molecules into their constituting atoms and their immediate neighbourhood; the respective calculations of the contribution of the atomic groups by means of the Gauss-Seidel fitting method is based on experimental data collected from literature. The feasibility of the calculations for each of the mentioned descriptors was verified by means of a 10-fold cross-validation procedure proving the good to high quality of the predicted values for the three mentioned enthalpies and for the entropy of fusion, whereas the predictive quality for the total phase-change entropy of liquid crystals was poor. The goodness of fit (Q²) and the standard deviation (σ) of the cross-validation calculations for the five descriptors was as follows: 0.9641 and 4.56 kJ/mol (N=3386 test molecules) for the enthalpy of vaporization, 0.8657 and 11.39 kJ/mol (N=1791) for the enthalpy of sublimation, 0.9546 and 4.34 kJ/mol (N=373) for the enthalpy of solvation, 0.8727 and 17.93 J/mol/K (N=2637) for the entropy of fusion and 0.5804 and 32.79 J/mol/K (N=2643) for the total phase-change entropy of liquid crystals. The large discrepancy between the results of the two closely related entropies is discussed in detail. Molecules, for which both the standard enthalpies of vaporization and sublimation were calculable, enabled the estimation of their standard enthalpy of fusion by simple subtraction of the former from the latter enthalpy. For 990 of them the experimental enthalpy-of-fusion values are also known, allowing their comparison with predictions, yielding a correlation coefficient R² of 0.6066.