Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “intergalactic medium”

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 37 records · Page 2

Optimal 1D Ly α forest power spectrum estimation – II. KODIAQ, SQUAD, and XQ-100

ABSTRACT We measure the 1D Ly α power spectrum P1D from Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ), The Spectral Quasar Absorption Database (SQUAD), and XQ-100 quasars using the optimal quadratic estimator. We combine KODIAQ and SQUAD at the spectrum level, but perform a separate XQ-100 estimation to control its large resolution corrections in check. Our final analysis measures P1D at scales k < 0.1 s km−1 between redshifts $z$ = 2.0–4.6 using 538 quasars. This sample provides the largest number of high-resolution, high-S/N observations; and combined with the power of optimal estimator it provides exceptional precision at small scales. These small-scale modes (k ≳ 0.02 s km−1), unavailable in Sloan Digital Sky Survey and Dark Energy Spectroscopic Instrument analyses, are sensitive to the thermal state and reionization history of the intergalactic medium, as well as the nature of dark matter. As an example, a simple Fisher forecast analysis estimates that our results can improve small-scale cut-off sensitivity by more than a factor of 2.

79 ASTRONOMY AND ASTROPHYSICS↗

Lyman-α constraints on cosmic heating from dark matter annihilation and decay

We derive new constraints on models of decaying and annihilating dark matter (DM) by requiring that the energy injected into the intergalactic medium (IGM) not overheat it at late times, when measurements of the Lyman-α forest constrain the IGM temperature. We improve upon previous analyses by using the recently developed DarkHistory code package, which self-consistently takes into account additional photoionization and photoheating processes due to reionization and DM sources. Our constraints are robust to the uncertainties of reionization and competitive with leading limits on sub-GeV DM that decays preferentially to electrons

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrogen reionization ends by z = 5.3: Lyman-α optical depth measured by the XQR-30 sample

ABSTRACT The presence of excess scatter in the Ly-α forest at z ∼ 5.5, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionization. However, low data quality and systematic uncertainties complicate the use of Ly-α transmission as a precision probe of reionization’s end stages. In this paper, we assemble a sample of 67 quasar sightlines at z > 5.5 with high signal-to-noise ratios of >10 per ≤15 km s−1 spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR > 20 per pixel) spectra of 30 quasars at z > 5.7. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-α systems. We present improved measurements of the mean Ly-α transmission over 4.9 < z < 6.1. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-α transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to z ≤ 5.2 (<1σ), and mild tension (∼2.5σ) at z = 5.3. Homogeneous UVB models are ruled out by excess Ly-α transmission scatter at z ≥ 5.4 with high confidence (>3.5σ). Our results indicate that reionization-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least z = 5.3 (t = 1.1 Gyr after the big bang). This is further evidence for a late end to reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

Recovering Density Fields inside Quasar Proximity Zones at z ~ 6

The matter density field at z ~ 6 is very challenging to probe. One of the traditional probes of the low density intergalactic medium that works successfully at lower redshifts is the Lyα forest in quasar spectra. However, at the end of reionization, the residual neutral hydrogen usually creates saturated absorption, thus much of the information about the gas density is lost. Luckily, in a quasar proximity zone, the ionizing radiation is exceptionally intense, thus creating a large region with nonzero transmitted flux. In this study we use the synthetic spectra from simulations to investigate how to recover the density fluctuations inside the quasar proximity zones. We show that, under ideal conditions, the density can be recovered accurately with a small scatter. We also discuss how systematics such as the quasar continuum fitting and reionization models affect the results. This study shows that by analyzing the absorption features inside quasar proximity zones we can potentially constrain quasar properties and the environments they reside in.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of varying inhomogeneous reionization histories on metrics of Ly α opacity

The epoch of hydrogen reionization is complete by z = 5⁠, but its progression at higher redshifts is uncertain. Measurements of Ly α forest opacity show large scatter at z < 6⁠, suggestive of spatial fluctuations in neutral fraction, temperature, or ionizing background, either individually or in combination. However, there are degeneracies in the impact of such fluctuations, necessitating careful modelling. We develop a framework for modelling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. We generate several reionization histories using seminumerical code AMBER, and implement them in the Nyx cosmological hydrodynamics code to examine the impact on the evolution of gas within the simulation and the associated metrics of the Ly α forest opacity. We find that the pressure smoothing scale within the intergalactic medium is strongly correlated with the adiabatic index of the temperature–density relation. We find that while models with 20 000 K photoheating at reionization are better able to reproduce the shape of the observed z = 5 1D flux power spectrum than colder ones, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at z < 6.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal 1D Ly α forest power spectrum estimation – I. DESI-lite spectra

ABSTRACT The 1D Ly α forest flux power spectrum P1D is sensitive to scales smaller than a typical galaxy survey, and hence ties to the intergalactic medium’s thermal state, suppression from neutrino masses, and new dark matter models. It has emerged as a competitive framework to study new physics, but also has come with various challenges and systematic errors in analysis. In this work, we revisit the optimal quadratic estimator for P1D, which is robust against the relevant problems such as pixel masking, time evolution within spectrum, and quasar continuum errors. We further improve the estimator by introducing a fiducial power spectrum, which enables us to extract more information by alleviating the discreteness of band powers. We meticulously apply our method to synthetic Dark Energy Spectroscopic Instrument (DESI) spectra and demonstrate how the estimator overcomes each challenge. We further apply an optimization scheme that approximates the Fisher matrix to three elements per row and reduces computation time by 60 per cent. We show that we can achieve per cent precision in P1D with 5-yr DESI data in the absence of systematics and provide forecasts for different spectral qualities.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on dark photon dark matter from Lyman- α forest simulations and an ultrahigh signal-to-noise quasar spectrum

The ultralight dark photon is a well-motivated, hypothetical dark matter candidate. In a dilute plasma, they can resonantly convert into photons, and heat up the intergalactic medium between galaxies. In this work, we explore the dark photon dark matter parameter space by comparing synthetic Lyman- α forest data from cosmological hydrodynamical simulations to observational data from VLT/UVES of the quasar HE0940-1050 ( z em = 3.09 ). We use a novel flux normalization technique that targets underdense gas, reshaping the flux probability distribution. Not only do we place robust constraints on the kinetic mixing parameter of dark photon dark matter, but notably our findings suggest that this model can still reconcile simulated and observed Doppler parameter distributions of z ∼ 0 Lyman- α lines, as seen by HST/COS. This work opens new pathways for the use of the Lyman- α forest to explore new physics, and can be extended to other scenarios such as primordial black hole evaporation, dark matter decay, and annihilation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reconstructing Ly$α$ Fields from Low-resolution Hydrodynamical Simulations with Deep Learning

Hydrodynamical cosmological simulations are a powerful tool for accurately predicting the properties of the intergalactic medium (IGM) and for producing mock skies that can be compared against observational data. However, the need to resolve density fluctuation in the IGM puts a stringent requirement on the resolution of such simulations, which in turn limits the volumes that can be modeled, even on the most powerful supercomputers. In this work, we present a novel modeling method that combines physics-driven simulations with data-driven generative neural networks to produce outputs that are qualitatively and statistically close to the outputs of hydrodynamical simulations employing eight times higher resolution. We show that the Ly$α$ flux field, as well as the underlying hydrodynamic fields, have greatly improved statistical fidelity over a low-resolution simulation. Importantly, the design of our neural network allows for sampling multiple realizations from a given input, enabling us to quantify the model uncertainty. Using test data, we demonstrate that this model uncertainty correlates well with the true error of the Ly$α$ flux prediction. Ultimately, our approach allows for training on small simulation volumes and applying it to much larger ones, opening the door to producing accurate Ly$α$ mock skies in volumes of Hubble size, as will be probed with DESI and future spectroscopic sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

TARDIS. I. A Constrained Reconstruction Approach to Modeling the z ~ 2.5 Cosmic Web Probed by Lyα Forest Tomography

Recent Lyα forest tomography measurements of the intergalactic medium (IGM) have revealed a wealth of cosmic structures at high redshift (z ~ 2.5). In this work, we present the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS), a new chronocosmographic analysis tool for understanding the formation and evolution of these observed structures. Here, we use maximum likelihood techniques with a fast nonlinear gravitational model to reconstruct the initial density field of the observed regions. We find that TARDIS allows accurate reconstruction of smaller-scale structures than standard Wiener-filtering techniques. Applying this technique to mock Lyα forest data sets that simulate ongoing and future surveys such as CLAMATO, Subaru PFS, or the ELTs, we are able to infer the underlying matter density field at observed redshift and classify the cosmic web structures. We find good agreement with the underlying truth in both the characteristic eigenvalues and eigenvectors of the pseudo-deformation tensor, with the eigenvalues inferred from 30 m class telescopes correlated at r = 0.95 relative to the truth. As an output of this method, we are able to further evolve the inferred structures to late time (z = 0) and also track the trajectories of coeval z = 2.5 galaxies to their z = 0 cosmic web environments.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconstructing large-scale temperature profiles around z ∼ 6 quasars

ABSTRACT High-redshift quasars ionize He ii into He iii around them, heating the intergalactic medium in the process and creating large regions with elevated temperature. In this work, we demonstrate a method based on a convolutional neural network (CNN) to recover the spatial profile for T0, the temperature at the mean cosmic density, in quasar proximity zones. We train the neural network with synthetic spectra drawn from a Cosmic Reionization on Computers simulation. We discover that the simple CNN is able to recover the temperature profile with an accuracy of ≈1400 K in an idealized case of negligible observational uncertainties. We test the robustness of the CNN and discover that it is robust against the uncertainties in quasar host halo mass, quasar continuum, and ionizing flux. We also find that the CNN has good generality with regard to the hardness of quasar spectra. This shows that with noiseless spectra, one could use a simple CNN to distinguish gas inside or outside the He iii region created by the quasar. Because the size of the He iii region is closely related to the total quasar lifetime, this method has great potential in constraining the quasar lifetime on ∼Myr time-scales. However, noise poses a big problem for accuracy and could downgrade the accuracy to ≈2340 K even for very high signal-to-noise (≳50) spectra. Future studies are needed to reduce the error associated with noise to constrain the lifetimes of reionization epoch quasars with currently available data.

79 ASTRONOMY AND ASTROPHYSICS↗

The Cosmic Evolution of C IV Absorbers at 1.4 < z < 4.5: Insights from 100,000 Systems in DESI Quasars

We present the largest catalog to date of triply ionized carbon (C IV ) absorbers detected in quasar spectra from the Dark Energy Spectroscopic Instrument. Using an automated matched-kernel convolution method with adaptive signal-to-noise thresholds, we identify 101,487 C IV systems in the redshift range 1.4 < z < 4.5 from 300,637 quasar spectra. Completeness is estimated via Monte Carlo simulations, and the catalog is 50% complete at EW C IV ≥ 0.4 Å. The differential equivalent width frequency distribution declines exponentially and shows weak redshift evolution. The absorber incidence per unit comoving path increases by a factor of 2–5 from z ≈ 4.5 to z ≈ 1.4, with stronger redshift evolution for strong systems. Using column densities derived from the apparent optical depth method, we constrain the cosmic mass density of C IV , Ω C IV , which increases by a factor of ∼3.8 from (0.82 ± 0.05) × 10 −8 at z ≈ 4.5 to (3.16 ± 0.2) × 10 −8 at z ≈ 1.4. From Ω C IV , we estimate a lower limit on intergalactic medium metallicity ${\mathrm{log}}({Z}_{{\rm{IGM}}}/{Z}_{\odot })\gtrsim -3.25$ at z ∼ 2.3, with a smooth decline at higher redshifts. These trends trace the cosmic star formation history and He II photoheating rate, suggesting a link between C IV enrichment, star formation, and UV background over ∼3 Gyr. The catalog also provides a critical resource for future studies connecting circumgalactic metals to galaxy evolution, especially near cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining IGM enrichment and metallicity with the C iv forest correlation function

ABSTRACT The distribution and abundance of metals in the diffuse intergalactic medium (IGM) have implications for galaxy formation and evolution models, and has been argued to be sensitive to the Universe’s reionization history. However, reduced sensitivity in the near-IR implies that probing IGM metals at z > 4 is currently out of reach with the traditional method of detecting individual absorbers. We present a new technique based on clustering analysis that enables the detection of these weak IGM absorbers. We investigate the two-point correlation function (2PCF) of the ${\rm C\, {\small IV}}$ forest as a probe of IGM metallicity and enrichment topology by simulating the z = 4.5 IGM with models of inhomogeneous metal distributions. The 2PCF of the ${\rm C\, {\small IV}}$ forest demonstrates a clear peak at a characteristic separation corresponding to the doublet separation of the ${\rm C\, {\small IV}}$ line.The peak amplitude scales quadratically with metallicity, while enrichment topology affects both the shape and amplitude of the 2PCF. For models consistent with the distribution of metals at z ∼ 3, we find that we can constrain [C/H] to within 0.2 dex, log$\, M_{\rm {min}}$ to within 0.4 dex, and R to within 15 per cent. We show that CGM absorbers can be reliably identified and masked, thus recovering the underlying IGM signal. The auto-correlation of the metal-line forest presents a compelling avenue to constrain the IGM metallicity and enrichment topology with high precision at z > 4, thereby pushing such measurements into the Epoch of Reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Convergence of small scale Ly α structure at high- z under different reionization scenarios

ABSTRACT The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc and 5 − 80 h−1 kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ~1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ( |$\langle f_{\rm esc}^{\rm sf} \rangle$| ) and AGN ( |$f_{\rm esc}^{\rm bh}$| ) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ( |$M_*\lt 10^9 \, {\rm \rm M_\odot }$| ) galaxies at z > 6 with the AGN contribution (driven by |$M_{bh}\gt 10^6 \, {\rm \rm M_\odot }$| black holes in |$M_* \gtrsim 10^9\, {\rm \rm M_\odot}$| galaxies) dominating at lower redshifts; (ii) AGN only contribute |$10-25{{\ \rm per\ cent}}$| to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of |$\langle f_{\rm esc}^{\rm sf} \rangle$| is shallower than |$f_{\rm esc}^{\rm bh}$| , at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ~ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.

79 ASTRONOMY AND ASTROPHYSICS↗

A framework to measure the properties of intergalactic metal systems with two-point flux statistics

ABSTRACT The abundance, temperature, and clustering of metals in the intergalactic medium are important parameters for understanding their cosmic evolution and quantifying their impact on cosmological analysis with the Ly α forest. The properties of these systems are typically measured from individual quasar spectra redward of the quasar’s Ly α emission line, yet that approach may provide biased results due to selection effects. We present an alternative approach to measure these properties in an unbiased manner with the two-point statistics commonly employed to quantify large-scale structure. Our model treats the observed flux of a large sample of quasar spectra as a continuous field and describes the one-dimensional, two-point statistics of this field with three parameters per ion: the abundance (column density distribution), temperature (Doppler parameter), and clustering (cloud–cloud correlation function). We demonstrate this approach on multiple ions (e.g. ${\rm C\, \small {\rm IV}}$ , ${\rm Si\, \small {\rm IV}}$ , and ${\rm Mg\, \small {\rm II}}$ ) with early data from the Dark Energy Spectroscopic Instrument (DESI) and high-resolution spectra from the literature. Our initial results show some evidence that the ${\rm C\, \small {\rm IV}}$ abundance is higher than previous measurements and evidence for abundance evolution over time. The first full year of DESI observations will have over an order of magnitude more quasar spectra than this study. In a future paper, we will use those data to measure the growth of clustering and its impact on the Ly α forest, as well as test other DESI analysis infrastructure such as the pipeline noise estimates and the resolution matrix.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal 1D Ly α forest power spectrum estimation – III. DESI early data

ABSTRACT The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and discuss the key sources of systematic error that we need to address with the upcoming DESI Year 1 analysis.

79 ASTRONOMY AND ASTROPHYSICS↗