Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic equilibrium”

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 217 records · Page 12

Frontier nonequilibrium materials science enabled by ultrafast electron methods

Most innovation in materials science and engineering resides in our ability to understand and control the intimate relationship between the structure of materials and their properties. Conventionally, the only route to discovering innovative new material properties has been to explore the structural and compositional phase space that is accessible at (or near) thermodynamic equilibrium. The characterization, manipulation and, ultimately, control of material properties far from equilibrium offers almost completely untapped possibilities for uncovering novel states and phases of materials. Investigations of materials far from equilibrium require the development of new techniques, which are capable of following dynamic processes in materials with extreme spatiotemporal resolution. Thus, ultrafast electron-based methods have become a major new frontier in materials science due to the capability of following dynamics on time scales as short as femtoseconds with the high spatial resolution and sensitivity afforded by electrons. The articles in this issue provide an exciting cross section of the novel research directions that have been enabled at this frontier. Imaging on ultrafast time scales carries critical information on phase transitions, fundamental processes involved in light-harvesting, magnetic, and plasmonic dynamics, the coupling of electronic and nuclear degrees of freedom in materials and has revealed previously hidden, nonequilibrium metastable states of matter that have no equilibrium analog.

36 MATERIALS SCIENCE↗

Effect of Impurity Atoms on Thermonuclear Yield

The presence of impurity atoms in an inertial confinement fusion (ICF) implosion can greatly affect the burning of the thermonuclear fuel. This is a serious concern for attaining ignition as high atomic number atoms, such as carbon, can mix into the hot spot and reach the deuterium–tritium (DT) fuel, absorb energy from the implosion, and radiate away that energy. The temperature of the DT fuel drops and reduces the reactivity of the fuel resulting in lower thermonuclear yield. Fuel impurities can also be a useful diagnostic as when higher-Z atoms, such as Ar, Kr, and Xe, are purposefully mixed into the fuel. From the spectra emitted by these atoms, the temperature and density of the fuel can be deduced. The long history of using such tracer atoms by the national ICF program shows the great advances in experimental technique, diagnostic instruments, and atomic-physics modeling over the past 50 years. A particular set of experiments, the “High-Z” campaign, is reviewed in depth as this was the most systematic study of tracer element composition and quantity ever reported. Difficulties in interpreting the results of those experiments due to the limited code capabilities at that time are discussed. A main concern for the High-Z campaign was the then current state of non-Local Thermodynamic Equilibrium (nLTE) atomic physics models. The present state of nLTE modeling and its implementation in Los Alamos’s multi-physics codes is discussed in this report. To further the understanding of the complex interaction of high-Z atoms in the gas, a survey of relevant existing experimental platforms at the National Ignition Facility (NIF) and the Omega Laser Facility are presented along with the current state of pertinent diagnostics at NIF and Omega. This review concludes with a discussion of the trade-off available between temperature and volume of the imploded ICF capsule with the goal of maximizing the total yield.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Opacity Table Generation

Radiation opacity is a key property of material systems which describes the absorption and scattering of photons. Such processes are an important part of energy transport, diagnostics, and astrophysical phenomena. Recent experiments in iron opacity and continuum lowering have sparked renewed interest in and careful vetting of opacity models. Even within the local thermodynamic equilibrium approximation, opacity varies over a large parameter space: material, temperature, density, and energy. Opacity tables that span wide regimes of this parameter space are of prime importance to radiation-hydrodynamics codes. Accurate opacity calculations cannot be achieved inline during a radiation hydrodynamics simulation since they involve accounting for all possible electron transitions of all important atomic states of a thermal ensemble. Livermore’s current framework utilizes precalculated tables at many densities, temperatures, and energies. Last year, the opacity theory team completed an L2 milestone concerning the development of a new code, Opus, which served the dual purpose of creating a modern code infrastructure and version control as well as the basis for training new members of that team, a critical need as key members near or are passed retirement. For the current L2 milestone, we exercised and validated Opus on three elements: boron, carbon, and nitrogen, which have few enough electrons to allow a careful convergence study over many input parameters. Here, we report on the production of new tables, show their convergence with respect to important parameters, state some of their limitations, and document a mostly automated production process. In this section, we address the milestone completion criteria, found on the title page and the subsection headings. In section 2, we validate Opus against Tycho and compare with TabOp and Atomic. Section 3 details some of the physics implemented in Opus and relevant to this report. We then cover details of the density-temperature grid in section 4 and do a spectral comparison between different parameter settings and different codes in section 5. After summarizing these results in section 6, we give instructions how to run the opacity table scripts in appendix A and document the scripts and input/output files in appendix B.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Novel Catalytic Membrane Reactor for DME Synthesis from Renewable Resources

Production of liquid fuels or chemicals from CO 2 (captured from the air or flue gases) and renewable hydrogen presents a new approach to producing clean fuels domestically. While significant progress has been made in the area of renewable electricity generation from solar and wind, a large gap remains with respect to the production of renewable liquid fuels/chemicals. Other processes for producing liquid fuels/chemicals from renewable electricity are constrained by thermodynamic limitations, making them prohibitively expensive and impractical. The team is overcoming these limitations and developing catalytic membrane reactor processes with high yields and low energy penalties. Supported by the Advanced Research Projects Agency-Energy (ARPA-E) of the US Department of Energy (DOE), GTI Energy and partners have been developing a technology for the production of renewable dimethyl ether (DME) from carbon dioxide (CO 2 ) and renewable hydrogen (H 2 ) using a novel catalytic membrane reactor and demonstration of this system at a scale of 1 kg/day. DME is a clean-burning, non-toxic fuel with a high cetane value (55-60), making it an excellent diesel alternative. DME can be stored as a liquid under moderate pressure, eliminating the need for the high-pressure containers used for CNG or cryogenics, as in the case of LNG. DME is also approved as a renewable fuel under the U.S. Environmental Protection Agency’s Renewable Fuels Standard (RFS), making it eligible for Renewable Identification Numbers (RINs) credits. By producing DME through the catalytic conversion of captured CO 2 and renewable H 2 , this process will produce renewable liquid transportation fuel and a means of large-scale utilization of captured CO 2 . In the DME synthesis process, CO 2 and H 2 are fed to a hollow fiber catalytic membrane reactor at 300-600 psig that contains a bi-functional catalyst that combines two reactions, methanol synthesis (CO 2 + 3H 2 → CH 3 OH + H 2 O) and methanol dehydration (2CH 3 OH → CH 3 OCH 3 + H 2 O), into a one-step process to produce DME. The bifunctional catalyst converts methanol to DME, enabling higher overall CO 2 conversion. A Cu/ZnO/ZrO 2 /Al 2 O 3 (CZZA) catalyst is used for methanol synthesis and is coupled with a zeolite catalyst H-ZSM-5 for dehydration. This one-step process intensifies a process that would otherwise require multiple reaction steps. However, combining these two reactions results in increased water production which inhibits catalytic activity. Here, the Na + -gated, water-transport membrane (Science, vol. 367, pp. 667, 2020), removes water in situ, shifting the thermodynamic equilibrium towards product formation while decreasing kinetic inhibition from water adsorption onto the catalyst surface. The Na + gated, water-transport nanochannel membrane showed H 2 O/CO 2 selectivity of 560 at 250 °C and 300 psig for H 2 O/CO 2 /CO/H 2 /MeOH gas mixtures. The selectivities of H 2 O/H 2 , H 2 O/CO, and H 2 O/MeOH were 190, 170, and 80, respectively. In a laboratory-scale membrane reactor, DME synthesis testing using this membrane, a DME production rate of 440 g DME /kg cat /h was achieved at 260 °C and 550 psig. Compared to the packed bed reactor, the CO 2 conversion and DME production rate in the membrane reactor were 80% and three times higher, respectively. A prototype test system (1 kg/day) was designed, constructed, and tested. A DME production rate of 1.31 kg/day and a DME productivity of 360 g/h/kg were achieved in the prototype membrane reactor. Good stability was demonstrated during 150-h continuous operation and multiple startups/shutdowns tests.

10 SYNTHETIC FUELS↗

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spectroscopic Investigations of Molecular Symmetry Breakdown

This project addressed precise measurements of the spectra of molecules to address fundamental questions in chemical physics in the realm of nuclear spin symmetries and their breakdown in molecules, and parity non-conservation caused by the weak nuclear force. Nuclear spin symmetry is intimately connected to parity in symmetrical molecules and plays a central role in the physical and chemical properties of symmetric molecules, their collisional relaxation, and the establishment of thermodynamic equilibrium in an ensemble. Such effects are present in every molecular system or process, but their small size means they are often only directly observable under special conditions. However, their measurement is of fundamental, and often practical importance, and strongly relates to the themes of coherence in light and matter in the recent Basic Energy Sciences Advisory Committee (BESAC) report “Challenges at the Frontiers of Matter and Energy: Transformative Opportunities for Discovery Science”. We begin by addressing lifetime broadening of spectroscopic lines caused by self-collisions in acetylene, where efficient resonant rotational energy transfer between a pair of molecules separated in energy by rotational quantum number ΔJ = ± 2, can result in a different pressure broadening for ortho- and para- nuclear spin state molecules. This work was prompted by a controversial report [Iwakuni et al. Phys. Rev Letts. 117 (2016) 143902] of the observation of such an effect in the v 1 +v 3 combination band transition of acetylene discussed in the original proposal. After delays caused by the COVID-19 pandemic, we have completed these measurements and a paper was published in the Journal of Chemical Physics describing the work. The measurements pointed out that the use of an approximate model in the analysis of the original work was the cause of the erroneous conclusion that ortho and para acetylene molecules had substantially different collisional lifetime broadening. In the process, much improved broadening measurements have been made that will permit more reliable modeling of the effects in the future. Subsequently, we have completed a new set of measurements of self- and nitrogen-collisional broadening at a series of low and intermediate temperatures to characterize the temperature dependence of the effects in the molecule. This work is currently being analyzed and will be the subject of a paper to be published in the relatively near future. In addition, some very preliminary measurements of CH 4 as described in the proposal, have been completed. If facilities and time permit, these will continue as an undergraduate research project.

74 ATOMIC AND MOLECULAR PHYSICS↗

Report on use of Inoculants in Missile Application Alloys

This report documents the status of current inoculant research relevant to missile application alloys and MTCR control language. The information is intended to provide data on current inoculants for us determining the current state of development and identifying potential research directions. Although there has been significant scientific research into the development and synthesis of inoculants, their current availability is limited is traditional powder inoculants employed during casting processes. However, research continues the development of complex oxides, ribbon materials, high entropy alloys, and other inoculant product forms, including the use of inoculants in the melt pools produced during additive manufacturing. Research to date has focused primarily on aluminum and steel alloys with emphasis on refining grain structures and evolving equiaxed morphologies while increasing strength and castability. The primary inoculants in steel and cast irons include TiN, SiC, FeSi75, and Ce which have increased strength properties. Chief inoculants for Al alloys often include TiC, SiC, Al3Sc(x) and TiB2 to aid in precipitation and refinement. Ti and Ni alloys have fewer research activities involving inoculants, although TiN, TiB, ZrN and LaB6 (for Ti alloys) and WC, Co3FeNb2, and CrFeNb (for Ni alloys) have been used. Sic, Al2O3, Mg and Ti are key inoculants for Mg alloys. Multiple cast alloys from each of the material classes demonstrated increased strength and performance properties using inoculants, with several approaching requirements applicable to missile service environments. The continued evolution of advanced manufacturing capabilities is making it easier to produce high temperature near net shape structural materials using inoculant powders. These shapes may include the geometric shapes addressed within the MTCR (tubes and limited wall thicknesses). The use of inoculants may enable further development of high temperature alloys into near net shapes traditionally produced via casting processes due to limited ductility. This may decrease material and manufacturing costs. In addition, inoculation provides controlled kinetics and achievable chemical segregation that enables potential for far-from equilibrium thermodynamic microstructures and chemistries that could provide new metastable alloy states and subsequent properties to address co-design engineering constraints, including needs for increased strength and ductility. It is recommended that specific material combinations within these alloy classes be carefully watched as the materials evolve, with controls aimed at those having material properties above current MTCR levels. This specifically includes the use of refractory inoculants in alloys, and the application of inoculants in high strength and high temperature alloys via additive manufacturing processes, with care to link capabilities to product forms similar to the current requirements on tube geometries and material feed stocks. The continued development of nanoparticle inoculants will increase strength and ductility of high strength castings and additive manufactured metallic components. For example, adding inoculants into the casting of maraging steels and other precipitation strengthened alloys may drastically elevate mechanical properties above the control limit of current regulations.

36 MATERIALS SCIENCE↗

Final Scientific/Technical Report Rapid Design and Manufacturing of High-Performance Materials for Turbine Blades

This research demonstrated the concept of carbide precipitation-strengthened refractory high entropy alloys (RHEA). The advantage of a precipitation strengthened alloy is all phases are in thermodynamic equilibrium promoting microstructural stability, and consequently retention of properties at elevated temperatures. Additionally, as with any precipitation strengthened – or age hardened alloy, components can be heat-treated after manufacturing to manipulate the microstructure and optimize properties for performance. This is an advantage of the precipitation strengthened alloys over composites and ceramics materials, where the microstructure and properties are essentially fixed upon the initial materials manufacturing stage. High throughput (HT), multi-scale computer modeling was used to identifying novel RHEA compositions with desired characteristics needed for precipitation strengthening. The results showed that carbides precipitated and consequently, the strength of the alloys (measured in compression) increased after heat-treatment, which is the desired effect. The project also demonstrated the feasibility of producing articles from RHEA by additive manufacturing (AM).

36 MATERIALS SCIENCE↗

Dehydration Membrane Reactor for Production of Valuable Chemicals from CO 2 and H 2

GTI Energy and partners have been developing a technology for production of liquefied petroleum gas (LPG) from carbon dioxide (CO 2 ) and hydrogen (H 2 ) using a novel catalytic membrane reactor. The reactor contains a bi-functional catalyst for methanol synthesis and LPG synthesis, resulting in LPG production in a single reactor. A transformational dehydration membrane is used to remove water in situ, shifting the thermodynamic equilibrium towards product formation. As a result, CO 2 conversion as high as 90.2% and LPG yield as high as 60.5% were obtained at 300°C and 20 bara in the membrane reactor, which significantly exceed the literature results of the traditional reactors

03 NATURAL GAS↗

Rapid Design and Manufacturing of High-Performance Materials for Turbine Blades

This research demonstrated the concept of carbide precipitation-strengthened refractory high entropy alloys (RHEA). The advantage of a precipitation strengthened alloy is all phases are in thermodynamic equilibrium promoting microstructural stability, and consequently retention of properties at elevated temperatures. Additionally, as with any precipitation strengthened – or age hardened - alloy, components can be heat-treated after manufacturing to manipulate the microstructure and optimize properties for performance. This is an advantage of the precipitation strengthened alloys over composites and ceramics materials, where the microstructure and properties are essentially fixed upon the initial materials manufacturing stage. High throughput (HT), multi-scale computer modeling was used to identifying novel RHEA compositions with desired characteristics needed for precipitation strengthening. Designs were validated by producing small ingots via are melting. The results showed that carbides precipitated and consequently, the strength of the alloys (measured in compression) increased after heat-treatment, which is the desired effect. The project also demonstrated the feasibility of producing articles from RHEA by additive manufacturing (AM). Electron beam melting (EBM) and laser direct energy deposition (L-DED) additive manufacturing methods, were explored with various processing parameters were interrogated for both methods. Sound (crack-free and dense) precipitation strengthened RHEA samples were produced via EDM AM, demonstrating the feasibility of the concept.

36 MATERIALS SCIENCE↗

Geochemical Modeling of the Confining Unit at a CO 2 Storage Site in the Williston Basin, North Dakota

Conference presentation for American Association of Petroleum Geologists (AAPG) Carbon Capture, Utilization, and Storage (CCUS) Conference, Houston, TX, April 25–27, 2023. This geochemical modeling study investigates the evolution of geochemical reactions among CO 2 , formation water, and rocks of the two confining units above and below the CO 2 storage unit at one of the permitted CO 2 storage sites in the Williston Basin. The model is built on PHREEQC-coded chemical equilibrium thermodynamics.

20 FOSSIL-FUELED POWER PLANTS↗

Increasing Compressed Gas Energy Storage Density Using CO2–N2 Gas Mixture

This paper demonstrates a new method by which the energy storage density of compressed air systems is increased by 56.8% by changing the composition of the compressed gas to include a condensable component. A higher storage density of 7.33 MJ/m3 is possible using a mixture of 88% CO2 and 12% N2 compared to 4.67 MJ/m3 using pure N2. This ratio of gases representing an optimum mixture was determined through computer simulations that considered a variety of different proportions from pure CO2 to pure N2. The computer simulations are based on a thermodynamic equilibrium model that predicts the mixture composition as a function of volume and pressure under progressive compression to ultimately identify the optimal mixture composition (88% CO2 + 12% N2). The model and simulations predict that the optimal gas mixture attains a higher energy storage density than using either of the pure gases.

25 ENERGY STORAGE↗

Microstructural Assessment of a Multiple-Intermetallic-Strengthened Aluminum Alloy Produced from Gas-Atomized Powder by Hot Extrusion and Friction Extrusion

An aluminum (Al) matrix with various transition metal (TM) additions is an effective alloying approach for developing high-specific-strength materials for use at elevated temperatures. Conventional fabrication processes such as casting or fusion-related methods are not capable of producing Al–TM alloys in bulk form. Solid phase processing techniques, such as extrusion, have been shown to maintain the microstructure of Al–TM alloys. In this study, extrusions are fabricated from gas-atomized aluminum powders (≈100–400 µm) that contain 12.4 wt % TM additives and an Al-based matrix reinforced by various Al–Fe–Cr–Ti intermetallic compounds (IMCs). Two different extrusion techniques, conventional hot extrusion and friction extrusion, are compared using fabricating rods. During extrusion, the strengthening IMC phases were extensively refined as a result of severe plastic deformation. Furthermore, the quasicrystal approximant IMC phase (70.4 wt % Al, 20.4 wt % Fe, 8.7 wt % Cr, 0.6 wt % Ti) observed in the powder precursor is replaced by new IMC phases such as Al3.2Fe and Al45Cr7-type IMCs. The Al3Ti-type IMC phase is partially dissolved into the Al matrix during extrusion. The combination of linear and rotational shear in the friction extrusion process caused severe deformation in the powders, which allowed for a higher extrusion ratio, eliminated linear voids, and resulted in higher ductility while maintaining strength comparable to that resulting from hot extrusion. Results from equilibrium thermodynamic calculations show that the strengthening IMC phases are stable at elevated temperatures (up to ≈ 600 °C), thus enhancing the high-temperature strength of the extrudates.

36 MATERIALS SCIENCE↗

Mechanical Activation and Cation Site Disorder in MgAl 2 O 4

The synthesis and crystallographic site occupancy were investigated for MgAl 2 O 4 with and without mechanical activation of the precursor powders. Heating to 1200 °C or higher resulted in the formation of a single spinel phase regardless of whether the powders were mechanically activated or not. Neutron diffraction analysis was used to determine cation site occupancy and revealed that mechanical activation resulted in a lower degree of cation site inversion compared to the nonactivated materials, which indicated that the powders were closer to thermodynamic equilibrium. This is the first study to characterize the effects of mechanical activation on crystallographic site occupancy in magnesium aluminate spinel using neutron diffraction.

36 MATERIALS SCIENCE↗

Determination of Sodium Abundance Ratio from Low-resolution Stellar Spectra and Its Applications

We present a method to determine sodium abundance ratios ([Na/Fe]) using the Na i D doublet lines in low-resolution (R ~ 2000) stellar spectra. As stellar Na i D lines are blended with those produced by the interstellar medium, we developed a technique for removing the interstellar Na I D lines using the relationship between extinction, which is proportional to E(B – V), and the equivalent width of the interstellar Na i D absorption lines. When measuring [Na/Fe], we also considered corrections for nonlocal thermodynamic equilibrium (NLTE) effects. Comparisons with data from high-resolution spectroscopic surveys suggest that the expected precision of [Na/Fe] from low-resolution spectra is better than 0.3 dex for stars with [Fe/H] > –3.0. We also present a simple application employing the estimated [Na/Fe] values for a large number of stellar spectra from the Sloan Digital Sky Survey (SDSS). After classifying the SDSS stars into Na-normal, Na-high, and Na-extreme, we explore their relation to stars in Galactic globular clusters (GCs). We find that while the Na-high SDSS stars exhibit a similar metallicity distribution function (MDF) to that of the GCs, indicating that the majority of such stars may have originated from GC debris, the MDF of the Na-normal SDSS stars follows that of typical disk and halo stars. As there is a high fraction of carbon-enhanced metal-poor stars among the Na-extreme stars, they may have a non-GC origin, perhaps due to mass-transfer events from evolved binary companions.

79 ASTRONOMY AND ASTROPHYSICS↗

Production of Very Light Elements and Strontium in the Early Ejecta of Neutron Star Mergers

Abstract We study the production of very light elements ( Z < 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10 −6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10 −6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r -process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗