Carbon-cement supercapacitors as a scalable bulk energy storage solution.

The large-scale implementation of renewable energy systems necessitates the development of energy storage solutions to effectively manage imbalances between energy supply and demand. Herein, we investigate such a scalable material solution for energy storage in supercapacitors constructed from readily available material precursors that can be locally sourced from virtually anywhere on the planet, namely cement, water, and carbon black. We characterize our carbon-cement electrodes by combining correlative EDS-Raman spectroscopy with capacitance measurements derived from cyclic voltammetry and galvanostatic charge-discharge experiments using integer and fractional derivatives to correct for rate and current intensity effects. Texture analysis reveals that the hydration reactions of cement in the presence of carbon generate a fractal-like electron-conducting carbon network that permeates the load-bearing cement-based matrix. The energy storage capacity of this space-filling carbon black network of the high specific surface area accessible to charge storage is shown to be an intensive quantity, whereas the high-rate capability of the carbon-cement electrodes exhibits self-similarity due to the hydration porosity available for charge transport. This intensive and self-similar nature of energy storage and rate capability represents an opportunity for mass scaling from electrode to structural scales. The availability, versatility, and scalability of these carbon-cement supercapacitors opens a horizon for the design of multifunctional structures that leverage high energy storage capacity, high-rate charge/discharge capabilities, and structural strength for sustainable residential and industrial applications ranging from energy autarkic shelters and self-charging roads for electric vehicles, to intermittent energy storage for wind turbines and tidal power stations.

The role of concrete in life cycle greenhouse gas emissions of US buildings and pavements.

Concrete is a critical component of deep decarbonization efforts because of both the scale of the industry and because of how its use impacts the building, transportation, and industrial sectors. We use a bottom-up model of current and future building and pavement stocks and construction in the United States to contextualize the role of concrete in greenhouse gas (GHG) reductions strategies under projected and ambitious scenarios, including embodied and use phases of the structures' life cycle. We show that projected improvements in the building sector result in a reduction of 49% of GHG emissions in 2050 relative to 2016 levels, whereas ambitious improvements result in a 57% reduction in 2050, which is 22.5 Gt cumulative saving. The pavements sector shows a larger difference between the two scenarios with a 14% reduction of GHG emissions for projected improvements and a 65% reduction under the ambitious scenario, which is ∼1.35 Gt. This reduction occurs despite the fact that concrete usage in 2050 in the ambitious scenario is over three times that of the projected scenario because of the ways in which concrete lowers use phase emissions. Over 70% of future emissions from new construction are from the use phase.

Multiscale poromechanics of wet cement paste.

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.

Time-Space-Resolved Chemical Deconvolution of Cementitious Colloidal Systems Using Raman Spectroscopy.

Concrete is one of the most used materials in the world, second only to water. One of the key advantages of this versatile material is its workability in the early stages before setting. Here, we use in situ underwater Raman microspectroscopy to investigate and visualize the early hydration kinetics of ordinary Portland cement (OPC) with submicron spatial and high temporal resolution. First, the spectral features of the C-S-H gel were analyzed in the hydroxyl stretching region to confirm the coexistence of Ca-OH and Si-OH bonds in a highly disordered C-S-H gel. Second, the disordered calcium hydroxide (Ca(OH)2) is experimentally identified for the first time in the mixture before setting, suggesting that Ca(OH)2 crystallization and growth are essential in the setting of cement paste. Finally, the phase transformations of clinker, C-S-H, and Ca(OH)2 are spatially and temporally resolved, and the hydration kinetics are studied by analyzing the spatial relationships of these phases using two-point correlation functions. The results quantitatively validate that the setting occurs as a percolation process, wherein the hydration products intersect and form an interconnected network. This time-space-resolved characterization method can map and quantitatively analyze the heterogeneous reaction of the cementitious colloidal system and thus provide potential application value in the field of cement chemistry and materials design more broadly.

Mesoscale structure, mechanics, and transport properties of source rocks' organic pore networks.

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

Effect of Confinement on Capillary Phase Transition in Granular Aggregates.

Using a 3D mean-field lattice-gas model, we analyze the effect of confinement on the nature of capillary phase transition in granular aggregates with varying disorder and their inverse porous structures obtained by interchanging particles and pores. Surprisingly, the confinement effects are found to be much less pronounced in granular aggregates as opposed to porous structures. We show that this discrepancy can be understood in terms of the surface-surface correlation length with a connected path through the fluid domain, suggesting that this length captures the true degree of confinement. We also find that the liquid-gas phase transition in these porous materials is of second order nature near capillary critical temperature, which is shown to represent a true critical temperature, i.e., independent of the degree of disorder and the nature of the solid matrix, discrete or continuous. The critical exponents estimated here from finite-size scaling analysis suggest that this transition belongs to the 3D random field Ising model universality class as hypothesized by F. Brochard and P.G. de Gennes, with the underlying random fields induced by local disorder in fluid-solid interactions.

Mesoscale texture of cement hydrates.

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.

Impact of Nanoporosity on Hydrocarbon Transport in Shales' Organic Matter.

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

Atomistic and mesoscale simulation of sodium and potassium adsorption in cement paste.

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.

Realistic molecular model of kerogen's nanostructure.

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.

Stress Transmission and Failure in Disordered Porous Media.

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.

Rigidity transition in materials: hardness is driven by weak atomic constraints.

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.

Optimized molecular reconstruction procedure combining hybrid reverse Monte Carlo and molecular dynamics.

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.

Nano-scale mechanics of colloidal C-S-H gels.

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.

Anomalous composition-dependent dynamics of nanoconfined water in the interlayer of disordered calcium-silicates.

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.

Molecular simulations of supercritical fluid permeation through disordered microporous carbons.

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.

Confined water dissociation in microporous defective silicates: mechanism, dipole distribution, and impact on substrate properties.

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.

Hydration of calcium oxide surface predicted by reactive force field molecular dynamics.

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 CO(2) 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.

Nanogranular origin of concrete creep.

Concrete, the solid that forms at room temperature from mixing Portland cement with water, sand, and aggregates, suffers from time-dependent deformation under load. This creep occurs at a rate that deteriorates the durability and truncates the lifespan of concrete structures. However, despite decades of research, the origin of concrete creep remains unknown. Here, we measure the in situ creep behavior of calcium-silicate-hydrates (C-S-H), the nano-meter sized particles that form the fundamental building block of Portland cement concrete. We show that C-S-H exhibits a logarithmic creep that depends only on the packing of 3 structurally distinct but compositionally similar C-S-H forms: low density, high density, ultra-high density. We demonstrate that the creep rate ( approximately 1/t) is likely due to the rearrangement of nanoscale particles around limit packing densities following the free-volume dynamics theory of granular physics. These findings could lead to a new basis for nanoengineering concrete materials and structures with minimal creep rates monitored by packing density distributions of nanoscale particles, and predicted by nanoscale creep measurements in some minute time, which are as exact as macroscopic creep tests carried out over years.

A realistic molecular model of cement hydrates.

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.