Search NASA⌕ Search

SEARCH · Search NASA

Results for “closed systems”

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 55 records · Page 3

Many-Body Physics in the NISQ Era: Quantum Programming a Discrete Time Crystal

Recent progress in the realm of noisy intermediate-scale quantum (NISQ) devices represents an exciting opportunity for many-body physics by introducing new laboratory platforms with unprecedented control and measurement capabilities. We explore the implications of NISQ platforms for many-body physics in a practical sense: we ask which physical phenomena, in the domain of quantum statistical mechanics, they may realize more readily than traditional experimental platforms. While a universal quantum computer can simulate any system, the eponymous noise inherent to NISQ devices practically favors certain simulation tasks over others in the near term. As a particularly well-suited target, we identify discrete time crystals (DTCs), novel nonequilibrium states of matter that break time translation symmetry. These can only be realized in the intrinsically out-of-equilibrium setting of periodically driven quantum systems stabilized by disorder-induced many-body localization. While promising precursors of the DTC have been observed across a variety of experimental platforms—ranging from trapped ions to nitrogen-vacancy centers to NMR crystals—none have all the necessary ingredients for realizing a fully fledged incarnation of this phase, and for detecting its signature long-range spatiotemporal order. We show that a new generation of quantum simulators can be programmed to realize the DTC phase and to experimentally detect its dynamical properties, a task requiring extensive capabilities for programmability, initialization, and readout. Specifically, the architecture of Google’s Sycamore processor is a remarkably close match for the task at hand. We also discuss the effects of environmental decoherence, and how they can be distinguished from ‘internal’ decoherence coming from closed-system thermalization dynamics. Already with existing technology and noise levels, we find that DTC spatiotemporal order would be observable over hundreds of periods, with parametric improvements to come as the hardware advances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Competing antiferromagnetic-ferromagnetic states in a d 7 Kitaev honeycomb magnet

The Kitaev model is a rare example of an analytically solvable and physically instantiable Hamiltonian yielding a topological quantum spin liquid ground state. Here we report signatures of Kitaev spin liquid physics in the honeycomb magnet Li 3 Co 2 SbO 6 , built of high-spin d 7 (Co 2+ ) ions, in contrast to the more typical low-spin d 5 electron configurations in the presence of large spin-orbit coupling. Neutron powder diffraction measurements, heat capacity, and magnetization studies support the development of a long-range antiferromagnetic order space group of C C 2/m, below T N = 11K at μ 0 H = 0T. The magnetic entropy recovered between T = 2 and 50 K is estimated to be 0.6Rln2, in good agreement with the value expected for systems close to a Kitaev quantum spin liquid state. The temperature-dependent magnetic order parameter demonstrates a β value of 0.19(3), consistent with XY anisotropy and in-plane ordering, with Ising-like interactions between layers. Further, we observe a spin-flop-driven crossover to ferromagnetic order with space group of C2/m under an applied magnetic field of μ 0 H ≈ 0.7T at T = 2K. Magnetic structure analysis demonstrates these magnetic states are competing at finite applied magnetic fields even below the spin-flop transition. Both the d 7 compass model, a quantitative comparison of the specific heat of Li 3 Co 2 SbO 6 , and related honeycomb cobaltates to the anisotropic Kitaev model further support proximity to a Kitaev spin liquid state. Furthermore, this material demonstrates the rich playground of high-spin d 7 systems for spin liquid candidates and complements known d 5 Ir- and Ru-based materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluating the state-of-the-art in remote volcanic eruption characterization Part I: Raikoke volcano, Kuril Islands

Raikoke, a small, unmonitored volcano in the Kuril Islands, erupted in June 2019. We integrate data from satellites (including Sentinel-2, TROPOMI, MODIS, Himawari-8), the International Monitoring System (IMS) infrasound network, and global lightning detection network (GLD360) with information from local authorities and social media to retrospectively characterize the eruptive sequence and improve understanding of the pre-, syn- and post- eruptive behavior. In this work, we observe six infrasound pulses beginning on 21 June at 17:49:55 UTC as well as the main Plinian phase on 21 June at 22:29 UTC. Each pulse is tracked in space and time using lightning and satellite imagery as the plumes drift eastward. Post-eruption visible satellite imagery shows expansion of the island's surface area, an increase in crater size, and a possibly-linked algal bloom south of the island. We use thermal satellite imagery and plume modeling to estimate plume height at 10–12 km asl and 1.5–2 × 10 6 kg/s mass eruption rate. Remote infrasound data provide insight into syn-eruptive changes in eruption intensity. Our analysis illustrates the value of interdisciplinary analyses of remote data to illuminate eruptive processes. However, our inability to identify deformation, pre-eruptive outgassing, and thermal signals, which may reflect the relatively short duration (~12 h) of the eruption and minimal land area around the volcano and/or the character of closed-system eruptions, highlights current limitations in the application of remote sensing for eruption detection and characterization.

58 GEOSCIENCES↗

The fate of carbon utilized by the subterranean termite Reticulitermes flavipes

Abstract Subterranean termites are ecosystem engineers that modulate the flow of carbon from dead wood to the atmosphere and soil, yet their contributions to the latter pool are largely unaccounted for in carbon cycling models. The fate of C from wood utilized by Reticulitermes flavipes (Kollar) was determined using a reductionist design in a closed system with δ 13 C labeled wood as a stable isotope tracer. The percentage of wood‐based carbon in termite respiratory gases, tissues, and organic deposits (frass and construction materials) was measured for five colonies to budget wood‐C mass distributed into metabolic and behavioral pathways during a 160‐h incubation period. We found that termites emitted 42% of the C from wood as gas (largely as carbon dioxide), returned 40% to the environment as organic deposits (frass and construction materials), and retained 18% in their tissues (whole alimentary tracts and de‐gutted bodies). Our findings affirm that termites are a source of greenhouse gases but are also ecosystem engineers that return approximately half the C from dead wood as organic deposits into their surrounding environment.

54 ENVIRONMENTAL SCIENCES↗

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimized V1G and V2G Electric Vehicle Fleet Management and Grid Transaction at Marine Corps Air Station Miramar in San Diego, CA

The overall technical goal of the project was to demonstrate an all-electric bi-directional non-tactical fleet at Marine Corps Air Station (MCAS) Miramar that was integrated and controlled with other distributed energy resources (DERs) (i.e., PV, stationary battery, and building loads) to provide resilience to critical electric loads in the event of grid outages, to minimize charging costs, and to provide economic energy resources to electricity markets. In this project, the specific, technical objectives were: 1. Demonstrate that bi-directional electric vehicles can provide critical complementary services to fixed storage batteries in microgrid applications while performing function as non-tactical vehicles. 2. Demonstrate participation of bi-directional (V2G) and unidirectional (V1G) PEVs for demand management and minimization of charging costs. 3. Demonstrate integration of multiple DERs for grid service participation. US Marine Corps Air Station (MCAS) Miramar in San Diego was the site of this electric vehicle-to-microgrid-utility grid test and demonstration project. Existing microgrid assets in this study included (1) a public works building; (2) a 30-kW rooftop photovoltaic (PV) system and (3) a separate 250 kW carport PV system. In this project, six bi-directional V2G vans were located at the MCAS Miramar’s showcase building-scale microgrid to develop and test technical capabilities that V2G can provide in microgrid applications (e.g., cost reduction and resiliency). These resources provided aggregated demand management and simulated participation in current retail DR programs. The vehicles used in this demonstration were selected because they provided functionality that MCAS Miramar needed, 15 passenger transport and facilities work cargo carrying capacity, and bi-directional charging capability that the research project required. All vehicles in this study were manufactured and distributed by VIA Motors, Inc. There were six vehicles total and each was VIA’s VTRUX eREV V2G model, a modified General Motors Chevrolet 2500 2WD van. Three of the vans were configured as passenger vans and the other three were configured as cargo vans. Each van had an on-board bi-direcrtional inverter/charger, Bel Power Solutions model 350INVCHGT150-120-240-8G nominally rated at +/-15 kW. The VIA van’s charging connector follows the J1772 charging protocol. The bi-directional EVSEs demonstrated in this study were manufactured by Coritech, Inc. Each VGI-80-AC charging station enabled enhanced V2G charging capability to a Clipper Creek CS-100 charging module. The enhanced capabilities included ethernet communication following the SEP2.0 protocol with a distributed energy resource function set and an operator screen displaying real-time SOC, voltage, and current. The VGI-80-AC charging stations are classified as level 2 with a maximum current output of 80 A or effectively 19 kW. The VIA van’s onboard charger limited the charging and discharging power to 15 kW in each direction. A control computer was installed in the EWOC and connected to an existing monitor. The V2G control communication network was a completely stand-alone closed system that did not have any connection to any other networks on the base. A cybersecure remote communication connection was created with a cellular modem, firewall hardware, and a virtual private network configuration.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mechanochemistry of Phosphorus and Arsenic Alloys for Visible and Infrared Photonics

Incorporating few‐layer 2D phosphorus–arsenic alloys (PAs) into optoelectronic devices requires a synthesis technique that allows control of the alloy composition while producing volumes of material suitable for application development such as photodetectors, solar cells, and lasers. With that goal in mind, high‐energy ball milling allows production of both orthorhombic (o‐PAs) and trigonal (t‐PAs) alloys by reacting red phosphorus and metallic arsenic powders. The synthesis follows a two‐step process in which arsenic rapidly reacts with red phosphorus to first produce t‐PAs followed by a slower phase transformation into o‐PAs; synthesis time and overall conversion rate are slightly enhanced by the presence of arsenic. Optical measurements on exfoliated alloys at the few‐layer atomic limit of the 2D PAs reveal emission spanning from the visible (≈1.9 eV) into the near‐infrared region, covering a broad application space. Rapid powder synthesis within a closed system for stochiometric control of the solid‐solution PAs alloys combined with solution‐based exfoliation opens up opportunities for a whole new class of optoelectronic devices based on PA nanomaterials.

Pedersen, Samuel V.↗

Thermal impacts on nondestructive analysis measurements of uranium hexafluoride

A series of nondestructive assay measurements were performed before, during, and after localized heating of an industrial standard 12B cylinder of uranium hexafluoride to replicate dose variations previously observed on commercial 30B cylinders. While cylinders represent a closed system, enrichment measurements may be impacted by altered uranium and daughter product distributions due to sublimation. A series of measurements were performed to examine localized heating impacts on dose rates, ratios of uranium and daughter product distributions, and measured 235 U enrichment as calculated by two different approaches. The less complex enrichment meter approaches were found to be more robust and less prone to deviation than isotope ratio approaches such as FRAM, but the impact can be minimized by using observable dose changes to avoid impacted areas.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Shifting modes of iron sulfidization at the onset of OAE-2 drive regional shifts in pyrite δ 34 S records

We report total reduced inorganic S isotope ratios (δ 34 S CRS ) shift toward more negative values across much of the southern North Atlantic just before the onset of the Cenomanian-Turonian Ocean Anoxic Event (OAE-2). At the same time, there is no parallel isotopic change in the significantly larger pool of kerogen (organic) S, which indicates that the distribution and S-isotope composition of sulfide in the environment likely did not drive the change in δ 34 S CRS . Here, we investigate possible explanations for the negative shift in δ 34 S CRS values and their divergence from organic S by isolating iron sulfides for morphological identification and grain-specific isotopic analysis using secondary ion mass spectrometry (SIMS). In pre- and syn-OAE-2 sedimentary rocks from Demerara Rise, we find four distinct morphologies of iron sulfides: pyrite framboids (1–20 μm diameter), irregular pyrite aggregates (1–38 μm diameter), large cemented pyrite aggregates (~60 μm diameter), and irregular and cemented aggregates of the pyrite polymorph marcasite (1–45 μm diameter). These different textural groups have distinct S-isotopic compositions that are largely consistent through the onset of OAE-2. As such, the secular change in bulk δ 34 S CRS values likely reflects the changing proportions of these phases stratigraphically across OAE-2. All textural groups feature resolvable intra-grain δ 34 S variability, suggesting that the environments in which they formed were characterized by dynamic sulfide δ 34 S values and/or by partial closed-system distillation. We use grain-specific δ 34 S distributions to rule out shoaling of the chemocline within the sediments as a mechanism for the observed decrease in δ 34 S CRS . Instead, we propose that changes in the reactivity of the iron species delivered to Demerara Rise over the ~200 kyr leading up to the onset of OAE-2 impacted the relative contributions of pyrite with S-isotope signatures reflecting the water column, shallow sediments, and deeper sediments to the bulk sedimentary δ 34 S CRS value. Specifically, the change in iron reactivity at the onset of OAE-2 favored the production of 34 S-depleted large, cemented aggregates and framboids at the expense of more 34 S-enriched irregular aggregates. Our results underscore that bulk δ 34 S CRS measurements integrate multiple reduced phases that form via distinct reaction mechanisms and potentially in different parts of the depositional environment. Grain-specific SIMS analyses dramatically enrich our ability to interpret pyrite isotopic patterns in the geologic record.

58 GEOSCIENCES↗

Ca-isotopes as a robust tracer of magmatic differentiation

The large mass difference (~10%) between the two most abundant isotopes of calcium, 40 Ca and 44 Ca, gives Ca great potential in tracking mass-dependent fractionation during magmatic processes. Resolvable Ca-isotope fractionation during fractional crystallization of magma, particularly by feldspar in evolved melts, has been theoretically inferred but not robustly tested in nature. To further explore the effects of magmatic differentiation on Ca-isotope systematics, we studied the late-Permian alkaline igneous suite of the Øyangen Caldera, Oslo Rift, Norway, consisting of volcanic and intrusive units ranging from basanitic to rhyolitic compositions. Major and trace element variations and modeling demonstrate that the main series of samples (N = 21), including basanites, ring-dyke syenites, and central-dome syenites, likely documents a co-genetic and closed-system fractional crystallization sequence. Our data show minimal δ 44/40 Ca variation (< 0.05 ‰) in the intermediate magma and a marked increase in δ 44/40 Ca in the felsic magma of the Øyangen Caldera (from 0.62 ± 0.02 ‰ to 1.15 ± 0.03 ‰ relative to Ca standard, SRM915a). The systematic increase is best explained by equilibrium isotopic fractionation dominated by alkali feldspar in the fractionating mineral assemblage. This is further supported by strong correlations between δ 44/40 Ca, CaO, and Eu/Eu* in the main-series samples. Implementing a Monte Carlo approach, isotopic modeling of the liquid line of descent using Rayleigh fractionation is highly consistent with the observed Ca-isotope evolution. For the first time, we confirm prominent Ca stable isotope fractionation in felsic-stage differentiation of alkaline magma and constrain the isotope fractionation factors of plagioclase and K-feldspar. Integrated with extant estimations on mineral fractionation factors from the literature, our results suggest increasing fractionation effects of rock-forming minerals with decreasing Ca content. Finally, the affirmation of significant Ca-isotope fractionation in alkaline magma by feldspar empowers the application of Ca as a versatile tracer of crustal evolution, allowing further tests in other magmatic conditions across various planetary objects.

58 GEOSCIENCES↗

A visco-plastic framework for interface processes in sedimentary reservoir rocks at HPHT conditions

As energy operations face the challenge of reservoirs at ever-increasing depths, modelling the response of reservoir rocks at high pressure and high temperature (HPHT) conditions is a crucial step for successfully unlocking new resources. At these conditions, rocks can experience thermal and pressure sensitivity, as well as admit internal phase changes at the pore-solid interfaces that determine their macroscopic response to external loading. It is therefore expected that constitutive laws of Geomechanics should be enriched to accommodate such effects. In this work, a multi-physics constitutive theory for sedimentary rocks is proposed within the general framework of viscoplasticity. The viscosity of the material is assumed to be a function of the temperature, pore-pressure and energy required to alter the inter-granular interfaces. This energy is expressed through the chemical potentials of the phases involved during processes like chemical dissolution/precipitation of the interfaces or mechanical debonding. The resulting flow law and corresponding stress equilibrium are coupled to the energy and mass conservation laws, constituting a closed system of equations, which is solved using the Finite Element simulator REDBACK. A series of numerical tests for performance and calibration is then successfully performed against different types of reservoir rocks, confining pressures and temperatures (from room temperature to over 800K). Here, we show that the mechanical response of sedimentary porous rocks at strains usually achieved in laboratory testing can be explained by accounting for the interface processes taking place at the cementitious material bonding the grains, a process that can also determine the brittle-to-ductile transition of sedimentary rocks.

58 GEOSCIENCES↗

On numerical energy conservation for an implicit particle-in-cell method coupled with a binary Monte-Carlo algorithm for Coulomb collisions

We report conventional particle-in-cell (PIC) methods suffer from enhanced numerical heating (explicit PIC) or cooling (semi-implicit PIC) when coupled with a binary Monte-Carlo algorithm for Coulomb collisions. In this work, a fully-implicit θ-PIC scheme (with adjustable time-biasing parameter 1/2 ≤ θ ≤) is considered. The discrete change in energy of a closed system after a time step for this scheme scales with (1/2 - θ)C θ , where C θ is a positive definite quantity that depends on the frequency spectrum of the energy in the fields. Collisions lead to additional energy in the field fluctuations associated with high-frequency light waves produced by a numerical Bremsstrahlung process, which can result in a large increase in the numerical cooling rate for θ > 1/2. However, for θ = 1/2, energy is exactly conserved. The energy in the field fluctuations on long time scales agrees with that calculated using the equipartition theorem for a classical system in thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detector requirements and simulation results for the EIC exclusive, diffractive and tagging physics program using the ECCE detector concept

This article presents a collection of simulation studies using the ECCE detector concept in the context of the EIC’s exclusive, diffractive, and tagging physics program, which aims to further explore the rich quark–gluon structure of nucleons and nuclei. To successfully execute the program, ECCE proposed to utilize the detector system close to the beamline to ensure exclusivity and tag ion beam/fragments for a particular reaction of interest. Preliminary studies confirm the proposed technology and design satisfy the requirements. Further, the projected physics impact results are based on the projected detector performance from the simulation at 10 or 100 fb -1 of integrated luminosity. Additionally, insights related to a potential second EIC detector are documented, which could serve as a guidepost for future development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Gas generation and intramolecular isotope study in laboratory pyrolysis of the Springfield coal from the Illinois Basin

Position-specific (PS) isotopes of propane have been proposed as a potential geochemical tool to decipher various geological processes (e.g., thermal cracking, biodegradation, H exchange) in natural reservoirs. The limited studies have been conducted on the PS isotopes of propane from the pyrolysis gases from marine shales, and natural gases sourced from lacustrine and marine kerogens, but little is known on gases produced from the humic kerogen. Here, this study investigated the PS δ 13 C of propane in the closed-system pyrolysis of the Springfield coal, Illinois Basin, Indiana, at 310 to 470 °C (Easy %Ro: 0.76 to 3.07). The C kinetic isotope effect (KIE) of CH 4 produced in both this study and previous low-temperature pyrolysis of the same coal indicates the cleavage of C—O bonds is the main generation pathway at the early kerogen cracking stage, followed by the breakdown of alkyl groups. At the wet-gas cracking stage, C 3 H 8 production from thermally stable compounds has a significant influence on the bulk and position-specific C KIE in the pyrolysis of marine Woodford kerogen and Springfield coal. According to the PS δ 13 C of propane, the central site is likely more enriched in 13 C and the δ 13 C of the terminal site is relatively heterogeneous within the propyl group attached to different functional groups of the gas-prone kerogen. Our findings based on the pyrolysis experiments and natural gas samples indicate thermal cracking and biodegradation appear to alter the δ 13 C cen values more significantly than the δ 13 C ter values of propane. The larger magnitude of ΔC c-t in the kMC simulations (Peterson et al., 2018) compared with those from the marine shale and coal possibly implies the non-random distributions of 13 C of propane precursors in the kerogens. As a new dimension of intramolecular isotopic information of propane, the PS δ 13 C values can contribute to fingerprinting the gas origins and identifying the various geological processes (e.g., kerogen cracking, wet-gas cracking, microbial activities) in sedimentary basins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Thermal Processes Using Alternating Electromagnetic Field for Methodical and Selective Release of Hydrogen Isotopes

The controlled and selective release of hydrogen isotopes from Lanthanum-Nickel-Aluminum (LANA) materials was achieved using an alternating electromagnetic field (AEMF). Here, in the presence of the AEMF, the rate of desorption for D 2 was faster than for H 2 . In a closed system, compared to control experiments with helium, hydrogen-loaded LANA samples showed a significant increase in pressure which illustrates the desorption of hydrogen from the applied AEMF. An exponential increase in the hydrogen release rate was observed with increased AEMF strength. Additionally, hydrogen release from LANA confined in a subzero environment (-78°C) using an AEMF was demonstrated. These results demonstrate that the LANA material can be heated directly using alternating electromagnetic fields rather than indirectly to rapidly and selectively desorb hydrogen isotopes. The rapid response and difference in desorption kinetics for different hydrogen isotopes in the presence of the alternating electromagnetic fields provide a promising pathway for applications requiring controlled release of hydrogen and separation of its isotopes.

07 ISOTOPE AND RADIATION SOURCES↗

Experimental Investigation of Barium Sources and Fluid–Rock Interaction in Unconventional Marcellus Shale Wells Using Ba Isotopes

Produced waters from unconventional Marcellus Shale gas wells have anomalously high barium (Ba) concentrations and yield some of the isotopically heaviest Ba measured to date. Experiments were conducted to constrain the source of Ba in these fluids and the controls on barite (BaSO 4 ) precipitation and dissolution in oil and gas wells. Experiments simulating the acidizing stage evaluated the solubility of pure barite and drilling mud in 2 M HCl at 80 °C for periods of 2, 6, and 48 h and resulted in <0.01% barite dissolution with no appreciable change in δ 138 Ba ( 138 Ba/ 134 Ba normalized to NIST standard 3104a). Static autoclave experiments conducted at 66 °C and 20.7 MPa with combinations of ground Marcellus Shale solids and/or barite-bearing drilling mud with synthetic low-Ba fracturing fluid resulted in decreased Ba concentrations in the fluid, with the largest decrease in the shale-only run. Fluid δ 138 Ba values increased by about 0.5‰ as Ba concentrations decreased, consistent with closed-system Rayleigh fractionation. Flow-through experiments in Marcellus Shale core conducted for 28 days resulted in effluent Ba concentrations an order of magnitude lower than the influent, while sulfate concentrations increased over time. Effluent δ 138 Ba values increased over the first 12 days and plateaued at about 1‰ higher than the influent. Modeling suggests a combination of the release of labile shale Ba and barite precipitation. This work indicates that the processes of Ba release from fluid–shale interactions and barite precipitation in fractures and the well bore, while capable of producing high δ 138 Ba fluids, are unlikely to generate fluids with high-Ba concentrations and δ 138 Ba values like those in Marcellus-produced waters. As a result, we find that the release of sulfate from shale pyrite oxidation rapidly catalyzes barite precipitation and that dissolution of drilling mud barite or natural barite in the shale is unlikely to be the major source of Ba in Marcellus-produced waters.

54 ENVIRONMENTAL SCIENCES↗

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗