Source author record

Benjamin Rotenberg

Benjamin Rotenberg 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

19works
6topics
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

19 published item(s)

preprint2022arXiv

Chemi-sorbed versus physi-sorbed surface charge and its impact on electrokinetic transport: carbon versus boron-nitride surface

The possibility to control electrokinetic transport through carbon and hexagonal boron nitride (hBN) nanotubes has recently opened new avenues for nanofluidic approaches to face outstanding challenges such as energy production and conversion or water desalination. The $p$H-dependence of experimental transport coefficients point to the sorption of hydroxide ions as the microscopic origin of the surface charge and recent ab initio calculations suggest that these ions behave differently on carbon and hBN, with only physisorption on the former and both physi- and chemisorption on the latter. Using classical non-equilibrium molecular dynamics simulations of interfaces between an aqueous electrolyte and three models of hBN and graphite surfaces, we demonstrate the major influence of the sorption mode of hydroxide ions on the interfacial transport properties. Physisorbed surface charge leads to a considerable enhancement of the surface conductivity as compared to its chemisorbed counterpart, while values of the $ζ$-potential are less affected . The analysis of the MD results for the surface conductivity and $ζ$-potential in the framework of Poisson-Boltzmann-Stokes theory, as is usually done to analyze experimental data, further confirms the importance of taking into account both the mobility of surface hydroxide ions and the decrease of the slip length with increasing titratable surface charge density.

preprint2022arXiv

Likelihood-based non-Markovian models from molecular dynamics

We introduce a new method to accurately and efficiently estimate the effective dynamics of collective variables in molecular simulations. Such reduced dynamics play an essential role in the study of a broad class of processes, ranging from chemical reactions in solution to conformational changes in biomolecules or phase transitions in condensed matter systems. The standard Markovian approximation often breaks down due to the lack of a proper separation of time scales and memory effects must be taken into account. Using a parametrization based on hidden auxiliary variables, we obtain a generalized Langevin equation by maximizing the statistical likelihood of the observed trajectories. Both the memory kernel and random noise are correctly recovered by this procedure. This data-driven approach provides a reduced dynamical model for multidimensional collective variables, enabling the accurate sampling of their long-time dynamical properties at a computational cost drastically reduced with respect to all-atom numerical simulations. The present strategy, based on the reproduction of the dynamics of trajectories rather than the memory kernel or the velocity-autocorrelation function, conveniently provides other observables beyond these two, including e.g. stationary currents in non-equilibrium situations, or the distribution of first passage times between metastable states.

preprint2020arXiv

Field-dependent ionic conductivities from generalized fluctuation-dissipation relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian and subsequently the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian and the conductivity increases with applied field. This nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady-states.

preprint2020arXiv

Sampling mobility profiles of confined fluids with equilibrium molecular dynamics simulations

We show how to evaluate mobility profiles, characterizing the transport of confined fluids under a perturbation, from equilibrium molecular simulations. The correlation functions derived with the Green-Kubo formalism are difficult to sample accurately and we consider two complementary strategies: improving the spacial sampling thanks to a new estimator of the local fluxes involving the forces acting on the particles in addition to their positions and velocities, and improving temporal sampling thanks to the Einstein-Helfand approach instead of the Green-Kubo one. We illustrate this method on the case of a binary mixture confined between parallel walls, under a pressure or chemical potential gradient. All equilibrium methods are compared to standard non-equilibrium molecular dynamics (NEMD) and provide the correct mobility profiles. We recover quantitatively fluid viscosity and diffusio-osmostic mobility in the bulk part of the pore. Interestingly, the matrix of mobility profiles for local fluxes is not symmetric, unlike the Onsager matrix for the total fluxes. Even the most computationally efficient equilibrium method (Einstein-Helfand combined with the force-based estimator) remains less efficient than NEMD to determine a single mobility profile. However, the equilibrium approach provides all responses to all perturbations simultaneously, whereas NEMD requires the simulation of several types of perturbations to determine the various responses, each with different magnitudes to check the validity of the linear regime. While NEMD seems more competitive for the present example, the balance should be different for more complex systems, in particular for electrolyte solutions for the responses to pressure, salt concentration and electric potential gradients.

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.

preprint2016arXiv

Molecular hydrodynamics from memory kernels

The memory kernel for a tagged particle in a fluid, computed from molecular dynamics simulations, decays algebraically as $t^{-3/2}$. We show how the hydrodynamic Basset-Boussinesq force naturally emerges from this long-time tail and generalize the concept of hydrodynamic added mass. This mass term is negative in the present case of a molecular solute, at odds with incompressible hydrodynamics predictions. We finally discuss the various contributions to the friction, the associated time scales and the cross-over between the molecular and hydrodynamic regimes upon increasing the solute radius.

preprint2016arXiv

Structural Transitions at Ionic Liquids Interfaces

Recent advances in experimental and computational techniques have allowed for an accurate description of the adsorption of ionic liquids on metallic electrodes. It is now well established that they adopt a multi-layered structure, and that the composition of the layers changes with the potential of the electrode. In some cases, potential-driven ordering transitions in the first adsorbed layer have been observed in experiments probing the interface on the molecular scale or by molecular simulations. This perspective gives an overview of the current understanding of such transitions and of their potential impact on the physical and (electro)chemical processes at the interface. In particular, peaks in the differential capacitance, slow dynamics at the interface and changes in the reactivity have been reported in electrochemical studies. Interfaces between ionic liquids and metallic electrodes are also highly relevant for their friction properties, the voltage-dependence of which opens the way to exciting applications.

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.

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

The electric double layer has a life of its own

Using molecular dynamics simulations with recently developed importance sampling methods, we show that the differential capacitance of a model ionic liquid based double-layer capacitor exhibits an anomalous dependence on the applied electrical potential. Such behavior is qualitatively incompatible with standard mean-field theories of the electrical double layer, but is consistent with observations made in experiment. The anomalous response results from structural changes induced in the interfacial region of the ionic liquid as it develops a charge density to screen the charge induced on the electrode surface. These structural changes are strongly influenced by the out-of-plane layering of the electrolyte and are multifaceted, including an abrupt local ordering of the ions adsorbed in the plane of the electrode surface, reorientation of molecular ions, and the spontaneous exchange of ions between different layers of the electrolyte close to the electrode surface. The local ordering exhibits signatures of a first-order phase transition, which would indicate a singular charge-density transition in a macroscopic limit.

preprint2013arXiv

Charge fluctuations in nano-scale capacitors

The fluctuations of the charge on an electrode contain information on the microscopic correlations within the adjacent fluid and their effect on the electronic properties of the interface. We investigate these fluctuations using molecular dynamics simulations in a constant-potential ensemble with histogram reweighting techniques. This approach offers in particular an efficient, accurate and physically insightful route to the differential capacitance that is broadly applicable. We demonstrate these methods with three different capacitors: pure water between platinum electrodes, and a pure as well as a solvent-based organic electrolyte each between graphite electrodes. The total charge distributions with the pure solvent and solvent-based electrolytes are remarkably Gaussian, while in the pure ionic liquid the total charge distribution displays distinct non-Gaussian features, suggesting significant potential-driven changes in the organization of the interfacial fluid.

preprint2013arXiv

Influence of solvation on the structural and capacitive properties of electrical double layer capacitors

We use molecular dynamics simulations to explore the impact of a non-ionic solvent on the structural and capacitive properties of supercapacitors based on an ionic liquid electrolyte and carbon electrodes. The study is focused on two pure ionic liquids, namely 1-butyl-3-methylimidazolium hexafluorophosphate and 1-butyl-3-methylimidazolium tetrafluoroborate, and their 1.5 M solutions in acetonitrile. The electrolytes, represented by coarse-grained models, are enclosed between graphite electrodes. We employ a constant potential methodology which allows us to gain insight into the influence of solvation on the polarization of the electrodes as well as the structural and capacitive properties of the electrolytes at the interface. We show that the interfacial characteristics, different for two distinct pure ionic liquids, become very similar upon mixing with acetonitrile.

preprint2012arXiv

A transferable ab-initio based force field for aqueous ions

We present a new polarizable force field for aqueous ions (Li+, Na+, K+, Rb+, Cs+, Mg2+, Ca2+, Sr2+ and Cl-) derived from condensed phase ab-initio calculations. We use Maximally Localized Wannier Functions together with a generalized force and dipole-matching procedure to determine the whole set of parameters. Experimental data is then used only for validation purposes and a good agreement is obtained for structural, dynamic and thermodynamic properties. The same procedure applied to crystalline phases allows to parametrize the interaction between cations and the chloride anion. Finally, we illustrate the good transferability of the force field to other thermodynamic conditions by investigating concentrated solutions.

preprint2012arXiv

Diffusion coefficient and shear viscosity of rigid water models

We report the diffusion coefficient and viscosity of popular rigid water models: Two non polarizable ones (SPC/E with 3 sites, and TIP4P/2005 with 4 sites) and a polarizable one (Dang-Chang, 4 sites). We exploit the dependence of the diffusion coefficient on the system size [Yeh and Hummer, J. Phys. Chem. B 108, 15873 (2004)] to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang water overestimate the diffusion coefficient and underestimate the viscosity.

preprint2012arXiv

Including many-body effects in models for ionic liquids

Realistic modeling of ionic systems necessitates taking explicitly account of many-body effects. In molecular dynamics simulations, it is possible to introduce explicitly these effects through the use of additional degrees of freedom. Here we present two models: The first one only includes dipole polarization effect, while the second also accounts for quadrupole polarization as well as the effects of compression and deformation of an ion by its immediate coordination environment. All the parameters involved in these models are extracted from first-principles density functional theory calculations. This step is routinely done through an extended force-matching procedure, which has proven to be very succesfull for molten oxides and molten fluorides. Recent developments based on the use of localized orbitals can be used to complement the force-matching procedure by allowing for the direct calculations of several parameters such as the individual polarizabilities.

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.

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.

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.

preprint2011arXiv

Molecular explanation for why talc surfaces can be both hydrophilic and hydrophobic

While individual water molecules adsorb strongly on a talc surface (hydrophilic behavior), a droplet of water beads up on the same surface (hydrophobic behavior). To rationalize this dichotomy, we investigate the influence of the microscopic structure of the surface and the strength of adhesive (surface-water) interactions on surface hydrophobicity. We show that at low relative humidity, the competition between adhesion and the favorable entropy of being in the vapor phase determines the surface coverage. However, at saturation, it is the competition between adhesion and cohesion (water-water interactions) that determines surface hydrophobicity. The adhesive interactions in talc are strong enough to overcome the unfavorable entropy, and water adsorbs strongly on talc surfaces. However, they are too weak to overcome the cohesive interactions, and water thus beads up on talc surfaces. Surprisingly, even (talc-like) surfaces that are highly adhesive, do not fully wet at saturation. Instead, a water droplet forms on top of a strongly adsorbed monolayer of water. Our results imply that the interior of hydrophobic zeolites suspended in water may contain adsorbed water molecules at pressures much smaller than the intrusion pressure.