Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Solvated Electron in Acetonitrile: Radiation Yield, Absorption Spectrum, and Equilibrium between Cavity- and Solvent-Localized States

The equilibrium between a solvent cavity-localized electron, ecav-, and a dimeric solvent anion, (CH 3 CN)$_{2}^{•-}$, which are the two lowest energy states of the solvated electron in acetonitrile, has been investigated by pulse radiolysis at 233-353 K. The enthalpy and entropy for the e$_{cav}^{-}$ to (CH 3 CN)$_{2}^{•-}$ conversion amount to -11.2 ± 0.3 kcal/mol and -39.3 ± 1.2 cal/(mol K), corresponding to a 0.44 ± 0.35 equilibrium constant at 25 °C. The radiation yield of the solvated electron has been quantified using a Co(II) macrocycle that scavenges electrons with a 1.55 × 10 11 M -1 s -1 rate constant. The apparent yield increases without saturation over the attainable scavenger concentration range, reaching 2.8 per 100 eV; this value represents the lower limit for the acetonitrile ionization yield in pulse radiolysis. The apparent molar absorption coefficient of (20.8 ± 1.5) × 10 3 M -1 cm -1 at 1450 nm and 20 °C for the solvated electron and individual vis-near-infrared (NIR) absorption spectra of e$_{cav}^{-}$ and (CH 3 CN)$_{2}^{•-}$ are derived from the data. Variances with previous reports are thoroughly discussed. All together, these results resolve several controversies concerning the solvated electron properties in acetonitrile and furnish requisite data for quantitative pulse radiolysis investigations in this commonly used solvent.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Quantifying Energetic and Entropic Pathways in Molecular Systems

When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Lastly, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Energetics and Structure of Nickel Atoms and Nanoparticles on MgO(100)

The growth morphology and interfacial energetics of vapor deposited Ni on the MgO(100) surface at 300 and 100 K have been studied using single crystal adsorption calorimetry (SCAC), He + low-energy ion scattering spectroscopy (LEIS), X-ray photoelectron spectroscopy (XPS), and low-energy electron diffraction (LEED). At 300 K, the Ni atoms grow as three-dimensional nanoparticles with a saturation number density of 5 × 10 16 particles/m 2 . The differential heat of adsorption at 300 K increases rapidly with coverage, from 276 (initially) to 311 kJ/mol by 0.4 ML. Thereafter, it slowly increases asymptotically to the sublimation enthalpy of bulk Ni (430 kJ/ml) by 9 ML. The Ni 2p 3/2 XPS peak binding energy at 300 K is initially (i.e., at 0.16 ML) 1.4 eV higher than that for bulk Ni(solid), but it decreases to that value at high coverage. The Ni atoms form a metastable phase at 300 K when in nanoparticles with diameter <2.5 nm, and the adhesion energy of such Ni nanoparticles to MgO(100) was found to be 3.05 J/m 2 . At 100 K, the Ni atoms form single adatoms and then 0.17 nm thick 2D islands at low coverage with fewer Ni-Ni bonds compared to the Ni nanoparticles formed at 300 K. Thus, the initial heat (i.e., for the first ~0.03 ML) is 148 kJ/mol at 100 K, 128 kJ/mol lower than at 300 K, and remains lower for the 2D islands. With increasing coverage at 100 K, the tiny 2D Ni islands grow in size to cover nearly the entire surface before thickening. Lastly, the XPS Ni 2p 3/2 peak binding energy for 0.21 ML Ni on MgO(100) at 100 K is 2.2 eV higher than that for bulk Ni(solid), suggesting charge transfer from Ni to MgO(100) and formation of Ni 2+ at very low coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Effect of Background Electrolyte Composition on the Interfacial Formation of Th(IV) Nanoparticles on the Muscovite (001) Basal Plane

Understanding the impact of actinide nanoparticle (NP) formation is important to assess radionuclide mobility in the environment. We combined surface X-ray diffraction (SXRD) and in situ AFM to investigate the previously reported unusual electrolyte effects on Th uptake on mica. At low [Th] (0.1 mM), interfacial structures show a broad Th electron density (~50 Å). A linear decrease of Th uptake with decreasing hydration enthalpy of the electrolyte cation (Li+, K+, NH4+, and Cs+) indicates a competitive effect between Th and the electrolyte cation. Na+ is a clear outlier from this trend. In situ AFM imaging confirms the results. Particles show a vertical size of ~1–2 nm and larger lateral dimensions of ~10–20 nm, which is typical for particles formed at interfaces (heterogeneous nucleation). At high [Th] = 1 and 3 mM, all investigated electrolytes (ACl, A = Li+, Na+, K+) show similar Th uptake, indicating a much smaller impact of electrolyte composition. The interfacial structures are dominated by a high Th loading at a distinct distance (~6.5 Å) from the surface. Therefore, the main retention mechanism at high [Th] is suggested to be the sorption of Th NPs aggregated from Th oligomers present in solution (homogeneous nucleation).

Neumann, J↗
DOE OSTI2021

Prediction of Transition-State Scaling Relationships and Universal Transition-State Vibrational and Entropic Correlations for Dehydrogenations

Linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) or transition-state scaling (TSS) relations aid with the prediction of electronic energies. However, temperature effects and pre-exponentials are often taken as constants across metal surfaces or a homologous series. Vibrational scaling relationships (VSRs) offer a way to determine such parameters. Transition-state VSRs (TSVSRs) between local minima and transition states of AH X (A = C, N, O) surface diffusions correlate with BEP relations and broaden to thermochemical property scaling. Using density functional theory, we extend TSVSRs to AHX dehydrogenation reactions on transition-metal surfaces, relating vibrational modes of local minima to transition states. We first predict the slopes of the TSS relations by incorporating bond angles using the Slater–Koster structure factors and hybridization through crystal orbital overlap population analysis and energy overlap integrals between adsorbates and metal surfaces. Additionally, we uncover universal thermochemical property scaling, enabling the estimation of entropies and temperature corrections to enthalpies across a homologous series. Here, we demonstrate both significant vibrational corrections in reactions with low intrinsic electronic barriers and considerable variation in the pre-exponential of a simple dehydrogenation reaction across metals and AH X adsorbates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Elucidating the Influence of Metal Surface Composition on Organic Adsorbate Binding Using Active Particle Dynamics

The adsorption strengths of organic compounds on metal surfaces are sensitive to the metal composition, and they play a central role in many catalytic reactions, helping to control the coverage of the reactant and altering the overall reaction rate. While adsorption energies are straightforward to measure and calculate in vacuum and gas-phase environments, adsorption energetics can be dramatically altered by the presence of solvent in liquid-phase reactions. However, the effects of metal composition on binding strengths in a liquid environment are less well understood, primarily due to the difficulty of accurate in situ measurements of organic binding on metal surfaces in the liquid phase. Here, we utilize the motion of active particles in water to probe the adsorption energies of an organic adsorbate (furfural) on a range of metal surfaces (pure Pd, pure Pt, and four PdAu alloy compositions) to elucidate the effect of metal composition. Janus particles with catalytic caps of particular metal compositions all exhibited active motion resulting from consumption of H 2 O 2 ; adsorbate binding was inferred through the decrease in the velocity of active motion and was modeled by a Langmuir adsorption isotherm. The measured adsorption affinities were used to extract the adsorption enthalpy of furfural on the different metals. The Pd surface was found to bind furfural more strongly than the Pt surface by some 10 kJ/mol. Furthermore, the adsorption of furfural on the alloys was found to increase monotonically in magnitude with Pd content. As a result, the data reported herein aid the development of accurate understanding of organic adsorption in the presence of solvent and the role of the metal surface in tuning adsorption strengths to optimize catalytic processes in the liquid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Polyethylene in Dead-End Silica Nanopores: Forces and Mobility from Non-Equilibrium Statistical Mechanics and Exchange Spectroscopy Nuclear Magnetic Resonance

Billions of tons of plastic have been produced, and only a small fraction of this has been recycled. Tennakoon et al. [Nature Catalysis 3, 893 (2020)] developed a catalyst that repeatedly cleaves C 10 –C 30 hydrocarbons from the end of a polyethylene chain. The reaction occurs at a Pt nanoparticle at the base of a cylindrical silica mesopore with a diameter of 2 nm and a length of 110 nm. Portions of the polymer situated inside the pore can be differentiated from those outside using 13 C nuclear magnetic resonance (NMR), allowing the dynamics and extent of polymer threading to be monitored using two-dimensional (2D) exchange spectroscopy NMR. Here we construct a Fokker–Planck equation for the polymer dynamics by assuming a reptation diffusivity and a graduated adsorption free energy that depends linearly on the depth of polymer penetration in the pore. The solutions allow us to predict the intensities of the 2D NMR resonances as a function of time. We use the solutions to extract a polymer diffusivity at each temperature and estimate the per-segment desorption free energy, enthalpy, and entropy. Random and systematic errors are examined to test key assumptions in the theory and interpretation of the experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ethylene Hydrogenation Molecular Mechanism on MoC y Nanoparticles

Ethylene hydrogenation catalyzed by MoC y nanoparticles has been studied by means of density functional theory methods and several models. These include MetCar (Mo 8 C 12 ), Nanocube (Mo 14 C 13 ), and Mo 12 C 12 nanoparticles as representatives of experimental MoC y nanostructures. The effect of hydrogen coverage has been studied in detail by considering low-, intermediate-, and high-hydrogen regimes. The calculated enthalpy and energy barriers show that ethylene hydrogenation is feasible on the MetCar, Mo 12 C 12 , and Nanocube but at low, medium, and high hydrogen coverages, respectively. An additional step, related to the H* migration from a Mo to a C site in the nanoparticle, has been found to be the key to establishing the best hydrogenation system. In most cases, the reactions are exothermic, featuring low hydrogenation energy barriers, especially for the Nanocube at high hydrogen coverage. In addition, the calculated adsorption Gibbs free energy shows that, for this system, the C 2 H 4 adsorption is feasible in the 300–400 K temperature range and pressures from 10 –10 to 2 atm. For the hydrogenation steps, calculated transition state theory rates show that the overall process is limited by the first hydrogenation step (C 2 H 4 → C 2 H 5 ) at temperatures of 330–400 K. However, at the lower temperatures of 300–320 K, the reaction rates are comparable for the two steps. The present results indicate that the Mo 14 C 13 Nanocube models of MoC y nanoparticles exhibit appropriate thermodynamic and kinetic features to catalyze ethylene hydrogenation at a high-hydrogen-coverage regime. The present findings provide a basis for understanding the chemistry of active MoC y catalysts, suggest appropriate working conditions for the reaction to proceed, and provide a basis for future experimental studies.

03 NATURAL GAS↗
DOE OSTI2023

Procedures for Assessing the Stability of Proposed Topological Materials

Here, we investigate the stability of MnPb 2 Bi 2 Te 6 (MPBT), which is predicted to be a magnetic topological insulator (TI), using density functional theory calculations. Our analysis includes various measures such as enthalpies of formation, Helmholtz free energies, defect formation energies, and dynamical stability. Our thermodynamic analysis shows that the phonon contribution to the energy gain from finite temperature is estimated to be less than 10 meV/atom, which may not be sufficient to stabilize MPBT at high temperatures, even with the most favorable reactions starting from binaries. While MPBT is generally robust against the formation of various defects, we find that the anti-site defect formation of MnPb is the most likely to occur, with a corresponding energy less than 60 meV. This can be attributed to the significant energy cost from compressive strain at the PbTe layer. Our findings suggest that MPBT is on the brink of stability in terms of thermodynamics and defect formation, underscoring the importance of conducting systematic analyses of the stability of proposed TIs, including MPBT, for their practical utilization. This study offers valuable insights into the design and synthesis of desirable magnetic TI materials with robust stabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Local Order in Liquid Gallium–Indium Alloys

Liquid metals such as eutectic Ga–In alloys have low melting points and low toxicity and are used in catalysis and micro-robotics. This study investigates the local atomic structure of liquid gallium-indium alloys by a combination of density measurements, diffraction data, and Monte-Carlo simulation via the empirical potential structure refinement approach. A high-Q shoulder observed in liquid Ga is related to structural rearrangements in the second coordination shell. Structure analysis found coordination environments close to a random distribution for eutectic Ga–In alloy, while electronic effects appear to dominate the mixing enthalpy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Anions on the Delithiation of [Li–Al] Layered Double Hydroxides: Thermodynamic Insights

The utilization of lithium aluminum [Li–Al] layered double hydroxides (LDHs) is explored for direct lithium extraction (DLE) from geothermal brines and minerals. Following extraction, the Li + ions need to be removed from the LDH structure and converted into LiOH or Li 2 CO 3 products, making them suitable for battery applications. The research investigates the delithiation of [Li–Al-X] LDHs (where X = Cl – , OH – , and SO 4 2– ), which were synthesized and dried under different conditions. Herein the study aims to understand how the choice of anions and drying conditions affects the delithiation process. To determine the stability of these [Li–Al] LDHs, high-temperature oxide melt solution calorimetry is employed. The results reveal that the stability of LDHs varies significantly based on postdrying conditions, anion choice, and the water content in the interlayer. The order of stability, as indicated by the enthalpy of formation values, is found to be OH-LDH-O > OH-LDH-A > SO 4 -LDH-O > SO4-LDH-A > Cl-LDH-O > Cl-LDH-A. This trend is attributed to the interactions among the interlayer species and the metal hydroxide layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Synthesis and Properties of the Helium Clathrate and Defect Perovskite [He 2– x $^\square$ x ][CaNb]F 6

The defect double perovskite [He 2– x $^\square$ x ][CaNb]F 6 , with helium on its A-site, can be prepared by the insertion of helium into ReO 3 -type CaNbF 6 at high pressure. Upon cooling from 300 to 100 K under 0.4 GPa helium, ~60% of the A-sites become occupied. Helium uptake was quantified by both neutron powder diffraction and gas insertion and release measurements. After the conversion of gauge pressure to fugacity, the uptake of helium by CaNbF 6 can be described by a Langmuir isotherm. The enthalpy of absorption for helium in [He 2– x $^\square$ x ][CaNb]F 6 is estimated to be ~+3(1) kJ mol –1 , implying that its formation is entropically favored. Helium is able to diffuse through the material on a time scale of minutes at temperatures down to ~150 K but is trapped at 100 K and below. The insertion of helium into CaNbF 6 reduces the magnitude of its negative thermal expansion, increases the bulk modulus, and modifies its phase behavior. On compressing pristine CaNbF 6 , at 50 and 100 K, a cubic (Fm$\overline{3}$m) to rhombohedral (R$\overline{3}$) phase transition was observed at <0.20 GPa. However, a helium-containing sample remained cubic at 0.4 GPa and 50 K. CaNbF 6 , compressed in helium at room temperature, remained cubic to >3.7 GPa, the limit of our X-ray diffraction measurements, in contrast to prior reports that upon compression in a nonpenetrating medium, a phase transition is detected at ~0.4 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Impact of Anions and Water Content on [Li–Al] Layered Double-Hydroxide Stability

[Li–Al] layered double hydroxides (LDHs) are compounds with potential as sorbents for lithium extraction from brine solutions. Here, in this work, heat capacities were measured from approximately 2.5 to 300 K for six [Li–Al] LDHs with differing anions (Cl – , OH – , and SO 4 2– ) and water content (denoted A for air-dried and O for oven-dried). These measurements were used to calculate the standard entropy at 298.15 K, and the results were combined with previously performed enthalpy measurements to calculate Gibbs energies of formation from the binary compounds. The calculated order of stability based on Gibbs energies of formation was Cl-LDH-O > OH-LDH-O > Cl-LDH-A > SO 4 -LDH-O > SO 4 -LDH-A > OH-LDH-A. Results support previous findings that higher water content generally raises the Gibbs energy of the LDH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Computational Optimization of Room Temperature Usable Capacity for Hydrogen Storage in MFU-4-Type Metal–Organic Frameworks via Pairwise Metal Substitutions

The efficient storage of hydrogen is a critical challenge in the quest for sustainable energy solutions. Current adsorbent-based methods achieve satisfactory storage densities predominantly under cryogenic temperatures and/or high pressures, which imposes problems with cost-efficient and safe implementation of this technology. Materials that can bind hydrogen gas reversibly at ambient temperatures and more moderate pressures could play a pivotal role in enabling hydrogen-powered technologies. In this study, we use reliable computational modeling to investigate two synthetically feasible paths for tuning the enthalpy of H2 binding in MFU-4-type metal–organic frameworks (MOFs), aiming to maximize usable capacity. This study examines MIM4 IICl3(bta)6 (bta– = benzotriazolate) Kuratowski-type clusters as a model for strong binding sites in MFU-4l frameworks. We systematically evaluate the impact of separately tuning the central MII metal ion (which plays a structural role) and the peripheral MI metal ion (which binds the substrate) on the energetics of H2 binding. Our computational study reveals that H2 binding at an MI site mostly follows the trend AgI < CuI < NiI < CoI < AuI while a larger central MII site generally weakens the H2 binding at a MI site. Importantly, we have identified three new combinations of MI and MII to achieve high fractional usable capacities of the total H2 adsorbed under a pressure swing from 5 to 100 bar at room temperature. Additionally, we examine the nature of the binding interaction between the peripheral metal atom and the hydrogen molecule. While charge transfer predominantly induces this interaction, for several atom combinations, a change in the polarization (associated with variations in the ionic radius of the MI binding atom) is another important factor for adjusting the strength of the interaction. We suggest that the proposed compositions of Kuratowski-type clusters are highly desirable synthetic targets for future laboratory study.

Tkachenko, Nikolay V↗
DOE OSTI2025

Adsorption Studies of Dilute Krypton and Xenon from Nitrogen on SBMOF-1 and Activated Charcoal for Applications in Isotope Harvesting

Adsorptive partitioning of dilute krypton (Kr) and xenon (Xe) onto Stony Brook Metal–Organic Framework (SBMOF-1) and activated charcoal (AC) from carrier nitrogen was experimentally measured at temperatures ranging from 195 to 293 K. For this purpose, a closed-loop system for gas adsorption experiments was developed. From the Kr adsorption measurements, the adsorption equilibrium constant for Kr on SBMOF-1 was calculated, yielding a value for the enthalpy of adsorption of −19 ± 1 kJ·mol –1 . The partition coefficients were utilized to estimate the extraction rates of 76 Kr, 77 Kr, and 122 Xe isotopes during isotope harvesting at the Facility for Rare Isotope Beams (FRIB). We conclude that both materials showed promising results for the extraction of noble gases from FRIB effluents using temperature swing adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Probing Conformational Evolution and Associated Dynamics of Mg(N(SO 2 CF 3 ) 2 ) 2 ·Dimethoxyethane Adduct Using Solid-State 19 F and 1 H NMR

Bis(trifluoromethanesulfonimide) or TFSI is widely used as a counter anion in electrolyte design due to its structural flexibility and chemical stability. Here we studied the conformational variations and associated dynamics of TFSI in adduct of Mg(TFSI) 2 with dimethoxyethane (DME), a solvate crystalline material using solid-state 1 H and 19 F NMR. TFSI molecular motion in this solvate structure falls within the timescale of the 19 F NMR experiment, yielding spectroscopic signatures for unique TFSI conformers under the coordination environment of Mg 2+ cation. Within the temperature range of -5 to 82 °C, we observe nine distinct TFSI sites in both crystalline and disordered regions using 19 F NMR, reflecting complexity of structural and dynamics of TFS1 anions within solvate structure. The four distinguishable sites in the disordered region for the two CF 3 groups of the same TFSI molecule are identified using chemical shift analysis. The exchange rate constants from site to site are calculated through variable temperature 19 F NMR and two-dimensional (2D) exchange spectroscopy (EXSY) experiments, along with respective activation enthalpies using Eyring's formulation. The flip rate of CF 3 around the S-C bond is estimated as ~ 15 s -1 at 8 °C with ΔH ≠ ~ 22 kJ/mol, but the rotation of the entire TFSI is 4.8 s -1 at 8 °C with a significantly greater ΔH ≠ = 98 ± 10 kJ/mol. Furthermore, the slow conversion of trans to cis conformers at a lower temperature (T ≤ 1 °C) in the crystalline region is monitored, with a conversion rate of ~ 2 x 10 -5 s -1 at -5 °C. Density functional theory (DFT)-based calculations were performed to support further the assignment of experimental chemical shifts, and the activation energy E a = 21.1 kJ/mol obtained for the cis to trans transition is consistent with experimental values. The combined set of 19 F and 1 H under both one-dimensional (1D) and 2D NMR methods demonstrated here can be further used for examining electrode-electrolyte interfaces to probe the motions of various constituents that can enable detailed studies of interfacial processes and dynamics. Ultimately, such studies will aid in the design and discovery of interfacial constructs in which directed defect chemistry, chemical moiety distribution, and nanostructure are employed to drive efficient charge transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗