Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Monte Carlo Potts model”

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.

1,030 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Heat Transfer Fluids as Co‐Diluents in Localized High‐Concentration Electrolytes for High‐Rate Lithium Metal Batteries With Enhanced Safety

Localized high-concentration electrolytes (LHCEs) have been identified as promising electrolyte formulations for lithium metal batteries, due to their effective interphase formation and promotion of compact Li deposition, yet their practical implementation is often limited by reduced ion transport kinetics. In this study, two industrially established fluorinated ethers are identified for the first time in battery research as effective co-diluents as they combine a broad electrochemical stability window with a low viscosity and intrinsic non-flammability. Incorporating these components, commonly used as heat transfer fluids, yields safer, less flammable electrolyte formulations with enhanced ion mobilities. In particular, the ternary co-diluent formulation shows improved ion mobility by reducing the electrolyte's viscosity while limiting excessive ion clustering. Based on the improved electrolyte transport kinetics, lower overvoltages and higher Coulombic efficiencies at current densities ≥ 1 mA cm −2 are achieved with the ternary co-diluent blend, resulting in markedly extended cycle life in an application-oriented zero-excess pouch cell compared with the baseline system. Complementary electrochemical and ex situ analysis of harvested electrodes at moderate current densities reveals no discernible differences in interphase morphology and composition, suggesting enhanced ion mobility as the primary cause of the improved high-rate performance.

electrolyte diluent↗
DOE OSTI2024

Development of Electrolytes under Lean Condition in Lithium–Sulfur Batteries

Lithium–sulfur (Li–S) batteries stand out as one of the promising candidates for next-generation electrochemical energy storage technologies. A key requirement to realize high-specific-energy Li–S batteries is to implement low amount of electrolyte, often characterized by the electrolyte/sulfur (E/S) ratio. Low E/S ratio aggravates the known challenges for Li–S batteries and introduces new ones originated from the high concentration of polysulfides in limited electrolyte reservoir. Here, in this review, the connections between the fundamental properties of electrolytes and the electrochemical/chemical reactions in Li–S batteries under lean electrolyte condition are elucidated. The emphasis is on how the solvating properties of the electrolyte affect the fate of polysulfides. Built upon the mechanistic analysis, different strategies to design lean electrolytes to improve the overall process of Li–S reactions and Li anode protection are discussed.

25 ENERGY STORAGE↗
DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2026

Amazonian fog harbors viable microbes

Fog formation over tropical forests remains poorly characterized, despite its potential role in bioaerosol dispersion and ecosystem processes. Here, we analyzed fog samples collected at the Amazon Tall Tower Observatory using flow cytometry and culture-based techniques to characterize viable microbial communities. Microbial cell concentrations varied over an order of magnitude across 13 fog events, reaching up to 8 × 104 cells per ml of fog water. Flow cytometry consistently detected metabolically active cells, while culturing and mass spectrometry-based identification yielded eight viable bacterial species and seven fungal taxa. The bacteria Serratia marcescens, Ralstonia pickettii and Sphingomonas paucimobilis exhibited seasonal variations in prevalence. The fungal species identified were primarily mesophilic saprophytes and endophytes, commonly associated with soil and plant surfaces. Our findings indicate that fog harbors viable microbes, including Serratia marcescens and Ralstonia pickettii, which may imply a relevance of fog for microbial dispersal, colonization and nutrient cycling in the Amazon rainforest.

Godoi, Ricardo H. (ORCID:0000000247744870)↗
NASA NTRS2026

Atacama Cosmology Telescope A Measurement of Galaxy Cluster Temperatures Through Relativistic Corrections to the Thermal Sunyaev-Zeldovich Effect

The high electron temperature in galaxy clusters (> 1 keV or > 107 K) leads to corrections at the level of a few percent in their thermal Sunyaev-Zeldovich effect signatures. Both the size and frequency dependence of these corrections, which are known as relativistic temperature corrections, depend upon the temperature of the objects. In this work we exploit this effect to measure the average temperature of a stack of Compton-y selected clusters. Specifically, we apply the “spectroscopic method” and search for the temperature that best fits the clusters’ signal measured at frequencies from 30 to 545 GHz by the Atacama Cosmology Telescope and Planck satellite. We measure the average temperature of clusters detected in the Atacama Cosmology Telescope maps to be 8.5 2.4 keV, with an additional systematic error of comparable amplitude dominated by passband uncertainty. Upcoming surveys, such as the Simons Observatory and CMB-S4, have the potential to dramatically improve upon these measurements and thereby enable precision studies of cluster temperatures with millimeter observations. The key challenge for future observations will be mitigating instrumental systematic effects, which already limit this analysis.

William R Coulton↗
DOE OSTI2026

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗
NASA NTRSDate not supplied

Evaluation of Low-Energy Hydrogen Separation Method Using Metal-Organic Frameworks (MOFs) for Closed-Loop ECLSS Air Revitalization (CLEAR)

The State-of-the-Art (SOA) air revitalization architecture onboard the International Space Station (ISS) recovered approximately 50% of the oxygen (O 2 ) from metabolic carbon dioxide (CO 2 ) via the Sabatier process from 2011 to 2017. O 2 recovery is currently constrained by the limited availability of reactant hydrogen (H 2 ) preventing complete conversion of CO 2 to H 2 O. Increasing O 2 recovery within Closed-Loop ECLSS is essential to reducing resupply mass for long-duration manned missions; specifically focusing on water (H 2 O) which supplies H 2 for Sabatier via water electrolysis. Past ground test endeavors at Marshall Space Flight Center (MSFC) have attempted to recover H 2 from Sabatier-produced CH 4 using technologies such as carbon vapor deposition (CVD) and plasma pyrolysis. The byproducts of these technologies can act as a catalyst poison or reactor deadload to the Sabatier reactor. Hydrogen separation techniques must be utilized to maintain the Sabatier catalyst during gas recycling and must be scalable, non-energy intensive, and safe to operate in a habitation setting. Research indicated that metal-organic frameworks (MOFs) could meet these criteria and were tested for their capability to capture the various carbon-based gaseous products of CVD and plasma pyrolysis such as acetylene (C 2 H 2 ), ethylene (C 2 H 4 ), ethane (C 2 H 6 ), and carbon monoxide (CO) which would purify the hydrogen gas stream passing through the MOF. A sub-scale adsorption column was developed by Marshall Space Flight Center to test three MOF candidates against a synthetic gas mixture comprised of process-relevant carbonous gases and hydrogen to evaluate the separation capability of the MOFs. The results of the hydrogen separation capability, isothermal desorption capability, and demonstrated cyclic reuse of the MOF are presented in this paper.

Kagen Crawford↗
DOE OSTI2026

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational modeling of coupled mechanical damage and electrochemistry in ternary oxide composite electrodes

Performance degradation of ternary layered oxide cathodes largely originates from their loss of structural integrity in cyclic usage. Mechanical damage, such as intergranular fracture of the active particles, is not only a mechanical cleavage process but also interferes with electrochemical kinetics such as infiltration of liquid electrolyte, surface corrosion of the constituent primary particles, and may eventually isolate the primary grains from the electron conducting network. Here, in this work, we develop a computational framework that integrates electrochemistry of a LiNi x Mn y Co 1−x−y O 2 (NMC) composite cathode with mechanical damage of the active particles. To fully examine the intricate chemomechanical behavior of the electrode, we evaluate the effects of the anisotropic material properties, the influence of mechanical potential on Li transport, and the concurrent intergranular fracture and electrolyte penetration along the grain boundaries upon multiple cycles. Electrolyte infiltration benefits capacity retention but aggravates further mechanical damage by corrosion. Structural failure mostly occurs in the first charging due to the anisotropic mechanical strain between the primary grains, while the resulting damage remains stable in the later few cycles. The results are consistent with experimental observations and the integration of electrochemistry and mechanical failure enables a step further understanding of the complex mechanism of battery degradation.

Battery degradation↗
DOE OSTI2026

Mitigation of polysulfide shuttle effect in Li-S batteries through catalytic disproportionation reaction

Polysulfides are poorly retained within porous cathodes and readily diffuse into the electrolyte over time, leading to the well-known shuttle effect that undermines the reversibility of Li-S batteries. Here, in this study, we demonstrate that catalytic disproportionation of polysulfides provides an effective pathway to suppress this process by rapidly converting dissolved species into solid sulfur and sulfides, thereby preventing their migration into the electrolyte. Fundamentally, the sluggish kinetics of sulfur redox reactions are responsible for the accumulation and redistribution of soluble polysulfides in the bulk electrolyte. By accelerating these kinetics, catalyzed disproportionation not only confines sulfur within the conductive cathode matrix but also promotes the homogeneous precipitation of Li₂S₂/Li₂S, which enhances electrochemical reversibility and cycling stability. Using nitrogen-doped carbon (NC800) as a model catalyst, we reveal its ability to drive a pseudo-16-electron reduction pathway, leading to a single dominant Li₂S product and uniform deposition within the porous framework. In contrast, a non-catalytic carbon (KB) yields multiple polysulfide intermediates and heterogeneous deposition. The mechanistic insights provided here highlight the pivotal role of catalytic disproportionation in reshaping sulfur redox pathways and offer a rational strategy for mitigating polysulfide shuttling in practical Li-S pouch cells.

25 ENERGY STORAGE↗
DOE OSTI2025

Advanced Model Development for Large Eddy Simulation of Oxy-Combustion and Supercritical Carbon Dioxide Power Cycles

A joint experimental and numerical study is performed to observe the characteristics of a supercritical carbon dioxide turbulent mixing layer in the presence of strong nonlinearities in the thermodynamic and transport properties. A bespoke experimental setup is designed and employed for this purpose and provides insight into macroscopic mixing behavior. The mixing is experimentally observed using two techniques: shadowgraphy and spontaneous Raman scattering. Qualitative and quantitative intensity fields obtained via these techniques yield instantaneous and mean density data. Spanwise temperature data is also collected using analogue resistance temperature detectors. These measurements are used to quantify the level of mixed material within the field. The experimental data are supplemented by a companion high-fidelity numerical study. The numerical results are obtained through fully resolved, three-dimensional direct numerical simulation. The numerical dataset permits observation of the near-field mixing characteristics, which are difficult to measure experimentally due to the rapid dynamics and sharp thermophysical gradients in this area. Qualitative field visualizations are presented, followed by quantitative mixed material results and observations regarding thermodynamic property trends at select locations within the field. One-dimensional spectra of the turbulent kinetic energy and solenoidal dissipation are provided to observe the spectral characteristics of the flow. Reynolds stress anisotropy is analyzed graphically through anisotropy invariance maps (Lumley triangles). The mixing quantification, spectral data and anisotropy analysis of a flow at these thermodynamic conditions represent the main outcomes of the work.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Electronic glasses from a broken gauge symmetry in disorder-free systems

Glass phases can be stabilized by quenched disorders, as in most spin-glass materials, or self-generated through kinetic freezing in disorder-free systems. A canonical example of the latter is structural glasses, which have been extensively studied for many decades. Yet, how the rugged energy landscape of a glass phase is spontaneously generated in disorder-free systems remains one of the key questions in glass physics. Here, in this work, we present a general electronic mechanism for the emergence of glassy phase using the example of itinerant electrons coupled to XY spins on a lattice. This model can also be viewed as the mean-field theory of a superconducting system with attractive density-density interactions. Intriguingly, the electron gauge symmetry in the strong pairing limit gives rise to a macroscopic degeneracy of XY spins. In the presence of electron hopping that breaks the gauge symmetry, the lifting of the extensive degeneracy leads to a glass phase with disordered pairings. Our findings highlight a scenario in which a glassy state originates from the breaking of quantum gauge symmetry without quenched disorders.

XY model↗
NASA NTRSDate not supplied

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES↗
DOE OSTI2026

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials↗
DOE OSTI2023

Electrolyte-Induced Restructuring of Acid-Stable Oxygen Evolution Catalysts

Crystalline metal oxide catalysts operating under oxygen evolution reaction (OER) conditions invariably restructure, resulting in active sites with hydroxo/oxo species in an amorphous environment. An increase in the population of terminal hydroxo/oxo species (i.e., edge sites) facilitates proton-coupled electron-transfer (PCET) kinetics for oxygen generation and thus improves catalyst competency. While amorphous films benefit from a greater density of active sites, they suffer from diminished charge transport as compared to that of extended crystalline lattices. Managing this amorphous–crystalline dichotomy is essential when designing OER catalysts, which we highlight with the examination of electrodeposited PbO x materials, which historically are very poor OER catalysts. Along these lines, the presence of phosphate during PbO x electrodeposition truncates the growth of an extended lattice owing to its strong bonding to oxide surfaces to afford an amorphous catalyst film (A-PbO x ) with significant charge-transfer resistance (138 ± 42 Ω) and poor OER kinetics (420 ± 105 mV dec –1 Tafel slope). Conversely, electrodeposition of Pb 2+ in the presence of less coordinating electrolytes such as nitrate affords crystalline β-PbO2 with improved charge-transfer resistance (42.6 ± 1.1 Ω), though still poor OER kinetics (134 ± 36 mV dec –1 Tafel slope). By operating amorphous A-PbOx in less coordinating electrolytes, however, a new partially crystalline material can be generated (μc-PbO x ) with further reduced charge-transfer resistance (33.0 ± 1.4 Ω) and improved OER kinetics (70 ± 15 mV dec –1 Tafel slope). The enhanced OER activity of μc-PbO x is the result of coupling the high edge-site population of an amorphous PbOx phase with crystalline-like charge transport properties. Finally, the ability to use an electrolyte to induce OER activity in an inactive amorphous form of PbO x highlights the benefits of optimizing the amorphous–crystalline phase compositions in the design of active OER catalysts.

catalysts↗
DOE OSTI2026

Single-Nucleon Transfer on Unstable 59 Cu Probes the NiCu Cycle in Astrophysical X-Ray Bursts

Recent models of the rapid proton (𝑟⁢𝑝) capture process indicate that a competition between the 59 Cu⁢(𝑝,𝛾) ⁢60 Zn and 59 Cu⁢(𝑝,𝛼) ⁢56 Ni reactions may result in the formation of a nickel-copper (NiCu) cycle that traps the flux of material between 56 Ni and 60 Zn . Here, in this work, we report the identification of 15 proton-unbound levels in 60 Zn , populated via 59 Cu⁢(𝑑,𝑛) transfer, which govern the rate of the 59 Cu⁢(𝑝,𝛾)⁢ 60 Zn reaction in XRBs. Precise excitation energies for levels in 60 Zn were obtained from observed 𝛾 decays, and spectroscopic factors were determined from angle-integrated cross sections. Incorporating these results into stellar-model calculations, we find that with experimentally constrained uncertainties a NiCu cycle in XRBs is indeed possible, though we limit its branching strength to less than 38%. While modest, such a branching has significant impact on the light curve, motivating further studies of the relevant rates. Our calculations also indicate that a significant NiCu cycle leads to an increase in the amount of odd-𝐴 nuclei in the burst ashes, which may affect Urca cooling processes in neutron star crusts.

direct reactions↗