Publications

Abstract not provided.

Risks associated with carbonation are a key limitation to greater replacement levels of ordinary portland cement (OPC) by supplementary cementitious materials (SCMs). The addition of pozzolanic SCMs in OPC alters the hydrate assemblage by forming phases like calcium-(alumina)-silicate-hydrate (C-(A)-S-H). The objective of the present study was to elucidate how such changes in hydrate assemblage influence the chemical mechanisms of carbonation in a realistic OPC system. Here, we show that synthetic zeolite Y (faujasite) is a highly reactive pozzolan in OPC that reduces the calcium content of hydration products via prompt consumption of calcium hydroxide from the evolving phase assemblage prior to CO2 exposure. Suppression of portlandite at moderate to high zeolite Y content led to a more damaging mechanism of carbonation by disrupting the formation of a passivating carbonate layer. Without this layer, carbonation depth and CO2 uptake are increased. Binders containing 12–18% zeolite Y by volume consumed all the calcium hydroxide from OPC during hydration and reduced the Ca/(Si+Al) ratio of the amorphous products to near 0.67. In these cases, higher carbonation depths were observed after exposure to ambient air with decalcification of C-(A)-S-H as the main source of CO2 buffering. Binders with either 0% or 4% zeolite Y contained calcium hydroxide in the hydrated microstructure, had higher Ca/(Si+Al) ratios, and formed a calcite-rich passivation layer that halted deep carbonation. Although the carbonated layer in the samples with 12% and 18% zeolite Y contained 70% and 76% less calcite than the OPC respectively, their higher carbonation depths resulted in total CO2 uptakes that were 12x greater than the OPC sample. Passivation layer formation in samples with calcium hydroxide explains this finding and was further supported by thermodynamic modeling. High Si/Al zeolite additives to OPC should be balanced with the calcium content for optimal carbonation resistance.

Porous liquids (PLs), which are solvent-based systems that contain permanent porosity due to the incorporation of a solid porous host, are of significant interest for the capture of greenhouse gases, including CO2. Type 3 PLs formed by using metal-organic frameworks (MOFs) as the nanoporous host provide a high degree of chemical turnability for gas capture. However, pore aperture fluctuation, such as gate-opening in zeolitic imidazole framework (ZIF) MOFs, complicates the ability to keep the MOF pores available for gas adsorption. Therefore, an understanding of the solvent molecular size required to ensure exclusion from MOFs in ZIF-based Type 3 PLs is needed. Through a combined computational and experimental approach, the solvent-pore accessibility of exemplar MOF ZIF-8 was examined. Density functional theory (DFT) calculations identified that the lowest-energy solvent-ZIF interaction occurred at the pore aperture. Experimental density measurements of ZIF-8 dispersed in various-sized solvents showed that ZIF-8 adsorbed solvent molecules up to 2 Å larger than the crystallographic pore aperture. Density analysis of ZIF dispersions was further applied to a series of possible ZIF-based PLs, including ZIF-67, −69, −71(RHO), and −71(SOD), to examine the structure-property relationships governing solvent exclusion, which identified eight new ZIF-based Type 3 PL compositions. Solvent exclusion was driven by pore aperture expansion across all ZIFs, and the degree of expansion, as well as water exclusion, was influenced by ligand functionalization. Using these results, a design principle was formulated to guide the formation of future ZIF-based Type 3 PLs that ensures solvent-free pores and availability for gas adsorption.

During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na2O-SiO2(NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na2O-SiO2(NS15) composition, the 25% Na2O-SiO2(NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na2O-SiO2 compositions. Coordination analysis of the Na+ ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na+ ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Calcite (CaCO3) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca-O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H+ and/or OH- species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17-0.01 OH/nm2 surface hydroxylation observed on the (104) surface at the end of the simulations. The limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO₂ binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO₂ energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO₂ binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO₂ molecules were unable to diffuse into the ZIF-8 pore, with CO₂ adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO₂ adsorption should be based on the strength of gas binding at the solvated particle surface.

Single-molecule stretching experiments are widely utilized within the fields of physics and chemistry to characterize the mechanics of individual bonds or molecules, as well as chemical reactions. Analytic relations describing these experiments are valuable, and these relations can be obtained through the statistical thermodynamics of idealized model systems representing the experiments. Since the specific thermodynamic ensembles manifested by the experiments affect the outcome, primarily for small molecules, the stretching device must be included in the idealized model system. Though the model for the stretched molecule might be exactly solvable, including the device in the model often prevents analytic solutions. In the limit of large or small device stiffness, the isometric or isotensional ensembles can provide effective approximations, but the device effects are missing. Here a dual set of asymptotically correct statistical thermodynamic theories are applied to develop accurate approximations for the full model system that includes both the molecule and the device. The asymptotic theories are first demonstrated to be accurate using the freely jointed chain model and then using molecular dynamics calculations of a single polyethylene chain.

Abstract not provided.

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Porous liquids (PLs) based on the zeolitic imidazole framework ZIF-8 are attractive systems for carbon capture since the hydrophobic ZIF framework can be solvated in aqueous solvent systems without porous host degradation. However, solid ZIF-8 is known to degrade when exposed to CO2 in wet environments, and therefore the long-term stability of ZIF-8-based PLs is unknown. Through aging experiments, the long-term stability of a ZIF-8 PL formed using the water, ethylene glycol, and 2-methylimidazole solvent system was systematically examined, and the mechanisms of degradation were elucidated. The PL was found to be stable for several weeks, with no ZIF framework degradation observed after aging in N2 or air. However, for PLs aged in a CO2 atmosphere, formation of a secondary phase occurred within 1 day from the degradation of the ZIF-8 framework. From the computational and structural evaluation of the effects of CO2 on the PL solvent mixture, it was identified that the basic environment of the PL caused ethylene glycol to react with CO2 forming carbonate species. These carbonate species further react within the PL to degrade ZIF-8. The mechanisms governing this process involves a multistep pathway for PL degradation and lays out a long-term evaluation strategy of PLs for carbon capture. Additionally, it clearly demonstrates the need to examine the reactivity and aging properties of all components in these complex PL systems in order to fully assess their stabilities and lifetimes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Decreasing cost of technologies for direct air capture of carbon can be achieved through the design of new materials with high CO2 selectivity that can be incorporated into existing industrial processes. An emerging class of materials for these applications are porous liquids (PLs). PLs are mixtures of porous hosts and solvents with intrinsic porosity due to steric exclusion of solvent from inside the porous host. It is currently unknown how solvent -porous host interactions affect porous host solubility in the bulk solvent. Here, density functional theory simulations were used to investigate interactions between nine solvents and a CC13 porous organic cage (POC). Calculations identified that solvent molecules were the most stable when placed either inside the CC13 POC or in the pore window compared to interfacial binding sites. Structural changes to the CC13 POC correlated with reported experimental solubilities, including expansion of the CC13 POC with solvent molecule infiltration and expansion or contraction of the pore window. Based on these results, new PL design guidelines should include compositions with (1) high concentrations of POCs with flexible cage structures that can expand when solvated and (2) solvent molecule-POC combinations that contract the pore window during solvent molecule-host binding.

Rare-earth terephthalic acid (BDC)-based metal-organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H2O, NO2, and SO2) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC-gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with −OH, −NH2, and −SH were investigated. The calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.

Abstract not provided.