Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cents”

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 505 records · Page 28

Light elements Na and Al in 58 bulge spheroid stars from APOGEE

ABSTRACT We identified a sample of 58 candidate stars with metallicity [Fe/H] ≲ −0.8 that likely belong to the old bulge spheroid stellar population, and analyse their Na and Al abundances from Apache Point Observatory Galactic Evolution Experiment (APOGEE) spectra. In a previous work, we inspected APOGEE-Stellar Parameter and Chemical Abundance Pipeline abundances of C, N, O, Mg, Al, Ca, Si, and Ce in this sample. Regarding Na lines, one of them appears very strong in about 20 per cent of the sample stars, but it is not confirmed by other Na lines, and can be explained by sky lines, which affect the reduced spectra of stars in a certain radial velocity range. The Na abundances for 15 more reliable cases were taken into account. Al lines in the H band instead appear to be very reliable. Na and Al exhibit a spread in abundances, whereas no spread in N abundances is found, and we found no correlation between them, indicating that these stars could not be identified as second-generation stars that originated in globular clusters. We carry out the study of the behaviour of Na and Al in our sample of bulge stars and literature data by comparing them with chemodynamical evolution model suitable for the Galactic bulge. The Na abundances show a large spread, and the chemodynamical models follow the main data, whereas for aluminum instead, the models reproduce very satisfactorily the nearly secondary-element behaviour of aluminum in the metallicity range below [Fe/H] ≲ −1.0. For the lower-metallicity end ([Fe/H < −2.5), hypernovae are assumed to be the main contributor to yields.

Astronomy & Astrophysics↗

Galaxy Zoo DESI: Detailed morphology measurements for 8.7M galaxies in the DESI Legacy Imaging Surveys

ABSTRACT We present detailed morphology measurements for 8.67 million galaxies in the DESI Legacy Imaging Surveys (DECaLS, MzLS, and BASS, plus DES). These are automated measurements made by deep learning models trained on Galaxy Zoo volunteer votes. Our models typically predict the fraction of volunteers selecting each answer to within 5–10 per cent for every answer to every GZ question. The models are trained on newly collected votes for DESI-LS DR8 images as well as historical votes from GZ DECaLS. We also release the newly collected votes. Extending our morphology measurements outside of the previously released DECaLS/SDSS intersection increases our sky coverage by a factor of 4 (5000–19 000 deg2) and allows for full overlap with complementary surveys including ALFALFA and MaNGA.

Astronomy & Astrophysics↗

Rates and properties of Type Ia supernovae in galaxy clusters within the dark energy survey

We identify 66 photometrically classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) that have occurred within red-sequence selected galaxy clusters. We compare light-curve and host galaxy properties of the cluster SNe to 1024 DES SNe Ia located in field galaxies, the largest comparison of two such samples at high redshift (z > 0.1). We find that cluster SN light curves decline faster than those in the field (97.7 per cent confidence). However, when limiting these samples to host galaxies of similar colour and mass, there is no significant difference in the SN light-curve properties. Motivated by previous detections of a higher-normalized SN Ia delay-time distribution in galaxy clusters, we measure the intrinsic rate of SNe Ia in cluster and field environments. We find the average ratio of the SN Ia rate per galaxy between high-mass (⁠|$10\le \log \mathrm{(\mathit{ M}_{*}/{\rm \mathit{ M}}_{\odot })} \le 11.25$|⁠) cluster and field galaxies to be 0.594 ± 0.068. This difference is mass-dependent, with the ratio declining with increasing mass, which suggests that the stellar populations in cluster hosts are older than those in field hosts. We show that the mass-normalized rate (or SNe per unit mass) in massive–passive galaxies is consistent between cluster and field environments. Additionally, both of these rates are consistent with rates previously measured in clusters at similar redshifts. We conclude that in massive–passive galaxies, which are the dominant hosts of cluster SNe, the cluster delay-time distribution is comparable to the field.

79 ASTRONOMY AND ASTROPHYSICS↗

The Calar Alto Legacy Integral Field Area survey: extended and remastered data release

ABSTRACT This paper describes the extended data release (eDR) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. It comprises science-grade quality data for 895 galaxies obtained with the Potsdam Multi Aperture Spectograph/PPak instrument at the 3.5-m telescope at the Calar Alto Observatory along the last 12 yr, using the V500 setup [3700–7500 Å, 6 Å/full-width at half-maximum (FWHM)] and the CALIFA observing strategy. It includes galaxies of any morphological type, star formation stage, a wide range of stellar masses (∼107–1012 M⊙), at an average redshift of ∼0.015 (90 per cent within 0.005 < z < 0.05). Primarily selected based on the projected size and apparent magnitude, we demonstrate that it can be volume corrected resulting in a statistically limited but representative sample of the population of galaxies in the nearby Universe. All the data were homogeneous re-reduced, introducing a set of modifications to the previous reduction. The most relevant is the development and implementation of a new cube-reconstruction algorithm that provides with an (almost) seeing-limited spatial resolution (FWHMPSF ∼ 1.0 arcsec). To illustrate the usability and quality of the data, we extracted two aperture spectra for each galaxy (central 1.5 arcsec and fully integrated), and analyse them using pyFIT3D. We obtain a set of observational and physical properties of both the stellar populations and the ionized gas, that have been compared for the two apertures, exploring their distributions as a function of the stellar masses and morphologies of the galaxies, comparing with recent results in the literature.

Sánchez, S. F. (ORCID:0000000164449307)↗

A near magnetic-to-kinetic energy equipartition flare from the relativistic jet in AO 0235 + 164 during 2013–2019

ABSTRACT We present the multiwavelength flaring activity of the blazar AO 0235 + 164 during its recent active period from 2013 to 2019. From a discrete correlation function analysis, we find a significant ($\gt 95~{{\ \rm per\ cent}}$) correlation between radio and gamma-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz Very Long Baseline Array data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of δvar = 28.5 ± 8.4, a bulk Lorentz factor of $\Gamma =16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $\theta _{\rm v}=1.42^{+1.07}_{-0.52}\textrm {~degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of 10−70 GHz and the peak flux density varying in the range of 0.7−4 Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}\textrm {~mG}$ at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}\textrm {~mG}$ calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235 + 164 suggests that it did not significantly deviate from equipartition during its recent flaring activity.

Astronomy & Astrophysics↗

Mixing bispectrum multipoles under geometric distortions

We derive general expressions for how the Alcock–Paczynski distortions affect the power spectrum and the bispectrum of cosmological fields. We compute explicit formulas for the mixing coefficients of bispectrum multipoles in the linear approximation. The leading-order effect for the bispectrum is the uniform dilation of all three wavevectors. The mixing coefficients depend on the shape of the bispectrum triplet. Our results for the bispectrum multipoles are framed in terms of the ‘natural’ basis of the lengths of three wavevectors but can be easily generalized for other bases and reduction schemes. Our validation tests confirm that the linear approximation is extremely accurate for all power spectrum multipoles. The linear approximation is accurate for the bispectrum monopole but results in sub-per cent level inaccuracies for the bispectrum quadrupole and fails for the bispectrum hexadecapole. Our results can be used to simplify the analysis of the bispectrum from galaxy surveys, especially the measurement of the baryon acoustic oscillation peak position. They can be used to replace numeric schemes with exact analytical formulae.

79 ASTRONOMY AND ASTROPHYSICS↗

Robust cosmological inference from non-linear scales with k -th nearest neighbour statistics

ABSTRACT We present the methodology for deriving accurate and reliable cosmological constraints from non-linear scales ($\lt 50\, h^{-1}$ Mpc) with k-th nearest neighbour (kNN) statistics. We detail our methods for choosing robust minimum scale cuts and validating galaxy–halo connection models. Using cross-validation, we identify the galaxy–halo model that ensures both good fits and unbiased predictions across diverse summary statistics. We demonstrate that we can model kNNs effectively down to transverse scales of $r_{\rm p}\sim 3\, h^{-1}$ Mpc and achieve precise and unbiased constraints on the matter density and clustering amplitude, leading to a 2 per cent constraint on σ8. Our simulation-based model pipeline is resilient to varied model systematics, spanning simulation codes, halo finding, and cosmology priors. We demonstrate the effectiveness of this approach through an application to the Beyond-2p mock challenge. We propose further explorations to test more complex galaxy–halo connection models and tackle potential observational systematics.

79 ASTRONOMY AND ASTROPHYSICS↗

The DESI One-Percent Survey: exploring a generalized SHAM for multiple tracers with the UNIT simulation

We perform SubHalo Abundance Matching (SHAM) studies on UNIT simulations with {σ, V ceil , v smear }-SHAM and {σ, V ceil , $f$ sat }-SHAM. They are designed to reproduce the clustering on 5–30 h -1 Mpc of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasi-stellar objects (QSOs) at 0.4 < z < 3.5 from DESI (Dark Energy Spectroscopic Instrument) One Percent Survey. V ceil is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. v smear models the clustering effect of redshift uncertainties, providing measurements consistent with those from repeat observations. A free satellite fraction $f$ sat is necessary to reproduce the clustering of ELGs. We find ELGs present a more complex galaxy–halo mass relation than LRGs reflected in their weak constraints on σ. LRGs, QSOs, and ELGs show increasing V ceil values, corresponding to the massive galaxy incompleteness of LRGs, the quenched star formation of ELGs and the quenched black hole accretion of QSOs. For LRGs, a Gaussian v smear presents a better profile for subsamples at redshift bins than a Lorentzian profile used for other tracers. The impact of the statistical redshift uncertainty on ELG clustering is negligible. The best-fitting satellite fraction for DESI ELGs is around 4 per cent, lower than previous estimations for ELGs. The mean halo mass log 10 ($\langle$M vir $\rangle$) in h -1 M ⊙ for LRGs, ELGs, and QSOs are 13.16 ± 0.01, 11.90 ± 0.06, and 12.66 ± 0.45, respectively. Our generalized SHAM algorithms facilitate the production of multitracer galaxy mocks for cosmological tests.

79 ASTRONOMY AND ASTROPHYSICS↗

Pulse profile modelling of thermonuclear burst oscillations – II. Handling variability

ABSTRACT Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii, and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 106 source counts, this yields uncertainties of approximately 10 per cent for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts.

Astronomy & Astrophysics↗

The impact of halo concentration on the Sunyaev Zel’dovich effect signal from massive galaxy clusters

ABSTRACT The Sunyaev Zel’dovich (SZ) effect is sensitive to the pressure of ionized gas inside galaxy clusters, which is in turn controlled largely by the gravitational potential of the cluster. Changing the concentration parameter describing the cluster mass distribution impacts the gravitational potential and thus the cluster SZ signal, with implications for cosmological and other analyses of SZ-selected clusters. We investigate the concentration–SZ relation in theory and simulations. We find that the impact of concentration on the inner SZ profile (R ≲ 0.75R200c) can be captured with standard polytropic gas models. However, we find that such models do a poor job of reproducing the outer SZ profiles (R ≳ 0.75R200c) and the relation between the integrated SZ signal, Y, and concentration. This disagreement results from a sharp truncation of the gas pressure profile near the splashback radius, likely caused by virial shocks. We develop a simple description of the truncation that leads to a good match with the simulated SZ profiles out to several R200c for clusters of varying mass and concentration, and that also accurately predicts the concentration–Y relationship. Finally, we determine how inference of the linear bias parameter and splashback radius for SZ-selected clusters can be biased by ignoring the concentration dependence of the SZ signal, finding that bias to the former is essentially negligible, while bias to the latter can be as much as 2 per cent.

Astronomy & 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↗

Hickson-like compact groups inhabiting different environments

ABSTRACT Although compact groups (CGs) of galaxies have been envisioned as isolated extremely dense structures in the Universe, it is accepted today that many of them could be not as isolated as thought. In this work, we study Hickson-like CGs identified in the Sloan Digital Sky Survey Data Release 16 to analyse these systems and their galaxies when embedded in different cosmological structures. To achieve this goal, we identify several cosmological structures where CGs can reside: nodes of filaments, loose groups, filaments, and cosmic voids. Our results indicate that 45 per cent of CGs do not reside in any of these structures, i.e. they can be considered non-embedded or isolated systems. Most of the embedded CGs are found inhabiting loose groups and nodes, while there are almost no CGs residing well inside cosmic voids. Some physical properties of CGs vary depending on the environment they inhabit. CGs in nodes show the largest velocity dispersions, the brightest absolute magnitude of the first-ranked galaxy, and the smallest crossing times, while the opposite occurs in non-embedded CGs. When comparing galaxies in all the environments and galaxies in CGs, CGs show the highest fractions of red/early-type galaxy members in most of the absolute magnitudes ranges. The variation between galaxies in CGs inhabiting one or another environment is not as significant as the differences caused by belonging or not to a CG. Our results suggest a plausible scenario for galaxy evolution in CGs in which both large-scale and local environments play essential roles.

Astronomy & Astrophysics↗

Late growth of early-type galaxies in low-z massive clusters

ABSTRACT We study a sample of 936 early-type galaxies (ETGs) located in 48 low-z regular galaxy clusters with M200 ≥ 1014 M⊙ at z < 0.1. We examine variations in the concentration index, radius, and colour gradient of ETGs as a function of their stellar mass and loci in the projected phase space (PPS) of the clusters. We aim to understand the environmental influence on the growth of ETGs according to the time since infall into their host clusters. Our analysis indicates a significant change in the behaviour of the concentration index C and colour gradient around $M_{\ast } \approx 2\times 10^{11} ~M_\odot \equiv \tilde{M}_{\ast }$. Objects less massive than $\tilde{M}_{\ast }$ present a slight growth of C with M*, with negative and approximately constant colour gradients in all regions of the PPS. Objects more massive than $\tilde{M}_{\ast }$ present a slight decrease of C with M*, with colour gradients becoming less negative and approaching zero. We also find that objects more massive than $\tilde{M}_{\ast }$, in all PPS regions, have smaller R90 for a given R50, suggesting a smaller external growth in these objects or even a shrinkage possibly due to tidal stripping. Finally, we estimate different dark matter fractions for galaxies in different regions of the PPS, with the ancient satellites having the largest fractions, fDM ≈ 65 per cent. These results favour a scenario where cluster ETGs experience environmental influence the longer they remain and the deeper into the gravitational potential they lie, indicating a combination of tidal stripping + harassment, which predominate during infall, followed by mergers + feedback effects affecting the late growth of ancient satellites and BCGs.

Astronomy & Astrophysics↗

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO’s and Advanced Virgo’s third observing run

ABSTRACT We describe a search for gravitational waves from compact binaries with at least one component with mass $0.2$–$1.0 \, \mathrm{M}_\odot$ and mass ratio q ≥ 0.1 in Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo data collected between 2019 November 1, 15:00 utc and 2020 March 27, 17:00 utc. No signals were detected. The most significant candidate has a false alarm rate of $0.2 \, \rm {yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO’s and Advanced Virgo’s third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim \, 0.6$ (at 90 per cent confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions, we are unable to rule out fPBH = 1. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound fDBH < 10−5 on the fraction of atomic dark matter collapsed into black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate model of the projected velocity distribution of galaxies in dark matter haloes

ABSTRACT We present a per cent-level accurate model of the line-of-sight velocity distribution of galaxies around dark matter haloes as a function of projected radius and halo mass. The model is developed and tested using synthetic galaxy catalogues generated with the UniverseMachine run on the Multi-Dark Planck 2 N-body simulations. The model decomposes the galaxies around a cluster into three kinematically distinct classes: orbiting, infalling, and interloping galaxies. We demonstrate that: (1) we can statistically distinguish between these three types of galaxies using only projected line-of-sight velocity information; (2) the halo edge radius inferred from the line-of-sight velocity dispersion is an excellent proxy for the three-dimensional halo edge radius; and (3) we can accurately recover the full velocity dispersion profile for each of the three populations of galaxies. Importantly, the velocity dispersion profiles of the orbiting and infalling galaxies contain five independent parameters – three distinct radial scales and two velocity dispersion amplitudes – each of which is correlated with mass. Thus, the velocity dispersion profile of galaxy clusters has inherent redundancies that allow us to perform non-trivial systematics checks from a single data set. We discuss several potential applications of our new model for detecting the edge radius and constraining cosmology and astrophysics using upcoming spectroscopic surveys.

Astronomy & Astrophysics↗

A better way to define dark matter haloes

ABSTRACT Dark matter haloes have long been recognized as one of the fundamental building blocks of large-scale structure formation models. Despite their importance – or perhaps because of it! – halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, effectively unique, and parameter-free: ‘A dark matter halo is comprised of the collection of particles orbiting in their own self-generated potential’. This definition is enabled by the fact that, even with as few as ≈300 particles per halo, nearly every particle in the vicinity of a halo can be uniquely classified as either orbiting or infalling based on its dynamical history. For brevity, we refer to haloes selected in this way as physical haloes. We demonstrate that (1) the mass function of physical haloes is Press–Schechter, provided the critical threshold for collapse is allowed to vary slowly with peak height; and (2) the peak-background split prediction of the clustering amplitude of physical haloes is statistically consistent with the simulation data, with accuracy no worse than ≈5 per cent.

Astronomy & Astrophysics↗

Forecasting constraints on the mean free path of ionizing photons at z ≥ 5.4 from the Lyman-α forest flux autocorrelation function

ABSTRACT Fluctuations in Lyman-α (Ly α) forest transmission towards high-z quasars are partially sourced from spatial fluctuations in the ultraviolet background, the level of which are set by the mean free path of ionizing photons (λmfp). The autocorrelation function of Ly α forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to λmfp. Recent measurements at z ∼ 6 suggest a rapid evolution of λmfp at z > 5.0 which would leave a signature in the evolution of the autocorrelation function. For this forecast, we model mock Ly α forest data with properties similar to the XQR-30 extended data set at 5.4 ≤ z ≤ 6.0. At each z, we investigate 100 mock data sets and an ideal case where mock data matches model values of the autocorrelation function. For ideal data with λmfp = 9.0 cMpc at z = 6.0, we recover $\lambda _{\text{mfp}}=12^{+6}_{-3}$ cMpc. This precision is comparable to direct measurements of λmfp from the stacking of quasar spectra beyond the Lyman limit. Hypothetical high-resolution data leads to a $\sim 40~{{\ \rm per\ cent}}$ reduction in the error bars over all z. The distribution of mock values of the autocorrelation function in this work is highly non-Gaussian for high-z, which should caution work with other statistics of the high-z Ly α forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements.

Astronomy & Astrophysics↗

DAHe white dwarfs from the DESI Survey

A new class of white dwarfs, dubbed DAHe, that present Zeeman-split Balmer lines in emission has recently emerged. However, the physical origin of these emission lines remains unclear. Here, we present here a sample of 21 newly identified DAHe systems and determine magnetic field strengths and (for a subset) periods that span the ranges of ≃6.5–147 MG and ≃0.4–36 h, respectively. All but four of these systems were identified from the Dark Energy Spectroscopic Instrument survey sample of more than 47 000 white dwarf candidates observed during its first year of observations. We present detailed analysis of the new DAHe WD J161634.36+541011.51 with a spin period of 95.3 min, which exhibits an anticorrelation between broad-band flux and Balmer line strength that is typically observed for this class of systems. All DAHe systems cluster closely on the Gaia Hertzsprung–Russell diagram where they represent ≃1 per cent of white dwarfs within that region. This grouping further solidifies their unexplained emergence at relatively late cooling times and we discuss this in context of current formation theories. Nine of the new DAHe systems are identifiable from Sloan Digital Sky Survey spectra of white dwarfs that had been previously classified as featureless DC-type systems. We suggest high-S/N (signal-to-noise ratios), unbiased observations of DCs as a possible route for discovering additional DAHe systems.

79 ASTRONOMY AND ASTROPHYSICS↗