Search NASA⌕ Search

SEARCH · Search NASA

Results for “density evolution”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comparative Study on the Photoelectrochemical Hydrogen Production from the Layout Optimization in a Unit Cell

Solar-driven water splitting is achieved by connecting photovoltaics with electrolysis. Recently, integrated approaches of the two configurations for photoelectrochemical reactions have shown great potential in a combined unit, but strategies to address further material and cell development face significant scientific challenges. Here, in this work, we present an experimental demonstration of photoelectrochemical (PEC) unit-cells with efficient separate dual compartments for oxygen evolution and hydrogen evolution reactions. The design and fabrication of cell architectures with various degrees of integration are investigated efficient and sustainable PEC processes. The effect of the nanostructured silicon (Si) photoelectrodes, including the adoption of layer architectures, surface protection layer deposition, and membrane electrode assembly is further studied to optimize the design parameters of the PEC unit-cells. This developed PEC cell with the Si photocathode achieved a photocurrent density of -4.8 mA cm –2 at -2.0 V cell and -2.87 mA cm –2 at -1.6 V cell . It clearly provides a critical milestone for unbiased solar water splitting. Practical solar to hydrogen (STH) can be directly estimated in investigating PEC performance in this unit-cell as assembled with a solar cell. It is a significantly meaningful step forward in practical solar fuel production.

08 HYDROGEN↗
DOE OSTI2024

High-order harmonics and the reverse of the squaring up process in the triangular-lattice magnet HoPdAl 4 ⁢Ge 2

We uncover high-order harmonics and the reverse of the squaring-up process, in terms of analyzing the evolution of magnetic orders in a centrosymmetric layered triangular-lattice magnet HoPdAl 4 ⁢Ge 2 , based on a detailed study of the crystal structure, magnetic susceptibility, magnetization, heat capacity, and magnetic structure. Temperature dependencies of magnetic susceptibility and heat capacity show two magnetic transitions at T N = 10.5 K and T t = 5.5 K. Below T N , Ho 3+ spins order antiferromagnetically as a transverse spin-density wave with the propagation vector k 1 = (001.5 – δ) with δ ≈ 0.18. Upon further cooling through T t , the high-order harmonics with k n = (001.5 – nδ), n = 3, 5, 7 develop, suggesting a squaring-up process. It is surprising that the squaring-up process does not continue down to 0 K but reverses the trend below 3 K. Magnetic-field-induced metastable transitions were observed in M⁡(H) curves with the fields applied both parallel and perpendicular to the triangular-lattice plane. Neutron-diffraction results suggest that the magnetization process in HoPdAl 4 ⁢Ge 2 involves the conversion of k n = (001.5 – n⁢δ) with n = 1, 3, 5, 7, to a ferromagnetic k F = (000) component. It is worth noting that the already weak seventh harmonic magnetic peak is enhanced by applying a small magnetic field in the ab plane at 1.5 K or warming up to 3 K, accompanied by a slight decrease of δ.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Effect of Specimen Thickness on Fracture Toughness and Plane Stress to Plane Strain Transition in Medium-Density Polyethylene

This study investigates the effect of specimen thickness on fracture toughness and the transition from plane stress to plane strain conditions in Medium-Density Polyethylene (MDPE) using Single Edge Notched Bend (SENB) specimens. Three thickness groups (7.5 mm, 9.0 mm, and 12.1 mm) were tested following ASTM D5045 protocol. Conditional stress intensity factors (KQ) increased from 2.2 MPavm to 2.8 MPavm with increasing thickness, demonstrating significant size dependency. Confocal microscopy revealed a 34% reduction in maximum crack tip opening displacement (CTOD) from 0.478 mm to 0.314 mm as thickness increased, with plastic zone lateral extent decreasing by 32%. This quantitative evidence validates Irwin's theoretical prediction of plastic zone size transition from plane stress r_y˜ (1/2p) (K_I/s_y )^2 to plane strain conditions r_y˜ (1/6p) (K_I/s_y )^2.Scanning electron microscopy confirmed progressive suppression of shear lips and evolution from ductile tearing with extensive polymer chain drawing in thin specimens to localized crazing in thick specimens. Despite all specimens satisfying the ASTM D5045 criterion (B = 2.5(KQ/s_y )^2 ) ,residual shear lips persisted even in the thickest specimens, demonstrating that nominal size requirements alone are insufficient for achieving complete plane strain conditions in highly ductile polymers. The findings emphasize the necessity of integrated analytical and morphological assessment for valid fracture toughness characterization, with critical implications for life assessment and integrity prediction in thick-section polymer components such as natural gas distribution pipelines.

Medium-Density Polyethylene↗
DOE OSTI2024

Thermal stability of 3D interface Cu/Nb nanolaminates

Nanocrystalline alloys are promising structural materials yet lack thermal stability in many cases. Recent work shows that interface structure has an outsize effect on the thermal behavior of nanostructured alloys. Here, this work focuses on the role of controlled heterophase interface structure in the thermal evolution of model Cu/Nb nanolaminates. We introduce 3D interfaces containing nanoscale heterogeneities in all spatial dimensions between Cu and Nb, forming 3D Cu/Nb. TEM, nanoindentation, and DSC are used in tandem to establish thermal stability and to identify shifts in microstructure as a function of static annealing temperature. 3D interfaces are shown to survive annealing to 300 °C for 1 hr., while 3D Cu/Nb microstructure evolves to form low-density and voided regions correlating to the onset of layer pinch-off between 500 and 600 °C annealing temperatures. A diffusivity- and vacancy energetics-based mechanism is developed to explain void formation driven by 3D interface degradation at elevated temperature.

3D interfaces↗
DOE OSTI2024

Unveiling the role of structural vacancies in Mn-based Prussian blue analogues for energy storage applications

Prussian blue analogues (PBAs) are promising cathode materials for monovalent- and multivalent-ion batteries due to their large framework structures. Nevertheless, the influence of lattice vacancies on their electrochemical performance has not been thoroughly clarified, hindering the further development of PBAs. Here we identify two types of functional vacancies, i.e., structural vacancies (SVs) and incidental vacancies (IVs), in manganese hexacyanoferrate (MnHCF) through synchrotron-based X-ray absorption spectroscopy and density functional theory calculations. Unlike structurally disordered IVs, the introduction of structurally ordered SVs promotes ion transport and reduces the interaction between host ions and the framework, enabling improved cycling and rate performance. The controllable adjustment of SVs in K-rich MnHCF is achieved through a co-reactant method. Furthermore, the partial introduction of SVs in K-rich MnHCF is demonstrated to favor both a milder structural evolution by alleviating the Jahn–Teller distortion of Mn 3+ and a stable dynamic process of interface reaction. As a result, this study unveils the potential importance of incorporating structural vacancies into MnHCF for advanced energy storage applications.

Gao, Chongwei↗
DOE OSTI2025

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Recent Progress in Solid‐State Lithium Batteries through Cathode Microstructure Engineering

A high‐performance cathode is necessary to realize the great potentials of solid‐state batteries such as high energy density and long cycle life. It is also needed to validate electrolyte performance, which is lacking. Currently, cathodes for solid‐state batteries are thinner and have lower cathode active material content than their lithium‐ion battery counterpart, resulting from insufficient conductivity and limiting the battery energy density. This review article provides an overview of recent development in cathode microstructures and their impact on battery properties, including compatibility between cathode and electrolyte, cathode architecture design, interface engineering, correlation between material properties, cathode processing approaches, and performance, as well as the advanced characterization methods used to understand the above correlations . Some perspectives on future development are shared including utilizing in situ and operando characterization tools to better understand dynamic evolution of the cathode/electrolyte interface, adapting artificial intelligence and machine learning to design and optimize cathode structures. The article is aimed to promote research interest on cathode development and advance solid‐state battery technologies.

cathode engineering↗
DOE OSTI2024

Assessment of effective elastic constants of U-10Mo fuel: A multiscale modeling and homogenization study

The significant microstructural changes that U-Mo fuel undergoes during operation degrades its mechanical properties and structural integrity. Microstructural evolution entails the formation, evolution, and redistribution of porosity in conjunction with grain refinement. In the present paper, we employ numerical approaches to assess the impact of the various microstructural features—grains, nanoscale intragranular fission gas bubbles, and mesoscale intergranular voids—on the degradation of elastic constants. Phase-field microstructure models are combined with the asymptotic expansion homogenization technique in order to derive the effective elastic constants as a function of porosity and fission density. Here the results are verified and compared against theoretical bounds. Using this approach, elastic degradation in operating nuclear fuels can be quantified when the distributions of microstructural features are known.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Preheat effects in laser-driven Rayleigh–Taylor instability experiments at intensities greater than $10^{15}$ $\textrm{W}$ $\textrm{cm}^{-2}$ at OMEGA EP and the NIF

The propagation of high-energy X-rays or hot electrons have the potential to alter the initial conditions in experimental target designs, especially at material interfaces, for laser-driven inertial confinement fusion (ICF) and high-energy density (HED) experimental platforms. Hot-electron preheat can drastically modify the initial conditions of experimental targets used to study the deceleration-stage Rayleigh–Taylor instability (RTI) both with and without applied magnetic fields. Therefore, it is necessary to understand and quantify the impact of hot-electron preheat. The hydrodynamic (HD) capabilities in the Ares code are used to study the effects varying levels of preheat can have on RTI evolution. The experimental and computational studies presented in this work demonstrate that at high laser intensities of around or greater than 10 15 W cm −2 , there is hot-electron generation from laser plasma instabilities which induces substantial preheat and impacts the morphology of RTI evolution and even inhibits the intended RTI growth such that it is not observable experimentally. The necessity of better quantifying hot-electron induced preheat and mitigating its impact on such high-intensity direct-drive laser experiments in the future is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tunable Collective Excitations in Epitaxial Perovskite Nickelates

The formation of plasmons through the collective excitation of charge density has generated intense discussions, offering insights to fundamental sciences and potential applications. While the underlying physical principles have been well-established, the effects of multibody interactions and orbital hybridization on plasmonic dynamics remain understudied. Here, in this work, we present the observation of conventional metallic and correlated plasmons in epitaxial La 1-x Sr x NiO 3 (LSNO) films with varying Sr doping concentrations (x = 0, 0.125, 0.25), unveiling their intriguing evolution. Unlike samples at other doping concentrations, the x = 0.125 intermediate doping sample does not exhibit the correlated plasmons despite showing high optical conductivity. Through experimental investigation using spectroscopic ellipsometry and X-ray absorption spectroscopy, that is further supported by theoretical calculations, the O2p-Ni3d orbital hybridization for x = 0.125 is found to be significantly enhanced, alongside a considerable weakening of its effective interaction comprising long-range Coulomb and variable interaction, U*. These factors account for the absence of correlated plasmons and the high optical conductivity observed in LSNO(0.125). Our findings highlight the significant impact of orbital hybridization on the electronic structures and the formation of quasiparticles in strongly correlated systems, opening new paths for plasmonic-based engineering research.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗
DOE OSTI2024

Influence of Annealing Temperature on the OER Activity of NiO(111) Nanosheets Prepared via Microwave and Solvothermal Synthesis Approaches

Earth-abundant transition metal oxides are promising alternatives to precious metal oxides as electrocatalysts for the oxygen evolution reaction (OER) and are intensively investigated for alkaline water electrolysis. OER electrocatalysis, like most other catalytic reactions, is surface-initiated, and the catalyst performance is fundamentally determined by the surface properties. Most transition metal oxide catalysts show OER activities that depend on the predominantly exposed crystal facets/surface structure. Therefore, the design of synthetic strategies to obtain the most active crystal facets is of significant research interest. In this work, rock salt NiO OER catalysts with (111) predominantly exposed facets were synthesized by a solvothermal (ST) method either heated under supercritical or microwave-assisted (MW) conditions. Particular emphasis was placed on the influence of the post annealing temperature on the structural configuration and OER activity to compare their catalytic performances. The as-prepared electrocatalysts are pure α-Ni hydroxides which were converted to rock salt NiO (111) nanosheets with hexagonal pores after heat treatment at different temperatures. The OER activity of the electrodes has been evaluated in 0.1 M KOH using geometric and intrinsic current densities via normalization by the disk area and BET area, respectively. The lowest overpotential at a geometric current density of 10 mA/cm 2 is found for samples pretreated by heating between 400 and 500 °C with a catalyst loading of 115 μg/cm 2 . Despite the very similar nature of the catalysts obtained from the two methods, the ST electrodes show a higher geometric and intrinsic current density for 500 °C pretreatment. The MW electrodes, however, achieve an optimal geometric current density for 400 °C pretreatment, while their intrinsic current density requires pretreatment over 600 °C. Interestingly, pretreated electrodes show consistently higher OER activity as compared to the poorly crystalline/less ordered hydroxide as-prepared electrocatalysts. Thus, our study highlights the importance of the synthesis method and pretreatment at an optimal temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rationalizing Acidic Oxygen Evolution Reaction over IrO 2 : Essential Role of Hydronium Cation

Abstract The development of active, stable, and more affordable electrocatalysts for acidic oxygen evolution reaction (OER) is of great importance for the practical application of electrolyzers and the advancement of renewable energy conversion technologies. Currently, IrO 2 is the only catalyst with high stability and activity, but a high cost. Further optimization of the catalyst is limited by the lack of understanding of catalytic behaviors at the acid‐IrO 2 interface. Here, in strong interaction with the experiment, we develop an explicit model based on grand‐canonical density function theory (GC‐DFT) calculations to describe acidic OER over IrO 2 . Compared to the explicit models reported previously, hydronium cations (H 3 O + ) are introduced at the electrochemical interface in the current model. As a result, a variation in stable IrO 2 surface configuration under the OER operating condition from previously proposed complete *O‐coverage to a mixture coverage of *OH and *O is revealed, which is well supported by in situ Raman measurements. In addition, the accuracy of predicted overpotential is increased in comparison with the experimentally measured. More importantly, an alteration of the potential limiting step from previously identified *O→*OOH to *OH→*O is observed, which opens new opportunities to advance the IrO 2 ‐based catalysts for acidic OER.

Mou, Tianyou↗
DOE OSTI2024

First Observation of Harmonics of Magnetosonic Waves in Martian Magnetosheath Region

Abstract The present study provides an evidence for the generation of harmonics of magnetosonic waves in the Martian magnetosheath region. The wave signatures are manifested in the magnetic field measurements recorded by the fluxgate magnetometer instrument onboard the Mars Atmosphere and Volatile Evolution missioN (MAVEN) spacecraft in the dawn sector around 5–10 LT at an altitude of 4,000–6,000 kms. The wave that is observed continuously from 19.1 to 20.7 UT below the proton cyclotron frequency (f ci ≈ 46 mHz) is identified as fundamental mode of the magnetosonic wave. Whereas harmonics of the magnetosonic wave are observed during 19.7–20.3 UT at frequencies that are multiple off ci . The ambient solar wind proton density and plasma flow velocity are found to vary with a fundamental mode frequency of 46 mHz. It is noticed that the fundamental mode is mainly associated with the left‐hand (LH), and higher frequency harmonics are associated with the right‐hand (RH) circular polarizations. A clear difference in the polarization and ellipticity is noticed during the time of occurrence of harmonics. The magnetosonic wave harmonics are found to propagate in the quasi‐perpendicular directions to the ambient magnetic field. The results of linear theory and Particle‐In‐Cell simulation performed here are in agreement with the observations. The present study provides a conclusive evidence for the occurrence of harmonics of magnetosonic wave in the close vicinity of the magnetosheath region of the unmagnetized planet Mars.

Astronomy & Astrophysics↗
DOE OSTI2025

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Axion Perturbations: A General Analytical Treatment

Cosmological data provides us two key constraints on dark matter (DM): it must have a particular abundance, and it must have an adiabatic spectrum of density perturbations in the early universe. Many different cosmological scenarios have been proposed that establish the abundance of axion DM in qualitatively different ways. In this paper we emphasize that, despite this variety of backgrounds, the perturbations in axion DM can be understood from universal principles. How does a feebly interacting axion field acquire perturbations proportional to those of photons? How do the isocurvature power spectrum and non-Gaussianity depend on the background evolution of the universe? We answer these questions for a completely general choice of cosmological background and temperature-dependent axion potential. We show that the most general solution to the axion field equation on super-horizon scales is entirely determined by the family of background solutions for different initial field values . This holds for both the component in the field perturbation solution contributing to the DM isocurvature perturbation (enhanced at late times by the sensitivity of the DM abundance to the initial condition, , which can be large for initial conditions near the hilltop), and the other component that contributes to the DM curvature perturbation. In particular, we explain that an unperturbed axion field in the early universe evolving into one with nontrivial adiabatic perturbations is guaranteed by Weinberg's theorem on adiabatic modes. These results have been derived before with various assumptions, such as a radiation dominated background or a quadratic potential. Our aim is to give a clear, simple derivation that is manifestly independent of those assumptions, and thus can be applied to any cosmological axion scenario.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗