BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Maximilien Levesque

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

ResearcherUnclaimed source record
physics.chem-phcond-mat.mtrl-scicond-mat.softcond-mat.stat-mechphysics.comp-phphysics.flu-dynArtificial IntelligenceQuantitative Methods

Catalog footprint

What is connected

26works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

26 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pressure Correction for Solvation Theories

Liquid state theories such as integral equations and classical density functional theory often overestimate the bulk pressure of fluids because they require closure relations or truncations of functionals. Consequently, the cost to create a molecular cavity in the fluid is no longer negligible and those theories predict wrong solvation free energies. We show how to correct them simply by computing an optimized Van der Walls volume of the solute and removing the undue free energy to create such volume in the fluid. Given this versatile correction, we demonstrate that state-of-the-art solvation theories can predict, within seconds, hydration free energies of a benchmark of small neutral drug-like molecules with the same accuracy as day-long molecular simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scaffold-constrained molecular generation

One of the major applications of generative models for drug Discovery targets the lead-optimization phase. During the optimization of a lead series, it is common to have scaffold constraints imposed on the structure of the molecules designed. Without enforcing such constraints, the probability of generating molecules with the required scaffold is extremely low and hinders the practicality of generative models for de-novo drug design. To tackle this issue, we introduce a new algorithm to perform scaffold-constrained in-silico molecular design. We build on the well-known SMILES-based Recurrent Neural Network (RNN) generative model, with a modified sampling procedure to achieve scaffold-constrained generation. We directly benefit from the associated reinforcement Learning methods, allowing to design molecules optimized for different properties while exploring only the relevant chemical space. We showcase the method's ability to perform scaffold-constrained generation on various tasks: designing novel molecules around scaffolds extracted from SureChEMBL chemical series, generating novel active molecules on the Dopamine Receptor D2 (DRD2) target, and, finally, designing predicted actives on the MMP-12 series, an industrial lead-optimization project.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Lattice Boltzmann Electrokinetics simulation of nanocapacitors

We propose a method to model metallic surfaces in Lattice Boltzmann Electrokinetics simulations (LBE), a lattice-based algorithm rooted in kinetic theory which captures the coupled solvent and ion dynamics in electrolyte solutions. This is achieved by a simple rule to impose electrostatic boundary conditions, in a consistent way with the location of the hydrodynamic interface for stick boundary conditions. The proposed method also provides the local charge induced on the electrode by the instantaneous distribution of ions under voltage. We validate it in the low voltage regime by comparison with analytical results in two model nanocapacitors: parallel plate and coaxial electrodes. We examine the steady-state ionic concentrations and electric potential profiles (and corresponding capacitance), the time-dependent response of the charge on the electrodes, as well as the steady-state electro-osmotic profiles in the presence of an additional, tangential electric field. The LBE method further provides the time-dependence of these quantities, as illustrated on the electro-osmotic response. While we do not consider this case in the present work, which focuses on the validation of the method, the latter readily applies to large voltages between the electrodes, as well as to time-dependent voltages. This work opens the way to the LBE simulation of more complex systems involving electrodes and metallic surfaces, such as sensing devices based on nanofluidic channels and nanotubes, or porous electrodes.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accounting for adsorption and desorption in lattice Boltzmann simulations

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. Associated to the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g. in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic, but also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a Lattice-Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account on the time-dependent diffusion coefficient in a more complex porous medium.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous Surface Segregation Profiles in Ferritic FeCr Stainless Steel

The iron-chromium alloy and its derivatives are widely used for their remarkable resistance to corrosion, which only occurs in a narrow concentration range around 9 to 13 atomic percent chromium. Although known to be due to chromium enrichment of a few atoms thick layer at the surfaces, the understanding of its complex atomistic origin has been a remaining challenge. We report an investigation of the thermodynamics of such surfaces at the atomic scale by means of Monte Carlo simulations. We use a Hamiltonian which provides a parameterization of previous ab initio results and successfully describes the alloy's unusual thermodynamics. We report a strong enrichment in Cr of the surfaces for low bulk concentrations, with a narrow optimum around 12 atomic percent chromium, beyond which the surface composition decreases drastically. This behavior is explained by a synergy between (i) the complex phase separation in the bulk alloy, (ii) local phase transitions that tune the layers closest to the surface to an iron-rich state and inhibit the bulk phase separation in this region, and (iii) its compensation by a strong and non-linear enrichment in Cr of the next few layers. Implications with respect to the design of prospective nanomaterials are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Démixtion et ségrégation superficielle dans les alliages fer-chrome

Ferritic steels possibly strengthened by oxide dispersion are candidates as structural materials for generation IV and fusion nuclear reactors. Their use is limited by incomplete knowledge of the iron-chromium phase diagram at low temperatures and of the phenomena inducing preferential segregation of one element at grain boundaries or at surfaces. In this context, this work contributes to the multi-scale study of the model iron-chromium alloy and their free surfaces by numerical simulations. This study begins with ab initio calculations of properties related to the mixture of atoms of iron and chromium. We highlight complex dependency of the magnetic moments of the chromium atoms on their local chemical environment. Surface properties are also proving sensitive to magnetism. This is the case of impurity segregation of chromium in iron and of their interactions near the surface. In a second step, we construct a simple energy model for high numerical efficiency. It is based on pair interactions on a rigid lattice to which are given local chemical environment and temperature dependencies. With this model, we reproduce the ab initio results at zero temperature and experimental results at high temperature. We also deduce the solubility limits at all intermediate temperatures with mean field approximations that we compare to Monte Carlo simulations. The last step of our work is to introduce free surfaces in our model. We then study the effect of ab initio calculated bulk and surface properties on surface segregation. Finally, we calculate segregation isotherms. We therefore propose an evolution model of surface composition of iron-chromium alloys as a function of bulk composition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of surface hydrogen on the anomalous surface segregation behavior of Cr in Fe-rich Fe-Cr alloys

The segregation behavior of Cr in dilute Fe-Cr alloys is known to be anomalous since the main barrier for surface segregation of Cr in these alloys arises not from the topmost surface layer but from the subsurface layer where the solution energy of Cr is much more endothermic as compared to the topmost surface layer. The Fe-Cr alloys are candidate structural materials for the new generation of nuclear reactors. The surfaces of these alloys will be exposed to hydrogen or its isotopes in these reactors, and although hydrogen is soluble neither in Fe nor in Fe-Cr alloys, it is known that the adsorption energy of hydrogen on the surface of iron is not only exothermic but relatively large. This clearly raises the question of the effect of the hydrogen adsorbed on the surface of iron on the segregation behavior of chromium towards the surface of iron. In this paper we show, on the basis of our ab initio density functional theory calculations, that the presence of hydrogen on the surface of iron leads to a considerably reduced barrier for Cr segregation to both the topmost surface layer and the subsurface layer, but the subsurface layer still controls the barrier for surface segregation. This reduction in the barrier for surface segregation is due to the nature of the Cr-H couple that acts in a complex and synergistic manner. The presence of Cr enhances the exothermic nature of hydrogen adsorption that in turn leads to a reduced barrier for surface segregation. These results should be included in the multiscale modeling of Fe-Cr alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electronic origin of the anomalous segregation behavior of Cr in Fe-rich Fe-Cr alloys

The energy of segregation of Cr from the bulk to the topmost surface layer in dilute Fe-Cr alloys is endothermic but small. On the other hand, the segregation energy from the bulk to the subsurface layer is not only endothermic but relatively large. Thus, the subsurface layer presents a huge barrier, not only for the segregation of Cr from the bulk to the surface, but also for its diffusion from the surface into the bulk. This means that the topmost layer alone cannot be used to determine the barrier for Cr segregation (the segregation energy), which is determined in these alloys by the subsurface layer. We present the results of our ab initio electronic structure calculations on the segregation of Cr as an impurity in Fe, which shed some light on the origin of this anomalous behavior. We find that the interaction of a Cr atom situated in the subsurface layer, in contrast to those in the other layers, is quite complex with its nearest neighbor Fe atoms due to the magnetism of both Fe and Cr, leading to a highly unfavorable electronic structure. This in turn leads to a large endothermic energy of solution of Cr in the subsurface layer. These results are important for a proper understanding of the properties of the technologically important class of Fe-Cr steels with low concentrations of chromium.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simple concentration-dependent pair interaction model for large-scale simulations of Fe-Cr alloys

This work is motivated by the need for large-scale simulations to extract physical information on the iron-chromium system that is a binary model alloy for ferritic steels used or proposed in many nuclear applications. From first-principles calculations and the experimental critical temperature we build a new energetic rigid lattice model based on pair interactions with concentration and temperature dependence. Density functional theory calculations in both norm-conserving and projector augmented-wave approaches have been performed. A thorough comparison of these two different ab initio techniques leads to a robust parametrization of the Fe-Cr Hamiltonian. Mean-field approximations and Monte Carlo calculations are then used to account for temperature effects. The predictions of the model are in agreement with the most recent phase diagram at all temperatures and compositions. The solubility of Cr in Fe below 700 K remains in the range of about 6 to 12%. It reproduces the transition between the ordering and demixing tendency and the spinodal decomposition limits are also in agreement with the values given in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Structure and dynamics in yttrium-based molten rare earth alkali fluorides

The transport properties of molten LiF-YF3 mixtures have been studied by pulsed field gradient nuclear magnetic resonance spectroscopy, potentiometric experiments, and molecular dynamics simulations. The calculated diffusion coefficients and electric conductivities compare very well with the measurements accross a wide composition range. We then extract static (radial distribution functions, coordination numbers distributions) and dynamic (cage correlation functions) quantities from the simulations. Then, we discuss the interplay between the microscopic structure of the molten salts and their dynamic properties. It is often considered that variations in the diffusion coefficient of the anions are mainly driven by the evolution of its coordination with the metallic ion (Y3+ here). We compare this system with fluorozirconate melts and demonstrate that the coordination number is a poor indicator of the evolution of the diffusion coefficient. Instead, we propose to use the ionic bonds lifetime. We show that the weak Y-F ionic bonds in LiF-YF3 do not induce the expected tendency of the fluoride diffusion coefficient to converge toward the one of yttrium cation when the content in YF3 increases. Implications on the validity of the Nernst-Einstein relation for estimating the electrical conductivity are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Molecular diffusion between walls with adsorption and desorption

The time dependency of the diffusion coefficient of particles in porous media is an efficient probe of their geometry. The analysis of this quantity, measured e.g. by nuclear magnetic resonance (PGSE-NMR), can provide rich information pertaining to porosity, pore size distribution, permeability and surface-to-volume ratio of porous materials. Nevertheless, in numerous if not all practical situations, transport is confined by walls where adsorption and desorption processes may occur. In this article, we derive explicitly the expression of the time-dependent diffusion coefficient between two confining walls in the presence of adsorption and desorption. We show that they strongly modify the time-dependency of the diffusion coefficient, even in this simple geometry. We finally propose several applications, from sorption rates measurements to the use as a reference for numerical implementations for more complex geometries.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Taylor Dispersion with Adsorption and Desorption

We use a stochastic approach to show how Taylor dispersion is affected by kinetic processes of adsorption and desorption onto surfaces. A general theory is developed, from which we derive explicitly the dispersion coefficients of canonical examples like Poiseuille flows in planar and cylindrical geometries, both in constant and sinusoidal velocity fields. These results open the way for the measurement of adsorption and desorption rate constants using stationary flows and molecular sorting using the stochastic resonance of the adsorption and desorption processes with the oscillatory velocity field.

0 citations0 reviewsNo code linkedOpen access