Photometric Bands Research Articles

Neural networks optimized by genetic algorithms in cosmology

The applications of artificial neural networks in the cosmological field have shone successfully during the past decade, this is due to their great ability of modeling large amounts of datasets and complex nonlinear functions. However, in some cases, their use still remains controversial because their ease of producing inaccurate results when the hyperparameters are not carefully selected. In this paper, to find the optimal combination of hyperparameters to artificial neural networks, we propose to take advantage of the genetic algorithms. As a proof of the concept, we analyze three different cosmological cases to test the performance of the architectures achieved with the genetic algorithms and compare them with the standard process, consisting of a grid with all possible configurations. First, we carry out a model-independent reconstruction of the distance modulus using a type Ia supernovae compilation. Second, the neural networks learn to infer the equation of state for the quintessence model, and finally with the data from a combined redshift catalog the neural networks predict the photometric redshift given six photometric bands (urgizy). We found that the genetic algorithms improve considerably the generation of the neural network architectures, which can ensure more confidence in their physical results because of the better performance in the metrics with respect to the grid method.

PHANGS–JWST First Results: Measuring Polycyclic Aromatic Hydrocarbon Properties across the Multiphase Interstellar Medium

Ratios of polycyclic aromatic hydrocarbon (PAH) vibrational bands are a promising tool for measuring the properties of the PAH population and their effect on star formation. The photometric bands of the MIRI and NIRCam instruments on JWST provide the opportunity to measure PAH emission features across entire galaxy disks at unprecedented resolution and sensitivity. Here we present the first results of this analysis in a sample of three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Based on the variations observed in the 3.3, 7.7, and 11.3 μm features, we infer changes to the average PAH size and ionization state across the different galaxy environments. High values of F335MPAH/F1130W and low values of F1130W/F770W are measured in H ii regions in all three galaxies. This suggests that these regions are populated by hotter PAHs, and/or that the PAH ionization fraction is larger. We see additional evidence of heating and/or changes in PAH size in regions with higher molecular gas content as well as increased ionization in regions with higher Hα intensity.

Data-driven Cosmology from Three-dimensional Light Cones

We present a data-driven technique to analyze multifrequency images from upcoming cosmological surveys mapping large sky area. Using full information from the data at the two-point level, our method can simultaneously constrain the large-scale structure (LSS), the spectra and redshift distribution of emitting sources, and the noise in the observed data without any prior assumptions beyond the homogeneity and isotropy of cosmological perturbations. In particular, the method does not rely on source detection or photometric or spectroscopic redshift estimates. Here, we present the formalism and demonstrate our technique with a mock observation from nine optical and near-infrared photometric bands. Our method can recover the input signal and noise without bias, and quantify the uncertainty on the constraints. Our technique provides a flexible framework to analyze the LSS observation traced by different types of sources, which has potential for wide application to current or future cosmological data sets such as SPHEREx, Rubin Observatory, Euclid, or the Nancy Grace Roman Space Telescope.

Near-infrared Polarization Charateristics of the Zodiacal Light Observed with DIRBE/COBE

We report near-infrared polarization of the zodiacal light (ZL) measured from space by the Diffuse Infrared Background Experiment (DIRBE) on board the Cosmic Background Explorer in photometric bands centered at 1.25, 2.2, and 3.5 μm. To constrain the physical properties of interplanetary dust, we use DIRBE Weekly Sky Maps to investigate the solar elongation (ϵ), ecliptic latitude (β), and wavelength (λ) dependence of ZL polarization. We find that the polarization of the ZL varies as a function of ϵ and β, consistent with observed polarization at λ = 550 nm. While the polarization dependence on wavelength at (ϵ, β) = (90°, 0°) is modest (increasing from 17.7% ± 0.2% at 1.25% μm to 21.0% ± 0.3% at 3.5 μm), the variation is more pronounced at the north ecliptic pole (23.1% ± 1.6%, 35.1% ± 2.0%, and 39.3% ± 2.1% at 1.25, 2.2, and 3.5 μm, respectively). The variation in ZL polarization with wavelength is not explained by either Rayleigh scattering or absorptive particles larger than 10 μm.

On the Posch ratio for irradiance in coastal waters and the high seas

The horizontal irradiance at the sea surface is an informative light pollution indicator to study Artificial Light at Night (ALAN) effects on marine biodiversity (e.g.: zooplankton diel vertical migration). The Posch ratio (PR) for the horizontal irradiance (that is, the ratio of the horizontal irradiance to the zenith radiance) is a useful tool for estimating the irradiance from easily available measurements of the zenith night sky brightness. The PR definition has already been generalized for any pair of linear radiance indicators in any pair of arbitrarily chosen photometric bands, and can also be applied to estimate e.g. the average sky radiance or the radiance at some elevation above the horizon as a function of the radiance in any other direction of the sky. The PR for a single light source depends on the distance from the source, its angular and spectral emission pattern, and the state of the atmosphere. The PR for any set of sources is a linear combination of the individual PRs that each one would produce separately, with weights that can be easily derived from the relative contribution of each source to the zenith radiance. Whereas in populated lands the ALAN PR varies relatively fast from one location to another, due to the particular spatial distribution of lights, in coastal waters and the high seas the light pollution PR is a smooth function of the distance to the shoreline, due to the progressive lack of neighboring sources and the absence of obstacles. In this work we present the fundamental equations of the model and an example of application for the waters surrounding the Iberian Peninsula, North Africa and the West Mediterranean islands.

Core-collapse supernovae in the Dark Energy Survey: luminosity functions and host galaxy demographics

ABSTRACT We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically confirmed supernovae over a redshift range 0.045 < z < 0.25. We fit the observed DES griz CCSN light curves and K-correct to produce rest-frame R-band light curves. We compare the sample with lower redshift CCSN samples from Zwicky Transient Facility (ZTF) and Lick Observatory Supernova Search (LOSS). Comparing luminosity functions, the DES and ZTF samples of SNe II are brighter than that of LOSS with significances of 3.0σ and 2.5σ, respectively. While this difference could be caused by redshift evolution in the luminosity function, simpler explanations such as differing levels of host extinction remain a possibility. We find that the host galaxies of SNe II in DES are on average bluer than in ZTF, despite having consistent stellar mass distributions. We consider a number of possibilities to explain this – including galaxy evolution with redshift, selection biases in either the DES or ZTF samples, and systematic differences due to the different photometric bands available – but find that none can easily reconcile the differences in host colour between the two samples and thus its cause remains uncertain.

Stellar ages, masses, extinctions, and orbital parameters based on spectroscopic parameters ofGaiaDR3

Context.Gaia’s third data release provides radial velocities for 33 million stars and spectroscopically derived atmospheric parameters for more than 5 million targets. When combined with the astrometric data, these allow us to derive orbital and stellar parameters that are key to understanding the stellar populations of the Milky Way and to perform Galactic archaeology.Aims.We used the calibrated atmospheric parameters, 2MASS andGaia-EDR3 photometry, and parallax-based distances to compute the ages, initial stellar masses, and reddenings for the stars with spectroscopic parameters. We also derived the orbits for all of the stars with measured radial velocities and astrometry, adopting two sets of line-of-sight distances from the literature.Methods.Four different sets of ages, masses, and absolute magnitudes in different photometric bands are obtained through an isochrone fitting method, considering different combinations of input parameters. The reddenings are obtained by comparing the observed colours with those obtained from the isochrone projection. Finally, the orbits are computed adopting an axisymmetric potential of the Galaxy.Results.Comparisons with reference catalogues of field and cluster stars suggest that reliable ages are obtained for stars younger than 9–10 Gyr when the estimated relative age uncertainty is < 50 per cent. For older stars, ages tend to be underestimated. The most reliable stellar type for age determination are turn-off stars, even when the input atmospheric parameters have large uncertainties. Ages for giants and main-sequence stars are retrieved with uncertainties of the order of 2 Gyr when extinction towards the star’s line of sight is less thanAV ≲ 2.5 mag.Conclusions.The catalogue of ages, initial stellar masses, reddenings, galactocentric positions and velocities, as well as the stellar orbital actions, eccentricities, apocentre, pericentre and maximum distance from the Galactic plane reached during their orbits, is made publicly available to be downloaded. With this catalogue, the full chemo-dynamical properties of the extended solar neighbourhood unfold and allow us to better identify the properties of the spiral arms, to parametrise the dynamical heating of the disc, and to thoroughly study the chemical enrichment of the Milky Way.

Photometric redshift estimation of strongly lensed galaxies

Context.Around 105strongly lensed galaxies are expected to be discovered with upcoming wide-field surveys such asEuclidand the LSST. Utilising these large samples to study the inner structure of lens galaxies requires source redshifts, which are needed to turn lens models into mass measurements. However, obtaining spectroscopic source redshifts for large samples of strong lenses is prohibitive with the current capacity of spectroscopic facilities.Aims.As an alternative to spectroscopy, we study the possibility of obtaining source photometric redshifts (photo-zs) for large samples of strong lenses. We pay particular attention to the problem of blending between the lens and the source light.Methods.Our strategy consists of deblending the source and lens light by simultaneously modelling the lens galaxy and the background source in all available photometric bands, and then feeding the derived source colours to a template-fitting photo-zalgorithm. We describe the lens and the source light with a Sérsic profile, and the lens mass with a singular isothermal ellipsoid. We first test our approach on a simulated sample of lenses. Then, we apply it to 23 real systems with broad-band photometry from the Hyper Suprime-Cam survey.Results.We identify the deviations of the lens light from a Sérsic profile and the contrast between the lens and source image as the main drivers of the source colour measurement error. Although the former is challenging to measure directly for real lenses, we find the latter to be sufficient for evaluating the accuracy of a measured source colour. We split the real sample based on the ratio Λ of the lens to source surface brightness measured at the image locations. In the Λ < 1 regime, the photo-zoutlier fraction is 20%, and the accuracy of photo-zestimation is limited by the performance of the template-fitting process. In the opposite regime, the photo-zoutlier fraction is 75%, and the errors from the source colour measurements dominate the photo-zuncertainty.Conclusions.Measuring source photo-zs for lenses with Λ < 1 poses no particular challenges compared to the isolated galaxy case. For systems with significant lens light contamination, however, improvements in the description of the surface brightness distribution of the lens are required: a single Sérsic model is not sufficiently accurate.

A Three-dimensional Map of the Milky Way Using 66,000 Mira Variable Stars

We study the three-dimensional structure of the Milky Way using 65,981 Mira variable stars discovered by the Optical Gravitational Lensing Experiment survey. The spatial distribution of the Mira stars is analyzed with a model containing three barred components that include the X-shaped boxy component in the Galactic center (GC) and an axisymmetric disk. We take into account the distance uncertainties by implementing the Bayesian hierarchical inference method. The distance to the GC is R 0 = 7.66 ± 0.01(stat.) ± 0.39(sys.) kpc, while the inclination of the major axis of the bulge to the Sun–GC line of sight is θ = 20.°2 ± 0.°6(stat.) ± 0.°7(sys.). We present, for the first time, a detailed three-dimensional map of the Milky Way composed of young and intermediate-age stellar populations. Our analysis provides independent evidence for both the X-shaped bulge component and the flaring disk (being plausibly warped). We provide the complete data set of properties of Miras that were used for calculations in this work. The table includes mean brightness and amplitudes in nine photometric bands (covering a range of wavelengths from 0.5 to 12 μm), photometric chemical type, estimated extinction, and calculated distance with its uncertainty for each Mira variable. The median distance accuracy to a Mira star is at the level of 6.6%.

Merger identification through photometric bands, colours, and their errors

Aims. We present the application of a fully connected neural network (NN) for galaxy merger identification using exclusively photometric information. Our purpose is not only to test the method’s efficiency, but also to understand what merger properties the NN can learn and what their physical interpretation is. Methods. We created a class-balanced training dataset of 5860 galaxies split into mergers and non-mergers. The galaxy observations came from SDSS DR6 and were visually identified in Galaxy Zoo. The 2930 mergers were selected from known SDSS mergers and the respective non-mergers were the closest match in both redshift and r magnitude. The NN architecture was built by testing a different number of layers with different sizes and variations of the dropout rate. We compared input spaces constructed using: the five SDSS filters: u, g, r, i, and z; combinations of bands, colours, and their errors; six magnitude types; and variations of input normalization. Results. We find that the fibre magnitude errors contribute the most to the training accuracy. Studying the parameters from which they are calculated, we show that the input space built from the sky error background in the five SDSS bands alone leads to 92.64 ± 0.15% training accuracy. We also find that the input normalization, that is to say, how the data are presented to the NN, has a significant effect on the training performance. Conclusions. We conclude that, from all the SDSS photometric information, the sky error background is the most sensitive to merging processes. This finding is supported by an analysis of its five-band feature space by means of data visualization. Moreover, studying the plane of the g and r sky error bands shows that a decision boundary line is enough to achieve an accuracy of 91.59%.

Classifying Astronomical Transients Using Only Host Galaxy Photometry

The Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory will discover tens of thousands of extragalactic transients each night. The high volume of alerts demands immediate classification of transient types in order to prioritize observational follow-ups before events fade away. We use host galaxy features to classify transients, thereby providing classification upon discovery. In contrast to past work that focused on distinguishing Type Ia and core-collapse supernovae (SNe) using host galaxy features that are not always accessible (e.g., morphology), we determine the relative likelihood across 12 transient classes based on only 19 host apparent magnitudes and colors from 10 optical and IR photometric bands. We develop both binary and multiclass classifiers, using kernel density estimation to estimate the underlying distribution of host galaxy properties for each transient class. Even in this pilot study, and ignoring relative differences in transient class frequencies, we distinguish eight transient classes at purities significantly above the 8.3% baseline (based on a classifier that assigns labels uniformly and at random): tidal disruption events (TDEs; 48% ± 27%, where ± indicates the 95% confidence limit), SNe Ia-91bg (32% ± 18%), SNe Ia-91T (23% ± 11%), SNe Ib (23% ± 13%), SNe II (17% ± 2%), SNe IIn (17% ± 6%), SNe II P (16% ± 4%), and SNe Ia (10% ± 1%). We demonstrate that our model is applicable to LSST and estimate that our approach can accurately classify 59% of LSST alerts expected each year for SNe Ia, Ia-91bg, II, Ibc, SLSN-I, and TDEs. Our code and data set are publicly available.

A Multiwavelength Study of the Ongoing Changing-look AGN AT2021fxu

We present multiwavelength analyses of an active optical transient AT2021fxu which shows the appearance of previously absent broad emission lines in a recent optical spectrum, suggesting a Changing-look Active Galactic Nucleus (AGN). During the spectral transition, the brightness in the individual photometric bands increased up to ≈0.6 in the optical bands and up to ≈1.1 mag in the UV bands. The brightening was accompanied by a blueward shift of the optical spectrum. AT2021fxu shows high X-ray (0.3–10 keV) flux variability before and after the outburst, with the average X-ray flux increasing by a factor of ≈2 post-outburst. However, the X-ray spectral shape remains roughly the same, with no significant change in the line-of-sight column density. AT2021fxu’s overall properties are consistent with an accretion-rate-driven transition from a Type-II to a Type-1 AGN.

Controlling the artificial radiance of the night sky: The Añora urban laboratory

We provide radiometric evidence of the relevance of nearby light sources on the artificial brightness of the night sky. To obtain the required data we developed a method based on the use of power-regulated urban lighting systems, which also provides relevant information on the propagation of light pollution at short distances from the sources. A controlled experiment was carried out in the Andalusian municipality of Añora (1526 inhab.) in which the radiant power of its regulated public outdoor lighting system was modified in a predefined way during the central part of the night while continuously recording the zenith sky brightness in the TESS-W photometric band from two separate locations inside the town. We determined that the power-regulated streetlight sources contribute a 60–62% to the total zenith sky brightness at these observing locations in clear and moonless nights when operated at their maximum power rating, being the remaining sky radiance due to constant local sources and sources from neighboring towns. These results, in combination with georeferenced information on the sources' location and properties, impose some constraints on the functional form of the effective point-spread functions (PSF) of the zenith sky brightness at short distances from the sources. For the conditions of our experiment we have found that the expected exponents of power-law PSFs are close to -1.

SFR estimations from z = 0 to z = 0.9

Context. The star formation rate (SFR) is a key ingredient for studying the formation and evolution of galaxies. Being able to obtain accurate estimations of the SFR, for a wide range of redshifts, is crucial for building and studying galaxy evolution paths over cosmic time. Aims. Based on a statistical sample of galaxies, the aim of this paper is to constrain a set of SFR calibrators that are able to work in a large redshift range, from z = 0 to z = 0.9. Those calibrators will help to homogenize SFR estimations of star-forming galaxies and to remove any possible biases from the study of galaxy evolution. Methods. Using the VIMOS Public Extragalactic Redshift Survey (VIPERS), we estimated a set of SFR based on photometric and spectroscopic data. We used, as estimators, photometric bands from ultraviolet (UV) to mid-infrared (mid-IR), and the spectral lines Hβ, [O II]λ3727, and [O III]λ5007. Assuming a reference SFR obtained from the spectral energy distribution reconstructed with Code Investigating GALaxy Emission (CIGALE), we estimated the reliability of each band as an SFR tracer. We used the GALEX-SDSS-WISE Legacy Catalog (GSWLC, z < 0.3) to trace the dependence of these SFR calibrators with redshift. Results. The far and near UV (FUV and NUV, respectively), u-band and 24 μm bands, as well as LTIR, are found to be good SFR tracers up to z ∼ 0.9 with a strong dependence on the attenuation prescription used for the bluest bands (scatter of SFR of 0.26, 0.14, 0.15, 0.23, and 0.24 dex for VIPERS, and 0.25, 0.24, 0.09, 0.12, and 0.12 dex for GSWLC). The 8 μm band provides only a rough estimate of the SFR as it depends on metallicity and polycyclic aromatic hydrocarbon properties (scatter of 0.23 dex for VIPERS). We estimated the scatter of rest-frame luminosity estimations from CIGALE to be 0.26, 0.14, 0.12, 0.15, and 0.20 dex for FUV, NUV, ugriz, Ks, and 8–24 μm-LTIR. At intermediate redshift, the Hβ line is a reliable SFR tracer (scatter of 0.19 dex) and the [O II]λ3727 line gives an equally good estimation when the metallicity from the R23 parameter is taken into account (0.17 for VIPERS and 0.20 dex for GSWLC). A calibration based on [O III] retrieves the SFR only when additional information such as the metallicity or the ionization parameter of galaxies are used (0.26 for VIPERS and 0.20 dex for GSWLC), diminishing its usability as a direct SFR tracer. Based on rest-frame luminosities estimated with CIGALE, we propose our own set of calibrations from FUV, NUV, u-band, 8, 24 μm, LTIR, Hβ, [O II], and [O III].

2D modelling of pulsating stars with rapid rotation

Rapid stellar rotation is an important phenomenon in stellar physics, particularly for massive and intermediate mass main-sequence stars. This affects all aspects of the star’s physics including its structure, evolution, and pulsations, and makes it necessary to use 2D numerical approaches rather than the 1D approaches typically used. In this contribution, we will review 2D numerical methods for modelling and interpreting pulsation modes in rapidly rotating stars. We will start by deriving the pulsation equations, both in an adiabatic and non-adiabatic setting, then provide a description of the 2D numerical implementation. We will then explain approximate implementations of the effects of rotation, namely first, second, and third order perturbative approaches, as well as the traditional approximation. This will then be followed by a description on how to calculate disk-integrated mode visibilities in various photometric bands, and how to apply this to mode identification in rapid rotators. Finally, we will review some of the recent works that interpret the pulsation spectra of various stars as viewed in either a single photometric band or in multiple bands, and including supplementary constraints from interferometry and spectroscopy.

RRab variables with identical light-curve shapes at different pulsation periods

ABSTRACT In this paper, we report on the detection of RRab stars with quasi-identical-shape light curves but period differences as large as 0.05–0.21 d using the Galactic bulge data of the OGLE-IV survey. We have examined stars with shorter periods than the Oosterhoff I ridge of the bulge. These stars generally have smaller amplitudes and larger Fourier phase-differences than the typical bulge RRab stars have at the same period. Many of these ‘anomalous’ stars have good-quality light curves without any sign of the Blazhko modulation. Examining their Fourier parameters revealed that several of these stars show very similar light curve to the typical bulge RR Lyrae. We found hundreds of quasi-identical-shape light-curve pairs with different periods between the ‘anomalous’- and the ‘normal’-position RRab stars based on the OGLE I-band data. The OGLE V-band, and the archive VVV and MACHO surveys Ks-, b- and r-band data of these stars were also checked for light-curve-shape similarity. Finally, 149 pairs with identical-shape light curves in each available photometric band were identified. Calculating the physical properties of the variables using empirical formulae, on average, −0.5 dex, −0.13 mag, 0.67, and 165 K differences between the [Fe/H], MV, R/R⊙, and Teff values of the members of the pairs are derived, being the short-period stars less metal-poor, fainter, smaller, and hotter than the long-period variables. To explain the existence of variables with different physical properties and pulsation periods but with identical-shape light curves is a challenging task for modelling.

Superhumps in the cataclysmic variable BG Triangulum

ABSTRACT We present a detailed photometric study of the bright cataclysmic variable, BG Triangulum (BG Tri), using ground-based observations mainly from the Rozhen Observatory, the All-Sky Automated Survey for Supernovae (ASAS-SN), the Transiting Exoplanet Survey Satellite (TESS) and the Wide Angle Search for Planets (WASP). We report on the discovery of a negative superhump with P−sh = 0.1515(2) d and a co-existing superorbital variation with P = 3.94(53) d in data from 2019 and 2020. A positive superhump with P+sh = 0.1727(14) d is also discovered in data from 2006. The obtained negative superhump deficit ε− = 0.044(1) and the positive superhump excess ε+ = 0.090(9) give us an independent photometric evaluation of the mass ratio (q) of the system, which we find to be q− = 0.37(2) and q+ = 0.40(5), respectively. We also present a study of the quasi-periodic oscillations and stochastic variability (flickering) in BG Tri. The light curves show a rich mixture of simultaneously overlapping quasi-periods ranging from 5 to 25 min. The multicolour (UBVRI) photometric observations from the Rozhen Observatory reveal the typical increase of the flickering amplitudes to the shorter wavelengths. The recently introduced A60 amplitude of the flickering light source in all studied photometric bands is systematically lower when the negative superhump is gone in season 2021.

Mira variables in the Milky Way’s nuclear stellar disc: discovery and classification

ABSTRACT The properties of the Milky Way’s nuclear stellar disc give crucial information on the epoch of bar formation. Mira variables are promising bright candidates to study the nuclear stellar disc, and through their period–age relation dissect its star formation history. We report on a sample of 1782 Mira variable candidates across the central $3\times 3\, \mathrm{deg}^2$ of the Galaxy using the multi-epoch infrared VISTA Variables in Via Lactea (VVV) survey. We describe the algorithms employed to select candidate variable stars and then model their light curves using periodogram and Gaussian process methods. By combining with WISE, 2MASS, and other archival photometry, we model the multiband light curves to refine the periods and inspect the amplitude variation between different photometric bands. The infrared brightness of the Mira variables means many are too bright and missed by VVV. However, our sample follows a well-defined selection function as expected from artificial star tests. The multiband photometry is modelled using stellar models with circumstellar dust that characterize the mass-loss rates. We demonstrate how ≳90 per cent of our sample is consistent with O-rich chemistry. Comparison to period–luminosity relations demonstrates that the bulk of the short period stars are situated at the Galactic Centre distance. Many of the longer period variables are very dusty, falling significantly under the O-rich Magellanic Cloud and solar neighbourhood period–luminosity relations and exhibit high mass-loss rates of $\sim 2.5\times 10^{-5}M_\odot \, \mathrm{yr}^{-1}$. The period distribution appears consistent with the nuclear stellar disc forming $\gtrsim 8\, \mathrm{Gyr}$ ago, although it is not possible to disentangle the relative contributions of the nuclear stellar disc and the contaminating bulge.

Photometric Metallicity Prediction of Fundamental-mode RR Lyrae Stars in the Gaia Optical and K s Infrared Wave Bands by Deep Learning

RR Lyrae stars are useful chemical tracers thanks to the empirical relationship between their heavy-element abundance and the shape of their light curves. However, the consistent and accurate calibration of this relation across multiple photometric wave bands has been lacking. We have devised a new method for the metallicity estimation of fundamental-mode RR Lyrae stars in the Gaia optical G and near-infrared VISTA K s wave bands by deep learning. First, an existing metallicity prediction method is applied to large photometric data sets, which are then used to train long short-term memory recurrent neural networks for the regression of the [Fe/H] to the light curves in other wave bands. This approach allows an unbiased transfer of our accurate, spectroscopically calibrated I-band formula to additional bands at the expense of minimal additional noise. We achieve a low mean absolute error of 0.1 dex and high R 2 regression performance of 0.84 and 0.93 for the K s and G bands, respectively, measured by cross-validation. The resulting predictive models are deployed on the Gaia DR2 and VVV inner bulge RR Lyrae catalogs. We estimate mean metallicities of −1.35 dex for the inner bulge and −1.7 dex for the halo, which are significantly less than the values obtained by earlier photometric prediction methods. Using our results, we establish a public catalog of photometric metallicities of over 60,000 Galactic RR Lyrae stars and provide an all-sky map of the resulting RR Lyrae metallicity distribution. The software code used for training and deploying our recurrent neural networks is made publicly available in the open-source domain.

New calibrated models for the tip of the red giant branch luminosity and a thorough analysis of theoretical uncertainties

ABSTRACT The luminosity of the tip of the red giant branch (TRGB) is instrumental for the construction of the distance ladder, and its accurate modelling is key for determining the local Hubble parameter. In this work, we present an extensive quantitative analysis of the TRGB luminosity, accounting for virtually all input physics that affect it, namely: chemical composition, opacity, diffusion, nuclear reaction rates, electron screening, neutrinos, convection efficiency, boundary conditions, and mass-loss, amongst others. Our analysis is based on a newly produced grid of ∼3 × 106 TRGB models, evolved from pre-main sequence up to the helium ignition at the TRGB, and covering a wide range of metallicity (Z = 0.0001–0.02) and initial mass (M = 0.8–1.4 M⊙). Through a Monte Carlo analysis, we study the systematic variation of the TRGB luminosity due to the combined effect of all above input physics, and show that a maximum theoretical uncertainty of about $1.6 {{\ \rm per\ cent}}$ is still present on the current generation of models, dominated by systematics of radiative opacity. Results are also provided in several photometric bands. As a by-product of our analysis, we demonstrate robust evidence for the linear response of the tip luminosity to individual changes of input physics, which can significantly simplify future analyses. A comparison of our results with other stellar evolution codes shows excellent agreement. Our full grid of models is available upon reasonable requests, while a partial grid is available online on Zenodo.