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.

This paper reports on a molecular simulation study of the thermodynamics, structure and dynamics of water confined at ambient temperature in hydroxylated silica nanopores of a width H = 10 and 20 Å. The adsorption isotherms for water in these nanopores resemble those observed for experimental samples; the adsorbed amount increases continuously in the multilayer adsorption regime until a jump occurs due to capillary condensation of the fluid within the pore. Strong layering of water in the vicinity of the silica surfaces is observed as marked density oscillations are observed up to 8 Å from the surface in the density profiles for confined water. Our results indicate that water molecules within the first adsorbed layer tend to adopt a H-down orientation with respect to the silica substrate. For all pore sizes and adsorbed amounts, the self-diffusivity of confined water is lower than the bulk, due to the hydrophilic interaction between the water molecules and the hydroxylated silica surface. Our results also suggest that the self-diffusivity of confined water is sensitive to the adsorbed amount.

Cement paste has a complex distribution of pores and molecular-scale spaces. This distribution controls the hysteresis of water sorption isotherms and associated bulk dimensional changes (shrinkage). We focus on two locations of evaporable water within the fine structure of pastes, each having unique properties, and we present applied physics models that capture the hysteresis by dividing drying and rewetting into two related regimes based on relative humidity (RH). We show that a continuum model, incorporating a poreblocking mechanism for desorption and equilibrium thermodynamics for adsorption, explains well the sorption hysteresis for a paste that remains above approximately 20% RH. In addition, we show with molecular models and experiments that water in spaces of ≲1 nm width evaporates below approximately 20% RH but reenters throughout the entire RH range. This water is responsible for a drying shrinkage hysteresis similar to that of clays but opposite in direction to typical mesoporous glass. Combining the models of these two regimes allows the entire drying and rewetting hysteresis to be reproduced accurately and provides parameters to predict the corresponding dimensional changes. The resulting model can improve the engineering predictions of long-term drying shrinkage accounting also for the history dependence of strain induced by hysteresis. Alternative strategies for quantitative analyses of the * Corresponding author. hmj@mit.edu PHYSICAL REVIEW APPLIED 3, 064009 (2015) 2331-7019=15=3(6)=064009 (17) 064009-1

No abstract

An optimization of the carbon force centers connecting the aromatic ring with methyl groups was performed to extend the new Anisotropic United Atoms (AUA) intermolecular potential recently developed for benzene molecules that include electrostatic interactions [Bonnaud, P.; Nieto-Draghi, C.; Ungerer, P. J. Phys. Chem. B 2007, 111, 3730−3741]. The new potential has been tested by computing thermodynamic and structural properties of methylbenzenes using Gibbs Ensemble, NPT Monte Carlo simulations. A comprehensive comparison of the new model is given with other intermolecular potentials taken from the literature. Overall thermodynamic and structural properties of our optimized model are in very good agreement with experimental data, represented by an absolute average deviation of <1% for liquid density, <2% for vaporization enthalpy, and <18% for vapor pressure. The new model also provides a coherent representation of the liquid structure, as revealed by three-dimensional spatial density functions and carbon−carbon radial distribution function. Preliminary calculations of shear viscosity by equilibrium molecular dynamics show a better agreement with experimental data as a result of introducing electrostatic charges in the AUA potential.

An optimization including electrostatic interactions has been performed for the parameters of an anisotropic united atoms intermolecular potential for benzene for thermodynamic and transport property prediction using Gibbs ensemble, isothermal-isobaric (NPT) Monte Carlo, and molecular dynamic simulations. The optimization procedure is based on the minimization of a dimensionless error criterion incorporating various thermodynamic data (saturation pressure, vaporization enthalpy, and liquid density) at ambient conditions and at 350 and 450 K. A comprehensive comparison of the new model is given with other intermolecular potentials taken from the literature. Overall thermodynamic, structural, reorientational, and translational dynamic properties of our optimized model are in very good agreement with experimental data. The new model also provides a good representation of the liquid structure, as revealed by three-dimensional spatial density functions and carbon-carbon radial distribution function. Shear viscosity variations with temperature and pressure are very well reproduced, revealing a significant improvement with respect to nonpolar models.

Calcium-silicate-hydrate (or C-S-H), an inosilicate, is the major binding phase in cement pastes and concretes and a porous hydrated material made up of a percolated and dense network of crystalline nanoparticles of a mean apparent spherical diameter of ∼5 nm that are each stacks of multiple C-S-H layers. Interaction forces between these nanoparticles are at the origin of C-S-H chemical, physical, and mechanical properties at the meso- and macroscales. These particle interactions and the resulting properties may be affected significantly by nanoparticle density and environmental conditions such as the temperature, relative humidity, or concentration of chemical species in the bulk solution. In this study, we combined grand canonical Monte Carlo simulations and an extension of the mean force integration method to derive the pair potentials. This approach enables realistic simulation of the physical environment surrounding the C-S-H particles. We thus constructed the pair potentials for C-S-H nanoparticles of defined chemical stoichiometry at 10% relative humidity (RH), varying the relative crystallographic orientations at a constant particle density of ρpart ∼ 2.21 mmol L(-1). We found that cohesion between nanoparticles is affected strongly by both the aspect ratio and the crystallographic misorientation of interacting particles. This method and the findings underscore the importance of accounting for relative dimensions and orientation among C-S-H nanoparticles in descriptions of physical and simulated multiparticle aggregates or mesoscale systems.

Materials enabling impact-energy absorption of high-velocity projectiles are of great interest for applications like aerospace. In such a frame, shear thickening fluids were found very useful. Here, we investigated nanorheological properties of neat and aqueous polyelectrolytes of low molecular weights containing poly­([2-(methacryloyloxy) ethyl] trimethyl ammonium) as polycations and poly­(acrylamide-co-acrylic acid) as polyanions. Results were compared with pure water. We employed nonequilibrium molecular dynamics with the SLLOD algorithm to compute the viscosity at various shear rates. Systems containing polyelectrolytes exhibit shear thickening. The analysis of molecular configurations revealed a strong disruption of the ionic structure and more clusters with smaller sizes on increasing the shear rate. Potential energies showed that shear thickening originates from an increase in intramolecular and van der Waals interactions resulting from the increasing difficultly of polyelectrolyte-based systems to relax at high shear rates. Our method and findings underscore the importance of accounting for the molecular scale in the design of materials absorbing the impact energy efficiently.