Despite decades of studies of calcium-silicate-hydrate (C-S-H), the structurally complex binder phase of concrete, the interplay between chemical composition and density remains essentially unexplored. Together these characteristics of C-S-H define and modulate the physical and mechanical properties of this ''liquid stone'' gel phase. With the recent determination of the calcium/silicon (C/S ‫؍‬ 1.7) ratio and the density of the C-S-H particle (2.6 g/cm 3 ) by neutron scattering measurements, there is new urgency to the challenge of explaining these essential properties. Here we propose a molecular model of C-S-H based on a bottom-up atomistic simulation approach that considers only the chemical specificity of the system as the overriding constraint. By allowing for short silica chains distributed as monomers, dimers, and pentamers, this C-S-H archetype of a molecular description of interacting CaO, SiO2, and H2O units provides not only realistic values of the C/S ratio and the density computed by grand canonical Monte Carlo simulation of water adsorption at 300 K. The model, with a chemical composition of (CaO)1.65(SiO2)(H2O)1.75, also predicts other essential structural features and fundamental physical properties amenable to experimental validation, which suggest that the C-S-H gel structure includes both glass-like short-range order and crystalline features of the mineral tobermorite. Additionally, we probe the mechanical stiffness, strength, and hydrolytic shear response of our molecular model, as compared to experimentally measured properties of C-S-H. The latter results illustrate the prospect of treating cement on equal footing with metals and ceramics in the current application of mechanism-based models and multiscale simulations to study inelastic deformation and cracking.atomistic simulation ͉ mechanical properties ͉ structural properties B y mixing water and cement, a complex hydrated oxide called calcium-silicate-hydrate (C-S-H) precipitates as nanoscale clusters of particles (1). Much of our knowledge of C-S-H has been obtained from structural comparisons with crystalline calcium silicate hydrates, based on HFW Taylor's postulate that real C-S-H was a structurally imperfect layered hybrid of two natural mineral analogs (2) (4)]. While this suggestion is plausible in morphological terms, this model is incompatible with two basic characteristics of real C-S-H; specifically the calcium-tosilicon ratio (C/S) and the density. Recently, small-angle neutron scattering measurements have fixed the C/S ratio at 1.7 and the density at 2.6 g/cm 3 (1), values that clearly cannot be obtained from either tobermorite (C/S ϭ 0.83, 2.18 g/cm 3 ) or jennite (C/S ϭ 1.5 and 2.27 g/cm 3 ). From the standpoint of constructing a molecular model of C-S-H, this means that these crystalline minerals are not strict structural analogs. Here we adopt the perspective that the chemical composition of C-S-H is the most essential property in formulating a realistic molecular description. We show that once the C/S ratio is described...

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 sp(2)/sp(3) hybridization ratio, entails a crossover from plastic-to-brittle rupture mechanisms.

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

We present a constrained reverse Monte Carlo method for structural modeling of porous carbons. As in the original reverse Monte Carlo method, the procedure is to stochastically change the atomic positions of a system of carbon atoms to minimize the differences between the simulated and the experimental pair correlation functions. However, applying the original reverse Monte Carlo method without further constraints yields nonunique structures for carbons, due to the presence of strong three-body forces. In this respect, the uniqueness theorem of statistical mechanics provides a helpful guide to the design of reverse Monte Carlo methods that give reliable structures. In our method, we constrain the bond angle distribution and the average carbon coordination number to describe the three-body correlations. Using this procedure, we have constructed structural models of two highly disordered porous carbons prepared by pyrolysis of saccharose at two different temperatures. The resulting pair correlation functions are in excellent agreement with those obtained by diffraction experiments. Simulated transmission electron microscopy (TEM) images of the resulting models are compared to experimental images. Many of the features observed in the experimental images are also observed in the simulations. The model carbons are characterized by determination of the porosity, pore size distribution, adsorbent−adsorbate potential energy distribution, and adsorption properties at zero coverage, using a model of nitrogen as the adsorbate. Grand canonical Monte Carlo simulations of nitrogen adsorption in the model materials are presented, and it is found that the results can be explained in terms of the adsorbent−adsorbate potential energy distribution but not in terms of the pore size distribution. For both models, the isosteric heat of adsorption is a decreasing function of coverage, in agreement with typical experimental results in other porous carbons.