@article {659, title = {Bulk Photodriven CO2 Conversion through TiO2@Si(HIPE) Monolithic Macrocellular Foams}, journal = {Advanced Functional Materials}, volume = {3}, year = {2019}, month = {Jan-07-2019}, pages = {1807767}, abstract = {

Operating photo-induced reactions exclusively on catalyst surfaces while not exploiting the full catalyst volume generates a major footprint penalty for the photocatalytic reactor and leads to an inefficient use of the catalytic material. Photonic investigations clearly show that the solid foams have a strongly multidiffusive character, with photons being significantly trapped within the sample cores while addressing a photon mean free path lt = 20.1 \± 1.3 \µm. This 3D process both greatly limits back-reactions and promotes outstanding selectivity toward methane (around 80\%) generation, and even ethane (around 18\%) through C-C coupling reaction, with residual carbon monoxide and dihydrogen contents (around 2\%). Silica\–titania TiO2@Si(HIPE) self-standing macrocellular catalysts lead to optimal efficient thicknesses up to 20 times those of powders, thereby enhancing the way for real 3D-photodriven catalytic processes above the millimeter scale and up to a 6 mm thickness. A rather simple Langmuir\–Hinshelwood based kinetic model is proposed which highlights the strong dependence of photocatalytic reaction rates on light scattering and the crucial role on oxidation back-reactions. In addition, a strong correlation between light attenuation coefficient and photon mean free path and median pore aperture diameter is demonstrated, offering thus a tool for photocatalytic behavior prediction.

}, doi = {10.1002/adfm.201807767}, url = {http://doi.wiley.com/10.1002/adfm.201807767}, author = {Bernadet, Sophie and Tavernier, Eug{\'e}nie and Ta, Duc-Minh and Vall{\'e}e, Renaud A. L. and Ravaine, Serge and F{\'e}cant, Antoine and Backov, Renal} } @article {594, title = {Investigation of mixed ionic/nonionic building blocks for the dual templating of macro-mesoporous silica}, journal = {Journal of Colloid and Interface Science}, volume = {533}, year = {2019}, month = {Jan-01-2019}, pages = {385 - 400}, abstract = {

Traditional porous monoliths Si(HIPE) (High Internal Phase Emulsion), prepared from the Tetradecyltrimethylammonium Bromide (TTAB)/dodecane/water system, offer high specific surface area, mainly due to microporosity. Aside, mesoporous materials SBA-15, prepared from Pluronic P123, have a high specific surface area, but are obtained as powder, which limits their applications. Starting from the mixed TTAB-P123 surfactant, it is expected to tune the mesoporosity of Si(HIPE), while keeping their monolithic character.

The ternary TTAB/P123/water phase diagram was established by varying the weight ratio between these two surfactants. The micellar structure as well as the structural parameters of the liquid crystal domains were determined by SAXS (Small Angle X-ray Scattering). The effect of dodecane solubilization was also investigated and concentrated emulsions were formulated from the (P123/TTAB)/dodecane/ water systems. After this soft matter dedicated study, the acquired knowledge was transferred toward the hierarchical porous silica generations, where the sol-gel process is involved.

Mixing P123 with TTAB, macro-mesoporous monolithic silica with an enhanced contribution of the specific surface area due to mesoporosity can be prepared. The variation of the TTAB/P123 weight ratio allows controlling the porosity at the mesoscale. Moreover, the macroporosity can be tuned by changing the preparation method, by mixing either the two micellar solutions or directly the two surfactants prior the emulsification process. (C) 2018 Elsevier Inc. All rights reserved.

}, issn = {00219797}, doi = {10.1016/j.jcis.2018.08.068}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0021979718309883}, author = {Roucher, Armand and Emo, Melanie and Vibert, Francois and Stebe, Marie-Jose and Schmitt, Veronique and Jonas, Florian and Backov, Renal and Blin, Jean-Luc} } @article {523, title = {Concrete material science: Past, present, and future innovations}, journal = {Cement and Concrete Research}, volume = {Volume 112}, year = {2018}, month = {Oct-01-2018}, pages = {5-24}, abstract = {

Concrete is flying off, but it is simultaneously facing tremendous challenges in terms of environmental impact, financial needs, societal acceptance and image. Based on an historical approach of the science of concrete and reinforced concrete in particular, this paper calls for the exploration of radical changes in three key aspects of concrete use: reinforcement, binder content, and implementation methods. More precisely, it is suggested that, in parallel to the introduction of robotic fabrication methods, digital technologies may be key for the introduction several innovations like (i) rebar-free reinforcement using non-convex granular media; (ii) compression-optimized concrete structures, using topology optimization, architectural geometry, and 3D-printing or origami-patterned formworks; (iii) truly digital concrete through the coupling of massive data collection and deep learning.

CCR DIGITAL CONCRETE.pdf

Fig. 1. Top: Comparative evolution of the post-WWII global cement, steel, and plastic{\textellipsis}

Fig. 2. Top: The four beams on this figure have the same load bearing capacity

Fig. 3. (a): A {\textquotedblleft}deterministic{\textquotedblright} or regular Apollonian packing of disks in a circular{\textellipsis}

Fig. 4. A surface area vs particle size plot for the different granular populations -{\textellipsis}

Fig. 5. Minimum porosity of 17 different concretes, plotted vs the ratio between the{\textellipsis}

Fig. 6. Schematic illustration of the difference between an Apollonian packing (a) and{\textellipsis}

Fig. 7. Two perturbations occurring in a dense assembly of particles with at least two{\textellipsis}

Fig. 8. SEM images of cement pastes hydrated without (a and c) or with (b and d){\textellipsis}

Fig. 9. Binder intensity index as a function of 28days strength for close to one{\textellipsis}

Fig. 10. Entangled assemblages of non-convex particles which could provide{\textellipsis}

Fig. 11. Free-form architecture with UHPFRC concrete: (A) Curved panels of the Louis{\textellipsis}

Fig. 12. (A) Aircraft hangar, Orvieto, architect-engineer Pier Luigi Nervi, 1935; (B) A{\textellipsis}

Fig. 13. Cartoon depicting a paradigm for transformative evolution of research on{\textellipsis}

}, issn = {00088846}, doi = {10.1016/j.cemconres.2018.05.002}, url = {https://www.sciencedirect.com/science/article/pii/S0008884618300802}, author = {Henri Van Damme} } @article {620, title = {Earth concrete. Stabilization revisited}, journal = {Cement and Concrete Research}, volume = {114}, year = {2018}, month = {Dec-2018}, pages = {90 - 102}, abstract = {

Not surprisingly, with the increased awareness of environmental issues, construction with raw (crude, unbaked) earth (subsoil) is gaining renewed interest. However, it suffers from a poor image and from the difficulty to meet modern productivity standards and to pass some durability tests designed for industrial materials. The recent trend is to overcome these drawbacks by \“stabilizing\” the material most often with Portland cement (PC). Here we show that stabilization with PC is in general neither technically nor environmentally advisable. It brings only moderate mechanical improvement at a high environmental cost. Rather than massively transforming crude earth into a low quality concrete, it would be more appropriate to adapt the architectural practice and/or to look for milder ways to improve properties. In this respect, the recent successful attempts to improve the workability and the strength of raw earth by controlling the dispersion of its fine fraction seem to be particularly promising.

}, issn = {00088846}, doi = {10.1016/j.cemconres.2017.02.035}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008884616308365}, author = {Henri Van Damme and Houben, Hugo} } @article {533, title = {Mechanical strength of wet particle agglomerates}, journal = {Mechanics Research Communications}, volume = {92}, year = {2018}, month = {Sep-2018}, pages = {1 - 7}, abstract = {

Using particle dynamics simulations, we investigate the strength and microstructure of agglomerates of wet frictional particles subjected to axial compression. The numerical model accounts for the cohesive and viscous effects of the binding liquid up to a debonding distance with the liquid assumed to be distributed homogeneously inside the agglomerate. We show that wet agglomerates undergo plastic deformation due to the rearrangements of primary particles during compression. The compressive strength is thus characterized by the plastic threshold before the onset of failure by the irreversible loss of wet contacts between primary particles. We find that the agglomerate plastic threshold is proportional to the characteristic cohesive stress defined from the liquid-vapor surface tension and the mean diameter of primary particles, with a prefactor that is a nearly linear function of the debonding distance and increases with size span. We analyze the agglomerate microstructure and, considering only the cohesive capillary forces at all bonds between primary particles, we propose an expression of the plastic strength as a function of the texture parameters such as the wet coordination number and packing fraction. This expression is shown to be consistent with our simulations up to a multiplicative factor reflecting the distribution of the capillary bridges. (C) 2018 Published by Elsevier Ltd.

Fig. 1. Schematic drawing of the forces acting on particle i by a contacting particle j{\textellipsis}

Fig. 2. Schematic representation of diametrical compression test (a), and snapshots of{\textellipsis}

Fig. 3. Snapshot of a granule under diametrical compression for α=1

Fig. 4. Vertical strength σzz normalized by the reference stress σc for different{\textellipsis}

Fig. 5. Normalized plastic strength of the granule for several values of the debonding{\textellipsis}

Fig. 6. Normalized plastic strength σp/σc of the granules for different values of the{\textellipsis}

Fig. 7. The initial wet coordination number Z0 of the granules for different values of{\textellipsis}

Fig. 8. The prefactor η in Eq

Fig. 9. The ZΦs in Eq

}, issn = {00936413}, doi = {10.1016/j.mechrescom.2018.07.003}, url = {https://www-sciencedirect-com.libproxy.mit.edu/science/article/pii/S0093641318301216}, author = {Vo, Thanh-Trung and Patrick Mutabaruka and Saeid Nezamabadi and Jean-Yves Delenne and Izard, Edouard and Roland Jean-Marc Pellenq and Farhang Radja{\"\i}} } @article {630, title = {Mesoscale structure, mechanics, and transport properties of source rocks{\textquoteright} organic pore networks}, journal = {Proceedings of the National Academy of Sciences}, volume = {115}, year = {2018}, month = {Dec-04-2018}, pages = {12365-12370}, abstract = {

Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores\&$\#$39; topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.

}, keywords = {porous media; electron tomography; mechanics; fluid transport; mesoscale}, issn = {0027-8424}, doi = {10.1073/pnas.1808402115}, url = {http://www.pnas.org/lookup/doi/10.1073/pnas.1808402115}, author = {Berthonneau, Jeremie and Ama{\"e}l Obliger and Valdenaire, Pierre-Louis and Grauby, Olivier and Ferry, Daniel and Chaudanson, Damien and Pierre E. Levitz and Kim, Jae Jin and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {570, title = {Poroelasticity of Methane-Loaded Mature and Immature Kerogen from Molecular Simulations}, journal = {Langmuir}, volume = {34}, year = {2018}, month = {Nov-13-2018}, pages = {13766 - 13780}, abstract = {

While hydrocarbon expulsion from kerogen is certainly the key step in shale oil/gas recovery, the poromechanical couplings governing this desorption process, taking place under a significant pressure gradient, are still poorly understood. Especially, most molecular simulation investigations of hydrocarbon adsorption and transport in kerogen have so far been performed under the rigid matrix approximation, implying that the pore space is independent of pressure, temperature, and fluid loading, or in other words, neglecting poromechanics. Here, using two hydrogenated porous carbon models as proxies for immature and overmature kerogen, that is, highly aliphatic hydrogen-rich vs highly aromatic hydrogen-poor models, we perform an extensive molecular-dynamics-based investigation of the evolution of the poroelastic properties of those matrices with respect to temperature, external pressure, and methane loading as a prototype alkane molecule. The rigid matrix approximation is shown to hold reasonably well for overmature kerogen even though accounting for flexibility has allowed us to observe the well-known small volume contraction at low fluid loading and temperature. Our results demonstrate that immature kerogen is highly deformable. Within the ranges of conditions considered in this work, its density can double and its accessible porosity (to a methane molecule) can increase from 0 to \∼30\%. We also show that these deformations are significantly nonaffine (i.e., nonhomogeneous), especially upon fluid adsorption or desorption.

}, keywords = {ADSORPTION-INDUCED DEFORMATION; ORGANIC-MATTER; COMPETITIVE ADSORPTION; NANOPOROUS MEDIA; POROUS CARBONS; GLASSY-CARBON; COAL; TRANSPORT; SORPTION; MODEL}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.8b02534}, url = {http://pubs.acs.org/doi/10.1021/acs.langmuir.8b02534}, author = {Obliger, Ama{\"e}l and Valdenaire, Pierre-Louis and Capit, Nicolas and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Jean-Marc Leyssale} } @article {525, title = {A reaction model for cement solidification: Evolving the C{\textendash}S{\textendash}H packing density at the micrometer-scale}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {118}, year = {2018}, month = {Sep-2018}, pages = {58 - 73}, abstract = {

Cement paste is a multiphase material of complex chemistry, of which 60\% by volume is typically composed of calcium\–silicate\–hydrates (C\–S\–H), the phase that lends the material its strength and stiffness. Moreover, it has been shown that the C\–S\–H phase is a dispersion of nanometer-sized particles that are characterized by an attractive\–repulsive potential and densify in course of the hydration reaction. Herein, we model the nucleation and growth of the nanoparticles as a continuous density field subject to a reaction equation. Using this phase-field approach, we aim to reduce the parameter space present in similar hydration models and create a vehicle to upscale mechanical information from the nanometer-scale to the micrometer-scale. Despite the apparent simplification of the physics at play, we readily reproduce the cement paste reaction kinetics and microtexture\—functions of temperature, coarseness of the calcium\–silicate source particles, and initial water-to-cement ratio\—, and vet them against experimental observations. Presenting results for two-dimensional simulations, we achieve excellent agreement with measurements of hydration heat curves, pore-chord-length and solid-chord-length density functions, distributions of low- and high-density C\–S\–H products, and elasticity.

Illustration of the dissolution and nucleation and growth processes driving{\textellipsis}

A comparison between the mean coordination number ⟨z⟩ of the particles modeled{\textellipsis}

Gibb{\textquoteright}s free energy that interpolates between the source phase ψ=1 and open pore{\textellipsis}

Gibb{\textquoteright}s energy of reaction and (b) reaction rate for nanoparticle{\textellipsis}

Evolution of the mean packing fraction of C{\textendash}S{\textendash}H (〈ϕ〉V'; opaque) and the{\textellipsis}

Diagram of the volume fractions of hydrating cement paste

}, issn = {00225096}, doi = {10.1016/j.jmps.2018.05.010}, url = {https://www-sciencedirect-com.libproxy.mit.edu/science/article/pii/S0022509617308876}, author = {Petersen, Thomas and Valdenaire, Pierre-Louis and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {600, title = {Cohesive strength of iron ore granules}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 08020}, abstract = {

We present an experimental and numerical investigation of the mechanical strength of crude iron ore (Hematite) granules in which capillary bonds between primary particles are the source of internal cohesion. The strength is measured by subjecting the granules to vertical compression between two plates. We show that the behavior of the granules is ductile with a well-defined plastic threshold which increases with the amount of water. It is found that the compressive strength scales with capillary cohesion with a pre-factor that is nearly independent of size polydispersity for the investigated range of parameters but increases with friction coefficient between primary particles. This weak dependence may be attributed to the class of fine particles which, due to their large number, behaves as a cohesive matrix that controls the strength of the granule.

}, doi = {10.1051/epjconf/201714008020}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714008020}, author = {Contreras, Rafael Jaimes and Berger, Nicolas and Izard, Edouard and Douce, Jean-Fran{\c c}ois and Koltsov, Alexey and Jean-Yves Delenne and Emilien Az{\'e}ma and Saeid Nezamabadi and van Loo, Fr{\'e}d{\'e}ric and Roland Jean-Marc Pellenq and Farhang Radja{\"\i}}, editor = {Luding, S.} } @article {588, title = {From cellulose to kerogen: molecular simulation of a geological process}, journal = {Chemical Science}, volume = {8}, year = {2017}, month = {Dec-01-2017}, pages = {8325-8335}, abstract = {

The process by which organic matter decomposes deep underground to form petroleum and its underlying kerogen matrix has so far remained a no man\&$\#$39;s land to theoreticians, largely because of the geological (Myears) timescale associated with the process. Using reactive molecular dynamics and an accelerated simulation framework, the replica exchange molecular dynamics method, we simulate the full transformation of cellulose into kerogen and its associated fluid phase under prevailing geological conditions. We observe in sequence the fragmentation of the cellulose crystal and production of water, the development of an unsaturated aliphatic macromolecular phase and its aromatization. The composition of the solid residue along the maturation pathway strictly follows what is observed for natural type III kerogen and for artificially matured samples under confined conditions. After expulsion of the fluid phase, the obtained microporous kerogen possesses the structure, texture, density, porosity and stiffness observed for mature type III kerogen and a microporous carbon obtained by saccharose pyrolysis at low temperature. As expected for this variety of precursor, the main resulting hydrocarbon is methane. The present work thus demonstrates that molecular simulations can now be used to assess, almost quantitatively, such complex chemical processes as petrogenesis in fossil reservoirs and, more generally, the possible conversion of any natural product into bio-sourced materials and/or fuel.

}, issn = {2041-6520}, doi = {10.1039/C7SC03466K}, url = {http://xlink.rsc.org/?DOI=C7SC03466K}, author = {Atmani, Lea and Christophe Bichara and Roland Jean-Marc Pellenq and Henri Van Damme and Adri CT Van Duin and Raza, Zamaan and Truflandier, Lionel A. and Ama{\"e}l Obliger and Kralert, Paul G. and Franz-Josef Ulm and Jean-Marc Leyssale} } @article {622, title = {The Potential of Mean Force concept for bridging (length and time) scales in the modeling of complex porous materials}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 01009}, abstract = {

We introduce the concept of Potential of Mean Force, PMF, as a way to implement upscaling modeling from the nano-scale to micron-scale. A PMF is a free energy function representing in an effective way the interactions between objects (cement hydrates, clay platelets, etc.) at thermodynamics conditions. The PMF is therefore the key piece of information allowing to coarse-grained Physical-Chemistry information in a meso-scale model formulation. The use of PMF offers a huge computational advantage as it allows a straight up-scaling to the meso-scale while keeping essential interactions information that are the hallmark of Physical-Chemistry processes. Such a coarse-grained modeling integrates atomistic response into inter-particle potentials that fully propagate molecular scale information all the way to the meso-scale.

}, doi = {10.1051/epjconf/201714001009}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714001009}, author = {Katerina Ioannidou and Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq}, editor = {Farhang Radja{\"\i} and Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @book {636, title = {Springer Series in Geomechanics and GeoengineeringAdvances in Laboratory Testing and Modelling of Soils and Shales (ATMSS)Measurement of Mechanical Properties of Thin Clay Films and Comparison with Molecular Simulations}, series = {Springer Series in Geomechanics and Geoengineering}, year = {2017}, pages = {78-84}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

Here, we focus on the hydromechanical behavior of self-standing clay films with a thickness of a few dozen microns. We measure their elastic and creep properties and how those properties depend on the interlayer cation and on the relative humidity (or water content). Those experimental results are compared with the elastic and creep behavior of nanometric clay particles, which we characterize by molecular simulations. Significant qualitative differences between the behavior of the clay films and that of the clay particles are observed, which suggests that the hydromechanical behavior of the clay films is significantly impacted by their mesostructure (i.e., by how the clay particles or tactoids are arranged in space). Upscaling the hydromechanical behavior of the clay films from that of the clay particles may be challenging.

}, isbn = {978-3-319-52772-7}, issn = {1866-8755}, doi = {10.1007/978-3-319-52773-410.1007/978-3-319-52773-4_8}, url = {http://link.springer.com/10.1007/978-3-319-52773-4}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Henri Van Damme}, editor = {Ferrari, Alessio and Laloui, Lyesse} } @article {599, title = {Strength of wet agglomerates of spherical particles: effects of friction and size distribution}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 08021}, abstract = {

We investigate the mechanical behavior of wet granular agglomerates composed of spherical particles by means of molecular dynamics simulations. The capillary cohesion force is modeled as an attraction force at the contact between two particles and expressed as an explicit function of the gap and volume of the liquid bridge. We are interested in the effect of the friction coefficient between primary particles. The agglomerates are subjected to diametrical compression tests. We find that the deformation is ductile involving particle rearrangements. However, a well-defined stress peak is observed and the peak stress is used as a measure of the compressive strength of the agglomerate. The strength increases with friction coefficient but levels off at friction coefficients above 0.4. Furthermore, the compressive strength is an increasing function of particle size span.

}, doi = {10.1051/epjconf/201714008021}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714008021}, author = {Vo, Thanh-Trung and Patrick Mutabaruka and Jean-Yves Delenne and Saeid Nezamabadi and Farhang Radja{\"\i}}, editor = {Luding, S.} } @article {609, title = {A time-dependent atomistic reconstruction of severe irradiation damage and associated property changes in nuclear graphite}, journal = {Carbon}, volume = {120}, year = {2017}, month = {Aug-2017}, pages = {111 - 120}, abstract = {

Wheat grains can be considered as a natural cemented granular material. They are milled under high forces to produce food products such as flour. The major part of the grain is the so-called starchy endosperm. It contains stiff starch granules, which show a multi-modal size distribution, and a softer protein matrix that surrounds the granules. Experimental milling studies and numerical simulations are going hand in hand to better understand the fragmentation behavior of this biological material and to improve milling performance. We present a numerical study of the effect of granule size distribution on the strength of such a cemented granular material. Samples of bi-modal starch granule size distribution were created and submitted to uniaxial tension, using a peridynamics method. We show that, when compared to the effects of starch-protein interface adhesion and voids, the granule size distribution has a limited effect on the samples\’ yield stress.

}, keywords = {TRANSMISSION ELECTRON-MICROSCOPY; MOLECULAR-DYNAMICS SIMULATIONS; HIGHLY ANISOTROPIC PYROCARBONS; THERMAL-CONDUCTIVITY; NEUTRON-IRRADIATION; CARBON; GRAPHENE; DEFECTS; ENERGY; NANOSTRUCTURE}, issn = {00086223}, doi = {10.1016/j.carbon.2017.05.009}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008622317304542}, author = {Baptiste Farbos and Freeman, Helen and Hardcastle, Trevor and Da Costa, Jean-Pierre and Brydson, Rik and Scott, Andrew J. and Weisbecker, Patrick and Germain Christian and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {637, title = {Adsorption in heterogeneous porous media: Hierarchical and composite solids}, journal = {Microporous and Mesoporous Materials}, volume = {229}, year = {2016}, month = {Jul-15-2016}, pages = {145 - 154}, abstract = {

Experiment and molecular simulation are used to investigate adsorption in heterogeneous porous media consisting of hierarchical solids (combining different porosity scales) or composite solids (such as silver nanoparticles adsorbed at the external surface of zeolite). It is shown that adsorption in such heterogeneous materials can be written as a linear combination of the adsorption isotherms in its different domains (i.e. porosity scales for the hierarchical sample and constituents for the composite sample). In the case of the composite material, we also show that the linear combination can be used with weighing parameters obtained for a different adsorbate. Such a superimposition principle, which is validated using well-characterized experimental samples, is of interest for characterization purpose as well as industrial applications as they can be used to determine accurately the amount of phases in a given sample (volume corresponding to a given porosity scale or constituent). In contrast, significant departure between the experimental adsorption isotherm and the linear combination can be used to detect coupling effects between the different domains or restrained access to a given domain type. Such a characterization strategy of complex heterogeneous media is complementary to other experiments, such as those probing capillary hysteresis shapes, scanning curves and subloops, which allow determining the distribution of domains within the framework of the independent domain theory. (C) 2016 Elsevier Inc. All rights reserved.

}, issn = {13871811}, doi = {10.1016/j.micromeso.2016.04.011}, url = {https://linkinghub.elsevier.com/retrieve/pii/S1387181116301032}, author = {Deliere, Ludovic and Fran{\c c}ois Villemot and David Farrusseng and Anne Galarneau and Topin, Sylvain and Benoit A. Coasne} } @article {217, title = {Effect of Water on Elastic and Creep Properties of Self-Standing Clay Films}, journal = {Langmuir}, volume = {32}, year = {2016}, month = {Feb-09-2016}, pages = {1370 - 1379}, abstract = {

We characterized experimentally the elastic and creep properties of thin self-standing clay films, and how their mechanical properties evolved with relative humidity and water content. The films were made of clay montmorillonite SWy-2, obtained by evaporation of a clay suspension. Three types of films were manufactured, which differed by their interlayer cation: sodium, calcium, or a mixture of sodium with calcium. The orientational order of the films was characterized by X-ray diffractometry. The films were mechanically solicited in tension, the resulting strains being measured by digital image correlation. We measured the Young\&$\#$39;s modulus and the creep over a variety of relative humidities, on a full cycle of adsorption-desorption for what concerns the Young\&$\#$39;s modulus. Increasing relative humidity made the films less stiff and made them creep more. Both the elastic and creep properties depended significantly on the interlayer cation. For the Young\&$\#$39;s modulus, this dependence must originate from a scale greater than the scale of the clay layer. Also, hysteresis disappeared when plotting the Young\&$\#$39;s modulus versus water content instead of relative humidity. Independent of interlayer cation and of relative humidity greater than 60\%, after a transient period, the creep of the films was always a logarithmic function of time. The experimental data gathered on these mesoscale systems can be of value for modelers who aim at predicting the mechanical behavior of clay-based materials (e.g., shales) at the engineering macroscopic scale from the one at the atomistic scale, for them to validate the first steps of their upscaling scheme. They provide also valuable reference data for bioinspired clay-based hybrid materials.

}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.5b03431}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Bornert, Michel and Ferrage, Eric and Hubert, Fabien and Henri Van Damme} } @article {247, title = {Evolution of organo-clay composites with respect to thermal maturity in type II organic-rich source rocks}, journal = {Geochimica et Cosmochimica Acta}, volume = {195}, year = {2016}, month = {Dec-15-2016}, pages = {68 - 83}, abstract = {

Among the proposed pathways of sequestration and progressive transformation of organic matter (OM) during burial, the sorptive protection mechanism explains the strong relationship between total organic carbon (TOC) and mineral surface area (MSA) noted in numerous black shales around the globe. The complete mechanistic framework of preservation and evolution of OM in organo-mineral complexes remains, however, an enigma. On one hand, smectite layers composing the clay fraction are known to have a major influence on available surface area. OM maturation occurs, on the other hand, concurrently with the smectite illitization that provokes the closure of the interlayer spaces. The potential of smectite layers in the sequestration and preservation of organic molecules and the fate of these molecules with respect to the smectite illitization were therefore addressed. Here, the mineralogy of three organic-rich source rocks of various maturities was characterized in regards with the geochemistry of their OM. A thorough examination of the clay minerals present in the clay matrices provided evidences of mixed layer minerals containing smectite and illite layers with an increasing illite component with respect to maturity. The comprehensive interpretation of the X-ray diffractograms and analytical electron microscopy results suggested the presence of organic molecules in the inter-particulate and possibly the interlayer spaces of the smectite-rich components in immature source rocks. This eventuality was further supported by the presence of intercalated clay-organic nanocomposites observed by transmitted electron microscopy coupled with energy dispersive spectroscopy. Textural observations also showed that the increased illite content found in the overmature sample led to the reorganization of the OM and the clay particles into nano-scale aggregates. These results clarify the geochemical mechanism beyond the reported relationship between TOC and MSA and allow generalizing it to various organic-rich source rocks. Eventually, the understanding of the close organo-clay association raises prospects on the prediction of the evolution of the mechanical properties of economically valuable source rocks with respect to maturity.

}, issn = {00167037}, doi = {10.1016/j.gca.2016.09.008}, author = {Berthonneau, Jeremie and Grauby, Olivier and Muhannad Abuhaikal and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Henri Van Damme} } @article {195, title = {Realistic molecular model of kerogen{\textquoteright}s nanostructure}, journal = {Nature Materials}, year = {2016}, month = {Feb-01-2016}, abstract = {

Despite kerogen\’s importance as the organic backbone for hydrocarbon production from source rocks such as gas shale, the interplay between kerogen\’s chemistry, morphology and mechanics remains unexplored. As the environmental impact of shale gas rises, identifying functional relations between its geochemical, transport, elastic and fracture properties from realistic molecular models of kerogens becomes all the more important. Here, by using a hybrid experimental\–simulation method, we propose a panel of realistic molecular models of mature and immature kerogens that provide a detailed picture of kerogen\’s nanostructure without considering the presence of clays and other minerals in shales. We probe the models\’ strengths and limitations, and show that they predict essential features amenable to experimental validation, including pore distribution, vibrational density of states and stiffness. We also show that kerogen\’s maturation, which manifests itself as an increase in the sp2/sp3 hybridization ratio, entails a crossover from plastic-to-brittle rupture mechanisms.

Kerogen in organic-rich shale formations.

Reconstruction of molecular models of kerogens.

Experimental validation of the molecular models.

Elastic properties of molecular models of kerogen.

}, issn = {1476-1122}, doi = {10.1038/nmat4541}, author = {Bousige, Colin and Ghimbeu, Cam{\'e}lia Matei and Vix-Guterl, Cathie and Pomerantz, Andrew E. and Suleimenova, Assiya and Vaughan, Gavin and Garbarino, Gaston and Feygenson, Mikhail and Wildgruber, Christoph and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @proceedings {333, title = {Creep of Clay: Numerical Results at the Scale of a Layer and Experimental Results at the Scale of Thin Self-Standing Films}, journal = {10th International Conference on Mechanics and Physics of Creep, Shrinkage, and Durability of Concrete and Concrete StructuresCONCREEP 10}, volume = {CONCREEP 10: MECHANICS AND PHYSICS OF CREEP, SHRINKAGE, AND DURABILITY OF CONCRETE AND CONCRETE STRUCTURES }, year = {2015}, month = {Sep-17-2015}, pages = {531-536}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA}, abstract = {

This work focuses on the creep of clay-based materials, which exhibit significant analogies with cement-based materials. Here, we studied the creep of clay at two scales and with two techniques: numerically (with molecular simulations) at the scale of a clay layer (nm), and experimentally at the scale of thin self-standing clay films (few dozen m). At the scale of the clay layer, numerical simulations showed that the shear rate was constant over time and an affine function of the shear stress. Creep experiments showed that, after a transient period, the creep function of our thin self-standing clay films was a logarithmic function of time. A comparison of the results obtained at the two scales shows that the origin of the logarithmic feature of clay creep must at least partly originate from a scale greater than that of an individual clay layer. By analogy, such result is likely to hold for cementitious materials, which are also known to creep logarithmically with respect to time in the long term: the origin of this logarithmic feature is likely to stem at least partly from a scale greater than the scale of an individual C-S-H layer.

}, doi = {10.1061/978078447934610.1061/9780784479346.064}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Henri Van Damme}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @article {163, title = {Diffusion of Interlayer Cations in Swelling Clays as a Function of Water Content: Case of Montmorillonites Saturated with Alkali Cations}, journal = {The Journal of Physical Chemistry C}, volume = {119}, year = {2015}, month = {May-14-2015}, pages = {10370 - 10378}, abstract = {
Abstract Image

Swelling clays act as ion-exchange membranes owing to the presence of extra-framework ions surrounded by water molecules. A better understanding of the motion of the ions present in the pore space is required to predict the diffusion properties in these solids. Since water molecules tend to adsorb both on ions and clay surfaces, they moderate the interactions between the ions and the framework. The hydration of ions therefore has an impact on their diffusion properties. In this paper, the MX-80 montmorillonite considered for nuclear waste disposal applications has been selected and studied using an approach combining measurements of complex impedance spectroscopy and water adsorption isotherms. The number of charge carriers has been estimated, and the diffusion coefficient for interlayer cations at various hydration states was determined. The evolution of the diffusion coefficient is subsequently correlated with the effect of the opening of the interlayer space and its hydration state. The resulting picture sheds light on the reversible or irreversible character of ion exchange. Finally, our results obtained from the difference between the interlayer diffusion coefficients extracted from conductivity measurements and those obtained at the macroscopic scale is discussed in terms of the textural parameters.

}, issn = {1932-7447}, doi = {10.1021/jp512986d}, author = {Salles, Fabrice and Douillard, Jean-Marc and Bildstein, Olivier and El Ghazi, Samira and Prelot, B{\'e}n{\'e}dicte and Zajac, Jerzy and Henri Van Damme} } @article {644, title = {First-principles prediction of kink-pair activation enthalpy on screw dislocations in bcc transition metals: V, Nb, Ta, Mo, W, and Fe}, journal = {Physical Review B}, volume = {91}, year = {2015}, month = {Mar-09-2015}, pages = { Article Number 094105 }, abstract = {

The atomistic study of kink pairs on screw dislocations in body-centered cubic (bcc) metals is challenging because interatomic potentials in bcc metals still lack accuracy and kink pairs require too many atoms to be modeled by first principles. Here, we circumvent this difficulty using a one-dimensional line tension model whose parameters, namely the line tension and Peierls barrier, are reachable to density functional theory calculations. The model parameterized in V, Nb, Ta, Mo, W, and Fe, is used to study the kink-pair activation enthalpy and spatial extension. Interestingly, we find that the atomistic line tension is more than twice the usual elastic estimates. The calculations also show interesting group tendencies with the line tension and kink-pair width larger in group V than in group VI elements. Finally, the present kink-pair activation energies are shown to compare qualitatively with experimental data and potential origins of quantitative discrepancies are discussed.

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.91.094105}, url = {https://link.aps.org/doi/10.1103/PhysRevB.91.094105}, author = {Dezerald, Lucile and Proville, L. and Ventelon, Lisa and Willaime, F. and Rodney, D.} } @article {162, title = {Fractal analysis of surface roughness of montmorillonite clay self-supported films: Effects of exchanged cations and of mechanical tensile stress}, journal = {Colloids and Surfaces A: Physicochemical and Engineering Aspects}, volume = {486}, year = {2015}, month = {Dec-05-2015}, pages = {38 - 44}, abstract = {

In this paper, surface roughness of self-supported montmorillonite clay films was studied by the use of fractal analysis to study its links with various inter-foliar cations, mechanical properties, and permeability. The Hurst exponent of fractional Brownian motion was estimated by maximum likelihood from optical images in reflection mode. Images corresponding to films with different inter-foliar exchanged cations visually show significant texture variations. It is confirmed by a continuous decrease of H parameter toward less correlated textures as one goes from the samples exchanged with monovalent ions to those exchanged with divalent and finally trivalent ions. The same films were analyzed while they were submitted to mechanical tensile stress. The films exchanged with the most charged cations or less hydrated ones do not show any evidence of texture change and the H parameter is almost constant. Concerning the permeability study, less permeable parent cakes \– with hydrated monovalent cations \– lead to relatively smooth and less disturbed surface films with high H parameter. In the opposite, very permeable cakes give birth to anti-persistent relief.

}, issn = {09277757}, doi = {10.1016/j.colsurfa.2015.09.001}, author = {Zabat, M. and Harba, R. and Henri Van Damme} } @article {146, title = {Investigating carbon materials nanostructure using image orientation statistics}, journal = {Carbon}, volume = {84}, year = {2015}, month = {Apr-2015}, pages = {160 - 173}, abstract = {

A new characterization method of the lattice fringe images of turbostratic carbons is proposed. This method is based on the computation of their orientation field without explicit detection of fringes. It allows meaningful insights into the material nanostructure and nanotexture at several scales, either qualitatively or quantitatively. The calculation of pairwise spatial statistics of the orientation field at short distance provides measurements of the coherence lengths along any direction, in particular along and orthogonally to the layers. These statistics also allow representing orientation coherence patterns typical of the observed nanostructure. At larger distances, the mean disorientation of the fringes is computed and information about the homogeneity of the sample is obtained. An experimental validation is carried out on various artificial images and an application to the characterization of four bulk turbostratic carbons is provided.

}, issn = {00086223}, doi = {10.1016/j.carbon.2014.11.048}, author = {J. Philippe Da-Costa and Weisbecker, P. and Baptiste Farbos and Jean-Marc Leyssale and G{\'e}rard L. Vignoles and Germain, C.} } @article {185, title = {The Maya blue nanostructured material concept applied to colouring geopolymers}, journal = {RSC Adv.}, volume = {5}, year = {2015}, month = {Nov-2015}, pages = {98834 - 98841}, abstract = {

Maya blue is an ancient nanostructured pigment synthetized by assembling indigo, a natural dye, with palygorskite, a microfibrous clay mineral. The novelty of our approach is to mimic \“pre-Columbian nanotechnology\” and to functionalize geopolymers with a sepiolite-based hybrid organic\–inorganic nanocomposite inspired from the Maya blue. It is acid- and UV-resistant, as confirmed by the stability of Maya mural paintings over time. We synthesized analogous pigments, using methylene blue (MB) and methyl red (MR) as organic dyes and sepiolite as fibrous clay mineral. We used an aqueous and a solid-state method, both leading to encapsulation of dye monomers into the clay micropores, as confirmed by UV-vis spectroscopy. This nanostructured pigment was then included into a geopolymer matrix at room temperature. The stability of the new material to UV and acid was tested. It was confirmed that it is the prior encapsulation of the dye into sepiolite that leads to the stability of the pigment in the geopolymer matrix. This first study opens the way to numerous possibilities for functionalizing inorganic binder materials with organic elements that would be otherwise sensitive to thermal treatment in conventional ceramic processing.

}, doi = {10.1039/C5RA14076E}, author = {Ouellet-Plamondon, Claudiane and Aranda, Pilar and Favier, Aur{\'e}lie and Habert, Guillaume and Henri Van Damme and Ruiz-Hitzky, Eduardo} } @article {145, title = {Mechanism of strength reduction along the graphenization pathway}, journal = {Science Advances}, volume = {1}, year = {2015}, month = {Nov-20-2015}, pages = {e1501009 - e1501009}, abstract = {

Even though polycrystalline graphene has shown a surprisingly high tensile strength, the influence of inherent grain boundaries on such property remains unclear. We study the fracture properties of a series of polycrystalline graphene models of increasing thermodynamic stability, as obtained from a long molecular dynamics simulation at an elevated temperature. All of the models show the typical and well-documented brittle fracture behavior of polycrystalline graphene; however, a clear decrease in all fracture properties is observed with increasing annealing time. The remarkably high fracture properties obtained for the most disordered (less annealed) structures arise from the formation of many nonpropagating prefracture cracks, significantly retarding failure. The stability of these reversible cracks is due to the nonlocal character of load transfer after a bond rupture in very disordered systems. It results in an insufficient strain level on neighboring bonds to promote fracture propagation. Although polycrystallinity seems to be an unavoidable feature of chemically synthesized graphenes, these results suggest that targeting highly disordered states might be a convenient way to obtain improved mechanical properties.

}, doi = {10.1126/sciadv.1501009}, author = {Antonio Gamboa and Baptiste Farbos and Aurel, P. and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {143, title = {Nanoscale elasticity of highly anisotropic pyrocarbons}, journal = {Carbon}, volume = {94}, year = {2015}, month = {Nov-2015}, pages = {285 - 294}, abstract = {

We report on the elastic properties of high-textured laminar pyrocarbons (PyCs) as obtained at the nanoscale using atomistic simulations on realistic models of the rough laminar (RL) PyC, the regenerative laminar (ReL) PyC, and the ReL PyC heat treated at temperatures up to 1700 \°C. The purely longitudinal ( and ) and transverse ( and ) elastic properties of such materials have values of around 50\–75\% of those found for graphite. Conversely, cross longitudinal-transverse properties ( and ) are much larger in PyCs than in graphite (up to around one order of magnitude for ). Our results also show that stiffness decreases with the hydrogen content, that longitudinal properties increase with the extent of graphene domains and decrease with the misorientation of the layers, and that transverse and cross properties significantly increase with the amount of interlayer cross-links present in the materials. Comparison to known experimental data indicates that texture effect at superior scales play an important role in the effective macroscopic elasticity constants. Finally, from a materials perspective and in agreement with the structural evolution, the elastic properties of the ReL PyC evolve with heat treatment towards those of the RL PyC.

}, issn = {00086223}, doi = {10.1016/j.carbon.2015.06.060}, author = {Baptiste Farbos and J. Philippe Da-Costa and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {178, title = {Polarized Water Wires under Confinement in Chiral Channels}, journal = {The Journal of Physical Chemistry B}, volume = {119}, year = {2015}, month = {Jul-16-2015}, pages = {8707 - 8717}, abstract = {
Abstract Image

The alignment of water molecules along chiral pores may activate proton/ion conduction along dipolar hydrophilic pathways. Here we show that a simple synthetic \“T-channel\” forms a directional pore with its carbonyl moieties solvated by chiral helical water wires. Atom-scale simulations and experimental crystallographic assays reveal a dynamical structure of water and electrolyte solutions (alkali metal chlorides) confined in these organic T-channels. Oscillations in the dipole orientation, which correspond to alternative ordering (dipole up\–dipole down) of the water molecules with a period of about 4.2 \Å (imposed by the distance between two successive carbonyl groups) are observed. When ions are added to the system, despite the strong Coulombic water/ion interaction, confined water remains significantly ordered in the T-channel and still exhibits surface-induced polarization. Cation permeation can be achieved through alternated hydration\–dehydration occurring along strongly oriented water wires. The T-channel, which exhibits chirality with strong water orientation, provides an opportunity to unravel novel water-channel systems that share many interesting properties of biomolecular systems.

}, issn = {1520-6106}, doi = {10.1021/acs.jpcb.5b03322}, author = {Barboiu, Mihail and Pierre-Andre Cazade and Le Duc, Yann and Yves-Marie Legrand and Arie Van Der Lee and Benoit A. Coasne} } @article {144, title = {On the prediction of graphene{\textquoteright}s elastic properties with reactive empirical bond order potentials}, journal = {Carbon}, volume = {89}, year = {2015}, month = {Aug-2015}, pages = {176 - 187}, abstract = {

The elastic properties of graphene as described by the reactive empirical bond order potential are studied through uniaxial tensile tests calculations at both zero temperature, with a conjugate gradient approach, and room temperature, with molecular dynamics simulations. A perfect linear elastic behavior is observed at 0\ K up to \≈0.1\% strain. The Young\’s modulus and Poisson\’s ratio obtained with this potential are of \≈730\ GPa and 0.39, respectively, with little chirality effects. These values differ significantly from former estimations, much closer to experimental values. We show that these former values have certainly been obtained by neglecting the effect of atomic relaxation, leading to a severe inaccuracy. At larger strains, an extended apparent linear domain is observed in the stress\–strain curves, which is relevant to Young\’s modulus calculations at finite temperature. Our molecular dynamics simulations at 300\ K have allowed obtaining the following, chirality dependent, apparent Young\’s moduli, 860 and 761\ GPa, and Poisson\’s ratios, 0.12 and 0.23, for armchair and zigzag loadings, respectively.

}, issn = {00086223}, doi = {10.1016/j.carbon.2015.03.035}, author = {Antonio Gamboa and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {108, title = {Adsorption and Dynamics in Hierarchical Metal-Organic Frameworks}, journal = {Journal of Physical Chemistry C}, volume = {118}, year = {2014}, month = {Apr-10-2014}, pages = {7423-7433}, type = {Article}, abstract = {

Adsorption and dynamics in hierarchical metal-organic frameworks are investigated by means of molecular simulation. The models of hierarchical porous solids are obtained by carving mesopores of different diameters out of a crystal of Cu-BTC (model A) or by inserting a microporous particle of Cu-BTC in an amorphous silica mesopore (model B). We show that the nitrogen adsorption isotherms at 77 K for the solids corresponding to model A can be described as a linear combination of reference adsorption isotherms for pure microporous and mesoporous solids. In contrast, the adsorption isotherms for model B cannot be described accurately as a sum of reference microporous and mesoporous adsorption isotherms. The inserted particle acts as a constriction which helps nucleate the liquid phase within the mesopore so that no capillary condensation hysteresis is observed. The dynamics of nitrogen adsorbed at 77 K inside the porosity of the hierarchical solids is also investigated. The Fickian regime is reached at long times which are not attainable with molecular dynamics simulations. At higher temperature, the faster self-diffusion makes it possible to obtain the diffusivity of the adsorbate. Nitrogen adsorbed in the microporosity of the hierarchical porous solids has a self-diffusion coefficient close to that of nitrogen adsorbed in pure Cu-BTC. In contrast, diffusion in the mesoporosity is faster than in the microporosity so that the overall diffusivity is faster than in pure Cu-BTC.

\<\>Abstract Image

}, issn = {1932-7447}, doi = {10.1021/jp500037z}, author = {Fran{\c c}ois Villemot and Anne Galarneau and Benoit A. Coasne} } @article {64, title = {Adsorption-based characterization of hierarchical metal{\textendash}organic frameworks}, journal = {Adsorption}, volume = {20}, year = {2014}, month = {Feb-2014}, pages = {349-357}, abstract = {

Nitrogen adsorption at 77 K on metal\–organic framework (MOF) is investigated by means of molecular simulations. We consider both regular Cu\–BTC crystal and a MOF-based hierarchical porous solid consisting of a mesopore carved out of a Cu\–BTC crystal. The t-plot method is applied to these solids by using a non-porous Cu\–BTC surface as the reference sample. The values of the mesoporous and external surface areas are determined from the t-plot, and the validity of the method for this type of hierarchical solid is discussed.

}, keywords = {Adsorption, Hierarchical porous materials, Molecular simulation, t-Plot}, issn = {0929-5607}, doi = {10.1007/s10450-013-9583-7}, author = {Fran{\c c}ois Villemot and Anne Galarneau and Benoit A. Coasne} } @article {107, title = {An artificial primitive mimic of the Gramicidin-A channel}, journal = {Nature Communications}, volume = {5}, year = {2014}, month = {Jun-26-2014}, type = {Article}, abstract = {

Gramicidin A (gA) is the simplest known natural channel, and important progress in improving conduction activity has previously been obtained with modified natural gAs. However, simple artificial systems mimicking the gA functions are unknown. Here we show that gA can be mimicked using a simple synthetic triazole or {\textquoteleft}T-channel\&$\#$39; forming compound (TCT), having similar constitutional functions as the natural gAs. As in gA channels, the carbonyl moieties of the TCT, which point toward the T-channel core and surround the transport direction, are solvated by water. The net-dipolar alignment of water molecules along the chiral pore surfaces influences the conduction of protons/ions, envisioned to diffuse along dipolar hydrophilic pathways. Theoretical simulations and experimental assays reveal that the conduction through the T-channel, similar to that in gA, presents proton/water conduction, cation/anion selectivity and large open channel-conductance states. T-channels-associating supramolecular chirality with dipolar water alignment-represent an artificial primitive mimic of gA.

Molecular structures and crystal packing.

Molecular simulations of water molecules and ions confined within T-channel.

Bilayer membrane transport experiments.

Single channel conductance of T-channels under {\textquoteleft}early stage{\textquoteright} conductance records.

}, issn = {2041-1723}, doi = {10.1038/ncomms5142}, author = {Barboiu, Mihail and Le Duc, Yann and Gilles, Arnaud and Pierre-Andre Cazade and Michau, Mathieu and Yves-Marie Legrand and Arie Van Der Lee and Benoit A. Coasne and Parvizi, Paria and Post, Joshua and Fyles, Thomas} } @article {72, title = {Combinatorial molecular optimization of cement hydrates}, journal = {Nat Commun}, volume = {5}, year = {2014}, month = {Sep-24-2014}, pages = {Article Number: 4960}, abstract = {

Despite its ubiquitous presence in the built environment, concrete\â\€$(1s (Bmolecular-level properties are only recently being explored using experimental and simulation studies. Increasing societal concerns about concrete\â\€$(1s (Benvironmental footprint have provided strong motivation to develop new concrete with greater specific stiffness or strength (for structures with less material). Herein, a combinatorial approach is described to optimize properties of cement hydrates. The method entails screening a computationally generated database of atomic structures of calcium-silicate-hydrate, the binding phase of concrete, against a set of three defect attributes: calcium-to-silicon ratio as compositional index and two correlation distances describing medium-range silicon-oxygen and calcium-oxygen environments. Although structural and mechanical properties correlate well with calcium-to-silicon ratio, the cross-correlation between all three defect attributes reveals an indentation modulus-to-hardness ratio extremum, analogous to identifying optimum network connectivity in glass rheology. We also comment on implications of the present findings for a novel route to optimize the nanoscale mechanical properties of cement hydrate.

Figure 1

}, doi = {10.1038/ncomms5960}, author = {Mohammad Javad Abdolhosseini Qomi and Konrad J. Krakowiak and Mathieu Bauchy and Stewart, K.L. and Rouzbeh Shahsavari and Jagannathan, D. and Brommer, D.B. and Alain Baronnet and Markus J Buehler and Sidney Yip and Franz-Josef Ulm and Krystyn J. Van Vliet and Roland Jean-Marc Pellenq} } @article {114, title = {Elastic Properties of Swelling Clay Particles at Finite Temperature upon Hydration}, journal = {Journal of Physical Chemistry C}, volume = {118}, year = {2014}, month = {May-01-2014}, pages = {8933-8943}, type = {Article}, abstract = {
Abstract Image

The knowledge of mechanical properties of swelling clays on the scale of the clay layer is of crucial importance to build predictive hydromechanical macroscopic constitutive laws in a bottom-up approach. In this work, we computed the elastic properties of particles of a swelling clay (i.e., montmorillonite) by molecular dynamics simulations. Because of the softness of the material in the direction orthogonal to the clay layer with increasing water content, the computation of the whole stiffness tensor was not trivial: we needed to implement the elastic bath method,1 which allows us to attenuate thermal strains and compute the stiffness via strain fluctuations in isothermal\–isobaric simulations. We investigated the effect of water content, temperature, and the interlayer cation on the elastic properties of swelling clays. In particular, we showed that the out-of-plane stiffness coefficients computed at 300 K differed significantly from the same coefficients computed at 0 K. The out-of-plane coefficients were very sensitive to both temperature and water content. In contrast, the dependence of the in-plane coefficients to temperature was slight. Moreover, the decrease of the in-plane coefficients with water content could be entirely explained by the change of geometry of the system due to the swelling of the interlayer space. The interlayer cation impacted the elastic properties of montmorillonite only in the driest states.

}, issn = {1932-7447}, doi = {10.1021/jp412160e}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Henri Van Damme} } @article {83, title = {Ionic liquid confined in silica nanopores: molecular dynamics in the isobaric-isothermal ensemble}, journal = {Molecular Physics}, volume = {112}, year = {2014}, month = {Avr-03-2014}, pages = {1350-1361}, type = {Article}, abstract = {

Molecular dynamics simulations in the isobaric-isothermal ensemble are used to investigate the structure and dynamics of an ionic liquid confined at ambient temperature and pressure in hydroxylated amorphous silica nanopores. The use of the isobaric-isothermal ensemble allows estimating the effect of confinement and surface chemistry on the density of the confined ionic liquid. The structure of the confined ionic liquid is investigated using density profiles and structural order parameters while its dynamics is assessed by determining the mobility and ionic conductivity of the confined phase. Despite the important screening of the electrostatic interactions (owing to the small Debye length in ionic liquids), the local structure of the confined ionic liquid is found to be mostly driven by electrostatic interactions. We show that both the structure and dynamics of the confined ionic liquid can be described as the sum of a surface contribution arising from the ions in contact with the surface and a bulk-like contribution arising from the ions located in the pore centre; as a result, most properties of the confined ionic liquid are a simple function of the surface-to-volume ratio of the host porous material. In contrast, the ionic conductivity of the confined ionic liquid, which is a collective dynamical property, is found to be similar to the bulk. This study sheds light on the complex behaviour of hybrid materials made up of ionic liquid confined in inorganic porous materials.

}, keywords = {amorphous silica, ionic liquids, molecular dynamics, surface}, issn = {0026-8976}, doi = {10.1080/00268976.2014.902138}, author = {Guido Ori and Fran{\c c}ois Villemot and Lydie Viau and Andre Vioux and Benoit A. Coasne} } @article {649, title = {Poroelastic Theory Applied to the Adsorption-Induced Deformation of Vitreous Silica}, journal = {The Journal of Physical Chemistry B}, volume = {118}, year = {2014}, month = {Dec-11-2014}, pages = {14519 - 14525}, abstract = {

When vitreous silica is submitted to high pressures under a helium atmosphere, the change in volume observed is much smaller than expected from its elastic properties. It results from helium penetration into the interstitial free volume of the glass network. We present here the results of concurrent spectroscopic experiments using either helium or neon and molecular simulations relating the amount of gas adsorbed to the strain of the network. We show that a generalized poromechanical approach, describing the elastic properties of microporous materials upon adsorption, can be applied successfully to silica glass in which the free volume exists only at the subnanometer scale. In that picture, the adsorption-induced deformation accounts for the small apparent compressibility of silica observed in experiments.

}, issn = {1520-6106}, doi = {10.1021/jp5094383}, url = {http://pubs.acs.org/doi/10.1021/jp5094383}, author = {Benoit A. Coasne and Weigel, Coralie and Polian, Alain and Kint, Mathieu and Rouquette, J{\'e}r{\^o}me and Haines, Julien and Foret, Marie and Vacher, Rene and Ruffle, Benoit} } @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 {74, title = {Ultrahigh interlayer friction in multiwalled boron nitride nanotubes}, journal = {Nat Mater}, volume = {13}, year = {2014}, month = {Jul-2014}, pages = {688 - 693}, abstract = {

Friction at the nanoscale has revealed a wealth of behaviours that depart strongly from the long-standing macroscopic laws of Amontons\–Coulomb1, 2. Here, by using a \‘Christmas cracker\’-type of system in which a multiwalled nanotube is torn apart between a quartz-tuning-fork-based atomic force microscope (TF\–AFM) and a nanomanipulator, we compare the mechanical response of multiwalled carbon nanotubes (CNTs) and multiwalled boron nitride nanotubes (BNNTs) during the fracture and telescopic sliding of the layers. We found that the interlayer friction for insulating BNNTs results in ultrahigh viscous-like dissipation that is proportional to the contact area, whereas for the semimetallic CNTs the sliding friction vanishes within experimental uncertainty. We ascribe this difference to the ionic character of the BN, which allows charge localization. The interlayer viscous friction of BNNTs suggests that BNNT membranes could serve as extremely efficient shock-absorbing\ surfaces.

Multiwalled nanotube Christmas-cracker experimental set-up.

Frequency shift measured by the tuning fork versus the displacement of the piezoscanner during complete tensile load experiments.

}, isbn = {1476-1122}, doi = {10.1038/nmat3985}, author = {Nigues, A. and A. Siria and Vincent, P. and Poncharal, P. and Lyd{\'e}ric Bocquet} } @article {92, title = {Validity of the t-plot Method to Assess Microporosity in Hierarchical Micro/Mesoporous Materials}, journal = {Langmuir}, volume = {30}, year = {2014}, month = {Nov-2014}, pages = {13266-13274}, type = {Article}, abstract = {
Abstract Image

The t-plot method is a well-known technique which allows determining the micro- and/or mesoporous volumes and the specific surface area of a sample by comparison with a reference adsorption isotherm of a nonporous material having the same surface chemistry. In this paper, the validity of the t-plot method is discussed in the case of hierarchical porous materials exhibiting both micro- and mesoporosities. Different hierarchical zeolites with MCM-41 type ordered mesoporosity are prepared using pseudomorphic transformation. For comparison, we also consider simple mechanical mixtures of microporous and mesoporous materials. We first show an intrinsic failure of the t-plot method; this method does not describe the fact that, for a given surface chemistry and pressure, the thickness of the film adsorbed in micropores or small mesopores (\< 10\σ, \σ being the diameter of the adsorbate) increases with decreasing the pore size (curvature effect). We further show that such an effect, which arises from the fact that the surface area and, hence, the free energy of the curved gas/liquid interface decreases with increasing the film thickness, is captured using the simple thermodynamical model by Derjaguin. The effect of such a drawback on the ability of the t-plot method to estimate the micro- and mesoporous volumes of hierarchical samples is then discussed, and an abacus is given to correct the underestimated microporous volume by the t-plot method.

}, issn = {0743-7463}, doi = {10.1021/la5026679}, author = {Anne Galarneau and Fran{\c c}ois Villemot and Rodriguez, Jeremy and Fajula, Fran{\c c}ois and Benoit A. Coasne} } @article {81, title = {Adsorption-based characterization of hierarchical metal-organic frameworks}, journal = {Adsorptio-Journal of the International Adsorption Society}, volume = {20}, year = {2013}, month = {Oct-06-2013}, pages = {349-357}, type = {Article}, abstract = {

Nitrogen adsorption at 77 K on metal-organic framework (MOF) is investigated by means of molecular simulations. We consider both regular Cu-BTC crystal and a MOF-based hierarchical porous solid consisting of a mesopore carved out of a Cu-BTC crystal. The t-plot method is applied to these solids by using a non-porous Cu-BTC surface as the reference sample. The values of the mesoporous and external surface areas are determined from the t-plot, and the validity of the method for this type of hierarchical solid is discussed.

}, keywords = {Adsorption, Hierarchical porous materials, Molecular simulation, t-Plot}, issn = {0929-5607}, doi = {10.1007/s10450-013-9583-7}, author = {Fran{\c c}ois Villemot and Anne Galarneau and Benoit A. Coasne} } @proceedings {367, title = {A benchmark for particle shape dependence}, journal = {7th International Conference on Micromechanics of Granular Media (Powders and Grains)}, volume = {Book Series: AIP Conference Proceedings POWDERS AND GRAINS 2013}, year = {2013}, month = {Jun-18-2013}, pages = {883-886}, publisher = {AIP}, address = {JUL 08-12 2013 Sydney, AUSTRALIA}, abstract = {

Particle shape is a major parameter for the space-filling and strength properties of granular materials. For a systematic investigation of shape effect, a numerical benchmark test was set up within a collaborative group using different numerical methods and particles of various shape characteristics such as elongation, angularity and nonconvexity. Extensive 2D shear simulations were performed in this framework and the shear strength and packing fraction were compared for different shapes. We show that the results may be analyzed in terms of a low-order shape parameter \η describing the degree of distortion from a perfectly circular shape. In particular, the shear strength is an increasing function of \η with nearly the same trend for all shapes, the differences being of second order compared to \η. We also observe a nontrivial behavior of packing fraction which, for all our simulated shapes, increases with \η from the random close packing fraction for disks, reaches a peak considerably higher than that for disks, and subsequently declines as \η is further increased. Finally, the analysis of contact forces for the same value of \η leads to very similar statistics regardless of our specific particle shapes.

}, doi = {10.1063/1.4812073}, author = {Gael Combe and C{\'e}cile Nouguier-Lehon and Emilien Az{\'e}ma and Krzysztof Szarf and Baptiste Saint-Cyr and Marie Chaze and Farhang Radja{\"\i} and Pascal Villard and Jean-Yves Delenne and Vincent Richefeu and Philippe Sornay and Charles Voivret and CEGEO Group}, editor = {Yu, A and Dong, K and Yang, R} } @article {46, title = {ESEM Study of the Humidity-Induced Swelling of Clay Film}, journal = {Langmuir}, volume = {29}, year = {2013}, month = {Oct-02-2013}, pages = {12823{\textendash}12833}, abstract = {
Abstract Image

We measured the humidity-induced swelling of thin self-standing films of montmorillonite clay by a combination of environmental scanning electron microscopy (ESEM) and digital image correlation (DIC). The films were about 40 \μm thick. They were prepared by depositing and evaporating a suspension of clay and peeling off the highly oriented deposits. The rationale for creating such original samples was to obtain mesoscopic samples that could be used to bridge experimentally the gap between the scale of the clay layer and the engineering scale of a macroscopic clay sample. Several montmorillonite samples were used: the reference clay Swy-2, the same clay homoionized with sodium or calcium ions, and a sodium-exchanged Cloisite. The edges of the clay films were observed by ESEM at various relative humidity values between 14\% and 95\%. The ESEM images were then analyzed by DIC to measure the swelling or the shrinkage of the films. We also measured the adsorption/desorption isotherms by weighing the film samples in a humidity-controlled environment. In order to analyze our results, we compared our swelling/shrinkage and adsorption/desorption data with previously published data on the interlayer spacing obtained by X-ray diffraction and with numerical estimates of the interlayer water obtained by molecular dynamics simulation. The swelling and the hysteresis of this swelling were found to be comparable for the overall macroscopic films and for the interlayer space. The same correspondence between film and interlayer space was observed for the amount of adsorbed water. This suggests that, in the range of relative humidities values explored, the films behave like freely swelling oriented stacks of clay layers, without any significant contribution from the mesoporosity. The relevance of this result for the behavior of clayey sedimentary rocks and the differences with the behavior of nonoriented samples (powders or compacted powders) are briefly discussed.

}, doi = {10.1021/la402781p}, author = {Benoit Carrier and Wang, Linlin and Matthieu Vandamme and Roland Jean-Marc Pellenq and Bornert, Michel and Tanguy, Alexandre and Henri Van Damme} } @article {34, title = {Molecular simulation of adsorption and transport in hierarchical porous materials}, journal = {Langmuir}, volume = {29}, year = {2013}, month = {May-29-2013}, pages = {7864{\textendash}7875}, abstract = {

\<\>Abstract Image

Adsorption and transport in hierarchical porous solids with micro- (\∼1 nm) and mesoporosities (\>2 nm) are investigated by molecular simulation. Two models of hierarchical solids are considered: microporous materials in which mesopores are carved out (model A) and mesoporous materials in which microporous nanoparticles are inserted (model B). Adsorption isotherms for model A can be described as a linear combination of the adsorption isotherms for pure mesoporous and microporous solids. In contrast, adsorption in model B departs from adsorption in pure microporous and mesoporous solids; the inserted microporous particles act as defects, which help nucleate the liquid phase within the mesopore and shift capillary condensation toward lower pressures. As far as transport under a pressure gradient is concerned, the flux in hierarchical materials consisting of microporous solids in which mesopores are carved out obeys the Navier\–Stokes equation so that Darcy\’s law is verified within the mesopore. Moreover, the flow in such materials is larger than in a single mesopore, due to the transfer between micropores and mesopores. This nonzero velocity at the mesopore surface implies that transport in such hierarchical materials involves slippage at the mesopore surface, although the adsorbate has a strong affinity for the surface. In contrast to model A, flux in model B is smaller than in a single mesopore, as the nanoparticles act as constrictions that hinder transport. By a subtle effect arising from fast transport in the mesopores, the presence of mesopores increases the number of molecules in the microporosity in hierarchical materials and, hence, decreases the flow in the micropores (due to mass conservation). As a result, we do not observe faster diffusion in the micropores of hierarchical materials upon flow but slower diffusion, which increases the contact time between the adsorbate and the surface of the microporosity.

}, doi = {10.1021/la401228k}, author = {Benoit A. Coasne and Anne Galarneau and Gerardin, Corine and Fajula, Fran{\c c}ois and Fran{\c c}ois Villemot} } @article {45, title = {Molecular simulations of supercritical fluid permeation through disordered microporous carbons}, journal = {Langmuir}, volume = {29}, year = {2013}, month = {Aug-13-2013}, pages = {9985{\textendash}9990}, abstract = {
Abstract Image

Fluid transport through microporous carbon-based materials is inherent in numerous applications, ranging from gas separation by carbon molecular sieves to natural gas production from coal seams and gas shales. The present study investigates the steady-state permeation of supercritical methane in response to a constant cross-membrane pressure drop. We performed dual control volume grand canonical molecular dynamics (DCV-GCMD) simulations to mimic the conditions of actual permeation experiments. To overcome arbitrary assumptions regarding the investigated porous structures, the membranes were modeled after the CS1000a and CS1000 molecular models, which are representative of real microporous carbon materials. When adsorption-induced molecular trapping (AIMT) mechanisms are negligible, we show that the permeability of the microporous material, although not significantly sensitive to the pressure gradient, monotonically decreases with temperature and reservoir pressures, consistent with diffusion theory. However, when AIMT occurs, the permeability increases with temperature in agreement with experimental data found in the literature.

}, doi = {10.1021/la402087r}, author = {Botan, Alexandru and Vermorel, Romain and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {31, title = {Nanoindentation investigation of creep properties of calcium silicate hydrates}, journal = {Cement and Concrete Research}, volume = {52}, year = {2013}, month = {Oct-2013}, pages = {38{\textendash}52}, abstract = {

The creep properties of calcium silicate hydrates (C-S-H) are assessed by means of nanoindentation creep experiments on a wide range of substoichiometric cement pastes. We observe that, after a few seconds, the measured creep compliance of C-S-H is very well captured by a logarithmic time function. The rate of the logarithmic creep is found to scale in a unique manner with indentation modulus, indentation hardness, and packing density, independent of processing, mix proportions, indenter geometry and load history. The comparison with macroscopic creep experiments on concrete shows that minutes-long nanoindentations enable a quantitative assessment of the long-term creep properties of cementitious materials, orders of magnitude faster than macroscopic testing. Finally, we show that a strong analogy exists between this logarithmic creep behavior of C-S-H and that of soils, which suggests a granular origin of creep of geomaterials.

Typical measurement of variation Δh(t) in penetration depth during the creep ...

Coefficient γ governing the linear part of the fitting function (16) versus (a) ...

}, doi = {10.1016/j.cemconres.2013.05.006}, author = {Matthieu Vandamme and Franz-Josef Ulm} } @article {94, title = {Comparison of computational water models for simulation of calcium-silicate-hydrate}, journal = {Computational Materials Science}, volume = {53}, year = {2012}, month = {Feb-2012}, pages = {234-240}, type = {Article}, abstract = {

Calcium silicate hydrate, or C-S-H, is the chief hydration product of Portland cement. The structure of the C-S-H phase within cement has been proposed and developed via molecular simulations. In such simulations, empirical interatomic potentials for water molecules within C-S-H are adopted to govern the position and relative motion of this key constituent. Initial simulations and force fields of C-S-H have assumed the simplest molecular model of H2O termed {\textquoteleft}{\textquoteleft}single point charge{\&$\#$39;\&$\#$39;} or SPC, but this choice has not been validated by comparison with other computational models of water that confer additional bond flexibility or charge distribution. To enable efficiently computational modeling of C-S-H and to explore the role that H2O plays in maintaining C-S-H structure and properties, the choice of an efficient and accurate water model is critical. Here, we consider five distinct, classical atomistic water models (SPC, TIP3P, TIP4P, TIP4P05, and TIP5P) to determine the effects of these computational simplifications on C-S-H properties. Quantitative comparison of all five water models shows that the appropriate water model depends on the C-S-H characteristics of interest. Among these models, both SPC and TIP5P models successfully predict key properties of the structure and elastic constants of C-S-H, as well as the dynamics of water molecules within C-S-H. (C) 2011 Elsevier B. V. All rights reserved.

Schematics of (a) five water models and (b) calcium{\textendash}silicate{\textendash}hydrate (C{\textendash}S{\textendash}H) ...

Comparison of C{\textendash}S{\textendash}H structural characteristics. (a) Unit cell dimensions, where ...

Comparison of C{\textendash}S{\textendash}H dynamic characteristics. (a) Infrared absorption spectrum ...

Comparison of C{\textendash}S{\textendash}H elastic properties. (a) Elastic constants and the relative ...

}, keywords = {C-S-H, Cement, Simulation, Water model}, issn = {0927-0256}, doi = {10.1016/j.commatsci.2011.08.024}, author = {Qing Ji and Roland Jean-Marc Pellenq and Krystyn J. Van Vliet} } @article {57, title = {Confined water dissociation in microporous defective silicates: Mechanism, dipole distribution, and impact on substrate properties}, journal = {Journal of the American Chemical Society}, volume = {134}, year = {2012}, month = {Feb-01-2012}, pages = {2208{\textendash}2215}, abstract = {
Abstract Image

Interest in microporous materials has risen in recent years, as they offer a confined environment that is optimal to enhance chemical reactions. Calcium silicate hydrate (C-S-H) gel, the main component of cement, presents a layered structure with sub-nanometer-size disordered pores filled with water and cations. The size of the pores and the hydrophilicity of the environment make C-S-H gel an excellent system to study the possibility of confined water reactions. To investigate it, we have performed molecular dynamics simulations using the ReaxFF force field. The results show that water does dissociate to form hydroxyl groups. We have analyzed the water dissociation mechanism, as well as the changes in the structure and water affinity of the C-S-H matrix and water polarization, comparing the results with the behavior of water in a defective zeolite. Finally, we establish a relationship between water dissociation in C-S-H gel and the increase of hardness due to a transformation from a two- to a three-dimensional structure.

}, doi = {10.1021/ja209152n}, author = {Hegoi Manzano and Moeini, Sina and Marinelli, Francis and Adri CT Van Duin and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {51, title = {Hydration of calcium oxide surface predicted by reactive force field molecular dynamics}, journal = {Langmuir}, volume = {28}, year = {2012}, month = {Mar-06-2012}, pages = {4187{\textendash}4197}, abstract = {
Abstract Image

In this work, we present the parametrization of Ca\–O/H interactions within the reactive force field ReaxFF, and its application to study the hydration of calcium oxide surface. The force field has been fitted using density functional theory calculations on gas phase calcium\–water clusters, calcium oxide bulk and surface properties, calcium hydroxide, bcc and fcc Ca, and proton transfer reactions in the presence of calcium. Then, the reactive force field has been used to study the hydration of the calcium oxide {001} surface with different water contents. Calcium oxide is used as a catalyzer in many applications such as CO2 sequestration and biodiesel production, and the degree of surface hydroxylation is a key factor in its catalytic performance. The results show that the water dissociates very fast on CaO {001} bare surfaces without any defect or vacancy. The surface structure is maintained up to a certain amount of water, after which the surface undergoes a structural rearrangement, becoming a disordered calcium hydroxyl layer. This transformation is the most probable reason for the CaO catalytic activity decrease.

}, doi = {10.1021/la204338m}, author = {Hegoi Manzano and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Markus J Buehler and Adri CT Van Duin} } @article {307, title = {Particle shape dependence in 2D granular media}, journal = {EPL (Europhysics Letters)}, volume = {98}, year = {2012}, month = {Apr-20-2012}, pages = {Article Number: 44008}, abstract = {

Particle shape is a key to the space-filling and strength properties of granular matter. We consider a shape parameter

describing the degree of distortion from a perfectly spherical shape. Encompassing most specific shape characteristics such as elongation, angularity and non-convexity, is a low-order but generic parameter that we used in a numerical benchmark test for a systematic investigation of shape dependence in sheared granular packings composed of particles of different shapes. We find that the shear strength is an increasing function of with nearly the same trend for all shapes, the differences appearing thus to be of second order compared to . We also observe a non-trivial behavior of packing fraction which, for all our simulated shapes, increases with from the random close packing fraction for disks, reaches a peak considerably higher than that for disks, and subsequently declines as is further increased. These findings suggest that a low-order description of particle shape accounts for the principal trends of packing fraction and shear strength. Hence, the effect of second-order shape parameters may be investigated by considering different shapes at the same level of .

}, issn = {0295-5075}, doi = {10.1209/0295-5075/98/44008}, author = {Baptiste Saint-Cyr and Krzysztof Szarf and Charles Voivret and Emilien Az{\'e}ma and Vincent Richefeu and Jean-Yves Delenne and Gael Combe and C{\'e}cile Nouguier-Lehon and Pascal Villard and Philippe Sornay and Marie Chaze and Farhang Radja{\"\i}} } @article {56, title = {Poromechanics of microporous media}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {60}, year = {2012}, month = {Apr-2012}, pages = {606{\textendash}622}, abstract = {

Microporous media, i.e., porous media made of pores with a nanometer size, are important for a variety of applications, for instance for sequestration of carbon dioxide in coal, or for storage of hydrogen in metal-organic frameworks. In a pore of nanometer size, fluid molecules are not in their bulk state anymore since they interact with the atoms of the solid: they are said to be in an adsorbed state. For such microporous media, conventional poromechanics breaks down.

In this work we derive poroelastic constitutive equations which are valid for a generic porous medium, i.e., even for a porous medium with pores of nanometer size. The complete determination of the poromechanical behavior of a microporous medium requires knowing how the amount of fluid adsorbed depends on both the fluid bulk pressure and the strain of the medium. The derived constitutive equations are validated with the help of molecular simulations on one-dimensional microporous media. Even when a microporous medium behaves linearly in the absence of any fluid (i.e., its bulk modulus does not depend on strain), we show that fluid adsorption can induce non-linear behavior (i.e., its drained bulk modulus can then depend significantly on strain). We also show that adsorption can lead to an apparent Biot coefficient of the microporous medium greater than unity or smaller than zero.

The poromechanical response of a microporous medium to adsorption significantly depends on the pore size distribution. Indeed, the commensurability (i.e., the ratio of the size of the pores to that of the fluid molecules) proves to play a major role. For a one-dimensional model of micropores with a variety of pore sizes, molecular simulations show that the amount of adsorbed fluid depends linearly on the strain of the medium. We derive linearized constitutive equations which are valid when such a linear dependence of the adsorbed amount of fluid on the strain is observed.

As an application, the case of methane and coal is considered. Molecular simulations of an adsorption of methane on a microporous realistic model for coal are performed with a flexible solid skeleton. The applicability of the set of linearized constitutive equations to this case is discussed and the results are shown to be consistent with swelling data measured during a classical adsorption experiment.

Example of molecular configuration. The fluid molecules can enter or leave the ...

Amount of fluid molecules adsorbed per unit length of the undeformed chain, in ...

Apparent tangent drained bulk modulus Kt of the crystalline chain.

}, doi = {10.1016/j.jmps.2012.01.001}, author = {Brochard, Laurent and Matthieu Vandamme and Roland Jean-Marc Pellenq} } @article {260, title = {Quantifying plasticity-independent creep compliance and relaxation of viscoelastoplastic materials under contact loading}, journal = {Journal of Materials Research}, volume = {27}, year = {2012}, month = {Jan-2012}, pages = {302 - 312}, abstract = {

Here we quantify the time-dependent mechanical properties of a linear viscoelastoplastic material under contact loading. For contact load relaxation, we showed that the relaxation modulus can be measured independently of concurrent plasticity exhibited during the loading phase. For indentation creep, we showed that the rate of change of the contact creep compliance can be measured independently of any plastic deformation exhibited during loading through , where a(t) is the contact radius, h(t) is the displacement of the contact probe, and Pmax is the constant applied load during the creep phase. These analytical relations were compared with numerical simulations of conical indentation creep for a viscoelastoplastic material and validated against sharp indentation creep experiments conducted on polystyrene. The derived relations enable extraction of viscoelastic material characteristics, even if sharp probes confer concurrent plasticity, applicable for a general axisymmetric contact probe geometry and a general time-independent plasticity.

}, issn = {0884-2914}, doi = {10.1557/jmr.2011.302}, author = {Matthieu Vandamme and Tweedie, Catherine A. and Constantinides, Georgios and Franz-Josef Ulm and Krystyn J. Van Vliet} } @article {53, title = {Set in stone? A perspective on the concrete sustainability challenge}, journal = {MRS bulletin}, volume = {37}, year = {2012}, month = {Apr-09-2012}, pages = {395{\textendash}402}, abstract = {

As the most abundant engineered material on Earth, concrete is essential to the physical infrastructure of all modern societies. There are no known materials that can replace concrete in terms of cost and availability. There are, however, environmental concerns, including the significant CO2 emissions associated with cement production, which create new incentives for university\–industry collaboration to address concrete sustainability. Herein, we examine one aspect of this challenge\—the translation of scientific understanding at the microscale into industrial innovation at the macroscale\—by seeking improvements in cement-paste processing, performance, and sustainability through control of the mechanisms that govern microstructure development. Specifically, we consider modeling, simulation, and experimental advances in fracture, dissolution, precipitation, and hydration of cement paste precursors, as well as properties of the hardened cement paste within concrete. The aim of such studies is to optimize the chemical reactivity, mechanical performance, and other physical properties of cement paste to enable more sustainable processing routes for this ubiquitous material.

}, doi = {10.1557/mrs.2012.55}, author = {Krystyn J. Van Vliet and Roland Jean-Marc Pellenq and Markus J Buehler and Jeffrey C Grossman and Hamlin M. Jennings and Franz-Josef Ulm and Sidney Yip} } @article {93, title = {Thermodynamics of Water Confined in Porous Calcium-Silicate-Hydrates}, journal = {Langmuir}, volume = {28}, year = {2012}, month = {Aug-07-2012}, pages = {11422-11432}, type = {Article}, abstract = {
Abstract Image

Water within pores of cementitious materials plays a crucial role in the damage processes of cement pastes, particularly in the binding material comprising calcium-silicate-hydrates (C\–S\–H). Here, we employed Grand Canonical Monte Carlo simulations to investigate the properties of water confined at ambient temperature within and between C\–S\–H nanoparticles or \“grains\” as a function of the relative humidity (\%RH). We address the effect of water on the cohesion of cement pastes by computing fluid internal pressures within and between grains as a function of \%RH and intergranular separation distance, from 1 to 10 \Å. We found that, within a C\–S\–H grain and between C\–S\–H grains, pores are completely filled with water for \%RH larger than 20\%. While the cohesion of the cement paste is mainly driven by the calcium ions in the C\–S\–H, water facilitates a disjoining behavior inside a C\–S\–H grain. Between C\–S\–H grains, confined water diminishes or enhances the cohesion of the material depending on the intergranular distance. At very low \%RH, the loss of water increases the cohesion within a C\–S\–H grain and reduces the cohesion between C\–S\–H grains. These findings provide insights into the behavior of C\–S\–H in dry or high-temperature environments, with a loss of cohesion between C\–S\–H grains due to the loss of water content. Such quantification provides the necessary baseline to understand cement paste damaging upon extreme thermal, mechanical, and salt-rich environments.

}, issn = {0743-7463}, doi = {10.1021/la301738p}, author = {Patrick A Bonnaud and Qing Ji and Benoit A. Coasne and Roland Jean-Marc Pellenq and Krystyn J. Van Vliet} } @article {59, title = {Adsorption-induced deformation of microporous materials: coal swelling induced by co2{\textendash}ch4 competitive adsorption}, journal = {Langmuir}, volume = {28}, year = {2011}, month = {Dec-20-2011}, pages = {2659{\textendash}2670}, abstract = {
Abstract Image

Carbon dioxide injection in coal seams is known to improve the methane production of the coal seam, while ensuring a safe and long-term carbon sequestration. This improvement is due to the preferential adsorption of CO2 in coal with respect to CH4: an injection of CO2 thus results in a desorption of CH4. However, this preferential adsorption is also known to cause a differential swelling of coal, which results in a significant decrease in the reservoir permeability during the injection process. Recent studies have shown that adsorption in coal micropores (few angstr\öms in size) is the main cause of the swelling. In this work, we focus on the competitive adsorption behavior of CO2 and CH4 in micropores. We perform molecular simulations of adsorption with a realistic atomistic model for coal. The competitive adsorption is studied at various temperatures and pressures representative of those in geological reservoirs. With the help of a poromechanical model, we then quantify the subsequent differential swelling induced by the computed adsorption behaviors. The differential swelling is almost insensitive to the geological temperatures and pressures considered here and is proportional to the CO2 mole fraction in the coal.

}, doi = {10.1021/la204072d}, author = {Brochard, Laurent and Matthieu Vandamme and Roland Jean-Marc Pellenq and Fen-Chong, Teddy} } @article {315, title = {Rheology of granular materials composed of nonconvex particles}, journal = {Physical Review E}, volume = {84}, year = {2011}, month = {Oct-10-2011}, pages = {Article Number: 041302 Part: 1 }, abstract = {

By means of contact dynamics simulations, we investigate the shear strength and internal structure of granular materials composed of two-dimensional nonconvex aggregates. We find that the packing fraction first grows as the nonconvexity is increased but declines at higher nonconvexity. This unmonotonic dependence reflects the competing effects of pore size reduction between convex borders of aggregates and gain in porosity at the nonconvex borders that are captured in a simple model fitting nicely the simulation data both in the isotropic and sheared packings. On the other hand, the internal angle of friction increases linearly with nonconvexity and saturates to a value independent of nonconvexity. We show that fabric anisotropy, force anisotropy, and friction mobilization, all enhanced by multiple contacts between aggregates, govern the observed increase of shear strength and its saturation with increasing nonconvexity. The main effect of interlocking is to dislocate frictional dissipation from the locked double and triple contacts between aggregates to the simple contacts between clusters of aggregates. This self-organization of particle motions allows the packing to keep a constant shear strength at high nonconvexity.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, issn = {1539-3755}, doi = {10.1103/PhysRevE.84.041302}, author = {Baptiste Saint-Cyr and Jean-Yves Delenne and Charles Voivret and Farhang Radja{\"\i} and Philippe Sornay} }