Search NASA⌕ Search

SEARCH · Search NASA

Results for “ASTRONOMY AND ASTROPHYSICS”

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 343 records · Page 19

Discovery and characterization of five new eclipsing AM CVn systems

ABSTRACT AM CVn systems are ultra-compact, hydrogen-depleted, and helium-rich, accreting binaries with degenerate or semidegenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 min. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) light curves of white dwarfs selected using Gaia parallaxes. We obtained phase-resolved spectroscopy to confirm that all systems are AM CVn binaries, and we obtained high-speed photometry to confirm the eclipse and characterize the systems. The spectra show double-peaked H e lines but also show metals, including K and Zn, elements that have never been detected in AM CVn systems before. By modelling the high-speed photometry, we measured the mass and radius of the donor star, potentially constraining the evolutionary channel that formed these AM CVn systems. We determined that the average mass of the accreting white dwarf is ≈0.8 M⊙, and that the white dwarfs in long-period systems are hotter than predicted by recently updated theoretical models. The donors have a high entropy and are a factor of ≈2 more massive compared to zero-entropy donors at the same orbital period. The large donor radius is most consistent with H e-star progenitors, although the observed spectral features seem to contradict this. The discovery of five new eclipsing AM CVn systems is consistent with the known observed AM CVn space density and estimated ZTF recovery efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗

The picasso map-making code: application to a simulation of the QUIJOTE northern sky survey

ABSTRACT Map-making is an important step for the data analysis of cosmic microwave background (CMB) experiments. It consists of converting the data, which are typically a long, complex, and noisy collection of measurements, into a map, which is an image of the observed sky. We present in this paper a new map-making code named picasso (Polarization and Intensity CArtographer for Scanned Sky Observations), which was implemented to construct intensity and polarization maps from the Multi Frequency Instrument (MFI) of the QUIJOTE (Q-U-I Joint TEnerife) CMB polarization experiment. picasso is based on the destriping algorithm, and is suited to address specific issues of ground-based microwave observations, with a technique that allows the fit of a template function in the time domain, during the map-making step. This paper describes the picasso code, validating it with simulations and assessing its performance. For this purpose, we produced realistic simulations of the QUIJOTE-MFI survey of the northern sky (approximately ∼20 000 deg2), and analysed the reconstructed maps with picasso, using real and harmonic space statistics. We show that, for this sky area, picasso is able to reconstruct, with high fidelity, the injected signal, recovering all the scales with ℓ > 10 in TT, EE, and BB. The signal error is better than 0.001 per cent at 20 < ℓ < 200. Finally, we validated some of the methods that will be applied to the real wide-survey data, like the detection of the CMB anisotropies via cross-correlation analyses. Despite that the implementation of picasso is specific for QUIJOTE-MFI data, it could be adapted to other experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

The galaxy–halo connection of emission-line galaxies in IllustrisTNG

ABSTRACT We employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the halo occupation distribution (HOD) and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at z ∼ 1. Such galaxies are key targets for current and upcoming cosmological surveys. We select model galaxies through cuts in colour–colour space allowing for a direct comparison with the Extended Baryon Oscillation Spectroscopic Survey and the Dark Energy Spectroscopic Instrument (DESI) surveys and then compare them with galaxies selected based on specific star formation rate (sSFR) and stellar mass. We demonstrate that the ELG populations are twice more likely to reside in lower density regions (sheets) compared with the mass-selected populations and twice less likely to occupy the densest regions of the cosmic web (knots). We also show that the colour-selected and sSFR-selected ELGs exhibit very similar occupation and clustering statistics, finding that the agreement is best for lower redshifts. In contrast with the mass-selected sample, the occupation of haloes by a central ELG peaks at ∼20 per cent. We furthermore explore the dependence of the HOD and the autocorrelation on environment, noticing that at fixed halo mass, galaxies in high-density regions cluster about 10 times more strongly than low-density ones. This result suggests that we should model carefully the galaxy–halo relation and implement assembly bias effects into our models (estimated at ∼4 per cent of the clustering of the DESI colour-selected sample at z = 0.8). Finally, we apply a simple mock recipe to recover the clustering on large scales (r ≳ 1 Mpc h−1) to within 1 per cent by augmenting the HOD model with an environment dependence, demonstrating the power of adopting flexible population models.

79 ASTRONOMY AND ASTROPHYSICS↗

Variable stars in Local Group galaxies – V. The fast and early evolution of the low-mass Eridanus II dSph galaxy

ABSTRACT We present a detailed study of the variable star population of Eridanus II (Eri II), an ultra-faint dwarf galaxy that lies close to the Milky Way virial radius. We analyse multi-epoch g, r, i ground-based data from Goodman and the Dark Energy Camera, plus F475W, F606W, F814W space data from the Advanced Camera for Surveys. We report the detection of 67 RR Lyrae (RRL) stars and 2 Anomalous Cepheids, most of them new discoveries. With the RRL stars, we measure the distance modulus of Eri II, μ0 = 22.84 ± 0.05 mag (D⊙ = 370 ± 9 kpc), and derive a metallicity spread of 0.3 dex (0.2 dex intrinsic). The colour distribution of the horizontal branch (HB) and the period distribution of the RRL stars can be nicely reproduced by a combination of two stellar models of [Fe/H] = (−2.62, −2.14). The overall low metallicity is consistent with the red giant branch bump location, 0.65 mag brighter than the HB. These results are in agreement with previous spectroscopic studies. The more metal-rich RRL and the RRab stars have greater central concentration than the more metal-poor RRL and the RRc stars that are mainly located outside ∼1 rh. This is similar to what is found in larger dwarf galaxies such as Sculptor, and in agreement with an outside-in galaxy formation scenario. This is remarkable in such a faint dwarf galaxy with an apparently single and extremely short (<1 Gyr) star formation burst. Finally, we have derived new and independent structural parameters for Eri II and its star cluster using our new data that are in very good agreement with previous estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗

Rotation curves and scaling relations of extremely massive spiral galaxies

We study the kinematics and scaling relations of a sample of 43 giant spiral galaxies that have stellar masses exceeding $10^{11} \, {\rm M}_\odot$ and optical discs up to 80 kpc in radius. We use a hybrid 3D–1D approach to fit 3D kinematic models to long-slit observations of the H α-$\rm{[N\, \small {II}]}$ emission lines and we obtain robust rotation curves of these massive systems. We find that all galaxies in our sample seem to reach a flat part of the rotation curve within the outermost optical radius. We use the derived kinematics to study the high-mass end of the two most important scaling relations for spiral galaxies: the stellar/baryonic mass Tully–Fisher relation and the Fall (mass-angular momentum) relation. All galaxies in our sample, with the possible exception of the two fastest rotators, lie comfortably on both these scaling relations determined at lower masses, without any evident break or bend at the high-mass regime. When we combine our high-mass sample with low-mass data from the Spitzer Photometry & Accurate Rotation Curves catalogue, we find a slope of α = 4.25 ± 0.19 for the stellar Tully–Fisher relation and a slope of γ = 0.64 ± 0.11 for the Fall relation. Our results indicate that most, if not all, of these rare, giant spiral galaxies are scaled up versions of less massive discs and that spiral galaxies are a self-similar population of objects up to the very high-mass end.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling type 1 quasar colours in the era of Rubin and Euclid

We construct a parametric SED model which is able to reproduce the average observed SDSS–UKIDSS–WISE quasar colours to within one-tenth of a magnitude across a wide range of redshift (0 < z < 5) and luminosity (–22 > M i > –29). This model is shown to provide accurate predictions for the colours of known quasars which are less luminous than those used to calibrate the model parameters, and also those at higher redshifts z > 5. Using a single parameter, the model encapsulates an up-to-date understanding of the intra-population variance in the rest-frame ultraviolet and optical emission lines of luminous quasars. At fixed redshift, there are systematic changes in the average quasar colours with apparent i-band magnitude, which we find to be well explained by the contribution from the host galaxy and our parametrization of the emission-line properties. By including redshift as an additional free parameter, the model could be used to provide photometric redshifts for individual objects. For the population as a whole we find that the average emission line and host-galaxy contributions can be well described by simple functions of luminosity which account for the observed changes in the average quasar colours across 18.1 < i AB < 21.5. Finally, we use these trends to provide predictions for quasar colours at the luminosities and redshifts which will be probed by the Rubin Observatory LSST and ESA-Euclid wide survey. The model code is applicable to a wide range of upcoming photometric and spectroscopic surveys, and is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessing the sources of reionization: a spectroscopic case study of a 30× lensed galaxy at z ~ 5 with Ly α , C IV, Mg ιι, and [Ne ιιι]

We present a detailed spectroscopic analysis of a galaxy at z ≃ 4.88 that is, by chance, magnified ~30× by gravitational lensing. Only three sources at z ≳ 5 are known with such high magnification. This particular source has been shown to exhibit widespread, high equivalent width ${{\rm C\, \small {IV}}}\ \lambda 1549\, \mathring{\rm A}$ emission, implying it is a unique example of a metal-poor galaxy with a hard radiation field, likely representing the galaxy population responsible for cosmic reionization. Using ultraviolet (UV) nebular line ratio diagnostics, Very Large Telescope (VLT)/X-shooter observations rule out strong active galactic nuclei (AGN) activity, indicating a stellar origin of the hard radiation field instead. We present a new detection of ${[{\rm Ne\, \small {III}}]}\ \lambda 3870\, \mathring{\rm A}$ and use the [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] line ratio to constrain the ionization parameter and gas-phase metallicity. Closely related to the commonly used [${\rm O\, \small {III}}$]/[${\rm O\, \small {II}}$] ratio, our [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] measurement shows this source is similar to local ‘Green Pea’ galaxies and Lyman-continuum leakers. It furthermore suggests this galaxy is more metal poor than expected from the fundamental metallicity relation, possibly as a consequence of excess gas accretion diluting the metallicity. Finally, we present the highest redshift detection of ${{\rm Mg\, \small {II}}}\ \lambda 2796\, \mathring{\rm A}$, observed at high equivalent width in emission, in contrast to more evolved systems predominantly exhibiting ${\rm Mg\, \small {II}}$ absorption. Strong ${\rm Mg\, \small {II}}$ emission has been observed in most z ~ 0 Lyman-continuum leakers known and has recently been proposed as an indirect tracer of escaping ionizing radiation. In conclusion, this strongly lensed galaxy, observed just $300\, \mathrm{Myr}$ after reionization ends, enables testing of observational diagnostics proposed to constrain the physical properties of distant galaxies in the James Webb Space Telescope (JWST)/Extremely Large Telescope (ELT) era.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining the Milky Way’s ultraviolet-to-infrared SED with Gaussian process regression

ABSTRACT Improving our knowledge of global Milky Way (MW) properties is critical for connecting the detailed measurements only possible from within our Galaxy to our understanding of the broader galaxy population. We here train Gaussian process regression (GPR) models on SDSS (Sloan Digital Sky Survey) galaxies to map from galaxy properties (stellar mass, apparent axial ratio, star formation rate, bulge-to-total ratio, disc scale length, and bar vote fraction) to ultraviolet (UV; GALEX FUV/NUV), optical (SDSS ugriz), and infrared (IR; 2MASS JHKs and WISE W1/W2/W3/W4) fluxes and uncertainties. With these models, we estimate the photometric properties of the MW, resulting in a full UV-to-IR spectral energy distribution (SED) as it would be measured externally, viewed face-on. We confirm that the MW lies in the green valley in optical diagnostic diagrams, but show for the first time that the MW is in the star-forming region in standard UV and IR diagnostics – characteristic of the population of red spiral galaxies. Although our GPR method predicts one band at a time, the resulting MW UV–IR SED is consistent with SEDs of local spirals with characteristics broadly similar to the MW, suggesting that these independent predictions can be combined reliably. Our UV–IR SED will be invaluable for reconstructing the MW’s star formation history using the same tools employed for external galaxies, allowing comparisons of results from in situ measurements to those from the methods used for extragalactic objects.

79 ASTRONOMY AND ASTROPHYSICS↗

A panoramic view of the Local Group dwarf galaxy NGC 6822

ABSTRACT We present a panoramic survey of the isolated Local Group dwarf irregular galaxy NGC 6822. Our photometry reaches ∼2–3 mag deeper than most previous studies and spans the widest area around the dwarf compared to any prior work. We observe no stellar overdensities in the outskirts of NGC 6822 to V ∼ 30 mag arcsec−2 and a projected radius of 16.5 kpc. This indicates that NGC 6822 has not experienced any recent interaction with a companion galaxy, despite previous suggestions to the contrary. Similarly, we find no evidence for any dwarf satellites of NGC 6822 to a limiting luminosity MV ≈ −5. NGC 6822 contains a disc of H i gas and young stars, oriented at ∼60○ to an extended spheroid composed of old stellar populations. We observe no correlation between the distribution of young stars and spheroid members. Our imaging allows us to trace the spheroid to nearly 11 kpc along its major axis, commensurate with the extent of the NGC 6822 globular cluster system. We find that the spheroid becomes increasingly flattened at larger radii, and its position angle twists by up to 40○. We use Gaia EDR3 astrometry to measure a proper motion for NGC 6822, and then sample its orbital parameter space. While this galaxy has spent the majority of its life in isolation, we find that it likely passed within the virial radius of the Milky Way ∼3–4 Gyr ago. This may explain the apparent flattening and twisting observed in the outskirts of its spheroid.

79 ASTRONOMY AND ASTROPHYSICS↗

A catalogue of white dwarfs in Gaia EDR3

ABSTRACT We present a catalogue of white dwarf candidates selected from Gaia Early Data Release 3 (EDR3). We applied several selection criteria in absolute magnitude, colour, and Gaia quality flags to remove objects with unreliable measurements while preserving most stars compatible with the white dwarf locus in the Hertzsprung–Russell diagram. We then used a sample of over 30 000 spectroscopically confirmed white dwarfs and contaminants from the Sloan Digital Sky Survey (SDSS) to map the distribution of these objects in the Gaia absolute magnitude–colour space. Finally, we adopt the same method presented in our previous work on Gaia Data Release 2 (DR2) to calculate a probability of being a white dwarf (PWD) for ≃1.3 million sources that passed our quality selection. The PWD values can be used to select a sample of ${\simeq} 359\,000$ high-confidence white dwarf candidates. We calculated stellar parameters (effective temperature, surface gravity, and mass) for all these stars by fitting Gaia astrometry and photometry with synthetic pure-H, pure-He, and mixed H–He atmospheric models. We estimate an upper limit of 93 per cent for the overall completeness of our catalogue for white dwarfs with G ≤ 20 mag and effective temperature (Teff) > 7000 K, at high Galactic latitudes (|b| > 20°). Alongside the main catalogue we include a reduced proper motion extension containing ${\simeq} 10\,200$ white dwarf candidates with unreliable parallax measurements that could, however, be identified on the basis of their proper motion. We also performed a cross-match of our catalogues with SDSS Data Release 16 (DR16) spectroscopy and provide spectral classification based on visual inspection for all resulting matches.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating highly eccentric common envelope jet supernova impostors

We conduct three-dimensional hydrodynamical simulations of eccentric common envelope jet supernova (CEJSN) impostors, i.e. a neutron star that crosses through the envelope of a red supergiant star on a highly eccentric orbit and launches jets as it accretes mass from the envelope. Because of numerical limitations, we apply a simple prescription where we inject the assumed jets’ power into two opposite conical regions inside the envelope. We find the outflow morphology to be very complicated, clumpy, and non-spherical, having a large-scale symmetry only about the equatorial plane. The outflow morphology can substantially differ between simulations that differ by their jets’ power. We estimate by simple means the light curve to be very bumpy, to have a rise time of one to a few months, and to slowly decay in about a year to several years. These eccentric CEJSN impostors will be classified as ‘gap’ objects, i.e. having a luminosity between those of classical novae and typical supernovae (termed also ILOTs for intermediate luminosity optical transients). We strengthen a previous conclusion that CEJSN impostors might account for some peculiar ILOTs, in particular those that might repeat over time-scales of months to years.

79 ASTRONOMY AND ASTROPHYSICS↗

The antesonic condition for the explosion of core-collapse supernovae – II. Rotation and turbulence

In the problem of steady free fall on to a standing shockwave around a central mass, the ‘antesonic’ condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\le 3/16$, where cT is the isothermal sound speed in the subsonic post-shock flow, and vesc is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and the shock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models of accretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a ~5 percent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by ~20 to 40 percent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multidimensional simulations explode more easily than their one-dimensional counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of void environment on AGN

ABSTRACT We study the population of active galaxies in void environment in the Sloan Digital Sky Survey. We use optical spectroscopic information to analyse characteristics of the emission lines of galaxies, accomplished by WHAN and BPT diagrams. Also, we study Wide-field Infrared Survey Explorer(WISE) mid-IR colours to assess active galactic nucleus (AGN) activity. We investigate these different AGN classification schemes, both optical and mid-IR, and their dependence on the spatial location with respect to the void centres. To this end, we define three regions: void, the spherical region defined by voidcentric distance relative to void radius (distance/rvoid) smaller than 0.8, comprising overdensities lesser than −0.9, an intermediate/transition shell region (namely, void-wall) 0.8 < distance/rvoid < 1.2, and a region sufficiently distant from voids, the field: distance/rvoid > 2. We find statistical evidence for a larger fraction of AGN and star-forming galaxies in the void region, regardless of the classification scheme addressed (either BPT, WHAN, or WISE). Moreover, we obtain a significantly stronger nuclear activity in voids compared to the field. We find an unusually large fraction of the most massive black holes undergoing strong accretion when their host galaxies reside in voids. Our results suggest a strong influence of the void environment on AGN mechanisms associated with galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Intrinsic alignments of galaxies around cosmic voids

ABSTRACT The intrinsic alignments of galaxies, i.e. the correlation between galaxy shapes and their environment, are a major source of contamination for weak gravitational lensing surveys. Most studies of intrinsic alignments have so far focused on measuring and modelling the correlations of luminous red galaxies with galaxy positions or the filaments of the cosmic web. In this work, we investigate alignments around cosmic voids. We measure the intrinsic alignments of luminous red galaxies detected by the Sloan Digital Sky Survey around a sample of voids constructed from those same tracers and with radii in the ranges: [20–30; 30–40; 40–50] h−1 Mpc and in the redshift range z = 0.4−0.8. We present fits to the measurements based on a linear model at large scales, and on a new model based on the void density profile inside the void and in its neighbourhood. We constrain the free scaling amplitude of our model at small scales, finding no significant alignment at 1σ for either sample. We observe a deviation from the null hypothesis, at large scales, of 2σ for voids with radii between 20 and 30 h−1 Mpc, and 1.5σ for voids with radii between 30 and 40 h−1 Mpc and constrain the amplitude of the model on these scales. We find no significant deviation at 1σ for larger voids. Our work is a first attempt at detecting intrinsic alignments of galaxy shapes around voids and provides a useful framework for their mitigation in future void lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗

The dark matter haloes of HI selected galaxies

We present the neutral hydrogen mass (M HI) function (HIMF) and velocity width (w 50 ) function (HIWF) based on a sample of 7857 galaxies from the 40 per cent data release of the ALFALFA survey (α.40). The low mass (velocity width) end of the HIMF (HIWF) is dominated by the blue population of galaxies whereas the red population dominates the HIMF (HIWF) at the high mass (velocity width) end. We use a deconvolution method to estimate the HI rotational velocity (V rot ) functions (HIVF) from the HIWF for the total, red, and blue samples. The HIWF and HIVF for the red and blue samples are well separated at the knee of the function compared to their HIMFs. We then use recent stacking results from the ALFALFA survey to constrain the halo mass (M h ) function of HI-selected galaxies. This allows us to obtain various scaling relations between M HI –w 50 –V rot –M h , which we present. The M HI –M h relation has a steep slope ~2.10 at small masses and flattens to ~0.34 at masses larger than a transition halo mass, $\log _{10}(M_{\rm{ht}}h_{70}^2/M_{\odot })=10.62$. Our scaling relation is robust and consistent with a volume-limited sample of α.40. The M HI –M h relation is qualitatively similar to the M star –M h relation but the transition halo mass is smaller by ~1.4 dex compared to that of the M star –M h relation. Our results suggest that baryonic processes like heating and feedback in larger mass haloes suppress HI gas on a shorter time-scale compared to star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

$\rm {H\,\small {I}}$ constraints from the cross-correlation of eBOSS galaxies and Green Bank Telescope intensity maps

We present the joint analysis of Neutral Hydrogen ($\rm {H\,\small {I}}$) Intensity Mapping observations with three galaxy samples: the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples from the eBOSS survey, and the WiggleZ Dark Energy Survey sample. The $\rm {H\,\small {I}}$ intensity maps are Green Bank Telescope observations of the redshifted $21\rm cm$ emission on $100 \, {\rm deg}^2$ covering the redshift range 0.6 < z < 1.0. We process the data by separating and removing the foregrounds present in the radio frequencies with FastI ICA. We verify the quality of the foreground separation with mock realizations, and construct a transfer function to correct for the effects of foreground removal on the $\rm {H\,\small {I}}$ signal. We cross-correlate the cleaned $\rm {H\,\small {I}}$ data with the galaxy samples and study the overall amplitude as well as the scale dependence of the power spectrum. We also qualitatively compare our findings with the predictions by a semianalytical galaxy evolution simulation. The cross-correlations constrain the quantity $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm opt}}$ at an effective scale k eff , where $\Omega _\rm {H\,\small {I}}$ is the $\rm {H\,\small {I}}$ density fraction, $b_\rm {H\,\small {I}}$ is the $\rm {H\,\small {I}}$ bias, and $r_{\rm {H\,\small {I}},{\rm opt}}$ the galaxy–hydrogen correlation coefficient, which is dependent on the $\rm {H\,\small {I}}$ content of the optical galaxy sample. At $k_{\rm eff}=0.31 \, h\,{\rm Mpc^{-1}}$ we find $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm Wig}} = [0.58 \pm 0.09 \, {\rm (stat) \pm 0.05 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-WiggleZ, $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm ELG}} = [0.40 \pm 0.09 \, {\rm (stat) \pm 0.04 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-ELG, and $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm LRG}} = [0.35 \pm 0.08 \, {\rm (stat) \pm 0.03 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-LRG, at z ≃ 0.8. We also report results at $k_{\rm eff}=0.24$ and $k_{\rm eff}=0.48 \, h\,{\rm Mpc^{-1}}$. With little information on $\rm {H\,\small {I}}$ parameters beyond our local Universe, these are amongst the most precise constraints on neutral hydrogen density fluctuations in an underexplored redshift range.

79 ASTRONOMY AND ASTROPHYSICS↗

A selection function toolbox for subsets of astronomical catalogues

ABSTRACT Large catalogues are ubiquitous throughout astronomy, but most scientific analyses are carried out on smaller samples selected from these catalogues by chosen cuts on catalogued quantities. The selection function of that scientific sample – the probability that a star in the catalogue will satisfy these cuts and so make it into the sample – is thus unique to each scientific analysis. We have created a general framework that can flexibly estimate the selection function of a sample drawn from a catalogue in terms of position, magnitude and colour. Our method is unique in using the binomial likelihood and accounting for correlations in the selection function across position, magnitude, and colour using Gaussian processes and spherical harmonics. We have created a new open-source python package selectionfunctiontoolbox that implements this framework and used it to make three different estimates of the APOGEE DR16 red giant sample selection function, as a subset of 2MASS, with each estimate using an increasing amount of technical knowledge of the APOGEE targeting. In a companion paper, we applied our methodology to derive estimates of the astrometric and spectroscopic selection functions of Gaia EDR3. Our framework will make it trivial for astrophysicists to estimate the selection function that they should be using with the custom sample of stars that they have chosen to answer their scientific question.

79 ASTRONOMY AND ASTROPHYSICS↗