Search

We present capivara, a fast and scalable multi-decomposition package designed to study astrophysical properties within distinct structural components of galaxies. Our spectro-decomposition code for analyzing integral field unit (IFU) data enables a more holistic approach, moving beyond conventional radial gradients and the bulge-plus-disk dichotomy. It facilitates comprehensive comparisons of integrated stellar ages and metallicities across various galactic structures. Our classification method naturally identifies outliers and organizes the different pixels based on their dominant spectral features. The algorithm leverages the scalability and GPU acceleration of Torch, outputting both a one-dimensional spectrum and a full data cube for each galaxy component, without relying on Voronoi binning. We demonstrate the capabilities of our approach using a sample of galaxies from the MaNGA survey, processing the resulting data cubes with the starlight spectral fitting code to derive both stellar population and ionized gas properties of the galaxy components. Our method effectively groups regions with similar spectral properties in both the continuum and emission lines. By aggregating the spectra of these regions, we enhance the signal-to-noise ratio of our analysis while significantly speeding up computations by reducing the number of spectra processed simultaneously. capivara will be freely available on GitHub., Comment: Submitted to MNRAS. Comments and suggestions are very welcome

Understanding faint dwarf galaxies is fundamental to the development of a robust theory of galaxy formation on small scales. Since the discovery of a population of ultra diffuse galaxies (UDGs) rich in globular clusters (GCs) in Coma, an increasing number of studies on low surface brightness dwarf galaxies (LSBds) have been published in recent years. The most massive LSBds have been observed predominantly in groups and clusters, with properties displaying dependence on the environment. In this work, we use deep DECam imaging to systematically identify LSBds and their GC populations around the low-density environment of NGC 3115. We carefully analyse the structure and morphology of 24 candidates, 18 of which are reported for the first time. Most candidates exhibit red colours suggesting a connection between their colour and distance to NGC 3115. We followed up with Gemini GMOS imaging 9 LSBds to properly identify their GC populations. We derive lower limits for the number of GCs associated with each galaxy. Our analysis reveals that they occur around of the same loci of Fornax LSB dwarf GC systems. The relationship between the number of GCs and total mass provides a tool in which, by counting the GCs in these galaxies, we estimate an upper limit for the total mass of these LSB dwarfs, obtaining the mean value of $\sim 3.3\times10^{10}$ M$_{\odot}$. Our results align with expectations for dwarf-sized galaxies, particularly regarding the distribution and specific frequency of their GC systems., Comment: 25 pages, 22 figures, accepted for publication in MNRAS

Strong gravitational lensing and stellar dynamics are independent and powerful methods to probe the total gravitational potential of galaxies, and thus, their total mass profile. However, inherent degeneracies in the individual models makes it difficult to obtain a full understanding of the distribution of baryons and dark matter (DM), although such degeneracies might be broken by the combination of these two tracers, leading to more reliable measurements of the mass distribution of the lens galaxy. We use mock data from IllustrisTNG50 to compare how dynamical-only, lens-only, and joint modelling can constrain the mass distribution of early-type galaxies (ETGs). The joint model consistently outperforms the other models, achivieng a $2\%$ accuracy in recovering the total mass within $2.5R_\text{eff}$. The Einstein radius is robustly recovered for both lens-only and joint models, with the first showing a median fractional error of $-5\%$ and the latter a fractional error consistent with zero. The stellar mass-to-light ratio and total mass density slope are well recovered by all models. In particular, the dynamical-only model achieves an accuracy of $1\%$ for the stellar mass-to-light ratio, while the accuracy of the mass density slope is typically of the order of $5\%$ for all models. However, all models struggle to constrain integrated quantities involving DM and the halo parameters. Nevertheless, imposing more restrictive assumptions on the DM halo, such as fixing the scale radius, could alleviate some of the issues. Finally, we verify that the number of kinematical constraints ($15, 35, 55$ bins) on the kinematical map does not impact the models outcomes., Comment: 42 pages, 13 Figures, and 2 Tables

We study 51 jellyfish galaxy candidates in the Fornax, Antlia, and Hydra clusters. These candidates are identified using the JClass scheme based on the visual classification of wide-field, twelve-band optical images obtained from the Southern Photometric Local Universe Survey. A comprehensive astrophysical analysis of the jellyfish (JClass > 0), non-jellyfish (JClass = 0), and independently organized control samples is undertaken. We develop a semi-automated pipeline using self-supervised learning and similarity search to detect jellyfish galaxies. The proposed framework is designed to assist visual classifiers by providing more reliable JClasses for galaxies. We find that jellyfish candidates exhibit a lower Gini coefficient, higher entropy, and a lower 2D S\'ersic index as the jellyfish features in these galaxies become more pronounced. Jellyfish candidates show elevated star formation rates (including contributions from the main body and tails) by $\sim$1.75 dex, suggesting a significant increase in the SFR caused by the ram-pressure stripping phenomenon. Galaxies in the Antlia and Fornax clusters preferentially fall towards the cluster's centre, whereas only a mild preference is observed for Hydra galaxies. Our self-supervised pipeline, applied in visually challenging cases, offers two main advantages: it reduces human visual biases and scales effectively for large datasets. This versatile framework promises substantial enhancements in morphology studies for future galaxy image surveys., Comment: 25 pages, 22 figures, accepted for publication in MNRAS

Most of the observed galaxies cannot be resolved into individual stars and are studied through their integrated spectrum using simple stellar populations (SSPs) models, with stellar libraries being a key ingredient in building them. Spectroscopic observations are increasingly being directed towards the near-infrared (NIR), where much is yet to be explored. SSPs in the NIR are still limited, and there are inconsistencies between different sets of models. One of the ways to minimize this problem is to have reliable NIR stellar libraries. The main goal of this work is to present SMARTY (mileS Moderate resolution neAr-infRared sTellar librarY) a ~0.9-2.4$\mu$m stellar spectral library composed of 31 stars observed with the Gemini Near-IR Spectrograph (GNIRS) at the 8.1m Gemini North telescope and make it available to the community. The stars were chosen from the SMARTY library, for which the atmospheric parameters are reliable (and well tested), to populate different regions of the Hertzsprung-Russell (HR) diagram. Furthermore, five of these stars have NIR spectra available that we use to assess the quality of SMARTY. The remaining 26 stars are presented for the first time in the NIR. We compared the observed SMARTY spectra with synthetic and interpolated spectra, finding a mean difference of ~20% in the equivalent widths and ~1% in the overall continuum shape in both sets of comparisons. We computed the spectrophotometric broadband magnitudes and colours and compared them with the 2MASS ones, resulting in mean differences up to 0.07 and 0.10mag in magnitudes and colours, respectively. In general, a small difference was noted between the SMARTY spectra corrected using the continuum from the interpolated and the theoretical stars.

We introduce galmoss, a python-based, torch-powered tool for two-dimensional fitting of galaxy profiles. By seamlessly enabling GPU parallelization, galmoss meets the high computational demands of large-scale galaxy surveys, placing galaxy profile fitting in the LSST-era. It incorporates widely used profiles such as the S\'ersic, Exponential disk, Ferrer, King, Gaussian, and Moffat profiles, and allows for the easy integration of more complex models. Tested on 8,289 galaxies from the Sloan Digital Sky Survey (SDSS) g-band with a single NVIDIA A100 GPU, galmoss completed classical S\'ersic profile fitting in about 10 minutes. Benchmark tests show that galmoss achieves computational speeds that are 6 $\times$ faster than those of default implementations., Comment: 12 pages, 8 figures, Accepted for publication in A&C

Nuclear star clusters (NSCs) are dense star clusters located at the centre of galaxies spanning a wide range of masses and morphologies. Analysing NSC occupation statistics in different environments provides an invaluable window into investigating early conditions of high-density star formation and mass assembly in clusters and group galaxies. We use HST/ACS deep imaging to obtain a catalogue of dwarf galaxies in two galaxy clusters in the Shapley Supercluster: the central cluster Abell 3558 and the northern Abell 1736a. The Shapley region is an ideal laboratory to study nucleation as it stands as the highest mass concentration in the nearby Universe. We investigate the NSC occurrence in quiescent dwarf galaxies as faint as $M_{I} = -10$ mag and compare it with all other environments where nucleation data is available. We use galaxy cluster/group halo mass as a proxy for the environment and employ a Bayesian logistic regression framework to model the nucleation fraction ($f_{n}$) as a function of galaxy luminosity and environment. We find a notably high $f_n$ in Abell 3558: at $M_{I} \approx -13.1$ mag, half the galaxies in the cluster host NSCs. This is higher than in the Virgo and Fornax clusters but comparable to the Coma Cluster. On the other hand, the $f_n$ in Abell 1736a is relatively lower, comparable to groups in the Local Volume. We find that the probability of nucleation varies with galaxy luminosity remarkably similarly in galaxy clusters. These results reinforce previous findings of the important role of the environment in NSC formation/growth., Comment: 20 pages, 10 figures, accepted for publication in the MNRAS

Relic galaxies are massive, compact, quiescent objects observed in the local Universe that have not experienced any significant interaction episodes or merger events since about $z = 2$, remaining relatively unaltered since their formation. On the other hand, massive and compact Early Type Galaxies (cETGs) in the local Universe appear to show similar properties to Relic galaxies, despite having substantial accretion history. Relic galaxies, with frozen history, can provide important clues about the intrinsic processes related to the evolutionary pathways of ETGs and the role that mergers play in their evolution. Using the high-resolution cosmological simulation TNG50-1 from the Illustris Project, we investigate the assembly history of a sample of massive, compact, old, and quiescent subhalos split by satellite accretion fraction. We compare the evolutionary pathways at three cosmic epochs: $z = 2$, $z = 1.5$, and $z = 0$, using the orbital decomposition numerical method to investigate the stellar dynamics of each galactic kinematical component and their environmental correlations. Our results point to a steady pathway across time that is not strongly dependent on the environment. Relics and cETGs do not show a clear preference for high or low-density environments within the volume explored at TNG50. However, progenitors of Relic galaxies are shown to be located in high density since $z = 2$. The merger history can be recovered from the hot inner stellar halo imprints in the local Universe. In the current scenario, the mergers that drive the growth of cETGs do not give rise to a new and distinct evolutionary pathway when compared to Relics. This is despite the reported effects on the age and metallicity of the kinematic components., Comment: 12 pages, 13 figures; accepted for publication in MNRAS

The Vera C. Rubin Observatory will undertake the Legacy Survey of Space and Time, providing an unprecedented, volume-limited catalog of star clusters in the Southern Sky, including Galactic and extragalactic star clusters. The Star Clusters subgroup of the Stars, Milky Way and Local Volume Working Group has identified key areas where Rubin Observatory will enable significant progress in star cluster research. This roadmap represents our science cases and preparation for studies of all kinds of star clusters from the Milky Way out to distances of tens of megaparsecs., Comment: accepted to PASP

We present a spatially resolved stellar population analysis of 61 jellyfish galaxies and 47 control galaxies observed with ESO/MUSE attempting to understand the general trends of the stellar populations as a function of the stripping intensity and mass. This is the public sample from the GASP programme, with $0.01 < z < 0.15$ and $8.9 <\log(M_{\star}/M_{\odot}) < 12.0$. We apply the spectral population synthesis code FADO to fit self-consistently both the stellar and nebular contributions to the spectra of the sources. We present 2D morphological maps for mean stellar ages, metallicities, gas-phase oxygen abundances, and star formation rates for the galaxies with Integrated Nested Laplace Approximation (INLA), which is efficient in reconstructing spatial data of extended sources. We find that "extreme stripping" and "stripping" galaxies are typically younger than the other types. Regarding stellar and nebular metallicities, the "stripping" and "control passive" galaxies are the most metal-poor. Based on the phase space for jellyfish cluster members we find trends in ages, metallicities, and abundances with different regions of the diagram. We also compute radial profiles for the same quantities. We find that both the stripping and the stellar masses seem to influence the profiles, and we see differences between various groups and distinct mass bins. The radial profiles for different mass bins present relations already shown in the literature for undisturbed galaxies, i.e., profiles of ages and metallicities tend to increase with mass. However, beyond $\sim0.75$ effective radius, the ages of the most massive galaxies become similar to or lower than the ages of the lower mass ones., Comment: Accepted by MNRAS. 71 pages (13 of main text and 58 of supplementay material). 123 figures

General Relativity (GR) has been successfully tested mainly at Solar system scales; however, galaxy-scale tests have become popular in the last few decades. In this work, we investigate the $\eta_\text{PPN}$ parameter, which is commonly defined by the ratio of two scalar potentials that appears in the cosmological linearly perturbed metric. Under the assumption of GR and a vanish anisotropic stress tensor, $\eta_\text{PPN}= 1$. Using ALMA, HST, and VLT/MUSE data, we combine mass measurements, using gravitational lensing and galactic dynamics, for the SDP.81 lens galaxy ($z = 0.299$) to constrain $\eta_\text{PPN}$. By using a flexible and self-consistent mass profile, our fiducial model takes into account the contribution of the stellar mass and a dark matter halo to reconstruct the lensed galaxy and the spatially-resolved stellar kinematics. We infer, after accounting for systematic uncertainties related to the mass model, cosmology and kinematics, $\eta_{\text{PPN}} = 1.13^{+0.03}_{-0.03}\pm0.20\,(\text{sys})$, which is in accordance with GR predictions. Better spectroscopy data are needed to push the systematics down and bring the uncertainty to the percentage level since our analysis shows that the main source of the systematics is related to kinematics, which heavily depends on the signal-to-noise ratio of the spectra., Comment: 18 pages, 8 figures. Accepted to MNRAS

The Ly$\alpha$ emission is an important tracer of neutral gas in a circum-galactic medium (CGM) around high-z QSOs. The origin of Lya emission around QSOs is still under debate which has significant implications for galaxy formation and evolution. In this paper, we study Ly$\alpha$ nebulae around two high redshift QSOs, SDSS J141935.58+525710.7 at $z=3.218$ (hereafter QSO1) and SDSS J141813.40+525240.4 at $z=3.287$ (hereafter QSO2), from the miniJPAS survey within the AEGIS field. Using the contiguous narrow-band (NB) images from the miniJPAS survey and SDSS spectra, we analyzed their morphology, nature, and origin. We report the serendipitous detection of double-core Ly\al\ morphology around two QSOs which is rarely seen among other QSOs. The separations of the two Ly\al~cores are 11.07 $\pm$ 2.26 kpcs (1.47 $\pm$ 0.3$^{\prime\prime}$) and 9.73 $\pm$ 1.55 kpcs (1.31 $\pm$ 0.21$^{\prime\prime}$) with Ly$\alpha$~line luminosities of $\sim$ 3.35 $\times 10^{44}$ erg s $^{-1} $ and $\sim$ 6.99 $\times$ 10$^{44}$ erg s $^{-1}$ for QSO1 and QSO2, respectively. The miniJPAS NB images show evidence of extended Ly$\alpha$ and CIV morphology for both QSOs and extended HeII morphology for QSO1. These two QSOs may be potential candidates for the new enormous Lyman alpha nebula (ELAN) found from the miniJPAS survey due to their extended morphology in the shallow depth and relatively high Ly$\alpha$ luminosities. We suggest that galactic outflows are the major powering mechanism for the double-core Ly$\alpha$ morphology. Considering the relatively shallow exposures of miniJPAS, the objects found here could be the tip of the iceberg of a promising number of such objects that will be uncovered in the upcoming full J-PAS survey and deep IFU observations with 8-10m telescopes will be essential for constraining the underlying physical mechanism that is responsible for the double-cored morphology., Comment: Five figures and two tables (Accepted in Astronomy and Astrophysics)

Galaxy morphological classification is a fundamental aspect of galaxy formation and evolution studies. Various machine learning tools have been developed for automated pipeline analysis of large-scale surveys, enabling a fast search for objects of interest. However, crowded regions in the image may pose a challenge as they can lead to bias in the learning algorithm. In this Research Note, we present galmask, an open-source package for unsupervised galaxy masking to isolate the central object of interest in the image. galmask is written in Python and can be installed from PyPI via the pip command., Comment: Submitted to RNAAS

Globular clusters (GCs) are proxies of the formation assemblies of their host galaxies. However, few studies exist targeting GC systems of spiral galaxies up to several effective radii. Through 12-band Javalambre Photometric Local Universe Survey (J-PLUS) imaging, we study the point sources around the M81/M82/NGC3077 triplet in search of new GC candidates. We develop a tailored classification scheme to search for GC candidates based on their similarity to known GCs via a principal components analysis (PCA) projection. Our method accounts for missing data and photometric errors. We report 642 new GC candidates in a region of 3.5 deg$^2$ around the triplet, ranked according to their Gaia astrometric proper motions when available. We find tantalising evidence for an overdensity of GC candidate sources forming a bridge connecting M81 and M82. Finally, the spatial distribution of the GC candidates $(g-i)$ colours is consistent with halo/intra-cluster GCs, i.e. it gets bluer as they get further from the closest galaxy in the field. We further employ a regression-tree based model to estimate the metallicity distribution of the GC candidates based on their J-PLUS bands. The metallicity distribution of the sample candidates is broad and displays a bump towards the metal-rich end. Our list increases the population of GC candidates around the triplet by 3-fold, stresses the usefulness of multi-band surveys in finding these objects, and provides a testbed for further studies analysing their spatial distribution around nearby (spirals) galaxies., Comment: 20 pages, 18 figures, MNRAS accepted

The presence of Multiple Stellar Populations (MSPs) in Galactic Globular Clusters (GCs) is a poorly understood phenomenon. By probing different spectral ranges that are affected by different absorption lines using the multi-band photometric survey S-PLUS, we study four GCs -- NGC 104, NGC 288, NGC 3201 and NGC 7089 -- that span a wide range in metallicities. With the combination of broad and narrow-band photometry in 12 different filters from 3485A (u) to 9114A (z), we identified MSPs along the rectified red-giant branch in colour-magnitude diagrams (CMDs) and separated them using a K-means clustering algorithm. Additionally, we take advantage of the large Field of View of the S-PLUS detector to investigate radial trends in our sample. We report on six colour combinations that can be used to successfully identify two stellar populations in all studied clusters and show that they can be characterized as Na-rich and Na-poor. For both NGC 288 and NGC 7089, their radial profiles show a clear concentration of 2P. This directly supports the formation theories that propose an enrichment of the intra-cluster medium and subsequent star formation in the more dense central regions. However, in the case of NGC 3201, the trend is reversed. The 1P is more centrally concentrated, in direct contradiction with previous literature studies. NGC 104 shows a well-mixed population. We also constructed radial profiles up to 1 half-light radius of the clusters with HST data to highlight that radial differences are lost in the inner regions of the GCs and that wide-field studies are essential when studying this., Comment: 9 pages, 11 figures

In this series of papers, we employ several machine learning (ML) methods to classify the point-like sources from the miniJPAS catalogue, and identify quasar candidates. Since no representative sample of spectroscopically confirmed sources exists at present to train these ML algorithms, we rely on mock catalogues. In this first paper we develop a pipeline to compute synthetic photometry of quasars, galaxies and stars using spectra of objects targeted as quasars in the Sloan Digital Sky Survey. To match the same depths and signal-to-noise ratio distributions in all bands expected for miniJPAS point sources in the range $17.5\leq r<24$, we augment our sample of available spectra by shifting the original $r$-band magnitude distributions towards the faint end, ensure that the relative incidence rates of the different objects are distributed according to their respective luminosity functions, and perform a thorough modeling of the noise distribution in each filter, by sampling the flux variance either from Gaussian realizations with given widths, or from combinations of Gaussian functions. Finally, we also add in the mocks the patterns of non-detections which are present in all real observations. Although the mock catalogues presented in this work are a first step towards simulated data sets that match the properties of the miniJPAS observations, these mocks can be adapted to serve the purposes of other photometric surveys., Comment: 20 pages, 18 figures, submitted to MNRAS

NGC 4382 is a merger-remnant galaxy that has been classified as morphological type E2, S0, and even Sa. In this work, we performed a photometric and spectroscopic analysis of the globular cluster (GC) system of this peculiar galaxy in order to provide additional information about its history. We used a combination of photometric data in different filters, and multi-object and long-slit spectroscopic data obtained using the Gemini/GMOS instrument. The photometric analysis of the GC system, using the Gaussian Mixture Model algorithm in the colour plane, reveals a complex colour distribution within $R_\mathrm{gal}<5$ arcmin (26.1 kpc), showing four different groups: the typical blue and red subpopulations, a group with intermediate colours, and the fourth group towards even redder colours. From the spectroscopic analysis of 47 GCs, confirmed members of NGC\,4382 based on radial velocities, we verified 3 of the 4 photometric groups from the analysis of their stellar populations using the ULySS code. NGC 4382 presents the classic blue ($10.4\pm2.8$ Gyr, $\mathrm{[Fe/H]}=-1.48\pm0.18$ dex) and red ($12.1\pm2.3$ Gyr, $\mathrm{[Fe/H]}=-0.64\pm0.26$ dex) GCs formed earlier in the lifetime of the galaxy, and a third group of young GCs ($2.2\pm0.9$ Gyr; $\mathrm{[Fe/H]}=-0.05\pm0.28$ dex). Finally, analysis of long-slit data of the galaxy reveals a luminosity-weighted mean age for the stellar population of $\sim$2.7 Gyr, and an increasing metallicity from [Fe/H]=$-0.1$ to $+0.2$ dex in $R_\mathrm{gal}<10$ arcsec (0.87 kpc). These values, and other morphological signatures in the galaxy, are in good agreement with the younger group of GCs, indicating a common origin as a result of a recent merger., Comment: 20 pages, 24 figures. Accepted to be published in MNRAS

Relic galaxies are massive compact quiescent galaxies that formed at high-redshift and remained almost unchanged since then. In this work, we search for analogues to relic galaxies in the TNG50 cosmological simulations to understand relic formation and test the ability of TNG50 to reproduce such rare objects. Using stellar mass, age, radius, quiescence and stellar assembly criteria, we find 5 subhalos in TNG50 that could be potential relic analogues. We compare their properties with other constraints imposed by a sample of 13 observed relic galaxies. We find one analogue in TNG50 that simultaneously satisfies most of the available observational constraints, such as metallicity and morphology. It also shows similarities to the confirmed relic NGC 1277, regarding environment and dark matter fraction. By taking into account a degree of relicness, a second relic analogue may also be considered. However, the central parts of the brightness and density profiles of the analogues are less steep than that of real relic galaxies, possibly due to limited numerical resolution. We identify two formation pathways of relic analogues in TNG50 depending on their environment: they either have their remaining gas stripped during the infall into a cluster at $z \lesssim 1.2$ or consume it before $z > 1.5$. They are then deprived of significant star formation, leaving their stellar populations almost unaltered during the last 9 Gyr. We also find that the analogue progenitors at $z \sim 4$ inhabit more massive halos than progenitors of quiescent galaxies with similar stellar mass at $z \sim 0$., Comment: 20 pages, 12 figures. Submitted to MNRAS, version after 1st referee report

In this paper, we describe a new method to select globular cluster (GC) candidates, including galaxy subtraction with unsharp masking, template fitting techniques and the inclusion of Gaia's proper motions. We report the use of the 12-band photometric system of S-PLUS to determine radial velocities and stellar populations of GCs around nearby galaxies. Specifically, we assess the effectiveness of identifying GCs around nearby and massive galaxies (D $< 20$ Mpc and $\sigma > 200$ km/s) in a multi-band survey such as S-PLUS by using spectroscopically confirmed GCs and literature GC candidate lists around the bright central galaxy in the Fornax cluster, NGC 1399 (D = 19 Mpc), and the isolated lenticular galaxy NGC 3115 (D = 9.4 Mpc). Despite the shallow survey depth, that limits the present work to $r < 21.3$ mag, we measure reliable photometry and perform robust SED fitting for a sample of 115 GCs around NGC 1399 and 42 GCs around NGC 3115, recovering radial velocities, ages, and metallicities for the GC populations., Comment: Accepted for publication in MNRAS. 11 pages, 5 figures, 2 tables