Design Principles Guiding Solvent Size Selection in ZIF-Based Type 3 Porous Liquids for Permanent Porosity

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.

Modeling single-molecule stretching experiments using statistical thermodynamics

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.

Experimental and Computational Mechanisms that Govern Long-Term Stability of CO2-Adsorbed ZIF-8-Based Porous Liquids

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 CO₂ 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 N₂ or air. However, for PLs aged in a CO₂ 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 CO₂ on the PL solvent mixture, it was identified that the basic environment of the PL caused ethylene glycol to react with CO₂ 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.

Effect of Linker Structure and Functionalization on Secondary Gas Formation in Metal-Organic Frameworks

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 (H₂O, NO₂, and SO₂) 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, -NH₂, 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.

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

Dramatic Enhancement of Rare-Earth Metal-Organic Framework Stability Via Metal Cluster Fluorination

Rare-earth polynuclear metal-organic frameworks (RE-MOFs) have demonstrated high durability for caustic acid gas adsorption and separation based on gas adsorption to the metal clusters. The metal clusters in the RE-MOFs traditionally contain RE metals bound by μ₃-OH groups connected via organic linkers. Recent studies have suggested that these hydroxyl groups could be replaced by fluorine atoms during synthesis that includes a fluorine-containing modulator. Here, a combined modeling and experimental study was undertaken to elucidate the role of metal cluster fluorination on the thermodynamic stability, structure, and gas adsorption properties of RE-MOFs. Through systematic density-functional theory calculations, fluorinated clusters were found to be thermodynamically more stable than hydroxylated clusters by up to 8-16 kJ/mol per atom for 100% fluorination. The extent of fluorination in the metal clusters was validated through a ¹⁹F NMR characterization of 2,5-dihydroxyterepthalic acid (Y-DOBDC) MOF synthesized with a fluorine-containing modulator. ¹⁹F magic-angle spinning NMR identified two primary peaks in the isotropic chemical shift (δ_iso) spectra located at -64.2 and -69.6 ppm, matching calculated ¹⁹F NMR δ_iso peaks at -63.0 and -70.0 ppm for fluorinated systems. Calculations also indicate that fluorination of the Y-DOBDC MOF had negligible effects on the acid gas (SO₂, NO₂, H₂O) binding energies, which decreased by only ∼4 kJ/mol for the 100% fluorinated structure relative to the hydroxylated structure. Additionally, fluorination did not change the relative gas binding strengths (SO₂ > H₂O > NO₂). Therefore, for the first time the presence of fluorine in the metal clusters was found to significantly stabilize RE-MOFs without changing their acid-gas adsorption properties.

Porous Liquids: Computational Design for Targeted Gas Adsorption

In this Perspective, we present the unique gas adsorption capabilities of porous liquids (PLs) and the value of complex computational methods in the design of PL compositions. Traditionally, liquids only contain transient pore space between molecules that limit long-term gas capture. However, PLs are stable fluids that that contain permanent porosity due to the combination of a rigid porous host structure and a solvent. PLs exhibit remarkable adsorption and separation properties, including increased solubility and selectivity. The unique gas adsorption properties of PLs are based on their structure, which exhibits multiple gas binding sites in the pore and on the cage surface, varying binding mechanisms including hydrogen-bonding and π-π interactions, and selective diffusion in the solvent. Tunable PL compositions will require fundamental investigations of competitive gas binding mechanisms, thermal effects on binding site stability, and the role of nanoconfinement on gas and solvent diffusion that can be accelerated through molecular modeling. With these new insights PLs promise to be an exceptional material class with tunable properties for targeted gas adsorption.

The Lanthanide Contraction Is a Variable

Upon examination of the bond distances of the recently reported series of [Ln(SST)₃(THF)₂] [Ln = lanthanides, SST = tris(trimethylsilyl)siloxide (OSi(SiMe₃)₃), and THF = tetrahydrofuran] compounds, it was found that over the Ln-series (La through Lu), the Ln-O(THF) bond changed by 0.257 Å, whereas the Ln-O(SST) bond varied by 0.164 Å. Examination of all similarly ligated Ln-O(THF) (Ln = La vs Lu) structures available in the Cambridge Structural Database (CSD) revealed that this previously unreported, increased Ln-contraction is pervasive. Further evaluations showed that this enhanced Ln-contraction also occurs for pyridine (py) in the [Ln(SST)₃(py)₂] family as well as the average Ln-N(py) (La vs Lu) structure distances recovered from the CSD. Additional ligands, such as halides (Cl and I) were found to display this enhanced Ln-contraction, while other species (i.e., cyclopentadienide, alkoxide, SST, and dimethyl sulfoxide) yielded a "normal" Ln-contraction (La-L vs Lu-L). Gas-phase electronic structure density functional theory calculations were carried out to evaluate the molecular orbital influence on the Ln-contraction between Ln-O(SST) and Ln-O(THF). The calculated [Ln(SST)₃(THF)₂] structures were found to demonstrate the same capricious Ln-contraction. Based on these studies, one can say that the variability of the Ln-contraction noted in the [Ln(SST)₃(THF)₂] experimental data is due to the different bonding types, ion-ion for the Ln-SST bond versus ion-dipole for the Ln-THF bond.

Design Elements for Enhanced Hydrogen Isotope Separations in Barely Porous Organic Cages

Barely porous organic cages (POCs) successfully separate hydrogen isotopes (H₂/D₂) at temperatures below 100 K. Identifying the mechanisms that control the separation process is key to the design of next-generation hydrogen separation materials. Here, ab initio molecular dynamics (AIMD) simulations are used to elucidate the mechanisms that control D₂ and H₂ separation in barely POCs with varying functionalization. The temperature and pore size dependence were identified, including the selective capture of D₂ in three different CC3 structures (RCC3, CC3-S, and 6ET-RCC3). The temperature versus capture trend was reversed for the 6ET-RCC3 structure, identifying that the D₂ and H₂ escape mechanisms are unique in highly functionalized systems. Analysis of calculated isotope velocities identified effective pore sizes that extend beyond the pore opening distances, resulting in increased capture in minimally functionalized CC3-S and RCC3. In a highly functionalized POC, 6ET-RCC3, higher velocities of the H isotopes were calculated moving through the restricted pore compared to the rest of the system, identifying a unique molecular behavior in the barely nanoporous pore openings. By using AIMD, mechanisms of H₂ and D₂ separation were identified, allowing for the targeted design of future novel materials for hydrogen isotope separation.

Kinetically Controlled Linker Binding in Rare Earth-2,5-Dihydroxyterepthalic Acid Metal-Organic Frameworks and Its Predicted Effects on Acid Gas Adsorption

In the pursuit of highly stable and selective metal-organic frameworks (MOFs) for the adsorption of caustic acid gas species, an entire series of rare earth MOFs have been explored. Each of the MOFs in this series (RE-DOBDC; RE = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu; DOBDC = 2,5-dihydroxyterepthalic acid) was synthesized in the tetragonal space group I4/m. Crystallized MOF samples, specifically Eu-DOBDC, were seen to have a combination of monodentate and bidentate binding when synthesized under typical reaction conditions, resulting in a contortion of the structure. However, extended crystallization times determined that this binding is kinetically controlled and that the monodentate binding option was crystallographically eliminated by extended reaction times at higher temperatures. Furthermore, this series allows for the direct study of the effect of the metal center on the structure of the of the MOF; herein, the lanthanide metal ionic radii contraction across the periodic table results in a reduction of the MOF pore size and lattice parameters. Scanning electron microscopy-energy-dispersive spectroscopy was used to investigate the stages of crystal growth for these RE-DOBDC MOFs. All MOFs, except Er-DOBDC had a minimum of two stages of growth. These analogues were demonstrated by analysis of neutron diffraction (PND) to exhibit a cooperative rotational distortion of the secondary building unit, resulting in two crystallographically distinct linker sublattices. Computational modeling efforts were used to show distinct differences on acid gas (NO₂ and SO₂) binding energies for RE-DOBDC MOFs when comparing the monodentate/bidentate combined linker with the bidentate-only linker crystal structures.

Temperature and Pressure Dependence of Salt-Brine Dihedral Angles in the Subsurface

Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na⁺ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.

Prediction of Reactive Nitrous Acid Formation in Rare-Earth MOFs via ab initio Molecular Dynamics

Reactive gas formation in pores of metal-organic frameworks (MOFs) is a known mechanism of framework destruction; understanding those mechanisms for future durability design is key to next generation adsorbents. Herein, an extensive set of ab initio molecular dynamics (AIMD) simulations are used for the first time to predict competitive adsorption of mixed acid gases (NO2 and H2 O) and the in-pore reaction mechanisms for a series of rare earth (RE)-DOBDC MOFs. Spontaneous formation of nitrous acid (HONO) is identified as a result of deprotonation of the MOF organic linker, DOBDC. The unique DOBDC coordination to the metal clusters allows for proton transfer from the linker to the NO2 without the presence of H2 O and may be a factor in DOBDC MOF durability. This is a previously unreported mechanisms of HONO formation in MOFs. With the presented methodology, prediction of future gas interactions in new nanoporous materials can be achieved.

Stability Evaluation of Candidate Precursors for Chemical Vapor Deposition of Hafnium Diboride (HfB2)

Alternative candidate precursors to [Hf(BH₄)₄] for low-temperature chemical vapor deposition of hafnium diboride (HfB₂) films were identified using density functional theory simulations of molecules with the composition [Hf(BH₄)₂L₂], where L = -OH, -OMe, -O-t-Bu, -NH₂, -N=C=O, -N(Me)₂, and -N(CH₂)₅NH₂ (1-piperidin-2-amine referred to as Pip2A). Disassociation energies (E _D), potential energy surface (PES) scans, ionization potentials, and electron affinities were all calculated to identify the strength of the Hf-L bond and the potential reactivity of the candidate precursor. Ultimately, the low E _D (2.07 eV) of the BH₄ ligand removal from the Hf atom in [Hf(BH₄)₄] was partially attributed to an intermediate state where [Hf(BH₄)₃(H)] and BH₃ is formed. Of the candidate precursors investigated, three exhibited a similar mechanism, but only -Pip2A had a PES scan that indicated binding competitive with [Hf(BH₄)₄], making it a viable candidate for further study.

Magnetic Tunability in RE-DOBDC MOFs via NOx Acid Gas Adsorption

The magnetic susceptibility of NO_x-loaded RE-DOBDC (rare earth (RE): Y, Eu, Tb, Yb; DOBDC: 2,5-dihydroxyterephthalic acid) metal-organic frameworks (MOFs) is unique to the MOF metal center. RE-DOBDC samples were synthesized, activated, and subsequently exposed to humid NO_x. Each NO_x-loaded MOF was characterized by powder X-ray diffraction, and the magnetic characteristics were probed by using a VersaLab vibrating sample magnetometer (VSM). Lanthanide-containing RE-DOBDC (Eu, Tb, Yb) are paramagnetic with a reduction in paramagnetism upon adsorption of NO_x. Y-DOBDC has a diamagnetic moment with a slight reduction upon adsorption of NO_x. The magnetic susceptibility of the MOF is determined by the magnetism imparted by the framework metal center. The electronic population of orbitals contributes to determining the extent of magnetism and change with NO_x (electron acceptor) adsorption. Eu-DOBDC results in the largest mass magnetization change upon adsorption of NO_x due to more available unpaired f electrons. Experimental changes in magnetic moment were supported by density functional theory (DFT) simulations of NO_x adsorbed in lanthanide Eu-DOBDC and transition metal Y-DOBDC MOFs.

Surface Energies and Structure of Salt-Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt-brine and salt-salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. Here, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt-brine interface to provide insight into the stability of salt systems. Mixed NaCl-KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt-brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2-3 times (5.0 M systems). Size differences in Na⁺ and K⁺ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt-brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

Computational and Experimental 1H-NMR Study of Hydrated Mg-Based Minerals

Magnesium oxide (MgO) can convert to different magnesium-containing compounds depending on exposure and environmental conditions. Many MgO-based phases contain hydrated species allowing ¹H-nuclear magnetic resonance (NMR) spectroscopy to be used in the characterization and quantification of proton-containing phases; however, surprisingly limited examples have been reported. Here, ¹H-magic angle spinning (MAS) NMR spectra of select Mg-based minerals are presented and assigned. These experimental results are combined with computational NMR density functional theory (DFT) periodic calculations to calibrate the predicted chemical shielding results. This correlation is then used to predict the NMR shielding for a series of different MgO hydroxide, magnesium chloride hydrate, magnesium perchlorate, and magnesium cement compounds to aid in the future assignment of ¹H-NMR spectra for complex Mg phases.

Tuned Hydrogen Bonding in Rare-Earth Metal-Organic Frameworks for Design of Optical and Electronic Properties: An Exemplar Study of Y-2,5-Dihydroxyterephthalic Acid

Organic linkers in metal-organic framework (MOF) materials exhibit differences in hydrogen bonding (H-bonding), which can alter the geometric, electronic, and optical properties of the MOF. Density functional theory (DFT) simulations were performed on a photoluminescent Y-2,5-dihydroxyterephthalic acid (DOBDC) MOF with H-bonding concentrations between 0 and 100%; the H-bonds were located on both bidentate- and monodentate-bound DOBDC linkers. At 0% H-bond concentration in the framework, the lattice parameters contracted, the density increased, and simulated X-ray diffraction patterns shifted. Comparison with published experimental data identified that Y-DOBDC MOF structures must have a degree of H-bond concentration. The concentration of H-bonds in the system shifted the calculated band gap energy from 2.25 eV at 100% to 3.00 eV at 0%. The band gap energies also indicate a distinction of H-bonds formed on bidentate-coordinated linkers compared to those on monodentate linkers. Additionally, when the calculated optical spectra are compared with experimental data, the ligand-to-ligand charge-transfer luminescence in Y-DOBDC MOFs is expected to result from an average of 20-40% H-bonding with at least 50% of the bidentate linkers containing H-bonding. Therefore, the type of H-bonding within the DOBDC linker determines the electronic structure and the optical absorption of the MOF framework structure. Tuning of the H-bonding in rare-earth MOFs provides an opportunity to control the specific optical and adsorption properties of the MOF framework on the basis of reactions between the linker and the environment.

NOx Adsorption and Optical Detection in Rare Earth Metal-Organic Frameworks

Acid gases (e.g., NO_x and SO_x), commonly found in complex chemical and petrochemical streams, require material development for their selective adsorption and removal. Here, we report the NO_x adsorption properties in a family of rare earth (RE) metal-organic frameworks (MOFs) materials. Fundamental understanding of the structure-property relationship of NO_x adsorption in the RE-DOBDC materials platform was sought via a combined experimental and molecular modeling study. No structural change was noted following humid NO_x exposure. Density functional theory (DFT) simulations indicated that H₂O has a stronger affinity to bind with the metal center than NO₂, while NO₂ preferentially binds with the DOBDC ligands. Further modeling results indicate no change in binding energy across the RE elements investigated. Also, stabilization of the NO₂ and H₂O molecules following adsorption was noted, predicted to be due to hydrogen bonding between the framework ligands and the molecules and nanoconfinement within the MOF structure. This interaction also caused distinct changes in emission spectra, identified experimentally. Calculations indicated that this is due to the adsorption of NO₂ molecules onto the DOBDC ligand altering the electronic transitions and the resulting photoluminescent properties, a feature that has potential applications in future sensing technologies.

Structure and electronic properties of rare earth DOBDC metal-organic-frameworks

Here, we apply density functional theory (DFT) to investigate rare-earth metal organic frameworks (RE-MOFs), RE12(μ3-OH)16(C8O6H4)8(C8O6H5)4 (RE = Y, Eu, Tb, Yb), and characterize the level of theory needed to accurately predict structural and electronic properties in MOF materials with 4f-electrons. A two-step calculation approach of geometry optimization with spin-restricted DFT and large core potential (LCPs), and detailed electronic structures with spin-unrestricted DFT with a full valence potential + Hubbard U correction is investigated. Spin-restricted DFT with LCPs resulted in good agreement between experimental lattice parameters and optimized geometries, while a full valence potential is necessary for accurate representation of the electronic structure. The electronic structure of Eu-DOBDC MOF indicated a strong dependence on the treatment of highly localized 4f-electrons and spin polarization, as well as variation within a range of Hubbard corrections (U = 1-9 eV). For Hubbard corrected spin-unrestricted calculations, a U value of 1-4 eV maintains the non-metallic character of the band gap with slight deviations in f-orbital energetics. When compared with experimentally reported results, the importance of the full valence calculation and the Hubbard correction in correctly predicting the electronic structure is highlighted.

An atomic-scale evaluation of the fracture toughness of silica glass

Using an atomistic technique consistent with continuum balance laws and drawing on classical fracture mechanics theory, we estimate the resistance to fracture propagation of amorphous silica. We discuss correspondence and deviations from classical linear elastic fracture mechanics theory including size dependence, rigid/floppy modes of deformation, and the effects of surface energy and stress.

Synthesis and characterization of thallium–salen derivatives for use as underground fluid flow tracers

[Tl₂(salo-Bu^t)] (1, shown) and [Tl₂(saloPh-Bu^t)] (2) were synthesized and characterized for use as monitors for deep subterranean fluid flows.

Surface Structure and Stability of Partially Hydroxylated Silica Surfaces

Surface energies of silicates influence crack propagation during brittle fracture and decrease with surface relaxation caused by annealing and hydroxylation. Molecular-level simulations are particularly suited for the investigation of surface processes. In this work, classical MD simulations of silica surfaces are performed with two force fields (ClayFF and ReaxFF) to investigate the effect of force field reactivity on surface structure and energy as a function of surface hydroxylation. An unhydroxylated fracture surface energy of 5.1 J/m² is calculated with the ClayFF force field, and 2.0 J/m² is calculated for the ReaxFF force field. The ClayFF surface energies are consistent with the experimental results from double cantilever beam fracture tests (4.5 J/m²), whereas ReaxFF underestimated these surface energies. Surface relaxation via annealing and hydroxylation was performed by creating a low-energy equilibrium surface. Annealing condensed neighboring siloxane bonds increased the surface connectivity, and decreased the surface energies by 0.2 J/m² for ClayFF and 0.8 J/m² for ReaxFF. Posthydroxylation surface energies decreased further to 4.6 J/m² with the ClayFF force field and to 0.2 J/m² with the ReaxFF force field. Experimental equilibrium surface energies are ∼0.35 J/m², consistent with the ReaxFF force field. Although neither force field was capable of replicating both the fracture and equilibrium surface energies reported from experiment, each was consistent with one of these conditions. Therefore, future computational investigations that rely on accurate surface energy values should consider the surface state of the system and select the appropriate force field.