Source author record

Daniel Borgis

Daniel Borgis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2020arXiv

Predicting hydration free energies of the FreeSolv database of druglike molecules with molecular density functional theory

We assess the performance of molecular densityfunctional theory (MDFT) to predict hydration freeenergies of the small drug-like molecules benchmark,FreeSolv. MDFT in the hyper-netted chain approx-imation (HNC) coupled with a pressure correctionpredicts experimental hydration free energies of theFreeSolv database within 1 kcal/mol with an averagecomputation time of two cpu.min per molecule. Thisis the same accuracy as for simulation based free en-ergy calculations that typically require hundreds ofcpu.h or tens of gpu.h per molecule.

preprint2020arXiv

Tackling solvent effect by coupling electronic and molecular Density Functional Theory

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.

preprint2016arXiv

Molecular density functional theory of water including density-polarization coupling

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.

preprint2015arXiv

Classical density functional theory to tackle solvation in molecular liquids

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.

preprint2015arXiv

Molecular Density Functional Theory for water with liquid-gas coexistence and correct pressure

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.

preprint2015arXiv

Pressure Correction in Classical Density Functional Theory: Hyper Netted Chain and Hard Sphere Bridge Functionals

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.

preprint2015arXiv

Solvation free-energy pressure corrections in the Three Dimensional Reference Interaction Site Model

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.

preprint2014arXiv

Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections

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.

preprint2014arXiv

Hydration of Clays at the Molecular Scale: The Promising Perspective of Classical Density Functional Theory

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.

preprint2014arXiv

Introduction to Classical Density Functional Theory by a Computational Experiment

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.

preprint2013arXiv

Molecular Density Functional Theory of Water

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.

preprint2013arXiv

Molecular Density Functional Theory of Water describing Hydrophobicity at Short and Long Length Scales

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.

preprint2012arXiv

Scalar fundamental measure theory for hard spheres in three dimensions. Application to hydrophobic solvation

Hard-sphere mixtures provide one a solvable reference system that can be used to improve the density functional theory of realistic molecular fluids. We show how the Kierlik-Rosinberg's scalar version of the fundamental measure density functional theory of hard spheres [Phys. Rev. A, {\bf 42}, 3382 (1990)], which presents computational advantages with respect to the original Rosenfeld's vectorial formulation or its extensions, can be implemented and minimized in three dimensions to describe fluid mixtures in complex environments. This implementation is used as a basis for defining a molecular density functional theory of water around molecular hydrophobic solutes of arbitrary shape.

preprint2012arXiv

Solvation of complex surfaces via molecular density functional theory

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.