Search

Combining classical density functional theory (cDFT) with quantum mechanics (QM) methods offers a computationally efficient alternative to traditional QM/molecular mechanics (MM) approaches for modeling mixed quantum-classical systems at finite temperatures. However, both QM/MM and QM/cDFT rely on approximations that introduce ambiguities. This series of papers aims to establish a rigorous theoretical framework to guide the development of a QM/cDFT formulation that clarifies the approximations involved. In this first paper, we establish a comprehensive density functional theory (DFT) framework for mixed quantum (QM)-classical (MM) systems within the canonical ensemble. We derive the equilibrium density matrix for a system comprising N qm QM and N mm MM particles and, based on this derivation, propose a variational formulation of the Helmholtz free energy in terms of the full QM-MM density matrix. By leveraging permutational symmetry and thanks to constrained-search methods, we reformulate the variational problem using functionals of objects having reduced dimensionalities. Three different formulations are presented, each using quantities with progressively lower dimensionalities. This ultimately enables computation of the Helmholtz free energy using only the QM and MM one-body densities, generalizing both cDFT and electronic DFT (eDFT) for QM-MM systems. In a subsequent paper, this formalism will be extended to the grand canonical ensemble, allowing classical particle number fluctuations. This extension, in which cDFT is most naturally formulated, provides a framework for studying solvation problems.

This paper presents the first implementation of a coupling between advanced wave function 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 density is obtained through the minimization of the functional associated to 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. Firstly, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Secondly, it preserves a molecular-level description lost in macroscopic continuum approaches. Excited states properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. Excitation energies and dipole moment 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., Comment: 10 pages, 1 figure

We re-examine the problem of the dielectric response of highly polar liquids such as water in confinement between two walls using a simple two-variable density functional theory involving number and polarisation densities. In the longitudinal polarisation case where a perturbing field is applied perpendicularly to the walls, we show that the notion of local dielectric constant, although ill-defined at a microscopic level, makes sense when a coarse-graining over the typical size of a particle is introduced. The approach makes it possible to study the effective dielectric response of thin liquid films of various thicknesses in connection to the recent experiments of [Fumagalli et al. , Science, 2018, 360, 1339-1342], and to discuss the notion interfacial dielectric constant. We argue that the observed properties as function of slab dimension, in particular the very low dielectric constants of the order of 2-3 measured for thin slabs of 1 nm thickness do not highlight any special property of water but can be recovered for a generic polar solvent having similar particle size and the same high dielectric constant. Regarding the transverse polarisation case where the perturbing field is parallel to the walls, the associated effective dielectric constant as a function of the slab dimension reaches bulk-like values at much shorter widths than in the longitudinal case., Comment: 10 pages, 9 figures

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic Density Functional Theory and Molecular Dynamics. Here we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference Quantum Monte Carlo calculations that allow us to benchmark various electronic Density Functional Theory functionals. We then fit an accurate carbon--lithium pair potential, which is used in molecular Density Functional Theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer., Comment: 10 pages 4 figures

Liquid states theories have emerged as a numerically efficient alternative to costly molecular dynamics simulations of electron transfer reactions in solution. In a recent paper [Chem. Sci., 2019, 10, 2130], we introduced the framework to compute energy gap, free energy profile and reorganization free energy using molecular density functional theory. However, this technique, as other molecular liquid state theories, overestimates the bulk pressure of the fluids. Because of the too high pressure, the predicted free energy is dramatically exaggerated. Several attempts were made to fix this issue, either based on simple a posteriori correction or by improving the description of the liquid introducing bridge terms. By studying two model half reactions in water, Cl -> Cl+ and Cl -> Cl-, we assess the correctness of these two types of corrections to study electron transfer reactions. We found that a posteriori corrections, because they violate the functional principle, lead to an inconsistency in the definition of the reorganization free energy and should not be used to study electron transfer reactions. The bridge approach, because it is theoretically well grounded, is perfectly suitable for this type of systems., Comment: 15 pages, 4 figures

Solvation effect might have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike previous work, both densities are minimized self consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effect in quantum calculation. It is computationally more effective than a direct quantum-mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway., Comment: 14 pages, 8 figures

Biredox ionic liquids are a new class of functionalized electrolytes that may play an important role in future capacitive energy storage devices. By allowing additional storage of electrons inside the liquids, they can improve device performance significantly. However current devices employ nanoporous carbons in which the diffusion of the liquid and the adsorption of the ions could be affected by the occurrence of electron-transfer reactions. It is therefore necessary to understand better the thermodynamics and the kinetics of such reactions in biredox ionic liquids. Here we perform ab initio molecular dynamics simulations of both the oxidized and reduced species of several redox-active ionic molecules (used in biredox ionic liquids) dissolved in acetonitrile solvent and compare them with the bare redox molecules. We show that in all the cases, it is necessary to introduce a two Gaussian state model to calculate the reaction free energies accurately. These reaction free energies are only slightly affected by the presence of the IL group on the molecule. We characterize the structure of the solvation shell around the redox active part of the molecules and show that in the case of TEMPO-based molecules strong reorientation effects occur during the oxidation reaction., Comment: 11 pages, 12 figures

A better understanding of interfacial mechanisms is needed to improve the performances of electrochemical devices. Yet, simulating an electrode surface at fixed electrolyte composition remains a challenge. Here we apply a finite electric field to a single electrode held at constant potential and in contact with an aqueous ionic solution, using classical molecular dynamics. The polarization yields two electrochemical interfaces on opposite sides of the same metal slab. While the net charge on one electrode surface is the opposite of the net charge on the other, maintaining overall charge neutrality of the metal. The electrode surface charges fluctuations are compensated by the adsorption of ions from the electrolyte, forming a pair of electric double layers with aligned dipoles. This opens the way towards the efficient simulation of electrochemical interfaces using any flavor of molecular dynamics, from classical to first principles-based methods., Comment: 6 pages, 7 figures (including SI)

Most of the performances of electrochemical devices are governed by molecular processes taking place at the solution-electrode interfaces and molecular simulation are the main way to study these processes. Aqueous electrochemical systems have often been studied using classical DFT but with too crude approximations to consider the system description to be realistic. We study the interface between graphene electrodes and liquid water at different applied voltage using molecular DFT, improving the state of the art by the following key points: 1) electrodes have a realistic atomic resolution, 2) classical DFT calculations are carried out at fixed imposed potential difference and 3) water is described by a molecular model. This allows to reveal the structural modification of water adsorbed at the graphene interface and the evolution of water dielectric permittivity when a voltage is applied. The computed capacitance of this device is in agreement with molecular dynamics simulations. This demonstrates the relevance of molecular DFT to study electrochemical systems at the molecular level., Comment: 6 pages, 5 figures, Supplementary infos

Beyond the dielectric continuum description initiated by Marcus theory, the nowadays standard theoretical approach to study electron transfer (ET) reactions in solution or at interfaces is to use classical force field or ab initio Molecular Dynamics simulations. We propose here an alternative method based on liquid-state theory, namely molecular density functional theory, which is numerically much more efficient than simulations while still retaining the molecular nature of the solvent. We begin by reformulating molecular ET theory in a density functional language and show how to compute the various observables characterizing ET reactions from an ensemble of density functional minimizations. In particular, we define in that formulation the relevant order parameter of the reaction, the so-called vertical energy gap, and determine the Marcus free energy curves of both reactant and product states along that coordinate. Important thermodynamic quantities such as the reaction free energy and the reorganization free energies follow. We assess the validity of the method by studying the model Cl$^0\rightarrow$ Cl$^+$ and Cl$^0\rightarrow$ Cl$^-$ ET reactions in bulk water for which molecular dynamics results are available. The anionic case is found to violate the standard Marcus theory. Finally, we take advantage of the computational efficiency of the method to study the influence of confinement on the ET, by investigating the evolution of the reorganization free energy of the Cl$^0\rightarrow$ Cl$^+$ reaction when the atom approaches an atomistically resolved wall., Comment: 30 pages, 10 figures, Chemical Science, 2019

Nanoconfinement generally leads to drastic effect on the physical and chemical properties of ionic liquids. Here we investigate how the electrochemical reactivity in such media may be impacted inside nanoporous carbon electrodes. To this end, we study a simple electron transfer reaction using molecular dynamics simulations. The electrodes are held at constant electric potential by allowing the atomic charges on the carbon atoms to fluctuate. We show that the $\mathrm{Fe^{3+}}/\mathrm{Fe^{2+}}$ couple dissolved in an ionic liquid exhibits a deviation with respect to Marcus theory. This behavior is rationalized by the stabilization of a solvation state of the Fe$^{3+}$ cation in the disordered nanoporous electrode that is not observed in the bulk. The simulation results are fitted with a recently proposed two solvation state model, which allows us to estimate the effect of such a deviation on the kinetics of electron transfer inside nanoporous electrodes., Comment: 8 pages, 10 figures

In this article we report a stochastic evaluation of the recently proposed LCC multireference perturbation theory [Sharma S., and Alavi A., J. Chem. Phys. 143, 102815, (2015)]. In this method both the zeroth order and first order wavefunctions are sampled stochastically by propagating simultaneously two populations of signed walkers. The sampling of the zeroth order wavefunction follows a set of stochastic processes identical to the one used in the FCIQMC method. To sample the first order wavefunction, the usual FCIQMC algorithm is augmented with a source term that spawns walkers in the sampled first order wavefunction from the zeroth order wavefunction. The second order energy is also computed stochastically but requires no additional overhead outside of the added cost of sampling the first order wavefunction. This fully stochastic method opens up the possibility of simultaneously treating large active spaces to account for static correlation and recovering the dynamical correlation using perturbation theory. This method is used to study a few benchmark systems including the carbon dimer and aromatic molecules. We have computed the singlet-triplet gaps of benzene and m-xylylene. For m-xylylene, which has proved difficult for standard CASSCF+PT, we find the singlet-triplet gap to be in good agreement with the experimental values., Comment: 22 pages, 4 figures

We extend the recently proposed heat-bath configuration interaction (HCI) method [Holmes, Tubman, Umrigar, J. Chem. Theory Comput. 12, 3674 (2016)], by introducing a semistochastic algorithm for performing multireference Epstein-Nesbet perturbation theory, in order to completely eliminate the severe memory bottleneck of the original method. The proposed algorithm has several attractive features. First, there is no sign problem that plagues several quantum Monte Carlo methods. Second, instead of using Metropolis-Hastings sampling, we use the Alias method to directly sample determinants from the reference wavefunction, thus avoiding correlations between consecutive samples. Third, in addition to removing the memory bottleneck, semistochastic HCI (SHCI) is faster than the deterministic variant for many systems if a stochastic error of 0.1 mHa is acceptable. Fourth, within the SHCI algorithm one can trade memory for a modest increase in computer time. Fifth, the perturbative calculation is embarrassingly parallel. The SHCI algorithm extends the range of applicability of the original algorithm, allowing us to calculate the correlation energy of very large active spaces. We demonstrate this by performing calculations on several first row dimers including F2 with an active space of (14e, 108o), Mn-Salen cluster with an active space of (28e, 22o), and Cr2 dimer with up to a quadruple-zeta basis set with an active space of (12e, 190o). For these systems we were able to obtain better than 1 mHa accuracy with a wall time of merely 55 seconds, 37 seconds, and 56 minutes on 1, 1, and 4 nodes, respectively.

We present a three-dimensional molecular density functional theory (MDFT) of water derived from first-principles that relies on the particle's density and multipolar polarization density and includes the density-polarization coupling. This brings two main benefits: ($i$) a scalar density and a vectorial multipolar polarization density fields are much more tractable and give more physical insight than the full position and orientation densities, and ($ii$) it includes the full density-polarization coupling of water, that is known to be non-vanishing but has never been taken into account. Furthermore, the theory requires only the partial charge distribution of a water molecule and three measurable bulk properties, namely the structure factor and the Fourier components of the longitudinal and transverse dielectric susceptibilities., Comment: 11 figures

Solvation free energies are efficiently predicted by molecular density functionnal theory (MDFT) if one corrects the overpressure introduced by the usual homogeneous reference fluid approximation. Sergiievskyi et al. [Sergiievskyi et al., JPCL, 2014, 5, 1935-1942] recently derived the rigorous compensation of this excess of pressure (PC) and proposed an empirical "ideal gas" supplementary correction (PC+) that further enhances the calculated solvation free energies. In a recent paper [Misin et al, JCP, 2015, 142, 091105], those corrections were applied to solvation free energy calculations using the three-dimensional reference interaction site model (3D-RISM). As for classical DFT, PC and PC+ corrections improve greatly the predictions of 3D-RISM, but PC+ is described as decreasing the accuracy. In this article, we first derive rigorously the PC and PC+ corrections for 3D-RISM. We show the reported discrepancy is then taken off by introducing the correct expression of the pressure in 3D-RISM. This provides a consistent way to correct the solvation free-energies calculated by 3D-RISM method.

Low accuracy of the Solvation Free Energy (SFE) calculation is a known problem of the numerical methods of the Integral Equation Theory of Liquids and the Classical Density Functional Theory (Classical DFT). Although functionals with empirical corrections can essentially improve the predictability of the methods, their universality is still a question. In our recent paper we connected the SFE calculation errors with the incorrect pressure in the Classical DFT and proposed the a posteriory correction to improve the results (J. Phys. Chem. Lett., 5, 1925-1942 ). This paper raised a discussion in the community. In particular, recently appeared a critical reply where pointed some thermodynamical inconsistencies of the derivations in our paper (J. Chem. Theory Comput., 11, 378-380). In the present work we re-derive the pressure correction in a more simple way and show that despite the inaccuracies during the derivation, the final form of the previously derived correction is correct. We also test the applicability of the proposed correction to the functionals which include a three- and many- body terms from the fundamental measure theory (FMT) for hard sphere fluid. We test all the functionals on a set of model systems and discuss the obtained results.

We present a brief review of the classical density functional theory of atomic and molecular fluids. We focus on the application of the theory to the determination of the solvation properties of arbitrary molecular solutes in arbitrary molecular solvent. This includes the prediction of the solvation free energies, as well as the characterization of the microscopic, three-dimensional solvent structure., Comment: 27 pages, 7 figures, book chapter; IAS Series Vol. 28, 2015, Computational Trends in Solvation and Transport in Liquids, edited by Godehard Sutmann, Johannes Grotendorst, Gerhard Gompper, Dominik Marx ISBN 978-3-95806-030-2

The solvation of hydrophobic solutes in water is special because liquid and gas are almost at coexistence. In the common hypernetted chain approximation to integral equations, or equivalently in the homogenous reference fluid of molecular density functional theory, coexistence is not taken into account. Hydration structures and energies of nanometer-scale hydrophobic solutes are thus incorrect. In this article, we propose a bridge functional that corrects this thermodynamic inconsistency by introducing a metastable gas phase for the homogeneous solvent. We show how this can be done by a third order expansion of the functional around the bulk liquid density that imposes the right pressure and the correct second order derivatives. Although this theory is not limited to water, we apply it to study hydrophobic solvation in water at room temperature and pressure and compare the results to all-atom simulations. With this correction, molecular density functional theory gives, at a modest computational cost, quantitative hydration free energies and structures of small molecular solutes like n-alkanes, and of hard sphere solutes whose radii range from angstroms to nanometers. The macroscopic liquid-gas surface tension predicted by the theory is comparable to experiments. This theory gives an alternative to the empirical hard sphere bridge correction used so far by several authors., Comment: 18 pages, 6 figures

A classical density functional theory is applied to study solvation of solutes in water. An approx- imate form of the excess functional is proposed for water. This functional requires the knowledge of pure solvent direct correlation functions. Those functions can be computed by using molecular simulations such as molecular dynamic or Monte Carlo. It is also possible to use functions that have been determined experimentally. The functional minimization gives access to the solvation free energy and to the equilibrium solvent density. Some correction to the functional are also proposed to get the proper tetrahedral order of solvent molecules around a charged solute and to reproduce the correct long range hydrophobic behavior of big apolar solutes. To proceed the numerical minimization of the functional, the theory has been discretized on two tridimensional grids, one for the space coordinates, the other for the angular coordinates, in a functional minimization code written in modern Fortran, mdft. This program is used to study the solvation in water of small solutes of several kind, atomic and molecular, charged or neutral. More complex solutes, a neutral clay and a small protein have also been studied by functional minimization. In each case the classical density functional theory is able to reproduce the exact results predicted by MD. The computational cost is at least three order of magnitude less than in explicit methods., Comment: PhD Manuscript 157 pages, written in French

Molecular Density Functional Theory (MDFT) offers an efficient implicit- solvent method to estimate molecule solvation free-energies whereas conserving a fully molecular representation of the solvent. Even within a second order ap- proximation for the free-energy functional, the so-called homogeneous reference uid approximation, we show that the hydration free-energies computed for a dataset of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by two to three orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canoni- cal to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justifcation to empirical partial molar volume corrections that have been proposed recently., Comment: Version with correct equation numbers is here: http://compchemmpi.wikispaces.com/file/view/sergiievskyi_et_al.pdf/513575848/sergiievskyi_et_al.pdf Supporting information available online at: http://compchemmpi.wikispaces.com/file/view/SuppInf_sergiievskyi_et_al_07-04-2014.pdf/513576008/SuppInf_sergiievskyi_et_al_07-04-2014.pdf

We report here how the hydration of complex surfaces can be efficiently studied thanks to recent advances in classical molecular density functional theory. This is illustrated on the example of the pyrophylite clay. After presenting the most recent advances, we show that the strength of this implicit method is that (i) it is in quantitative or semi-quantitative agreement with reference all-atoms simulations (molecular dynamics here) for both the solvation structure and energetics, and that (ii) the computational cost is two to three orders of magnitude less than in explicit methods. The method remains imperfect, in that it locally overestimates the polarization of water close to hydrophylic sites of the clay. The high numerical efficiency of the method is illustrated and exploited to carry a systematic study of the electrostatic and van der Waals components of the surface-solvant interactions within the most popular force field for clays, CLAYFF. Hydration structure and energetics are found to weakly depend upon the electrostatics. We conclude on the consequences of such findings in future force-field development., Comment: 24 pages, 8 figures. Molecular Physics (2014)

We propose an in-silico experiment to introduce classical density functional theory (cDFT). Den- sity functional theories, whether quantum or classical, rely on abstract concepts that are non- intuitive. However, they are at the heart of powerful tools and active fields of research in both physics and chemistry. They led to the 1998 Nobel Prize in chemistry. DFT is illustrated here in its most simple and yet physically relevant form: the classical density functional theory of an ideal fluid of classical particles. For illustration purpose, it is applied to the prediction of the molecular structure of liquid neon. The numerical experiment proposed therein is built around the writing of a cDFT code by students in Mathematica. Students thus have to deal with (i) the cDFT theory, (ii) some basic concepts of statistical mechanics of simple fluids, (iii) functional minimization, and (iv) a useful functional programming language. This computational experiment is proposed during a molecular simulation class, but may also be of interest in a quantum chemistry class to illustrate electronic DFT, if one highlights the analogies between the quantum and classical DFTs., Comment: 4 pages, 2 figures, submitted to: The Journal of Chemical Education

We present an extension of our recently introduced molecular density functional theory of water [G. Jeanmairet et al., J. Phys. Chem. Lett. 4, 619, 2013] to the solvation of hydrophobic solutes of various sizes, going from angstroms to nanometers. The theory is based on the quadratic expansion of the excess free energy in terms of two classical density fields, the particle density and the multipolar polarization density. Its implementation requires as input a molecular model of water and three measurable bulk properties, namely the structure factor and the k-dependent longitudinal and transverse dielectric susceptibilities. The fine three-dimensional water structure around small hydrophobic molecules is found to be well reproduced. In contrast the computed solvation free-energies appear overestimated and do not exhibit the correct qualitative behavior when the hydrophobic solute is grown in size. These shortcomings are corrected, in the spirit of the Lum-Chandler-Weeks theory, by complementing the functional with a truncated hard-sphere functional acting beyond quadratic order in density. It makes the resulting functional compatible with the Van-der-Waals theory of liquid-vapor coexistence at long range. Compared to available molecular simulations, the approach yields reasonable solvation structure and free energy of hard or soft spheres of increasing size, with a correct qualitative transition from a volume-driven to a surface-driven regime at the nanometer scale., Comment: 24 pages, 8 figures

Three dimensional implementations of liquid state theories offer an efficient alternative to computer simulations for the atomic-level description of aqueous solutions in complex environments. In this context, we present a (classical) molecular density functional theory (MDFT) of water that is derived from first principles and is based on two classical density fields, a scalar one, the particle density, and a vectorial one, the multipolar polarization density. Its implementation requires as input the partial charge distribution of a water molecule and three measurable bulk properties, namely the structure factor and the k-dependent longitudinal and transverse dielectric constants. It has to be complemented by a solute-solvent three-body term that reinforces tetrahedral order at short range. The approach is shown to provide the correct three-dimensional microscopic solvation profile around various molecular solutes, possibly possessing H-bonding sites, at a computer cost two-three orders of magnitude lower than with explicit simulations., Comment: 19 pages, 4 figures

We show that classical molecular density functional theory (MDFT), here in the homogeneous reference fluid approximation in which the functional is inferred from the properties of the bulk solvent, is a powerful new tool to study, at a fully molecular level, the solvation of complex surfaces and interfaces by polar solvents. This implicit solvent method allows for the determination of structural, orientational and energetic solvation properties that are on a par with all-atom molecular simulations performed for the same system, while reducing the computer time by two orders of magnitude. This is illustrated by the study of an atomistically-resolved clay surface composed of over a thousand atoms wetted by a molecular dipolar solvent. The high numerical efficiency of the method is exploited to carry a systematic analysis of the electrostatic and non-electrostatic components of the surface-solvent interaction within the popular CLAYFF force field. Solvent energetics and structure are found to depend weakly upon the atomic charges distribution of the clay surface, even for a rather polar solvent. We conclude on the consequences of such findings for force-field development., Comment: accepted by J. Chem. Phys., nov. 2012