@article {175, title = {Solubility of Gases in Water Confined in Nanoporous Materials: ZSM-5, MCM-41, and MIL-100}, journal = {The Journal of Physical Chemistry C}, volume = {119}, year = {2015}, month = {Sep-17-2015}, pages = {21547 - 21554}, abstract = {
Abstract Image

Oversolubility, which corresponds to large apparent gas solubilities in liquids confined in nanoporous solids, has been proposed as a means to develop novel adsorption, phase separation, or catalytic processes. We report a molecular simulation study to help design such hybrid adsorbents consisting of a solvent confined in host nanoporous materials. Water is selected as the confined solvent because of its ubiquity in industrial applications, while N2, CH4, and CO2 are selected as they allow probing the effect of specific molecular interaction and polarity. For each system, we consider a zeolite, an ordered mesoporous silica, and a metal\–organic framework as the confining host as they present different morphologies, porosities, and surface chemistries. We show that oversolubility is always observed because the apparent solubility in these materials surpasses the bulk solubility. In all cases, oversolubility is an enhanced bulklike solubility in which solubility is favored in the regions of low water density formed by the layering of the solvent. Such an oversolubility mechanism, which arises from the fact that the gas\–solid interactions are weaker than the solvent\–solid interactions, leads to large uptakes as high as a few hundred times that expected from bulk solubility.

}, issn = {1932-7447}, doi = {10.1021/acs.jpcc.5b06660}, author = {Linh Ngoc Ho and Schuurman, Yves and David Farrusseng and Benoit A. Coasne} } @book {648, title = {Springer Series in Materials ScienceMolecular Dynamics Simulations of Disordered MaterialsMolecular Modeling of Glassy Surfaces}, volume = {215}, year = {2015}, pages = {345 - 365}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

Progress in computational materials science has allowed the development of realistic models for a wide range of materials including both crystalline and glassy solids. In recent years, with the growing interest in nanoparticles and porous materials, more attention has been devoted to the design of realistic models of glassy surfaces and finely divided materials. The structural disorder in glassy surfaces, however, poses a major challenge which consists of describing such surfaces using computer simulations. In this paper, we show how atomic-scale simulations can be used to develop and investigate the properties of glassy surfaces. We illustrate how both first principles calculations and classical molecular mechanics can be used to follow the trajectory at finite temperature of these systems, and obtain statistical thermodynamic averages to compare against available experiments. Both glassy oxide (silica) and non-oxide (chalcogenide) surfaces are considered.

}, isbn = {978-3-319-15674-3}, issn = {0933-033X}, doi = {10.1007/978-3-319-15675-010.1007/978-3-319-15675-0_13}, url = {http://link.springer.com/10.1007/978-3-319-15675-0}, author = {Guido Ori and Carlo Massobrio and Bouzid, A. and Benoit A. Coasne}, editor = {Carlo Massobrio and Du, Jincheng and Bernasconi, Marco and Salmon, Philip S.} } @article {180, title = {Structure and Dynamics of Ionic Liquids Confined in Amorphous Porous Chalcogenides}, journal = {Langmuir}, volume = {31}, year = {2015}, month = {Jun-23-2015}, pages = {6742 - 6751}, abstract = {
Abstract Image

Besides the abundant literature on ionic liquids in porous silica and carbon, the confinement of such intriguing liquids in porous chalcogenides has received very little attention. Here, molecular simulation is employed to study the structural and dynamical properties of a typical ionic liquid confined in a realistic molecular model of amorphous chalcogenide with various pore sizes and surface chemistries. Using molecular dynamics in the isobaric\–isothermal (NPT) ensemble, we consider confinement conditions relevant to real samples. Both the structure and self-dynamics of the confined phase are found to depend on the surface-to-volume ratio of the host confining material. Consequently, most properties of the confined ionic liquid can be written as a linear combination of surface and bulk-like contributions, arising from the ions in contact with the surface and the ions in the pore center, respectively. On the other hand, collective dynamical properties such as the ionic conductivity remain close to their bulk counterpart and almost insensitive to pore size and surface chemistry. These results, which are in fair agreement with available experimental data, provide a basis for the development of novel applications using hybrid organic\–inorganic solids consisting of ionic liquids confined in porous chalcogenides.

}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.5b00982}, author = {Guido Ori and Carlo Massobrio and Pradel, Annie and Ribes, Michel and Benoit A. Coasne} } @article {132, title = {Subcontinuum mass transport of condensed hydrocarbons in nanoporous media.}, journal = {Nat Commun}, volume = {6}, year = {2015}, month = {Apr-22-2015}, pages = {Article Number: 6949}, abstract = {

Although hydrocarbon production from unconventional reservoirs, the so-called shale gas, has exploded recently, reliable predictions of resource availability and extraction are missing because conventional tools fail to account for their ultra-low permeability and complexity. Here, we use molecular simulation and statistical mechanics to show that continuum description--Darcy\&$\#$39;s law--fails to predict transport in shales nanoporous matrix (kerogen). The non-Darcy behaviour arises from strong adsorption in kerogen and the breakdown of hydrodynamics at the nanoscale, which contradict the assumption of viscous flow. Despite this complexity, all permeances collapse on a master curve with an unexpected dependence on alkane length. We rationalize this non-hydrodynamic behaviour using a molecular description capturing the scaling of permeance with alkane length and density. These results, which stress the need for a change of paradigm from classical descriptions to nanofluidic transport, have implications for shale gas but more generally for transport in nanoporous media.

Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

}, issn = {2041-1723}, doi = {10.1038/ncomms7949}, author = {Kerstyn I. Falk and Benoit A. Coasne and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Lyd{\'e}ric Bocquet} } @article {63, title = {Structure and Dynamics of an Electrolyte Confined in Charged Nanopores}, journal = {The Journal of Physical Chemistry C}, volume = {118}, year = {2014}, month = {Feb-2014}, pages = {5061{\textendash}5072}, abstract = {
Abstract Image

Molecular Dynamics simulations are used to investigate the structure and dynamics of an aqueous electrolyte (NaCl) confined within a nanomembrane, which consists of a nanopore with a diameter 3 nm having a negatively charged surface. Both nanomembranes with a diffuse charge and with local charges are considered (in both cases, two surface charge densities are considered, \−0.9 e/nm2 and \−1.8 e/nm2). For all nanomembranes, significant layering of water and ions in the vicinity of the nanomembrane surface is observed. While the distribution of water and chloride ions is nearly insensitive to the nanomembrane charge and type, the arrangement of sodium cations within the nanomembrane depends on the system being considered. The water and ion density profiles in the nanomembranes are compared with the predictions of a modified Poisson\–Boltzmann equation in which charge image, solvation effects, and dispersion interactions with the surface are taken into account [Huang et al. Langmuir, 2008, 24, 1442]. The self-diffusion coefficient for a given species is smaller than its bulk counterpart and is at most 75\% of the bulk value. While the self-diffusion coefficients for water and sodium cations decrease with decreasing the overall negative charge of the nanomembrane, the self-diffusion coefficient for the chloride anions is nearly independent of the nanomembrane type and charge. We also estimate the dynamics of the confined aqueous electrolyte by calculating time correlation functions which allow estimating solvation, ion pairing, and residence times.

}, isbn = {1932-7447}, doi = {10.1021/jp4098638}, author = {Pierre-Andre Cazade and Hartkamp, Remco and Benoit A. Coasne} } @article {90, title = {Structure and transport of aqueous electrolytes: From simple halides to radionuclide ions}, journal = {Journal of Chemical Physics}, volume = {141}, year = {2014}, month = {Sept-28-2014}, pages = {124508}, type = {Article}, abstract = {

Molecular simulations are used to compare the structure and dynamics of conventional and radioactive aqueous electrolytes: chloride solutions with sodium, potassium, cesium, calcium, and strontium. The study of Cs+ and Sr2+ is important because these radioactive ions can be extremely harmful and are often confused by living organisms for K+ and Ca2+, respectively. Na+, Ca2+, and Sr2+ are strongly bonded to their hydration shell because of their large charge density. We find that the water molecules in the first hydration shell around Na+ form hydrogen bonds between each other, whereas molecules in the first hydration shell around Ca2+ and Sr2+ predominantly form hydrogen bonds with water molecules in the second shell. In contrast to these three ions, K+ and Cs+ have low charge densities so that they are weakly bonded to their hydration shell. Overall, the structural differences between Ca2+ and Sr2+ are small, but the difference between their coordination numbers relative to their surface areas could potentially be used to separate these ions. Moreover, the different decays of the velocity-autocorrelation functions corresponding to these ions indicates that the difference in mass could be used to separate these cations. In this work, we also propose a new definition of the pairing time that is easy to calculate and of physical significance regardless of the problem at hand. (C) 2014 AIP Publishing LLC.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
Image of FIG. 4.
Image of FIG. 5.
Image of FIG. 6.
}, issn = {0021-9606}, doi = {10.1063/1.4896380}, author = {Hartkamp, Remco and Benoit A. Coasne} } @article {79, title = {Structure-property relationships of water adsorption in metal-organic frameworks}, journal = {New Journal of Chemistry}, volume = {38}, year = {2014}, month = {Apr-16-2014}, pages = {3102-3111}, type = {Article}, abstract = {

A set of 15 metal\–organic frameworks (MIL-53, MIL-68, MIL-125, UiO-66, ZIF) exhibiting different pore size, morphology, and surface chemistry is used to unravel the numerous behaviors of water adsorption at room temperature in this class of materials. Outstanding \“S\”-shaped (type V) adsorption isotherms are observed for MIL-68 type solids. We show that the underlying mechanism of water adsorption can be rationalized using a simple set of three parameters: the Henry constant (i.e. the slope of the adsorption pressure in the low pressure range), the pressure at which pore filling occurs, and the maximum water adsorption capacity. While the Henry constant and pore filling pressure mostly depend on the affinity of water for the surface chemistry and on pore size, respectively, these two parameters are correlated as they both reflect different aspects of the hydrophobicity\–hydrophilicity of the material. For a given type of porous structure, the functionalization of the material by hydrophilic moieties such as hydrogen bonding groups (amine or aldehyde) systematically leads to an increase in the Henry constant concomitantly with a decrease in the pore filling pressure. As for the adsorption mechanism, we show that, for a given temperature, there is a critical diameter (Dc \∼ 20 \Å for water at room temperature) above which pore filling occurs through irreversible capillary condensation accompanied by capillary hysteresis loops. Below this critical diameter, pore filling is continuous and reversible unless the material exhibits some adsorption-induced flexibility.

Graphical abstract: Structure{\textendash}property relationships of water adsorption in metal{\textendash}organic frameworks
}, issn = {1144-0546}, doi = {10.1039/c4nj00076e}, author = {Canivet, Jerome and Bonnefoy, Jonathan and Daniel, Cecile and Legrand, Alexandre and Benoit A. Coasne and David Farrusseng} } @article {86, title = {Surface of glassy GeS2: A model based on a first-principles approach}, journal = {Physical Review B}, volume = {90}, year = {2014}, month = {Jul-2014}, type = {Article}, chapter = {045423}, abstract = {

First-principles calculations within the framework of the density functional theory are used to construct realistic models for the surface of glassy GeS2(g\−GeS2). Both calculations at T=0 K and at finite temperature (T=300 K) are considered. This allows for a comparison between the structural and electronic properties of surface and bulk g\−GeS2. Although the g\−GeS2 surface recovers the main tetrahedral structural motif of bulk g\−GeS2, the number of fourfold coordinated Ge atoms and twofold coordinated S atoms is smaller than in the bulk. On the contrary, the surface system features a larger content of overcoordinated S atoms and threefold coordinated Ge atoms. This effect is more important for the g\−GeS2 surface relaxed at 0 K. Maximally localized Wannier functions (WF) are used to inspect the nature of the chemical bonds of the structural units present at the g\−GeS2 surface. We compare the ability of several charge derivation methods to capture the atomic charge variations induced by a coordination change. Our estimate for the charges allows exploiting the first-principles results as a data base to construct a reliable interatomic force field.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.90.045423}, author = {Guido Ori and Carlo Massobrio and Bouzid, A. and Boero, M. and Benoit A. Coasne} } @article {82, title = {Surfactant Behavior of Ionic Liquids Involving a Drug: From Molecular Interactions to Self-Assembly}, journal = {Langmuir}, volume = {30}, year = {2014}, month = {Feb-2014}, pages = {1229-1238}, type = {Article}, abstract = {
Abstract Image

Aggregates formed in an aqueous medium by three ionic liquids CnMImIbu made up of 1-alkyl-3-methyl-imidazolium cation (n = 4, 6, 8) and ibuprofenate anion are investigated. Dynamic light scattering (DLS), cryogenic transmission electron microscopy (cryo-TEM), 1H nuclear magnetic resonance measurements, and atom-scale molecular dynamics simulations are used to shed light on the main interactions governing the formation of the aggregates and their composition. At high concentration, mixed micelles are formed with a composition that depends on the imidazolium alkyl chain length. For the shortest alkyl chain, micelles are mainly composed of ibuprofenate anions with some imidazolium cations intercalated between the anions. Upon increasing the alkyl chain length, the composition of the aggregates gets enriched in imidazolium cations and aggregates of stoichiometric composition are obtained. Attractive interactions between these aggregates led to the formation of larger aggregates. As suggested by molecular simulations, these larger aggregates might constitute the early stage of phase separation. Transitions from micelles to vesicles or ribbons are observed due to dilution effects and changes in the chemical composition of the aggregates. We also show that aggregation can be probed using simple microscopic quantities such as radial distribution functions and average solvation numbers.

}, issn = {0743-7463}, doi = {10.1021/la404166y}, author = {Tourne-Peteilh, Corine and Benoit A. Coasne and In, Martin and Brevet, David and Devoisselle, Jean-Marie and Andre Vioux and Lydie Viau} } @article {49, title = {Solvated calcium ions in charged silica nanopores}, journal = {The Journal of Chemical Physics}, volume = {137}, year = {2012}, month = {Aug-14-2012}, pages = {Article Number: 064706}, abstract = {

Hydroxyl surface density in porous silica drops down to nearly zero when the pH of the confined aqueous solution is greater than 10.5. To study such extreme conditions, we developed a model of slit silicananopores where all the hydrogen atoms of the hydroxylated surface are removed and the negative charge of the resulting oxygen dangling bonds is compensated by Ca2+ counterions. We employed grand canonical Monte Carlo and molecular dynamics simulations to address how the Ca2+ counterions affect the thermodynamics, structure, and dynamics of confined water. While most of the Ca2+ counterions arrange themselves according to the so-called \“Stern layer,\” no diffuse layer is observed. The presence of Ca2+ counterions affects the pore filling for strong confinement where the surfaceeffects are large. At full loading, no significant changes are observed in the layering of the first two adsorbed water layers compared to nanopores with fully hydroxylated surfaces. However, the water structure and water orientational ordering with respect to the surface is much more disturbed. Due to the super hydrophilicity of the Ca2+-silica nanopores,waterdynamics is slowed down and vicinal water molecules stick to the pore surface over longer times than in the case of hydroxylated silicasurfaces. These findings, which suggest the breakdown of the linear Poisson\–Boltzmann theory, provide important information about the properties of nanoconfined electrolytes upon extreme conditions where the surface charge and ion concentration are large.

Image of FIG. 1.

Image of FIG. 2.

Image of FIG. 3.

Image of FIG. 4.

Image of FIG. 5.

Image of FIG. 6.

}, doi = {10.1063/1.4742854}, author = {Patrick A Bonnaud and Benoit A. Coasne and Roland Jean-Marc Pellenq} }