Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “big”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Impact of vertical and seasonal variation in leaf traits on simulating soybean canopy photosynthesis via 1D and 3D modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax 25 ), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax 25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax 25 profile achieved comparable performance to the scheme with the measured Vcmax 25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R 2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R 2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effective direct steam regeneration of bis-iminoguanidine solid sorbent used for carbon dioxide capture

A cost-effective, energy-efficient sorbent regeneration process for phase-changing guanidines used for CO 2 capture was developed based on direct-steam stripping. This approach enhances the regeneration rate, simplifies the overall CO 2 capture process, and reduces the energy cost compared to conventional conductive thermal regeneration. A direct-steam sorbent regeneration reactor was developed, demonstrating that aqueous bis(iminoguanidines) (BIG) sorbents, e.g., methylglyoxal-bis(iminoguanidine) (MGBIG) and glyoxal-bis(iminoguanidine) (GBIG), could be efficiently regenerated with up to ~ 99 % CO 2 recovery through direct-steam stripping. Using low-temperature steam at 100 °C, a 4.5 times faster regeneration rate for GBIG carbonate sorbent (e.g., 30 min for 10 g) was demonstrated compared to conductive-heating (e.g., 135 min for 10 g) at 130 °C. Additionally, fully regenerated MGBIG converts into an aqueous MGBIG solution when the steam condenses onto the sorbent surface. Condensed steam with the guanidine can be easily recycled as an aqueous solution into the gas–liquid contactor to achieve a continuous-flow CO 2 -capture process. Molecular dynamics simulation was employed to provide a better understanding of the process. Higher heat transfer rates from steam to guanidine carbonate, compared to air heating, were attributed to the vibration resonance of water molecules within MGBIG with that of vapor molecules and the effective transfer of kinetic energy from vapor to solid. Technoeconomic analysis demonstrated that direct-steam stripping significantly decreases the CO 2 capture cost by 50 % compared to traditional conductive heating methods. Further, enhanced mass transfer facilitated by low-temperature steam and subsequent condensation effectively heats up the H 2 O-containing BIG-carbonate crystals, facilitating the desorption of CO 2 from the solid crystals, thereby leading to fast, effective, and energy-efficient sorbent regeneration.

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

Dark matter is the new BBN

Measurements of the primordial element abundances provide us with an important probe of our universe’s early thermal history, allowing us to constrain the expansion rate and composition of our universe as early as ~ 1 s after the Big Bang. Prior to this time, we have essentially no empirical information on which to base any such claims. In this paper, we imagine a future time in which we have not only detected the particles that make up the dark matter, but have measured their mass and annihilation cross section with reasonable precision. In analogy to the light element abundances, the dark matter abundance in this scenario could be used to study and constrain the expansion rate and composition of our universe at the time of dark matter freeze out, which for a standard thermal relic occurs at T f ~ m χ /20, corresponding to t ~ 4 x 10 -10 s x (TeV/m χ ) 2 , many orders of magnitude prior to the onset of Big Bang nucleosynthesis. Here we consider how such measurements could be used to constrain scenarios which feature exotic forms of radiation or matter, a ultralight scalar, or modifications to gravity, each of which have the potential to be much more powerfully probed with dark matter than with the light element abundances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Unsupported and carbon-supported silver catalysts for oxygen reduction reaction in alkaline media

Quick and easy Ag catalysts preparation via wet chemical synthesis method using only reducing agent (pure-Ag); reducing agent and citric acid as the capping agent (Ag-CA); and carbon support (KetjenBlack 600J), capping agent, and the reducing agent (Ag/C) is demonstrated. The Ag-based electrocatalysts are characterized by high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) with energy-dispersive X-ray spectroscopy (EDS), scanning electron microscopy (SEM), X-ray diffraction (XRD) analysis, and X-ray photoelectron spectroscopy (XPS). The electrocatalytic activity of Ag catalysts for O 2 reduction reaction (ORR) in 1 M KOH is evaluated using the rotating (ring)-disc electrode method. SEM and HAADF-STEM results show that the unsupported pure-Ag and Ag-CA catalysts consist mainly of big agglomerates, and Ag/C has the smallest agglomerates and some sub-3 nm Ag nanoparticles. The XPS results reveal that Ag in all the catalysts is in the metallic form (Ag 0 ). Despite consisting of big agglomerates, the Ag-CA catalyst exhibits similar ORR electrocatalytic activity to that of Ag/C. Ag-CA (unsupported) shows the lowest hydrogen peroxide yield. These results are of great importance for the development of Ag-based catalysts that can be prepared in a fast, simple and easily up scalable fashion, for anion exchange membrane fuel cells.

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

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H 0 , the σ 8 – S 8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0 σ tension between the Planck CMB estimate of the Hubble constant H 0 and the SH0ES collaboration measurements. After showing the H 0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω m , and the amplitude or rate of the growth of structure ( σ 8 , f σ 8 ). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H 0 – S 8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Metal Halide Perovskites in Quantum Dot Solar Cells: Progress and Prospects

“Go big or go home” could never be a truer statement, especially when it comes to energy. The world population is increasing, the energy use per person is growing more rapidly, and the total electricity use per person is growing even more quickly than that. To handle this demand, energy production must be increased, and it is vital for renewable sources to be used. Currently three quarters of a terawatt of power generated from photovoltaics (PVs) has been installed. At a somewhat optimistic average power-conversion efficiency of 18%, the distributed surface area of high purity optoelectronic grade semiconductor photovoltaic panels pointed at the sun is substantially larger than the area of the whole state of Rhode Island or nearly double the land mass of Luxembourg. While photovoltaic production continues to ramp to immense proportions, perovskite semiconductors are poised to greatly complement existing technology. On the other end of the size spectrum, the study of semiconductor nanocrystals or quantum dots (QDs) has led to advanced structures and deeper understanding within halide perovskite semiconductors. In this article, we show how the development of nanoscale metal halide perovskite semiconductors have gained prominence surpassing all other QD materials in terms of efficiency, and are becoming a platform for further improving technology to solve big energy challenges.

14 SOLAR ENERGY↗
DOE OSTI2022

Scalable algorithms for physics-informed neural and graph networks

Physics-informed machine learning (PIML) has emerged as a promising new approach for simulating complex physical and biological systems that are governed by complex multiscale processes for which some data are also available. In some instances, the objective is to discover part of the hidden physics from the available data, and PIML has been shown to be particularly effective for such problems for which conventional methods may fail. Unlike commercial machine learning where training of deep neural networks requires big data, in PIML big data are not available. Instead, we can train such networks from additional information obtained by employing the physical laws and evaluating them at random points in the space–time domain. Such PIML integrates multimodality and multifidelity data with mathematical models, and implements them using neural networks or graph networks. Here, we review some of the prevailing trends in embedding physics into machine learning, using physics-informed neural networks (PINNs) based primarily on feed-forward neural networks and automatic differentiation. For more complex systems or systems of systems and unstructured data, graph neural networks (GNNs) present some distinct advantages, and here we review how physics-informed learning can be accomplished with GNNs based on graph exterior calculus to construct differential operators; we refer to these architectures as physics-informed graph networks (PIGNs). We present representative examples for both forward and inverse problems and discuss what advances are needed to scale up PINNs, PIGNs and more broadly GNNs for large-scale engineering problems.

42 ENGINEERING↗
DOE OSTI2023

Near Quantitative Removal of Selenate and Sulfate Anions from Wastewaters by Cocrystallization with Chelating Hydrogen-Bonding Guanidinium Ligands

Selenium (Se) has become an environmental contaminant of aquatic ecosystems as a result of human activities, particularly mining, fossil fuel combustion, and agricultural activities. By leveraging the high sulfate concentrations relative to Se oxyanions (i.e., SeO n 2– , n = 3, 4) present in some wastewaters, we have developed an efficient approach to Se-oxyanion removal by cocrystallization with bisiminoguanidinium (BIG) ligands that form crystalline sulfate/selenate solid solutions. The crystallization of the sulfate, selenate and selenite, oxyanions and of sulfate/selenate mixtures with five candidate BIG ligands are reported along with the thermodynamics of crystallization and aqueous solubilities. Oxyanion removal experiments with the top two performing candidate ligands show a near quantitative removal (>99%) of sulfate or selenate from solution. When both sulfate and selenate are present, there is near quantitative removal (>99%) of selenate, down to sub-ppb Se levels, with no discrimination between the two oxyanions during cocrystallization. Reducing the selenate concentrations by 3 orders of magnitude or more relative to sulfate, as found in many wastewaters, led to no measurable loss in Se removal efficiencies. This work offers a simple and effective alternative to selective separation of trace amounts of highly toxic selenate oxyanions from wastewaters, to meet stringent regulatory discharge limits.

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

A Vertically Resolved Canopy Improves Chemical Transport Model Predictions of Ozone Deposition to North Temperate Forests

Abstract Dry deposition is the second largest tropospheric ozone (O 3 ) sink and occurs through stomatal and nonstomatal pathways. Current O 3 uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O 3 fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O 3 fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O 3 deposition and its drivers. We show that GEOS‐Chem overpredicts observed O 3 fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O 3 deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O 3 flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O 3 deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O 3 fluxes; and (c) the CTM big‐leaf treatment overestimates O 3 ‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O 3 predictions.

Vermeuel, Michael P. [Department of Soil, Water, a↗
DOE OSTI2024

Disturbance of hibernating bats due to researchers entering caves to conduct hibernacula surveys

Estimating population changes of bats is important for their conservation. Population estimates of hibernating bats are often calculated by researchers entering hibernacula to count bats; however, the disturbance caused by these surveys can cause bats to arouse unnaturally, fly, and lose body mass. We conducted 17 hibernacula surveys in 9 caves from 2013 to 2018 and used acoustic detectors to document cave-exiting bats the night following our surveys. We predicted that cave-exiting flights (i.e., bats flying out and then back into caves) of Townsend’s big-eared bats (Corynorhinus townsendii) and western small-footed myotis (Myotis ciliolabrum) would be higher the night following hibernacula surveys than on nights following no surveys. Those two species, however, did not fly out of caves more than predicted the night following 82% of surveys. Nonetheless, the activity of bats flying out of caves following surveys was related to a disturbance factor (i.e., number of researchers × total time in a cave). We produced a parsimonious model for predicting the probability of Townsend’s big-eared bats flying out of caves as a function of disturbance factor and ambient temperature. That model can be used to help biologists plan for the number of researchers, and the length of time those individuals are in a cave to minimize disturbing bats.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The structures of ordered defects in thiocyanate analogues of Prussian Blue

We report the structures of six new divalent transition metal hexathiocyanatobismuthate frameworks with the generic formula $\Big\{ M[Bi(SCN)_6]_\frac{2+x}{3} \Big\}$ $\{guest\}^{x+}$, M = Mn, Co, Ni and Zn. These frameworks are defective analogues of the perovskite-derived trivalent transition metal hexathiocyanatobismuthates M III [Bi(SCN) 6 ]. The defects in these new thiocyanate frameworks order and produce complex superstructures due to the low symmetry of the parent structure, in contrast to the related and more well-studied cyanide Prussian Blue analogues. Despite the close similarities in the chemistries of these four transition metal cations, we find that each framework contains a different mechanism for accommodating the lowered transition metal charge, making use of some combination of Bi(SCN) 6 3- vacancies, M Bi antisite defects, water substitution for thiocyanate, adventitious extra-framework cations and reduced metal coordination number. These materials provide an unusually clear view of defects in molecular framework materials and their variety suggests that similar richness may be waiting to be uncovered in other hybrid perovskite frameworks.

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

Rapid early gas accretion for the inner Galactic disc: A case for short accretion timescale

Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion timescale, this has not been directly constrained using observations. Using data for the inner regions of the Milky Way recently produced by APOGEE and Gaia and of unprecedented quantity and quality, we aim to derive strong constraints on the infall history of the inner (< 6 kpc) Galaxy (with a focus on stars between 4 and 6 kpc, which we show is an appropriate proxy for the entire inner disc). We implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt–Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time–[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius of between 4 and 6 kpc from APOGEE. Our results point to a fast, early gas accretion, and an upper limit on the accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2–3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Non-oscillatory no-scale inflation

We propose a non-oscillatory no-scale supergravity model of inflation (NO-NO inflation) in which the inflaton does not oscillate at the end of the inflationary era. Instead, the Universe is then dominated by the inflaton kinetic energy density (kination). During the transition from inflation to kination, the Universe preheats instantly through a coupling to Higgs-like fields. These rapidly annihilate and scatter into ultra-relativistic matter particles, which subsequently dominate the energy density, and reheating occurs at a temperature far above that of Big Bang Nucleosynthesis. After the electroweak transition, the inflaton enters a tracking phase as in some models of quintessential inflation. The model predictions for cosmic microwave background observables are consistent with Planck 2018 data, and the density of gravitational waves is below the upper bound from Big Bang Nucleosynthesis. Furthermore, we also find that the density of supersymmetric cold dark matter produced by gravitino decay is consistent with Planck 2018 data over the expected range of supersymmetric particle masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Time variability of ultra fast BAL outflows using SALT: C IV equivalent width analysis

We study the time variability (over ≤7.3 yr) of ultra fast outflows (UFOs) detected in a sample of 64 C IV broad absorption line (BAL) quasars (with 80 distinct BAL components) monitored using the Southern African Large Telescope. By comparing the properties of the quasar in our sample with those of a control sample of non-BAL quasars, we show that the distributions of black hole mass are different and the bolometric luminosities and optical photometric variations of UFO BAL quasars are slightly smaller compared to that of non-BAL quasars. The detection fraction of C IV equivalent width (W) variability (∼95 per cent), the fractional variability amplitude $\left(\frac{\Delta W}{W}\right)$ and the fraction of ‘highly variable’ BAL (i.e. $\big|$$\frac{\Delta W}{W}$$\big|$ > 0.67) components (∼33 per cent) are higher in our sample compared to the general BAL population. The scatter in $\frac{\Delta W}{W}$ and the fraction of ‘highly variable’ BALs increase with the time-scale probed. The $\frac{\Delta W}{W}$ distribution is asymmetric at large time scales. We attribute this to the BAL strengthening time-scales being shorter than the weakening time-scales. The BAL variability amplitude correlates strongly with the BAL properties compared to the quasar properties. BALs with low W, high-velocity, shallow profiles, and low-velocity width tend to show more variability. When multiple BAL components are present, a correlated variability is seen between low- and high-velocity components with the latter showing a larger amplitude variations. We find an anticorrelation between the fractional variations in the continuum flux and W. While this suggests photoionization induced variability, the scatter in continuum flux is much smaller than that of W.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Clustering and runaway merging in a primordial black hole dominated universe

If primordial black holes (PBHs) are present in the early Universe, their contribution to the energy budget grows relative to that of radiation and generically becomes dominant unless the initial abundance is exponentially small. This black hole domination scenario is largely unconstrained for PBHs with masses ≲ 10 9 g, which evaporate prior to big bang nucleosynthesis. However, if the era of PBH domination is sufficiently long, the PBHs form clusters and can merge appreciably within these objects. We calculate the population statistics of these clusters within the Press-Schechter formalism and find that, for a wide range of PBH masses and Hubble rates at the onset of PBH domination, the mergers within PBH clusters can exhibit runaway behavior, where the majority of the cluster will eventually form a single black hole with a mass much greater than the original PBH mass. Finally, these mergers can dramatically alter the PBH mass distribution and leave behind merged relic black holes that evaporate after big bang nucleosynthesis and yield various observational signatures, excluding parameter choices previously thought to be viable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Evolution of decaying particles and decay products in various scenarios for the future expansion of the Universe

We examine nonrelativistic particles that decay into relativistic products in big rip, little rip, and pseudorip models for the future evolution of the Universe. In contrast to decays that occur in standard ΛCDM, the evolution of the ratio r of the energy density of the relativistic decay products to the energy density of the initially decaying particles can decrease with time in all of these models. In big rip and little rip models, r always goes to zero asymptotically, while this ratio evolves to infinity or a constant in pseudorip models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Signatures of primordial energy injection from axion strings

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultralight axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a subdominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants f a ≲ 10 12 GeV , up to model dependence on the ultraviolet completion, for axion masses m a ≲ 10 − 29 eV ; future CMB surveys could find striking evidence of axion strings with lower decay constants. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulations of multicomponent relativistic thermalization

Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of the particles and fluids in such a mixture has been a challenge due to the fact that such simulations are computationally expensive in three spatial dimensions. Here, we report on the development and application of a multidimensional relativistic Monte Carlo code to explore the thermalization process in a relativistic multicomponent environment in a computationally inexpensive way. As an illustration we simulate the fully relativistic three-dimensional Brownian-motion-like solution to the thermalization of a high-mass particle (proton) in a bath of relativistic low-mass particles (electrons). Here, we follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during big-bang nucleosynthesis. We also simulate the thermalization of energetic particles injected into the plasma as can occur, for example, by the decay of massive unstable particles during the big bang.

79 ASTRONOMY AND ASTROPHYSICS↗