Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Explosive”

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 1,351 records · Page 75

Flavor triangle of the diffuse supernova neutrino background

Although Galactic core-collapse supernovae (SNe) only happen a few times per century, every hour a vast number of explosions happen in the whole universe, emitting energy in the form of neutrinos, resulting in the diffuse supernova neutrino background (DSNB). The DSNB has not yet been detected, but Super-Kamiokande doped with gadolinium is expected to yield the first statistically significant observation within the next several years. Since the neutrinos produced at the core collapse undergo mixing during their propagation to Earth, the flavor content at detection is a test of oscillation physics. In this paper, we estimate the expected DSNB data at the DUNE, Hyper-K and JUNO experiments which when combined are sensitive to all different neutrino flavors. Here, we determine how well the flavor content of the DSNB will be reconstructed in the future, for a Mikheyev-Smirnov-Wolfenstein (MSW) scenario as well as a neutrino decay scenario. A large fraction of the flavor space will be excluded, but the heavy-lepton neutrino flux remains a challenge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗

Successful νp -process in neutrino-driven outflows in core-collapse supernovae

The origin of the solar system abundances of several proton-rich isotopes, especially 92,94 Mo and 96,98 Ru, has been an enduring mystery in nuclear astrophysics. An attractive proposal to solve this problem is the νp-process, which can operate in neutrino-driven outflows in a core-collapse supernova after the shock is launched. Years of detailed studies, however, have cast doubt over the ability of this process to generate sufficiently high absolute and relative amounts of various p-nuclei. The νp-process is also thought to be excluded by arguments based on the long-lived radionuclide 92 Nb. Here, we present explicit calculations, in which both the abundance ratios and the absolute yields of the p-nuclei up to A ≲ 105 are successfully reproduced, even when using the modern (medium enhanced) triple-α reaction rates. The process is also shown to produce the necessary amounts of 92 Nb. The models are characterized by subsonic outflows and by the protoneutron star masses in the ≳ 1.7 M ⊙ range. In conclusion, this suggests that the Mo and Ru p-nuclides observed in the Solar System were made in CCSN explosions characterized by an extended accretion stage.

79 ASTRONOMY AND ASTROPHYSICS↗

Keck Infrared Transient Survey. I. Survey Description and Data Release 1

Abstract We present the Keck Infrared Transient Survey, a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roman SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. Here, we describe the first release of data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients <17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift z < 0.01 ( D ≈ 50 Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to z = 0.005. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events, luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using PypeIt , which requires minimal human interaction to ensure reproducibility.

Tinyanont, S. (ORCID:0000000214814676)↗

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of runaway electron impact on instrumented sacrificial limiters on DIII-D

Instrumented sacrificial limiter heads, both domed (proud) and flat (flush) are used in DIII-D runaway electron (RE) wall strikes to study the wall impact dynamics with high spatial and time resolution. The approximate structure of the RE wetted area and heating depth on the domed limiter heads were predicted qualitatively using orbit-tracking simulations, although a strong left–right asymmetry (about the magnetic field direction) was not captured well by the simulations. It is hypothesized that this difference is perhaps due to the local 3D magnetic field perturbation of the dome limiter head. The average kinetic energy K and pitch angle θ of REs striking the limiter head were estimated from the spatial distribution of local HXR emission and were estimated to be roughly K ≈ 4 MeV and θ ≈ 0.2 rad. These values are roughly consistent with in-plasma values estimated before the loss event, indicating that RE kinetic energy and pitch angle are not drastically altered when transporting to the wall. Large shot–shot variations (1–10 kJ) in energy deposition into the limiter head were observed and were explained by shot–shot variations in locked magneto-hydrodynamics mode toroidal phase. For the largest deposited energies (10 kJ), graphite material failure and explosive dust release was observed, and the depth of material failure at higher energy deposition was successfully reproduced using modelling of volumetric energy deposition and brittle failure. The presence of energetic (keV) level ion impact during the RE wall strike was confirmed by three different surface analysis techniques. The ratio of energetic ion to RE flux appears to be larger on flat surfaces within the RE wetted area, although the energetic ion flux and total energy flux due to energetic ions have not yet been quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear saturation of ballooning modes in stellarators

Ballooning mode saturation is investigated in realistic stellarator configurations using the flux tube approach of Ham et al (2018 Plasma Phys. Control. Fusion 60 075017), Ham et al (2016 Phys. Rev. Lett. 116 235001). The method is adapted to account for the lack of exact force balance in stellarator equilibrium solvers that assume existence of nested flux surfaces. A variational approach for calculating flux tube energy is developed to overcome this force error problem in stellarator numerical equilibria. Saturated (equilibrium) flux tube states that cross 10%–20% of the plasma minor radius are shown to exist for linearly ballooning unstable profiles. It is shown that several features of the displaced flux tube structure in a full nonlinear MHD simulation of Wendelstein 7X are reproduced by our model. Saturated states are found in a compact stellarator equilibrium close but below the marginal ballooning linear instability, i.e. the unperturbed equilibrium is metastable. This suggests that edge-localized-mode-like explosive MHD behavior may be possible in stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a high-resolution Rayleigh-Taylor experiment with the Crystal Backlighter Imager on the National Ignition Facility

The Rayleigh-Taylor (RT) instability affects a vast range of High Energy Density (HED) length scales, spanning from supernova explosions (10 13 m) to inertial confinement fusion (10 -6 m). In inertial confinement fusion, the RT instability is known to induce mixing or turbulent transition, which in turn cools the hot spot and hinders ignition. The fine-scale features of the RT instability, which are difficult to image in HED physics, may help determine if the system is mixing or is transitioning to turbulence. Earlier diagnostics lacked the spatial and temporal resolution necessary to diagnose the dynamics that occur along the RT structure. A recently developed diagnostic, the Crystal Backlighter Imager (CBI), [1,2] can now produce an x-ray radiograph capable of resolving the fine-scale features expected in these RT unstable systems. This paper describes an experimental design that adapts a well-characterized National Ignition Facility (NIF) platform to accommodate the CBI diagnostic. Simulations and synthetic radiographs highlight the resolution capabilities of the CBI in comparison to previous diagnostics. Furthermore, the improved resolution of the system can provide new observations to study the RT instability's involvement in mixing and the transition to turbulence in the HED regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of a liquid Ga target for Laser Ion Source

We experimentally proved the feasibility of a liquid-based target for Laser Ion Source (LIS) application. The target consists of melted metal gallium contained in a heated crucible. Ions current resulting from laser irradiation. Moreover, given the explosive ablation mechanism involved, each laser shot has been found to induce fluctuations in the surface level. It takes about 0.7 sec to recover the steady condition fully. Despite the maximum fluctuation being 2.7 mm (top-bottom maximum displacement), these fluctuations have shown no significant influence on total ion current and are independent of the temperature of the sample within the tested repetition rate. This study provides valuable insights into the potential of employing such a system for LIS.

43 PARTICLE ACCELERATORS↗

PRODeepSyn: predicting anticancer synergistic drug combinations by embedding cell lines with protein–protein interaction network

Abstract Although drug combinations in cancer treatment appear to be a promising therapeutic strategy with respect to monotherapy, it is arduous to discover new synergistic drug combinations due to the combinatorial explosion. Deep learning technology holds immense promise for better prediction of in vitro synergistic drug combinations for certain cell lines. In methods applying such technology, omics data are widely adopted to construct cell line features. However, biological network data are rarely considered yet, which is worthy of in-depth study. In this study, we propose a novel deep learning method, termed PRODeepSyn, for predicting anticancer synergistic drug combinations. By leveraging the Graph Convolutional Network, PRODeepSyn integrates the protein–protein interaction (PPI) network with omics data to construct low-dimensional dense embeddings for cell lines. PRODeepSyn then builds a deep neural network with the Batch Normalization mechanism to predict synergy scores using the cell line embeddings and drug features. PRODeepSyn achieves the lowest root mean square error of 15.08 and the highest Pearson correlation coefficient of 0.75, outperforming two deep learning methods and four machine learning methods. On the classification task, PRODeepSyn achieves an area under the receiver operator characteristics curve of 0.90, an area under the precision–recall curve of 0.63 and a Cohen’s Kappa of 0.53. In the ablation study, we find that using the multi-omics data and the integrated PPI network’s information both can improve the prediction results. Additionally, the case study demonstrates the consistency between PRODeepSyn and previous studies.

Wang, Xiaowen↗

Standardized naming of microbiome samples in Genomes OnLine Database

The power of next-generation sequencing has resulted in an explosive growth in the number of projects aiming to understand the metagenomic diversity of complex microbial environments. The interdisciplinary nature of this microbiome research community, along with the absence of reporting standards for microbiome data and samples, poses a significant challenge for follow-up studies. Commonly used names of metagenomes and metatranscriptomes in public databases currently lack the essential information necessary to accurately describe and classify the underlying samples, which makes a comparative analysis difficult to conduct and often results in misclassified sequences in data repositories. The Genomes OnLine Database (GOLD) (https:// gold.jgi.doe.gov/) at the Department of Energy Joint Genome Institute has been at the forefront of addressing this challenge by developing a standardized nomenclature system for naming microbiome samples. GOLD, currently in its twenty-fifth anniversary, continues to enrich the research community with hundreds of thousands of metagenomes and metatranscriptomes with well-curated and easy-to-understand names. Through this manuscript, we describe the overall naming process that can be easily adopted by researchers worldwide. Additionally, we propose the use of this naming system as a best practice for the scientific community to facilitate better interoperability and reusability of microbiome data.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluating the location capabilities of a regional infrasonic network in Utah, US, using both ray tracing-derived and empirical-derived celerity-range and backazimuth models

SUMMARY More realistic models for infrasound signal propagation across a region can be used to improve the precision and accuracy of spatial and temporal source localization estimates. Motivated by incomplete infrasound event bulletins in the Western US, the location capabilities of a regional infrasonic network of stations located between 84–458 km from the Utah Test and Training Range, Utah, USA, is assessed using a series of near-surface explosive events with complementary ground truth (GT) information. Signal arrival times and backazimuth estimates are determined with an automatic F-statistic based signal detector and manually refined by an analyst. This study represents the first application of three distinct celerity-range and backazimuth models to an extensive suite of realistic signal detections for event location purposes. A singular celerity and backazimuth deviation model was previously constructed using ray tracing analysis based on an extensive archive of historical atmospheric specifications and is applied within this study to test location capabilities. Similarly, a set of multivariate, season and location specific models for celerity and backazimuth are compared to an empirical model that depends on the observations across the infrasound network and the GT events, which accounts for atmospheric propagation variations from source to receiver. Discrepancies between observed and predicted signal celerities result in locations with poor accuracy. Application of the empirical model improves both spatial localization precision and accuracy; all but one location estimates retain the true GT location within the 90 per cent confidence bounds. Average mislocation of the events is 15.49 km and average 90 per cent error ellipse areas are 4141 km2. The empirical model additionally reduces origin time residuals; origin time residuals from the other location models are in excess of 160 s while residuals produced with the empirical model are within 30 s of the true origin time. We demonstrate that event location accuracy is driven by a combination of signal propagation model and the azimuthal gap of detecting stations. A direct relationship between mislocation, error ellipse area and increased station azimuthal gaps indicate that for sparse networks, detection backazimuths may drive location biases over traveltime estimates.

58 GEOSCIENCES↗

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗

A low-luminosity core-collapse supernova very similar to SN 2005cs

ABSTRACT We present observations and analysis of PSN J17292918+7542390, a low-luminosity Type II-P supernova (LL SN IIP). The observed sample of such events is still low, and their nature is still under debate. Such SNe are similar to SN 2005cs, a well-observed LL Type II-P event, having low expansion velocities, and small ejected 56Ni mass. We have developed a robust and relatively fast Monte Carlo code that fits semi-analytic models to light curves of core-collapse SNe. This allows the estimation of the most important physical parameters, like the radius of the progenitor star, the mass of the ejected envelope, the mass of the radioactive nickel synthesized during the explosion, among others. PSN J17292918+7542390 has $R_0 = 91_{-70}^{+119} \times 10^{11} \, \text{cm}$, $M_\text{ej} = 9.89_{-1.00}^{+2.10} \, \mathrm{ M}_{\odot }$, $E_{\mbox{kin}} = 0.65_{-0.18}^{+0.19} \, \text{foe}$, and $v_{\mbox{exp}} = 3332_{-347}^{+216}$ km s−1, for its progenitor radius, ejecta mass, kinetic energy, and expansion velocity, respectively. The initial nickel mass of the PSN J17292918+7542390 turned out to be $1.55_{-0.70}^{+0.75} \times 10^{-3} \, \mathrm{M}_{\odot }$. The measured photospheric velocity at the earliest observed phase is 7000 km s−1. As far as we can tell based on the small population of observed LL SNe IIP, the determined values are typical for these events.

Jäger, Jr, Zoltán↗

The host galaxies of 106 rapidly evolving transients discovered by the Dark Energy Survey

Rapidly evolving transients (RETs), also termed fast blue optical transients, are a recently discovered group of astrophysical events that display rapid luminosity evolution. RETs typically rise to peak in less than 10 d and fade within 30, a time-scale unlikely to be compatible with the decay of Nickel-56 that drives conventional supernovae (SNe). Their peak luminosity spans a range of -15 < Mg < -22.5, with some events observed at redshifts greater than 1. Their evolution on fast time-scales has hindered high-quality follow-up observations, and thus their origin and explosion/emission mechanism remains unexplained. In this paper, we present the largest sample of RETs to date, comprising 106 objects discovered by the Dark Energy Survey, and perform the most comprehensive analysis of RET host galaxies. Using deep-stacked photometry and emission lines from OzDES spectroscopy, we derive stellar masses and star formation rates (SFRs) for 49 host galaxies, and metallicities ([O/H]) for 37. Here, we find that RETs explode exclusively in star-forming galaxies and are thus likely associated with massive stars. Comparing RET hosts to samples of host galaxies of other explosive transients as well as field galaxies, we find that RETs prefer galaxies with high specific SFRs ($\langle$log(sSFR)$\rangle$ ~ -9.6), indicating a link to young stellar populations, similar to stripped-envelope SNe. RET hosts appear to show a lack of chemical enrichment, their metallicities akin to long-duration gamma-ray bursts and superluminous SN host galaxies ($\langle$12+ log (O/H)$\rangle$ ~ 9.4). There are no clear relationships between mass or SFR of the host galaxies and the peak magnitudes or decline rates of the transients themselves.

79 ASTRONOMY AND ASTROPHYSICS↗