SEARCH · Engineering Papers

Results for “outflows”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2020

NuSTAR Observations of G11.2–0.3

In this work, we present the hard X-ray view of the pulsar wind nebula in G11.2-0.3 and its central pulsar powered pulsar J1811-1925 as seen by NuSTAR. We complement the data with Chandra for a more complete picture and confirm the existence of a hard, power-law component in the shell with photon index ${\rm{\Gamma }}=2.1\pm 0.1$, which we attribute to synchrotron emission. Our imaging observations of the shell show a slightly smaller radius at higher energies, consistent with Chandra results, and we find shrinkage as a function of increased energy along the jet direction, indicating that the electron outflow in the PWN may be simpler than that seen in other young PWNe. Combining NuSTAR with INTEGRAL, we find that the pulsar spectrum can be fit by a power law with ${\rm{\Gamma }}=1.32\pm 0.07$ up to 300 keV without evidence of curvature.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

THEMIS: A Parameter Estimation Framework for the Event Horizon Telescope

The Event Horizon Telescope (EHT) provides the unprecedented ability to directly resolve the structure and dynamics of black hole emission regions on scales smaller than their horizons. This has the potential to critically probe the mechanisms by which black holes accrete and launch outflows, and the structure of supermassive black hole spacetimes. However, accessing this information is a formidable analysis challenge for two reasons. First, the EHT natively produces a variety of data types that encode information about the image structure in nontrivial ways; these are subject to a variety of systematic effects associated with very long baseline interferometry and are supplemented by a wide variety of auxiliary data on the primary EHT targets from decades of other observations. Second, models of the emission regions and their interaction with the black hole are complex, highly uncertain, and computationally expensive to construct. As a result, the scientific utilization of EHT observations requires a flexible, extensible, and powerful analysis framework. We present such a framework, Themis, which defines a set of interfaces between models, data, and sampling algorithms that facilitates future development. We describe the design and currently existing components of Themis, how Themis has been validated thus far, and present additional analyses made possible by Themis that illustrate its capabilities. Importantly, we demonstrate that Themis is able to reproduce prior EHT analyses, extend these, and do so in a computationally efficient manner that can efficiently exploit modern high-performance computing facilities. Themis has already been used extensively in the scientific analysis and interpretation of the first EHT observations of M87.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Origin of Large Meteoritic SiC Stardust Grains in Metal-rich AGB Stars

Stardust grains that originated in ancient stars and supernovae are recovered from meteorites and carry the detailed composition of their astronomical sites of origin. We present evidence that the majority of large (μm-sized) meteoritic silicon carbide (SiC) grains formed in C-rich asymptotic giant branch (AGB) stars that were more metal-rich than the Sun. In the framework of the slow neutron captures (the s process) that occur in AGB stars, the lower-than-solar $^{88}$Sr/$^{86}$Sr isotopic ratios measured in the large SiC grains can only be accompanied by Ce/Y elemental ratios that are also lower than solar and predominately observed in metal-rich barium stars—the binary companions of AGB stars. Such an origin suggests that these large grains represent the material from high-metallicity AGB stars needed to explain the s-process nucleosynthesis variations observed in bulk meteorites. In the outflows of metal-rich, C-rich AGB stars, SiC grains are predicted to be small (≃0.2 μm); large (≃μm-sized) SiC grains can grow if the number of dust seeds is 2–3 orders of magnitude lower than the standard value of 10 -13 times the number of H atoms. We therefore predict that with increasing metallicity, the number of dust seeds might decrease, resulting in the production of larger SiC grains.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

The Rise and Fall of ASASSN-18pg: Following a TDE from Early to Late Times

Here, we present nearly 500 days of observations of the tidal disruption event (TDE) ASASSN-18pg, spanning from 54 days before peak light to 441 days after peak light. Our data set includes X-ray, UV, and optical photometry, optical spectroscopy, radio observations, and the first published spectropolarimetric observations of a TDE. ASASSN-18pg was discovered on 2018 July 11 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $\textit{d}$ = 78.6 Mpc; with a peak UV magnitude of $\textit{m}$ ≃ 14, it is both one of the nearest and brightest TDEs discovered to-date. The photometric data allow us to track both the rise to peak and the long-term evolution of the TDE. ASASSN-18pg peaked at a luminosity of $\textit{L}$ ≃ 2.4 × 10 44 erg s -1 , and its late-time evolution is shallower than a flux ∝$\textit{t}^{–5/3}$ power-law model, similar to what has been seen in other TDEs. ASASSN-18pg exhibited Balmer lines and spectroscopic features consistent with Bowen fluorescence prior to peak, which remained detectable for roughly 225 days after peak. Analysis of the two-component Hα profile indicates that, if they are the result of reprocessing of emission from the accretion disk, the different spectroscopic lines may be coming from regions between ~10 and ~60 lt-days from the black hole. No X-ray emission is detected from the TDE, and there is no evidence of a jet or strong outflow detected in the radio. Our spectropolarimetric observations indicate that the projected emission region is likely not significantly aspherical, with the projected emission region having an axis ratio of ≳0.65.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

IceCube Search for High-energy Neutrino Emission from TeV Pulsar Wind Nebulae

Pulsar wind nebulae (PWNe) are the main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit nonthermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma-rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 yr of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. Additionally, in the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Kinetic Simulations of Radiative Magnetic Reconnection in the Coronae of Accreting Black Holes

We perform 2D and 3D particle-in-cell simulations of reconnection in magnetically dominated e ± plasmas subject to strong Compton cooling. Magnetic reconnection under such conditions can operate in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically dominated plasmoids with a small thermal pressure. Compton drag clears cavities inside plasmoids and also affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in nonradiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two distinct components: (1) A mildly relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the dissipated reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives ~20% of the dissipated reconnection power. It is populated by particles accelerated impulsively at X-points or "picked up" by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support magnetic reconnection as a mechanism powering hard X-ray emission from magnetically dominated regions of accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. I. A Systematic Study

Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal-dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell (PIC) simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical evolution of magnetic field during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands (ultraviolet to X-ray) tend to flare earlier than lower-frequency bands (infrared to optical) in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the variation in temporal polarization depends strongly on the guide field strength. Specifically, reconnection with a significant guide field component leads to a very high polarization degree that contradicts typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly antiparallel magnetic field lines. These rotations are at least 90° and can extend to >180°, and they may be in either direction. Furthermore, these results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly antiparallel magnetic field morphology.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above ${Y}_{{\rm{e}}}\sim 0.3$ and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. Finally, we compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The Physics of the Coronal-line Region for Galaxies in Mapping Galaxies at Apache Point Observatory

The fundamental nature and extent of the coronal-line region (CLR), which may serve as a vital tracer for active galactic nucleus (AGN) activity, remain unresolved. Previous studies suggest that the CLR is produced by AGNdriven outflows and occupies a distinct region between the broad-line region and the narrow-line region, which places it tens to hundreds of parsecs from the galactic center. Here, we investigate 10 coronal line (CL; ionization potential ≳100 eV) emitting galaxies from the Sloan Digital Sky Survey IV Mapping Galaxies at Apache Point Observatory (MaNGA) catalog with emission from one or more CLs detected at ≥5$σ$ above the continuum in at least 10 spaxels—the largest such MaNGA catalog. We find that the CLR is far more extended, reaching out to 1.3–23 kpc from the galactic center. We crossmatch our sample of 10 CL galaxies with the largest existing MaNGA AGN catalog and identify seven in it; two of the remaining three are galaxy mergers and the final one is an AGN candidate. Further, we measure the average CLR electron temperatures as ranging between 12,331 and 22,530 K, slightly above the typical threshold for pure AGN photoionization (~20,000 K) and indicative of shocks (e.g., merger induced or from supernova remnants) in the CLR. In conclusion, we reason that ionizing photons emitted by the central continuum source (i.e., AGN photoionization) primarily generate the CLs, and that energetic shocks are an additional ionization mechanism that likely produce the most extended CLRs we measure.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The Challenges Ahead for Multimessenger Analyses of Gravitational Waves and Kilonova: A Case Study on GW190425

In recent years, there have been significant advances in multimessenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, in this work, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so-called kilonova, and particularly on black hole–neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole or a neutron star. We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. II. Connection with γ-Rays

Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LAT γ -ray bands as well as closely correlated optical and γ -ray flares. The swings in optical polarization angle are also accompanied by γ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast γ -ray flares or orphan γ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the kglobal kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma β, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Radio and X-Ray Observations of the Luminous Fast Blue Optical Transient AT 2020xnd

We present deep X-ray and radio observations of the fast blue optical transient (FBOT) AT 2020xnd/ZTF 20acigmel at z = 0.2433 from 13 days to 269 days after explosion. AT 2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT 2018cow. AT 2020xnd shows luminous radio emission reaching L ν ≈ 8 × 10 29 erg s -1 Hz -1 at 20 GHz and 75 days post-explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at L X ≈ 6 × 10 42 erg s -1 . Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment, we find that AT 2020xnd launched a high-velocity outflow (v ~ 0.1c–0.2c) propagating into a dense circumstellar medium (effective $\dot{M}\approx {10}^{-3}\,{M}_{\odot }$ yr -1 for an assumed wind velocity of v w = 1000 km s -1 ). Similar to AT 2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT 2020xnd a high-redshift analog to AT 2018cow, and establish AT 2020xnd as the fourth member of the class of optically luminous FBOTs with luminous multiwavelength counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Hydrogen Burning of 29 Si and Its Impact on Presolar Stardust Grains from Classical Novae

Abstract Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small N ( 12 C)/ N ( 13 C) and N ( 14 N)/ N ( 15 N) number abundance ratios compared to solar values, and abundance excesses in 30 Si relative to 29 Si, as previously predicted by models of classical novae. We report on a direct measurement of the 29 Si(p, γ ) 30 P reaction, which strongly impacts simulated δ 29 Si values from classical novae. Our new experimental 29 Si(p, γ ) 30 P thermonuclear reaction rate differs from previous results by up to 50% in the classical nova temperature range ( T = 100–400 MK), while the rate uncertainty is reduced by up to a factor of 3. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate δ 29 Si values with much reduced uncertainties. Our results establish δ 29 Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that δ 30 Si values from nova simulations are currently not a useful diagnostic tool unless the large uncertainty of the 30 P(p, γ ) 31 S reaction rate can be significantly reduced.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The Nascent Milliquasar VT J154843.06+220812.6: Tidal Disruption Event or Extreme Accretion State Change?

We present a detailed multiwavelength follow-up of the nuclear radio flare VT J154843.06+220812.6, hereafter VT J1548. VT J1548 was selected as a ~1 mJy radio flare in 3 GHz observations from the Very Large Array Sky Survey. It is located in the nucleus of a low-mass ($\mathrm{log}{M}_{\mathrm{BH}}/{M}_{\odot }\sim 6$) host galaxy with weak or no past active galactic nuclei (AGN) activity. VT J1548 is associated with a slow rising (multiple year), bright mid-IR flare in the Wide-field Infrared Survey Explorer survey, peaking at ~10%L edd. . No associated optical transient is detected, although we cannot rule out a short, early optical flare given the limited data available. Constant late-time (~3 yr post-flare) X-ray emission is detected at ~10 42 erg s –1 . The radio spectral energy distribution is consistent with synchrotron emission from an outflow incident on an asymmetric medium. A follow-up, optical spectrum shows transient, bright, high-ionization coronal line emission ([Fe xi] λ6375, [Fe xi] λ7894, [S xii] λ7612). Transient broad Hα is also detected but without corresponding broad Hβ emission, suggesting high nuclear extinction. We interpret this event as either a tidal disruption event or an extreme flare of an AGN, in both cases obscured by a dusty torus. Although these individual properties have been observed in previous transients, the combination is unprecedented. This event highlights the importance of searches across all wave bands for assembling a sample of nuclear flares that spans the range of observable properties and possible triggers.

79 ASTRONOMY AND ASTROPHYSICS↗