The Overlooked Potential of Generalized Linear Models in Astronomy-III: Bayesian Negative Binomial Regression and Globular Cluster Populations

In this paper, the third in a series illustrating the power of generalized linear models (GLMs) for the astronomical community, we elucidate the potential of the class of GLMs which handles count data. The size of a galaxy's globular cluster population $N_{\rm GC}$ is a prolonged puzzle in the astronomical literature. It falls in the category of count data analysis, yet it is usually modelled as if it were a continuous response variable. We have developed a Bayesian negative binomial regression model to study the connection between $N_{\rm GC}$ and the following galaxy properties: central black hole mass, dynamical bulge mass, bulge velocity dispersion, and absolute visual magnitude. The methodology introduced herein naturally accounts for heteroscedasticity, intrinsic scatter, errors in measurements in both axes (either discrete or continuous), and allows modelling the population of globular clusters on their natural scale as a non-negative integer variable. Prediction intervals of 99% around the trend for expected $N_{\rm GC}$comfortably envelope the data, notably including the Milky Way, which has hitherto been considered a problematic outlier. Finally, we demonstrate how random intercept models can incorporate information of each particular galaxy morphological type. Bayesian variable selection methodology allows for automatically identifying galaxy types with different productions of GCs, suggesting that on average S0 galaxies have a GC population 35% smaller than other types with similar brightness.Comment: 14 pages, 12 figures. Accepted for publication in MNRA

J-PLUS: Identification of low-metallicity stars with artificial neural networks using SPHINX

We present a new methodology for the estimation of stellar atmospheric parameters from narrow- and intermediate-band photometry of the Javalambre Photometric Local Universe Survey (J-PLUS), and propose a method for target pre-selection of low-metallicity stars for follow-up spectroscopic studies. Photometric metallicity estimates for stars in the globular cluster M15 are determined using this method. By development of a neural-network-based photometry pipeline, we aim to produce estimates of effective temperature, $T_{\rm eff}$, and metallicity, [Fe/H], for a large subset of stars in the J-PLUS footprint. The Stellar Photometric Index Network Explorer, SPHINX, is developed to produce estimates of $T_{\rm eff}$ and [Fe/H], after training on a combination of J-PLUS photometric inputs and synthetic magnitudes computed for medium-resolution (R ~ 2000) spectra of the Sloan Digital Sky Survey. This methodology is applied to J-PLUS photometry of the globular cluster M15. Effective temperature estimates made with J-PLUS Early Data Release photometry exhibit low scatter, \sigma($T_{\rm eff}$) = 91 K, over the temperature range 4500 < $T_{\rm eff}$ (K) < 8500. For stars from the J-PLUS First Data Release with 4500 < $T_{\rm eff}$ (K) < 6200, 85 $\pm$ 3% of stars known to have [Fe/H] <-2.0 are recovered by SPHINX. A mean metallicity of [Fe/H]=-2.32 $\pm$ 0.01, with a residual spread of 0.3 dex, is determined for M15 using J-PLUS photometry of 664 likely cluster members. We confirm the performance of SPHINX within the ranges specified, and verify its utility as a stand-alone tool for photometric estimation of effective temperature and metallicity, and for pre-selection of metal-poor spectroscopic targets.Comment: 18 pages, 12 figure

J-PLUS: analysis of the intracluster light in the Coma cluster

The intracluster light (ICL) is a luminous component of galaxy clusters composed of stars that are gravitationally bound to the cluster potential but do not belong to the individual galaxies. Previous studies of the ICL have shown that its formation and evolution are intimately linked to the evolutionary stage of the cluster. Thus, the analysis of the ICL in the Coma cluster will give insights into the main processes driving the dynamics in this highly complex system. Using a recently developed technique, we measure the ICL fraction in Coma at several wavelengths, using the J-PLUS unique filter system. The combination of narrow- and broadband filters provides valuable information on the dynamical state of the cluster, the ICL stellar types, and the morphology of the diffuse light. We use the Chebyshev-Fourier Intracluster Light Estimator (CICLE) to disentangle the ICL from the light of the galaxies, and to robustly measure the ICL fraction in seven J-PLUS filters. We obtain the ICL fraction distribution of the Coma cluster at different optical wavelengths, which varies from $\sim 7\%-21\%$, showing the highest values in the narrowband filters J0395, J0410, and J0430. This ICL fraction excess is distinctive pattern recently observed in dynamically active clusters (mergers), indicating a higher amount of bluer stars in the ICL compared to the cluster galaxies. Both the high ICL fractions and the excess in the bluer filters are indicative of a merging state. The presence of younger/lower-metallicity stars the ICL suggests that the main mechanism of ICL formation for the Coma cluster is the stripping of the stars in the outskirts of infalling galaxies and, possibly, the disruption of dwarf galaxies during past/ongoing mergers.Comment: 10 pages, 3 figures, 1 table. Accepted for publication in A&

Rules, Standards, and the Internal Point of View

Large scale structure and cosmolog

J-PLUS: Detecting and studying extragalactic globular clusters -- the case of NGC 1023

Extragalactic globular clusters (GCs) are key objects for studying the formation and evolution of galaxies. The arrival of wide-field surveys such as the Javalambre Photometric Local Universe Survey (J-PLUS) offers new possibilities for the study of GCs. Nevertheless, GCs are not detected a priori by the data reduction pipeline of J-PLUS and, due to its pixel scale, the standard techniques of GCs detection are challenged. To fill this gap, we develop a semi-automatic pipeline to detect GCs in J-PLUS that can also be adapted to similar surveys. As a case study, we use data from the S0 galaxy NGC 1023 and we also study the stellar population content of GC candidates in the galaxy. To detect GCs, our methodology is based on Source Extractor and does not require a previous filtering or modelling of the host galaxy. We study colors and perform spectral energy distribution (SED) analysis on our final GC candidate catalog to obtain stellar population parameters. In NGC 1023, GCFinder identifies 523 GC candidates. We observe evidence of color bimodality in a few broad-band colors but not on narrow-band colors. The SED analysis reveals a clear metallicity bimodality and we observe that narrow-band filters are very useful to constrain metallicities. We also identified a broad age-metallicity relation as well as a wide metallicity distribution that are evidence that NGC 1023 experienced accretion events in the past. It is the first time this kind of study is performed with J-PLUS data. By detecting GC candidates in wide-field images without modeling the light of the galaxy, GCFinder becomes considerably faster, at a marginal loss of centrally-located GC candidates of about 7 percent. As GCFinder is entirely based on Source Extractor, it could be easily incorporated into automated software handling wide-field surveys.Comment: 21 pages, 19 figures, submitted to A&

The Globular Cluster System of NGC 5846 Revisited: Colours, Sizes and X-Ray Counterparts

NGC 5846 is a giant elliptical galaxy with a previously well studied globular cluster system (GCS), known to have a bimodal colour distribution with a remarkably high red fraction. Here we revisit the central galaxy regions searching for new globular cluster (GC) candidates, and measuring, magnitudes, colours and sizes for them. We also search for their X-ray counterparts. We use archival Hubble Space Telescope WFPC2 images, from which we modelled and subtracted the host light distribution and increased the available sample of GCs. We performed photometry on the central objects, and measured sizes and equatorial coordinates for the entire system known in this galaxy. We detect two dozen previously unknown GC candidates in the central regions. Reliable sizes are obtained for about 60 GCs; their typical effective radii are in the range 3-5 pc. The largest clusters are located in the central regions. We find 7 X-ray counterparts to globular clusters, most of them in the central region. They are among the most luminous X-ray sources in NGC 5846. They are also optically luminous, compact and belong to the red subpopulation.Comment: 8 pages, 9 figures, A&A accepte

Molecular Genetic Variability, Population Structure and Mating System in Tropical Forages

Microsatellite (SSR) markers were developed for the following tropical forage species, using accessions available from the plant genetic resources (PGR) collections held by EMBRAPA (Brazilian Agricultural Research Corporation): Brachiaria brizantha, B. humidicola, Panicum maximum, Paspalum spp., Stylosanthes capitata, S. guianensis, S. macrocephala, Calopogonium mucunoides and Centrosema spp. The markers were used to analyse population structure and genetic diversity, evolution and origin of the genetic variability in the centre of origin, mating systems and genetic resources in EMBRAPA’s germplasm bank. The results shed light on the amount of genetic variation within and between populations, revealed the need in some cases for further plant collection to adequately represent the species in PGR collections, allowed us to assemble core collections (subsets of the total collections) that should contain most of the available diversity and (in the case of the legumes) showed the need to avoid unwanted outcrossing when regenerating conserved material. The data will allow plant breeders to better select accessions for hybrid production, discriminate between genotypes and use marker-assisted selection in breeding programs. Our results will also underpin the construction of genetic maps, mapping of genes of agronomic interest and numerous other studies on genetic variability, population structure, gene flow and reproductive systems for the tropical forage species studied in this work