Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “X-ray binaries”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

X-Ray Polarization of Z-Type Neutron Star Low-Mass X-Ray Binaries II. Spectropolarimetric Analysis

IXPE has provided for the first time detailed energy- and time-resolved X-ray polarimetry of Z-type neutron star low-mass X-ray binaries (NS-LMXBs) as they move along their color-color diagrams (CCDs). These sources can reach the highest polarization observed for NS-LXMBs in the 2–8 keV range when they move along the horizontal branch. In a previous paper, we characterized the spectral state of a sample of Z-sources using the CCD and estimated the polarization with model-independent analysis. Here, we present detailed spectropolarimetric analysis for each source on each branch using data from IXPE, NICER, and NuSTAR. The continuum X-ray emission of all the sources is well described with a combination of thermal accretion disk emission plus a harder Comptonized component. In addition, reflection features, in particular the relativistically broadened Fe Kα line, are observed for our sources, except GX 5–1. For most of the sources and branches, the main contribution to the X-ray emission and polarization is due to Comptonization: moving from the horizontal branch (HB) to the normal branch (NB), the polarization degree (PD) in the 2–8 keV band varies from about 6% to 3–4%, while the PD is loosely constrained in the flaring branch (FB), due to the shorter exposures. These PD values are significantly higher than theoretical expectations for typical spreading or boundary layer configurations. The polarization of the disk is generally lower (below 3%) but still higher than predictions for an electron scattering-dominated, plane-parallel atmosphere above the disk observed at the corresponding inclination. Moreover, the polarization angle (PA) of the disk seems to be significantly misaligned and not perpendicular to that of Comptonization. We find no correlation between the polarization signal and the inclination, nor with the contribution of reflected photons throughout the Z-track.

accretion↗
DOE OSTI2026

Energy shift of Fe-K fluorescence lines due to low ionization demonstrated with XRISM in Centaurus X-3

The Fe K$\alpha$ fluorescence line at 6.4 keV is a powerful probe of cold matter surrounding X-ray sources and has been widely used in various astrophysical contexts. The X-ray microcalorimeter spectrometer onboard XRISM can measure line shifts with unprecedented precision of ${\sim }$0.2 eV, equivalent to a line-of-sight velocity of ${\sim }$10 km s$^{-1}$. At this level of accuracy, however, several factors that influence the line energy must be carefully considered prior to astrophysical interpretation. One such important factor is the ionization degree, Fe$^{q+}$. The K$\alpha$ line shifts redward by ${\sim }$4 eV as q increases from 0 (neutral) to 8 (Ar-like). Additionally, the accompanying Fe K$\beta$ line at 7.06 keV shifts blueward by ${\sim }$30 eV from $q=0$ to 8. We demonstrate that this effect is actually observable in the XRISM data of the high-mass X-ray binary Centaurus X-3 (Cen X-3). We advocate that the differential energy shift between the K$\alpha$ and K$\beta$ line provides a robust estimate of q by decoupling from other effects that shift the two lines in the same direction. We derived $q \sim 5$ (Sc-like) for the fluorescing matter by comparing the observation with atomic structure calculations of our own and in the literature. By accounting for the derived charge state and the corresponding shift in the rest-frame line energy, we made corrections for this effect and reached a consistent residual shift among the K$\alpha$, K$\beta$, and the optical measurement attributable to the systemic velocity of the system. Consequently, we obtained a new constraint on the location of the cold matter. This ionization effect needs to be assessed in all use cases of the Fe K$\alpha$ line shift beyond Cen X-3, and the proposed metric is generally applicable to all of them.

Astronomy and AstroPhysics↗
DOE OSTI2026

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Orbital and short-term variability of iron absorption

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe K band on second timescales. Using data from the Performance Verification phase in 2024 April, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser, and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve’s capability to probe plasma environments in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
NASA NTRSDate not supplied

High EneRgy Observatory for Imaging X-rays (HEROIX) Mission Concept

The High EneRgy Observatory for Imaging X-rays (HEROIX) is a next-generation broadband X-ray space telescope concept under development as a Medium-Class Explorer mission. Building on more than three decades of replicated NiCo shell technology at NASA’s Marshall Space Flight Center, HEROIX will deliver 5 arcsecond half-power-diameter imaging resolution in the hard X-ray band. The observatory will incorporate four co-aligned telescopes with advanced multilayer coatings, extending the bandpass into the hard X-ray regime and providing an integrated effective area of 380 cm² at 30 keV. A suite of focal plane detectors will enable imaging capability throughout its 1–80 keV bandpass. This combination of angular resolution and high-energy coverage will open new windows into non-thermal processes in extreme astrophysical environments. HEROIXwill deliver transformative science, advancing the astrophysics community’s understanding of supernova physics, the X-ray binary population of the Galactic Center, accretion-driven galaxy evolution, and the origin of the Cosmic X-ray Background.

Nick Thomas↗
DOE OSTI2026

4th TDAMM Workshop White Paper

Time-Domain and Multi-Messenger Astrophysics (TDAMM) is entering a fundamentally new phase characterized by an unprecedented increase in the rate and diversity of astrophysical transient detections. The community is transitioning from a discovery-limited to a follow-up-limited era, driven by major investments across electromagnetic, gravitational-wave, and neutrino observatories. Upcoming facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and wide-field survey instruments will produce a deluge of time-domain alerts, reaching millions of events per night. Simultaneously, upgrades to the gravitational-wave network (LVK O5 and beyond) and neutrino observatories (IceCube Gen2) will significantly increase the detection rates of non-electromagnetic messengers. New high-energy missions and expansions of the InterPlanetary Network (IPN) will further enhance discovery capabilities across the gamma-ray and X-ray regimes. This convergence of capabilities represents a transformative opportunity: for the first time, the community will routinely detect rare and high-impact events across multiple messengers. However, the scientific return from these discoveries will depend critically on the ability to rapidly identify, prioritize, and coordinate follow-up observations across a heterogeneous and globally distributed set of facilities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2026

X-Ray Polarization From the Atoll 4U 1735−44 Suggests a Low Inclination

X-ray polarimetry is a new tool capable of probing the geometry of accretion onto weakly magnetized neutron stars. Here we present the first X-ray spectro-polarimetric results from coordinated observations of the atoll source 4U 1735−44, conducted with the Imaging X-ray Polarimetry Explorer (IXPE), NICER, and NuSTAR . Over the 2–8 keV energy range, we obtained a marginal detection of polarization with the polarization degree of 1.4% ± 0.7% and polarization angle of −29◦ ± 14◦, corresponding to a 3σ upper limit on the polarization degree of 3.5%. The best-fit model to describe the spectrum comprises a thermal component associated with the accretion disk, a Comptonized black body component, and a relativistic reflection component. From the reflection model, we infer a disk inclination of ∼40◦. The spectroscopic and polarimetric properties of 4U 1735−44 are consistent with those observed in other atoll sources studied by IXPE, with its low polarization likely due to its low inclination

X-rays: binaries↗
DOE OSTI2026

High-spectral-resolution X-Ray Observations of the Evolved Supermassive Stellar Binary System η Carinae: Fe K α Band Profile Revealed with XRISM

The supermassive binary system, η Carinae, is experiencing enormous wind-driven mass loss at a rate unparalleled in the rest of the Galaxy. Their wind–wind collision (WWC) continuously produces shock heated, X-ray-emitting plasmas. The XRISM X-ray observatory observed the system in 2023 and 2024 when the X-ray emission began to increase toward periastron passage in 2025. This paper reports unprecedentedly high-resolution X-ray spectra in the Fe Kα band between 6.2 and 7.1 keV, obtained with the Resolve X-ray microcalorimeter. The hydrogen-like (Lyα) and helium-like (Heα) lines reveal three velocity components. Two of them are broadened with maximum velocities of 2000–3000 km s −1 , likely originating from the postshock companion wind. The other is relatively narrow, with a Gaussian broadening of only ∼290 km s −1 in 1σ, which may originate from the postshock companion wind at the WWC stagnation point or penetrating the primary wind. The Fe fluorescent lines exhibit a moderate blueshift and broadening with velocities at 100–200 km s −1 , consistent with the primary wind’s velocity field. The spectra also confirm a Compton shoulder of the Heα line complex for the first time. Both fluorescing and scattering spectral profiles indicate that the binary system is seen from the companion side during these observations. The flux ratio of the Compton-scattering emission to the fluorescent line suggests substantial hydrogen depletion of the primary wind, expected from CNO-cycled hydrogen nuclear fusion gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Harmonic suppression gratings for soft X-ray monochromators

We describe an approach to harmonic suppression in soft X-ray monochromators by engineering the reflection grating’s diffraction pattern to approximate a sinusoidal amplitude. At synchrotron and free-electron laser sources, X-ray beamlines powered by insertion devices produce a spectrum containing harmonic photon energies that can couple unwanted light into experiments. Beamlines in the soft X-ray energy range (100 eV to 2 keV) commonly employ energy-filtering elements to suppress these harmonics. Available approaches tend to be inefficient, significantly reducing the transmitted power. We show that with pseudo-grayscale binary halftone patterns, gratings can approximate a sinusoidal amplitude and suppress higher diffraction orders. Prototype demonstrations of lithographically fabricated gratings were conducted on a soft X-ray beamline with photon energies of 110 eV and 330 eV. Relative to a square-wave amplitude grating, the third-harmonic intensity was reduced by a factor of 9.0 with a first-order efficiency reduction of 38%.

Goldberg, Kenneth A. [Lawrence Berkeley National L↗
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Metal oxide-promoted calcium cuprate catalysts for diol oxidative dehydrocyclization to lactones

Here, this work investigates structure-function relationships in electronically tunable, redox-active, basic Cu-Ca mixed metal oxide catalysts for oxidative dehydrocyclization of liquid diols to lactones. Compositional screening identified Ni 2+ and Zn 2+ as effective promoters that increase the surface Cu 2+ population by ∼1.7× and Cu-normalized activity for liquid 1,4-butanediol conversion to γ-butyrolactone by ∼3–4×. In situ Raman spectroscopy, in situ X-ray absorption spectroscopy (XAS), in situ diffuse-reflectance Fourier transform infrared spectroscopy (DRIFTS), ex situ X-ray diffraction (XRD), and H 2 -temperature-programmed reduction (H 2 -TPR) show that Ni 2+ or Zn 2+ incorporation promotes the formation of Ca 0.82 Cu 1.00 O 2 nanoparticles under mild calcination conditions. This cuprate phase features stronger and shorter Cu–O bonds (1.90 Å) than inactive bulk CuO (1.95 Å) and square-planar Cu 2+ O 4 sites with enhanced d z2 electrophilicity, strengthening alkoxy adsorption. Pyridine-DRIFTS confirms the purely basic nature of the catalyst surface, while methanol-DRIFTS indicates Cu 2+ surface enrichment with Ni or Zn promotion, where Cu–O(Ca)–Cu sites can exist as amorphous domains or a truncation layer on crystalline nanoparticles.

09 BIOMASS FUELS↗
DOE OSTI2026

Exploring the connection between compact object mergers and fast X-ray transients

Context. The connection between compact object mergers and some extragalactic fast X-ray transients (FXRTs) has long been hypothesised but never ultimately established. Aims. In this work, we investigate two FXRTs, the LEIA X-ray Transient LXT 240402A and the Einstein Probe EP 250207b, whose precise positions lie close to nearby (z ≲ 0.1) quiescent galaxies with a negligible probability of chance coincidence, identifying them as particularly promising cases of merger-driven explosions in the local Universe. Methods. We used Chandra to derive accurate localisations for both events and secure otherwise ambiguous associations with their optical counterparts. Deep optical and near-infrared observations with VLT, GTC, and LBT were performed to characterise the surrounding environment and search for kilonova emission, the hallmark of neutron star mergers. Complementary early-time X-ray monitoring with Swift and Einstein Probe was used to constrain the non-thermal afterglow. Results. We find that both FXRTs remain compatible with a compact binary merger progenitor, which produced low-mass ejecta and kilonova emission subdominant to the afterglow. However, alternative explanations such as a distant (z ≳ 1) core-collapse supernova cannot be conclusively ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗
DOE OSTI2026

Connecting Current and Future Dual Active Galactic Nucleus Searches to LISA and Pulsar Timing Array Gravitational-wave Detections

Abstract Dual active galactic nuclei (DAGN) mark the observable stage of massive black hole (MBH) pairing during galaxy mergers and are the progenitors of the MBH binaries that generate low-frequency gravitational waves. Using the large-volume ASTRID cosmological simulation, we construct mock DAGN catalogs tailored to the selection functions of current (COSMOS-Web and DESI) and forthcoming (Nancy Grace Roman Space Telescope (Roman) and Lynx X-ray Observatory (Lynx)) surveys, enabling direct comparisons between simulations and observations. With realistic observational selections, ASTRID reproduces the observed dual fractions, projected separations, and host-galaxy properties across redshifts. We predict a substantial population of small-separation (<5 kpc) duals that remain inaccessible to current surveys, demonstrating that the apparent paucity of subkiloparsec systems in COSMOS-Web is primarily a consequence of observational selection rather than an intrinsic absence. Following each simulated dual to coalescence, we show that DAGN are effective tracers of MBH mergers: ∼30%–70% merge within ≲1 Gyr, and 20%–60% of these mergers produce gravitational-wave signals detectable by the Laser Interferometer Space Antenna (LISA). Duals observable with Roman and Lynx are the progenitors of ∼10%–50% of low-redshift LISA sources and contribute ∼30% of the PTA-band stochastic gravitational-wave background. We further identify massive green-valley galaxies hosting moderate-luminosity active galactic nuclei (AGN), together with massive star-forming galaxies containing bright quasars at z > 1, as the environments most likely to host imminent MBH binaries. These results establish a unified cosmological framework connecting DAGN demographics, MBH binary evolution, and gravitational-wave sources, while identifying high-priority targets for coordinated electromagnetic and multimessenger observations in the coming decade.

Chen, Nianyi [Max-Planck-Institut für Astrophysik;↗
DOE OSTI2026

Conventional superconductivity in single-crystalline BiPt

Binary Bi-Pd/Pt systems have attracted a lot of interest because of their topologically nontrivial nature along with superconductivity. We report the structural and superconducting properties of high-quality single-crystalline BiPt using a comprehensive range of experimental techniques, including X-ray diffraction, electron microscopy, muon spin rotation/relaxation (𝜇⁢SR), magnetization, resistivity, and heat capacity. Our findings establish that BiPt is a weak type-II superconductor with a transition temperature (𝑇 𝑐 ) of 1.2 K which exhibits pronounced anisotropic superconducting characteristics attributed to its hexagonal crystal structure. Magnetization and electronic transport studies reveal that BiPt lies within the dirty limit, while 𝜇⁢SR and heat capacity data indicate conventional 𝑠-wave superconductivity that maintains time-reversal symmetry. Here, this work provides valuable insights into the pairing symmetry and superconducting mechanism of topologically trivial BiPt, a sound comparison system for other Bi-based topologically nontrivial superconductors.

Demagnetization↗
DOE OSTI2026

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗
DOE OSTI2026

A high transmission tender X-ray monochromator employing a matched pair of multilayer grating and mirror

It is now well established that multilayer coated gratings may offer high diffraction efficiencies over the tender X-ray range, from about 1 keV to 5 keV, covering the gap between single layer coated grating monochromators and crystal monochromators. Nevertheless, few beamlines in the world are using such gratings in their monochromator. The successful implementation of a multilayer grating monochromator requires producing a matched pair of a multilayer grating and a multilayer mirror, and this matching is not straightforward as it must account for different incidence angles and refraction effects on each element. Here we review the realization of the multilayer grating monochromator of the SIRIUS beamline which has been successfully in service for several years. We show how, by alternating computer simulation with our diffraction code and measurements, we could optimize the monochromator transmission on a very wide energy range. After the grating was coated, it was found that the angle of optimal efficiency versus photon energy was significantly different from what was predicted by a simple conformal model of binary layers. Layer interdiffusion and profile smoothing during the deposition process must be added to the multilayer model to reproduce the measured data. The critical adjustment of the mirror multilayer period is achieved by the lateral translation of the mirror, which was given a small transverse period gradient. The monochromator is thus providing high transmission efficiency in the 1 to 5 keV energy range, more than 30% over 2.5 keV and up to 46% at 4.6 keV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Interfacial Chemistry Involved in Selective Separation of NMC/LMO and LCO/LMO Binary Cathode Materials by Froth Flotation Using Oleic Acid

The variability in cathode compositions within recycled lithium-ion battery (LIB) feedstocks poses a significant challenge to efficient downstream refining processes. This study demonstrates the feasibility of using froth flotation with oleic acid as a collector to selectively separate lithium nickel-manganese-cobalt oxide (NMC) and lithium cobalt oxide (LCO) from lithium manganese oxide (LMO) materials. Laboratory-scale flotation tests achieved an 80% separation efficiency in a single stage, producing a froth product with >90% purity of NMC/LCO at approximately 90% yield. Concurrently, the LMO materials were enriched in the sink product with ∼90% purity and ∼90% yield. This approach was further validated using recycled cathode materials, confirming its applicability to realistic feedstocks. The underlying mechanism governing the selective separation of NMC/LCO from LMO was investigated using ζ-potential measurements, contact angle measurements, bubble-particle attachment experiments, and X-ray photoelectron spectroscopy (XPS) analysis. Both contact angle and bubble-particle attachment results confirmed that oleic acid adsorption rendered NMC and LCO surfaces hydrophobic, thereby enhancing flotation recovery. At pH 5, oleic acid adsorbed preferentially onto NMC and LCO surfaces via electrostatic interactions, while exhibiting minimal adsorption on LMO surfaces. However, separation efficiency deteriorated at higher pH, which was attributed to the co-flotation of LMO materials caused by oleate chemisorption on MnOH + species. This work establishes froth flotation as a viable cathode/cathode separation strategy, providing a low-cost, scalable pathway to preconcentrate and enrich nickel-rich and cobalt-rich cathode active materials from incompatible cathode chemistries for direct recycling or hydrometallurgical processing. Furthermore, this study reveals, for the first time, the mechanism of oleate adsorption on the surface of different cathode materials.

PH↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗