Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probabilistic learning on manifold”

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.

2,087 records · Page 103
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Investigation Into LiCl-LiF+Li 2 O as a Low Temperature Electrolyte in Direct Electrolytic Reduction of UO 2

Direct electrolytic reduction (DER) of UO 2 in a binary LiCl-Li2O (1–3 wt%) at 650 °C has been widely reported. The process temperature can be reduced to 550 °C in theory by using a near eutectic ternary LiCl-LiF+Li 2O (∼2 wt%) electrolyte, which will result in substantial reduction in rate of salt vaporization. The feasibility of using the ternary electrolyte for the process was tested via material compatibility testing, electrochemical polarization testing, and controlled potential or current DER experiments. MgO, Pt, and 625 nickel samples were each contacted with the ternary salt for 168 h. Negligible changes to dimensions and mass of the samples were measured, and concentrations of corrosion products in the salt were 102 ppm or less. O2−oxidation was observed to occur at a lower potential than Cl-or F-on a Pt anode. DER of UO2was performed in six trials. Process variables included type of reference electrode, electrochemical control method, type of cathode basket, and electrolyte composition. UO 2conversion ranged from 28.2 to 99.9 % , and current efficiency ranged from 38 to 52%. Overall, results indicate that the use of the ternary salt is feasible for DER.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Weak Lensing Mass Calibration of the ACT DR5 Galaxy Clusters with the DES Year 3 Weak Lensing Data

We use weak gravitational lensing measurements from Year 3 Dark Energy Survey data to calibrate the masses of 443 galaxy clusters selected via the Sunyaev-Zel'dovich effect from Atacama Cosmology Telescope Data Release 5 maps of the cosmic microwave background. We incorporate redshift and SZ measurements for individual clusters into a hierarchical model for the stacked lensing signals and perform Bayesian analyses to constrain the hydrostatic mass bias of the clusters. Our treatment of systematic uncertainties includes a prescription for measuring and accounting for the weak lensing boost factor, consideration of a miscentering effect, as well as marginalization over uncertainties in the source galaxy photometric redshift distributions and shear calibration. The resultant constraints on the normalization of the mass-observable relation have a precision of approximately 7%, with the mean WL halo mass of M $_{500c}$ = 5.4 × 10$^{14}$ M $_{⊙}$. We measure the bias between the true cluster mass and the mass estimated from the SZ signal based on an X-ray-calibrated scaling relation assuming hydrostatic equilibrium, to be 1 - b = 0.74$^{+0.06}$ $_{-0.05}$ over the full sample. When splitting the clusters into high (z = 0.43-0.70) and low (z = 0.15-0.43) redshift bins, we measure 1 - b = 0.58$^{+0.06}$ $_{-0.05}$ and 0.81$^{+0.08}$ $_{-0.06}$, respectively. When introducing additional freedom in redshift and mass to the hydrostatic bias model, we find that 1 - b decreases with redshift (with the power law of -1.8$^{+0.5}$ $_{-0.6}$, 99.95% confidence), consistent with findings from other recent studies, while we do not find any significant trend in mass. We also demonstrate that our result is robust against various systematics such as a scale cut, priors on baryonic and miscentering parameters, and degree of scatter in mass-observable relation. The weak-lensing mass calibration presented in this study will be a useful tool for using the ACT clusters as probes of astrophysics, and as a step towards using their abundance as a cosmological probe.

Shin, T. [Carnegie Mellon U.] (ORCID:0000000263895↗
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

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
NASA NTRSDate not supplied

Thermal-Fluid Analysis of a Liquid-Cooled Battery Module for Electrified Aircraft

The development of safe, energy-dense batteries is critical to advancing hybrid electric and fully electrified aircraft propulsion. Achieving this capability requires a thermal management system that can maintain battery performance and safety under demanding operational conditions. The objective of this project is to support the maturation of next-generation lithium-ion batteries for electrified aircraft by conducting performance testing on integrated battery modules, specifically a 2-cell series configuration module housed within an aluminum enclosure. Designed to operate at a nominal 7.2 V with discharge rates up to 2.5C, the module will eventually be used to power an electric motor and DC-DC converter, generating substantial thermal loads that must be effectively managed to increase the usable energy and power density of electrified aircraft. To address these thermal challenges, this study presents the development and thermal-fluid analysis of a liquid-cooled thermal management system. While the full aircraft architecture utilizes an eight-string configuration, the present work evaluates a representative single-string water coolant loop to characterize baseline performance. The active cooling loop circulates water through a reservoir, pump, the battery module, and a variable area flow meter. Key performance metrics including component-level temperatures, mass flow rates, and pressure drops are quantified across the loop. Across discharge rates ranging from 0.5C to 2.5C, the active thermal loop consistently and effectively removed heat from the module, validating the design approach and confirming readiness for further development. The validated thermal performance indicates a path toward scalable battery modules that could enable energy and power dense systems for hybrid electric aircraft.

Electrified Aircraft↗
DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Heavy long-lived dark vector via a gluonic portal

We study a dark gauge boson Z′ that exclusively couples to the QCD gluons through higher dimensional operators. These operators are generated from integrating out of heavy ultraviolet resonances carrying both QCD and dark gauge charges. With SU(3)C gauge invariance, charge and parity symmetries preserved, we find that the leading effective operators are restricted to have the form of Z′GGG and Z′Z′GG at dimension-eight, which can naturally render the Z′ particle long-lived, and serve as a viable dark matter candidate. We investigate the phenomenology of these operators with both collider experiments and cosmological observation, without and with the assumption that this dark gauge boson plays the role of the dominant dark matter component. For an unstable Z′, we show that depending on its lifetime, it can be probed by various observables up to ultraviolet physics scale around 10 9 GeV. For Z′ being dark matter, we find that $m_{Z'}$ ≳ 1 TeV is consistent with the thermal freeze-out scenario. In contrast, in the freeze-in scenario, the extremely small couplings leave the relevant parameter space largely unconstrained by current experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Smart building HVAC control challenge: experience and solutions from the ADRENALIN project

A smart building HVAC control competition crowdsourced and compared algorithms on fair and equal ground using the standardized BOPTEST framework. The competition attracted 138 participants, but only 9% submitted valid solutions for the final stage, highlighting the complexity of advanced HVAC control design. The winning solutions showed significant potential to reduce energy use and cost by shifting demand, without compromising occupant comfort. Across scenarios, thermal energy cost reductions of 36–76% relative to a baseline, were achieved. In peak heat periods, the cost reduction leveraged limited energy use reduction (0–15%), but more significant energy price reduction (34–62%). This shows smart controls' ability to avoid as much as possible consumption during the morning peak hours, when spot prices are tendentially the highest. Hosting the competition has highlighted challenges in creating competitions that both are fair and promotes solutions that are transferable to real life implementation.

BOPTEST↗
DOE OSTI2026

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc↗
DOE OSTI2026

Peptoid-Based Nanosheets Exhibiting Broad Antiviral Activity Against Enveloped RNA Viruses

Enveloped RNA viruses, such as Influenza A (H1N1) and Sindbis virus, pose persistent global health threats due to their high mutation rates, efficient transmission, and frequent drug resistance. By mimicking host cell membrane receptors, multivalent virus inhibitors can block viral attachment, making them promising broad-spectrum antiviral agents. However, most of existing antivirals are often limited by strain specificity, short-lived efficacy, and toxicity. Here, we introduce a broad-spectrum antiviral platform based on highly tunable and biocompatible two-dimensional nanomembranes (2DNMs) self-assembled from amphiphilic peptoids, operating via a non-genomic, mutation-insensitive mechanism. By varying peptoid sequence, we design and synthesize over twenty different 2DNMs with various surface charge and high density of viral-attachment ligands (VALs). The self-assembled architecture of these stable 2DNMs provides cooperative noncovalent multivalent binding to virus particles that result in effective inhibition of viral infection. Screening of variants identified three leads that potently suppressed Influenza A (H1N1) and Sindbis virus infection across median tissue culture infectious dose (TCID50), plaque, RT–qPCR, and immunofluorescence assays, while maintaining >90% cell viability. These nanosheets significantly reduced infectious titers, viral RNA replication, and intracellular viral protein expression, indicating inhibition at early stages of viral entry and propagation. The sequence programmability, chemical robustness, and mutation-insensitive antiviral activity distinguish 2DNMs from traditional antivirals and positions them as a versatile materials platform for antiviral coatings, protective barriers, and prophylactic biomedical applications.

Influenza A virus↗
DOE OSTI2026

Preliminary Assessment for the Electric Load Shifting Potential of Integrating Thermal Energy Storage with Heat Pumps in Residential Buildings in Texas of United States

The widespread adoption of electric-driven heat pumps for heating and cooling is expected to significantly increase electric demand. Cooling electric demand will result in an increase in peak hours electric loads, placing additional strain on the grid during peak hours. This challenge, combined with the current rapidly increasing demand from data centers, exacerbates the electric demand duck curve problem. Integrating thermal energy storage (TES) with heat pump can shift electric use for heating and cooling from peak to off-peak hours of the electric grid, which can help flatten the daily electric demand profile of a building. This paper presents a novel design that integrates heat pump with TES (HP-TES), which uses phase change materials (PCM). Heat pump charges TES by melting or freezing PCM during off-peak hours when there is no thermal demand from the building. TES is then discharged (i.e., by freezing or melting PCM) during peak hours to provide a more favorable heat source or heat sink for the heat pump to meet the thermal demand of the building with lower electricity use than conventional air-source heat pumps. Computer simulations were developed to predict the performance of HP-TES applied to a typical single-family house in the US. The building-level simulation results were scaled up to preliminarily assess the aggregated impacts of deploying HP-TES across all single-family houses in Texas of the United States, including reduction of peak demand of the electric grid.

Anees, Fady [ORNL]↗
DOE OSTI2026

From Bricks to Clicks: Mapping the White Space in Building Innovation

It is a critical national imperative to transform the buildings sector, yet innovation is impeded by deployment failures that leave promising technologies stranded. Conventional market reports and techno-economic analysis provide an insufficient understanding of markets and resource allocation for emerging building technologies. They omit crucial commercialization factors such as ecosystem maturity and adoption friction, where the coordinated participation of a network of suppliers, contractors, financiers, regulators, and integrators is required to scale solutions. This study addresses these gaps by introducing an evaluation framework grounded in front-line data from six years of the DOE's IMPEL incubator, comprising experience from 300 building-sector innovators and the adjacent, complex ecosystem. Our methodology synthesizes top-down market analysis with bottom-up, practitioner-level data across five megatrends: (M1) Affordable materials and industrialized construction; (M2) Healthy and efficient mechanical systems; (M3) Intelligent building operations; (M4) Buildings as grid assets; and (M5) High-density power and cooling for data centers and therein identify twelve "white space" technology opportunities. Next, we develop a multi-criteria scoring rubric to rank these opportunities based on parameters, i.e., Affordability, Quality of Life, Reliability, and Security, yielding composite ‘Demand’ and ‘Maturity’ indices. Our results indicate that the most significant white spaces may not be incremental products but a new class of ‘Ecosystem Enablers’, such as logistics platforms, orchestration layers, and automated compliance software that solve structural deployment gaps. This paper summarizes this transparent, evidence-based, practitioner-informed evaluation framework for policymakers and investors to re-evaluate policy and resource allocation and unlock scalable market transformation.

Singh, Reshma↗
DOE OSTI2026

Fast Active-Set Thresholding Method for Nonnegative Least Squares

Nonnegative Least Squares (NNLS) is a fundamental constrained optimization problem encountered in many applications such as image deblurring, signal processing, nonnegative matrix factorization, magnetic microscopy, and hyperspectral imaging. Active-set based methods are a common class of algorithms for solving NNLS which identify the optimal variable set of the NNLS solution. They do so by iteratively solving a series of unconstrained least squares problems, identifying which variables violate the nonnegativity constraints, and then swapping variables in/out of consideration until the optimal set of variables is found. Several variations improving upon this method exist in the literature. In this work, we propose an active-set swap heuristic which further improves upon existing active-set based methods for NNLS. Our optimizations are based upon adding multiple variables to the passive set within a threshold of the smallest gradient value and removing variables within a similar threshold of the closest boundary constraint. We leverage these optimizations to yield a Fast Active-Set Thresholding NNLS (FAST-NNLS) algorithm which significantly outperforms the existing state-of-the-art NNLS algorithms for a wide range of problems. Rigorous convergence guarantees are proven for the proposed method. We demonstrate the effectiveness of our proposed method on multiple synthetic datasets and two realworld text analysis applications. In doing so, we present the most comprehensive NNLS solver comparison in the literature to date.

Cobb, Benjamin [Georgia Institute of Technology]↗
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

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

Carbon source–driven metabolic and regulatory remodeling defines phenomic states in Lipomyces starkeyi

Lipomyces is a genus of oleaginous yeasts with potential for contributing to reliable biomanufacturing supply chains. However, progress in advanced strain designs and engineering efforts are still constrained by a lack of understanding of the underlying molecular drivers of Lipomyces phenotypes. To address this gap, we collected a suite of multi-omic data to dissect how carbon source availability reshapes the metabolic network, lipid allocation, and regulatory architecture of Lipomyces starkeyi. We observed that glucose promotes biosynthetic and proliferative processes supported by abundant energy and carbon intermediates, xylose enhances redox-balancing mechanisms centered on the pentose phosphate pathway, and glycerol activates respiratory metabolism, ß-oxidation, and the glyoxylate cycle. Lipid species distributions remained consistent in both nitrogen replete and depleted conditions across the carbon sources, indicating robust production mechanisms. Regulatory protein identification and network analysis revealed glycerol-driven respiratory growth favors regulatory programs integrating stress tolerance, redox balance, and lipid-associated metabolism, whereas xylose growth activates compensatory transcriptional responses aimed at maintaining mitochondrial function. Nitrogen limitation modulates the strength of these responses but does not fundamentally alter their direction, reinforcing carbon source as the dominant driver of regulatory architecture. Taken together, this data enhances the understanding of Lipomyces molecular rearrangements and provides a foundation for further development of predictive phenotypic tools in this genus.

Biotechnology↗
DOE OSTI2026

Boron-Based Neutron Scintillator Screen Characterization with X-Rays and Neutrons

Recent work on boron-based neutron scintillator screens suggests these screens can offer superior performance when compared to commonly used screens. Borated neutron scintillator screens perform well in terms of light output (5-6 times greater than a standard Gadox screen) and detection effi-ciency (larger than standard LiF+ZnS screens). However, previously manu-factured boron-based screens have exhibited non-uniform surface coating and a poor mixture between phosphor and converter particles. The objective of this work was to evaluate newly fabricated scintillator screens to deter-mine if enhanced fabrication methods produced a more homogeneous distribution between neutron converter and scintillation phosphor particles. Uniformity of scintillator material deposition was also inspected. This new iteration of screens appeared more uniform than previous generations with the new coating method improving surface chemistry and scintillator material homogeneity. Additionally, a new methodology for screen characterization, involving the correlation of a neutron image taken with a borated scintillator screen to X-ray computed tomography of that same screen, was demonstrated to elucidate a relationship between scintillator screen thickness and relative light output of the screen under neutron exposure. This method suggested that the ideal thickness of scintillator material was ~150 µm to maximize light output of the screen.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Laser-induced graphene gas sensors for environmental monitoring

Artemesia tridentatais a foundational plant taxon in western North America and an important medicinal plant threatened by climate change. Low-cost fabrication of sensors is critical for developing large-area sensor networks for understanding and monitoring a range of environmental conditions. However, the availability of materials and manufacturing processes is still in the early stages, limiting the capacity to develop cost-effective sensors at a large scale. In this study, we demonstrate the fabrication of low-cost flexible sensors using laser-induced graphene (LIG); a graphitic material synthesized using a 450-nm wavelength bench top laser patterned onto polyimide substrates. We demonstrate the effect of the intensity and focus of the incident beam on the morphology and electrical properties of the synthesized material. Raman analyses of the synthesized LIG show a defect-rich graphene with a crystallite size in the tens of nanometers. This shows that the high level of disorder within the LIG structure, along with the porous nature of the material provide a good surface for gas adsorption. The initial characterization of the material has shown an analyte response represented by a change in resistance of up to 5% in the presence of volatile organic compounds (VOCs) that are emitted and detected byArtemisiaspecies. Bend testing up to 100 cycles provides evidence that these sensors will remain resilient when deployed across the landscapes to assess VOC signaling in plant communities. The versatile low-cost laser writing technique highlights the promise of low-cost and scalable fabrication of LIG sensors for gas sensor monitoring.

Chemistry↗
DOE OSTI2024

Understanding the origin of early-type dwarfs: the spectrophotometric study of CGCG014−074

ABSTRACT Early-type dwarf galaxies constitute a prevalent population in the central regions of rich groups and clusters in the local Universe. These low-luminosity and low-mass stellar systems play a fundamental role in the assembly of the luminous galaxies observed today, according to the Lambda cold dark matter hierarchical theory. The origin of early-type dwarfs has been linked to the transformation of disc galaxies interacting with the intracluster medium, especially in dense environments. However, the existence of low-luminosity early-type galaxies in low-density environments presents a challenge to this scenario. This study presents a comprehensive photometric and spectroscopic analysis of the early-type dwarf galaxy CGCG014−074 using deep Gemini GMOS (Gemini Multi-Object Spectrograph) data, focusing on its peculiarities and evolutionary implications. CGCG014−074 exhibits distinct features, including a rotating inner disc, an extended stellar formation with a quiescent phase since about 2 Gyr ago, and the presence of boxy isophotes. From the kinematic analysis, we confirm CGCG014−074 as a nucleated early-type dwarf galaxy with embedded disc. The study of its stellar population parameters using different methods provides significant insights into the galaxy’s evolutionary history. These results show an old and metal-poor nucleus (${\sim}9.3$ Gyr and $\mathrm{[Z/H]}\sim -0.84$ dex), while the stellar disc is younger (${\sim}4.4$ Gyr) with a higher metallicity ($\mathrm{[Z/H]}\sim -0.40$ dex). These distinctive features collectively position CGCG014−074 as a likely building block galaxy that has evolved passively throughout its history.

Astronomy & Astrophysics↗