BAYES-LOSVD: a bayesian framework for non-parametric extraction of the line-of-sight velocity distribution of galaxies

We introduce BAYES-LOSVD, a novel implementation of the non-parametric extraction of line-of-sight velocity distributions (LOSVDs) in galaxies. We employ bayesian inference to obtain robust LOSVDs and associated uncertainties. Our method relies on principal component analysis to reduce the dimensionality of the base of templates required for the extraction and thus increase the performance of the code. In addition, we implement several options to regularise the output solutions. Our tests, conducted on mock spectra, confirm the ability of our approach to model a wide range of LOSVD shapes, overcoming limitations of the most widely used parametric methods (e.g. Gauss-Hermite expansion). We present examples of LOSVD extractions for real galaxies with known peculiar LOSVD shapes, i.e. NGC4371, IC0719 and NGC4550, using MUSE and SAURON integral-field unit (IFU) data. Our implementation can also handle data from other popular IFU surveys (e.g. ATLAS3D, CALIFA, MaNGA, SAMI). Details of the code and relevant documentation are freely available to the community in the dedicated repositories.Comment: 13 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics. Public repository with the code can be found at: https://github.com/jfalconbarroso/BAYES-LOSV

The loss of the intracluster medium in globular clusters

Stars in globular clusters (GCs) lose a non negligible amount of mass during their post-main sequence evolution. This material is then expected to build up a substantial intra-cluster medium (ICM) within the GC. However, the observed gas content in GCs is a couple of orders of magnitude below these expectations. Here we follow the evolution of this stellar wind material through hydrodynamical simulations to attempt to reconcile theoretical predictions with observations. We test different mechanisms proposed in the literature to clear out the gas such as ram-pressure stripping by the motion of the GC in the Galactic halo medium and ionisation by UV sources. We use the code ramses to run 3D hydrodynamical simulations to study for the first time the ICM evolution within discretised multi-mass GC models including stellar winds and full radiative transfer. We find that the inclusion of both ram-pressure and ionisation is mandatory to explain why only a very low amount of ionised gas is observed in the core of GCs. The same mechanisms operating in ancient GCs that clear the gas could also be efficient at younger ages, meaning that young GCs would not be able to retain gas and form multiple generations of stars as assumed in many models to explain "multiple populations". However, this rapid clearing of gas is consistent with observations of young massive clusters

How to bend galaxy disc profiles - II. Stars surfing the bar in Type-III discs

The radial profiles of stars in disc galaxies are observed to be either purely exponential (Type-I), truncated (Type-II) or antitruncated (Type-III) exponentials. Controlled formation simulations of isolated galaxies can reproduce all of these profile types by varying a single parameter, the initial halo spin. In this paper, we examine these simulations in more detail in an effort to identify the physical mechanism that leads to the formation of Type-III profiles. The stars in the antitruncated outskirts of such discs are now on eccentric orbits, but were born on near-circular orbits at much smaller radii. We show that, and explain how, they were driven to the outskirts via non-linear interactions with a strong and long-lived central bar, which greatly boosted their semimajor axis but also their eccentricity. While bars have been known to cause radial heating and outward migration to stellar orbits, we link this effect to the formation of Type-III profiles. This predicts that the antitruncated parts of galaxies have unusual kinematics for disc-like stellar configurations: high radial velocity dispersions and slow net rotation. Whether such discs exist in nature, can be tested by future observations

Triggering of merger-induced starbursts by the tidal field of galaxy groups and clusters

Star formation in galaxies is for a part driven by galaxy mergers. At low redshift, star formation activity is low in high-density environments like groups and clusters, and the star formation activity of galaxies increases with their isolation. This star formation-density relation is observed to be reversed at z~1, which is not explained by theoretical models so far. We study the influence of the tidal field of a galaxy group or cluster on the star formation activity of merging galaxies, using N-body simulations including gas dynamics and star formation. We find that the merger-driven star formation is significantly more active in the vicinity of such cosmological structures compared to mergers in the field. The large-scale tidal field can thus enhance the activity of galaxies in dense cosmic structures, and should be particularly efficient at high redshift before quenching processes take effect in the densest regions.Comment: MNRAS Letters, in pres

BAYES-LOSVD: A Bayesian framework for non-parametric extraction of the line-of-sight velocity distribution of galaxies

We introduce BAYES-LOSVD, a novel implementation of the non-parametric extraction of line-of-sight velocity distributions (LOSVDs) in galaxies. We employed Bayesian inference to obtain robust LOSVDs and associated uncertainties. Our method relies on a principal component analysis to reduce the dimensionality on the set of templates required for the extraction and thus increase the performance of the code. In addition, we implemented several options to regularise the output solutions. Our tests, conducted on mock spectra, confirm the ability of our approach to model a wide range of LOSVD shapes, overcoming limitations of the most widely used parametric methods (e.g., Gauss-Hermite expansion). We present examples of LOSVD extractions for real galaxies with known peculiar LOSVD shapes, including NGC 4371, IC 0719, and NGC 4550, using MUSE and SAURON integral-field unit (IFU) data. Our implementation can also handle data from other popular IFU surveys (e.g., ATLAS3D, CALIFA, MaNGA, SAMI)

The [α/Fe]-[Fe/H] relation in the E-MOSAICS simulations: its connection to the birth place of globular clusters and the fraction of globular cluster field stars in the bulge

The {\alpha}-element abundances of the globular cluster (GC) and field star populations of galaxies encode information about the formation of each of these components. We use the E-MOSAICS cosmological simulations of ~L* galaxies and their GCs to investigate the [{\alpha}/Fe]-[Fe/H] distribution of field stars and GCs in 25 Milky Way-mass galaxies. The [{\alpha}/Fe]-[Fe/H] distribution go GCs largely follows that of the field stars and can also therefore be used as tracers of the [{\alpha}/Fe]-[Fe/H] evolution of the galaxy. Due to the difference in their star formation histories, GCs associated with stellar streams (i.e. which have recently been accreted) have systematically lower [{\alpha}/Fe] at fixed [Fe/H]. Therefore, if a GC is observed to have low [{\alpha}/Fe] for its [Fe/H] there is an increased probability that this GC was accreted recently alongside a dwarf galaxy. There is a wide range of shapes for the field star [{\alpha}/Fe]-[Fe/H] distribution, with a notable subset of galaxies exhibiting bimodal distributions, in which the high [{\alpha}/Fe] sequence is mostly comprised of stars in the bulge, a high fraction of which are from disrupted GCs. We calculate the contribution of disrupted GCs to the bulge component of the 25 simulated galaxies and find values between 0.3-14 per cent, where this fraction correlates with the galaxy's formation time. The upper range of these fractions is compatible with observationally-inferred measurements for the Milky Way, suggesting that in this respect the Milky Way is not typical of L* galaxies, having experienced a phase of unusually rapid growth at early times

Local variations of the Stellar Velocity Ellipsoid-I: the disc of galaxies in the Auriga simulations

DW-M and JF-B acknowledge support through the RAVET project by the grant PID2019-107427GB-C32 from the Spanish Ministry of Science, Innovation and Universities (MCIU), and through the IAC project TRACES which is partially supported through the state budget and the regional budget of the Consejer ' ia de Econom ' ia, Industria, Comercio y Conocimiento of the Canary Islands Autonomous Community.The connection between the Stellar Velocity Ellipsoid (SVE) and the dynamical evolution of galaxies has been a matter of debate in the last years and there is no clear consensus whether different heating agents (e.g. spiral arms, giant molecular clouds, bars and mergers) leave clear detectable signatures in the present day kinematics. Most of these results are based on a single and global SVE and have not taken into account that these agents do not necessarily equally affect all regions of the stellar disc. We study the two-dimensional (2D) spatial distribution of the SVE across the stellar discs of Auriga galaxies, a set of high resolution magnetohydrodynamical cosmological zoom-in simulations, to unveil the connection between local and global kinematic properties in the disc region.We find very similar, global, σz/σr = 0.80 ± 0.08 values for galaxies of different Hubble types. This shows that the global properties of the SVE at z = 0 are not a good indicator of the heating and cooling events experienced by galaxies.We also find that similar σz/σr radial profiles are obtained through different combinations of σz and σr trends: at a local level, the vertical and radial components can evolve differently, leading to similar σz/σr profiles at z = 0. By contrast, the 2D spatial distribution of the SVE varies a lot more from galaxy to galaxy. Present day features in the SVE spatial distribution may be associated with specific interactions such as fly-by encounters or the accretion of low mass satellites even in the cases when the global SVE is not affected. The stellar populations decomposition reveals that young stellar populations present colder and less isotropic SVEs and more complex 2D distributions than their older and hotter counterparts.Spanish Ministry of Science, Innovation and Universities (MCIU) PID2019-107427GB-C32Consejería de Economía, Industria, Comercio y Conocimiento of the Canary Islands Autonomous Communit

Field ecology of sylvatic Rhodnius populations (Heteroptera, Triatominae): risk factors for palm tree infestation in western Ecuador.

Most Rhodnius species (Triatominae) are primarily associated with palm trees. They maintain enzootic Trypanosoma cruzi transmission and are responsible for human infection (causing Chagas disease) through the Neotropics. Assessing whether individual palm traits (ecological and/or botanical) may increase the risk of palm infestation by triatomines is relevant in areas where bugs invade houses flying from peridomestic palms. We developed a novel fieldwork approach with that objective, and applied it to study infestation by sylvatic Rhodnius ecuadoriensis in 110 tagua palms (Phytelephas aequatorialis). Palm infestation (23% overall) was non-randomly distributed in our sample. Palms located in anthropic landscapes were frequently infested (>27%, n=92), whereas no bugs were collected from palms surveyed within forest remnants (n=18; P=0.01). The presence of abundant decaying vegetable matter (P=0.001) and (to a lesser extent) epiphytic plants (P=0.049) on palm crowns and stems increased the probability of infestation and was positively correlated with the apparent density of bug colonies (R2=0.68). A trend towards higher infestation rates in male palms (34% vs. 18%) could relate to female palm management (removal of infrutescences and vegetable debris) in areas where palm seeds are harvested. An outline of 'risk palm ecotopes' and environmental management-based strategies for the control of peridomestic, palm tree-living vector populations are proposed

The Relationship Between Mono-Abundance and Mono-Age Stellar Populations in the Milky Way Disk

Studying the Milky Way disk structure using stars in narrow bins of [Fe/H] and [α/Fe] has recently been proposed as a powerful method to understand the Galactic thick and thin disk formation. It has been assumed so far that these mono-abundance populations (MAPs) are also coeval, or mono-age, populations. Here we study this relationship for a Milky Way chemodynamical model and show that equivalence between MAPs and mono-age populations exists only for the high-[α/Fe] tail, where the chemical evolution curves of different Galactic radii are far apart. At lower [α/Fe]-values an MAP is composed of stars with a range in ages, even for small observational uncertainties and a small MAP bin size. Due to the disk inside-out formation, for these MAPs younger stars are typically located at larger radii, which results in negative radial age gradients that can be as large as 2 Gyr kpc−1. Positive radial age gradients can result for MAPs at the lowest [α/Fe] and highest [Fe/H] end. Such variations with age prevent the simple interpretation of observations for which accurate ages are not available. Studying the variation with radius of the stellar surface density and scale height in our model, we find good agreement to recent analyses of the APOGEE red-clump (RC) sample when 1–4 Gyr old stars dominate (as expected for the RC). Our results suggest that the APOGEE data are consistent with a Milky Way model for which mono-age populations flare for all ages. We propose observational tests for the validity of our predictions and argue that using accurate age measurements, such as from asteroseismology, is crucial for putting constraints on Galactic formation and evolution