Search NASA⌕ Search

SEARCH · Search NASA

Results for “matric potential”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bottomonium properties in the quark-gluon plasma from a lattice-QCD informed 𝑇-matrix approach

Recent computations of bottomonium correlation functions with extended sources in lattice-discretized quantum chromodynamics (lQCD) provide new insights into heavy-quark dynamics at distance scales which are of the order of the inverse temperature. We analyze these results employing the thermodynamic 𝑇-matrix framework, in a continued effort to interpret lQCD data for quarkonium correlation functions in a nonperturbative and self-consistently solved quantum-many-body approach to a strongly coupled quark-gluon plasma (QGP). Its key inputs are the in-medium driving kernel (potential) of the scattering equation and an interference function which implements 3-body effects in the quarkonium coupling to the thermal medium. A simultaneous description of lQCD results for the bottomonium correlators with extended operators and the previously analyzed Wilson line correlators only requires minor refinements of the potential but calls for stronger interference effects at larger separation of the bottom quark and antiquark. We then analyze the poles of the self-consistent 𝑇 matrices on the real axis to assess the survival of the various bound states. Here, we estimate the pertinent temperatures where the poles disappear for the various bottomonium states and discuss the relation to the corresponding peaks in the bottomonium spectral functions. We also recalculate the spatial diffusion coefficient of the QGP and find it to be similar to that in our previous study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Spectroscopic Studies of Charge-Transfer Character and Photoresponses of F 4 TCNQ-Based Donor–Acceptor Complexes

F 4 TCNQ (2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane) is used widely as a hole-doping agent in photoresponsive organic semiconducting materials, yet relatively little is known about the photoresponses of the F 4 TCNQ· - anion generated via doping. Furthermore, there is still relatively little systematic exploration of how the properties of the local material or chemical environment impacts the driving force for generating these charge-transfer complexes. Here we present spectroscopic and photophysical studies of F 4 TCNQ in charge-transfer complexes (CTCs) with the electron donor N,N’-Diphenyl-N-N’-di-p-tolylbenzene-1,4-diamine (MPDA) both in dichloroethane solution and polystyrene matrices. Integer charge transfer (ICT) between donor and acceptor occurs readily in dichloroethane solvent to form F 4 TCNQ· - :MPDA + CTCs, due to a ~150 mV difference in MPDA + /MPDA and F 4 TCNQ/F 4 TCNQ· - reduction potentials. Ultrafast spectroscopic studies of the CTC as well as electrochemically generated F 4 TCNQ· - and MDPA + reveal that the photoresponses of these CTCs are dominated by that of the dopant anion, including rapid deactivation (800 fs) after excitation to the anion D 1 excited state, followed by slower (~10 ps) vibrational cooling in the anion D 0 state. Excitation to the higher-lying D 2 state results in a rapid relaxation to the D 1 state, in contrast to direct D 2 -D 0 relaxation previously observed for F 4 TCNQ· - in gas phase. CTCs embedded in polystyrene (PS) matrices are observed to lose their integer charge-transfer character upon evaporation of solvent, as evidenced by changes to electronic and vibrational absorption features associated with F 4 TCNQ· - . This change is attributed to the loss of solvent stabilization of the ion pair formed through the charge-transfer reaction. Ultrafast spectral measurements reveal that the photoresponses of the partial charge-transfer (PCT) species embedded in PS are still highly similar to those of the ICT species and unlike that of neutral F 4 TCNQ, implying the electronic properties of the PCT state are likewise dominated by properties of the reduced acceptor molecule. We conclude that excitation of ICT or PCT states introduces optical losses for photoresponses of doped organic semiconductor materials due to the large anion absorption cross section and its rapid, dissipative deactivation dynamics.

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

DOME: Directional medical embedding vectors from Electronic Health Records

Motivation: The increasing availability of Electronic Health Record (EHR) systems has created enormous potential for translational research. Recent developments in representation learning techniques have led to effective large-scale representations of EHR concepts along with knowledge graphs that empower downstream EHR studies. However, most existing methods require training with patient-level data, limiting their abilities to expand the training with multi-institutional EHR data. On the other hand, scalable approaches that only require summary-level data do not incorporate temporal dependencies between concepts. Methods: We introduce a DirectiOnal Medical Embedding (DOME) algorithm to encode temporally directional relationships between medical concepts, using summary-level EHR data. Specifically, DOME first aggregates patient-level EHR data into an asymmetric co-occurrence matrix. Then it computes two Positive Pointwise Mutual Information (PPMI) matrices to correspondingly encode the pairwise prior and posterior dependencies between medical concepts. Following that, a joint matrix factorization is performed on the two PPMI matrices, which results in three vectors for each concept: a semantic embedding and two directional context embeddings. They collectively provide a comprehensive depiction of the temporal relationship between EHR concepts. Results: We highlight the advantages and translational potential of DOME through three sets of validation studies. First, DOME consistently improves existing direction-agnostic embedding vectors for disease risk prediction in several diseases, for example achieving a relative gain of 5.5% in the area under the receiver operating characteristic (AUROC) for lung cancer. Second, DOME excels in directional drug-disease relationship inference by successfully differentiating between drug side effects and indications, correspondingly achieving relative AUROC gain over the state-of-the-art methods by 10.8% and 6.6%. Finally, DOME effectively constructs directional knowledge graphs, which distinguish disease risk factors from comorbidities, thereby revealing disease progression trajectories. The source codes are provided at https://github.com/celehs/Directional-EHRembedding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Perturbation theory approach to predict the covariance matrices of the galaxy power spectrum and bispectrum in redshift space

In this paper, we predict the covariance matrices of both the power spectrum and the bispectrum, including full non-Gaussian contributions, redshift space distortions, linear bias effects, and shot-noise corrections, using perturbation theory (PT). To quantify the redshift-space distortion effect, we focus mainly on the monopole and quadrupole components of both the power and bispectra. We, for the first time, compute the 5- and 6-point spectra to predict the cross-covariance between the power and bispectra, and the autocovariance of the bispectrum in redshift space. Further, we test the validity of our calculations by comparing them with the covariance matrices measured from the MultiDark-Patchy mock catalogues that are designed to reproduce the galaxy clustering measured from the Baryon Oscillation Spectroscopic Survey Data Release 12. We argue that the simple, leading-order PT works because the shot-noise corrections for the Patchy mocks are more dominant than other higher order terms we ignore. In the meantime, we confirm some discrepancies in the comparison, especially of the cross-covariance. We discuss potential sources of such discrepancies. We also show that our PT model reproduces well the cumulative signal-to-noise ratio of the power spectrum and the bispectrum as a function of maximum wavenumber, implying that our PT model captures successfully essential contributions to the covariance matrices.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Role of Mineralogy in Controlling Fracture Formation

The presence of fractures in caprocks can pose increased risks in subsurface energy systems and processes like CO 2 sequestration by introducing high-permeability leakage paths. Fracture apertures and permeability can be altered through mineral dissolution and precipitation reactions, but the reactive evolution of fractures is not well understood. In fractures, minerals that are otherwise inaccessible to reactive fluids can become exposed, resulting in mineral reactions unpredicted by bulk formation data. This work seeks to understand the relationship between mineralogy and fracture formation to enhance our understanding of reactive fracture evolution and CO 2 leakage potential. Here, the mineral compositions of mechanically induced fracture surfaces in samples of the Mancos and Marcellus shales have been quantified and compared to those of the near-fracture matrices using imaging and bulk X-ray diffraction (XRD) data. In the Mancos shale, the concentrations of clay minerals are enhanced along fracture surfaces with respect to the bulk, and the fracture is most likely to form at kaolinite–kaolinite interfaces. Further, evaluation of the mineralogical spatial variability through cross-correlation analysis of the surrounding matrix in images of samples cut perpendicular to the fracture shows that clay is 16.7 times more likely to be present than carbonate minerals near the fracture surface. The high correlation persists roughly 200 μm into the surrounding matrix for the Mancos sample and implies that the fracture formed within a defined clay-rich lithofacies.

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

Microbial potentiometric sensor technology for real-time detecting and monitoring of toxic metals in aquatic matrices

Considering that toxic metals can affect metabolic processes in microorganisms adversely, it can be hypothesized that these metals in water matrices would induce a decrease in metabolic activity of the biofilm microorganisms populating the surface of a sensing electrode, which could be registered as a change in the open-circuit potential (OCP) generated by the biofilm microorganisms. The goal of this study was to test this hypothesis and demonstrate the underlying principle that microbial potentiometric sensor (MPS) technology could be used for long-term and real-time monitoring and detection of rapid changes in metal concentrations in realistic aquatic environments. To address the goal, four objective were addressed: (1) a batch reactor with three graphite-based MPS electrodes was fabricated; (2) a set of single-ion solutions and one multiple ion solution were prepared reflecting realistic concentrations of metals found in electroplating wastewaters; (3) the responses of the MPS to the simultaneous presence of multiple toxic metal ions in a single solution were measured; and (4) the changes of the MPS signals to the presence of individual metal ion solutions were examined. While the hypothesis was validated, the study also revealed that the MPS was sufficiently sensitive to not only detect, but also quantify, toxic metal ion concentrations in aqueous solutions. The coefficients of determination, which were R 2 > 0.995, and responsiveness of <1 μmol/L for some toxic metal cations, strongly support the performance of MPS technology in the echelons of expensive analytical tools capable detecting and measuring trace elements.The magnitude of the MPS response was toxic metal specific. When the molar concertation normalizes the inhibition portion of the signal area, the assessed sensitivity order was: Se>Cd>Pb>Ag>Ni> Zn. The study provides valuable information for enforcement agents, environmental professionals, and wastewater treatment operators, so toxic metal pollution and its detrimental impacts can be prevented and mitigated.

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

Implications of Microstructure in Helium-Implanted Nanocrystalline Metals

Helium bubbles are known to form in nuclear reactor structural components when displacement damage occurs in conjunction with helium exposure and/or transmutation. If left unchecked, bubble production can cause swelling, blistering, and embrittlement, all of which substantially degrade materials and—moreover—diminish mechanical properties. On the mission to produce more robust materials, nanocrystalline (NC) metals show great potential and are postulated to exhibit superior radiation resistance due to their high defect and particle sink densities; however, much is still unknown about the mechanisms of defect evolution in these systems under extreme conditions. Here, the performances of NC nickel (Ni) and iron (Fe) are investigated under helium bombardment via transmission electron microscopy (TEM). Bubble density statistics are measured as a function of grain size in specimens implanted under similar conditions. While the overall trends revealed an increase in bubble density up to saturation in both samples, bubble density in Fe was over 300% greater than in Ni. To interrogate the kinetics of helium diffusion and trapping, a rate theory model is developed that substantiates that helium is more readily captured within grains in helium-vacancy complexes in NC Fe, whereas helium is more prone to traversing the grain matrices and migrating to GBs in NC Ni. Our results suggest that (1) grain boundaries can affect bubble swelling in grain matrices significantly and can have a dominant effect over crystal structure, and (2) an NC-Ni-based material can yield superior resistance to irradiation-induced bubble growth compared to an NC-Fe-based material and exhibits high potential for use in extreme environments where swelling due to He bubble formation is of significant concern.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Green synthesis of sulfur‐doped g‐C 3 N 4 nanosheets for enhanced removal of oxytetracycline under visible‐light irradiation and reduction of its N ‐nitrosodimethylamine formation potential

Abstract BACKGROUND In this study, sulfur‐doped g‐C 3 N 4 (SCN) nanosheets were synthesized using thiourea as a precursor, and then used to catalyze the degradation of oxytetracycline (OTC) under visible‐light irradiation. OTC was chosen as target contaminant because of it being an important precursor of N ‐nitrosodimethylamine (NDMA). RESULTS The removal of OTC was influenced by its initial concentration, initial pH, SCN dosage and water matrices. OTC (≤10 mg L −1 ) was almost completely removed in 40 min under the conditions of an initial pH of 7.0 and a SCN dose of 1.0 g L −1 , and its degradation was well described by the pseudo‐first‐order kinetics model. Water matrices would also influence the degradation of OTC in the system. Electron paramagnetic resonance analysis and trapping experiments confirmed that the function of the radicals for OTC removal followed the order of • O 2 − > • OH > h + . Eight transformation products were determined during the degradation of OTC, it being removed through four pathways: methyl oxidization, demethylation, decarbonylation and secondary alcohol oxidation. SCN was stable in the system, as the removal efficiency of OTC was not significantly reduced even after three recycles. In addition, 55.5% of NDMA formation potential of OTC was reduced after SCN catalytic oxidation. CONCLUSIONS The results imply that the visible‐light‐driven SCN catalytic oxidation is an effective technology not only for the removal of OTC but also for the reduction of NDMA formation potential. © 2021 Society of Chemical Industry

Dou, Yicheng↗
DOE OSTI2023

Search for low mass dark matter in DarkSide-50: the bayesian network approach

Abstract We present a novel approach for the search of dark matter in the DarkSide-50 experiment, relying on Bayesian Networks. This method incorporates the detector response model into the likelihood function, explicitly maintaining the connection with the quantity of interest. No assumptions about the linearity of the problem or the shape of the probability distribution functions are required, and there is no need to morph signal and background spectra as a function of nuisance parameters. By expressing the problem in terms of Bayesian Networks, we have developed an inference algorithm based on a Markov Chain Monte Carlo to calculate the posterior probability. A clever description of the detector response model in terms of parametric matrices allows us to study the impact of systematic variations of any parameter on the final results. Our approach not only provides the desired information on the parameter of interest, but also potential constraints on the response model. Our results are consistent with recent published analyses and further refine the parameters of the detector response model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Control of iodine and tellurium for use in domestic molybdenum-99 production efforts

Three initial proof-of-concept experiments were performed to test the removal of I and Te from solutions containing Mo for potential abatement during chemical processing to purify 99 Mo. All three proof-of-concept experiments produced results indicating that further optimization could provide methods for I and/or Te abatement from multiple matrices. Decamp and Happel (2013) describe the use of CL resin from Eichrom Technologies loaded with Ag + combined with Amberlite XAD-4 in the chromatographic separation of irradiated fission products in 1 M HNO 3 to separate iodine species where iodide and iodate react with Ag + and I2 adheres to the Amberlite resin. The method worked well under processing conditions in a medical isotope facility although many separation methods use higher concentrations of HNO 3 which could adversely impact the Ag + loading. The experiment performed here used a roughly 2 mL column of CL resin loaded with Ag + and Amberlite XAD-4 to separate µg quantities of I, Te, and Mo in 3 M HNO 3 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Quantum simulation of quantum phase transitions using the convex geometry of reduced density matrices

Transitions of many-particle quantum systems between distinct phases at absolute-zero temperature, known as quantum phase transitions, require an exacting treatment of particle correlations. Here in this work, we present a general quantum-computing approach to quantum phase transitions that exploits the geometric structure of reduced density matrices. While typical approaches to quantum phase transitions examine discontinuities in the order parameters, the origin of phase transitions—their order parameters and symmetry breaking—can be understood geometrically in terms of the set of two-particle reduced density matrices (2-RDMs). The convex set of 2-RDMs provides a comprehensive map of the quantum system including its distinct phases as well as the transitions connecting these phases. Because 2-RDMs can potentially be computed on quantum computers at nonexponential cost, even when the quantum system is strongly correlated, they are ideally suited for a quantum-computing approach to quantum phase transitions. We compute the convex set of 2-RDMs for a Lipkin-Meshkov-Glick spin model on IBM superconducting-qubit quantum processors. Even though computations are limited to few-particle models due to device noise, comparisons with a classically solvable 1000-particle model reveal that the finite-particle quantum solutions capture the key features of the phase transitions including the strong correlation and the symmetry breaking.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Competing ionization and dissociation: Extension of the energy-dependent frame transformation to the gerade symmetry of H 2

This article solves two major tasks that frequently arise in the theory of electron collisions with a target molecular cation. First, it extends the energy-dependent frame transformation (EDFT) treatment, which is needed to map fixed-nuclei electron-molecule scattering matrices into an energy-dependent laboratory-frame scattering matrix with vibrational channel indices. The EDFT mapping can now be carried out even when the target molecule possesses multiple low-energy potential curves, significantly transcending previous applications. Second, it implements a method to extract the rest of the full laboratory-frame scattering matrix, i.e., the columns and rows describing input and/or output dissociation channels. The treatment is benchmarked in this article against the essentially exact solution of a refined two-dimensional model of the singlet gerade Σ symmetry of H 2 . Our tests demonstrate that the theory accurately maps fixed-nuclei scattering information, of the type provided by existing electron-molecule computer codes, into a laboratory-frame scattering matrix that includes both ionization and dissociation. Furthermore, this treatment can provide a general framework applicable to a broad class of electron collision processes involving diatomic target ions, suitable for an accurate description of challenging processes such as dissociative recombination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Actinium purification from elemental components of stainless steel and aluminum

Radioisotopes of actinium are valuable because of their potential use in the medical industry. Actinium-225 shows promise for treating disease—like cancer—via targeted alpha therapy (TAT), and the longer-lived actinium-227 is the parent of two radionuclides with potential application in TAT radiopharmaceuticals, thorium-227 and radium-223. Continuing progress in the development of these medical applications requires robust and diverse methods for extracting and purifying actinium from a wide range of matrices. Further, to define the strengths and limitations of separation methods commonly employed in actinium processing, we characterized the performance of ion exchange and extraction chromatographic steps for removing contaminants found in stainless steel (chromium, iron, nickel, and silicon) and aluminum from actinium-227. While AG-MP1 anion exchange resin with HCl (aq) successfully removed iron, cation exchanger AG 50W-X8 removed most other contaminants. The most persistent contaminant was aluminum, which was removed using a DGA Normal extraction resin. These results are presented within the context of applying that methodology toward actinium purification strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Layered Heterostructure Ionogel Electrolytes for High‐Performance Solid‐State Lithium‐Ion Batteries

Abstract Ionogel electrolytes based on ionic liquids and gelling matrices offer several advantages for solid‐state lithium‐ion batteries, including nonflammability, wide processing compatibility, and favorable electrochemical and thermal properties. However, the absence of ionic liquids that are concurrently stable at low and high potentials constrains the electrochemical windows of ionogel electrolytes and thus their high‐energy‐density applications. Here, ionogel electrolytes with a layered heterostructure are introduced, combining high‐potential (anodic stability: >5 V vs Li/Li + ) and low‐potential (cathodic stability: <0 V vs Li/Li + ) imidazolium ionic liquids in a hexagonal boron nitride nanoplatelet matrix. These layered heterostructure ionogel electrolytes lead to extended electrochemical windows, while preserving high ionic conductivity (>1 mS cm −1 at room temperature). Using the layered heterostructure ionogel electrolytes, full‐cell solid‐state lithium‐ion batteries with a nickel manganese cobalt oxide cathode and a graphite anode are demonstrated, exhibiting voltages that are unachievable with either the high‐potential or low‐potential ionic liquid alone. Compared to ionogel electrolytes based on mixed ionic liquids, the layered heterostructure ionogel electrolytes enable higher stability operation of full‐cell lithium‐ion batteries, resulting in significantly enhanced cycling performance.

Hyun, Woo Jin↗
DOE OSTI2021

CASTELO: clustered atom subtypes aided lead optimization—a combined machine learning and molecular modeling method

Background: Drug discovery is a multi-stage process that comprises two costly major steps: pre-clinical research and clinical trials. Among its stages, lead optimization easily consumes more than half of the pre-clinical budget. We propose a combined machine learning and molecular modeling approach that partially automates lead optimization workflow in silico, providing suggestions for modification hot spots. Results: The initial data collection is achieved with physics-based molecular dynamics simulation. Contact matrices are calculated as the preliminary features extracted from the simulations. To take advantage of the temporal information from the simulations, we enhanced contact matrices data with temporal dynamism representation, which are then modeled with unsupervised convolutional variational autoencoder (CVAE). Finally, conventional and CVAE-based clustering methods are compared with metrics to rank the submolecular structures and propose potential candidates for lead optimization. Conclusion: With no need for extensive structure-activity data, our method provides new hints for drug modification hotspots which can be used to improve drug potency and reduce the lead optimization time. It can potentially become a valuable tool for medicinal chemists.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An extended numerical manifold method for unsaturated soil-water interaction analysis at micro-scale

To investigate unsaturated soil-water interaction at micro-scale, this work extends the numerical manifold method (NMM) by incorporating a soil-water coupling model considering specific capillary water distribution and capillary force calculation. The soil skeleton is constructed by a soil skeleton generation algorithm with random polygons. To more realistically capture the interaction between soil grains and capillary water, a capillary mechanics-based geometric algorithm is proposed to iteratively calculate the capillary water distribution. The capillary forces corresponding to the capillary water distribution are calculated based on the Young-Laplace equation. The proposed capillary water solving framework is first verified by reproducing the soil-water characteristic curve and the capillary water distribution of an ideal contact-disk model against analytical solutions. To further validate the ability of the capillary water solving framework to predict hydraulic behavior of the real soil, a laboratory test on the Toyoura sand is reproduced numerically. Then an ideal direct shear test is performed to further validate the two-way soil-water coupling procedure, in which a comparison between the numerical and analytical results regarding the shear strength and matric suction is presented. Finally, microscopic hydraulic and compression tests are conducted on two soil specimens with the same porosity and mean grain diameter but different uniformity coefficients. The results elucidate that the extended method is a potential tool to explore unsaturated soil behaviors at micro-scale.

58 GEOSCIENCES↗
DOE OSTI2023

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

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

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗