Classical Pauli repulsion: An anisotropic, atomic multipole model

Quantification of Anisotropy in Exchange and Dispersion Interactions: A Simple Model for Physics-Based Force Fields

In some compounds, exchange repulsion is orientation dependent. However, in contrast to quantum chemical methods that treat exchange explicitly, empirical models assume exchange to be spherically symmetric, yielding an average description only. Here we quantify the anisotropy of exchange and dispersion energy for hydrogen halides and water by probing these compounds with a helium atom using the symmetry-adapted perturbation theory (SAPT). The exchange interaction is reduced by up to 33% due to the σ-hole in hydrogen iodide, depending on the location of the probe. We demonstrate how this anisotropy can be modeled in empirical force fields either using an angle-dependent potential or by introducing virtual sites, reducing the error in the empirical model by a factor of 5 compared to isotropic atoms. Lone-pairs on water, positioned close to perpendicular to the plane of the molecule, on a line with the oxygen atom, and, surprisingly, σ-holes on water both modulate the exchange interaction strongly. Both lone-pairs and σ-holes can be modeled by virtual sites, leading to an 80% reduced error.

Completely Multipolar Model as a General Framework for Many-Body Interactions as Illustrated for Water

We introduce a general framework for many-body force fields, the Completely Multipolar Model (CMM), that utilizes multipolar electrical moments modulated by exponential decay of electron density as a common functional form for all terms of an energy decomposition analysis of intermolecular interactions. With this common functional form, the CMM model establishes well-formulated damped tensors that reach the correct asymptotes at both long- and short-range while formally ensuring no short-range catastrophes. CMM describes the separable EDA terms of dispersion, exchange polarization, and Pauli repulsion with short-ranged anisotropy, polarization as intramolecular charge fluctuations and induced dipoles, while charge transfer describes explicit movement of charge between molecules, and naturally describes many-body charge transfer by coupling into the polarization equations. We also utilize a new one-body potential that accounts for intramolecular polarization by including an electric field-dependent correction to the Morse potential to ensure that CMM reproduces all physically relevant monomer properties including the dipole moment, molecular polarizability, and dipole and polarizability derivatives. The quality of CMM is illustrated through agreement of individual terms of the EDA and excellent extrapolation to energies and geometries of an extensive validation set of water cluster data.

Understanding the Electronic Effects of Lewis Acid Catalysts in Accelerating Polar Diels-Alder Reactions

The electronic effects of Lewis acid (LA) catalysts in reducing the activation energies of polar Diels-Alder (P-DA) reactions have been studied within Molecular Electron Density Theory. To this end, a quantum topological energy partitioning scheme, namely, the Relative Interacting Atomic Energy (RIAE) analysis, is applied to the transition state structures (TSs) and the ground state of the reagents of two different LA-catalyzed P-DA reactions. Analyses of the ξEtotalX total energies of the two interacting frameworks f(X) show that the electronic energy stabilization of the electrophilic frameworks, resulting from the global electron density transfer (GEDT), is the cause of an effective decrease of the activation energies. On the other hand, an in-depth analysis of the ξEintraA intra-atomic energies of the atoms belonging to the electrophilic ethylenic framework in the LA-catalyzed P-DA reactions of cyclopentadiene with acrolein indicates that the strong electronic stabilization of the carbonyl carbon, resulting from the GEDT taking place at the TSs, is the main factor responsible for the decrease of the activation energies in these LA-catalyzed P-DA reactions. Finally, the increase in GEDT at the TSs of these P-DA reactions causes an increase in the larger C-C distance, resulting from the stabilization of the electrophilic framework, thereby decreasing the suggested Pauli repulsion.

A Simple Model for the Pauli Repulsion with Possible Utility in QM, MM and Chemical Education

The Pauli repulsion is the intermolecular force responsible for the volume and low compressibility of condensed-phase matter at normal conditions. A simple model for this force is presented, wherein per-atom electron densities are represented as spherical charge distributions that are prevented from significantly overlapping. In the example of two noble gas atoms approaching one another beyond their van der Waals radii, the distance between the centers of the electronic charge distributions becomes larger than the distance between the nuclei, giving rise to an unfavorable electrostatic interaction. For the purpose of calculating this interaction, the model is further simplified by representing the per-atom electron density as a negative point charge, loosely inspired by the classical Drude oscillator. The dispersion interaction is simplified to an R-6 term, centered on the aforementioned point charges. Despite the gross simplicity of the resulting formalism, near-quantitative agreement with high-level QM interaction energies of noble gas dimers is achieved. Accordingly, the present model is thought to have utility in force fields, in post-HF and post-DFT dispersion corrections, and in chemical education.

PHAHST Potential: Modeling Sorption in a Dispersion-Dominated Environment

PHAHST (potentials with high accuracy, high speed, and transferability) is a recently developed force field that utilizes exponential repulsion, multiple dispersion terms, explicit many-body polarization, and many-body van der Waals interactions. The result is a systematic approach to force field development that is computationally practical. Here, PHAHST is employed in the simulation for rare gas uptake of krypton and xenon in the metal-organic material, HKUST-1. This material has shown promise in use as an adsorptive separating agent and presents a challenge to model due to the presence of heterogeneous interaction sorption surfaces, which include pores with readily accessible, open-metal sites that compete with dispersion-dominated pores. Such environments are difficult to simulate with commonly used empirical force fields, such as the Lennard-Jones (LJ) potential, which perform better when electrostatics are dominant in determining the nature of sorption and alone are incapable of modeling interactions with open-metal sites. The effectiveness of PHAHST is compared to the LJ potential in a series of mixed Kr-Xe gas simulations. It has been demonstrated that PHAHST compares favorably with experimental results, and the LJ potential is inadequate. Overall, we establish that force fields with physically grounded repulsion/dispersion terms are required in order to accurately model sorption, as these interactions are an important component of the energy. Furthermore, it is shown that the simple mixing rules work nearly quantitatively for the true pair potentials, while they are not transferable for effective potentials like LJ.

Coupling molecular density functional theory with converged selected configuration interaction methods to study excited states in aqueous solution

This paper presents the first implementation of a coupling between advanced wavefunction theories and molecular density functional theory (MDFT). This method enables the modeling of solvent effect into quantum mechanical (QM) calculations by incorporating an electrostatic potential generated by solvent charges into the electronic Hamiltonian. Solvent charges are deduced from the spatially and angularly dependent solvent particle density. Such a density is obtained through the minimization of the functional associated with the molecular mechanics (MM) Hamiltonian describing the interaction between the fluid particles. The introduced QM/MDFT framework belongs to QM/MM family of methods, but its originality lies in the use of MDFT as the MM solver, offering two main advantages. First, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Second, it preserves a molecular-level description lost in macroscopic continuum approaches. The excited state properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. The excitation energies and dipole moments have been compared with experimental data and previous theoretical work. A key finding is that using the Hartree-Fock method to describe the solute allows for predicting the solvent charge around the ground state with sufficient precision for the subsequent SCI calculations of excited states. This significantly reduces the computational cost of the described procedure, paving the way for the study of more complex molecules.

Beyond isotropic repulsion: Classical anisotropic repulsion by inclusion of p orbitals

Accurate modeling of intermolecular repulsion is an integral component in force field development. Although repulsion can be explicitly calculated by applying the Pauli exclusion principle, this approach is computationally viable only for systems of limited sizes. Instead, it has previously been shown that repulsion can be reformulated in a "classical" picture: the Pauli exclusion principle prohibits electrons from occupying the same state, leading to a depletion of electronic charge between atoms, giving rise to an enhanced nuclear-nuclear electrostatic repulsion. This classical picture is called the isotropic S2/R approximation, where S is the overlap and R is the interatomic distance. This approximation accurately captures the repulsion of isotropic atoms such as noble gas dimers; however, a key deficiency is that it fails to capture the angular dependence of the repulsion of anisotropic molecules. To include directionality, the wave function must at least be a linear combination of s and p orbitals. We derive a new anisotropic S2/R repulsion model through the inclusion of the anisotropic p orbital term in the total wave function. Because repulsion is pairwise and decays rapidly, it can be truncated at a short range, making it amenable for efficient calculation of energy and forces in complex biomolecular systems. We present a parameterization of the S101 dimer database against the ab initio benchmark symmetry-adapted perturbation theory, which yields an rms error of only 0.9 kcal/mol. The importance of the anisotropic term is demonstrated through angular scans of water-water dimers and dimers involving halobenzene. Simulation of liquid water shows that the model can be computed efficiently for realistic system sizes.

Single-Molecule Exchange inside a Nanocage Provides Insights into the Origin of π-π Interactions

The attractive interactions between aromatic rings, also known as π-π interactions, have been widely used for decades. However, the origin of π-π interactions remains controversial due to the difficulties in experimentally measuring the weak interactions between π-systems. Here, we construct an elaborate system to accurately compare the strength of the π-π interactions between phenylalanine derivatives via molecular exchange processes inside a protein nanopore. Based on quantitative comparison of binding strength, we find that in most cases, the π-π interaction is primarily driven by dispersive attraction, with the electrostatic interaction playing a secondary role and tending to be repulsive. However, in cases where electronic effects are particularly strong, electrostatic induction may exceed dispersion forces to become the primary driving force for interactions between π-systems. The results of this study not only deepen our understanding of π-stacking but also have potential implications in areas where π-π interactions play a crucial role.

Anisotropic Metal-Metal Pauli Repulsion in Polynuclear d10 Metal Clusters

Metallophilicity has been widely considered to be the driving force for self-assembly of closed-shell d10 metal complexes, but this view has been challenged by recent studies showing that metallophilicity in linear d10-d10 dimers is repulsive. This is due to strong metal-metal (M-M') Pauli repulsion (Wan, Q., Proc. Natl. Acad. Sci. U. S. A. 2021, 118, e2019265118). Here, we study M-M' Pauli repulsion in d10 metal clusters. Our results show that M-M' Pauli repulsion in d10 polynuclear clusters is 6-52% weaker than in similar linear d10 complexes due to the anisotropic shape of (n+1)s-nd hybridized orbitals. The overall M-M' interactions in closed-shell d10 polynuclear metal clusters remain repulsive. The effects of coordination geometry, relativistic effects, and the ligand's electronegativity on M-M' Pauli repulsion in polynuclear d10 clusters have been explored. These findings provide valuable guidance for the design and development of ligands and coordination geometries that alleviate M-M' Pauli repulsion in d10 metal cluster systems.

A Machine Learning Force Field for Bio-Macromolecular Modeling Based on Quantum Chemistry-Calculated Interaction Energy Datasets

Accurate energy data from noncovalent interactions are essential for constructing force fields for molecular dynamics simulations of bio-macromolecular systems. There are two important practical issues in the construction of a reliable force field with the hope of balancing the desired chemical accuracy and working efficiency. One is to determine a suitable quantum chemistry level of theory for calculating interaction energies. The other is to use a suitable continuous energy function to model the quantum chemical energy data. For the first issue, we have recently calculated the intermolecular interaction energies using the SAPT0 level of theory, and we have systematically organized these energies into the ab initio SOFG-31 (homodimer) and SOFG-31-heterodimer datasets. In this work, we re-calculate these interaction energies by using the more advanced SAPT2 level of theory with a wider series of basis sets. Our purpose is to determine the SAPT level of theory proper for interaction energies with respect to the CCSD(T)/CBS benchmark chemical accuracy. Next, to utilize these energy datasets, we employ one of the well-developed machine learning techniques, called the CLIFF scheme, to construct a general-purpose force field for biomolecular dynamics simulations. Here we use the SOFG-31 dataset and the SOFG-31-heterodimer dataset as the training and test sets, respectively. Our results demonstrate that using the CLIFF scheme can reproduce a diverse range of dimeric interaction energy patterns with only a small training set. The overall errors for each SAPT energy component, as well as the SAPT total energy, are all well below the desired chemical accuracy of ~1 kcal/mol.

Intermolecular Non-Bonded Interactions from Machine Learning Datasets

Accurate determination of intermolecular non-covalent-bonded or non-bonded interactions is the key to potentially useful molecular dynamics simulations of polymer systems. However, it is challenging to balance both the accuracy and computational cost in force field modelling. One of the main difficulties is properly representing the calculated energy data as a continuous force function. In this paper, we employ well-developed machine learning techniques to construct a general purpose intermolecular non-bonded interaction force field for organic polymers. The original ab initio dataset SOFG-31 was calculated by us and has been well documented, and here we use it as our training set. The CLIFF kernel type machine learning scheme is used for predicting the interaction energies of heterodimers selected from the SOFG-31 dataset. Our test results show that the overall errors are well below the chemical accuracy of about 1 kcal/mol, thus demonstrating the promising feasibility of machine learning techniques in force field modelling.

Direct Observation of Twin van der Waals Molecular Chains

The van der Waals (vdW) assemblies are the most common structures of materials. However, direct mapping of intermolecular electron clouds of a vdW assembly has never been obtained, even though the intramolecular electron clouds were visualized by atomic-resolution techniques. In this report, we unprecedentedly mapped the intermolecular electron cloud of the assemblies of ethanol molecules via ethyl groups with high-resolution atomic force microscopy and scanning tunneling microscopy at 5 K, leading to the first visualization of vdW molecular chains, in which ethanol molecules assemble into twin vdW molecular chains in a reverse parallel configuration on the Ag(111) plane. Furthermore, spontaneous order-disorder transitions in the chain were dynamically observed, suggesting its unusual properties different from those of 2D vdW materials. These findings provide an "eye" to see the atomic world of vdW materials.

Accuracy of One- and Two-Photon Intensities with the Extended Polarizable Density Embedding Model

The recently developed extended polarizable density embedding (PDE-X) model is evaluated for the spectroscopic properties of organic chromophores solvated in water, including both one- and two-photon absorption properties. The PDE-X embedding model systematically improves vertical excitation energies over the preceding polarizable density embedding model (PDE). PDE-X shows more modest improvements over existing embedding models for oscillator strengths and two-photon absorption cross-sections, which are more sensitive properties. We argue that the origin of these discrepancies is related to the description of polarization effects, suggesting directions for future development of the embedding model.

Hybrid classical/machine-learning force fields for the accurate description of molecular condensed-phase systems

Electronic structure methods offer in principle accurate predictions of molecular properties, however, their applicability is limited by computational costs. Empirical methods are cheaper, but come with inherent approximations and are dependent on the quality and quantity of training data. The rise of machine learning (ML) force fields (FFs) exacerbates limitations related to training data even further, especially for condensed-phase systems for which the generation of large and high-quality training datasets is difficult. Here, we propose a hybrid ML/classical FF model that is parametrized exclusively on high-quality ab initio data of dimers and monomers in vacuum but is transferable to condensed-phase systems. The proposed hybrid model combines our previous ML-parametrized classical model with ML corrections for situations where classical approximations break down, thus combining the robustness and efficiency of classical FFs with the flexibility of ML. Extensive validation on benchmarking datasets and experimental condensed-phase data, including organic liquids and small-molecule crystal structures, showcases how the proposed approach may promote FF development and unlock the full potential of classical FFs.

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.

General Charge Transfer Dipole Model for AMOEBA-Like Force Fields

The development of highly accurate force fields is always an importance aspect in molecular modeling. In this work, we introduce a general damping-based charge transfer dipole (D-CTD) model to describe the charge transfer energy and the corresponding charge flow for H, C, N, O, P, S, F, Cl, and Br elements in common bio-organic systems. Then, two effective schemes to evaluate the charge flow from the corresponding induced dipole moment between the interacting molecules were also proposed and discussed. The potential applicability of the D-CTD model in ion-containing systems was also demonstrated in a series of ion-water complexes including Li⁺, Na⁺, K⁺, Mg²⁺, Ca²⁺, Fe²⁺, Zn²⁺, Pt²⁺, F^-, Cl^-, Br^-, and I^- ions. In general, the D-CTD model demonstrated good accuracy and good transferability in both charge transfer energy and the corresponding charge flow for a wide range of model systems. By distinguishing the intermolecular charge redistribution (charge transfer) under the influence of an external electric field from the accompanying intramolecular charge redistribution (polarization), the D-CTD model is theoretically consistent with current induced dipole-based polarizable dipole models and hence can be easily implemented and parameterized. Along with our previous work in charge penetration-corrected electrostatics, a bottom-up approach constructed water model was also proposed and demonstrated. The structure-maker and structure-breaker roles of cations and anions were also correctly reproduced using Na⁺, K⁺, Cl^-, and I^- ions in the new water model, respectively. This work demonstrates a cost-effective approach to describe the charge transfer phenomena. The water and ion models also show the feasibility of a modulated development approach for future force fields.

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

Embedding Beyond Electrostatics: The Extended Polarizable Density Embedding Model

The polarizable density embedding (PDE) model is a focused QM/QM fragment-based embedding model designed to model solvation effects on molecular properties. We extend the PDE model to include exchange and nonadditive exchange-correlation (for DFT) in the embedding potential in addition to the existing electrostatic, polarization, and nonelectrostatic effects already present. The resulting model, termed PDE-X, yields localized electronic excitation energies that accurately capture the range dependence of the solvent interaction and gives close agreement with full quantum mechanical (QM) results, even when using minimal QM regions. We show that the PDE-X embedding description consistently improves the accuracy of excitation energies for a diverse set of organic chromophores. The improved embedding description leads to systematic solvent effects that do not average out when applying configurational sampling.

The Stereoselective Total Synthesis of the Elusive Cephalosporolide F

Here we report a seven-step protecting-group-free stereoselective total synthesis of the elusive (+)-cephalosporolide F from d-glucose. A microwave-assisted reaction between the Meldrum's acid and the d-glucose to the respective octono-1,4-lactone derivative, and a low temperature visible-light photoredox spirocyclization of a chiral N-alkoxyphthalimide to ceph F, are the two key chemical reactions that allowed the accomplishment of this unprecedented feat under an environmentally friendly processes.

Charge Density Evolution Governing Interfacial Friction

It is well-known that the electron nature of a solid in contact plays a predominant role in determining the many properties of the contact systems, but the general rules of electron coupling that govern interfacial friction remain an open issue for the surface/interface community. Here, density functional theory calculations were used to investigate the physical origins of friction of solid interfaces. It was found that interfacial friction can be inherently traced back to the electronic barrier to the change in the contact configuration of the joints in slip due to the resistance of energy level rearrangement leading to electron transfer, which applies for various interface types ranging from van der Waals, metallic, and ionic to covalent joints. The variation of the electron density accompanying contact conformation changes along the sliding pathways is defined to track the frictional energy dissipation process occurring in slip. The results demonstrate that the frictional energy landscapes evolve synchronously with responding charge density evolution along sliding pathways, yielding an explicitly linear dependence of frictional dissipation on electronic evolution. The correlation coefficient enables us to interpret the fundamental concept of shear strength. The present charge evolution model thereby provides insights into the classic hypothesis that the friction force scales with the real contact area. This may shed light on the intrinsic origin of friction at the electronic level, opening the way to the rational design of nanomechanical devices as well as the understanding of the natural faults.

Modeling transport and filtration of nanoparticle suspensions in porous media

Recently membrane filters have gained in significance due to the need to provide protection against airborne pollution. A question of importance, and some controversy, is the efficiency of filters for small nanoparticles with diameters below 100 nm as these are considered particularly dangerous due to possible penetration into the lungs. The efficiency is measured by the number of particles blocked by the pore structure after passing though the filter. To study the penetration into pores by nanoparticles suspended in a fluid, a stochastic transport theory based on an atomistic model is used to calculate particle density and flow within the pores, resulting pressure gradient, and filter efficiency. The importance of pore size relative to particle diameter and of the parameters of the pore wall interactions are investigated. The theory is applied to aerosols in fibrous filters and found to reproduce common trends in measurements. As particles enter the initially empty pores on relaxation to the steady state the small penetration measured at the onset of filtration increases faster in time the smaller the nanoparticle diameter. Control of pollution by filtration is achieved by strong repulsion of pore walls for particle diameters greater than twice the effective pore width. For smaller nanoparticles the steady-state efficiency decreases as the pore wall interactions weaken. Effective efficiency is increased when the suspended nanoparticles inside the pores combine into clusters of sizes greater than the filter channel width.

Development and Implementation of Atomically Anisotropic First-Principles Force Fields: A Benzene Case Study

π-interactions are an important motif in chemical and biochemical systems. However, due to their anisotropic electron densities and complex balance of intermolecular interactions, aromatic molecules represent an ongoing challenge for accurate and transferable force field development. Historically, ab initio force fields for aromatics have not exhibited good accuracy with respect to bulk properties or have only been used to study gas-phase dimers. Using benzene as a proof of concept, herein we show how our own ab initio MASTIFF force field incorporates an atomically anisotropic description of intermolecular interactions to yield an accurate and robust model for aromatic interactions irrespective of phase. Compared to existing models, the MASTIFF benzene force field not only is accurate for liquid phase properties but also offers transferability to the gas and solid phases. Additionally, we introduce a computationally efficient OpenMM plugin which enables customizable anisotropic intermolecular functional forms and which can be generically used in any MD simulation where a model for nonspherical atomic features is required. Overall, our results demonstrate the importance of atomic-level anisotropy in enabling next-generation ab initio force field development.

Why much of Chemistry may be indisputably non-bonded?

In this compendium, the wide scope of all intermolecular interactions ever known has been revisited, in particular giving emphasis the capability of much of the elements of the periodic table to form non-covalent contacts. Either hydrogen bonds, dihydrogen bonds, halogen bonds, pnictogen bonds, chalcogen bonds, triel bonds, tetrel bonds, regium bonds, spodium bonds or even the aerogen bond interactions may be cited. Obviously that experimental techniques have been used in some works, but it was through the theoretical methods that these interactions were validate, wherein the QTAIM integrations and SAPT energy partitions have been useful in this regard. Therefore, the great goal concerns to elucidate the interaction strength and if the intermolecular system shall be total, partial or non-covalently bonded, wherein this last one encompasses the most majority of the intermolecular interactions what leading to affirm that chemistry is debatably non-bonded.

Regularized by Physics: Graph Neural Network Parametrized Potentials for the Description of Intermolecular Interactions

Simulations of molecular systems using electronic structure methods are still not feasible for many systems of biological importance. As a result, empirical methods such as force fields (FF) have become an established tool for the simulation of large and complex molecular systems. The parametrization of FF is, however, time-consuming and has traditionally been based on experimental data. Recent years have therefore seen increasing efforts to automatize FF parametrization or to replace FF with machine-learning (ML) based potentials. Here, we propose an alternative strategy to parametrize FF, which makes use of ML and gradient-descent based optimization while retaining a functional form founded in physics. Using a predefined functional form is shown to enable interpretability, robustness, and efficient simulations of large systems over long time scales. To demonstrate the strength of the proposed method, a fixed-charge and a polarizable model are trained on ab initio potential-energy surfaces. Given only information about the constituting elements, the molecular topology, and reference potential energies, the models successfully learn to assign atom types and corresponding FF parameters from scratch. The resulting models and parameters are validated on a wide range of experimentally and computationally derived properties of systems including dimers, pure liquids, and molecular crystals.

Reactivity and Enantioselectivity in NHC Organocatalysis Provide Evidence for the Complex Role of Modifications at the Secondary Sphere

Secondary-sphere interactions are often harnessed to control reactivity and selectivity in organometallic and enzymatic catalysis. Yet, such strategies have only recently been explicitly applied in the context of organocatalytic systems. Although increased stability, reproducibility, and selectivity were obtained in previous work using this approach, the precise mechanistic pathway promoted by secondary-sphere modification in organocatalysis remained unclear. Herein, we report a comprehensive mechanistic study on the origin of the unique reactivity patterns and stereocontrol observed with boronic acids (BAs) as secondary-sphere modifiers of N-heterocyclic carbene (NHC) organocatalysts. Kinetic experiments revealed partial order in catalyst upon the addition of BA and unusual preactivation behavior, indicating the presence of stable off-cycle catalyst aggregation and BA-base adducts. These hypotheses were supported both by computations and by a series of NMR and nonlinear effect experiments. Furthermore, computations indicated a rate-limiting, water-assisted hydrogen atom transfer mechanism. This finding led to a considerable enhancement in the experimental reaction rate while maintaining excellent enantioselectivity by adding catalytic amounts of water. Finally, computations and racemization experiments uncovered an uncommon Curtin-Hammett-controlled enantioselectivity in the presence of secondary-sphere modifiers.

Classical Exchange Polarization: An Anisotropic Variable Polarizability Model

Polarizability, or the tendency of the electron distribution to distort under an electric field, often depends on the local chemical environment. For example, the polarizability of a chloride ion is larger in gas phase compared to a chloride ion solvated in water. This effect is due to the restriction the Pauli exclusion principle places on the allowed electron states. Because no two electrons can occupy the same state, when a highly polarizable atom comes in close contact with other atoms or molecules, the space of allowed states can dramatically decrease. This constraint suggests that an accurate molecular mechanics polarizability model should depend on the radial distance between neighboring atoms. This paper introduces a variable polarizability model within the framework of the HIPPO (Hydrogen-like Intermolecular Polarizable Potential) force field, by damping the polarizability as a function of the orbital overlap of two atoms. This effectively captures the quantum mechanical exchange polarization effects, without explicit utilization of antisymmetrized wave functions. We show that the variable polarizability model remarkably improves the two-body polarization energies and three-body energies of ion-ion and ion-water systems. Under this model, no manual tuning of atomic polarizabilities for monatomic ions is required; the gas-phase polarizability can be used because an appropriate damping function is able to correct the polarizability at short range.

Comparison of approximate intermolecular potentials for ab initio fragment calculations on medium sized N-heterocycles

The ground state intermolecular potential of bimolecular complexes of N-heterocycles is analyzed for the impact of individual terms in the interaction energy as provided by various, conceptually different theories. Novel combinations with several formulations of the electrostatic, Pauli repulsion, and dispersion contributions are tested at both short- and long-distance sides of the potential energy surface, for various alignments of the pyrrole dimer as well as the cytosine-uracil complex. The integration of a DFT/CCSD density embedding scheme, with dispersion terms from the effective fragment potential (EFP) method is found to provide good agreement with a reference CCSD(T) potential overall; simultaneously, a quantum mechanics/molecular mechanics approach using CHELPG atomic point charges for the electrostatic interaction, augmented by EFP dispersion and Pauli repulsion, comes also close to the reference result. Both schemes have the advantage of not relying on predefined force fields; rather, the interaction parameters can be determined for the system under study, thus being excellent candidates for ab initio modeling.

Using atomic charges to model molecular polarization

We review different models for introducing electric polarization in force fields, with special focus on methods where polarization is modelled at the atomic charge level. While electric polarization has been included in several force fields, the common approach has been to focus on atomic dipole polarizability. Several approaches allow modelling electric polarization by using charge-flow between charge sites instead, but this has been less exploited, despite that atomic charges and charge-flow is expected to be more important than atomic dipoles and dipole polarizability. A number of challenges are required to be solved for charge-flow models to be incorporated into polarizable force fields, for example how to parameterize the models and how to make them computational efficient.

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by -95.2 kJ mol^-1 at the EEF strength 0.005 a.u.

Unravelling the Interactions of Magnetic Ionic Liquids by Energy Decomposition Schemes: Towards a Transferable Polarizable Force Field

This work aims at unravelling the interactions in magnetic ionic liquids (MILs) by applying Symmetry-Adapted Perturbation Theory (SAPT) calculations, as well as based on those to set-up a polarisable force field model for these liquids. The targeted MILs comprise two different cations, namely: 1-butyl-3-methylimidazolium ([Bmim]⁺) and 1-ethyl-3-methylimidazolium ([Emim]⁺), along with several metal halides anions such as [FeCl₄]⁻, [FeBr₄]⁻, [ZnCl₃]⁻ and [SnCl₄]²⁻ To begin with, DFT geometry optimisations of such MILs were performed, which in turn revealed that the metallic anions prefer to stay close to the region of the carbon atom between the nitrogen atoms in the imidazolium fragment. Then, a SAPT study was carried out to find the optimal separation of the monomers and the different contributions for their interaction energy. It was found that the main contribution to the interaction energy is the electrostatic interaction component, followed by the dispersion one in most of the cases. The SAPT results were compared with those obtained by employing the local energy decomposition scheme based on the DLPNO-CCSD(T) method, the latter showing slightly lower values for the interaction energy as well as an increase of the distance between the minima centres of mass. Finally, the calculated SAPT interaction energies were found to correlate well with the melting points experimentally measured for these MILs.

The Catalysis Effect of Na and Point Defect on NO Heterogeneous Adsorption on Carbon during High-Sodium Zhundong Coal Reburning: Structures, Interactions and Thermodynamic Characteristics

The reburning process in a furnace, a key way to reduce NOx emissions, is a heterogeneous reaction during coal combustion, in which the heterogeneous adsorption is dominant. Zhundong coal with a high content of alkali metal can enhance the reburning process. In this paper, the influence of sodium and a defect on NO heterogeneous adsorption was studied by the density functional theory, and the thermodynamic characteristic was also analyzed. The results indicate that the binding energy for NO adsorption on the pristine graphene surface (graphene-NO), Na-decorated pristine graphene surface (graphene-Na-NO), defect graphene surface (gsv-NO) and Na-decorated defect graphene (gsv-Na-NO) is −5.86, −137.12, −48.94 and −74.85 kJ/mol, respectively, and that the heterogeneous adsorption is an exothermic reaction. Furthermore, except for covalent bonds of C and N, C and O for gsv-NO, other interactions are a closed-shell one, based on the analysis of AIM, ELF and IGM. The area of electron localization for NO is graphene-Na-NO > gsv-Na-NO > gsv-NO > graphene-NO. The dispersion interaction is the main interaction force between NO and the pristine graphene surface. The δg index for the atom pairs about N–C and O–C on the pristine graphene surface is also the smallest. The density of spikes at graphene-Na-NO is bigger than that at gsv-Na-NO. Moreover, the thermodynamics characteristic showed that the reaction equilibrium constant of graphene-NO is less than those on the other surfaces under the same temperature. Thus, NO on the pristine graphene surface is the most difficult to adsorb, but the presence of sodium and a defect structure can promote its adsorption.

Modeling ligand electrochemical parameters by repulsion-corrected eigenvalues

Ligand electrochemical parameters, E_L , more commonly known as Lever parameters, have played a major research role in understanding redox processes involved in inorganic electrochemistry, enzymatic reactions, catalysis, solar cells, biochemistry, and materials science. Despite their broad usefulness, Lever parameters are not well understood at a first-principles level. Using density functional theory, we demonstrate in this contribution that a ligand's Lever parameter is fundamentally related to the ligand's ability to alter the eigenvalue of the electroactive spin-orbital in an octahedral transition metal complex. Our analysis furthers a first-principles understanding of the nature of Lever parameters.

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.

Development and Validation of a DFT-Based Force Field for a Hydrated Homoalanine Polypeptide

A new force field has been created for simulating hydrated alanine polypeptides using the adaptive force matching (AFM) method. Only density functional theory calculations using the Perdew-Burke-Ernzerhof exchange-correlation functional and the D3 dispersion correction were used to fit the force field. The new force field, AFM2020, predicts NMR scalar coupling constants for hydrated homopolymeric alanine in better agreements with experimental data than several other models including those fitted directly to such data. For Ala₇, the new force field shows about 15% helical conformations, 20% conformation in the β basin, and 65% polyproline II. The predicted helical population of short hydrated alanine is higher than previous estimates based on the same experimental data. Gas-phase simulations indicate that the force field developed by AFM solution-phase data is likely to produce a reasonable conformation distribution when hydration water is no longer present, such as the interior of a protein.

Evolution of Conformation, Nanomechanics, and Infrared Nanospectroscopy of Single Amyloid Fibrils Converting into Microcrystals

Nanomechanical properties of amyloid fibrils and nanocrystals depend on their secondary and quaternary structure, and the geometry of intermolecular hydrogen bonds. Advanced imaging methods based on atomic force microscopy (AFM) have unravelled the morphological and mechanical heterogeneity of amyloids, however a full understanding has been hampered by the limited resolution of conventional spectroscopic methods. Here, it is shown that single molecule nanomechanical mapping and infrared nanospectroscopy (AFM-IR) in combination with atomistic modelling enable unravelling at the single aggregate scale of the morphological, nanomechanical, chemical, and structural transition from amyloid fibrils to amyloid microcrystals in the hexapeptides, ILQINS, IFQINS, and TFQINS. Different morphologies have different Young's moduli, within 2-6 GPa, with amyloid fibrils exhibiting lower Young's moduli compared to amyloid microcrystals. The origins of this stiffening are unravelled and related to the increased content of intermolecular β-sheet and the increased lengthscale of cooperativity following the transition from twisted fibril to flat nanocrystal. Increased stiffness in Young's moduli is correlated with increased density of intermolecular hydrogen bonding and parallel β-sheet structure, which energetically stabilize crystals over the other polymorphs. These results offer additional evidence for the position of amyloid crystals in the minimum of the protein folding and aggregation landscape.

AMOEBA binding free energies for the SAMPL7 TrimerTrip host-guest challenge

As part of the SAMPL7 host-guest binding challenge, the AMOEBA force field was applied to calculate the absolute binding free energy for 16 charged organic ammonium guests to the TrimerTrip host, a recently reported acyclic cucurbituril-derived clip host structure with triptycene moieties at its termini. Here we report binding free energy calculations for this system using the AMOEBA polarizable atomic multipole force field and double annihilation free energy methodology. Conformational analysis of the host suggests three families of conformations that do not interconvert in solution on a time scale available to nanosecond molecular dynamics (MD) simulations. Two of these host conformers, referred to as the "indent" and "overlap" structures, are capable of binding guest molecules. As a result, the free energies of all 16 guests binding to both conformations were computed separately, and combined to produce values for comparison with experiment. Initial ranked results submitted as part of the SAMPL7 exercise had a mean unsigned error (MUE) from experimental binding data of 2.14 kcal/mol. Subsequently, a rigorous umbrella sampling reference calculation was used to better determine the free energy difference between unligated "indent" and "overlap" host conformations. Revised binding values for the 16 guests pegged to this umbrella sampling reference reduced the MUE to 1.41 kcal/mol, with a correlation coefficient (Pearson R) between calculated and experimental binding values of 0.832 and a rank correlation (Kendall τ) of 0.65. Overall, the AMOEBA results demonstrate no significant systematic error, suggesting the force field provides an accurate energetic description of the TrimerTrip host, and an appropriate balance of solvation and desolvation effects associated with guest binding.

Next-Generation Accurate, Transferable, and Polarizable Potentials for Material Simulations

PHAHST (potentials with high accuracy, high speed, and transferability) intermolecular potential energy functions have been developed from first principles for H₂, N₂, the noble gases, and a metal-organic material, HKUST-1. The potentials are designed from the outset to be transferable to heterogeneous environments including porous materials, interfaces, and material simulations. This is accomplished by theoretically justified choices for all functional forms, parameters, and mixing rules, including explicit polarization in every environment and fitting to high quality electronic structure calculations using methods that are tractable for real systems. The models have been validated in neat systems by comparison to second virial coefficients and bulk pressure-density isotherms. For inhomogeneous applications, our main target, comparisons are presented to previously published experimental studies on the metal-organic material HKUST-1 including adsorption, isosteric heats of adsorption, binding site locations, and binding site energies. A systematic prescription is provided for developing compatible potentials for additional small molecules and materials. The resulting models are recommended for use in complex heterogeneous simulations where existing potentials may be inadequate.

Protobranching as repulsion-induced attraction: a prototype for geminal stabilization

Noncovalent interactions are traditionally defined within the context of their attractive components, such as electrostatics and dispersion. Sources of molecular strain are derived through the destabilization of Coulombic and exchange repulsion. Due to this binary designation, the underlying origin of geminal stability with respect to alkanes (referred to as protobranching) has been an active subject for debate between these competing perspectives. We recast this stabilization as a complementary (Gestalt) interaction between dispersion and exchange repulsion, each impacting the other. We use triplet hydrogen and argon dimer as foundational van der Waals adducts to develop a procedure for the visualization and quantification of both exchange repulsion, Δρ_SCF, and medium-range correlation, ΔΔρ, as perturbations in electron density. We use the framework of the DFT-D3 correction to reproduce the shape of the dispersion potential at medium range and successfully model the trend in stability for the eighteen isomers of octane with a diverse series of functionals: BLYP, B3LYP, BP86, PBE, and PBE0. Collectively, our findings show that protobranching is a manifestation of steric repulsion-reduction in vibrational enthalpy and medium-range electron correlation.

Fragment Exchange Potential for Realizing Pauli Deformation of Interfragment Interactions

In fragment-based methods, the lack of explicit short-range exchange interactions between monomers can result in unphysical deformation in charge density. In this study, we describe a fragment exchange potential (XFP) to explicitly account for interfragmental Pauli deformation. In our implementation, a Kohn-Sham exchange potential is adopted along with the Yukawa potential. The method has been validated by comparison of the computed exchange energies using the XFP potential with results obtained from antisymmetrized fragmental orbitals on the S66×8 data set containing 528 bimolecular interactions of equilibrium and arbitrary geometries. It was also found that it is only necessary to deploy numerical grids on atoms within their van der Waals contacts, significantly reducing the small, albeit extra, computational cost. We anticipate that the XFP presented here may be applied to molecular dynamics simulations of macromolecules using a fragment-based quantum mechanical potential with improved SCF convergence and computational accuracy.

SAPT codes for calculations of intermolecular interaction energies

Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.

Psi4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

Accurate Biomolecular Simulations Account for Electronic Polarization

In this perspective, we discuss where and how accounting for electronic many-body polarization affects the accuracy of classical molecular dynamics simulations of biomolecules. While the effects of electronic polarization are highly pronounced for molecules with an opposite total charge, they are also non-negligible for interactions with overall neutral molecules. For instance, neglecting these effects in important biomolecules like amino acids and phospholipids affects the structure of proteins and membranes having a large impact on interpreting experimental data as well as building coarse grained models. With the combined advances in theory, algorithms and computational power it is currently realistic to perform simulations with explicit polarizable dipoles on systems with relevant sizes and complexity. Alternatively, the effects of electronic polarization can also be included at zero additional computational cost compared to standard fixed-charge force fields using the electronic continuum correction, as was recently demonstrated for several classes of biomolecules.

Molecular and Supramolecular Interactions in Systems with Nitroxide-Based Radicals

Nitroxide-based radicals, having the advantage of firm chemical stability, are usable as probes in the detection of nanoscale details in the chemical environment of various multi-component systems, based on subtle variations in their electron paramagnetic resonance spectra. We propose a systematic walk through the vast area of problems and inquires that are implied by the rationalization of solvent effects on the spectral parameters, by first-principle methods of structural chemistry. Our approach consists of using state-of-the-art procedures, like Density Functional Theory (DFT), on properly designed systems, kept at the border of idealization and chemical realism. Thus, we investigate the case of real solvent molecules intervening in different configurations between two radical molecules, in comparison with radicals taken in vacuum or having the solvent that is treated by surrogate models, such as polarization continuum approximation. In this work, we selected the dichloromethane as solvent and the prototype radicals abbreviated TEMPO ((2,2,6,6-Tetramethylpiperidin-1-yl) oxyl). In another branch of the work, we check the interaction of radicals with large toroidal molecules, β-cyclodextrin, and cucurbit[6]uril, modeling the interaction energy profile at encapsulation. The drawn synoptic view offers valuable rationales for understanding spectroscopy and energetics of nitroxide radicals in various environments, which are specific to soft chemistry.