SM6: A Density Functional Theory Continuum Solvation Model for Calculating Aqueous Solvation Free Energies of Neutrals, Ions, and Solute-Water Clusters

Quantum chemical package Jaguar: A survey of recent developments and unique features

This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.

PW-SMD: A Plane-Wave Implicit Solvation Model Based on Electron Density for Surface Chemistry and Crystalline Systems in Aqueous Solution

Electron density-based implicit solvation models are a class of techniques for quantifying solvation effects and calculating free energies of solvation without an explicit representation of solvent molecules. Integral to the accuracy of solvation modeling is the proper definition of the solvation shell separating the solute molecule from the solvent environment, allowing for a physical partitioning of the free energies of solvation. Unlike state-of-the-art implicit solvation models for molecular quantum chemistry calculations, e.g., the solvation model based on solute electron density (SMD), solvation models for systems under periodic boundary conditions with plane-wave (PW) basis sets have been limited in their accuracy. Furthermore, a unified implicit solvation model with both homogeneous solution-phase and heterogeneous interfacial structures treated on equal footing is needed. In order to address this challenge, we developed a high-accuracy solvation model for periodic PW calculations that is applicable to molecular, ionic, interfacial, and bulk-phase chemistry. Our model, PW-SMD, is an extension of the SMD molecular solvation model to periodic systems in water. The free energy of solvation is partitioned into the electrostatic and cavity-dispersion-solvent structure (CDS) contributions. The electrostatic contributions of the solvation shell surrounding solute structures are parametrized based on their geometric and physical properties. In addition, the nonelectrostatic contribution to the solvation energy is accounted for by extending the CDS formalism of SMD to incorporate periodic boundary conditions. We validate the accuracy and robustness of our solvation model by comparing predicted solvation free energies against experimental data for molecular and ionic systems, carved-cluster composite energetic models of solvated reaction energies and barriers on surface systems, and deep-learning-accelerated ab initio molecular dynamics (AIMD). Our developed periodic implicit solvation model shows significantly improved accuracy compared to previous work (namely, solvation models in aqueous solution) and can be applied to simulate solvent effects in a wide range of surface and crystalline materials.

The Chlorido-Bismuth Dication: A Potent Lewis Acid Captured in a Hepta-Coordinate Species with a Stereochemically Active Lone Pair

The stabilization of simple, highly reactive cationic species in molecular complexes represents an important strategy to isolate and characterize compounds with uncommon or even unprecedented structural motifs and properties. Here we report the synthesis, isolation, and full characterization of chlorido-bismuth dications, stabilized only by monodentate dimethylsulfoxide (dmso) ligands: [BiCl(dmso)6][BF4]2 (1) and [BiCl(μ2-dmso)(dmso)4]2[BF4]4 (2). These compounds show unusual distorted pentagonal bipyramidal coordination geometries along with high Lewis acidities and have been analyzed by multinuclear NMR spectroscopy, elemental analysis, IR spectroscopy, single-crystal X-ray diffraction, and density functional theory calculations. Attempts to generate the bromido- and iodido-analogs gave dmso-stabilized tricationic bismuth species.

Lambda-ABF: Simplified, Portable, Accurate, and Cost-Effective Alchemical Free-Energy Computation

We introduce the lambda-Adaptive Biasing Force (lambda-ABF) method for the computation of alchemical free-energy differences. We propose a software implementation and showcase it on biomolecular systems. The method arises from coupling multiple-walker adaptive biasing force with λ-dynamics. The sampling of the alchemical variable is continuous and converges toward a uniform distribution, making manual optimization of the λ schedule unnecessary. Contrary to most other approaches, alchemical free-energy estimates are obtained immediately without any postprocessing. Free diffusion of λ improves orthogonal relaxation compared to fixed-λ thermodynamic integration or free-energy perturbation. Furthermore, multiple walkers provide generic orthogonal space coverage with minimal user input and negligible computational overhead. We show that our high-performance implementations coupling the Colvars library with NAMD and Tinker-HP can address real-world cases including ligand-receptor binding with both fixed-charge and polarizable models, with a demonstrably richer sampling than fixed-λ methods. The implementation is fully open-source, publicly available, and readily usable by practitioners of current alchemical methods. Thanks to the portable Colvars library, lambda-ABF presents a unified user interface regardless of the back-end (NAMD, Tinker-HP, or any software to be interfaced in the future), sparing users the effort of learning multiple interfaces. Finally, the Colvars Dashboard extension of the visual molecular dynamics (VMD) software provides an interactive monitoring and diagnostic tool for lambda-ABF simulations.

Designing solvent systems using self-evolving solubility databases and graph neural networks

Designing solvent systems is key to achieving the facile synthesis and separation of desired products from chemical processes, so many machine learning models have been developed to predict solubilities. However, breakthroughs are needed to address deficiencies in the model's predictive accuracy and generalizability; this can be addressed by expanding and integrating experimental and computational solubility databases. To maximize predictive accuracy, these two databases should not be trained separately, and they should not be simply combined without reconciling the discrepancies from different magnitudes of errors and uncertainties. Here, we introduce self-evolving solubility databases and graph neural networks developed through semi-supervised self-training approaches. Solubilities from quantum-mechanical calculations are referred to during semi-supervised learning, but they are not directly added to the experimental database. Dataset augmentation is performed from 11 637 experimental solubilities to >900 000 data points in the integrated database, while correcting for the discrepancies between experiment and computation. Our model was successfully applied to study solvent selection in organic reactions and separation processes. The accuracy (mean absolute error around 0.2 kcal mol-1 for the test set) is quantitatively useful in exploring Linear Free Energy Relationships between reaction rates and solvation free energies for 11 organic reactions. Our model also accurately predicted the partition coefficients of lignin-derived monomers and drug-like molecules. While there is room for expanding solubility predictions to transition states, radicals, charged species, and organometallic complexes, this approach will be attractive to predictive chemistry areas where experimental, computational, and other heterogeneous data should be combined.

An implicit electrolyte model for plane wave density functional theory exhibiting nonlinear response and a nonlocal cavity definition

We have developed and implemented an implicit electrolyte model in the Vienna Ab initio Simulation Package (VASP) that includes nonlinear dielectric and ionic responses as well as a nonlocal definition of the cavities defining the spatial regions where these responses can occur. The implementation into the existing VASPsol code is numerically efficient and exhibits robust convergence, requiring computational effort only slightly higher than the original linear polarizable continuum model. The nonlinear + nonlocal model is able to reproduce the characteristic "double hump" shape observed experimentally for the differential capacitance of an electrified metal interface while preventing "leakage" of the electrolyte into regions of space too small to contain a single water molecule or solvated ion. The model also gives a reasonable prediction of molecular solvation free energies as well as the self-ionization free energy of water and the absolute electron chemical potential of the standard hydrogen electrode. All of this, combined with the additional ability to run constant potential density functional theory calculations, should enable the routine computation of activation barriers for electrocatalytic processes.

Experimental Compilation and Computation of Hydration Free Energies for Ionic Solutes

Although charged solutes are common in many chemical systems, traditional solvation models perform poorly in calculating solvation energies of ions. One major obstacle is the scarcity of experimental data for solvated ions. In this study, we release an experiment-based aqueous ionic solvation energy data set, IonSolv-Aq, that contains hydration free energies for 118 anions and 155 cations, more than 2 times larger than the set of hydration free energies for singly charged ions contained in the 2012 Minnesota Solvation Database commonly used in benchmarking studies. We discuss sources of systematic uncertainty in the data set and use the data to examine the accuracy of popular implicit solvation models COSMO-RS and SMD for predicting solvation free energies of singly charged ionic solutes in water. Our results indicate that most SMD and COSMO-RS modeling errors for ionic solutes are systematic and correctable with empirical parameters. We discuss two systematic offsets: one across all ions and one that depends on the functional group of the ionization site. After correcting for these offsets, solvation energies of singly charged ions are predicted using COSMO-RS to 3.1 kcal mol-1 MAE against a challenging test set and 1.7 kcal mol-1 MAE (about 3% relative error) with a filtered test set. The performance of SMD is similar, with MAE against those same test sets of 2.7 and 1.7 kcal mol-1. These results underscore the importance of compiling larger experimental data sets to improve solvation model parametrization and fairly assess performance.

Solvent-controlled formation of alkali and alkali-earth-secured cucurbituril/guest trimers

Cucurbit[7]uril (CB[7]) encapsulates adamantyl and trimethylsilyl substituents of positively charged guests in dimethyl sulfoxide (DMSO). Unlike in water or deuterium oxide, addition of a selection of alkali and alkali-earth cations with van der Waals radii between 1.0 and 1.4 Å (Na+, K+, Ca2+, Sr2+, Ba2+ and Eu3+) to the CB[7]/guest complexes triggers their cation-mediated trimerization, a process that is very slow on the nuclear magnetic resonance (NMR) time scale. Smaller (Li+, Mg2+) or larger cations (Rb+, Cs+ or NH4+) are inert. The trimers display extensive CH-O interactions between the equatorial and pseudo-equatorial hydrogens of CB[7] and the carbonyl rim of the neighboring CB[7] unit in the trimer, and a deeply nested cation between the three interacting carbonylated CB[7] rims; a counteranion is likely perched in the shallow cavity formed by the three outer walls of CB[7] in the trimer. Remarkably, a guest must occupy the cavity of CB[7] for trimerization to take place. Using a combination of semi-empirical and density functional theory techniques in conjunction with continuum solvation models, we showed that trimerization is favored in DMSO, and not in water, because the penalty for the partial desolvation of three of the six CB[7] portals upon aggregation into a trimer is less unfavorable in DMSO compared to water.

Theoretical Insights on the Sensing Performance for Newly-synthesized Two-photon Fluorescent N2H4 Probes Based on Spirobifluorence

The development of fluorescent probe for hydrazine (N2H4) detection has attracted much attention due to the important role of N2H4 plays in the fields of medicine, agriculture, biology and environments. In this paper, the optical properties and water solubility of two novel two-photon fluorescent molecular probes (Probe1 and Probe2) before and after the reaction with N2H4 are studied by using the density function theory. The results show that electronic distribution and transition dipole moment of the probes are obviously changed after the reaction with N2H4, thus the optical properties of the molecules are influenced and the detection of N2H4 are realized. In addition, photoinduced electron transfer processes for Probe1 and Probe2 in the presence of N2H4 are theoretically characterized, which explains the experimental observations from the microscopic mechanism. Special attention has been paid on the analysis of the two-photon absorption for the probes with the absence and presence of N2H4 by the response theory method. Both probes with good water solubility show large variation on the two-photon absorption cross section when reacts with N2H4. In particular, the two-photon absorption response of Probe2 is more obvious, so it possesses preferable two-photon fluorescence microscopic imaging ability. More importantly, the receptor effect on the sensing performances of the probes are demonstrated, providing a theoretical reference for the design and synthesis on more efficient two-photon fluorescence N2H4 probes. Our study provides necessary information on the response mechanism of the studied chemosensors and helps to establish the relationship between the structure and optical properties of two-photon fluorescence N2H4 probes.

A generalized Kirkwood implicit solvent for the polarizable AMOEBA protein model

Computational simulation of biomolecules can provide important insights into protein design, protein-ligand binding interactions, and ab initio biomolecular folding, among other applications. Accurate treatment of the solvent environment is essential in such applications, but the use of explicit solvents can add considerable cost. Implicit treatment of solvent effects using a dielectric continuum model is an attractive alternative to explicit solvation since it is able to describe solvation effects without the inclusion of solvent degrees of freedom. Previously, we described the development and parameterization of implicit solvent models for small molecules. Here, we extend the parameterization of the generalized Kirkwood (GK) implicit solvent model for use with biomolecules described by the AMOEBA force field via the addition of corrections to the calculation of effective radii that account for interstitial spaces that arise within biomolecules. These include element-specific pairwise descreening scale factors, a short-range neck contribution to describe the solvent-excluded space between pairs of nearby atoms, and finally tanh-based rescaling of the overall descreening integral. We then apply the AMOEBA/GK implicit solvent to a set of ten proteins and achieve an average coordinate root mean square deviation for the experimental structures of 2.0 Å across 500 ns simulations. Overall, the continued development of implicit solvent models will help facilitate the simulation of biomolecules on mechanistically relevant timescales.

Evolution of vibrational bands upon gradual protonation/deprotonation of arsinic acid H2As(O)OH in media of different polarity

This computational work is devoted to the investigation (MP2/def2-TZVP) of the geometry and IR parameters of arsinic acid H₂AsOOH and its hydrogen-bonded complexes under vacuum and in media with different polarity. The medium effects were accounted for in two ways: (1) implicitly, using the IEFPCM model, varying the dielectric permittivity (ε) and (2) explicitly, by considering hydrogen-bonded complexes of H₂As(O)OH with various hydrogen bond donors (41 complexes) or acceptors (38 complexes), imitating a gradual transition to the As(OH)₂⁺ or AsO₂^- moiety, respectively. It was shown that the transition from vacuum to a medium with ε > 1 causes the As(O)OH fragment to lose its flatness. The solvent polar medium introduces significant changes in the geometry and IR spectral parameters of hydrogen-bonded complexes too: as the polarity of a medium increases, weak hydrogen bonds become weaker, and strong and medium hydrogen bonds become stronger; in the case of a complex with two hydrogen bonds cooperativity effects are observed. In almost all cases the driving force of these changes appears to be preferential solvation of charge-separated structures. In the limiting case of complete deprotonation (or conversely complete protonation) the vibrational frequencies of ν_AsO and ν_As-O turn into ν_As-O(asym) and ν_As-O(sym), respectively. In the intermediate cases the distance between ν_AsO and ν_As-O is sensitive to both implicit solvation and explicit solvation and the systematic changes of this distance can be used for estimation of the degree of proton transfer within the hydrogen bond.

Metal-Assisted Complexation of Fluorogenic Dyes by Cucurbit[7]uril and Cucurbit[8]uril: A DFT Evaluation of the Key Factors Governing the Host-Guest Recognition

With the emergence of host-guest systems, a novel branch of complexation chemistry has found wide application in industries such as food, pharmacy, medicine, environmental protection and cosmetics. Along with the extensively studied cyclodextrins and calixarenes, the innovative cucurbiturils (CB) have enjoyed increased popularity among the scientific community as they possess even better qualities as cavitands as compared to the former molecules. Moreover, their complexation abilities could further be enhanced with the assistance of metal cations, which can interestingly exert a dual effect on the complexation process: either by competitively binding to the host entity or cooperatively associating with the CB@guest structures. In our previous work, two metal species (Mg²⁺ and Ga³⁺) have been found to bind to CB molecules in the strongest fashion upon the formation of host-guest complexes. The current study focuses on their role in the complex formation with three dye molecules: thiazole orange, neutral red, and thioflavin T. Various key factors influencing the process have been recognized, such as pH and the dielectric constant of the medium, the cavity size of the host, Mⁿ⁺ charge, and the presence/absence of hydration shell around the metal cation. A well-calibrated DFT methodology, solidly based and validated and presented in the literature experimental data, is applied. The obtained results shed new light on several aspects of the cucurbituril complexation chemistry.

DFT calculations in solution systems: solvation energy, dispersion energy and entropy

DFT calculations of reaction mechanisms in solution have always been a hot topic, especially for transition-metal-catalyzed reactions. The calculation of solvation energy is performed using either the polarizable continuum model (PCM) or the universal solvation model SMD. The PCM calculation is very sensitive to the choice of atomic radii to form a cavity, where the self-consistent isodensity PCM (SCI-PCM) has been recognized as the best choice and our IDSCRF radii can provide a similar cavity. Moving from a gas-phase case to a solution case, dispersion energy and entropy should be carefully treated. The solvent-solute dispersion is also important in solution systems, and it should be calculated together with the solute dispersion. Only half of the solvent-solute dispersion energy from the PCM calculation belongs to the solute molecules to maintain a thermal equilibrium between a solute molecule and its cavity, similar to the treatment of electrostatic energy. Relative solute dispersion energy should also be shared equally with the newly formed cavity. The entropy change from a gas phase to a liquid phase is quite large, but the modern quantum chemistry programs can only calculate the gas-phase translational entropy based on the idea-gas equation. In this review, we will provide an operable method to calculate the solution translational entropy, which has been coded in our THERMO program.

Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach

The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H₂O-···-H₂O), the free energy to transform an H₃O⁺ cation into an H₅O₂⁺ cation, and the hydrogen bond free energy of protonated water clusters (H₃O⁺-···-H₂O and H₅O₂⁺-···-H₂O). The computed equilibrium hydrogen bond free energy of H₂O-···-H₂O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H₃O⁺ by roughly 1-2 k_BT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H⁺ cation and H₂O in solution were found to be stronger than between two H₂O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.

Thermodynamics of π–π Interactions of Benzene and Phenol in Water

The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.

Towards a transferable nonelectrostatic model for continuum solvation: The electrostatic and nonelectrostatic energy correction model

In this work, we introduce an electrostatic and non-electrostatic (ENE) correction to the solvation energy based on the Solvent-Accessible Surface Area (SASA) of the solute and the solvent static dielectric constant. The proposed correction was developed for neutral solutes in non-aqueous solvents, considering three different implicit solvation models based on a Self-Consistent Reaction Field treatment of solute-solvent mutual polarization using an Apparent Surface Charge formalism, namely the Integral Equation Formalism of the Polarizable Continuum Model using a continuous surface charge scheme (PCM), the Solvation Model based on solute electron density (SMD), and the generalized Finite-Difference Poisson-Boltzmann (FDPB) model. The proposed correction was parametrized on a diverse training set of 4980 solvation data from the Solv@tum database of experimental solvation energies, and validated on the non-aqueous subset of the MNSOL database comprising 2140 solvation energies. The performances of the proposed ENE models with minimal and extended parameters formulations have been analyzed and the latter variant has been further compared to the widely used Cavity, Dispersion, and Solvent structural effects (CDS) non-electrostatic model originally developed for the SMx family of implicit solvation models. Overall, a very good agreement between the computed solvation energies with the ENE correction and the reference experimental data has been found on both the training and test sets for all continuum solvation models considered. Furthermore, results for the ENE correction are on par with the reference CDS non-electrostatic model for both SMD and FDPB electrostatics, but with the advantage of using a lower number of parameters and thus an improved transferability between different electrostatics treatments.

Theoretical Insight into the Phosphate-Targeted Silver's Antibacterial Action: Differentiation between Gram (+) and Gram (-) Bacteria

Although silver is one of the first metals finding broad applications in everyday life, specific key points of the intimate mechanism of its bacteriostatic/bactericidal activity lack explanation. It is widely accepted that the antimicrobial potential of the silver cation depends on the composition and thickness of the bacterial external envelope: the outer membrane in Gram-negative bacteria is more prone to Ag⁺ attack than the cell wall in Gram-positive bacteria. The major cellular components able to interact strongly with Ag⁺ (teichoic acids, phospholipids, and lipopolysaccharides) contain mono/diesterified phosphate moieties. By applying a reliable DFT/SMD methodology, we modeled the reactions between the aforementioned constituents in typical Gram-positive and Gram-negative bacteria and hydrated Ag⁺ species, thus disclosing the factors that govern the process of metal-model ligand complexation. The conducted research indicates thermodynamically possible reactions in all cases but still a greater preference of the Ag⁺ toward the constituents in Gram-negative bacteria in comparison with their counterparts in Gram-positive bacteria. The observed tendencies shed light on the specific interactions of the silver cation with the modeled phosphate-containing units at the atomic level.

Can a Wanzlick-like equilibrium exist between dicoordinate borylenes and diborenes?

Boron chemistry has experienced tremendous progress in the last few decades, resulting in the isolation of a variety of compounds with remarkable electronic structures and properties. Some examples are the singly Lewis-base-stabilised borylenes, wherein boron has a formal oxidation state of +I, and their dimers featuring a boron-boron double bond, namely diborenes. However, no evidence of a Wanzlick-type equilibrium between borylenes and diborenes, which would open a valuable route to the latter compounds, has been found. In this work, we combine DFT, coupled-cluster, multireference methods, and natural bond orbital/natural resonance theory analyses to investigate the electronic, structural, and kinetic factors controlling the reactivity of the transient CAAC-stabilised cyanoborylene, which spontaneously cyclotetramerises into a butterfly-type, twelve-membered (BCN)4 ring, and the reasons why its dimerisation through the boron atoms is hampered. The computations are also extended to the NHC-stabilised borylene counterparts. We reveal that the borylene ground state multiplicity dictates the preference for self-stabilising cyclooligomerisation over boron-boron dimerisation. Our comparison between NHC- vs. CAAC-stabilised borylenes provides a convincing rationale for why the reduction of the former always gives diborenes while a range of other products is found for the latter. Our findings provide a theoretical background for the rational design of base-stabilised borylenes, which could pave the way for novel synthetic routes to diborenes or alternatively non-dimerising systems for small-molecule activation.

Predicting the Solubility of Nonelectrolyte Solids Using a Combination of Molecular Simulation with the Solubility Parameter Method MOSCED: Application to the Wastewater Contaminants Monuron, Diuron, Atrazine and Atenolol

Methods to predict the equilibrium solubility of nonelectrolyte solids are indispensable for early-stage process development, design, and feasibility studies. Conventional analytic methods typically require reference data to regress parameters, which may not be available or limited for novel systems. Molecular simulation is a promising alternative, but is computationally intensive. Here, we demonstrate the ability to use a small number of molecular simulation free energy calculations to generate reference data to regress model parameters for the analytical MOSCED (modified separation of cohesive energy density) model. The result is an efficient analytical method to predict the equilibrium solubility of nonelectrolyte solids. The method is demonstrated for the wastewater contaminants monuron, diuron, atrazine and atenolol. Predictions for monuron, diuron and atrazine are in reasonable agreement with MOSCED parameters regressed using experimental solubility data. Predictions for atenolol are inferior, suggesting a potential limitation in the adopted molecular models, or the solvents selected to generate the necessary reference data.

Role of Explicit Hydration in Predicting the Aqueous Standard Reduction Potential of Sulfate Radical Anion by DFT and Insight into the Influence of pH on the Reduction Potential

Sulfate radical anion (SO₄^•-) is a potent oxidant capable of destroying recalcitrant environmental contaminants such as perfluoroalkyl carboxylic acids. In addition, it is thought to participate in important atmospheric reactions. Its standard reduction potential (E°) is fundamental to its reactivity. Using theoretical methods to accurately predict the aqueous phase E° requires solvation with explicit water molecules. Herein, using density functional theory, we calculated the aqueous E° of SO₄^•- and evaluated sensitivity to explicit water count. The E° increased considerably with more waters until ca. 24 were included, after which change in E° was small. When a proton was added to these systems, the E° was similar regardless of the explicit water count and this value was similar to the E° for systems with a large number of explicit waters but no proton. This result agrees with literature evidence that the E° is pH independent. Natural Bond Orbital natural population analysis indicated that in the case of both SO₄^2- and SO₄^•-, considerable charge was donated from the SO₄ center to the explicit solvation shells.

Accelerating Solvent Selection for Type II Porous Liquids

Type II porous liquids, comprising intrinsically porous molecules dissolved in a liquid solvent, potentially combine the adsorption properties of porous adsorbents with the handling advantages of liquids. Previously, discovery of appropriate solvents to make porous liquids had been limited to direct experimental tests. We demonstrate an efficient screening approach for this task that uses COSMO-RS calculations, predictions of solvent pK_a values from a machine-learning model, and several other features and apply this approach to select solvents from a library of more than 11,000 compounds. This method is shown to give qualitative agreement with experimental observations for two molecular cages, CC13 and TG-TFB-CHEDA, identifying solvents with higher solubility for these molecules than had previously been known. Ultimately, the algorithm streamlines the downselection of suitable solvents for porous organic cages to enable more rapid discovery of Type II porous liquids.

Acid Dissociation Constants of the Benzimidazole Unit in the Polybenzimidazole Chain: Configuration Effects

The acid dissociation constant of three benzimidazoles, namely 2,2′-bibenzo[d]imidazole, 2,5′-bibenzo[d]imidazole, and 5,5′-bibenzo[d]imidazole, have been investigated by means of density functional theory calculations in gas phase and in aqueous solution. The theoretical approach was validated by the comparing of predicted and experimentally determined pK_a values in imidazole, benzimidazole, and 2-phenylbenzimidazole. From the studied compounds, 2,2′-bibenzo[d]imidazole was found to be the most acidic, which made it a valuable candidate as a material for polymer electrolyte membrane fuel cells.

Unveiling the [3+2] cycloaddition between difluoromethyl diazomethane and 3-ylideneoxindole from the perspective of molecular electron density theory

Evolution of some key ELF basins along the IRC of the most favorable ortho/endo reaction path.

Piecewise Multipole-Expansion Implicit Solvation for Arbitrarily Shaped Molecular Solutes

The multipole-expansion (MPE) model is an implicit solvation model used to efficiently incorporate solvent effects in quantum chemistry. Even within the recent direct approach, the multipole basis used in MPE to express the dielectric response still solves the electrostatic problem inefficiently or not at all for solutes larger than approximately ten non-hydrogen atoms. In existing MPE parametrizations, the resulting systematic underestimation of the electrostatic solute-solvent interaction is presently compensated for by a systematic overestimation of nonelectrostatic attractive interactions. Even though the MPE model can thus reproduce experimental free energies of solvation of small molecules remarkably well, the inherent error cancellation makes it hard to assign physical meaning to the individual free-energy terms in the model, raising concerns about transferability. Here we resolve this issue by solving the electrostatic problem piecewise in 3D regions centered around all non-hydrogen nuclei of the solute, ensuring reliable convergence of the multipole series. The resulting method thus allows for a much improved reproduction of the dielectric response of a medium to a solute. Employing a reduced nonelectrostatic model with a single free parameter, in addition to the density isovalue defining the solvation cavity, our method yields free energies of solvation of neutral, anionic, and cationic solutes in water in good agreement with experiment.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Theoretical Study of a Two-Photon Fluorescent Probe Based on Nile Red Derivatives with Controllable Fluorescence Wavelength and Water Solubility

Hypochloric acid (HOCl) plays a vital role in the natural defense system, but abnormal levels of it can cause cell damage, accelerated human aging, and various diseases. It is of great significance to develop new probes for detecting HOCl in biosystems nondestructively and noninvasively. The purpose of this work is to explore new chemical modification strategies of two-photon excitation fluorescence (TPEF) probes to improve the poor water solubility and low efficiency in imaging applications. Nil-OH-6 has a two-photon absorption cross-section value as high as 243 GM and attains a good quantum yield of 0.49. In addition, the modification of terminal groups with different azetidine-heterospirocycles or N,N-dialkyl fused amino groups to Nile Red can effectively improve the fluorescence efficiency as well as increase the solubility to some extent. This study provides some strategies to simultaneously improve the fluorescence performance and solubility of these two-photon probes and, hence, reliable guidance and a foundation for the subsequent synthesis of TPEF probes based on Nile Red.

Corrections in the CHARMM36 Parametrization of Chloride Interactions with Proteins, Lipids, and Alkali Cations, and Extension to Other Halide Anions

The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl^- with proteins and lipids. Accordingly, when examined in solution, little to no Cl^-association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl^--specific transport protein CLC-ec1 shows irreversible Cl^- dissociation from the so-called S_cen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl^- complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F^- complexes being the least stable. The NBFIX corrections also preserve the Cl^- occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.

The Promoter Role of Amines in the Condensation of Silicic Acid: A First-Principles Investigation

Though well-recognized, the molecular-level understanding of the multifunctional roles of amines in the condensation of polysilicic acids, which is one of the key processes in hydrothermal synthesis of zeolites, is still limited. Taking ethylamine as a prototype, we investigated the mechanism of polysilicic acid condensation in the existence of organic amines in aqueous solution with extensive first-principles-based calculations. Because of the high proton affinity, ethylamine exists as amine silicates and alters the subsequent condensation mechanisms from a 1-step lateral attack mechanism accompanied with simultaneous intermolecular proton transfer in neutral aqueous solution to a 2-step S_N2-like mechanism. Specifically, the 5-coordinated Si species that were not observed on pathways of condensation in neutral solution are effectively stabilized by the ethylamine cations as intermediates, and the barriers for condensation of ortho-silicic acid are significantly reduced from 133 kJ/mol in neutral solution to 58 and 63 kJ/mol for formation of the 5-coordinated Si intermediate and proton transfer for water release, respectively. Similar variations of mechanisms and barriers for condensation were also observed in the formation of cyclic trimers as well as linear and cyclic tetramers of ortho-silicic acids. Based on these, it was proposed that apart from acting as structure-directing agents, pore fillers, and pH adjusters, organic amines can also function as promoters in the condensation of polysilicic acids.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

H2 activation by hydrogenase-inspired NiFe catalyst using frustrated Lewis pair: effect of buffer and halide ion in the heterolytic H-H bond cleavage

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production. We investigate the reaction mechanism of H₂ activation in an aqueous solution by the recently developed NiFe complex (Ogo et al. Sci. Adv. 2020, 6, eaaz8181) using density functional theory (DFT) calculation. Our computational results showed that H₂ is activated using frustrated Lewis pair. That is, H₂ binds to the Fe site of the NiFe complex, acting as a Lewis acid, while the added buffer, acting as Lewis base, abstracts protons to form a hydride complex. Furthermore, the higher basicity in the proton abstraction reaction characterises reaction more exergonic and lowers the reaction barrier. In addition, in the proton abstraction by the water molecule, the reaction barrier was lowered when anion such as Cl⁻ is in the vicinity of the water. Understanding the chemical species that contribute to the catalytic process in cooperation with the metal catalyst at the atomic level should help to maximise the function of the catalyst.

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production.

Rational Development of a Metal-Free Bifunctional System for the C-H Activation of Methane: A Density Functional Theory Investigation

The activation or heterolytic splitting of methane, a challenging substrate usually restricted to transition metals, has so far proven elusive in experimental frustrated Lewis pair (FLP) chemistry. In this article, we demonstrate, using density functional theory (DFT), that 1-aza-9-boratriptycene is a conceptually simple intramolecular FLP for the activation of methane. Systematic comparison with other FLP systems allows to gain insight into their reactivity with methane. The thermodynamics and kinetics of methane activation are interpreted by referring to the analysis of the natural charges and by employing the distortion-interaction/activation strain (DIAS) model. These showed that the nature of the Lewis base influences the selectivity over the reaction pathway, with N Lewis bases favoring the deprotonation mechanism and P bases the hydride abstraction one. The lower barrier of activation for 1-aza-9-boratriptycene and the higher products stability are due to a better interaction energy than its counterparts, itself due to electrostatic interactions with the methane moiety, favorable orbital overlaps allowed by the side-attack, and space proximity between the B and N atoms.

Robust and Efficient Implicit Solvation Model for Fast Semiempirical Methods

We present a robust and efficient method to implicitly account for solvation effects in modern semiempirical quantum mechanics and force fields. A computationally efficient yet accurate solvation model based on the analytical linearized Poisson-Boltzmann (ALPB) model is parameterized for the extended tight binding (xTB) and density functional tight binding (DFTB) methods as well as for the recently proposed GFN-FF general force field. The proposed methods perform well over a broad range of systems and applications, from conformational energies over transition-metal complexes to large supramolecular association reactions of charged species. For hydration free energies of small molecules, GFN1-xTB(ALPB) is reaching the accuracy of sophisticated explicitly solvated approaches, with a mean absolute deviation of only 1.4 kcal/mol compared to the experiment. Logarithmic octanol-water partition coefficients (log K_ow) are computed with a mean absolute deviation of about 0.65 using GFN2-xTB(ALPB) compared to experimental values indicating a consistent description of differential solvent effects. Overall, more than twenty solvents for each of the six semiempirical methods are parameterized and tested. They are readily available in the xtb and dftb+ programs for diverse computational applications.

W-RESP: Well-Restrained Electrostatic Potential-Derived Charges. Revisiting the Charge Derivation Model

Representation of electrostatic interactions by a Coulombic pairwise potential between atom-centered partial charges is a fundamental and crucial part of empirical force fields used in classical molecular dynamics simulations. The broad success of the AMBER force-field family originates mainly from the restrained electrostatic potential (RESP) charge model, which derives partial charges to reproduce the electrostatic field around the molecules. However, the description of the electrostatic potential around molecules by standard RESP may be biased for some types of molecules. In this study, we modified the RESP charge derivation model to improve its description of the electrostatic potential around molecules and thus electrostatic interactions in the force field. In particular, we reoptimized the atomic radii for definition of the grid points around the molecule, redesigned the restraining scheme, and included extra point (EP) charges. The RESP fitting was significantly improved for aromatic heterocyclic molecules. Thus, the suggested W-RESP(-EP) charge derivation model shows some potential for improving the performance of the nucleic acid force fields, for which the poor description of nonbonded interactions, such as the underestimated stability of base pairing, is well-established. We also report some preliminary simulation tests (around 1 ms of simulation data) on A-RNA duplexes, tetranucleotides, and tetraloops. The simulations reveal no adverse effects, while the description of base-pairing interactions might be improved. The new charges can thus be used in future attempts to improve the nucleic acid simulation force fields, in combination with reparametrization of the other terms.

Is natural fraxin an overlooked radical scavenger?

Fraxin (FX) (7-hydroxy-6-methoxycoumarin 8-glucoside) is a typical natural product of the coumarin family. This compound was shown to protect endothelial cells from oxidative stress; however, the nature of its antioxidant properties is still ambiguous. In this study, we report on a systematic evaluation of the radical scavenging activity of FX using a two-tier protocol based on thermodynamic and kinetic calculations. The results show that FX has moderate activity in the aqueous physiological environment against a range of radicals including HO˙, CCl3O˙, CCl3OO˙, NO2, , and HOO˙. The latter was examined in detail due to the prevalence of HOO˙ as a source of oxidative stress in biological systems. HOO˙ scavenging activity was promising in the gas phase but low in physiological environments with k overall = 1.57 × 106, 3.13 × 102 and 2.68 × 103 M-1 s-1 in the gas phase, pentyl ethanoate and water solvents, respectively. The formal hydrogen transfer mechanism at the O7-H bond dominates the hydroperoxyl radical scavenging of FX in the nonpolar media, whereas, in the polar environment, the activity is exerted by the single electron transfer mechanism of the anion state. This activity falls behind typical antioxidants such as Trolox, ascorbic acid, and trans-resveratrol under the studied conditions. Thus FX may have multiple health benefits, but it is not an outstanding natural antioxidant.

Implicit Solvents for the Polarizable Atomic Multipole AMOEBA Force Field

Computational protein design, ab initio protein/RNA folding, and protein-ligand screening can be too computationally demanding for explicit treatment of solvent. For these applications, implicit solvent offers a compelling alternative, which we describe here for the polarizable atomic multipole AMOEBA force field based on three treatments of continuum electrostatics: numerical solutions to the nonlinear and linearized versions of the Poisson-Boltzmann equation (PBE), the domain-decomposition conductor-like screening model (ddCOSMO) approximation to the PBE, and the analytic generalized Kirkwood (GK) approximation. The continuum electrostatics models are combined with a nonpolar estimator based on novel cavitation and dispersion terms. Electrostatic model parameters are numerically optimized using a least-squares style target function based on a library of 103 small-molecule solvation free energy differences. Mean signed errors for the adaptive Poisson-Boltzmann solver (APBS), ddCOSMO, and GK models are 0.05, 0.00, and 0.00 kcal/mol, respectively, while the mean unsigned errors are 0.70, 0.63, and 0.58 kcal/mol, respectively. Validation of the electrostatic response of the resulting implicit solvents, which are available in the Tinker (or Tinker-HP), OpenMM, and Force Field X software packages, is based on comparisons to explicit solvent simulations for a series of proteins and nucleic acids. Overall, the emergence of performative implicit solvent models for polarizable force fields opens the door to their use for folding and design applications.

Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges?

Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.

Splitting of multiple hydrogen molecules by bioinspired diniobium metal complexes: a DFT study

Splitting of molecular hydrogen (H₂) into bridging and terminal hydrides is a common step in transition metal chemistry. Herein, we propose a novel organometallic platform for cleavage of multiple H₂ molecules, which combines metal centers capable of stabilizing multiple oxidation states, and ligands bearing positioned pendant basic groups. Using quantum chemical modeling, we show that low-valent, early transition metal diniobium(ii) complexes with diphosphine ligands featuring pendant amines can favorably uptake up to 8 hydrogen atoms, and that the energetics are favored by the formation of intramolecular dihydrogen bonds. This result suggests new possible strategies for the development of hydrogen scavenger molecules that are able to perform reversible splitting of multiple H₂ molecules.

Towards a converged strategy for including microsolvation in reaction mechanism calculations

A major part of chemical conversions is carried out in the fluid phase, where an accurate modeling of the involved reactions requires to also take into account solvation effects. Implicit solvation models often cover these effects with sufficient accuracy but can fail drastically when specific solvent-solute interactions are important. In those cases, microsolvation, i.e., the explicit inclusion of one or more solvent molecules, is a commonly used strategy. Nevertheless, microsolvation also introduces new challenges-a consistent workflow as well as strategies how to systematically improve prediction performance are not evident. For the COSMO and COSMO-RS solvation models, this work proposes a simple protocol to decide if microsolvation is needed and how the corresponding molecular model has to look like. To demonstrate the improved accuracy of the approach, specific application examples are presented and discussed, i.e., the computation of aqueous pK_a values and a mechanistic study of the methanol mediated Morita-Baylis-Hillman reaction.