Search NASA⌕ Search

SEARCH · Search NASA

Results for “compact binaries”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dust supply to close binary systems

Binary systems can be born surrounded by circumbinary discs. The gaseous disc around either of the two stellar companions can have its life extended by the supply of mass arriving from the circumbinary disc. The objective of this study is to investigate the gravitational interactions exerted by a compact and eccentric binary system on the circumbinary and circumprimary discs, and the resulting transport of gas and solids between the disc components. We assume that the gas in the system behaves as a fluid, and we model its evolution by means of high-resolution hydrodynamical simulations. Dust grains are modelled as Lagrangian particles that interact with the gas and the stars. Our models indicate that significant fluxes of gas and dust proceed from the circumbinary disc towards the circumprimary disc. For the applied system parameters, grains of certain sizes are segregated outside the tidal gap generated by the stars. Consequently, the size distribution of the transported dust is not continuous, but presents a gap in the millimetre size range. In close binaries, the lifetime of an isolated circumprimary disc is found to be short, ~10 5 years, because of its low mass. However, because of the influx of gas from beyond the tidal gap, the disc around the primary star can survive much longer, ~10 6 years, as long as gas accretion from the circumbinary disc continues. The supply of solids and the extended lifetime of a circumbinary disc also aids in the possible formation of giant planets. Compared to close binary systems without a circumbinary disc, we expect a higher frequency of singleplanet or multiple-planet systems. Additionally, a planetesimal or debris belt can form in the proximity of the truncation radius of the circumprimary disc and/or around the location of the exterior edge of the tidal gap.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Exotic compact objects: The dark white dwarf

Several dark matter models allow for the intriguing possibility of exotic compact object formation. These objects might have unique characteristics that set them apart from their baryonic counterparts. Furthermore, gravitational wave observations of their mergers may provide the only direct window on a potentially entirely hidden sector. Here, we discuss dark white dwarfs, starting with an overview of the microphysical model and analytic scaling relations of macroscopic properties derived from the nonrelativistic limit. We use the full relativistic formalism to confirm these scaling relations and demonstrate that dark white dwarfs, if they exist, would have masses and tidal deformabilities that are very different from those of baryonic compact objects. Further, and most importantly, we demonstrate that dark white dwarf mergers would be detectable by current or planned gravitational observatories across several orders of magnitude in the particle-mass parameter space. Lastly, we find universal relations analogous to the compactness-Love and binary Love relations in neutron star literature. Using these results, we show that gravitational wave observations would constrain the properties of the dark matter particles constituting these objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Electrical conductivity of porous binary powder mixtures

Simultaneous data of the quasi-static compaction and electrical conductivity of porous, binary powder mixtures have been collected as a function of bulk density. The powder mixtures consist of a metal conductor, either titanium or iron, an insulator, and pores filled with ambient air. The data show a dependency of the conductivity in terms of relative bulk density and metal volume fraction on conductor type and conductor particle characteristics of size and shape. Finally, finite element models using particle domains generated by discrete element method are used to simulate the bulk conductivity near its threshold while the general effective media equation is used to model the conductivity across the compression regime.

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

XRISM spectroscopy on orbital modulation of Fe Ly α lines in Cygnus X-3

To understand physical processes such as mass transfer and binary evolution in X-ray binaries, the orbital parameters of the system are fundamental and crucial information. Cygnus X-3 is a high-mass X-ray binary composed of a compact object of unknown nature and a Wolf–Rayet star, which is of great interest in the context of wind-fed mass accretion and binary evolution. Here we present XRISM/Resolve high-resolution spectroscopy, focusing on the Fe Ly$\alpha$ lines when the source was in its hypersoft state. We perform an orbital phase-resolved spectral analysis of the lines to study the orbital modulation of the emission and absorption lines. It is found that the emission lines reflect the orbital motion of the compact object whose estimated velocity amplitude is $430^{+150}_{-140}$ km s$^{-1}$, while the absorption lines show a variation that can be interpreted as originating from the stellar wind. We discuss possible mass ranges for the binary components using the mass function with the estimated value of the velocity amplitude in this work, combined with the relation between the mass-loss rate and the orbital period derivative and the empirical mass and mass-loss rate relation for Galactic Wolf–Rayet stars. They are constrained to be (1.3–5.1) and (9.3–12) M$_\odot$, respectively, for the assumed inclination angle of $i =25^\circ$, and become more relaxed, (1.3–24) and (9.3–16) M$_\odot$, for $i = 35^\circ$. Thus, it remains unclear whether the system harbors a black hole or a neutron star.

X-rays: binaries↗
DOE OSTI2025

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗
DOE OSTI2022

The role of supernova convection for the lower mass gap and the isolated binary formation of gravitational wave sources

Understanding astrophysical phenomena involving compact objects requires an insight about the engine behind core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of detected black hole-black hole, black holeneutron star and neutron star-neutron star mergers. The timescale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants (Fryer et al. 2022) and check how they impact population of double compact object (DCO) mergers formed via isolated binary evolution. The new formulas give one ability to test wide spectrum of assumptions on the convection growth time. In particular, different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low mass black holes. We present distribution of masses, mass ratios and the local merger rate densities of DCO for different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that mass distribution of DCO mergers is sensitive to adopted assumption on SN convection growth timescale up to $m_1$ + $m_2$ ≲ 35$M_⊙$. Between two extreme tested variants the probability of compact object formation within the mass gap may differ up to ~ 2 orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Catalog of Planetary Nebulae Detected by GALEX and Corollary Optical Surveys

Planetary nebulae (PNs) consist of an ionized envelope surrounding a hot central star (CSPN) that emits mostly at ultraviolet (UV) wavelengths. UV observations, therefore, provide important information on both the CSPN and the nebula. We have matched the PNs in the Hong Kong/Australian Astronomical Observatory/Strasbourg Hα catalog with the Galaxy Evolution Explorer (GALEX) UV sky surveys, the Sloan Digital Sky Survey (SDSS) data release 16, and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) PS1 second release. A total of 671 PNs were observed by GALEX with the far-UV (FUV; 1344 Å–1786 Å) and/or the near-UV (NUV; 1771 Å–2831 Å) detectors on (GUVPNcat); 83 were observed by SDSS (PNcatxSDSSDR16) and 1819 by Pan-STARRS (PNcatxPS1MDS). We merged a distilled version of these matched catalogs into GUVPNcatxSDSSDR16xPS1MDS, which contains a total of 375 PNs with both UV and optical photometry over a total spectral coverage of ~1540 Å–9610 Å. We analyzed separately 170 PNs resolved in GALEX images and determined their UV radius by applying a flux profile analysis. The CSPN flux could be extracted separately from the PN emission for 8 and 50 objects with SDSS and Pan-STARRS counterparts, respectively. The multiband photometry was used to distinguish between compact and extended PNs and CSPNs (binary CSPNs) by color–color diagram analysis. We found that compact PN candidates could be identified by using the r - i < - 0.4 and -1 < FUV-NUV < 1 colors, whereas binary CSPN candidates in given T eff ranges (all with color r-i > -0.4) can be identified in the color region (FUV-NUV) ≤ 6(r-i)+1.3, -0.8 < FUV-NUV < 0.4, and r-i < 0.75.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Radius and equation of state constraints from massive neutron stars and GW190814

Motivated by the unknown nature of the 2.50–2.67M⊙ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. Here we demonstrate that a modification of the highly uncertain suprasaturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of 2.6M⊙ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of ≈ 2.0M⊙ neutron stars (but not necessarily the radii of 1.4 M⊙ neutron stars), which further motivates future NICER radius measurements of PSR J1614 – 2230 and PSR J0740+6620.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Assessing equation of state-independent relations for neutron stars with nonparametric models

Relations between neutron star properties that do not depend on the nuclear equation of state offer insights on neutron star physics and have practical applications in data analysis. Such relations are obtained by fitting to a range of phenomenological or nuclear physics equation of state models, each of which may have varying degrees of accuracy. In this study we revisit commonly used relations and reassess them with a very flexible set of phenomenological nonparametric equation of state models that are based on Gaussian processes. Our models correspond to two sets: equations of state which mimic hadronic models, and equations of state with rapidly changing behavior that resemble phase transitions. Here we quantify the accuracy of relations under both sets and discuss their applicability with respect to expected upcoming statistical uncertainties of astrophysical observations. We further propose a goodness-of-fit metric which provides an estimate for the systematic error introduced by using the relation to model a certain equation-of-state set. Overall, the nonparametric distribution is more poorly fit with existing relations, with the I–Love–Q relations retaining the highest degree of universality. Fits degrade for relations involving the tidal deformability, such as the binary-Love and compactness-Love relations, and when introducing phase transition phenomenology. For most relations, systematic errors are comparable to current statistical uncertainties under the nonparametric equation of state distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The role of supernova convection for the lower mass gap in the isolated binary formation of gravitational wave sources

ABSTRACT Understanding the astrophysical phenomena involving compact objects requires an insight about the engine behind the core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of the detected black hole–black hole, black hole–neutron star, and neutron star–neutron star mergers. The time-scale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants and check how they impact the population of double compact object (DCO) mergers formed via the isolated binary evolution. The new formulas give one the ability to test a wide spectrum of assumptions on the convection growth time. In particular, the different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low-mass black holes. We present the distribution of masses, mass ratios, and the local merger rate densities of DCO for the different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that the mass distribution of DCO mergers is sensitive to the adopted assumption on the SN convection growth time-scale up to $m_1+m_2 \lesssim 35 \, \mathrm{M}_{\odot }$. Between the two extreme tested variants the probability of compact object formation within the mass gap may differ by up to approximately two orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

ARMADA. I. Triple Companions Detected in B-type Binaries α Del and ν Gem

Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ∼10 μas precision. This precision makes it possible to detect “wobbles” in the binary motion due to the gravitational pull from additional short-period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the Michigan Infra-Red Combiner (MIRC)/MIRC-X instrument at the Center for High Angular Resoloution Astronomy (CHARA) array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of microarcseconds level for <02 binaries. We test our instrument performance on a known triple system, κ Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30 day companion to α Del, and a 50 day companion to ν Gem. We also collected radial velocity data for α Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published radial velocity orbit for the inner pair of ν Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Peeking Between the Pulses: The Far-UV Spectrum of the Previously Unseen White Dwarf in AR Scorpii

The compact object in the interacting binary AR Sco has widely been presumed to be a rapidly rotating, magnetized white dwarf (WD), but it has never been detected directly. Isolating its spectrum has proven difficult because the spin-down of the WD generates pulsed synchrotron radiation that far outshines the WD's photosphere. As a result, a previous study of AR Sco was unable to detect the WD in the averaged far-ultraviolet spectrum from a Hubble Space Telescope (HST) observation. In an effort to unveil the WD's spectrum, we reanalyze these HST observations by calculating the average spectrum in the troughs between synchrotron pulses. We identify weak spectral features from the previously unseen WD and estimate its surface temperature to be 11,500 ± 500 K. Additionally, during the synchrotron pulses, we detect broad Lyα absorption consistent with hot WD spectral models. We infer the presence of a pair of hotspots, with temperatures between 23,000 and 28,000 K, near the magnetic poles of the WD. As the WD is not expected to be accreting from its companion, we describe two possible mechanisms for heating the magnetic poles. Here, the Lyα absorption of the hotspots appears relatively undistorted by Zeeman splitting, constraining the WD's field strength to be ≲100 MG, but the data are insufficient to search for the subtle Zeeman splits expected at lower field strengths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Origin of spin–orbit misalignments: the microblazar V4641 Sgr

ABSTRACT Of the known microquasars, V4641 Sgr boasts the most severe lower limit (>52°) on the misalignment angle between the relativistic jet axis and the binary orbital angular momentum. Assuming the jet and black hole spin axes coincide, we attempt to explain the origin of this extreme spin–orbit misalignment with a natal kick model, whereby an aligned binary system becomes misaligned by a supernova kick imparted to the newborn black hole. The model inputs are the kick velocity distribution, which we measure customized to V4641 Sgr, and the immediate pre/post-supernova binary system parameters. Using a grid of binary stellar evolution models, we determine post-supernova configurations that evolve to become consistent with V4641 Sgr today and obtain the corresponding pre-supernova configurations by using standard prescriptions for common envelope evolution. Using each of these potential progenitor system parameter sets as inputs, we find that a natal kick struggles to explain the origin of the V4641 Sgr spin–orbit misalignment. Consequently, we conclude that evolutionary pathways involving a standard common envelope phase followed by a supernova kick are highly unlikely for V4641 Sgr. An alternative interpretation is that the jet axis does not reliably trace the black hole spin axis. Our results raise concerns about compact object merger statistics gleaned from binary population synthesis models, which rely on unverified prescriptions for common envelope evolution and natal kicks. We also challenge the spin–orbit alignment assumption routinely invoked to measure black hole spin magnitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamic Excitonic Beam Switching with Atomically‐Thin Binary Blazed Gratings

Abstract Beam steering metasurfaces are ultra‐compact optical coatings that offer on‐demand redirection of optical power to specific diffraction orders. To achieve this, spatial gradients are commonly introduced in the phase of light scattered by plasmon or Mie resonant nanoparticles within the metasurface grating's unit cell. However, these phase gradients are oftentimes difficult to tune post‐fabrication. Recently, excitons in monolayer 2D semiconductors have emerged as a new metasurface building block, due to their strong and electrically‐tunable resonant light‐matter interaction. These 2D excitonic metasurfaces offer the tantalizing prospect of beam switching within a single monolayer. Here, it is demonstrated how the 2D analog of binary blazed gratings enables such beam switching by mere nanopatterning of a large monolayer WS 2 , even though nanoscale ribbons of WS 2 do not support geometrical resonances. By introducing a gradient in the nanoribbon width within the metasurface unit cell, an amplitude gradient combined with a small phase gradient in the scattered fields results in asymmetric diffraction efficiencies. Using a scattered‐field analysis, it is shown that these gradients can be further engineered via interference effects with the substrate reflection. Finally, the electrical tunability of the exciton resonance is leveraged to achieve selective and dynamic beam switching with an atomically‐thin metasurface.

Guarneri, Ludovica [Van der Waals‐Zeeman Institute↗
DOE OSTI2021

Orbital Evolution of Binary Black Holes in Active Galactic Nucleus Disks: A Disk Channel for Binary Black Hole Mergers?

We perform a series of high-resolution 2D hydrodynamical simulations of equal-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) accretion disks to study whether these binaries can be driven to merger by the surrounding gas. We find that the gravitational softening adopted for the BBH has a profound impact on this result. When the softening is less than 10% of the binary separation, we show that, in agreement with recent simulations of isolated equal-mass binaries, prograde BBHs expand in time rather than contract. Eventually, however, the binary separation becomes large enough that the tidal force of the central AGN disrupts them. Only when the softening is relatively large do we find that prograde BBHs harden. We determine through detailed analysis of the binary torque, that this dichotomy is due to a loss of spiral structure in the circum-single disks orbiting each black hole when the softening is a significant fraction of the binary separation. Properly resolving these spirals—both with high resolution and small softening—results in a significant source of binary angular momentum. Only for retrograde BBHs do we find consistent hardening, regardless of softening, as these BBHs lack the important spiral structure in their circum-single disks. This suggests that the gas-driven inspiral of retrograde binaries can produce a population of compact BBHs in the gravitational-wave-emitting regime in AGN disks, which may contribute a large fraction to the observed BBH mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Los Alamos National Laboratory: SULI Appointment

The progenitors of gamma-ray bursts - the brightest explosions in our universe - are still poorly understood. It has been suggested that massive stars collapsing in interacting binary systems (i.e. with a closely orbiting compact object companion) can provide the necessary angular momentum and dense surrounding medium to explain observations of radio-loud gamma-ray bursts. We explore this progenitor system using the population synthesis code COSMIC. We explore the dependence of the properties of an interacting binary systems - including the remnant types and their masses - on metallicity and mass raction of the binary system. We also explore how the rates of populations of these systems depend on metallicity. Our findings help constrain what properties lead to viable GRB progenitors, especially radio loud GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

On the Anticorrelation between Duration and Redshift in Gamma-Ray Bursts

Abstract For gamma-ray bursts (GRBs) with durations greater than two seconds (so-called long GRBs), the intrinsic prompt gamma-ray emission appears, on average, to last longer for bursts at lower redshifts. We explore the nature of this duration–redshift anticorrelation, describing systems and conditions in which this cosmological evolution could arise. In particular, we explore its dependence on the metallicity of a massive star progenitor, because we can securely count on the average stellar metallicity to increase with decreasing redshift. Although stars with higher metallicity/lower redshift lose mass and angular momentum through line-driven winds, in some cases these stars are able to form more extended accretion disks when they collapse, potentially leading to longer-duration GRBs. We also examine how this duration–redshift trend may show up in interacting binary models composed of a massive star and compact object companion, recently suggested to be the progenitors of radio-bright GRBs. Under certain conditions, mass loss and equation-of-state effects from massive stars with higher metallicity and lower redshift can decrease the binary separation. This can then lead to spin-up of the massive star and allow for a longer-duration GRB upon the massive star’s collapse. Finally, the duration–redshift trend may also be supported by a relatively larger population of small-separation binaries born in situ at low redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Occurrence of Gravitational Collapse in the Accreting Neutron Stars of Binary-driven Hypernovae

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon–oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB's observational properties and taxonomy. We perform three-dimensional, smoothed particle hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state, with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, t col . We include a nonzero angular momentum of the NSs and find that t col ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about 5 minute orbital period), promptly form a BH, and release ≳10 52 erg of energy. They form NS–BH binaries with tens of kiloyears merger timescales by gravitational-wave emission. BdHNe II and III do not form BHs, and release ~10 50 –10 52 erg and ≲10 50 erg of energy, respectively. They form NS–NS binaries with a range of merger timescales larger than for NS–BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a nonrotating NS. Magnetic braking by a 10 13 G field can delay BH formation, leading to BH–BH or NS–BH with tens of kiloyears merger timescales.

79 ASTRONOMY AND ASTROPHYSICS↗