An efficient Monte Carlo method for calculating ab initio transition state theory reaction rates in solution

In this article, we propose an efficient method for sampling the relevant state space in condensed phase reactions. In the present method, the reaction is described by solving the electronic Schr\"{o}dinger equation for the solute atoms in the presence of explicit solvent molecules. The sampling algorithm uses a molecular mechanics guiding potential in combination with simulated tempering ideas and allows thorough exploration of the solvent state space in the context of an ab initio calculation even when the dielectric relaxation time of the solvent is long. The method is applied to the study of the double proton transfer reaction that takes place between a molecule of acetic acid and a molecule of methanol in tetrahydrofuran. It is demonstrated that calculations of rates of chemical transformations occurring in solvents of medium polarity can be performed with an increase in the cpu time of factors ranging from 4 to 15 with respect to gas-phase calculations.Comment: 15 pages, 9 figures. To appear in J. Chem. Phy

Analysis of the conformational profiles of fenamates shows route towards novel, higher accuracy, force-fields for pharmaceuticals

In traditional molecular mechanics force fields, intramolecular non-bonded interactions are modelled as intermolecular interactions, and the form of the torsion potential is based on the conformational profiles of small organic molecules. We investigate how a separate model for the intramolecular forces in pharmaceuticals could be more realistic by analysing the low barrier to rotation of the phenyl ring in the fenamates (substituted N-phenyl-aminobenzoic acids), that results in a wide range of observed angles in the numerous fenamate crystal structures. Although the conformational energy changes by significantly less than 10 kJmol-1 for a complete rotation of the phenyl ring for fenamic acid, the barrier is only small because of small correlated changes in the other bond and torsion angles. The maxima for conformations where the two aromatic rings approach coplanarity arise from steric repulsion, but the maxima when the two rings are approximately perpendicular arise from a combination of an electronic effect and intramolecular dispersion. Representing the ab initio conformational energy profiles as a cosine series alone is ineffective; however, combining a cos2ξ term to represent the electronic barrier with an intramolecular atom-atom exp-6 term for all atom pairs separated by three or more bonds (1-4 interactions) provides a very effective representation. Thus we propose a new, physically motivated, generic analytical model of conformational energy, which could be combined with an intermolecular model to form more accurate force-fields for modelling the condensed phases of pharmaceutical-like organic molecules

The Monomer Electron Density Force Field (MEDFF) : a physically inspired model for noncovalent interactions

We propose a methodology to derive pairwise-additive noncovalent force fields from monomer electron densities without any empirical input. Energy expressions are based on the symmetry-adapted perturbation theory (SAPT) decomposition of interaction energies. This ensures a physically motivated force field featuring an electrostatic, exchange repulsion, dispersion, and induction contribution, which contains two types of parameters. First, each contribution depends on several fixed atomic parameters, resulting from a partitioning of the monomer electron density. Second, each of the last three contributions (exchange-repulsion, dispersion, and induction) contains exactly one linear fitting parameter. These three so-called interaction parameters in the model are initially estimated separately using SAPT reference calculations for the S66x8 database of noncovalent dimers. In a second step, the three interaction parameters are further refined simultaneously to reproduce CCSD(T)/CBS interaction energies for the same database. The limited number of parameters that are fitted to dimer interaction energies (only three) avoids ill-conditioned fits that plague conventional parameter optimizations. For the exchange repulsion and dispersion component, good results are obtained for all dimers in the S66x8 database using one single value for the associated interaction parameters. The values of those parameters can be considered universal and can also be used for dimers not present in the original database used for fitting. For the induction component such an approach is only viable for the dispersion dominated dimers in the S66x8 database. For other dimers (such as hydrogen-bonded complexes), we show that our methodology remains applicable. However, the interaction parameter needs to be determined on a case-specific basis. As an external validation:, the force field predicts interaction energies in good agreement with CCSD(T)/CBS values for dispersion dominated dimers extracted from an HIV-II protease crystal structure with a bound ligand (indinavir). Furthermore, experimental second virial coefficients of small alkanes and alkenes are well reproduced

Range separation: The divide between local structures and field theories

This work presents parallel histories of the development of two modern theories of condensed matter: the theory of electron structure in quantum mechanics, and the theory of liquid structure in statistical mechanics. Comparison shows that key revelations in both are not only remarkably similar, but even follow along a common thread of controversy that marks progress from antiquity through to the present. This theme appears as a creative tension between two competing philosophies, that of short range structure (atomistic models) on the one hand, and long range structure (continuum or density functional models) on the other. The timeline and technical content are designed to build up a set of key relations as guideposts for using density functional theories together with atomistic simulation.Comment: Expanded version of a 30 minute talk delivered at the 2018 TSRC workshop on Ions in Solution, to appear in the March, 2019 issue of Substantia (https://riviste.fupress.net/index.php/subs/index

How Water's Properties Are Encoded in Its Molecular Structure and Energies.

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties

The genesis of the quantum theory of the chemical bond

An historical overview is given of the relevant steps that allowed the genesis of the quantum theory of the chemical bond, starting from the appearance of the new quantum mechanics and following later developments till approximately 1931. General ideas and some important details are discussed concerning molecular spectroscopy, as well as quantum computations for simple molecular systems performed within perturbative and variational approaches, for which the Born-Oppenheimer method provided a quantitative theory accounting for rotational, vibrational and electronic states. The novel concepts introduced by the Heitler-London theory, complemented by those underlying the method of the molecular orbitals, are critically analyzed along with some of their relevant applications. Further improvements in the understanding of the nature of the chemical bond are also considered, including the ideas of one-electron and three-electron bonds introduced by Pauling, as well as the generalizations of the Heitler-London theory firstly performed by Majorana, which allowed the presence of ionic structures into homopolar compounds and provided the theoretical proof of the stability of the helium molecular ion. The study of intermolecular interactions, as developed by London, is finally examined.Comment: amsart, 34 pages, 2 figure