@article {679, title = {Multiscale poromechanics of wet cement paste}, journal = {Proceedings of the National Academy of Sciences}, volume = {2}, year = {2019}, month = {May-09-2019}, pages = {201901160}, abstract = {

Capillary effects, such as imbibition drying cycles, impact the mechanics of granular systems over time. A multiscale poromechanics framework was applied to cement paste, which is the most common building material, experiencing broad humidity variations over the lifetime of infrastructure. First, the liquid density distribution at intermediate to high relative humidity is obtained using a lattice gas density functional method together with a realistic nanogranular model of cement hydrates. The calculated adsorption/desorption isotherms and pore size distributions are discussed and compare well with nitrogen and water experiments. The standard method for pore size distribution determination from desorption data is evaluated. Second, the integration of the Korteweg liquid stress field around each cement hydrate particle provided the capillary forces at the nanoscale. The cement mesoscale structure was relaxed under the action of the capillary forces. Local irreversible deformations of the cement nanograins assembly were identified due to liquid\–solid interactions. The spatial correlations of the nonaffine displacements extend to a few tens of nanometers. Third, the Love\–Weber method provided the homogenized liquid stress at the micrometer scale. The homogenization length coincided with the spatial correlation length of nonaffine displacements. Our results on the solid response to capillary stress field suggest that the micrometer-scale texture is not affected by mild drying, while nanoscale irreversible deformations still occur. These results pave the way for understanding capillary phenomena-induced stresses in heterogeneous porous media ranging from construction materials to hydrogels and living systems.

}, issn = {0027-8424}, doi = {10.1073/pnas.1901160116}, url = {http://www.pnas.org/lookup/doi/10.1073/pnas.1901160116}, author = {Zhou, Tingtao and Katerina Ioannidou and Franz-Josef Ulm and Bazant, Martin Z. and Roland Jean-Marc Pellenq} } @article {534, title = {Atomistic and mesoscale simulation of sodium and potassium adsorption in cement paste}, journal = {The Journal of Chemical Physics}, volume = {14963}, year = {2018}, month = {Aug-21-2018}, pages = {074705}, abstract = {

An atomistic and mesoscopic assessment of the effect of alkali uptake in cement paste is performed. Semi-grand canonical Monte Carlo simulations indicate that Na and K not only adsorb at the pore surface of calcium silicate hydrates (C-S-H) but also adsorb in the C-S-H hydrated interlayer up to concentrations of the order of 0.05 and 0.1 mol/kg, respectively. Sorption of alkali is favored as the Ca/Si ratio of C-S-H is reduced. Long timescale simulations using the Activation Relaxation Technique indicate that characteristic diffusion times of Na and K in the C-S-H interlayer are of the order of a few hours. At the level of individual grains, Na and K adsorption leads to a reduction of roughly 5\% of the elastic moduli and to volume expansion of about 0.25\%. Simulations using the so-called primitive model indicate that adsorption of alkali ions at the pore surface can reduce the binding between C-S-H grains by up to 6\%. Using a mesoscopic model of cement paste, the combination of individual grain swelling and changes in inter-granular cohesion was estimated to lead to overall expansive pressures of up to 4 MPa\—and typically of less than 1 MPa\—for typical alkali concentrations observed at the proximity of gel veins caused by the alkali-silica reaction.

PDF FULL TEXT

}, issn = {0021-9606}, doi = {10.1063/1.5042755}, url = {https://aip.scitation.org/doi/10.1063/1.5042755}, author = {Dufresne, Alice and Arayro, Jack and Zhou, Tingtao and Katerina Ioannidou and Franz-Josef Ulm and Roland Jean-Marc Pellenq and B{\'e}land, Laurent Karim} } @article {571, title = {Griffith{\textquoteright}s postulate: Grand Canonical Monte Carlo approach for fracture mechanics of solids}, journal = {Engineering Fracture Mechanics}, volume = {199}, year = {2018}, month = {Aug-2018}, pages = {544 - 554}, abstract = {

A Grand Canonical Monte Carlo Approach (GCMC) is proposed for the fracture analysis of solids discretized as mass points and bond interactions. In contrast to classical load-driven fracture processes, the GCMC approach introduces an auxiliary field, the bond rupture potential , to which the system is subjected; in addition to changes in volume V and temperature T. In this VT-ensemble, bond isotherms that link the average number of bonds to the bond rupture potential (Nk-) are obtained that carry critical information for fracture analysis. Specifically, the slope of the bond isotherm reflects bond fluctuations, permitting identification of (1) a fluctuation-based damage variable, and (2) the competition in energy fluctuations between the redistribution of strain energy induced by bond rupture, and the dissipation of the groundstate energy. Based on these fluctuations, it is shown that the GCMC-approach allows the identification of a critical bond energy release rate of material samples, when strain energy fluctuations equal groundstate energy fluctuations \– much akin to Griffith\’s 1921 stationarity postulate to \“predicting the breaking loads of elastic solids\”. This is illustrated by means of thermodynamic integration of bond isotherms to determine force-displacement curves, for both notched and unnotched homogeneous samples discretized by regular 2-D lattices with bonds exhibiting harmonic potentials.

Fig. 1. Bond-Isotherms Nk-μfor three systems: (a) unnotched; (b) notched; (c) broken

Fig. 2. Fluctuation-based damage in function of the prescribed bond rupture potential,{\textellipsis}

Fig. 3. Energy dissipation due to fluctuations for (a) the unnotched and (b) the{\textellipsis}

Fig. 4. Thermodynamic integration of the bond isotherms of (a) an unnotched sample; and{\textellipsis}

}, keywords = {Grand Canonical Monte Carlo Simulations; Bond rupture; Bond rupture potential; Bond isotherms; Damage; Fluctuations}, issn = {00137944}, doi = {10.1016/j.engfracmech.2018.06.001}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0013794418301747}, author = {Al-Mulla, Talal and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {576, title = {Impact of Nanoporosity on Hydrocarbon Transport in Shales{\textquoteright} Organic Matter}, journal = {Nano Letters}, volume = {18}, year = {2018}, month = {Jan-16-2018}, pages = {832 - 837}, abstract = {

In a context of growing attention for shale gas, the precise impact of organic matter (kerogen) on hydrocarbon recovery from unconventional reservoirs still has to be assessed. Kerogen\&$\#$39;s microstructure is characterized by a very disordered pore network that greatly affects hydrocarbon transport. The specific structure and texture of this organic matter at the nanoscale is highly dependent on its origin. In this study, by the use of statistical physics and molecular dynamics, we shed some new lights on hydrocarbon transport through realistic molecular models of kerogen at different level of maturity [Bousige et al. Nat. Mater. 2016, 15, 576]. Despite the apparent complexity, severe confinement effects controlled by the porosity of the various kerogens allow linear alkanes (from methane to dodecane) transport to be studied only via the self-diffusion coefficients of the species. The decrease of the transport coefficients with the amount of adsorbed fluid can be described by a free volume theory. Ultimately, the transport coefficients of hydrocarbons can be expressed simply as a function of the porosity (volume fraction of void) of the microstructure, thus paving the way for shale gas recovery predictions.

https://pubs-acs-org.libproxy.mit.edu/appl/literatum/publisher/achs/journals/content/nalefd/2018/nalefd.2018.18.issue-2/acs.nanolett.7b04079/20180208/images/medium/nl-2017-04079p_0004.gif

}, keywords = {Kerogen; microporosity; diffusion; hydrocarbons; shales}, issn = {1530-6984}, doi = {10.1021/acs.nanolett.7b04079}, url = {http://pubs.acs.org/doi/10.1021/acs.nanolett.7b04079}, author = {Ama{\"e}l Obliger and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {580, title = {Improving the practicality and safety of artificial corneas: Pre-assembly and gamma-rays sterilization of the Boston Keratoprosthesis}, journal = {The Ocular Surface}, volume = {16}, year = {2018}, month = {Jul-2018}, pages = {322 - 330}, abstract = {

Purpose: To make the Boston keratoprosthesis (B-KPro), together with its carrier corneal graft, more easily procured, transported and stored, as well as less expensive, easier for the surgeon to implant and safer for the patient, it is proposed that the B-KPro-graft combination be pre-assembled by an expert technician, followed by sterilization with gamma ray irradiation (GI) allowing long-term storage at room temperature. For this to be possible, it must be shown that the B-KPro itself (not only the graft) remains unharmed by the irradiation.

Methods: Polymethyl methacrylate (PMMA) discs and B-KPros were submitted to either ethylene oxide sterilization or different doses of GI. Cell biocompatibility, mechanical strength and optical quality were evaluated. The feasibility of assembling the B-KPro to a corneal graft, and gamma-radiate afterwards, was also assessed.

Results: There were no differences in cell biocompatibility between the samples. The optical evaluation showed high levels of transparency for all the groups. The absorbance of ultraviolet was higher for the groups treated with GI. The mechanical evaluation by nanoindentation showed no alterations of the PMMA discs after GI. The flexure test revealed a similar mechanical behavior. Technically, pre-assembly and GI of the B-KPro revealed no problems.

Conclusions: Sterilization of B-KPro using GI has no detrimental influence on the device. The pre-assembly of B-KPro to a donor cornea, followed by gamma sterilization, emerges as an efficient and safe procedure. (C) 2018 Elsevier Inc. All rights reserved.

}, keywords = {Artificial corneas; Boston keratoprosthesis; Gamma radiation; Gamma rays; Preassembly; Sterilization; Polymethyl methacrylate}, issn = {15420124}, doi = {10.1016/j.jtos.2018.04.002}, url = {https://linkinghub.elsevier.com/retrieve/pii/S1542012418300442}, author = {Gonzalez-Andrades, Miguel and Sharifi, Roholah and Islam, Mohammad Mirazul and Divoux, Thibaut and Haist, Michael and Paschalis, Eleftherios I. and Gelfand, Larisa and Mamodaly, Shamina and Di Cecilia, Luca and Cruzat, Andrea and Franz-Josef Ulm and Chodosh, James and Delori, Francois and Dohlman, Claes H.} } @article {575, title = {Le Ch{\^a}telier{\textquoteright}s conjecture: Measurement of colloidal eigenstresses in chemically reactive materials}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {112}, year = {2018}, month = {Mar-2018}, pages = {334 - 344}, abstract = {

Volume changes in chemically reactive materials, such as hydrating cement, play a critical role in many engineering applications that require precise estimates of stress and pressure developments. But a means to determine bulk volume changes in the absence of other deformation mechanisms related to thermal, pressure and load variations, is still missing. Herein, we present such a measuring devise, and a hybrid experimental-theoretical technique that permits the determination of colloidal eigenstresses. Applied to cementitious materials, it is found that bulk volume changes in saturated cement pastes at constant pressure and temperature conditions result from a competition of repulsive and attractive phenomena that originate from the relative distance of the solid particles- much as Henry Louis Le Chatelier, the father of modern cement science, had conjectured in the late 19th century. Precipitation of hydration products in confined spaces entails a repulsion, whereas the concurrent reduction in interparticle distance entails activation of attractive forces in charged colloidal particles. This cross-over from repulsion to attraction can be viewed as a phase transition between a liquid state (below the solid percolation) and the limit packing of hard spheres, separated by an energy barrier that defines the temperature-dependent eigenstress magnitude. (C) 2017 Elsevier Ltd. All rights reserved.

Fig. 1. Pressure vessel measurement device for volume changes under constant pressure{\textellipsis}

Fig. 2. Bulk Volume Change measurement results under drained and fluid saturated{\textellipsis}

Fig. 3. Eigenstress Scaling: (a) Normalized eigenstress (|σ*|/max|σ*|) vs

Fig. 4. Eigenstress interpretation: (a) volume fractions vs

}, keywords = {Poromechanics; Colloidal stresses; Cement-based materials; Expansion; Shrinkage}, issn = {00225096}, doi = {10.1016/j.jmps.2017.12.012}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0022509617309985}, author = {Muhannad Abuhaikal and Katerina Ioannidou and Petersen, Thomas and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @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 {658, title = {A methodology to calibrate and to validate effective solid potentials of heterogeneous porous media from computed tomography scans and laboratory-measured nanoindentation data}, journal = {Acta Geotechnica}, volume = {13}, year = {2018}, month = {Jan-12-2018}, pages = {1369 - 1394}, abstract = {

Built on the framework of effective interaction potentials using lattice element method, a methodology to calibrate and to validate the elasticity of solid constituents in heterogeneous porous media from experimentally measured nanoindentation moduli and imported scans from advanced imaging techniques is presented. Applied to computed tomography (CT) scans of two organic-rich shales, spatial variations of effective interaction potentials prove instrumental in capturing the effective elastic behavior of highly heterogeneous materials via the first two cumulants of experimentally measured distributions of nanoindentation moduli. After calibration and validation steps while implicitly accounting for mesoscale texture effects via CT scans, Biot poroelastic coefficients are simulated. Analysis of stress percolation suggests contrasting pathways for load transmission, a reflection of microtextural differences in the studied cases. This methodology to calibrate elastic energy content of real materials from advanced imaging techniques and experimental measurements paves the way to study other phenomena such as wave propagation and fracture while providing a platform to fine-tune effective behavior of materials given advancements in additive manufacturing and machine learning algorithms .

}, issn = {1861-1125}, doi = {10.1007/s11440-018-0687-9}, url = {http://link.springer.com/10.1007/s11440-018-0687-9}, author = {Monfared, Siavash and Hadrien Laubie and Farhang Radja{\"\i} and Hubler, Mija H. and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {610, title = {A minimal power model for human running performance}, journal = {PLOS ONE}, volume = {13}, year = {2018}, month = {Nov-16-2018}, pages = {e0206645}, abstract = {

Models for human running performances of various complexities and underlying principles have been proposed, often combining data from world record performances and bio-energetic facts of human physiology. The purpose of this work is to develop a novel, minimal and universal model for human running performance that employs a relative metabolic power scale. The main component is a self-consistency relation for the time dependent maximal power output. The analytic approach presented here is the first to derive the observed logarithmic scaling between world (and other) record running speeds and times from basic principles of metabolic power supply. Our hypothesis is that various female and male record performances (world, national) and also personal best performances of individual runners for distances from 800m to the marathon are excellently described by this model. Indeed, we confirm this hypothesis with mean errors of (often much) less than 1\%. The model defines endurance in a way that demonstrates symmetry between long and short racing events that are separated by a characteristic time scale comparable to the time over which a runner can sustain maximal oxygen uptake. As an application of our model, we derive personalized characteristic race speeds for different durations and distances.

}, keywords = {ENERGY-COST; OXYGEN-UPTAKE; ATHLETIC RECORDS; ECONOMY; RUNNERS; RUN; ENDURANCE; MARATHON; EXERCISE; EXHAUSTION}, doi = {10.1371/journal.pone.0206645}, url = {http://dx.plos.org/10.1371/journal.pone.0206645}, author = {Mulligan, Matthew and Adam, Guillaume and Emig, Thorsten}, editor = {Piknova, Barbora} } @article {627, title = {Multiscale modeling for bioresources and bioproducts}, journal = {Innovative Food Science \& Emerging Technologies}, volume = {46}, year = {2018}, month = {Avr-2018}, pages = {41 - 53}, abstract = {

Designing and processing complex matter and materials are key objectives of bioresource and bioproduct research. Modeling approaches targeting such systems have to account for their two main sources of complexity: their intrinsic multi-scale nature; and the variability and heterogeneity inherent to all living systems. Here we provide insight into methods developed at the Food \& Bioproduct Engineering division (CEPIA) of the French National Institute of Agricultural Research (INRA). This brief survey focuses on innovative research lines that tackle complexity by mobilizing different approaches with complementary objectives. On one hand cognitive approaches aim to uncover the basic mechanisms and laws underlying the emerging collective properties and macroscopic behavior of soft-matter and granular systems, using numerical and experimental methods borrowed from physics and mechanics. The corresponding case studies are dedicated to the structuring and phase behavior of biopolymers, powders and granular materials, and to the evolution of these structures caused by external constraints. On the other hand machine learning approaches can deal with process optimizations and outcome predictions by extracting useful information and correlations from huge datasets built from experiments at different length scales and in heterogeneous conditions. These predictive methods are illustrated in the context of cheese ripening, grape maturity prediction and bacterial production.

}, issn = {14668564}, doi = {10.1016/j.ifset.2017.09.015}, url = {https://linkinghub.elsevier.com/retrieve/pii/S1466856417302230}, author = {Barnabe, M. and Blanc, Nicolas and Chabin, T. and Jean-Yves Delenne and Duri, A. and Frank, Xavier and Hugouvieux, V. and Lutton, E. and Mabille, F. and Saeid Nezamabadi and Perrot, N. and Farhang Radja{\"\i} and Ruiz, T. and Tonda, A.} } @article {532, title = {Phase separation of stable colloidal clusters}, journal = {Physical Review Materials}, volume = {2}, year = {2018}, month = {Sep-25-2018}, abstract = {

This article presents a nonequilibrium thermodynamic theory for the mean-field precipitation, aggregation, and pattern formation of colloidal clusters. A variable gradient energy coefficient and the arrest of particle diffusion upon \“jamming\” of cluster aggregates in the spinodal region predicts observable gel patterns that, at high intercluster attraction, form system-spanning, out-of-equilibrium networks with glasslike, quasistatic structural relaxation. For reactive systems, we incorporate the free energy landscape of stable prenucleation clusters into the Allen-Cahn reaction equation. We show that pattern formation is dominantly controlled by the Damk\öhler number and the stability of the clusters, which modifies the autocatalytic rate of precipitation. As clusters individually become more stable, bulk phase separation is suppressed.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, doi = {10.1103/PhysRevMaterials.2.095602}, url = {https://journals.aps.org/prmaterials/abstract/10.1103/PhysRevMaterials.2.095602}, author = {Petersen, Thomas and Bazant, Martin Z. and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @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 {572, title = {Potential-of-Mean-Force Approach for Molecular Dynamics{\textendash}Based Resilience Assessment of Structures}, journal = {Journal of Engineering Mechanics}, volume = {144}, year = {2018}, month = {Aug-2018}, pages = {04018066}, abstract = {

A molecular dynamics (MD)-based structural mechanics approach is proposed for the assessment of resilience of buildings. At the core of the approach, potentials of mean force (PMFs) suitable for structural members for both two-body (stretch) and three-body (bending) interactions are derived to define the energy states between mass points discretizing structural members. An original potential parameter calibration procedure is proposed: for close-to-equilibrium potential parameters, the procedure is based on matching measured frequency of a structure with the frequency of the molecular model. In turn, for bond-rupture parameters, it is shown that classical interatomic potential expressions, such as Morse potential, can be used to calibrate the energy content of many structural members and connections. By way of example, the MD-based structural mechanics approach is applied to a large-scale structure. Compared with classical continuum-based approaches, the added value of the method thus proposed is a rational means of determining the progressive structural collapse load of structures based on thermodynamic integration. By redefining structural mechanics within the context of statistical physics, molecular simulations, and potentials of mean force, the approach provides a powerful means of determining fragility curves required for the assessment of resilience of buildings.

}, keywords = {Molecular dynamics; Structural mechanics; Potential of mean force; Morse potential; Progressive structural collapse; Fragility curves}, issn = {0733-9399}, doi = {10.1061/(ASCE)EM.1943-7889.0001491}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29EM.1943-7889.0001491}, author = {Keremides, Konstantinos and Abdolhosseini Qomi, Mohammad Javad and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @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 {574, title = {Role of City Texture in Urban Heat Islands at Nighttime}, journal = {Physical Review Letters}, volume = {120}, year = {2018}, month = {Mar-09-2018}, pages = {Article Number: 108701 }, abstract = {

An urban heat island (UHI) is a climate phenomenon that results in an increased air temperature in cities when compared to their rural surroundings. In this Letter, the dependence of an UHI on urban geometry is studied. Multiyear urban-rural temperature differences and building footprints data combined with a heat radiation scaling model are used to demonstrate for more than 50 cities worldwide that city texture-measured by a building distribution function and the sky view factor-explains city-to-city variations in nocturnal UHIs. Our results show a strong correlation between nocturnal UHIs and the city texture.

Figure

Figure

Figure

Figure

}, keywords = {IDEAL CONDITIONS; CLIMATE-CHANGE; CITIES; SIMULATION; GEOMETRY; ENERGY; HEALTH; IMPACT; MODEL}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.120.108701}, url = {https://link.aps.org/doi/10.1103/PhysRevLett.120.108701}, author = {Sobstyl, Jake M. and Emig, Thorsten and Qomi, M. J. Abdolhosseini and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {652, title = {Synthesis, characterization, and catalytic application in aldehyde hydrosilylation of half-sandwich nickel complexes bearing (κ1-C)- and hemilabile (κ 2-C, S)-thioether-functionalised NHC ligands}, journal = {Dalton Transactions}, volume = {47}, year = {2018}, month = {Dec-21-2018}, pages = {17134-17145}, abstract = {

Neutral nickel-N-heterocyclic carbene complexes, ((1)-C)-[NiCpBr{R-NHC-(CH2)(2)SR}] [Cp = (5)-C5H5; R-NHC-(CH2)(2)SR = 1-mesityl-3-[2-(tert-butylthio)ethyl]- (1a), 1-mesityl-3-[2-(phenylthio)ethyl]- (1b), 1-benzyl-3-[2-(tert-butylthio)ethyl]- (1c), 1-benzyl-3-[2-(phenylthio)ethyl]-imidazol-2-ylidene (1d)], which bear a N-bound thioether side arm, were prepared by the reaction of nickelocene with the corresponding imidazolium bromides [R-NHC-(CH2)(2)SRHBr] (a-d), via conventional or microwave heating. The H-1 NMR spectra of the benzyl-substituted species 1c and 1d showed signals for diastereotopic NCH2CH2S protons at room temperature. However, structural studies established the absence of coordination of the sulphur atom in the solid state, and solvent DFT calculations showed that bromide displacement by sulphur is an unfavourable process (G = +13.5 kcal mol(-1) for 1d), thereby suggesting that the observed disatereotopicity is more likely due to significant steric congestion rather than to a possible C,S-chelation in solution. Treatment of these complexes with KPF6 in tetrahydrofuran (THF) led to bromide abstraction to afford the cationic complexes [NiCp{R-NHC-(CH2)(2)SR}](PF6) (2a-c). Alternatively, 2a-c could also be prepared by the direct reaction of nickelocene with the corresponding imidazolium hexafluorophosphate salts [R-NHC-(CH2)(2)SRHPF6]. Inversely to the neutral species, whereas X-ray crystallography established C,S-chelation in the solid state, the H-1 NMR spectra (CDCl3, CD2Cl2, or thf-d(8)) at room temperature showed no diastereotopic NCH2CH2S protons, thus suggesting the possible displacement of the sulphur atom by the respective solvents and/or very fast sulphur inversion. DFT calculations established a low energy inversion process in all cases (+9 G(double dagger) +13 kcal mol(-1)) as well as a favourable solvent coordination process (G(double dagger) approximate to +11 kcal mol(-1); G approximate to -7 kcal mol(-1)) with a solvent such as THF, thus suggesting that sulphur inversion and/or solvent coordination can both account for the absence of diastereotopy at room temperature, depending on the solvent. While all complexes catalysed the hydrosilylation of benzaldehyde in the absence of any additive, the cationic C,S-chelated complexes 2 proved more active than the sterically constrained neutral species 1. In particular, 2c proved to be the most active pre-catalyst and its catalytic charge could be lowered down to 2 mol\% with PhSiH3 as the hydrogen source.

}, issn = {1477-9226}, doi = {10.1039/C8DT03882A}, url = {http://xlink.rsc.org/?DOI=C8DT03882A}, author = {Franz-Josef Ulm and Poblador-Bahamonde, Amalia I. and Choppin, Sabine and Bellemin-Laponnaz, St{\'e}phane and Chetcuti, Michael J. and Achard, Thierry and Ritleng, Vincent} } @article {535, title = {Thermalizing and Damping in Structural Dynamics}, journal = {Journal of Applied Mechanics}, volume = {85}, year = {2018}, month = {Aug-2018}, abstract = {

Structural damping, that is the presence of a velocity dependent dissipative term in the equation of motion, is rationalized as a thermalization process between a structure (here a beam) and an outside bath (understood in a broad sense as a system property). This is achieved via the introduction of the kinetic temperature of structures and formalized by means of an extended Lagrangian formulation of a structure in contact with an outside bath at a given temperature. Using the Nose-Hoover thermostat, the heat exchange rate between structure and bath is identified as a mass damping coefficient, which evolves in time in function of the kinetic energy/temperature history exhibited by the structure. By way of application to a simple beam structure subjected to eigen-vibrations and dynamic buckling, commonality and differences of the Nose-Hoover beam theory with constant mass damping models are shown, which permit a handshake between classical damping models and statistical mechanics-based thermalization models. The solid foundation of these thermalization models in statistical physics provides new insights into stability and instability for engineering structures. Specifically, since two systems are considered in (thermodynamic) equilibrium when they have the same temperature, we show in the case of dynamic buckling that a persistent steady-state difference in kinetic temperature between structure and bath is but indicative of the instability of the system. This shows that the kinetic temperature can serve as a structural order parameter to identify and comprehend failure of structures, possibly well beyond the elastic stability considered here.

PDF FULL TEXT

jam_085_08_081001_f001.png

jam_085_08_081001_f002.png

jam_085_08_081001_f003.png

}, keywords = {structural dynamics; damping; Nose-Hoover Bath; kinetic temperature; dynamic buckling}, issn = {0021-8936}, doi = {10.1115/1.4040080}, url = {http://appliedmechanics.asmedigitalcollection.asme.org.libproxy.mit.edu/article.aspx?articleid=2680040}, author = {Arghavan Louhghalam and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {573, title = {Thermodynamics, kinetics, and mechanics of cesium sorption in cement paste: A multiscale assessment}, journal = {Physical Review Materials}, volume = {2}, year = {2018}, month = {May-30-2018}, chapter = {Article Number: 053608 }, abstract = {

Cesium-137 is a common radioactive byproduct found in nuclear spent fuel. Given its 30 year half life, its interactions with potential storage materials-such as cement paste-is of crucial importance. In this paper, simulations are used to establish the interaction of calcium silicate hydrates (C-S-H)-the mam binding phase of cement paste-with Cs at the nano- and mesoscale. Different C-S-H compositions are explored, including a range of Ca/Si ratios from 1.0 to 2.0. These calculations are based on a set of 150 atomistic models, which qualitatively and quantitatively reproduce a number of experimentally measured features of C-S-H-within limits intrinsic to the approximations imposed by classical molecular dynamics and the steps followed when building the models. A procedure where hydrated Ca2+ ions are swapped for Cs1+ ions shows that Cs adsorption in the C-S-H interlayer is preferred to Cs adsorption at the nanopore surface when Cs concentrations are lower than 0.19 Mol/kg. Interlayer sorption decreases as the Ca/Si ratio increases. The activation relaxation technique nouveau is used to access timescales out of the reach of traditional molecular dynamics (MD). It indicates that characteristic diffusion time for Cs1+ in the C-S-H interlayer is on the order of a few hours. Cs uptake in the interlayer has little impact on the elastic response of C-S-H. It leads to swelling of the C-S-H grams, but mesoscale calculations that access length scales out of the range of MD indicate that this leads to practically negligible expansive pressures for Cs concentrations relevant to nuclear waste repositories.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, keywords = {CALCIUM-SILICATE-HYDRATE; C-S-H; MOLECULAR-DYNAMICS; CRYSTAL-STRUCTURE; MONTE-CARLO; IMMOBILIZATION; RADIONUCLIDES; SIMULATIONS; RELAXATION; COLLOIDS}, doi = {10.1103/PhysRevMaterials.2.053608}, url = {https://link.aps.org/doi/10.1103/PhysRevMaterials.2.053608}, author = {Arayro, Jack and Dufresne, Alice and Zhou, Tingtao and Katerina Ioannidou and Franz-Josef Ulm and Roland Jean-Marc Pellenq and B{\'e}land, Laurent Karim} } @article {633, title = {Two-dimensional numerical simulation of chimney fluidization in a granular medium using a combination of discrete element and lattice Boltzmann methods}, journal = {Physical Review E}, volume = {97}, year = {2018}, month = {May-10-2018}, pages = { Article Number 052902 }, abstract = {

We present here a numerical study dedicated to the fluidization of a submerged granular medium induced by a localized fluid injection. To this end, a two-dimensional (2D) model is used, coupling the lattice Boltzmann method (LBM) with the discrete element method (DEM) for a relevant description of fluid-grains interaction. An extensive investigation has been carried out to analyze the respective influences of the different parameters of our configuration, both geometrical (bed height, grain diameter, injection width) and physical (fluid viscosity, buoyancy). Compared to previous experimental works, the same qualitative features are recovered as regards the general phenomenology including transitory phase, stationary states, and hysteretic behavior. We also present quantitative findings about transient fluidization, for which several dimensionless quantities and scaling laws are proposed, and about the influence of the injection width, from localized to homogeneous fluidization. Finally, the impact of the present 2D geometry is discussed, by comparison to the real three-dimensional (3D) experiments, as well as the crucial role of the prevailing hydrodynamic regime within the expanding cavity, quantified through a cavity Reynolds number, that can presumably explain some substantial differences observed regarding upward expansion process of the fluidized zone when the fluid viscosity is changed.

}, issn = {2470-0045}, doi = {10.1103/PhysRevE.97.052902}, url = {https://link.aps.org/doi/10.1103/PhysRevE.97.052902}, author = {Jeff Ngoma and Pierre Philippe and St{\'e}phane Bonelli and Farhang Radja{\"\i} and Jean-Yves Delenne} } @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 {244, title = {Crystal-chemistry control of the mechanical properties of 2:1 clay minerals}, journal = {Applied Clay Science}, volume = {143}, year = {2017}, month = {Jul-2017}, pages = {387 - 398}, abstract = {

Clay minerals are the main constituents of the clay matrix of a wide variety of sedimentary deposits. When subjected to burial, some of these minerals undergo phase transitions accompanied by atomic substitutions which ultimately impact the cohesive interactions between their constitutive layers. The most common among such transitions is smectite illitization, which is also highly relevant for oil and gas exploration and production from source rocks. The impact of this transition on the mechanical properties of clay minerals as well as clay bearing rocks remains, however, to be properly addressed. To this end, a set of macroscopic single 2:1 clay minerals (pyrophyllite, talc, vermiculite, phlogopite, muscovite, and clintonite) representative of the two octahedral fillings and the variation in surface charge densities was investigated. A hybrid experimental-modeling approach is proposed, which combines nanoindentation in orthogonal directions (x1 and x3-directions) with analytical derivations of the cohesive energy (Ui) using XRD and TEM-EDS measurements. The results highlight that the interlayer energy defines the elasticity (modulus M), strength (hardness H) and ductility behavior (M/H) of these materials; and not the bond energy stored into the constituent layers. This finding permits the derivation of predictive stiffness and strength functional relations as derivatives of the interlayer energy that account for the arrangement of nanoscale layers and the interlayer composition. These relations suggest that as the cohesion increases with the coulombic interactions through progressive smectite illitization, the system loses its capacity to dissipate applied energy by dislocation mechanisms in between the layers, entailing an increase in brittleness of the clay particles with burial. By way of conclusion, the geophysics and geochemistry implications of these results for predicting the macroscopic mechanical performance of clay bearing geomaterials, such as economically valuable source rocks, are discussed.

}, issn = {01691317}, doi = {10.1016/j.clay.2017.04.010}, author = {Berthonneau, Jeremie and Christian G. Hoover and Grauby, Olivier and Alain Baronnet and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {239, title = {Disorder-induced stiffness degradation of highly disordered porous materials}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {106}, year = {2017}, month = {Sep-2017}, chapter = {207-228}, abstract = {

The effective mechanical behavior of multiphase solid materials is generally modeled by means of homogenization techniques that account for phase volume fractions and elastic moduli without considering the spatial distribution of the different phases. By means of extensive numerical simulations of randomly generated porous materials using the lattice element method, the role of local textural properties on the effective elastic properties of disordered porous materials is investigated and compared with different continuum micromechanics-based models. It is found that the pronounced disorder-induced stiffness degradation originates from stress concentrations around pore clusters in highly disordered porous materials. We identify a single disorder parameter, \φsa, which combines a measure of the spatial disorder of pores (the clustering index, sa) with the pore volume fraction (the porosity, \φ) to scale the disorder-induced stiffness degradation. Thus, we conclude that the classical continuum micromechanics models with one spherical pore phase, due to their underlying homogeneity assumption fall short of addressing the clustering effect, unless additional texture information is introduced, e.g. in form of the shift of the percolation threshold with disorder, or other functional relations between volume fractions and spatial disorder; as illustrated herein for a differential scheme model representative of a two-phase (solid\–pore) composite model material.

Fig. 1. Two-dimensional porous media: (a) ordered system, (b) type 1 disorder (pores{\textellipsis}

Fig. 2. Two-point probability functions and associated microstructures of four systems{\textellipsis}

Fig. 3. Probability Density Function (PDF) of the local porosity (φa) and associated{\textellipsis}

Fig. 4. (a) Degrees of freedom of a link element between points i and j, (b) D3Q18 unit{\textellipsis}

Fig. 5. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of the porosity{\textellipsis}

Fig. 6. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of the porosity

Fig. 7. Ordered system under uniaxial strain, (a) schematic stress map, the darker the{\textellipsis}

Fig. 8. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of the porosity

Fig. 9. Probability Density Function (PDF) of the normalized stress (σxx/〈σxx〉) in the{\textellipsis}

Fig. 10. Dimensionless disordered Young{\textquoteright}s modulus: F(φsa)=Eeff(φ,sa)/Eordered(φ) as{\textellipsis}

Fig. 11. Two-phase periodic porous solid

Fig. 12. Dimensionless disordered Young{\textquoteright}s modulus: FMT(φsa)=Eeff(φ,sa)/EMT(φ) as{\textellipsis}

Fig. 13. Integration path for the lower bound differential scheme

Fig. 14. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of porosity the

Fig. 15. Dimensionless disordered Young{\textquoteright}s modulus: FMT(φsa)=Eeff(φ,sa)/EMT(φ) as{\textellipsis}

Fig. B.16. (a) Porosity as a function of the number of pores

Fig. D.17. High porosity limit geometry

Fig. E.18. Calculation of the local porosity in a square observation window

}, doi = {10.1016/j.jmps.2017.05.008}, author = {Hadrien Laubie and Monfared, Siavash and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {601, title = {Effect of particle size distribution on 3D packings of spherical particles}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 02030}, abstract = {

We use molecular dynamics simulations of frictionless spherical particles to investigate a class of polydisperse granular materials in which the particle size distribution is uniform in particle volumes. The particles are assembled in a box by uniaxial compaction under the action of a constant stress. Due to the absence of friction and the nature of size distribution, the generated packings have the highest packing fraction at a given size span, defined as the ratio of the largest size to the smallest size. We find that, up to = 5, the packing fraction is a nearly linear function of . While the coordination number is nearly constant due to the isostatic nature of the packings, we show that the connectivity of the particles evolves with . In particular, the proportion of particles with 4 contacts represents the largest proportion of particles mostly of small size. We argue that this particular class of particles occurs as a result of the high stability of local configurations in which a small particle is stuck by four larger particles.

}, doi = {10.1051/epjconf/201714002030}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714002030}, author = {Taiebat, Mahdi and Patrick Mutabaruka and Roland Jean-Marc Pellenq and Farhang Radja{\"\i}}, editor = {Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {591, title = {Effective Potentials and Elastic Properties in the Lattice-Element Method: Isotropy and Transverse Isotropy}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Sep-2017}, pages = { Article Number: UNSP 04017007 }, abstract = {

Lattice approaches have emerged as a powerful tool to capture the effective mechanical behavior of heterogeneous materials using harmonic interactions inspired from beam-type stretch and rotational interactions between a discrete number of mass points. In this paper, the lattice element method (LEM) is reformulated within the conceptual framework of empirical force fields employed at the lattice scale. Within this framework, because classical harmonic formulations are but a Taylor expansion of nonharmonic potential expressions, they can be used to model both the linear and the nonlinear response of discretized material systems. Specifically, closed-form calibration procedures for such interaction potentials are derived for both the isotropic and the transverse isotropic elastic cases on cubic lattices, in the form of linear relations between effective elasticity properties and energy parameters that define the interactions. The relevance of the approach is shown by an application to the classical Griffith crack problem. In particular, it is shown that continuum-scale quantities of linear-elastic fracture mechanics, such as stress intensity factors (SIFs), are well captured by the method, which by its very discrete nature removes geometric discontinuities that provoke stress singularities in the continuum case. With its strengths and limitations thus defined, the proposed LEM is well suited for the study of multiphase materials whose microtextural information is obtained by, e.g., X-ray micro-computed tomography. (c) 2017 American Society of Civil Engineers.

}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000125}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000125}, author = {Hadrien Laubie and Monfared, Siavash and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @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 {245, title = {Inhomogeneity in Cement Hydrates: Linking Local Packing to Local Pressure}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Mar-03-2017}, pages = {Article Number: UNSP 04017003}, abstract = {

Nanoscale structural heterogeneities were recently revealed in computational and experimental studies of calcium silicate hydrates in hardened cement pastes. In this work their consequences for the mechanics are analyzed by computing local pressures in model samples at different overall densities, corresponding to different initial water-to-cement ratios. The correlations between pore size distributions, local density, local cohesive energy, and local pressure clearly show how in these materials structural heterogeneities may be the origin of significant mechanical heterogeneities. The results indicate that even at high density pressure, heterogeneities develop during the densification of cement hydrates and result in the coexistence of regions of high positive pressure with regions of negative pressure, in spite of the overall mechanical stability of the samples. Furthermore, the regions of negative pressure, prone to mechanical instabilities and local plastic processes, tend to be localized close to the surface of large mesopores and hence to be more significant at higher initial water content.
Inhomogeneity in Cement Hydrates: Linking Local Packing to Local Pressure. Available from: https://www.researchgate.net/publication/314222493_Inhomogeneity_in_Cement_Hydrates_Linking_Local_Packing_to_Local_Pressure [accessed Oct 2, 2017].

}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000120}, author = {Katerina Ioannidou and Emanuela Del Gado and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {613, title = {Inorganic, Hybridized and Living Macrocellular Foams: {\textquotedblleft}Out of the Box{\textquotedblright} Heterogeneous Catalysis}, journal = {The Chemical Record}, volume = {18}, year = {2017}, month = {Dec-01-2017}, pages = {776-787}, abstract = {

With this personal account we show how the Integrative Chemistry, when combining the sol-gel process and concentrated emulsions, allows to trigger inorganic, hybrid or living materials when dedicated toward heterogeneous catalysis applications. In here we focus on 3D-macrocellular monolithic foams bearing hierarchical porosities and applications thereof toward heterogeneous catalysis where both activities and mass transport are enhanced. We thereby first depict the general background of emulsions, focusing on concentrated ones, acting as soft templates for the design of solid (HIPE) foams, HIPE being the acronym for High Internal Phase Emulsions while encompassing both sol-gel and polymer chemistry. Secondly we extend this approach toward the design of inorganic cellular materials labeled Si(HIPE) and hybrid organic-inorganic foams, labeled Organo-Si(HIPE), where heterogeneous catalysis applications are addressed considering acidic, metallic, enzymatic and bacterial-based modified Si-HIPE. Along, we will show how the fluid hydrodynamic within the macrocellular foams is offering advanced \"out of the box\" heterogeneous catalytic capabilities.

}, keywords = {catalysts; emulsion-templating; integrative chemistry; porous matter; sol-gel process}, doi = {10.1002/tcr.v18.7-810.1002/tcr.201700075}, url = {http://doi.wiley.com/10.1002/tcr.v18.7-8}, author = {Roucher, Armand and Depardieu, Martin and Pekin, Deniz and Morvan, Micka{\"e}l and Backov, Renal} } @article {410, title = {Magnetism as indirect tool for carbon content assessment in nickel nanoparticles}, journal = {Journal of Applied Physics}, volume = {122}, year = {2017}, month = {Jul-12-2017}, pages = {213902}, abstract = {

We report a combined experimental and theoretical study to ascertain carbon solubility in nickel nanoparticles embedded into a carbon matrix via the one-pot method. This original approach is based on the experimental characterization of the magnetic properties of Ni at room temperature and Monte Carlo simulations used to calculate the magnetization as a function of C content in Ni nanoparticles. Other commonly used experimental methods fail to accurately determine the chemical analysis of these types of nanoparticles. Thus, we could assess the C content within Ni nanoparticles and it decreases from 8 to around 4\ at. \% with increasing temperature during the synthesis. This behavior could be related to the catalytic transformation of dissolved C in the Ni particles into graphite layers surrounding the particles at high temperature. The proposed approach is original and easy to implement experimentally since only magnetization measurements at room temperature are needed. Moreover, it can be extended to other types of magnetic nanoparticles dissolving carbon.

}, issn = {0021-8979}, doi = {10.1063/1.5006138}, url = {http://aip.scitation.org/doi/10.1063/1.5006138http://aip.scitation.org/doi/pdf/10.1063/1.5006138}, author = {Oumellal, Y. and Yann Magnin and {\'\i}nez de Yuso, A. and Aguiar Hualde, J. M. and Amara, H. and Paul-Boncour, V. and Matei Ghimbeu, C. and Malouche, A. and Christophe Bichara and Roland Jean-Marc Pellenq and Zlotea, C.} } @article {586, title = {Mesoscale Poroelasticity of Heterogeneous Media}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Oct-19-2017}, pages = {Article Number: UNSP 04017016 }, abstract = {

The poromechanics of heterogeneous media is reformulated in a discrete framework using the lattice element method (LEM) that accounts for the presence of interfaces as well as local microtextural and elastic variations. The exchange of mechanical information between pore and solid(s) is captured by means of force field potentials for these domains, which eliminate the requirement of scale separability of continuum-based poromechanics approaches. In congruence with mu VT and NPT ensembles of statistical mechanics, discrete expressions for Biot poroelastic coefficients are derived. Considering harmonic-type interaction potentials for each link, analytical expressions for both isotropic and transversely isotropic effective elasticity are presented. The theory is validated against continuum-based expressions of Biot poroelastic coefficients for porous media with isotropic and transversely isotropic elastic solid behavior. (C) 2017 American Society of Civil Engineers.

}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000136}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000136}, author = {Monfared, Siavash and Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {587, title = {Methodology for Estimation of Nanoscale Hardness via Atomistic Simulations}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Dec-2017}, pages = { Article Number: UNSP 04017011 }, abstract = {

Statistical mechanics has provided powerful techniques to measure mechanical properties of materials at the nanoscale and paved the way for bottom-up computational materials design. The introduction of such techniques in civil engineering applications, namely construction and geotechnical materials, remains limited to the elastic and fracture properties. This paper presents an atomistic approach to calculate the nanoscale cohesion, friction angle, and hardness. This method is based on the application of biaxial external deformation, or stress, in the weakest crystallographic direction in the material. The onset of the failure is characterized by investigating the unloading paths from several points on the stress-strain curve. Such calculations of the failure stress along different deformation paths provide multiple failure Mohr circles in the normal-shear stress space, which is found to provide a failure envelope akin to the Mohr-Coulomb failure criterion that is widely used for the plastic analysis of granular geomaterials. The failure envelope characterizes the nanoscale cohesion and friction angle, which in conjunction with continuum mechanics can be utilized to estimate the nanoscale hardness of layered materials. Application of this method to tobermorite and Na-montmorillonite crystals yields values that are close to the experimental measurements obtained using nanoindentation and atomic force microscopy techniques. (C) 2017 American Society of Civil Engineers.

}, keywords = {Atomistic simulation; Friction; Cohesion; Hardness; Tobermorite; Montmorillonite}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000127}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000127}, author = {Qomi, M. J. Abdolhosseini and Ebrahimi, Davoud and Mathieu Bauchy and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {241, title = {Model representations of kerogen structures: An insight from density functional theory calculations and spectroscopic measurements}, journal = {Scientific Reports}, volume = {7}, year = {2017}, month = {Aug-01-2017}, pages = {Article Number: 7068}, abstract = {

Molecular structures of kerogen control hydrocarbon production in unconventional reservoirs. Significant progress has been made in developing model representations of various kerogen structures. These models have been widely used for the prediction of gas adsorption and migration in shale matrix. However, using density functional perturbation theory (DFPT) calculations and vibrational spectroscopic measurements, we here show that a large gap may still remain between the existing model representations and actual kerogen structures, therefore calling for new model development. Using DFPT, we calculated Fourier transform infrared (FTIR) spectra for six most widely used kerogen structure models. The computed spectra were then systematically compared to the FTIR absorption spectra collected for kerogen samples isolated from Mancos, Woodford and Marcellus formations representing a wide range of kerogen origin and maturation conditions. Limited agreement between the model predictions and the measurements highlights that the existing kerogen models may still miss some key features in structural representation. A combination of DFPT calculations with spectroscopic measurements may provide a useful diagnostic tool for assessing the adequacy of a proposed structural model as well as for future model development. This approach may eventually help develop comprehensive infrared (IR)-fingerprints for tracing kerogen evolution.

}, doi = {10.1038/s41598-017-07310-9}, author = {Weck, Philippe F. and Kim, Eunja and Wang, Yifeng and Kruichak, Jessica N. and Mills, Melissa M. and Matteo, Edward N. and Roland Jean-Marc Pellenq} } @article {246, title = {Molecular Modeling and Adsorption Properties of Ordered Silica-Templated CMK Mesoporous Carbons}, journal = {Langmuir}, volume = {33}, year = {2017}, month = {Feb-06-2017}, pages = {2109 - 2121}, abstract = {

Realistic molecular models of silica-templated CMK-1, CMK-3, and CMK-5 carbon materials have been developed by using carbon rods and carbon pipes that were obtained by adsorbing carbon in a model MCM-41 pore. The interactions between the carbon atoms with the silica matrix were described using the PN-Traz potential, and the interaction between the carbon atoms was calculated by the reactive empirical bond order (REBO) potential. Carbon rods and pipes with different thicknesses were obtained by changing the silica\–carbon interaction strength, the temperature, and the chemical potential of carbon vapor adsorption. These equilibrium structures were further used to obtain the atomic models of CMK-1, CMK-3, and CMK-5 materials using the same symmetry as found in TEM pictures. These models are further refined and made more realistic by adding interconnections between the carbon rods and carbon pipes. We calculated the geometric pore size distribution of the different models of CMK-5 and found that the presence of interconnections results in some new features in the pore size distribution. Argon adsorption properties were investigated using GCMC simulations to characterize these materials at 77 K. We found that the presence of interconnection results greatly improves the agreement with available experimental data by shifting the capillary condensation to lower pressures. Adding interconnections also induces smoother adsorption/condensation isotherms, and desorption/evaporation curves show a sharp jump. These features reflex the complexity of the nanovoids in CMKs in terms of their pore morphology and topology.

Abstract Image

}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.6b04169}, author = {Jain, Surendra Kumar and Roland Jean-Marc Pellenq and Keith E. Gubbins and Peng, Xuan} } @article {623, title = {Nano-granular texture of cement hydrates}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 15027}, abstract = {

Mechanical behavior of concrete crucially depends on cement hydrates, the \“glue\” of cement. The design of high performance and more environmentally friendly cements demands a deeper understanding of the formation of the multiscale structure of cement hydrates, when they precipitate and densify. We investigate the precipitation and setting of nano-grains of cement hydrates using a combination of Monte Carlo and Molecular Dynamics numerical simulations and study their texture from nano up to the micron scale. We characterize the texture of cement hydrates using the local volume fraction distribution, the pore size distribution, the scattering intensity and the chord length distribution and we compare them with experiments. Our nano-granular model provides cement structure with realistic texture and mechanics and can be further used to investigate degradation mechanisms.

}, doi = {10.1051/epjconf/201714015027}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714015027}, author = {Katerina Ioannidou and Franz-Josef Ulm and Pierre E. Levitz and Emanuela Del Gado and Roland Jean-Marc Pellenq}, editor = {Farhang Radja{\"\i} and Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {604, title = {Numerical insight into the micromechanics of jet erosion of a cohesive granular material}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 15017}, abstract = {

Here we investigate the physical mechanisms behind the surface erosion of a cohesive granular soil induced by an impinging jet by means of numerical simulations coupling fluid and grains at the microscale. The 2D numerical model combines the Discrete Element and Lattice Boltzmann methods (DEM-LBM) and accounts for the granular cohesion with a contact model featuring a paraboloidal yield surface. Here we review first the hydrodynamical conditions imposed by the fluid jet on a solid granular packing, turning then the attention to the impact of cohesion on the erosion kinetics. Finally, the use of an additional subcritical debonding damage model based on the work of Silvani and co-workers provides a novel insight into the internal solicitation of the cohesive granular sample by the impinging jet.

}, doi = {10.1051/epjconf/201714015017}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714015017}, author = {Cuellar, Pablo and Benseghier, Zeyd and Luu, Li-Hua and St{\'e}phane Bonelli and Jean-Yves Delenne and Farhang Radja{\"\i} and Pierre Philippe}, editor = {Saeid Nezamabadi and Luding, S.} } @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} } @article {242, title = {A potential-of-mean-force approach for fracture mechanics of heterogeneous materials using the lattice element method}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {105}, year = {2017}, month = {Aug-2017}, pages = {116 - 130}, abstract = {

Fracture of heterogeneous materials has emerged as a critical issue in many engineering applications, ranging from subsurface energy to biomedical applications, and requires a rational framework that allows linking local fracture processes with global fracture descriptors such as the energy release rate, fracture energy and fracture toughness. This is achieved here by means of a local and a global potential-of-mean-force (PMF) inspired Lattice Element Method (LEM) approach. In the local approach, fracture-strength criteria derived from the effective interaction potentials between mass points are shown to exhibit a scaling commensurable with the energy dissipation of fracture processes. In the global PMF-approach, fracture is considered as a sequence of equilibrium states associated with minimum potential energy states analogous to Griffith\’s approach. It is found that this global approach has much in common with a Grand Canonical Monte Carlo (GCMC) approach, in which mass points are randomly removed following a maximum dissipation criterion until the energy release rate reaches the fracture energy. The duality of the two approaches is illustrated through the application of the PMF-inspired LEM for fracture propagation in a homogeneous linear elastic solid using different means of evaluating the energy release rate. Finally, by application of the method to a textbook example of fracture propagation in a heterogeneous material, it is shown that the proposed PMF-inspired LEM approach captures some well-known toughening mechanisms related to fracture energy contrast, elasticity contrast and crack deflection in the considered two-phase layered composite material.

}, issn = {00225096}, doi = {10.1016/j.jmps.2017.05.006}, author = {Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {590, title = {Production of H 2 by water radiolysis in cement paste under electron irradiation: A joint experimental and theoretical study}, journal = {Cement and Concrete Research}, volume = {100}, year = {2017}, month = {Oct-2017}, pages = {110 - 118}, abstract = {

Long-term confinement of nuclear waste is one of the main challenges faced by the nuclear industry. Fission products such as Sr-90 and Cs-137, both beta(-) emitters known to induce serious health hazards, represent the largest fraction of nuclear waste. Cement is a good candidate to store them, provided it can resist the effects of irradiation over time. Here, we have investigated the effects of beta(-) decay on cement by performing electron irradiation experiments on different samples. We show that H-2 production in cement, the main effect of water radiolysis, depends strongly on composition and relative humidity. First-principles calculations indicate that the water-rich interlayer regions with Ca2+ ions act as electron traps that promote the formation of H-2. They also show that holes localize in water-rich regions in low Ca content samples and are then able to participate in H-2 production. This work provides new understanding of radiolysis effects in cements.

}, issn = {00088846}, doi = {10.1016/j.cemconres.2017.05.022}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008884617302065}, author = {Le Caer, Sophie and Dezerald, Lucile and Boukari, Khaoula and Laine, Maxime and Taupin, s{\'e}bastien and Kavanagh, Ryan M. and Johnston, Conrad S.N. and Foy, Eddy and Charpentier, Thibault and Konrad J. Krakowiak and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Tribello, Gareth A. and Kohanoff, Jorge J. and Andres Sa{\`u}l} } @article {589, title = {Role of Interfaces in Elasticity and Failure of Clay{\textendash}Organic Nanocomposites: Toughening upon Interface Weakening?}, journal = {Langmuir}, volume = {33}, year = {2017}, month = {Oct-24-2017}, pages = {11457 - 11466}, abstract = {

Synthetic organic-inorganic composites constitute a new class of engineering materials finding applications in an increasing range of fields. The interface between the constituting phases plays a pivotal role in the enhancement of mechanical properties. In exfoliated clay-organic nanocomposites, individual, high aspect ratio clay sheets are dispersed in the organic matrix providing large interfaces and hence efficient stress transfer. In this study, we aim at elucidating molecular-scale reinforcing mechanisms in a series of model clay-organic composite systems by means of reactive molecular simulations. In our models, two possible locations of failure initiation are present: one is the interlayer space of the clay platelet, and the other one is the clay-organic interface. We systematically modify the cohesiveness of the interface and assess how the failure mechanism changes when the different model composites are subjected to a tensile test. Besides a change in the failure mechanism, an increase in the released energy at the interface (meaning an increased overall toughness) are observed upon weakening the interface by bond removal. We propose a theoretical analysis of these results by considering a cohesive law that captures the effect of the interface on the composite mechanics. We suggest an atomistic interpretation of this cohesive law, in particular, how it relates to the degree of bonding at the interface. In a broader perspective, this work sheds light on the importance of the orthogonal behavior of interfaces to nanocomposites.

https://pubs-acs-org.libproxy.mit.edu/appl/literatum/publisher/achs/journals/content/langd5/2017/langd5.2017.33.issue-42/acs.langmuir.7b01071/20171025/images/medium/la-2017-01071k_0008.gif

}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.7b01071}, url = {http://pubs.acs.org/doi/10.1021/acs.langmuir.7b01071}, author = {Gy{\"o}rgy Hantal and Brochard, Laurent and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Benoit A. Coasne} } @article {607, title = {Scaling behavior of immersed granular flows}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 09044}, abstract = {

The shear behavior of granular materials immersed in a viscous fluid depends on fluid properties (viscosity, density), particle properties (size, density) and boundary conditions (shear rate, confining pressure). Using computational fluid dynamics simulations coupled with molecular dynamics for granular flow, and exploring a broad range of the values of parameters, we show that the parameter space can be reduced to a single parameter that controls the packing fraction and effective friction coefficient. This control parameter is a modified inertial number that incorporates viscous effects.

}, doi = {10.1051/epjconf/201714009044}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714009044}, author = {L. Amarsid and Jean-Yves Delenne and Patrick Mutabaruka and Yann Monerie and Perales, F. and Farhang Radja{\"\i}}, editor = {Saeid Nezamabadi and Luding, S.} } @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 {240, title = {Stress Transmission and Failure in Disordered Porous Media}, journal = {PHYSICAL REVIEW LETTERS}, volume = {119}, year = {2017}, month = {Aug-14-2017}, pages = {Article Number: 075501}, abstract = {

By means of extensive lattice-element simulations, we investigate stress transmission and its relation with failure properties in increasingly disordered porous systems. We observe a non-Gaussian broadening of stress probability density functions under tensile loading with increasing porosity and disorder, revealing a gradual transition from a state governed by single-pore stress concentration to a state controlled by multipore interactions and metric disorder. This effect is captured by the excess kurtosis of stress distributions and shown to be nicely correlated with the second moment of local porosity fluctuations, which appears thus as a (dis)order parameter for the system. By generating statistical ensembles of porous textures with varying porosity and disorder, we derive a general expression for the fracture stress as a decreasing function of porosity and disorder. Focusing on critical sites where the local stress is above the global fracture threshold, we also analyze the transition to failure in terms of a coarse-graining length. These findings provide a general framework which can also be more generally applied to multiphase and structural heterogeneous materials.

Figure

Figure

Figure

Figure

}, doi = {10.1103/PhysRevLett.119.075501}, author = {Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and et al} } @article {593, title = {Stress Transmission and Failure in Disordered Porous Media}, journal = {Physical Review Letters}, volume = {119}, year = {2017}, month = {Aug-14-2017}, pages = { Article Number: 075501 }, abstract = {

By means of extensive lattice-element simulations, we investigate stress transmission and its relation with failure properties in increasingly disordered porous systems. We observe a non-Gaussian broadening of stress probability density functions under tensile loading with increasing porosity and disorder, revealing a gradual transition from a state governed by single-pore stress concentration to a state controlled by multipore interactions and metric disorder. This effect is captured by the excess kurtosis of stress distributions and shown to be nicely correlated with the second moment of local porosity fluctuations, which appears thus as a (dis)order parameter for the system. By generating statistical ensembles of porous textures with varying porosity and disorder, we derive a general expression for the fracture stress as a decreasing function of porosity and disorder. Focusing on critical sites where the local stress is above the global fracture threshold, we also analyze the transition to failure in terms of a coarse-graining length. These findings provide a general framework which can also be more generally applied to multiphase and structural heterogeneous materials.

}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.119.075501}, url = {https://link.aps.org/doi/10.1103/PhysRevLett.119.075501}, author = {Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {243, title = {Topological Control on the Structural Relaxation of Atomic Networks under Stress}, journal = {Physical Review Letters}, volume = {119}, year = {2017}, month = {Jul-21-2017}, pages = {Article Number: 035502}, abstract = {

Upon loading, atomic networks can feature delayed irreversible relaxation. However, the effect of composition and structure on relaxation remains poorly understood. Herein, relying on accelerated molecular dynamics simulations and topological constraint theory, we investigate the relationship between atomic topology and stress-induced structural relaxation, by taking the example of creep deformations in calcium silicate hydrates (C-S-H), the binding phase of concrete. Under constant shear stress, C-S-H is found to feature delayed logarithmic shear deformations. We demonstrate that the propensity for relaxation is minimum for isostatic atomic networks, which are characterized by the simultaneous absence of floppy internal modes of relaxation and eigenstress. This suggests that topological nanoengineering could lead to the discovery of nonaging materials.

Figure

Figure

Figure

}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.119.035502}, author = {Mathieu Bauchy and Wang, Mengyi and Yu, Yingtian and Wang, Bu and Krishnan, M. Anoop and Enrico Masoero and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {263, title = {Viscoinertial regime of immersed granular flows}, journal = {Physical Review E}, volume = {96}, year = {2017}, month = {Oct-10-2017}, pages = {Article Number: 012901}, abstract = {

By means of extensive coupled molecular dynamics\–lattice Boltzmann simulations, accounting for grain dynamics and subparticle resolution of the fluid phase, we analyze steady inertial granular flows sheared by a viscous fluid. We show that, for a broad range of system parameters (shear rate, confining stress, fluid viscosity, and relative fluid-grain density), the frictional strength and packing fraction can be described by a modified inertial number incorporating the fluid effect. In a dual viscous description, the effective viscosity diverges as the inverse square of the difference between the packing fraction and its jamming value, as observed in experiments. We also find that the fabric and force anisotropies extracted from the contact network are well described by the modified inertial number, thus providing clear evidence for the role of these key structural parameters in dense suspensions.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, issn = {2470-0045}, doi = {10.1103/PhysRevE.96.012901}, author = {L. Amarsid and Jean-Yves Delenne and Patrick Mutabaruka and Yann Monerie and Perales, F. and Farhang Radja{\"\i}} } @article {268, title = {Binary mixtures of disks and elongated particles: Texture and mechanical properties}, journal = {Physical Review E}, volume = {94}, year = {2016}, month = {Oct-14-2016}, pages = {Article Number: 042901}, abstract = {

We analyze the shear strength and microstructure of binary granular mixtures consisting of disks and elongated particles by varying systematically both the mixture ratio and degree of homogeneity (from homogeneous to fully segregated). The contact dynamics method is used for numerical simulations with rigid particles interacting by frictional contacts. A counterintuitive finding of this work is that the shear strength, packing fraction, and, at the microscopic scale, the fabric, force, and friction anisotropies of the contact network are all nearly independent of the degree of homogeneity. In other words, homogeneous mixtures have the same strength properties as segregated packings of the two particle shapes. In contrast, the shear strength increases with the proportion of elongated particles correlatively with the increase of the corresponding force and fabric anisotropies. By a detailed analysis of the contact network topology, we show that various contact types contribute differently to force transmission and friction mobilization.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, issn = {2470-0045}, doi = {10.1103/PhysRevE.94.042901}, author = {Emilien Az{\'e}ma and Preechawuttipong, Itthichai and Farhang Radja{\"\i}} } @article {643, title = {Cement As a Waste Form for Nuclear Fission Products: The Case of 90Sr and Its Daughters}, journal = {Environmental Science \& Technology}, volume = {49}, year = {2016}, month = {Nov-17-2016}, pages = {13676-13683}, abstract = {

One of the main challenges faced by the nuclear industry is the long-term confinement of nuclear waste. Because it is inexpensive and easy to manufacture, cement is the material of choice to store large volumes of radioactive materials, in particular the low-level medium-lived fission products. It is therefore of utmost importance to assess the chemical and structural stability of cement containing radioactive species. Here, we use ab initio calculations based on density functional theory (DFT) to study the effects of Sr-90 insertion and decay in C-S-H (calcium-silicate-hydrate) in order to test the ability of cement to trap and hold this radioactive fission product and to investigate the consequences of its beta-decay on the cement paste structure. We show that Sr-90 is stable when it substitutes the Ca2+ ions in C-S-H, and so is its daughter nucleus Y-90 after beta-decay. Interestingly, Zr-90, daughter of Y-90 and final product in the decay sequence, is found to be unstable compared to the bulk phase of the element at zero K but stable when compared to the solvated ion in water. Therefore, cement appears as a suitable waste form for Sr-90 storage.

}, issn = {0013-936X}, doi = {10.1021/acs.est.5b02609}, url = {http://pubs.acs.org/doi/10.1021/acs.est.5b02609}, author = {Dezerald, Lucile and Kohanoff, Jorge J. and Alfredo A. Correa and Caro, Alfredo and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Andres Sa{\`u}l} } @article {252, title = {Data analytics for simplifying thermal efficiency planning in cities}, journal = {Journal of The Royal Society Interface}, volume = {13}, year = {2016}, month = {Apr-2016}, pages = {Article Number: 20150971}, abstract = {

More than 44\% of building energy consumption in the USA is used for space heating and cooling, and this accounts for 20\% of national CO2emissions. This prompts the need to identify among the 130 million households in the USA those with the greatest energy-saving potential and the associated costs of the path to reach that goal. Whereas current solutions address this problem by analysing each building in detail, we herein reduce the dimensionality of the problem by simplifying the calculations of energy losses in buildings. We present a novel inference method that can be used via a ranking algorithm that allows us to estimate the potential energy saving for heating purposes. To that end, we only need consumption from records of gas bills integrated with a building\&$\#$39;s footprint. The method entails a statistical screening of the intricate interplay between weather, infrastructural and residents\&$\#$39; choice variables to determine building gas consumption and potential savings at a city scale. We derive a general statistical pattern of consumption in an urban settlement, reducing it to a set of the most influential buildings\&$\#$39; parameters that operate locally. By way of example, the implications are explored using records of a set of (N= 6200) buildings in Cambridge, MA, USA, which indicate that retrofitting only 16\% of buildings entails a 40\% reduction in gas consumption of the whole building stock. We find that the inferred heat loss rate of buildings exhibits a power-law data distribution akin to Zipf\&$\#$39;s law, which provides a means to map an optimum path for gas savings per retrofit at a city scale. These findings have implications for improving the thermal efficiency of cities\&$\#$39; building stock, as outlined by current policy efforts seeking to reduce home heating and cooling energy consumption and lower associated greenhouse gas emissions.

}, issn = {1742-5689}, doi = {10.1098/rsif.2015.0971}, author = {Mohammad Javad Abdolhosseini Qomi and Noshadravan, Arash and Sobstyl, Jake M. and Toole, Jameson and Ferreira, Joseph and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Gonzalez, Marta C.} } @article {321, title = {Early-Age Stress and Pressure Developments in a Wellbore Cement Liner: Application to Eccentric Geometries}, journal = {Journal of Applied Mechanics}, volume = {83}, year = {2016}, month = {Sep-2016}, pages = {Article Number: 091012}, abstract = {

This paper introduces a predictive model for the stress and pressure evolutions in a wellbore cement liner at early ages. A pressure state equation is derived that observes the coupling of the elastic changes of the solid matrix, the eigenstress developments in the solid and porespaces, and the mass consumption of water in course of the reaction. Here, the transient constitution of the solid volume necessitates advancing the mechanical state of the poroelastic cement skeleton incrementally and at constant hydration degree. Next, analytic function theory is employed to assess the localization of stresses along the steel\–cement (SC) and rock\–cement (RC) interfaces by placing the casing eccentrically with respect to the wellbore hole. Though the energy release rate due to complete debonding of either interface is only marginally influenced by the eccentricity, the risk of evolving a microcrack along the thick portion of the sheath is substantially increased. Additionally, it is observed that the risk of microannulus formation is principally affected by the pressure rebound, which is engendered by the slowing reaction rate and amplified for rock boundaries with low permeability.

}, issn = {0021-8936}, doi = {10.1115/1.4034013}, author = {Thomas Alexander Petersen and Franz-Josef Ulm} } @article {250, title = {Effect of Polydispersity of Clay Platelets on the Aggregation And Mechanical Properties of Clay at the Mesoscale}, journal = {Clays and Clay Minerals}, volume = {64}, year = {2016}, month = {Aug-2016}, pages = {425 - 437}, abstract = {

The results from mesoscale simulations of the formation and evolution of microstructure for assemblies of Na-smectite particles based on assumed size distributions of individual clay platelets are presented here. The analyses predicted particle arrangements and aggregation (i.e. platelets linked in face-face configurations) and are used to link geometric properties of the microstructure and mechanical properties of the particle assemblies. Interactions between individual ellipsoidal clay platelets are represented using the Gay-Berne potential based on atomistic simulations of the free energy between two Na-smectite clay-platelets in liquid water, following a novel coarse-graining method developed previously. The current study describes the geometric (aggregate thickness, orientation, and porosity) and elastic properties in the \‘jammed states\’ from the mesoscale simulations for selected ranges of clay particle sizes and confining pressures. The thickness of clay aggregates for monodisperse assemblies increases (with average stack thickness consisting of n = 3-8 platelets) with the diameter of the individual clay platelets and with the level of confining pressure. Aggregates break down at high confining pressures (50-300 atm) due to slippage between the platelets. Polydisperse simulations generate smaller aggregates (n = 2) and show much smaller effects of confining pressure. All assemblies show increased order with confining pressure, implying more anisotropic microstructure. The mesoscale simulations are also in good agreement with macroscopic compression behavior measured in conventional 1-D laboratory compression tests. The mesoscale assemblies exhibit cubic symmetry in elastic properties. The results for larger platelets (D = 1000 \Å) are in good agreement with Nano-indentation measurements on natural clays and shale samples.

}, issn = {0009-8604}, doi = {10.1346/CCMN.2016.0640407}, author = {Ebrahimi, Davoud and Andrew J. Whittle and Roland Jean-Marc Pellenq} } @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 {248, title = {Fracture toughness anomalies: Viewpoint of topological constraint theory}, journal = {Acta Materialia}, volume = {121}, year = {2016}, month = {Dec-2016}, pages = {234 - 239}, abstract = {

The relationship between composition, structure, and resistance to fracture remains poorly understood. Here, based on molecular dynamics simulations, we report that sodium silicate glasses and calcium\–silicate\–hydrates feature an anomalous maximum in fracture toughness. In the framework of topological constraint theory, this anomaly is correlated to a flexible-to-rigid transition, driven by pressure or composition for sodium silicate and calcium\–silicate\–hydrates, respectively. This topological transition, observed for an isostatic network, is also shown to correspond to a ductile-to-brittle transition. At this state, the network is rigid but free of eigen-stress and features stress relaxation through crack blunting, resulting in optimal resistance to fracture. Our topological approach could therefore enable the computational design of tough inorganic solids, which has long been a \“holy grail\” within the non-metallic materials chemistry community.

Image 1

}, issn = {13596454}, doi = {10.1016/j.actamat.2016.09.004}, author = {Mathieu Bauchy and Wang, Bu and Wang, Mengyi and Yu, Yingtian and Mohammad Javad Abdolhosseini Qomi and Smedskjaer, Morten M. and Christophe Bichara and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {249, title = {Free Volume Theory of Hydrocarbon Mixture Transport in Nanoporous Materials}, journal = {The Journal of Physical Chemistry Letters}, volume = {7}, year = {2016}, month = {Jun-10-2016}, pages = {3712 - 3717}, abstract = {
Despite recent focus on shale gas, hydrocarbon recovery from the ultraconfining and disordered porosity of organic matter in shales (kerogen) remains poorly understood. Key aspects such as the breakdown of hydrodynamics at the nanoscale and strong adsorption effects lead to unexplained non-Darcy behaviors. Here, molecular dynamics and statistical mechanics are used to elucidate hydrocarbon mixture transport through a realistic molecular model of kerogen [Bousige, C.; et al. Nat. Mater. 2016, 15, 576]. Owing to strong adsorption effects, velocity cross-correlations between the mixture components and between molecules of the same species are shown to be negligible. This allows estimation of each component permeance from its self-diffusivity, which can be obtained from single-component data. These permeances are found to scale with the reciprocal of the alkane length and decrease with the number of adsorbed molecules following a simple free volume theory, therefore allowing mixture transport prediction as a function of the amount of trapped fluid.
Abstract Image
}, issn = {1948-7185}, doi = {10.1021/acs.jpclett.6b01684}, author = {Ama{\"e}l Obliger and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Benoit A. Coasne} } @article {196, title = {Mesoscale texture of cement hydrates}, journal = {Proceedings of the National Academy of Sciences}, volume = {113}, year = {2016}, month = {Feb-23-2016}, pages = {2029 - 2034}, abstract = {

Strength and other mechanical properties of cement and concrete rely upon the formation of calcium\–silicate\–hydrates (C\–S\–H) during cement hydration. Controlling structure and properties of the C\–S\–H phase is a challenge, due to the complexity of this hydration product and of the mechanisms that drive its precipitation from the ionic solution upon dissolution of cement grains in water. Departing from traditional models mostly focused on length scales above the micrometer, recent research addressed the molecular structure of C\–S\–H. However, small-angle neutron scattering, electron-microscopy imaging, and nanoindentation experiments suggest that its mesoscale organization, extending over hundreds of nanometers, may be more important. Here we unveil the C\–S\–H mesoscale texture, a crucial step to connect the fundamental scales to the macroscale of engineering properties. We use simulations that combine information of the nanoscale building units of C\–S\–H and their effective interactions, obtained from atomistic simulations and experiments, into a statistical physics framework for aggregating nanoparticles. We compute small-angle scattering intensities, pore size distributions, specific surface area, local densities, indentation modulus, and hardness of the material, providing quantitative understanding of different experimental investigations. Our results provide insight into how the heterogeneities developed during the early stages of hydration persist in the structure of C\–S\–H and impact the mechanical performance of the hardened cement paste. Unraveling such links in cement hydrates can be groundbreaking and controlling them can be the key to smarter mix designs of cementitious materials.

Fig. 1.

}, issn = {0027-8424}, doi = {10.1073/pnas.1520487113}, author = {Katerina Ioannidou and Konrad J. Krakowiak and Mathieu Bauchy and Christian G. Hoover and Enrico Masoero and Sidney Yip and Franz-Josef Ulm and Pierre E. Levitz and Roland Jean-Marc Pellenq and Emanuela Del Gado} } @article {640, title = {Nanoporous chalcogenides for adsorption and gas separation}, journal = {Physical Chemistry Chemical Physics}, volume = {18}, year = {2016}, month = {May-21-2016}, pages = {13449-13458}, abstract = {

The adsorption and gas separation properties of amorphous porous chalcogenides such as GeS2 are investigated using statistical mechanics molecular simulation. Using a realistic molecular model of such amorphous adsorbents, we show that they can be used efficiently to separate different gases relevant to environmental and energy applications (H-2, CO2, CH4, N-2). In addition to shedding light on the microscopic adsorption mechanisms, we show that coadsorption in this novel class of porous materials can be described using the ideal adsorbed solution theory (IAST). Such a simple thermodynamic model, which allows avoiding complex coadsorption measurements, describes the adsorption of mixture from pure component adsorption isotherms. Our results, which are found to be in good agreement with available experimental data, paves the way for the design of gas separation membranes using the large family of porous chalcogenides.

}, issn = {1463-9076}, doi = {10.1039/C6CP00467A}, url = {http://xlink.rsc.org/?DOI=C6CP00467A}, author = {Guido Ori and Carlo Massobrio and Pradel, Annie and Ribes, Miche and Benoit A. Coasne} } @article {324, title = {Radial fracture in a three-phase composite: Application to wellbore cement liners at early ages}, journal = {Engineering Fracture Mechanics}, volume = {154}, year = {2016}, month = {Mar-2016}, pages = {272 - 287}, abstract = {

Little understanding exists between the early-age stress developments in a wellbore cement sheath and its risk of impairment. During hydration, the cement morphology and pore-pressure changes induce eigenstresses in the solid and pore volumes. Utilizing these stresses as the driving mechanism of fracture, this paper formalizes the inspection of a radial crack in an elastic cement sheath constrained by an inner steel casing and an outer rock formation. The solution is constructed in the framework of analytic function theory and seeks the Green\’s function for an edge dislocation in the intermediate cement phase. A dislocation pile-up along the line of fracture constructs a singular integral equation for the crack opening displacement derivative, from which the energy release rate is readily deduced.

Under the uniform development of eigenstresses, the stiffness ratios of steel-to-cement and rock-to-cement generally predict the crack to initiate along the steel-cement interface. Here, the impacts of (i) a rigid bond and (ii) a sliding interface with no shear are assessed. This leads to the primary result of the paper: the potential for radial fracture is substantially mitigated by ensuring the shear connection between the steel casing and the cement sheath.

Diagram of a wellbore cement system with a radial crack emanating from SC.

Simulated evolution of the shear modulus G, bulk modulus K, and Poisson{\textquoteright}s ratio ...

Simulated evolution of the effective hoop stress Σθθ=σθθ+p along SC (blue) and ...

Continuation of the annular region of the cement domain across SC and RC.

}, issn = {00137944}, doi = {10.1016/j.engfracmech.2016.01.005}, author = {Thomas Alexander Petersen and Franz-Josef Ulm} } @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} } @article {179, title = {Advances in design and modeling of porous materials}, journal = {The European Physical Journal Special Topics}, volume = {224}, year = {2015}, month = {Jul-2015}, pages = {1653 - 1653}, abstract = {

This special issue of the European Physical Journal Special Topics is dedicated to selected papers from the symposium \“High surface area porous and granular materials\” organized in the frame of the conference \“Mat\ériaux 2014\”, held on November 24\–28, 2014 in Montpellier, France.

Porous materials and granular materials gather a wide variety of heterogeneous, isotropic or anisotropic media made of inorganic, organic or hybrid solid skeletons, with open or closed porosity, and pore sizes ranging from the centimeter scale to the sub-nanometer scale. Their technological and industrial applications cover numerous areas from building and civil engineering to microelectronics, including also metallurgy, chemistry, health, waste water and gas effluent treatment. Many emerging processes related to environmental protection and sustainable development also rely on this class of materials. Their functional properties are related to specific transfer mechanisms (matter, heat, radiation, electrical charge), to pore surface chemistry (exchange, adsorption, heterogeneous catalysis) and to retention inside confined volumes (storage, separation, exchange, controlled release). The development of innovative synthesis, shaping, characterization and modeling approaches enables the design of advanced materials with enhanced functional performance. The papers collected in this special issue offer a good overview of the state-of-the-art and science of these complex media.

We would like to thank all the speakers and participants for their contribution to the success of the symposium. We also express our gratitude to the organization committee of \“Mat\´eriaux 2014\”. We finally thank the reviewers and the staff of the European Physical Journal Special Topics who made the publication of this special issue possible.

}, issn = {1951-6355}, doi = {10.1140/epjst/e2015-02487-6}, author = {Ayral, Andre and Calas-Etienne, Sylvie and Benoit A. Coasne and Deratani, Andr{\'e} and Evstratov, Alexis and Anne Galarneau and Grande, Daniel and Hureau, Matthieu and Jobic, Herv{\'e} and Morlay, Catherine and Parmentier, Julien and Prelot, B{\'e}n{\'e}dicte and Rossignol, Sylvie and Simon-Masseron, Ang{\'e}lique and Thibault-Starzyk, Fr{\'e}d{\'e}ric} } @article {128, title = {Bottom-up model of adsorption and transport in multiscale porous media.}, journal = {Phys Rev E Stat Nonlin Soft Matter Phys}, volume = {91}, year = {2015}, month = {Mar-20-2015}, pages = {Article Number: 032133 }, abstract = {

We develop a model of transport in multiscale porous media which accounts for adsorption in the different porosity scales. This model employs statistical mechanics to upscale molecular simulation and describe adsorption and transport at larger time and length scales. Using atom-scale simulations, which capture the changes in adsorption and transport with temperature, pressure, pore size, etc., this approach does not assume any adsorption or flow type. Moreover, by relating the local chemical potential \μ(r) and density \ρ(r), the present model accounts for adsorption effects and possible changes in the confined fluid state upon transport. This model constitutes a bottom-up framework of adsorption and transport in multiscale materials as it (1) describes the adsorption-transport interplay, (2) accounts for the hydrodynamics breakdown at the nm scale, and (3) is multiscale.

}, issn = {1550-2376}, doi = {10.1103/PhysRevE.91.032133}, author = {Botan, Alexandru and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @article {136, title = {Capturing material toughness by molecular simulation: accounting for large yielding effects and limits}, journal = {International Journal of Fracture}, volume = {194}, year = {2015}, month = {Aug-2015}, pages = {149 - 167}, abstract = {

The inherent computational cost of molecular simulations limits their use to the study of nanometric systems with potentially strong size effects. In the case of fracture mechanics, size effects due to yielding at the crack tip can affect strongly the mechanical response of small systems. In this paper we consider two examples: a silica crystal for which yielding is limited to a few atoms at the crack tip, and a nanoporous polymer for which the process zone is about one order of magnitude larger. We perform molecular simulations of fracture of those materials and investigate in particular the system and crack size effects. The simulated systems are periodic with an initial crack. Quasi-static loading is achieved by increasing the system size in the direction orthogonal to the crack while maintaining a constant temperature. As expected, the behaviors of the two materials are significantly different. We show that the behavior of the silica crystal is reasonably well described by the classical framework of linear elastic fracture mechanics (LEFM). Therefore, one can easily upscale engineering fracture properties from molecular simulation results. In contrast, LEFM fails capturing the behavior of the polymer and we propose an alternative analysis based on cohesive crack zone models. We show that with a linear decreasing cohesive law, this alternative approach captures well the behavior of the polymer. Using this cohesive law, one can anticipate the mechanical behavior at larger scale and assess engineering fracture properties. Thus, despite the large yielding of the polymer at the scale of the molecular simulation, the cohesive zone analysis offers a proper upscaling methodology.

}, issn = {0376-9429}, doi = {10.1007/s10704-015-0045-y}, author = {Brochard, Laurent and Gy{\"o}rgy Hantal and Hadrien Laubie and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @proceedings {98, title = {Is cement a glassy material?}, journal = {Euro-C Conference}, volume = {COMPUTATIONAL MODELLING OF CONCRETE STRUCTURES, VOL 1}, year = {2015}, month = {Jun-22-2015}, pages = {169-176}, publisher = {CRC PRESS-TAYLOR \& FRANCIS GROUP}, address = {MAR 24-27, 2014, St Anton am Alberg, AUSTRIA}, abstract = {

The nature of Calcium-Silicate-Hydrate (C-S-H), the binding phase of cement, remains a controversial question. In particular, contrary to the former crystalline model, it was recently proposed that its nanoscale structure was actually amorphous. To elucidate this issue, we analyzed the structure of a realistic simulation of C-S-H, and compared the latter to crystalline tobermorite, a natural analogue to cement, and to an artificial ideal glass. Results clearly support that C-S-H is amorphous. However, its structure shows an intermediate degree of order, retaining some characteristics of the crystal while acquiring an overall glass-like disorder. Thanks to a detailed quantification of order and disorder, we show that its amorphous state mainly arises from its hydration.

}, isbn = {978-1-138-02641-4; 978-1-315-76203-6}, doi = {10.1201/b16645-19}, author = {Mathieu Bauchy and Mohammad Javad Abdolhosseini Qomi and Roland Jean-Marc Pellenq and Franz-Josef Ulm}, editor = {Bicanic, N and Mang, H and Meschke, G and DeBorst, R} } @article {133, title = {Cement As a Waste Form for Nuclear Fission Products: The Case of 90 Sr and Its Daughters}, journal = {Environ Sci Technol}, volume = {49}, year = {2015}, month = {Oct-2015}, pages = {13676-83}, abstract = {

One of the main challenges faced by the nuclear industry is the long-term confinement of nuclear waste. Because it is inexpensive and easy to manufacture, cement is the material of choice to store large volumes of radioactive materials, in particular the low-level medium-lived fission products. It is therefore of utmost importance to assess the chemical and structural stability of cement containing radioactive species. Here, we use ab initio calculations based on density functional theory (DFT) to study the effects of (90)Sr insertion and decay in C-S-H (calcium-silicate-hydrate) in order to test the ability of cement to trap and hold this radioactive fission product and to investigate the consequences of its \β-decay on the cement paste structure. We show that (90)Sr is stable when it substitutes the Ca(2+) ions in C-S-H, and so is its daughter nucleus (90)Y after \β-decay. Interestingly, (90)Zr, daughter of (90)Y and final product in the decay sequence, is found to be unstable compared to the bulk phase of the element at zero K but stable when compared to the solvated ion in water. Therefore, cement appears as a suitable waste form for (90)Sr storage.

}, issn = {1520-5851}, doi = {10.1021/acs.est.5b02609}, author = {Dezerald, Lucile and Kohanoff, Jorge J and Alfredo A. Correa and Caro, Alfredo and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Andres Sa{\`u}l} } @proceedings {332, title = {Creep of Bulk C-S-H: Insights from Molecular Dynamics Simulations}, 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 = {511-516}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA}, abstract = {

Understanding the physical origin of creep in calcium\–silicate\–hydrate (C\–S\–H) is of primary importance, both for fundamental and practical interest. Here, we present a new method, based on molecular dynamics simulation, allowing us to simulate the long-term visco-elastic deformations of C\–S\–H. Under a given shear stress, C\–S\–H features a gradually increasing shear strain, which follows a logarithmic law. The computed creep modulus is found to be independent of the shear stress applied and is in excellent agreement with nanoindentation measurements, as extrapolated to zero porosity.

}, doi = {10.1061/978078447934610.1061/9780784479346.061}, author = {Mathieu Bauchy and Enrico Masoero and Franz-Josef Ulm and Roland Jean-Marc Pellenq}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @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} } @proceedings {331, title = {C-S-H across Length Scales: From Nano to Micron}, 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 = {39-48}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA}, abstract = {

Despite their impact on longevity, serviceability, and environmental footprint of our built infrastructure, the chemo-physical origins of nanoscale properties of cementitious materials, and their link to macroscale properties still remain rather obscure. Here, we discuss a multi-scale approach that describes different aspects of physical properties of C-S-H at the nano- and meso-scales. These include dynamics of water, thermal properties and mechanical behavior of C-S-H and its effect on properties of cement paste at different scales.

}, doi = {10.1061/978078447934610.1061/9780784479346.006}, url = {http://ascelibrary.org/doi/book/10.1061/9780784479346http://ascelibrary.org/doi/pdf/10.1061/9780784479346http://ascelibrary.org/doi/10.1061/9780784479346.006http://ascelibrary.org/doi/pdf/10.1061/9780784479346.006}, author = {Mohammad Javad Abdolhosseini Qomi and Enrico Masoero and Mathieu Bauchy and Franz-Josef Ulm and Emanuela Del Gado and Roland Jean-Marc Pellenq}, 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 {135, title = {Effect of Chain Length and Pore Accessibility on Alkane Adsorption in Kerogen}, journal = {Energy \& Fuels}, volume = {29}, year = {2015}, month = {Oct-28-2015}, pages = {7889 - 7896}, abstract = {
Abstract Image

Configurational-biased Grand Canonical Monte Carlo and Molecular Dynamics simulations were used to investigate the adsorption and structure of n-alkanes in realistic models of kerogen (the organic phase in gas shales). Both the effects of the n-alkane length (from methane to dodecane) and of the adsorbent porosity/density and connectivity were considered. For all n-alkanes, due to the subnano pore size of kerogen, adsorption follows a Langmuir-type adsorption isotherm as the adsorbed amount increases continuously with pressure until maximum loading is reached. While all n-alkanes get adsorbed very efficiently in low density kerogen, size exclusion in high density kerogen lead to low adsorbed amounts for molecules larger than hexane. Using Molecular Dynamics simulations reflecting the setup of actual adsorption experiments, we also probed the effect of pore accessibility on n-alkane adsorption in kerogen. All pores are connected to the surface and, hence, accessible for the low density kerogen. In contrast, most pores are isolated and inaccessible in the kerogen with the high density. We also address the effect of mesoporosity and its connection to nanoporosity on alkane adsorption in kerogen.

}, issn = {0887-0624}, doi = {10.1021/acs.energyfuels.5b02015}, author = {Kerstyn I. Falk and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Benoit A. Coasne} } @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 {137, title = {Fracture toughness of calcium{\textendash}silicate{\textendash}hydrate from molecular dynamics simulations}, journal = {Journal of Non-Crystalline Solids}, volume = {419}, year = {2015}, month = {Jul-01-2015}, pages = {58 - 64}, abstract = {

Concrete is the most widely manufactured material in the world. Its binding phase, calcium\–silicate\–hydrate (C\–S\–H), is responsible for its mechanical properties and has an atomic structure fairly similar to that of usual calcium silicate glasses, which makes it appealing to study this material with tools and theories traditionally used for non-crystalline solids. Here, following this idea, we use molecular dynamics simulations to evaluate the fracture toughness of C\–S\–H, inaccessible experimentally. This allows us to discuss the brittleness of the material at the atomic scale. We show that, at this scale, C\–S\–H breaks in a ductile way, which prevents one from using methods based on linear elastic fracture mechanics. Knowledge of the fracture properties of C\–S\–H at the atomic scale opens the way for an upscaling approach to the design of tougher cement paste, which would allow for the design of slender environment-friendly infrastructures, requiring less material.

}, issn = {00223093}, doi = {10.1016/j.jnoncrysol.2015.03.031}, author = {Mathieu Bauchy and Hadrien Laubie and Mohammad Javad Abdolhosseini Qomi and Christian G. Hoover and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @proceedings {335, title = {Hydration Kinetics and Gel Morphology of C-S-H}, 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 = {565-573}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA}, abstract = {

Calcium-silicate hydrate (C-S-H) is the main binder in cement and concrete. It starts forming from the early stages of cement hydration and it progressively densifies as cement sets. C-S-H nanoscale building blocks form a cohesive gel, whose structure and mechanics are still poorly understood, in spite of its practical importance. Here we review a statistical physics approach recently developed, which allows us to investigate the C-S-H gel formation under the out-of-equilibrium conditions typical of cement hydration. Our approach is based on colloidal particles, precipitating in the pore solution and interacting with effective forces associated to the ionic environment. We present the evolution of the space filling of C-S-H with different particle interactions and compare them with experimental data at different lime concentrations. Moreover, we discuss the structural features of C-S-H in the mesoscale in terms of the scattering intensity. The comparison of our early stage C-S-H structures with small angle neutron scattering (SANS) experiments shows that long range spatial correlations and structural heterogeneties that develop in that early stages of hydration persist also in the hardened paste.

}, doi = {10.1061/978078447934610.1061/9780784479346.069}, author = {Katerina Ioannidou and Enrico Masoero and Pierre E. Levitz and Roland Jean-Marc Pellenq and Emanuela Del Gado}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @proceedings {334, title = {Kinetic Simulations of Cement Creep: Mechanisms from Shear Deformations of Glasses}, 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 = {555-564}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA}, abstract = {

The logarithmic deviatoric creep of cement paste is a technical and scientific challenge. Transition State Theory (TST) indicates that some nanoscale mechanisms of shear deformation, associated with a specific kind of strain hardening, can explain the type of deviatoric creep observed experimentally in mature cement pastes. To test this possible explanation, we simulate the shear deformations of a colloidal model of cement hydrates at the nanoscale. Results from quasi-static simulations indicate a strain hardening analogous to that postulated by the TST approach. Additional results from oscillatory shear (fatigue) simulations show an increase of deformation with number of loading cycles that is consistent with the observed creep. These findings indicate that nanoscale simulations can improve our current understanding of the mechanisms underlying creep, with potential to go beyond the logarithmic creep and explore the onset of failure during tertiary creep.

}, doi = {10.1061/978078447934610.1061/9780784479346.068}, author = {Enrico Masoero and Mathieu Bauchy and Emanuela Del Gado and Hegoi Manzano and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Sidney Yip}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @article {183, title = {Mechanism of H2O Insertion and Chemical Bond Formation in AlPO 4-54{\textperiodcentered}x H2O at High Pressure}, journal = {Journal of the American Chemical Society}, volume = {137}, year = {2015}, month = {Jan-21-2015}, pages = {584 - 587}, abstract = {

The insertion of H2O in AlPO4-54\•xH2O at high pressure was investigated by single-crystal x-ray diffraction and Monte Carlo molecular simulation. H2O molecules are concentrated, in particular near the pore walls. Upon insertion, the additional water is highly disordered. Insertion of H2O (superhydration) is found to impede pore collapse in the material, thereby strongly modifying its mechanical behavior. However, instead of stabilizing the structure with respect to amorphization, the results provide evidence for the early stages of chemical bond formation between H2O molecules and tetrahedrally coordinated aluminum, which is at the origin of the amorphization/reaction process.

}, issn = {0002-7863}, doi = {10.1021/ja511153n}, author = {Alabarse, Frederico G. and Rouquette, J{\'e}r{\^o}me and Benoit A. Coasne and Haidoux, Abel and Paulmann, Carsten and Cambon, Olivier and Haines, Julien} } @proceedings {330, title = {The Meso-Scale Texture of Cement Hydrate Gels: Out-of-Equilibrium Evolution and Thermodynamic Driving}, journal = {10th International Conference on Mechanics and Physics of Creep, Shrinkage, and Durability of Concrete and Concrete Structures CONCREEP 10}, volume = {CONCREEP 10: MECHANICS AND PHYSICS OF CREEP, SHRINKAGE, AND DURABILITY OF CONCRETE AND CONCRETE STRUCTURES}, year = {2015}, month = {Sep-17-2015}, pages = {34-38}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015 Vienna, AustriaReston, VA}, abstract = {

By the end of cement hydration calcium-silicate-hydrate (C-S-H) gels extends over tens and hundreds of nanometers. Their complex texture affects directly, and to a large extent, the macroscopic hygrothermal and mechanical behavior of cement. Here we review a statistical physics approach recently developed, which allows us to investigate the gel formation under the out-of-equilibrium conditions typical of cement hydration and the role of the nano-scale structure in C\—S\—H mechanics upon hardening. Our investigations have unveiled the role, in the C-S-H gels, of nano-scale structural and mechanical heterogeneities that develop due to the the far-from-equilibrium physico-chemical environment in which the material forms. A subtle interplay between the out-of-equilibrium evolution and the effective interactions emerging between the nano-scale units of the gels at different stages of the hydration process ultimately determines the mesoscale texture of cement hydrates and their material properties.

}, doi = {10.1061/978078447934610.1061/9780784479346.005}, author = {Emanuela Del Gado and Katerina Ioannidou and Enrico Masoero and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Sidney Yip}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @proceedings {358, title = {Micromechanical analysis of the surface erosion of a cohesive soil by means of a coupled LBM-DEM model}, journal = {IV International Conference on Particle-based Methods (PARTICLES 2015)}, volume = {PARTICLE-BASED METHODS IV-FUNDAMENTALS AND APPLICATIONS}, year = {2015}, month = {sept-30-2015}, pages = {519-528}, publisher = {International Center for Numerical Methods in Engineering (CIMNE)}, address = {SEP 28-30 2015 Barcelona, SPAIN}, abstract = {

The elementary mechanisms driving the ubiquitous surface erosion of cohesive geomaterials can be analysed from a micromechanical perspective combining well-known numerical techniques. Here, a coupled model combining the Discrete Element and Lattice Boltzmann methods (DEM-LBM) provides an insight into the solid-fluid interaction during the transient erosion caused by a vertical fluid jet impinging on the surface of a granular assembly. The brittle cementation providing cohesion between the solid grains is described here by means of a simple bond model with a single-parameter yield surface. The initial topology of the surface erosion tends to mimic the profile of fluid velocity directly above the soil surface. We find that both the rate of erosion and the magnitude of eroded mass depend directly on the micromechanical strength of the single solid bonds.

}, keywords = {Cohesion, DEM, LBM, Surface erosion}, url = {https://hal.archives-ouvertes.fr/hal-01269324}, author = {Cuellar, Pablo and Pierre Philippe and St{\'e}phane Bonelli and Benahmed, Nadia and Brunier-Coulin, Florian and Jeff Ngoma and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Onate, E and Bischoff, M and Owen, DRJ} } @article {142, title = {Optimized molecular reconstruction procedure combining hybrid reverse Monte Carlo and molecular dynamics}, journal = {The Journal of Chemical Physics}, volume = {142}, year = {2015}, month = {Mar-21-2015}, pages = {Article Number: 114112}, abstract = {

We report an efficient atom-scale reconstruction method that consists of combining the Hybrid Reverse Monte Carlo algorithm (HRMC) with Molecular Dynamics (MD) in the framework of a simulated annealing technique. In the spirit of the experimentally constrained molecular relaxation technique [Biswas et al., Phys. Rev. B 69, 195207 (2004)], this modified procedure offers a refined strategy in the field of reconstruction techniques, with special interest for heterogeneous and disordered solids such as amorphous porous materials. While the HRMC method generates physical structures, thanks to the use of energy penalties, the combination with MD makes the method at least one order of magnitude faster than HRMC simulations to obtain structures of similar quality. Furthermore, in order to ensure the transferability of this technique, we provide rational arguments to select the various input parameters such as the relative weight of the energy penalty with respect to the structure optimization. By applying the method to disordered porous carbons, we show that adsorption properties provide data to test the global texture of the reconstructed sample but are only weakly sensitive to the presence of defects. In contrast, the vibrational properties such as the phonon density of states are found to be very sensitive to the local structure of the sample.

Image of FIG. 1.

Image of FIG. 2.

Image of FIG. 3.

Image of FIG. 4.

Image of FIG. 5.

Image of FIG. 6.

Image of FIG. 7.

}, issn = {0021-9606}, doi = {10.1063/1.4914921}, author = {Bousige, Colin and Botan, Alexandru and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @article {141, title = {Organic{\textendash}Clay Interfacial Chemical Bonds Probed by ab Initio Calculations}, journal = {The Journal of Physical Chemistry C}, volume = {119}, year = {2015}, month = {Mar-26-2015}, pages = {6511 - 6517}, abstract = {
Abstract Image

Evaluating the nature of binding of the organic phase (called kerogen) with clays and other minerals in gas and oil shale is of critical importance for developing optimal processes for hydrocarbon extraction and recovery. Here we report the results of an ab initio reactive study of a gaseous oil/clay interface, which consists of grafting fragments of sp3 methane and sp2 benzene onto the basal surface of illite. Methane and benzene fragments were selected as they correspond to the simplest monomers of the immature (aliphatic sp3 carbon atoms) and mature (aromatic sp2 carbon atoms) stages of the kerogen macromolecule while illite is encountered in many shale formations. We find that the methyl or phenyl radicals always bond upright to one of the outmost Si atoms of the clay phase, while the H atom attaches to the inner Al\–OH group forming a water molecule. The covalent attachment of the radicals to the oxygen atoms of the clay surface is always less favorable than to the silicon atoms. The energetics of these grafted adducts are discussed along with the maturation state of the shale. An unanticipated signature of the covalent Si\–C binding is the formation of hypervalent Si complexes, characterized by elongated bonds at the oil/clay interface. We complement our atomistic modeling by simulating IR spectra of the two covalently grafted interfaces (methyl\–illite and phenyl\–illite) which fingerprint elongated Si\–C bonds in the frequency range 1100\–1250 cm\–1. This novel finding should have a large impact on the understanding of fracking conditions in shales through this tight- and loose-bonded organic/inorganic interface.

}, issn = {1932-7447}, doi = {10.1021/jp510885f}, author = {Bocquet, Marie-Laure and Benoit A. Coasne and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {138, title = {Physical Origins of Thermal Properties of Cement Paste}, journal = {Physical Review Applied}, volume = {3}, year = {2015}, month = {Jun-17-2015}, pages = {Article Number: 064010}, abstract = {

Despite the ever-increasing interest in multiscale porous materials, the chemophysical origin of their thermal properties at the nanoscale and its connection to the macroscale properties still remain rather obscure. In this paper, we link the atomic- and macroscopic-level thermal properties by combining tools of statistical physics and mean-field homogenization theory. We begin with analyzing the vibrational density of states of several calcium-silicate materials in the cement paste. Unlike crystalline phases, we indicate that calcium silicate hydrates (CSH) exhibit extra vibrational states at low frequencies (\<2\ \ THz) compared to the vibrational states predicted by the Debye model. This anomaly is commonly referred to as the boson peak in glass physics. In addition, the specific-heat capacity of CSH in both dry and saturated states scales linearly with the calcium-to-silicon ratio. We show that the nanoscale-confining environment of CSH decreases the apparent heat capacity of water by a factor of 4. Furthermore, full thermal conductivity tensors for all phases are calculated via the Green-Kubo formalism. We estimate the mean free path of phonons in calcium silicates to be on the order of interatomic bonds. This satisfies the scale separability condition and justifies the use of mean-field homogenization theories for upscaling purposes. Upscaling schemes yield a good estimate of the macroscopic specific-heat capacity and thermal conductivity of cement paste during the hydration process, independent of fitting parameters.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, doi = {10.1103/PhysRevApplied.3.064010}, author = {Mohammad Javad Abdolhosseini Qomi and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {131, title = {Rigidity transition in materials: hardness is driven by weak atomic constraints.}, journal = {Phys Rev Lett}, volume = {114}, year = {2015}, month = {Mar-23-2015}, pages = {Article Number: 125502}, abstract = {

Understanding the composition dependence of the hardness in materials is of primary importance for infrastructures and handled devices. Stimulated by the need for stronger protective screens, topological constraint theory has recently been used to predict the hardness in glasses. Herein, we report that the concept of rigidity transition can be extended to a broader range of materials than just glass. We show that hardness depends linearly on the number of angular constraints, which, compared to radial interactions, constitute the weaker ones acting between the atoms. This leads to a predictive model for hardness, generally applicable to any crystalline or glassy material.

Figure

Figure

Figure

Figure

}, issn = {1079-7114}, doi = {10.1103/PhysRevLett.114.125502}, author = {Mathieu Bauchy and Mohammad Javad Abdolhosseini Qomi and Christophe Bichara and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @proceedings {329, title = {Shrinkage Due to Colloidal Force Interactions}, journal = {10th International Conference on Mechanics and Physics of Creep, Shrinkage, and Durability of Concrete and Concrete Structures CONCREEP 10}, volume = {CONCREEP 10: MECHANICS AND PHYSICS OF CREEP, SHRINKAGE, AND DURABILITY OF CONCRETE AND CONCRETE STRUCTURES}, year = {2015}, month = {Sep-17-2015}, pages = {13-16}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA }, abstract = {

It has long been argued that shrinkage in cementitious materials can occur even in the absence of drying phenomena. Herein we argue that such shrinkage is attributable to colloidal forces that develop between C-S-H particles during the hydration in an out-of-equilibrium process. We base our argument on meso-scale simulation results, and show via microporomechanics theory that such colloidal forces can lead to the build-up of solid eigenstresses causing macroscopic shrinkage deformation.

}, doi = {10.1061/978078447934610.1061/9780784479346.002}, author = {Franz-Josef Ulm and Roland Jean-Marc Pellenq}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @article {180, title = {Structure and Dynamics of Ionic Liquids Confined in Amorphous Porous Chalcogenides}, journal = {Langmuir}, volume = {31}, year = {2015}, month = {Jun-23-2015}, pages = {6742 - 6751}, abstract = {
Abstract Image

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

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

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

Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

}, issn = {2041-1723}, doi = {10.1038/ncomms7949}, author = {Kerstyn I. Falk and Benoit A. Coasne and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Lyd{\'e}ric Bocquet} } @proceedings {359, title = {Transient regime to fluidized chimney within a granular bed by means of a 2D DEM/LBM modeling}, journal = {IV International Conference on Particle-based Methods (PARTICLES 2015)}, volume = {PARTICLE-BASED METHODS IV-FUNDAMENTALS AND APPLICATIONS}, year = {2015}, month = {sept-30-2015}, pages = {558-566}, publisher = {International Center for Numerical Methods in Engineering (CIMNE)}, address = {SEP 28-30 2015 Barcelona, SPAIN}, abstract = {

Beyond a given threshold, an upward fluid flow at constant flowrate, injected through a small size section, is able to generate a fluidization along a vertical chimney over the entire height of a granular assembly. Fluidization is first initiated in the immediate vicinity of the injection hole and then the fluidized zone grows gradually until reaching the upper surface of the granular packing. In this work, we present numerical results on the kinetics of chimney fluidization in an immersed granular bed produced with two-dimensional simulations coupling the Discrete Element and Lattice Boltzmann Methods (DEM-LBM). A parametric study is carried out with 11 different sets of physical parameters and analyzed based on spatio-temporal diagrams. Then a dimensional analysis allows finding general scaling laws for both threshold and growth rate of the fluidized zone by use of two dimensionless numbers, namely Reynolds and Archimedes numbers, while quite simple empirical relationships can also be proposed.

}, keywords = {DEM, Fluidization, Granular Materials, LBM}, url = {https://hal.archives-ouvertes.fr/hal-01269325}, author = {Jeff Ngoma and Pierre Philippe and St{\'e}phane Bonelli and Cuellar, Pablo and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Onate, E and Bischoff, M and Owen, DRJ} } @article {102, title = {Anomalous composition-dependent dynamics of nanoconfined water in the interlayer of disordered calcium-silicates}, journal = {Journal of Chemical Physics}, volume = {140}, year = {2014}, month = {Feb-07-2014}, pages = {Article Number: 054515}, type = {Article}, abstract = {

With shear interest in nanoporous materials, the ultraconfining interlayer spacing of calcium-silicate-hydrate (C-S-H) provides an excellent medium to study reactivity, structure, and dynamic properties of water. In this paper, we present how substrate composition affects chemo-physical properties of water in ultraconfined hydrophilic media. This is achieved by performing molecular dynamics simulation on a set of 150 realistic models with different compositions of calcium and silicon contents. It is demonstrated that the substrate chemistry directly affects the structural properties of water molecules. The motion of confined water shows a multi-stage dynamics which is characteristic of supercooled liquids and glassy phases. Inhomogeneity in that dynamics is used to differentiate between mobile and immobile water molecules. Furthermore, it is shown that the mobility of water molecules is composition-dependent. Similar to the pressure-driven self-diffusivity anomaly observed in bulk water, we report the first study on composition-driven diffusion anomaly, the self diffusivity increases with increasing confined water density in C-S-H. Such anomalous behavior is explained by the decrease in the typical activation energy required for a water molecule to escape its dynamical cage. (C) 2014 AIP Publishing LLC.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
}, issn = {0021-9606}, doi = {10.1063/1.4864118}, author = {Mohammad Javad Abdolhosseini Qomi and Mathieu Bauchy and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @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 {84, title = {Atomic-scale modelling of elastic and failure properties of clays}, journal = {Molecular Physics}, volume = {112}, year = {2014}, month = {May-19-2014}, pages = {1294-1305}, type = {Article}, abstract = {

The elastic and failure properties of a typical clay, illite, are investigated using molecular simulation. We employ a reactive (ReaxFF) and a non-reactive (ClayFF) force field to assess the elastic properties of the clay. As far as failure is concerned, ReaxFF was used throughout the study; however, some calculations were also performed with ClayFF. A crack parallel to the clay layers is found to have low fracture resistance when submitted to a tensile loading perpendicular to the crack. The mechanism of both yield and fracture failures is decohesion in the interlayer space. In contrast, under shear loading, the nanoscale failure mechanism is a stick-slip between clay layers. No fracture propagation is observed as the clay layers slide on top of each other. The low fracture resistance in mode I and the stick-slip failure in mode II are both the consequence of the lack of chemical bonds between clay layers where the cohesion is provided by non-covalent interactions. This work, which provides a description of the failure of clays at the microscopic scale, is the first step towards describing the failure of clays at a larger scale where the polycrystalline distribution of clay grains must be taken into account.

}, keywords = {clay, elastic properties, fracture, reactive molecular simulation}, issn = {0026-8976}, doi = {10.1080/00268976.2014.897393}, author = {Gy{\"o}rgy Hantal and Brochard, Laurent and Hadrien Laubie and Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Benoit A. Coasne} } @article {39, title = {An atomistic modeling of the xenon bubble behavior in the UO2 matrix}, journal = {Journal of Nuclear Materials}, volume = {444}, year = {2014}, month = {Jan-2014}, pages = {153{\textendash}160}, abstract = {

The behavior of the xenon nanoinclusions/bubbles in the uranium dioxide (UO2) matrix and their influence on its swelling were investigated through atomistic simulation techniques. The pressure in bubbles of less than 2\ nm in diameter, calculated using a virial equation that takes into account the xenon/matrix interactions, is larger than the pressure calculated in simulations of the equivalent density and temperature of super critical bulk xenon. The radial distribution function of confined xenon is characteristic of a dense (\ \>\ 4\ g/cm3) glassy phase. The swelling of the UO2 induced by the intragranular bubbles is proportional to the Xe/U ratio but independent of the temperature.

The UO2 elementary cell. It contains 4 uranium atoms (large spheres) forming a ...

Comparison between the density versus pressure ρ(P) curves calculated at 298K ...

Comparison between the bulk modulus versus pressure K(P) curves calculated at ...

}, doi = {10.1016/j.jnucmat.2013.09.041}, author = {Andrei Jelea and Roland Jean-Marc Pellenq and F. Ribeiro} } @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 {40, title = {Controlling local packing and growth in calcium-silicate-hydrate gels}, journal = {Soft Matter}, volume = {10}, year = {2014}, month = {Oct-2014}, pages = {1121-1133}, abstract = {

We investigate the development of gels under out-of-equilibrium conditions, such as calcium\–silicate\–hydrate (C\–S\–H) gels that form during cement hydration and are the major factor responsible for cement mechanical strength. We propose a new model and numerical approach to follow the gel formation upon precipitation and aggregation of nano-scale colloidal hydrates, whose effective interactions are consistent with forces measured in experiments at fixed lime concentrations. We use Grand Canonical Monte Carlo to mimic precipitation events during Molecular Dynamics simulations, with their rate corresponding to the hydrate production rate set by the chemical environment. Our results display hydrate precipitation curves that indeed reproduce the acceleration and deceleration regime typically observed in experiments and we are able to correctly capture the effect of lime concentration on the hydration kinetics and the gel morphology. Our analysis of the evolution of the gel morphology indicates that the acceleration is related to the formation of an optimal local crystalline packing that allows for large, elongated aggregates to grow and that is controlled by the underlying thermodynamics. The defects produced during precipitation favor branching and gelation that end up controlling the deceleration. The effects on the mechanical properties of C\–S\–H gels are also discussed.

Graphical abstract: Controlling local packing and growth in calcium{\textendash}silicate{\textendash}hydrate gels
}, doi = {10.1039/C3SM52232F}, author = {Katerina Ioannidou and Roland Jean-Marc Pellenq and Emanuela Del Gado} } @article {37, title = {Docking 90 Sr radionuclide in cement: An atomistic modeling study}, journal = {Physics and Chemistry of the Earth, Parts A/B/C}, volume = {70{\textendash}71}, year = {2014}, month = {2014}, pages = {39{\textendash}44}, abstract = {

Cementitious materials are considered to be a waste form for the ultimate disposal of radioactive materials in geological repositories. We investigated by means of atomistic simulations the encapsulation of strontium-90, an important radionuclide, in calcium\–silicate\–hydrate (C\–S\–H) and its crystalline analog, the 9\ \Å-tobermorite. C\–S\–H is the major binding phase of cement. Strontium was shown to energetically favor substituting calcium in the interlayer sites in C\–S\–H and 9\ \Å-tobermorite with the trend more pronounced in the latter. The integrity of the silicate chains in both cementitious waste forms were not affected by strontium substitution within the time span of molecular dynamics simulation. Finally, we observed a limited degradation of the mechanical properties in the strontium-containing cementitious waste form with the increasing strontium concentration. These results suggest the cement hydrate as a good candidate for immobilizing radioactive strontium.

Image for unlabelled figure

}, doi = {10.1016/j.pce.2013.11.007}, author = {Youssef, Mostafa and Roland Jean-Marc Pellenq and Bilge Yildiz} } @article {97, title = {Docking Sr-90 radionuclide in cement: An atomistic modeling study}, journal = {Physics and Chemistry of the Earth}, volume = {70-71}, year = {2014}, month = {2014}, pages = {39-44}, type = {Article}, abstract = {

Cementitious materials are considered to be a waste form for the ultimate disposal of radioactive materials in geological repositories. We investigated by means of atomistic simulations the encapsulation of strontium-90, an important radionuclide, in calcium silicate hydrate (C-S-H) and its crystalline analog, the 9 angstrom-tobermorite. C-S-H is the major binding phase of cement. Strontium was shown to energetically favor substituting calcium in the interlayer sites in C-S-H and 9 angstrom-tobermorite with the trend more pronounced in the latter. The integrity of the silicate chains in both cementitious waste forms were not affected by strontium substitution within the time span of molecular dynamics simulation. Finally, we observed a limited degradation of the mechanical properties in the strontium-containing cementitious waste form with the increasing strontium concentration. These results suggest the cement hydrate as a good candidate for immobilizing radioactive strontium. (C) 2013 Elsevier Ltd. All rights reserved.

}, keywords = {Cement, mechanical properties, Molecular simulation, Nuclear waste storage}, issn = {1474-7065}, doi = {10.1016/j.pce.2013.11.007}, author = {Youssef, Mostafa and Roland Jean-Marc Pellenq and Bilge Yildiz} } @article {101, title = {Effect of temperature on adsorption of mixtures in porous materials}, journal = {Molecular Simulation}, volume = {40}, year = {2014}, month = {Jan-02-2014}, pages = {45-51}, type = {Article}, abstract = {

The effect of temperature on the adsorption of a simple mixture (Ar/Kr) in disordered porous materials is investigated by means of molecular simulation. In the larger mesopores of porous silica glasses, capillary condensation occurs upon decreasing the temperature. At temperatures above the capillary condensation temperature, Kr is preferentially adsorbed at the pore surface and Ar adsorption occurs in regions of low Kr density. For temperatures below the capillary condensation temperature, Ar density surprisingly increases as temperature increases, the behaviour that is consistent with an over-solubility effect. In contrast, in the disordered sub-nanoporous carbon, filling of the pores occurs in a reversible and continuous way upon decreasing the temperature, owing to the small size and amorphous shape of the pores. These results show that the crossover between capillary condensation and continuous reversible filling observed for pure fluids in pores also exists for mixtures. We also show that the Kr selectivity exhibits a minimum in the disordered porous silica that is located at the capillary condensation temperature. In contrast, in the disordered porous carbon where no capillary condensation occurs, the selectivity decreases monotonically with increasing the temperature. These results shed light on low-temperature adsorption of mixtures confined in porous materials and provide a guide to design efficient phase separation processes.

}, keywords = {Adsorption, grand canonical Monte Carlo simulations, Mixtures, porous materials}, issn = {0892-7022}, doi = {10.1080/08927022.2013.852192}, author = {Kerstyn I. Falk and Benoit A. Coasne 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 {85, title = {Enhanced H-2 Uptake of n-Alkanes Confined in Mesoporous Materials}, journal = {Journal of Physical Chemistry C }, volume = {118}, year = {2014}, month = {Apr-30-2014}, pages = {10720-10727}, type = {Article}, abstract = {
Abstract Image

The hydrogen uptake in hybrid sorbents consisting of n-alkane solvents confined in mesoporous silica aerogel is measured at different temperatures from 273 to 313 K and pressures up to 40 bar. An apparent \“oversolubility\” effect is observed as the H2 uptake in the hybrid sorbents is much larger than that in bulk solvents. The H2 uptake in the hybrid sorbents is found to increase with increasing temperature, which suggests that the flexibility and conformation of n-alkane molecules confined in the aerogel play a crucial role; high-entropy (disordered) alkane configurations lead to the creation of numerous cavities which make it possible to solubilize a larger number of H2 molecules. This departs from adsorption-driven solubility effects for which the number of solubilized molecules decreases with increasing temperature. For a given temperature and pressure, it is found that the number of solubilized H2 molecules per unit volume increases with decreasing alkane chain length. Such an effect, which is observed for both the bulk alkanes and the alkanes confined in the silica aerogel, can be rationalized by considering the number density of CHx (x = 2 or 3) groups; for a given temperature, the latter number density decreases with decreasing alkane chain length so that the free volume available to solubilize H2 molecules increases.

}, issn = {1932-7447}, doi = {10.1021/jp411526f}, author = {St{\'e}phanie Clauzier and Linh Ngoc Ho and Pera-Titus, M. and David Farrusseng and Benoit A. Coasne} } @article {328, title = {Erratum to {\textquotedblleft}Static states of cohesive granular media{\textquotedblright}}, journal = {Journal of Mechanical Science and Technology}, volume = {28}, year = {2014}, month = {Oct-2014}, pages = {4345 - 4345}, abstract = {

vol 21, pg 1957, 2007

This paper investigates the influence of cohesion on the behavior of 2D granular media simulated by using a molecular dynamics method, involving a simple contact law with adhesion. The study considers an adhesion index which is non intrinsic but significant for interpreting the mechanical effects of applied pressure. A static state system of 4000 particles is numerically simulated by oedometrical compression under a constant force, without gravity. The results show how the geometrical texture and the network of contact forces change according to the level of adhesion. This effect is explained essentially in term of the increase in the number of the tensile contacts and by the development of the internal self-stress structure. As in the case of non cohesive granular media, a high spatial heterogeneity of the contact forces is observed.

}, issn = {1738-494X}, doi = {10.1007/s12206-014-0951-3}, author = {Preechawuttipong, I. and Peyroux, R. and Farhang Radja{\"\i} and Rangsri, W.} } @proceedings {346, title = {Geomechanics from Micro to MacroInteraction between two localized fluidization cavities in granular media : Experiments and numerical simulation}, journal = {3rd International Symposium on Geomechanics from Micro to Macro}, volume = {Geomechanics from Micro to Macro}, year = {2014}, month = {Sep-2014}, pages = {1571 - 1576}, publisher = {CRC Press}, address = {SEP 01-03-2014 Univ Cambridge, Cambridge, ENGLAND}, abstract = {

In this work, we present experimental and numerical results on the interaction between two localized fluidization cavities in an immersed granular packing. According to the gap between the two locally injected upward fluid flows, each fluidized cavity will evolve independently of the another, or conversely, the two cavities can interact with each other: they merge and eventually behave like a single cavity. Combined optical techniques are used for visualization of particle motion inside the granular media (Refractive Index-Matching between liquid and beads and Planar Laser-Induced Fluorescence), the experimental results are compared to a two-dimensional simulation based on coupled Discrete Element and Lattice Boltzmann Methods (DEM-LBM).

}, isbn = {978-1-138-02707-7}, doi = {10.1201/b17395-286}, author = {Jeff Ngoma and Pierre Philippe and St{\'e}phane Bonelli and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Kenichi Soga and Krishna Kumar and Giovanna Biscontin and Kuo, Matthew} } @article {113, title = {Mesoscale properties of clay aggregates from potential of mean force representation of interactions between nanoplatelets}, journal = {Journal of Chemical Physics }, volume = {140}, year = {2014}, month = {Apr-21-2014}, pages = {Article Number: 154309}, type = {Article}, abstract = {

Face-to-face and edge-to-edge free energy interactions of Wyoming Na-montmorillonite platelets were studied by calculating potential of mean force along their center to center reaction coordinate using explicit solvent (i.e., water) molecular dynamics and free energy perturbation methods. Using a series of configurations, the Gay-Berne potential was parametrized and used to examine the mesoscale aggregation and properties of platelets that are initially random oriented under isothermal-isobaric conditions. Aggregates of clay were defined by geometrical analysis of face-to-face proximity of platelets with size distribution described by a log-normal function. The isotropy of the microstructure was assessed by computing a scalar order parameter. The number of platelets per aggregate and anisotropy of the microstructure both increases with platelet plan area. The system becomes more ordered and aggregate size increases with increasing pressure until maximum ordered state at confining pressure of 50 atm. Further increase of pressure slides platelets relative to each other leading to smaller aggregate size. The results show aggregate size of (3-8) platelets for sodium-smectite in agreement with experiments (3-10). The geometrical arrangement of aggregates affects mechanical properties of the system. The elastic properties of the meso-scale aggregate assembly are reported and compared with nanoindentation experiments. It is found that the elastic properties at this scale are close to the cubic systems. The elastic stiffness and anisotropy of the assembly increases with the size of the platelets and the level of external pressure. (C) 2014 AIP Publishing LLC.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
Image of FIG. 4.
}, issn = {0021-9606}, doi = {10.1063/1.4870932}, author = {Ebrahimi, Davoud and Andrew J. Whittle and Roland Jean-Marc Pellenq} } @proceedings {251, title = {Mesoscale simulation of clay aggregate formation and mechanical properties}, journal = {3rd International Symposium on Geomechanics from Micro to Macro}, volume = {Geomechanics from Micro to Macro: Granular Matter}, year = {2014}, month = {Sep-2014}, pages = {539-544}, address = {SEP 01-03 2014 Univ Cambridge, Cambridge, ENGLAND}, abstract = {

This paper proposes a novel methodology for understanding the meso-scale aggregation of clay platelets in water. We use Molecular Dynamics simulations using the CLAYFF force fields to represent the interactions between two layers of Wyoming montmorillonite (Na-smectite) in bulk water. The analyses are used to establish the potential of mean force at different spacings between the layers for edge-to-edge and face-to-face interactions. This is accomplished by finding the change in free energy as a function of the separation distance between the platelets using thermodynamic perturbation theory with a simple overlap sampling method. These nanoscale results are then used to calibrate the Gay\–Berne (GB) potential that represents each platelet as a single-site ellipsoidal body. A coarse-graining upscaling approach then uses the GB potentials and molecular dynamics to represent the meso-scale aggregation of clay platelets (at submicron length scale). Results from meso-scale simulations obtain the equilibrium/jamming configurations for mono-disperse clay platelets. The results show aggregation for a range of clay platelets dimensions and pressures with mean stack size ranging from 3 to 8 platelets. The particle assemblies become more ordered and exhibit more pronounced elastic anisotropy at higher confining pressures. The results are in good agreement with previously measured nano-indentation moduli over a wide range of clay packing densities.

}, issn = {1434-5021}, doi = {10.1007/s10035-016-0655-8}, author = {Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Andrew J. Whittle}, editor = {Kenichi Soga and Krishna Kumar and Giovanna Biscontin} } @proceedings {99, title = {Modelling cement at fundamental scales: From atoms to engineering strength and durability}, journal = {Euro-C Conference}, volume = {COMPUTATIONAL MODELLING OF CONCRETE STRUCTURES, VOL 1}, year = {2014}, month = {2014}, pages = {139-148}, publisher = {CRC Press-Taylor \& Francis Group}, address = {MAR 24-27, 2014, St Anton am Alberg, AUSTRIA}, abstract = {

The mechanics and rheology of cement pastes are largely determined by the Calcium-Silicate-Hydrate (C-S-H) gel, which is normally the main binding phase. The C-S-H gel is an amorphous material with a nanoscale pore network. Here we consider a model structure of the gel, where a continuous molecular model of non-porous C-S-H is combined with a colloidal model that incorporates the porosity of the gel. Numerical simulations predict mechanical properties and microscopic processes associated with yielding under shear. The proposed model can now be used to extract further properties and to contribute to a broader effort aimed at describing and engineering the multi-scale structure and mechanical behaviour of the cement paste. We finally discuss the opportunities offered by our model with respect to the scientific and technological challenges of modelling and controlling ageing and creep over long time scales.

}, isbn = {978-1-138-02641-4; 978-1-315-76203-6}, doi = {10.1201/b16645-16}, author = {Enrico Masoero and Hamlin M. Jennings and Franz-Josef Ulm and Emanuela Del Gado and Hegoi Manzano and Roland Jean-Marc Pellenq and Sidney Yip}, editor = {Bicanic, N and Mang, H and Meschke, G and DeBorst, R} } @article {96, title = {Nano-scale mechanics of colloidal C-S-H gels}, journal = {Soft Matter}, volume = {10}, year = {2014}, month = {Jan-21-2014}, pages = {491-499}, type = {Article}, abstract = {

Gels of calcium-silicate-hydrates (C-S-H) are the glue that is largely responsible for the mechanical properties of cement. Despite their practical relevance, their nano-scale structure and mechanics are still mainly unexplored, because of the difficulties in characterizing them in a complex material like cement. We propose a colloidal model to investigate the gel mechanics emerging in the critical range of length-scales from several tens to hundreds of nanometers. We show that the size polydispersity of the hydrates and size-dependent effective interactions can explain the mechanical heterogeneities detected in nano-indentation experiments. We also show how these features control the arising of irreversible structural rearrangements under deformation, which are good candidates as nano-scale mechanisms underlying mechanical aging and slow structural relaxation in the gels.

Graphical abstract: Nano-scale mechanics of colloidal C{\textendash}S{\textendash}H gels
}, issn = {1744-683X}, doi = {10.1039/c3sm51815a}, author = {Enrico Masoero and Emanuela Del Gado and Roland Jean-Marc Pellenq and Sidney Yip and Franz-Josef Ulm} } @article {115, title = {Nanoscale Structure of Cement: Viewpoint of Rigidity Theory}, journal = {Journal of Physical Chemistry C}, volume = {118}, year = {2014}, month = {Jun-12-2014}, pages = {12485-12493}, type = {Article}, abstract = {
Abstract Image

Rigidity theory is a powerful tool to predict the properties of glasses with respect to composition. By reducing such molecular networks to simple mechanical trusses, topological constraint theory filters out all the unnecessary details that ultimately do not affect macroscopic properties. However, the usual constraint enumeration is restricted to networks that are amorphous, homogeneous, and fully connected. On the contrary, calcium\–silicate\–hydrate (C\–S\–H), the binding phase of cement, is partially crystalline and heterogeneous and shows some isolated water molecules. Here, we report how rigidity theory can be used to describe the nanoscale structure of this material by relying on molecular dynamics simulations. The distinction between intact and broken constraints is clearly defined at the atomic scale, thus allowing a precise enumeration of the topological constraints. We show that the rigidity of the C\–S\–H network can be increased by decreasing the Ca/Si molar ratio, which, as predicted by rigidity theory, allows improvement of the hardness of the material. This study suggests that rigidity theory could be applied with great rewards to a broader range of materials than glasses.

}, issn = {1932-7447}, doi = {10.1021/jp502550z}, author = {Mathieu Bauchy and Mohammad Javad Abdolhosseini Qomi and Christophe Bichara and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {262, title = {Order and disorder in calcium{\textendash}silicate{\textendash}hydrate}, journal = {The Journal of Chemical Physics}, volume = {140}, year = {2014}, month = {Jun-07-2014}, pages = {Article Number: 214503}, abstract = {

Despite advances in the characterization and modeling of cement hydrates, the atomic order in Calcium\–Silicate\–Hydrate (C\–S\–H), the binding phase of cement, remains an open question. Indeed, in contrast to the former crystalline model, recent molecular models suggest that the nanoscale structure of C\–S\–H is amorphous. To elucidate this issue, we analyzed the structure of a realistic simulated model of C\–S\–H, and compared the latter to crystalline tobermorite, a natural analogue of C\–S\–H, and to an artificial ideal glass. The results clearly indicate that C\–S\–H appears as amorphous, when averaged on all atoms. However, an analysis of the order around each atomic species reveals that its structure shows an intermediate degree of order, retaining some characteristics of the crystal while acquiring an overall glass-like disorder. Thanks to a detailed quantification of order and disorder, we show that, while C\–S\–H retains some signatures of a tobermorite-like layered structure, hydrated species are completely amorphous.

}, issn = {0021-9606}, doi = {10.1063/1.4878656}, author = {Mathieu Bauchy and Mohammad Javad Abdolhosseini Qomi and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @proceedings {100, title = {Poro-chemo-fracture-mechanics ... bottom-up: Application to risk of fracture design of oil and gas cement sheath at early ages}, journal = {Euro-C Conference}, volume = {COMPUTATIONAL MODELLING OF CONCRETE STRUCTURES, VOL 1 }, year = {2014}, month = {Jan-2014}, pages = {61-67}, publisher = {CRC PRESS-TAYLOR \& FRANCIS GROUP}, address = {MAR 24-27, 2014, St Anton am Alberg, AUSTRIA}, abstract = {

With ever more challenging (T, p) environments for cementing applications in oil and gas wells, there is a need to change the classical paradigm of constitutive modeling of the early-age behavior of cementitious materials in oil-and gas well applications. Herein, we propose a bottom-up approach, which starts at the molecular scale of calcium-silicate-hydrates (C-S-H), and propagates this understanding via incremental poroelastic upscaling methods to the meso-and macroscale, and ultimately to the scale of engineering applications. With a clear focus on engineering risk-of-fracture evaluations, we show how to integrate this bottom-up approach with fracture mechanics.

}, isbn = {978-1-138-02641-4; 978-1-315-76203-6}, doi = {10.1201/b16645-8}, author = {Franz-Josef Ulm and Muhannad Abuhaikal and Thomas Alexander Petersen and Roland Jean-Marc Pellenq}, editor = {Bicanic, N and Mang, H and Meschke, G and DeBorst, R} } @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 {261, title = {A soft matter in construction {\textendash} Statistical physics approach to formation and mechanics of C{\textendash}S{\textendash}H gels in cement}, journal = {The European Physical Journal Special Topics}, volume = {223}, year = {2014}, month = {Oct-2014}, pages = {2285 - 2295}, abstract = {

Calcium-silicate hydrate (C-S-H) is the main binding agent in cement and concrete. It forms at the beginning of cement hydration, it progressively densifies as cement hardens and is ultimately responsible of concrete performances. This hydration product is a cohesive nano-scale gel, whose structure and mechanics are still poorly understood, in spite of its practical importance. Here we review some of the open questions for this fascinating material and a statistical physics approach recently developed, which allows us to investigate the gel formation under the out-of-equilibrium conditions typical of cement hydration and the role of the nano-scale structure in C-S-H mechanics upon hardening. Our approach unveils how some distinctive features of the kinetics of cement hydration can be related to changes in the morphology of the gels and elucidates the role of nano-scale mechanical heterogeneities in the hardened C-S-H.

}, issn = {1951-6355}, doi = {10.1140/epjst/e2014-02264-1}, author = {Emanuela Del Gado and Katerina Ioannidou and Enrico Masoero and Alain Baronnet and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Sidney Yip} } @article {32, title = {Surface Chemistry and Atomic-Scale Reconstruction of Kerogen-Silica Composites}, journal = {The Journal of Physical Chemistry C}, volume = {118}, year = {2014}, month = {Feb-6-2014}, pages = {2429-2438}, abstract = {
Abstract Image

Interest in gas shale, a novel source rock of natural gas, has increased tremendously in recent years. Better understanding of the kerogen\–rock interaction is of crucial importance for efficient gas extraction and, hence, asset management. In this study, we explore the possible chemical bonds between kerogen and silica, one of the most predominant mineral constituents of gas shale, by means of quantum chemistry. Energetically favorable bond formation reactions are found between alcoholic hydroxyl, carboxylate, and aldehyde groups, as well as aliphatic double bonds of kerogen and the silica surface. The performance of a reactive force field was also assessed in a representative set of chemical reactions and found to be satisfactory. The potential impact of bond formation reactions between the two phases on the actual kerogen\–silica interface is discussed as a function of the kerogen type, maturity, and density. Finally, a methodology aiming to reconstruct realistic kerogen\–silica interfaces is presented.

}, doi = {10.1021/jp406329n}, author = {Gy{\"o}rgy Hantal and Brochard, Laurent and Dias Soeiro Cordeiro, M Natalia and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @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 {95, title = {Adsorption, intrusion and freezing in porous silica: the view from the nanoscale}, journal = {Chemical Society Reviews}, volume = {42}, year = {2013}, month = {Jan-24-2013}, pages = {4141-4171}, type = {Review}, abstract = {

This review presents the state of the art of molecular simulation and theory of adsorption, intrusion and freezing in porous silica. Both silica pores of a simple geometry and disordered porous silicas which exhibit morphological and topological disorders are considered. We provide a brief description of the numerical models of porous silicas available in the literature and present the most common molecular simulation and theoretical methods. Adsorption in regular and irregular pores is discussed in the light of classical theories of adsorption and capillary condensation in pores. We also present the different evaporation mechanisms for disordered systems: pore blocking and cavitation. The criticality of fluids confined in pores, which is still the matter of debate, is then discussed. We review theoretical results for intrusion/extrusion and freezing in silica pores and discuss the validity of classical approaches such as the Washburn-Laplace equation and Gibbs-Thomson equation to describe the thermodynamics of intrusion and in-pore freezing. The validity of the most widely used characterization techniques is then discussed. We report some concluding remarks and suggest directions for future work.

Graphical abstract: Adsorption, intrusion and freezing in porous silica: the view from the nanoscale

}, issn = {0306-0012}, doi = {10.1039/c2cs35384a}, author = {Benoit A. Coasne and Anne Galarneau and Roland Jean-Marc Pellenq and Di Renzo, Francesco} } @article {35, title = {Adsorption of Aromatic Hydrocarbon Molecules at the Surface of Ice, As Seen by Grand Canonical Monte Carlo Simulation}, journal = {The Journal of Physical Chemistry C}, volume = {117}, year = {2013}, month = {Apr-2013}, pages = {6719{\textendash}6729}, abstract = {
Abstract Image

The adsorption of four aromatic hydrocarbon compounds, benzene, naphthalene, anthracene, and phenanthrene, at the surface of Ih ice is investigated by grand canonical Monte Carlo (GCMC) computer simulation under tropospheric conditions at 200 K. By systematic variation of the value of adsorbate chemical potential in the simulations, the adsorption isotherms are determined. It is found that adsorption follows the Langmuir mechanism only up to a rather low relative pressure value in every case. In this range specific surface sites, called \α sites, to which adsorbate molecules can be bound particularly strongly in specific orientation, are occupied. In these \α sites, presumably the dangling OH bonds of the ice surface form O\–H....\π-type hydrogen bonds with the delocalized \π electrons of the adsorbed aromatic molecule lying parallel with the ice surface. Once these \α sites are saturated, lateral interactions become increasingly important, leading to large fluctuations of the lateral density of the adsorption layer and an increasing deviation of the adsorption isotherm from the Langmuir shape. The adsorption layer is found to be strictly monomolecular and even unsaturated in every case, as condensation well precedes the saturation of this monolayer for all four aromatic adsorbates considered in this study.

}, doi = {10.1021/jp401532x}, author = {M{\'e}sz{\'a}r, Zsuzsanna E and Gy{\"o}rgy Hantal and Picaud, Sylvain and Jedlovszky, P{\'a}l} } @article {24, title = {Adsorption of Volatile Organic Compounds in pure silica CHA, *BEA, MFI and STT-type zeolites}, journal = {Microp. Mesop. Mater.}, volume = {173}, year = {2013}, month = {Jun-2013}, pages = {147{\textendash}154}, abstract = {

Four pure silica zeolites, chabazite (CHA-structure type) and SSZ-23 (STT-structure type) with cage-like structure and silicalite-1 (MFI-structure type) and beta (\∗BEA-structure type) with channel structure, were synthesized and fully characterized. Their sorption properties were examined using gravimetric method combined with Grand Canonical Monte Carlo simulations (GCMC). Of particular interest is the large difference in the adsorption rates of n-hexane, p-xylene and acetone observed for these zeosils at 25, 75 and 150\ \°C. As expected, in most cases, a decrease in the adsorption capacity is observed with increasing the molar volume of the used probe molecule and the temperature. An exception is observed for the pure silica CHA-type zeolite due to its small pore size which prevents molecules from entering its porosity. However, at higher temperature, the window size widens slightly and allows n-hexane and acetone, which have kinetic diameter close to the pore opening, to enter easily in chabazite. As a result, for this zeosil, an increase in the adsorption capacity compared to that obtained at 25\ \°C is observed. All these zeosils are promising for technological uses in car exhaust gas decontamination.

Image for unlabelled figure

}, doi = {10.1016/j.micromeso.2013.02.009}, author = {A. F. Cosseron and T. J. Daou and L. Tzanis and H. Nouali and I. Deroche and Benoit A. Coasne and V. Tchaber and J. Patarin} } @proceedings {366, title = {AIP Conference ProceedingsTumbling sandpiles in a fluid}, 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 = {Jul-2013}, pages = {662-665}, publisher = {AIP}, address = {JUL 08-12 2013 Sydney, AUSTRALIA}, abstract = {

By means of contact dynamics simulations interfaced with computational fluid dynamics, we analyze the effect of a suspending fluid on the dynamics of collapse and spread of a granular column. We find that the runout distance increases as a power law with the aspect ratio of the column and, for a given aspect ratio, it may be the same in the grain-inertial and fluid-inertial regimes but with considerably longer duration in the latter case. We show that, in both viscous and fluid-inertial regimes, this behavior results from compensation between two effects of the fluid: 1) reduction of the kinetic energy during collapse and 2) enhancement of the flow by lubrication during spread.

}, doi = {10.1063/1.4812018}, author = {Farhang Radja{\"\i} and Vincent Topin and Perales, F. and Yann Monerie}, 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} } @proceedings {370, title = {GRANULAR SLUMPING IN A FLUID : FOCUS ON RUNOUT DISTANCES}, journal = {3rd International Conference on Particle-based Methods}, volume = {PARTICLE-BASED METHODS III: FUNDAMENTALS AND APPLICATIONS}, year = {2013}, month = {Sept-2013}, pages = {115-123}, address = {SEP 18-20 2013 Stuttgart, GERMANY}, author = {Vincent Topin and Yann Monerie and Perales, F.}, editor = {Bischoff, M and Ramm, E and Onate, E} } @article {284, title = {An interdisciplinary approach towards improved understanding of soil deformation during compaction}, journal = {Soil and Tillage Research}, volume = {128}, year = {2013}, month = {Apr-2013}, pages = {61 - 80}, abstract = {

Soil compaction not only reduces available pore volume in which fluids are stored, but it alters the arrangement of soil constituents and pore geometry, thereby adversely impacting fluid transport and a range of soil ecological functions. Quantitative understanding of stress transmission and deformation processes in arable soils remains limited. Yet such knowledge is essential for better predictions of effects of soil management practices such as agricultural field traffic on soil functioning. Concepts and theory used in agricultural soil mechanics (soil compaction and soil tillage) are often adopted from conventional soil mechanics (e.g. foundation engineering). However, in contrast with standard geotechnical applications, undesired stresses applied by agricultural tyres/tracks are highly dynamic and last for very short times. Moreover, arable soils are typically unsaturated and contain important secondary structures (e.g. aggregates), factors important for affecting their soil mechanical behaviour. Mechanical processes in porous media are not only of concern in soil mechanics, but also in other fields including geophysics and granular material science. Despite similarity of basic mechanical processes, theoretical frameworks often differ and reflect disciplinary focus. We review concepts from different but complementary fields concerned with porous media mechanics and highlight opportunities for synergistic advances in understanding deformation and compaction of arable soils. We highlight the important role of technological advances in non-destructive measurement methods at pore (X-ray tomography) and soil profile (seismic) scales that not only offer new insights into soil architecture and enable visualization of soil deformation, but are becoming instrumental in the development and validation of new soil compaction models. The integration of concepts underlying dynamic processes that modify soil pore spaces and bulk properties will improve the understanding of how soil management affect vital soil mechanical, hydraulic and ecological functions supporting plant growth.

}, issn = {01671987}, doi = {10.1016/j.still.2012.10.004}, author = {Keller, T. and Lamand, M. and Peth, S. and Berli, M. and Jean-Yves Delenne and Baumgarten, W. and Rabbel, W. and Farhang Radja{\"\i} and Rajchenbach, J. and Selvadurai, A.P.S. and Or, D.} } @article {25, title = {Large pressure enhancements in confined phases: effect of surface chemistry and roughness}, journal = {J. Chem. Phys.}, volume = {139}, year = {2013}, month = {Oct-10-2013}, chapter = {144701}, abstract = {

Experiments and simulations both suggest that the pressure experienced by an adsorbed phase confined within a carbon nanoporous material can be several orders of magnitude larger than the bulk phase pressure in equilibrium with the system. To investigate this pressure enhancement, we report a molecular-simulation study of the pressure tensor of argon confined in slit-shaped nanopores with walls of various models, including carbon and silica materials. We show that the pressure is strongly enhanced by confinement, arising from the effect of strongly attractive wall forces; confinement within purely repulsive walls does not lead to such enhanced pressures. Simulations with both the Lennard-Jones and Barker-Fisher-Watts intermolecular potentials for argon-argon interactions give rise to similar results. We also show that an increase in the wall roughness significantly decreases the in-pore pressure due to its influence on the structure of the adsorbate. Finally, we demonstrate that the pressures calculated from the mechanical (direct pressure tensor calculations) and the thermodynamic (volume perturbation method) routes yield almost identical results, suggesting that both methods can be used to calculate the local pressure tensor components in the case of these planar geometries.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
Image of FIG. 4.
}, doi = {10.1063/1.4824125}, author = {Yun Long and Jeremy C. Palmer and Benoit A. Coasne and Malgorzata Sliwinska-Bartkowiak and George Jackson and Erich A. Mueller and Keith E. Gubbins} } @article {23, title = {Mesoscopic Monte Carlo simulations of microstructure and conductivity of zirconia-based ceramics}, journal = {J. Eur. Ceram. Soc.}, volume = {61}, year = {2013}, month = {Aug-2013}, pages = {5018-5025}, abstract = {

Zirconia-based refractories are ceramics that serve as furnaces for the production of high-quality speciality glasses. These materials are obtained by cooling a melt composed of different oxides and oxide precursors. The microstructure of these refractories consists of monoclinic zirconia embedded in a three-dimensional interconnected amorphous phase. Using mesoscopic Monte Carlo simulations, the present paper addresses the effect of the amount of amorphous phase and temperature on the microstructure of the material. The simulated microstructures resemble those obtained using tomography experiments. A theoretical model describing the electrical behavior of the material as a function of temperature and composition is also reported. In this model, the conductivity of the ZrO2 particles is assumed to be constant, but the overall conductivity of the sample depends on its tortuosity s. Comparison with experimental data suggests that the present model provides a realistic picture of the electrical behavior of zirconia-based ceramics and leads to quantitative predictions.

Image for unlabelled figure

}, doi = {10.1016/j.actamat.2013.04.044}, author = {P. Epicoco and Benoit A. Coasne and A. Gioia and P. Papet and I. Cabodi and M. Gaubil} } @article {76, title = {Molecular intermittent dynamics of interfacial water: probing adsorption and bulk confinement}, journal = {Soft Matter}, volume = {9}, year = {2013}, month = {Jul-26-2013}, pages = {8654-8663}, type = {Article}}, abstract = {

Numerous natural and manufactured systems such as colloidal suspensions, geological pore networks, catalysts, and nanofluidic devices develop a large and sometimes complex interface strongly influencing the dynamics of the fluid entrapped inside these materials. A coarse grain picture of this molecular dynamics can be considered as an intermittence of adsorption steps and bulk relocations from one point to another point of the interface. Adsorption statistics such as the adsorption time distribution and its first moment reflect the degree of interaction of the molecule with the colloidal interface. Relocation statistics strongly depend on the shape of the pore, the surface forces and the bulk confinement. In this paper, a theoretical analysis of this intermittent dynamics is presented. A direct comparison with molecular dynamics simulations is proposed in the case of liquid water confined inside a hydrophilic substrate (silica slit pore with hydroxylated surfaces) or inside a hydrophobic substrate (carbon nanotube). Analysis of this intermittent dynamics allows quantification of the level of interaction of the vicinal water with the solid interface inside an independent adsorption region in exchange with the confined bulk fluid. The possibility of experimentally probing this dynamics using NMR relaxometry is emphasized.

}, issn = {1744-683X}, doi = {10.1039/c3sm51940f}, author = {Pierre E. Levitz and Patrick A Bonnaud and Pierre-Andre Cazade and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @article {103, title = {On the molecular origin of high-pressure effects in nanoconfinement: The role of surface chemistry and roughness}, journal = {Journal of Chemical Physics}, volume = {139}, year = {2013}, month = {Oct-14-2013}, type = {Article}, abstract = {

Experiments and simulations both suggest that the pressure experienced by an adsorbed phase confined within a carbon nanoporous material can be several orders of magnitude larger than the bulk phase pressure in equilibrium with the system. To investigate this pressure enhancement, we report a molecular-simulation study of the pressure tensor of argon confined in slit-shaped nanopores with walls of various models, including carbon and silica materials. We show that the pressure is strongly enhanced by confinement, arising from the effect of strongly attractive wall forces; confinement within purely repulsive walls does not lead to such enhanced pressures. Simulations with both the Lennard-Jones and Barker-Fisher-Watts intermolecular potentials for argon-argon interactions give rise to similar results. We also show that an increase in the wall roughness significantly decreases the in-pore pressure due to its influence on the structure of the adsorbate. Finally, we demonstrate that the pressures calculated from the mechanical (direct pressure tensor calculations) and the thermodynamic (volume perturbation method) routes yield almost identical results, suggesting that both methods can be used to calculate the local pressure tensor components in the case of these planar geometries. (C) 2013 AIP Publishing LLC.

}, issn = {0021-9606}, doi = {10.1063/1.4824125}, author = {Yun Long and Jeremy C. Palmer and Benoit A. Coasne and Malgorzata Sliwinska-Bartkowiak and George Jackson and Erich A. Mueller and Keith E. Gubbins} } @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 {30, title = {Nanofluidic Osmotic Diodes: Theory and Molecular Dynamics Simulations}, journal = {Physical Review Letters}, volume = {111}, year = {2013}, month = {Dec-2013}, pages = {244501}, abstract = {

Osmosis describes the flow of water across semipermeable membranes powered by the chemical free energy extracted from salinity gradients. While osmosis can be expressed in simple terms via the van \’t Hoff ideal gas formula for the osmotic pressure, it is a complex phenomenon taking its roots in the subtle interactions occurring at the scale of the membrane nanopores. Here we use new opportunities offered by nanofluidic systems to create an osmotic diode exhibiting asymmetric water flow under reversal of osmotic driving. We show that a surface charge asymmetry built on a nanochannel surface leads to nonlinear couplings between water flow and the ion dynamics, which are capable of water flow rectification. This phenomenon opens new opportunities for water purification and complex flow control in nanochannels.

Figure

Figure

Figure

Figure

}, doi = {10.1103/PhysRevLett.111.244501}, author = {Picallo, Clara B and Simon Gravelle and Joly, Laurent and Charlaix, Elisabeth and Lyd{\'e}ric Bocquet} } @conference {36, title = {Water Isotherms, Shrinkage and Creep of Cement Paste: Hypotheses, Models and Experiments}, booktitle = {Mechanics and Physics of Creep, Shrinkage, and Durability of Concrete: A Tribute to Zdenek P. Ba{\.z}ant}, year = {2013}, month = {Sep-25-2013}, publisher = {ASCE Publications}, organization = {ASCE Publications}, abstract = {

Cement paste has a complex mesoscale structure, and small changes in its pore network potentially cause large variation in measurements such as the water isotherm (also nitrogen). We deconvolute the water isotherm with the help of advanced computational techniques, hypotheses, and a re-examination of published data. The pore system is divided into four different categories, each containing water with its own physical properties. By viewing the highly interdependent roles of water in each of the pore categories as a system, new insights are gained regarding possible mechanisms that control drying, shrinkage, and creep, and experimental strategies for verification.

}, doi = {10.1061/9780784413111.015}, author = {Hamlin M. Jennings and Enrico Masoero and Pinson, Matthew B and Strekalova, Elena G and Patrick A Bonnaud and Hegoi Manzano and Qing Ji and Jeffrey J. Thomas and Roland Jean-Marc Pellenq and Franz-Josef Ulm and others} } @article {294, title = {Collapse Dynamics and Runout of Dense Granular Materials in a Fluid}, journal = {Physical Review Letters}, volume = {109}, year = {2012}, month = {Nov-02-2012}, pages = {Article Number: 188001}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.109.188001}, author = {Vincent Topin and Yann Monerie and Perales, F. and Farhang Radja{\"\i}} } @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 {52, title = {Enhanced CO2 solubility in hybrid adsorbents: optimization of solid support and solvent properties for CO2 capture}, journal = {The Journal of Physical Chemistry C}, volume = {116}, year = {2012}, month = {Feb-9-2012}, pages = {3600{\textendash}3607}, abstract = {
Abstract Image

In this study, we optimize the CO2 adsorption performance of hybrid adsorbents prepared by confining physical solvents in porous solid supports. A number of prospective solid supports and physical solvents are chosen to prepare hybrid adsorbents, and are subsequently evaluated in CO2 adsorption experiments. Generally, all the hybrid adsorbents show an enhancement of CO2 solubility compared to the bulk physical solvent. However, not all the adsorbents positively display an improvement in the CO2 adsorption performance as compared with the original solids after confining the physical solvent into the solids\’ pore. The micropore blocking effect is observed in the impregnated forms of zeolite, activated carbon, silicagel, and cecagel. Furthermore, we have obtained certain requisites for a good solid support, as efficient structures should be mesoporous with large surface area. In addition, there is an optimized solvent\’s size to achieve an optimized enhanced solubility. As a result, among the candidates, N-methyl-2-pyrrolidone confined in MCM-41 and alumina are identified as the most suitable hybrid adsorbents for an effective CO2-removal application.

}, doi = {10.1021/jp2099625}, author = {Linh Ngoc Ho and Perez-Pellitero, Javier and Porcheron, Fabien and Roland Jean-Marc Pellenq} } @article {54, title = {Evidence on the Dual Nature of Aluminum in the Calcium-Silicate-Hydrates Based on Atomistic Simulations}, journal = {Journal of the American Ceramic Society}, volume = {95}, year = {2012}, month = {Jan-23-2012}, pages = {1128{\textendash}1137}, abstract = {

Hydration of tri-calcium silicate (C3S) and di-calcium silicate (C2S) precipitates calcium-silicate-hydrate (CSH) which is the bonding phase responsible for the strength of cementitious materials. Substitution of part of C3S and C2S with aluminum-containing additives alters the chemical composition of hydration products by precipitating calcium-aluminate-silicate-hydrate (CASH). Incorporation of aluminum in the molecular building blocks of CSH entails structural and chemo-mechanical consequences. These alterations can be measured through solid state nuclear magnetic resonance (NMR) experiments. By conducting a wide spectrum of atomistic simulation methods on thousands of aluminum-containing molecular CASH structures, an overall molecular approach for determination of CASH nanostructure is presented. Through detailed analysis of different order parameters, it is found that aluminum can exhibit a tetra-/penta-/octahedral behavior which is fully consistent with the recent NMR observations. This corresponds to the formation of a class of complex three-dimensional alumino-silicate skeletons with partial healing effect in the CASH nanostructure potentially increasing durability and strength of hydration products. We explored the variation of mechanical observables by increasing aluminum content in CASH structures of varying calcium to silicon ratio. Finally, deformation of CSHs and CASHs of different chemical formula in a multi-scale fashion unravels the effect of chemical composition on the strength and kinematics of deformation in this particular type of composites.

}, doi = {10.1111/j.1551-2916.2011.05058.x}, author = {Mohammad Javad Abdolhosseini Qomi 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} } @proceedings {354, title = {Microscopic Origins of Shear Strength in Packings Composed of Elongated Particles}, journal = {9th International Workshop on Bifurcation and Degradation in Geomaterials (IWBDG 2011)}, volume = {Springer Series in Geomechanics and Geoengineering - ADVANCES IN BIFURCATION AND DEGRADATION IN GEOMATERIALS}, year = {2012}, month = {Apr-11-2012}, pages = {21-27}, address = {MAY 23-26 2011 Porquerolles, FRANCE}, abstract = {

We investigate the rheology, force transmission and texture of granular materials composed of elongated particles by means of contact dynamics simulations. The particles have a rounded-cap rectangular (RCR) shape described by a single elongation parameter varying from 0 for a circular particle to 1 for an infinitely thin or long particle. We study the quasi-static behavior, structural and force anisotropies as a function of the elongation parameter for packings submitted to biaxial compression. The shear strength is found to increase linearly with this parameter whereas the solid fraction both at the initial isotropic state and in the critical state is nonmonotonous. We show that for these elongated particles a harmonic decomposition of the stress tensor provides a fairly good approximation of the internal state. Our data suggest that the increase of shear strength with reflects both enhanced friction mobilization and anisotropic particle orientation as the elongation of the particles increases.

}, keywords = {Elongated particles, Fabric properties, force transmission, Harmonic decomposition}, url = {https://hal.archives-ouvertes.fr/hal-00686728}, author = {Emilien Az{\'e}ma and Farhang Radja{\"\i}}, editor = {St{\'e}phane Bonelli and Dascalu, Cristian and Fran{\c c}ois Nicot} } @article {47, title = {Nanoscale Elastic Properties of Montmorillonite upon Water Adsorption}, journal = {Langmuir}, volume = {28}, year = {2012}, month = {Dec-11-2012}, pages = {16855{\textendash}16863}, abstract = {
Abstract Image

Smectites are an important group of clay minerals that experience swelling upon water adsorption. This paper uses molecular dynamics with the CLAYFF force field to simulate isothermal isobaric water adsorption of interlayer Wyoming Na-montmorillonite, a member of the smectite group. Nanoscale elastic properties of the clay\–interlayer water system are calculated from the potential energy of the model system. The transverse isotropic symmetry of the elastic constant matrix was assessed by calculating Euclidean and Riemannian distance metrics. Simulated elastic constants of the clay mineral are compared with available results from acoustic and nanoindentation measurements.

}, doi = {10.1021/la302997g}, author = {Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Andrew J. Whittle} } @article {48, title = {Nanostructure and nanomechanics of cement: polydisperse colloidal packing}, journal = {Physical review letters}, volume = {109}, year = {2012}, month = {Oct-10-2012}, pages = { Article Number: 155503}, abstract = {

Cement setting and cohesion are governed by the precipitation and growth of calcium-silicate-hydrate, through a complex evolution of microstructure. A colloidal model to describe nucleation, packing, and rigidity of calcium-silicate-hydrate aggregates is proposed. Polydispersity and particle size dependent cohesion strength combine to produce a spectrum of packing fractions and of corresponding elastic properties that can be tested against nanoindentation experiments. Implications regarding plastic deformations and reconciling current structural characterizations are discussed.

Figure

Figure

Figure

}, doi = {10.1103/PhysRevLett.109.155503}, author = {Enrico Masoero and Emanuela Del Gado and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Sidney Yip} } @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 {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 {49, title = {Solvated calcium ions in charged silica nanopores}, journal = {The Journal of Chemical Physics}, volume = {137}, year = {2012}, month = {Aug-14-2012}, pages = {Article Number: 064706}, abstract = {

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

Image of FIG. 1.

Image of FIG. 2.

Image of FIG. 3.

Image of FIG. 4.

Image of FIG. 5.

Image of FIG. 6.

}, doi = {10.1063/1.4742854}, author = {Patrick A Bonnaud and Benoit A. Coasne and Roland Jean-Marc Pellenq} } @article {58, title = {Structure and properties of nanoscale materials: theory and atomistic computer simulation}, journal = {International Journal of Nanotechnology}, volume = {9}, year = {2012}, month = {2012}, pages = {576{\textendash}604}, abstract = {

We present a review of a few research topics developed within the \"Theory and Atomistic Computer Simulation\" Department at CINaM. The bottom line of the scientific activity is to use up\–to\–date theoretical and computer simulation techniques to address physics and materials science problems, often at the nanometric scale, in close contact with experimental groups. It ranges from the study of the structure and properties of molecular systems for organic electronics to metallic clusters and alloys, magnetic oxides, nuclear fuels and carbon\–based nanostructures. These studies are motivated by fundamental research questions as well as more applied goals including environmental and energy issues, or information technologies. This broad spectrum of activities requires a large range of techniques, from theory and ab initio calculations to semi\–empirical models incorporated in Monte Carlo or molecular dynamics simulations.

}, doi = {10.1504/IJNT.2012.045335 }, author = {Christophe Bichara and Marsal, P and Mottet, C and Roland Jean-Marc Pellenq and F. Ribeiro and Andres Sa{\`u}l and G. Tr{\'e}glia and Weissker, H{\textendash}Ch} } @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 {55, title = {Understanding and Controlling the Reactivity of the Calcium Silicate phases from First Principles}, journal = {Chemistry of Materials}, volume = {24}, year = {2012}, month = {Apr-10-2012}, pages = {1262{\textendash}1267}, abstract = {
Abstract Image

First principles calculations are employed to provide a fundamental understanding of the relationship between the reactivity of synthetic calcium silicate phases and their electronic structure. Our aim is to shed light on the wide range of hydration kinetics observed in different phases of calcium silicate. For example, while the dicalcium silicate (Ca2SiO4) phase slowly reacts with water, the tricalcium silicate (Ca3SiO5) shows much faster hydration kinetics. We show that the high reactivity of Ca3SiO5 is mainly related to the reactive sites around its more ionic oxygen atoms. Ca2SiO4 does not contain these types of oxygen atoms, although experiments suggest that impurities may play a role in changing the reactivity of these materials. We analyze the electronic structure of a wide range of possible substitutions in both Ca3SiO5 and Ca2SiO4 and show that while the influence of different types of impurities on structural properties is similar, their effect on reactivity is very different. Our calculations suggest that the variation of electronic structure is mainly related to the formation of new hybridized orbitals and the charge exchange between the impurity atoms and the bulk material. The charge localization upon introducing impurities is quantified to predict candidate substitutions that could increase the reactivity of Ca2SiO4, which would broaden the applicability of this lower temperature and thus less costly and energetically less demanding phase.

}, doi = {10.1021/cm203127m}, author = {Engin Durgun and Hegoi Manzano and Roland Jean-Marc Pellenq and Jeffrey C Grossman} } @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 {289, title = {An atomistic modelling of the porosity impact on UO2 matrix macroscopic properties}, journal = {Journal of Nuclear Materials}, volume = {415}, year = {2011}, month = {Aug-15-2011}, pages = {210 - 216}, abstract = {

The porosity impact on the UO2 matrix thermomechanical properties was investigated using atomistic simulation techniques. The porosity modifies the thermal expansion coefficient and this is attributed to pore surface effects. The elastic moduli at 0\ K and at finite temperature decrease with porosity, this variation being well approximated using affine functions. These results agree with other mesoscale model predictions and experimental data, showing the ability of the semiempirical potential atomistic simulations to give an overall good description of the porous UO2. However, the surface effects are incompletely described.

Fig. 1. The UO2 elementary cell contains 12 atoms: 4 uranium atoms (large spheres){\textellipsis}

fig.=

Fig. 3. Example of 6{\texttimes}6{\texttimes}6 one pore relaxed supercells of different pore radius (4, 6 and{\textellipsis}

Fig. 4. One layer of UO2 as considered in the present work: top {\textendash} the (100){\textellipsis}

Fig. 5. Variation of the UO2 thermal expansion coefficient versus temperature: top {\textendash}{\textellipsis}

Fig. 6. Variation of the shear modulus (G) and Young{\textquoteright}s modulus (E) versus porosity at 0K

}, issn = {00223115}, doi = {10.1016/j.jnucmat.2011.06.006}, author = {Andrei Jelea and Colbert, M. and F. Ribeiro and G. Tr{\'e}glia and Roland Jean-Marc Pellenq} } @proceedings {356, title = {Compressive strength of an unsaturated granular material during cementation}, journal = {Symposium on Science and Technology of Powders and Sintered Materials (STPMF 2009)}, volume = {POWDER TECHNOLOGY}, year = {2011}, month = {Mar-25-2011}, pages = {308 - 311}, address = {MAY 25-27 2009 Montpellier, FRANCE}, abstract = {

The cohesive behaviour of unsaturated granular materials is due to the presence of cohesive bonds between grains. These bonds can have various physico-chemical characteristics and may evolve with environmental conditions. We study the case of a granular material partially saturated by an aqueous solution. The bonds are thus initially of capillary type and the mechanical strength is weak. At low relative humidity, the phase change of water involves crystallization of the solute at the contact points between grains, generating thus solid bonds. The mechanical strength of the material is then enhanced. An experimental study of the evolution of the mechanical strength during crystallization of the solute shows clearly the transition from capillary regime to cemented regime. This transition is not correlated with the mass of the crystallized solute, but rather with the residual degree of saturation. This behavior is analyzed here in the light of discrete element simulations. We introduce a local cohesion law that accounts for transition from capillary to cemented bonding. This law is formulated in terms of the degree of crystallization as a result of the evaporation of water at the boundary of the sample. The cohesion of the packing is initially of capillary type. A crystallization front then spreads from the sample boundaries to the center of the sample, and the strength increases as a result. Uniaxial compression allows us to determine the strength at different times. The numerical strength agrees well with the experimental data, and reveals strength enhancement as the solute crystallizes, as well as the transition from capillary to cementation regime.

}, issn = {00325910}, doi = {10.1016/j.powtec.2010.08.021}, author = {Jean-Yves Delenne and Souli{\'e}, Fabien and Moulay Sa{\"\i}d El Youssoufi and Farhang Radja{\"\i}} } @article {290, title = {Enhanced CO2 Solubility in Hybrid MCM-41: Molecular Simulations and Experiments}, journal = {Langmuir}, volume = {27}, year = {2011}, month = {Jul-05-2011}, pages = {8187 - 8197}, abstract = {
Abstract Image

Grand canonical Monte Carlo simulations are performed in a hybrid adsorbent model in order to interpret the CO2 solubility behavior. The hybrid adsorbent is prepared by confining a physical solvent (OMCTS) into the pores of a mimetic MCM-41 solid support. As a result, simulated adsorption isotherms of CO2 nicely match the experimental data for three distinctive systems: bulk solvent, raw MCM-41, and hybrid MCM-41. The microscopic mechanisms underlying the apparition of enhanced solubility are then clearly identified. In fact, the presence of solvent molecules favors the layering of CO2 molecules within the pores; therefore, the CO2 solubility in the hybrid adsorbent markedly increases in comparison to that found in the raw adsorbent as well as in the bulk solvent. In addition, a good understanding of confined solvents\’ properties and solid surface structures is essential to fully evaluate the efficiency of hybrid adsorbents in capturing CO2. The sorbent\–solid interactions along with the solvent molecular size\’s impact on CO2 solubility are therefore investigated in this study. We found that an ideal hybrid system should possess a weak solvent\–solid interaction but a strong solvent\–CO2 interaction. Besides, an optimal solvent size is obtained for the enhanced CO2 solubility in the hybrid system. According to the simulation results, the solvent layer builds pseudomicropores inside the mesoporous MCM-41, enabling more CO2 molecules to be absorbed under the greater influence of spatial confinement and surface interaction. In addition, the molecular sieving effect is clearly observed in the case of larger solvent molecular sizes.

}, issn = {0743-7463}, doi = {10.1021/la2012765}, author = {Linh Ngoc Ho and Perez Pellitero, Javier and Porcheron, Fabien and Roland Jean-Marc Pellenq} } @article {292, title = {Glassy Nature of Water in an Ultraconfining Disordered Material: The Case of Calcium-Silicate-Hydrate}, journal = {Journal of the American Chemical Society}, volume = {133}, year = {2011}, month = {Mar-02-2011}, pages = {2499 - 2510}, abstract = {
Abstract Image

We present the structural and dynamic nature of water ultraconfined in the quasi-two-dimensional nanopores of the highly disordered calcium\−silicate\−hydrate (C-S-H), the major binding phase in cement. Our approach is based on classical molecular simulations. We demonstrate that the C-S-H nanopore space is hydrophilic, particularly because of the nonbridging oxygen atoms on the disordered silicate chains which serve as hydrogen-bond acceptor sites, directionally orienting the hydrogen atoms of the interfacial water molecules toward the calcium\−silicate layers. The water in this interlayer space adopts a unique multirange structure: a distorted tetrahedral coordination at short range up to 2.7 \Å, a disordered structure similar to that of dense fluids and supercooled phases at intermediate range up to 4.2 \Å, and persisting spatial correlations through dipole\−dipole interactions up to 10 \Å. A three-stage dynamics governs the mean square displacement (MSD) of water molecules, with a clear cage stage characteristic of the dynamics in supercooled liquids and glasses, consistent with its intermediate-range structure identified here. At the intermediate time scales corresponding to the \β-relaxation of glassy materials, coincident with the cage stage in MSD, the non-Gaussian parameter indicates a significant heterogeneity in the translational dynamics. This dynamic heterogeneity is induced primarily because of the heterogeneity in the distribution of hydrogen bond strengths. The strongly attractive interactions of water molecules with the calcium silicate walls serve to constrain their motion. Our findings have important implications on describing the cohesion and mechanical behavior of cement from its setting to its aging.

}, issn = {0002-7863}, doi = {10.1021/ja107003a}, author = {Youssef, Mostafa and Roland Jean-Marc Pellenq and Bilge Yildiz} } @article {291, title = {Impact of Chemical Impurities on the Crystalline Cement Clinker Phases Determined by Atomistic Simulations}, journal = {Crystal Growth \& Design}, volume = {11}, year = {2011}, month = {Jul-2011}, pages = {2964 - 2972}, abstract = {
Abstract Image

The presence of chemical substitutions is believed to play a crucial role in the hydration reactions, structure, and elastic properties of cement clinker phases. Hence, substitutions are of great technological interest, as more efficient production of cement clinkers would result in a reduction of CO2 emissions, as well as possible economic benefits. Here we use a combination of classical and quantum mechanical simulation methods to study the detailed physicochemical changes of the clinker phases alite (Ca3SiO5) and belite (Ca2SiO4) when Mg2+, Al3+ and Fe3+ guest ions are incorporated into their structure. Using classical force field methods, we considered random substitutions among possible sites and different compositions in order to identify the preferential substitution sites on the crystalline structures. Then, the resulting structural changes that take place to accommodate the guest ions are investigated and discussed in detail. Using quantum mechanical density functional theory calculations the electronic structure of representative configurations has been computed to determine the potential impact of impurities on the reactivity.

}, issn = {1528-7483}, doi = {10.1021/cg200212c}, author = {Hegoi Manzano and Engin Durgun and Mohammad Javad Abdolhosseini Qomi and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Jeffrey C Grossman} } @proceedings {352, title = {Modeling Porous Granular Aggregates}, journal = {9th International Workshop on Buifurcation and Degradation in Geomaterials (IWBDG 2011)}, volume = {Springer Series in Geomechanics and Geoengineering - ADVANCES IN BIFURCATION AND DEGRADATION IN GEOMATERIALS}, year = {2011}, month = {May-28-2011}, pages = {249 - 254}, publisher = {Springer Netherlands}, address = {MAY 23-26 2011 Porquerolles, FRANCE}, abstract = {

We rely on 3D simulations based on the Lattice Element Method (LEM) to analyze the failure of porous granular aggregates under tensile loading. We investigate crack growth by considering the number of broken bonds in the particle phase as a function of the matrix volume fraction and particle-matrix adhesion. Three regimes are evidenced, corresponding to no particle damage, particle abrasion and particle fragmentation, respectively. We also show that the probability density of strong stresses falls off exponentially at high particle volume fractions where a percolating network of jammed particles occurs. Decreasing the matrix volume fraction leads to increasingly broader stress distribution and hence a higher stress concentration. Our findings are in agreement with 2D results previously reported in the literature.

}, isbn = {978-94-007-1420-5}, issn = {1866-8755}, doi = {10.1007/978-94-007-1421-2_32}, author = {Rafik Aff{\`e}s and Vincent Topin and Jean-Yves Delenne and Yann Monerie}, editor = {St{\'e}phane Bonelli and Dascalu, Cristian and Fran{\c c}ois Nicot} } @proceedings {353, title = {Onset of Immersed Granular Avalanches by DEM-LBM Approach}, journal = {9th International Workshop on Buifurcation and Degradation in Geomaterials (IWBDG 2011)}, volume = {Springer Series in Geomechanics and Geoengineering - ADVANCES IN BIFURCATION AND DEGRADATION IN GEOMATERIALS}, year = {2011}, month = {May-28-2011}, pages = {109 - 115}, publisher = {Springer Netherlands}, address = {MAY 23-26 2011 Porquerolles, FRANCE}, abstract = {

We present 3D grain-fluid simulations based on the discrete element method interfaced with the lattice Boltzmann method and applied to investigate the initiation of underwater granular flows. We prepare granular beds of 800 spherical grains with different values of the initial solid fraction in a biperiodic rectangular box. In order to trigger an avalanche, the bed is instantaneously tilted to a finite slope angle above the maximum angle of stability. We simulate the dynamics of the transient flow for different solid fractions. In agreement with the experimental work of Iverson (Water Resour Res 36(7):1897\–1910, 2000) and Pailha et\ al. (Phys Fluids 20(11):111701, 2008), we find that the flow onset is controlled by the initial solid fraction.

}, isbn = {978-94-007-1420-5}, issn = {1866-8755}, doi = {10.1007/978-94-007-1421-210.1007/978-94-007-1421-2_14}, author = {Jean-Yves Delenne and Mansouri, M. and Farhang Radja{\"\i} and Moulay Sa{\"\i}d El Youssoufi and Seridi, A.}, editor = {St{\'e}phane Bonelli and Dascalu, Cristian and Fran{\c c}ois Nicot} } @article {287, title = {Reply to {\textquotedblleft}Comment on {\textquoteleft}Determination of the bulk melting temperature of nickel using Monte Carlo simulations: Inaccuracy of extrapolation from cluster melting temperatures{\textquoteright} {\textquotedblright}}, journal = {Physical Review B}, volume = {84}, year = {2011}, month = {Sep-26-2011}, pages = { Article Number: 096103}, abstract = {

The thermodynamic integration (TI) method for determining the melting temperature has to be performed with great care to ensure accuracy. Since the slopes of the Gibbs free energy curves as a function of temperature (or pressure) are usually rather close, small shifts in the free energy curves can lead to a large change in the calculated melting temperature. In this Reply to a Comment by Harvey and Gheribi [ Phys. Rev. B 84 096102 (2011)] on our paper [ Phys. Rev. B 81 064112 (2010)], we show that the error in the calculated melting temperature of nickel induced by some technical imperfections in our application of the TI method, as noticed and described in the Comment, is small compared to the error bar given in our paper. We also clarify supposed doubts about the accuracy of the TI method based on Johnson\&$\#$39;s free energy expression and parametrization [ Johnson et al. Mol. Phys. 78 591 (1993)] of the Lennard-Jones fluid.

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.84.096103}, author = {Los, J. H. and Roland Jean-Marc Pellenq} } @proceedings {364, title = {Springer Series in Geomechanics and Geoengineering Multiscale and Multiphysics Processes in Geomechanics SIMULATION OF HYDRATION AND ELASTIC PROPERTIES OF MONTMORILLONITE USING MOLECULAR DYNAMICS}, journal = {International Workshop on Multiscale and Multiphysics Processes in Geomechanics}, volume = {Springer Series in Geomechanics and Geoengineering}, year = {2011}, month = {2011}, pages = {105 - 108}, publisher = {Springer Berlin Heidelberg}, address = {JUN 23-25 2010 Stanford Univ, Stanford, CA}, abstract = {

This article describes molecular-level simulations of the adsorption isotherm and elastic properties for a Wyoming Na-montmorillonite using the General Utility Lattice Program. The selected clay mineral includes isomorphous substitutions in both the silica tetrahedral and octahedral gibbsite sheets. The resulting negative charges are balanced by Na+ cations. The authors have investigated the importance of different interatomic potentials in understanding the adsorption of water molecules (at ambient temperature, 300K). A unique combination of NPT and Grand Canonical ensemble simulations using Core\–Shell potential shows intra-layer adsorption, while the same approach based on empirical CLAYFF potential produces only interlayer water adsorption. There is a correspondingly large difference in the theoretical elastic modulus normal to the clay sheets using the two different sets of force potentials. This article discusses the underlying reasons for these differences and reviews the experimental evidence supporting intra-layer water absorption in montmorillonite.

}, isbn = {978-3-642-19629-4}, issn = {1866-8755}, doi = {10.1007/978-3-642-19630-0_27}, author = {Roland Jean-Marc Pellenq and Andrew J. Whittle and Ebrahimi, Davoud}, editor = {Borja, Ronaldo I. and Wu, Wei} } @proceedings {355, title = {Stress fields in granular solids: Effect of composition}, journal = {Symposium on Science and Technology of Powders and Sintered Materials (STPMF 2009)}, volume = {POWDER TECHNOLOGY}, year = {2011}, month = {Mar-25-2011}, pages = {568 - 573}, address = {MAY 25-27 2009 Montpellier, FRANCE}, abstract = {

We use the lattice element method to investigate stress fields at the sub-particle scale in granular solids composed of particles embedded in a cementing matrix. The stress distributions are found to be similar in 2D and 3D samples subjected to vertical loading with free lateral boundaries. We find that the number of strong forces falls off exponentially at high particle volume fractions where a percolating network of jammed particles occurs. The influence of the matrix volume fraction and particle/matrix stiffness ratio with respect to stress distribution is analyzed in 2D and 3D. We show that both decreasing the matrix volume fraction and increasing the stiffness ratio lead to increasingly broader distributions within a limit beyond which the distribution is independent of one or both of these parameters.

}, issn = {00325910}, doi = {10.1016/j.powtec.2010.08.060}, author = {Vincent Topin and Farhang Radja{\"\i} and Jean-Yves Delenne} } @article {288, title = {Tight binding within the fourth moment approximation: Efficient implementation and application to liquid Ni droplet diffusion on graphene}, journal = {Physical Review B}, volume = {84}, year = {2011}, month = {Aug-31-2011}, pages = { Article Number: 085455}, abstract = {

Application of the fourth moment approximation (FMA) to the local density of states within a tight binding description to build a reactive, interatomic interaction potential for use in large scale molecular simulations, is a logical and significant step forward to improve the second moment approximation, standing at the basis of several, widely used (semi-)empirical interatomic interaction models. In this paper we present a sufficiently detailed description of the FMA and its technical implications, containing the essential elements for an efficient implementation in a simulation code. Using a recent, existing FMA-based model for C-Ni systems, we investigated the size dependence of the diffusion of a liquid Ni cluster on a graphene sheet and find a power law dependence of the diffusion constant on the cluster size (number of cluster atoms) with an exponent very close to \−2/3, equal to a previously found exponent for the relatively fast diffusion of solid clusters on a substrate with incommensurate lattice matching. The cluster diffusion exponent gives rise to a specific contribution to the cluster growth law, which is due to cluster coalescence. This is confirmed by a simulation for Ni cluster growth on graphene, which shows that cluster coalescence dominates the initial stage of growth, overruling Oswald ripening.

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.84.085455}, author = {Los, J. H. and Christophe Bichara and Roland Jean-Marc Pellenq} } @article {293, title = {Empirical force fields for complex hydrated calcio-silicate layered materials}, journal = {Phys. Chem. Chem. Phys.}, volume = {13}, year = {2010}, month = {Nov-10-2010}, pages = {1002 - 1011}, abstract = {

The use of empirical force fields is now a standard approach in predicting the properties of hydrated oxides which are omnipresent in both natural and engineering applications. Transferability of force fields to analogous hydrated oxides without rigorous investigations may result in misleading property predictions. Herein, we focus on two common empirical force fields, the simple point charge ClayFF potential and the core-shell potential to study tobermorite minerals, the most prominent family of Calcium-Silicate-Hydrates that are complex hydrated oxides. We benchmark the predictive capabilities of these force fields against first principles results. While the structural information seem to be in close agreement with DFT results, we find that for higher order properties such as elastic constants, the core-shell potential quantitatively improves upon the simple point charge model, and shows a larger degree of transferability to complex materials. In return, to remedy the deficiencies of the simple point charge potential for hydrated calcio-silicates, we suggest using both structural data and elasticity data for potential calibration, a new force field potential, CSH-FF. This re-parameterized version of ClayFF is then applied to simulating an atomistic model of cement (Pellenq et al., PNAS, 2009). We demonstrate that this force field improves the predictive capabilities of ClayFF, being considerably less computational intensive than the core-shell model.

}, issn = {1463-9076}, doi = {10.1039/C0CP00516A}, author = {Rouzbeh Shahsavari and Roland Jean-Marc Pellenq and Franz-Josef Ulm} }