Search

Ultra metal-poor stars ([Fe/H] < -4.0) are very rare, and finding them is a challenging task. Both narrow-band photometry and low-resolution spectroscopy have been useful tools for identifying candidates, and in this work we combine both approaches. We cross-matched metallicity-sensitive photometry from the Pristine survey with the low-resolution spectroscopic LAMOST database, and re-analysed all LAMOST spectra with [Fe/H]_Pristine < -2.5. We find that ~1/3rd of this sample (selected without [Fe/H]_Pristine quality cuts) also have spectroscopic [Fe/H] < -2.5. From this sample, containing many low signal-to-noise (S/N) spectra, we selected eleven stars potentially having [Fe/H] < -4.0 or [Fe/H] < -3.0 with very high carbon abundances, and we performed higher S/N medium-resolution spectroscopic follow-up with OSIRIS on the 10.4m Gran Telescopio Canarias (GTC). We confirm their extremely low metallicities, with a mean of [Fe/H] = -3.4 and the most metal-poor star having [Fe/H]= -3.8. Three of these are clearly carbon-enhanced metal-poor (CEMP) stars with +1.65 < [C/Fe] < +2.45. The two most carbon-rich stars are either among the most metal-poor CEMP-s stars or the most carbon-rich CEMP-no stars known, the third is likely a CEMP-no star. We derived orbital properties for the OSIRIS sample and find that only one of our targets can be confidently associated with known substructures/accretion events, and that three out of four inner halo stars have prograde orbits. Large spectroscopic surveys may contain many hidden extremely and ultra metal-poor stars, and adding additional information from e.g. photometry as in this work can uncover them more efficiently and confidently., 13 pages, accepted for publication in MNRAS

Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate/geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to study the nature (i.e. geometry, rates) of mass-loss episodes. In this context, yellow hypergiants are great targets. Methods. We analyse a large set of optical/near-infrared data, in spectroscopic and photometric (X-shooter/VLT), spectropolarimetric (ISIS/WHT), and interferometric GRAVITY-AMBER/VLTI) modes, toward the yellow hypergiant IRAS 17163-3907. We present the first model-independent reconstructed images of IRAS 17163-3907 at these wavelengths at milli-arcsecond scales. Lastly, we apply a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published VISIR images at 8.59 {\mu}m, 11.85 {\mu}m and 12.81 {\mu}m simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 {\mu}m show that the Br{\gamma} emission is more extended and asymmetric than the Na i and the continuum emission. In addition to the two known shells surrounding IRAS 17163-3907 we report on the existence of a third hot inner shell with a maximum dynamical age of only 30 yr. Conclusions. The interpretation of the presence of Na i emission at closer distances to the star compared to Br{\gamma} has been a challenge in various studies. We argue that the presence of a pseudophotosphere is not needed, but it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform the theories on mass-loss mechanisms, which is a topic still under debate. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms., Comment: 26 pages, 22 Figures, 6 tables, accepted for publication in A&A. Original abstract will be available after publication

One of the major goals for astronomy in the next decades is the remote search for biosignatures (i.e.\ the spectroscopic evidence of biological activity) in exoplanets. Here, we adopt a Bayesian statistical framework to discuss the implications of such future searches, both in the case when life is detected, and when no definite evidence is found. We show that even a single detection of biosignatures in the vicinity of our stellar system, in a survey of similar size to what will be obtainable in the next two decades, would affect significantly our prior belief on the frequency of life in the universe, even starting from a neutral or pessimistic stance. In particular, after such discovery, an initially agnostic observer would be led to conclude that there are more than $10^5$ inhabited planets in the galaxy with a probability exceeding $95$\%. However, this conclusion would be somewhat weakened by the viability of transfer of biological material over interstellar distances, as in panspermia scenarios. Conversely, the lack of significant evidence of biosignatures would have little effect, leaving the assessment of the abundance of life in the galaxy still largely undetermined., Comment: Published on PNAS

We present DARKFLUX, a software tool designed to analyze indirect-detection signatures for next-generation models of dark matter (DM) with multiple annihilation channels. Version 1.0 of this tool accepts user-generated models with 2 -> 2 tree-level dark matter annihilation to pairs of Standard Model (SM) particles and analyzes DM annihilation to gamma rays. The tool consists of three modules, which can be run in a loop in order to scan over DM mass if desired, (I) The annihilation fraction module calls an internal installation of MADDM, a dark matter phenomenology plugin for the Monte Carlo event generator MADGRAPH5_AMC@NLO, to compute the thermally averaged cross section (i) for each annihilation channel chi chi((chi) over bar, chi(dagger)) -> i is an element of {SM, SM}. The module then computes the fractional annihilation rate (annihilation fraction) into each channel., (II) The flux module combines the flux spectrum from each annihilation channel, weighted by the appropriate annihilation fractions, to compute the total flux at Earth due to DM annihilation. In DARKFLUX v1.0, this module specifically computes the gamma-ray flux for each channel using the publicly available PPPC4DMID tables., (III) The analysis module compares the total flux to observational data and computes the upper limit at 95% confidence level (CL) on the total thermally averaged DM annihilation cross section. In DARKFLUX v1.0, this module compares the total gamma-ray flux to a joint-likelihood analysis of fifteen dwarf spheroidal galaxies (dSphs) analyzed by the Fermi-LAT collaboration., DARKFLUX v1.0 automatically provides data tables and can plot the output of these three modules. In this manual, we briefly motivate this indirect-detection computer tool and review the essential DM physics. We then describe the several modules of DARKFLUX in greater detail. Finally, we show how to install and run DARKFLUX and provide two worked examples demonstrating its capabilities. DARKFLUX is available on GitHub at https://github.com/carpenterphysics/DarkFlux. (C) 2022 Elsevier B.V. All rights reserved.

There is now ample evidence that dust is already present in abundance at high z. However, given the faintness of distant galaxies in the optical and the NIR, datasets are still limited and how the dust affects the emerging radiation of galaxies at very high redshift is not yet fully understood. Using the ALPINE survey, our objective is to quantify the dust attenuation properties in galaxies at z=4.4-5.5, and in particular the shape of their attenuation curve. Using the CIGALE code, we model the stellar populations and their interaction with the dust in order to measure some of the physical properties of a subsample of 23 main-sequence ALPINE galaxies. We find that the attenuation curves span a broad range of properties, from curves that are much steeper than the SMC extinction curve, to shallower than the starburst attenuation curve. The shape of the attenuation curves strongly depends on the V-band attenuation. Galaxies with the lowest attenuation also present the steepest curves. The steepness of such curves is probably the consequence of the combination of the intrinsic physical properties of the dust, the relative distribution of stars and dust in the interstellar medium, and the differential reddening. The broad range of attenuation curves found at z~5 shows that no single attenuation curve is appropriate for main sequence galaxies and that assuming a fixed curve can lead to large errors, for instance in the interpretation and use of the IRX-beta diagram, if SED modeling is not feasible. Great caution should be exercised when correcting high redshift galaxies for the presence of dust using the UV slope beta as it can affect the estimation of both SFR and stellar mass even at low V-band attenuation due to the steepness of the attenuation curve. However, when SED modeling can be used, the impact of the choice of the attenuation curve on the SFR and the stellar mass is limited., Comment: 20 pages, 14 figures. Accepted for publication in A&A

The Cepheid period-luminosity (PL) relation (or Leavitt law) has served as the first rung of the most widely used extragalactic distance ladder and is central to the determination of the local value of the Hubble constant ($H_0$). We investigate the influence of metallicity on Cepheid brightness, a term that significantly improves the overall fit of the distance ladder, to better define its wavelength dependence. To this aim, we compare the PL relations obtained for three Cepheid samples having distinct chemical composition (in the Milky Way and Magellanic Clouds) and focusing on the use of improved and recent data while covering a metallicity range of about 1 dex. We estimate the metallicity effect (hereafter $\gamma$) in 15 filters from mid-IR to optical wavelengths, including five Wesenheit indices, and we derive a significant metallicity term in all filters, in agreement with recent empirical studies and models, in the sense of metal-rich Cepheids being brighter than metal-poor ones. We describe the contribution of various systematic effects in the determination of the $\gamma$ term. We find no evidence of $\gamma$ changing over the wavelength range $0.5-4.5 \, \rm \mu m$, indicating that the main influence of metallicity on Cepheids is in their luminosity rather than color. Finally, we identify factors that sharpen the empirical constraints on the metallicity term over past studies, including corrections for the depth of the Magellanic Clouds, better-calibrated Cepheid photometry, improved Milky Way extinction estimates, and revised and expanded metallicity measurements in the LMC., Comment: Accepted to ApJ

We have conducted a search for open clusters in the vicinity of classical Galactic Cepheids based on high-quality astrometry from the third data release (DR3) of the ESA mission Gaia to improve the calibration of the Leavitt law (LL). Our approach requires no prior knowledge of existing clusters, allowing us to both detect new host clusters and cross-check previously reported associations. Our Gold sample consists of 34 Cepheids residing in 28 open clusters, including 27 fundamental mode and 7 overtone Cepheids. Three new bona fide cluster Cepheids are reported (V0378 Cen, ST Tau, and GH Lup) and the host cluster identifications for three others (VW Cru, IQ Nor, and SX Vel) are corrected. The fraction of Cepheids occurring in open clusters within 2 kpc of the Sun is $f_{CC,2kpc} = 0.088^{+0.029}_{-0.019}$. By combining cluster and field Cepheids, we calibrate the LL for several individual photometric passbands, together with reddening-free Wesenheit magnitudes based on Gaia and HST photometry, while solving for the residual offset applicable to Cepheid parallaxes, $\Delta \varpi_{\mathrm{Cep}}$. The most direct comparison of our results with the SH0ES distance ladder yields excellent ($0.3\sigma$) agreement for both the absolute magnitude of a 10d solar metallicity Cepheid in the near-IR HST Wesenheit magnitudes, $M_{H,1}^W=-5.914\pm 0.017$ mag, and the residual parallax offset, $\Delta \varpi_{\mathrm{Cep}}=-13 \pm 5\,\mu$as. Using the larger sample of 26 Gold cluster Cepheids and $225$ MW Cepheids with recent Gaia DR3 astrometry and photometry, we determine at solar metallicity $M_{G,1}^W = -6.004 \pm 0.019$\,mag and $\Delta \varpi_{\mathrm{Cep}}=-19 \pm 3\,\mu$as. These results mark the currently most accurate absolute calibrations of the Cepheid luminosity scale based purely on observations of MW Cepheids., Comment: 24 pages, Version accepted in A&A

Context. The [Mg/Fe] abundance ratios are a fundamental fossil signature to trace the chemical evolution of the disc. Despite of the huge observational and theoretical effort, discrepancies between models and data are still present and several explanations have been put forward to explain the [$\alpha$/Fe] bimodality. Aims. In this work, we take advantage of a new AMBRE:HARPS dataset, which provides new and more precise [Mg/Fe] estimations, as well as reliable stellar ages for a subsample of stars, to study the evolution of the solar neighbourhood. Methods. The above data are compared with detailed chemical evolution models for the Milky Way, exploring the most used prescriptions for stellar yields and different formation scenarios for the Galactic disc, i.e. the delayed two-infall and the parallel model, also including prescriptions for stellar radial migration. Results. We see that most of the stellar yields struggle to reproduce the observed trend of the data and that semi-empirical yields are still the best to describe the [Mg/Fe] evolution in the thick and thin discs. In particular, most of the yields still predict a steeper decrease of the [Mg/Fe] ratio at high metallicity than what is shown by the data. The bulk of the data are well reproduced by the parallel and two-infall scenarios, but both scenarios have problems in explaining the most metal-rich and metal-poor tails of the low-$\alpha$ data. These tails can be explained in light of radial migration from inner and outer disc regions, respectively. Conclusions. Despite of the evidence of stellar migration, it is difficult to estimate the actual contribution of stars from other parts of the disc to the solar vicinity. However, the comparison between data and models suggests that peculiar histories of star formation, such as that of the two-infall model, are still needed to reproduce the observed distribution of stars., Comment: 16 pages, 16 figures. Accepted for publication in Astronomy and Astrophysics (A&A)

We report the identification of a set of old super metal-rich dwarf stars with orbits of low eccentricity that reach a maximum height from the Galactic plane between ~0.5-1.5 kpc. We discuss their properties to understand their origins. We use data from the internal data release 6 of the Gaia-ESO Survey. We selected stars observed at high resolution with abundances of 21 species of 18 individual elements. We apply hierarchical clustering to group the stars with similar chemical abundances within the complete chemical abundance space. According to their chemical properties, this set of super metal-rich stars can be arranged into five subgroups. Four seem to follow a chemical enrichment flow, where nearly all abundances increase in lockstep with Fe. The fifth subgroup shows different chemical characteristics. All subgroups have the following features: median ages of the order of 7-9 Gyr, Solar or sub-Solar [Mg/Fe] ratios, maximum height between 0.5-1.5 kpc, low eccentricities, and a detachment from the expected metallicity gradient with guiding radius. The high metallicity of our stars is incompatible with a formation in the Solar neighbourhood. Their dynamic properties agree with theoretical expectations that these stars travelled from the inner Galaxy due to blurring and, most importantly, to churning. We suggest that most of this population's stars originated in the Milky Way's inner regions (inner disc and/or the bulge) and later migrated to the Solar neighbourhood. The region from where the stars originated had a complex chemical enrichment history, with contributions from supernovae types Ia and II and possibly asymptotic giant branch stars., Comment: 18 pages, 13 figures. Accepted for publication in A&A. Version after language proofs. The catalogue remains to be released. Abridged abstract to fit arxiv's requirements

Numerical simulations of dwarf galaxies have so far failed to reproduce the observed metallicity-luminosity relation, down to the regime of ultra-faint dwarfs (UFDs). We address this issue by exploring how the first generations of metal-free stars (Pop III) could help increase the mean metallicity ([Fe/H]) of those small and faint galaxies. We ran zoom-in chemo-dynamical simulations of 19 halos extracted from a Λ Cold Dark Matter (CDM) cosmological box and followed their evolution down to redshift z = 0. Models were validated not only on the basis of galaxy global properties, but also on the detailed investigation of the stellar abundance ratios ([α/Fe]). We identified the necessary conditions for the formation of the first stars in mini-halos and derived constraints on the metal ejection schemes. The impact of Pop III stars on the final metallicity of UFDs was evaluated by considering different stellar mass ranges for their initial mass function (IMF), the influence of pair-instability supernovae (PISNe), and their energetic feedback, as well as the metallicity threshold that marks the transition from the first massive stars to the formation of low-mass long-lived stars. The inclusion of Pop III stars with masses below 140 M⊙, and a standard IMF slope of −1.3 does increase the global metallicity of UFDs, although these are insufficient to resolve the tension with observations. The PISNe with progenitor masses above 140 M⊙ do allow the metal content of UFDs to further increase. However, as PISNe are very rare and sometimes absent in the faintest UFDs, they have a limited impact on the global faint end of the metallicity-luminosity relation. Despite a limited number of spectroscopically confirmed members in UFDs, which make the stellar metallicity distribution of some UFDs uncertain, our analysis reveals that this is essentially the metal-rich tail that is missing in the models. The remaining challenges are thus both observational and numerical: (i) to extend high-resolution spectroscopy data samples and confirm the mean metallicity of the faintest UFDs; and (ii) to explain the presence of chemically enriched stars in galaxies with very short star formation histories.

The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemodynamical mapping of the Milky Way. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) “Ancillary Science Programs” competitively awarded to Sloan Digital Sky Survey IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5 yr expansion of the survey, known as the Bright Time Extension (BTX), made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The BTX permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new data sets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, F. Santana et al. (submitted; AAS29036), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

We combine JWST observations with ALMA CO and VLT-MUSE H$\alpha$ data to examine off-spiral arm star formation in the face-on, grand-design spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3-4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21$\mu$m and MUSE H$\alpha$ of around 40$^\circ$ (CO-rich) and 55$^\circ$ (CO-poor) from the spiral arm. The star formation rate is higher in the regions of the spurs near to spiral arms, but the star formation efficiency appears relatively constant. Given the spiral pattern speed and rotation curve of this galaxy and assuming material exiting the arms undergoes purely circular motion, these offsets would be reached in 100-150 Myr, significantly longer than the 21$\mu$m and H$\alpha$ star formation timescales (both, Comment: 11 pages, 5 Figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ. Updated to include missing author and published paper references

(abridged) Star formation in the outer Galaxy, i.e., outside of the Solar circle, has been lightly studied in part due to low CO brightness of molecular clouds linked with the negative metallicity gradient. Recent infrared surveys provide an overview of dust emission in large sections of the Galaxy, but suffer from cloud confusion and poor spatial resolution at far-infrared wavelengths. We aim to develop a methodology to identify and classify young stellar objects (YSOs) in star-forming regions in the outer Galaxy, and use it to solve a long-standing confusion with the distance and evolutionary status of IRAS 22147+5948. We use Support Vector Machine learning algorithm to complement standard color-color and color-magnitude diagrams in search for YSOs in the IRAS 22147 region using publicly available data from the `Spitzer Mapping of the Outer Galaxy' survey. The agglomerative hierarchical clustering algorithm is used to identify clusters, along with the Robitaille et al. (2017) code to calculate physical properties of individual YSOs. We identify 13 Class I and 13 Class II YSO candidates using the color-color diagrams, and additional 2 and 21 sources, respectively, using the machine learning techniques. Spectral energy distributions of 23 sources are modelled with a star and a passive disk, corresponding to Class II objects. Models of 3 sources include envelopes typical for Class I objects. The objects are grouped in 2 clusters located at the distance of ~2.2 kpc, and 5 clusters at ~5.6 kpc. The spatial extent of CO, radio continuum, and dust emission confirms the origin of YSOs in two distinct star-forming regions along a similar line-of-sight. The outer Galaxy might serve as a unique laboratory of star formation across environments on condition that complementary methods and ancillary data are used to properly account for cloud confusion and distance uncertainties., 29 pages, 17 figures, 9 tables

In this fourth article on weighing the Galactic disk using the shape of the phase-space spiral, we have tested our method on a billion particle three-dimensional N-body simulation, comprised of a Milky Way like host galaxy and a merging dwarf satellite. The main purpose of this work was to test the validity of our model’s fundamental assumptions that the spiral inhabits a locally static and vertically separable gravitational potential. These assumptions might be compromised in the complex kinematic system of a disturbed three-dimensional disk galaxy; in fact, the statistical uncertainty and any potential biases related to these assumptions are expected to be amplified for this simulation, which differs from the Milky Way in that it is more strongly perturbed and has a phase-space spiral that inhabits higher vertical energies. We constructed 44 separate data samples from different spatial locations in the simulated host galaxy. Our method produced accurate results for the vertical gravitational potential of these 44 data samples, with an unbiased distribution of errors with a standard deviation of 7%. We also tested our method under severe and unknown spatially dependent selection effects, also with robust results; this sets it apart from traditional dynamical mass measurements that are based on the assumption of a steady state and which are highly sensitive to unknown or poorly modelled incompleteness. Hence, we will be able to make localised mass measurements of distant regions in the Milky Way disk, which would otherwise be compromised by complex and poorly understood selection effects., AW acknowledges support from the Carlsberg Foundation via a Semper Ardens grant (CF15-0384). JH is supported by a Flatiron Research Fellowship at the Flatiron institute, which is supported by the Simons Foundation. CL acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 852839). GM acknowledges funding from the Agence Nationale de la Recherche (ANR project ANR-18-CE31-0006 and ANR-19-CE31-0017) and from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 834148). This work made use of an HPC facility funded by a grant from VILLUM FONDEN (projectnumber 16599). This research utilised the following open-source Python packages: MATPLOTLIB (Hunter 2007), NUMPY (Harris et al. 2020), SCIPY (Virtanen et al. 2020), PANDAS (McKinney et al. 2010), TENSORFLOW (Abadi et al. 2015).

We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities., 29 pages, 22 figures, accepted for publication in A&A

We combine the power of two stream-searching tools, STREAMFINDER and StarGO applied to the Gaia EDR3 data, to detect stellar debris belonging to the Cetus stream system that forms a complex, nearly polar structure around the Milky Way. In this work, we find the southern extensions of the northern Cetus stream as the Palca stream and a new southern stream, which overlap on the sky but have different distances. These two stream wraps extend over more than $\sim100\deg$ on the sky ($-60\deg, Comment: accepted to ApJ. The N-body model of the Cetus system is published at https://zenodo.org/record/5771585

Aims. We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of different aspects of the Milky Way structure and evolution and we provide, at the same time, a description of several practical aspects of the data and examples of their usage., Methods. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. In this direction, the Gaia astrometric data alone enable the calculation of the vertical and azimuthal velocities; also, the extinction is relatively low compared to other directions in the Galactic plane. We then explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1., Results. With the improved astrometry and photometry of EDR3, we find that: (i) the dynamics of the Galactic disc are very complex with oscillations in the median rotation and vertical velocities as a function of radius, vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; (ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; (iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; (iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and (v) the open clusters Berkeley 29 and Saurer 1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits., Conclusions. Even with our simple preliminary exploration of the Gaia EDR3, we demonstrate how, once again, these data from the European Space Agency are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history.

It is still poorly constrained how the densest phase of the interstellar medium varies across galactic environment. A large observing time is required to recover significant emission from dense molecular gas at high spatial resolution, and to cover a large dynamic range of extragalactic disc environments. We present new NOrthern Extended Millimeter Array (NOEMA) observations of a range of high critical density molecular tracers (HCN, HNC, HCO+) and CO isotopologues (13CO, C18O) towards the nearby (11.3 Mpc), strongly barred galaxy NGC 3627. These observations represent the current highest angular resolution (1.85"; 100 pc) map of dense gas tracers across a disc of a nearby spiral galaxy, which we use here to assess the properties of the dense molecular gas, and their variation as a function of galactocentric radius, molecular gas, and star formation. We find that the HCN(1-0)/CO(2-1) integrated intensity ratio does not correlate with the amount of recent star formation. Instead, the HCN(1-0)/CO(2-1) ratio depends on the galactic environment, with differences between the galaxy centre, bar, and bar end regions. The dense gas in the central 600 pc appears to produce stars less efficiently despite containing a higher fraction of dense molecular gas than the bar ends where the star formation is enhanced. In assessing the dynamics of the dense gas, we find the HCN(1-0) and HCO+(1-0) emission lines showing multiple components towards regions in the bar ends that correspond to previously identified features in CO emission. These features are co-spatial with peaks of Halpha emission, which highlights that the complex dynamics of this bar end region could be linked to local enhancements in the star formation., Comment: 25 pages, 12 figures; accepted for publication in MNRAS

Elemental abundances of the most metal-poor stars reflect the conditions in the early Galaxy and the properties of the first stars. We present a spectroscopic follow-up of two ultra metal-poor stars ([Fe/H], Comment: 17 pages, 11 figures, accepted for publication in MNRAS

We present a new method for inferring the gravitational potential of the Galactic disk, using the time-varying structure of a phase-space spiral in the $(z,w)$-plane (where $z$ and $w$ represent vertical position and vertical velocity). Our method of inference extracts information from the shape of the spiral and disregards the bulk density distribution that is usually used to perform dynamical mass measurements. In this manner, it is complementary to traditional methods that are based on the assumption of a steady state. Our method consists of fitting an analytical model for the phase-space spiral to data, where the spiral is seen as a perturbation of the stellar number density in the $(z,w)$-plane. We tested our method on one-dimensional simulations, which were initiated in a steady state and then perturbed by an external force similar to that of a passing satellite. We were able to retrieve the true gravitational potentials of the simulations with high accuracy. The gravitational potential at 400-500 parsec distances from the disk mid-plane was inferred with an error of only a few percent. This is the first paper of a series in which we plan to test and refine our method on more complex simulations, as well as apply our method to Gaia data., Comment: 16 pages, 12 figures

We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more stars than current spectroscopic MDF determinations. We find that Eri II has a mean metallicity of [Fe/H]$=$-2.50$^{+0.07}_{-0.07}$ and a dispersion of $��_{\mbox{[Fe/H]}}=0.42^{+0.06}_{-0.06}$, which are consistent with spectroscopic MDFs, though more precisely constrained owing to a larger sample. We identify a handful of extremely metal-poor star candidates (EMP; [Fe/H] $< -3$) that are marginally bright enough for spectroscopic follow up. Eri II's MDF appears well-described by a leaky box chemical evolution model. We also compute an updated orbital history for Eri II using Gaia eDR3 proper motions, and find that it is likely on first infall into the Milky Way. Our findings suggest that Eri II underwent an evolutionary history similar to that of an isolated galaxy. Compared to MDFs for select cosmological simulations of similar mass galaxies, we find that Eri II has a lower fraction of stars with [Fe/H] $< -3$, though such comparisons should currently be treated with caution due to a paucity of simulations, selection effects, and known limitations of Ca H&K for EMPs. This study demonstrates the power of deep HST CaHK imaging for measuring the MDFs of UFDs., Resubmitted to ApJ with revisions following a positive referee report of the initial draft; 29 pages, 19 figures, and 3 tables, including appendix. Comments welcome!

We present a catalogue of 362 million stellar parameters, distances, and extinctions derived from Gaia's early third data release (EDR3) cross-matched with the photometric catalogues of Pan-STARRS1, SkyMapper, 2MASS, and AllWISE. The higher precision of the Gaia EDR3 data, combined with the broad wavelength coverage of the additional photometric surveys and the new stellar-density priors of the {\tt StarHorse} code allow us to substantially improve the accuracy and precision over previous photo-astrometric stellar-parameter estimates. At magnitude $G=14\, (17)$, our typical precisions amount to 3% (15%) in distance, 0.13 mag (0.15 mag) in $V$-band extinction, and 140 K (180 K) in effective temperature. Our results are validated by comparisons with open clusters, as well as with asteroseismic and spectroscopic measurements, indicating systematic errors smaller than the nominal uncertainties for the vast majority of objects. We also provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps, and extensive stellar density maps that reveal detailed substructures in the Milky Way and beyond. The new density maps now probe a much greater volume, extending to regions beyond the Galactic bar and to Local Group galaxies, with a larger total number density. We publish our results through an ADQL query interface ({\tt gaia.aip.de}) as well as via tables containing approximations of the full posterior distributions. Our multi-wavelength approach and the deep magnitude limit make our results useful also beyond the next Gaia release, DR3., Comment: A&A, accepted, 27 pages, 25 figures. Fig. 5 has changed with respect to v1 (post-processing bug corrected). For various options to access the data see https://data.aip.de/projects/starhorse2021.html

The ultraviolet (UV) to sub-millimetre (submm) spectral energy distribution of galaxies can be roughly divided into two sections: the stellar emission (attenuated by dust) at UV to near-infrared wavelengths and dust emission at longer wavelengths. In Dobbels et al. (2020), we show that these two sections are strongly related, and we can predict the global dust properties from the integrated UV to mid-infrared emission with the help of machine learning techniques. We investigate if these machine learning techniques can also be extended to resolved scales. Our aim is to predict resolved maps of the specific dust luminosity, specific dust mass, and dust temperature starting from a set of surface brightness images from UV to mid-infrared wavelengths. We used a selection of nearby galaxies retrieved from the DustPedia sample, in addition to M31 and M33. These were convolved and resampled to a range of pixel sizes, ranging from 150 pc to 3 kpc. We trained a random forest model which considers each pixel individually. We find that the predictions work well on resolved scales, with the dust mass and temperature having a similar root mean square error as on global scales (0.32 dex and 3.15 K on 18" scales respectively), and the dust luminosity being noticeably better (0.11 dex). We find no significant dependence on the pixel scale. Predictions on individual galaxies can be biased, and we find that about two-thirds of the scatter can be attributed to scatter between galaxies (rather than within galaxies). A machine learning approach can be used to create dust maps, with its resolution being only limited to the input bands, thus achieving a higher resolution than Herschel. These dust maps can be used to improve global estimates of dust properties, they can lead to a better estimate of dust attenuation, and they can be used as a constraint on cosmological simulations that trace dust., Comment: 16 pages, 13 figures, accepted by A&A

Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions., Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar systembarycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution., Methods. Theeffect of the acceleration was obtained as a part of the general expansion of the vector field of proper motions in vector spherical harmonics (VSH). Various versions of the VSH fit and various subsets of the sources were tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution were used to get a better idea of the possible systematic errors in the estimate., Results. Our best estimate of the acceleration based on Gaia EDR3 is (2.32 +/- 0.16) x 10(-10) m s(-2) (or 7.33 +/- 0.51 km s(-1) Myr-1) towards alpha = 269.1 degrees +/- 5.4 degrees, delta = -31.6 degrees +/- 4.1 degrees, corresponding to a proper motion amplitude of 5.05 +/- 0.35 mu as yr(-1). This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 mu as yr(-1).

RR Lyrae stars are an important and widely used tracer of the most ancient populations of our Galaxy, mainly due to their standard candle nature. The availability of large scale surveys of variable stars is allowing us to trace the structure of our entire Galaxy, even in previously inaccessible areas like the Galactic disc. In this work we aim to provide an empirical assessment of the completeness of the three largest RR Lyrae catalogues available: Gaia DR2, PanSTARRS-1 and ASAS-SN-II. Using a joint probabilistic analysis of the three surveys we compute 2D and 3D completeness maps in each survey's full magnitude range. At the bright end (G20deg); ASAS-SN-II has the best completeness at low latitude for RRab and at all latitudes for RRc. At the faint end (G>13), Gaia DR2 is the most complete catalogue for both RR Lyrae types, at any latitude, with median completeness rates of 95% (RRab) and >85% (RRc) outside the ecliptic plane (|��|>25deg). We confirm a high and uniform completeness of PanSTARRS-1 RR Lyrae at 91% (RRab) and 82% (RRc) down to G~18, and provide the first estimate of its completeness at low galactic latitude (|b, 15 pages, 14 figures, 4+10 tables. Accepted for publication at MNRAS. Completeness maps and Python utils to plot and/or recompute them, available at https://github.com/cmateu/rrl_completeness

In Meidt et al. (2018), we showed that gas kinematics on the scale of individual molecular clouds are not dominated by self-gravity but also track a component that originates with orbital motion in the potential of the host galaxy. This agrees with observed cloud line widths, which show systematic variations from virial motions with environment, pointing at the influence of the galaxy potential. In this paper, we hypothesize that these motions act to slow down the collapse of gas and so help regulate star formation. Extending the results of Meidt et al. (2018), we derive a dynamical collapse timescale that approaches the free-fall time only once the gas has fully decoupled from the galactic potential. Using this timescale we make predictions for how the fraction of free-falling, strongly self-gravitating gas varies throughout the disks of star-forming galaxies. We also use this collapse timescale to predict variations in the molecular gas star formation efficiency, which is lowered from a maximum, feedback-regulated level in the presence of strong coupling to the galactic potential. Our model implies that gas can only decouple from the galaxy to collapse and efficiently form stars deep within clouds. We show that this naturally explains the observed drop in star formation rate per unit gas mass in the Milky Way's CMZ and other galaxy centers. The model for a galactic bottleneck to star formation also agrees well with resolved observations of dense gas and star formation in galaxy disks and the properties of local clouds., Comment: Accepted for publication in ApJ, 32 pages, 11 figures

Square Kilometre Array: et al., arXiv:1810.02680v3, The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy, and spacetime., A. Racanelli has received funding from the People Programme (Marie Curie Actions) of the European Union H2020 Programme under REA grant agreement number 706896 (COSMOFLAGS). Funding for this work was partially provided by the Spanish MINECO under MDM-2014-0369 of ICCUB (Unidad de Excelencia ‘Maria de Maeztu’).

We present the second dust continuum data release in the Census of High- and Medium-mass Protostars (CHaMP), expanding the methodology trialed in Pitts et al. 2019 to the entire CHaMP survey area ($280^{\circ}, Comment: 33 pages, 13 figures, 2 tables, 51 pages of appendices; accepted to ApJS (in print)

Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_{\phi}\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers., Comment: 6 pages, 5 figures, submitted to A&A

In our HST photometric survey, we have been searching for multiple stellar populations (MPs) in Magellanic Clouds (MCs) massive star clusters which span a significant range of ages ($\sim 1.5-11$ Gyr). In the previous papers of the series, we have shown that the age of the cluster represents one of the key factors in shaping the origin of the chemical anomalies. Here we present the analysis of four additional clusters in the MCs, namely Lindsay 38, Lindsay 113, NGC 2121 and NGC 2155, for which we recently obtained new UV HST observations. These clusters are more massive than $\sim 10^4M_{\odot}$ and have ages between $\sim 2.5-6$ Gyr, i.e. located in a previously unexplored region of the cluster age/mass diagram. We found chemical anomalies, in the form of N spreads, in three out of four clusters in the sample, namely in NGC 2121, NGC 2155 and Lindsay 113. By combining data from our survey and HST photometry for 3 additional clusters in the Milky Way (namely 47 Tuc, M15 and NGC 2419), we show that the extent of the MPs in the form of N spread is a strong function of age, with older clusters having larger N spreads with respect to the younger ones. Hence, we confirm that cluster age plays a significant role in the onset of MPs., Accepted for publication in MNRAS. 12 pages, 9 figures

We utilise the final catalogue from the Pan-Andromeda Archaeological Survey to investigate the links between the globular cluster system and field halo in M31 at projected radii $R_p=25-150$ kpc. In this region the cluster radial density profile exhibits a power-law decline with index $\Gamma=-2.37\pm0.17$, matching that for the stellar halo component with [Fe/H] $, Comment: Accepted for publication in MNRAS; this is the accepted version

In the previous papers in this series, we have measured the stellar and \hi content in a sample of edge-on galaxies. In the present paper, we perform a simultaneous rotation curve and vertical force field gradient decomposition for five of these edge-on galaxies. The rotation curve decomposition provides a measure of the radial dark matter potential, while the vertical force field gradient provide a measure of the vertical dark matter potential. We fit dark matter halo models to these potentials. Using our \hi self-absorption results, we find that a typical dark matter halo has a less dense core ($0.094\pm0.230$\,M$_\odot$/pc$^3$) compared to an optically thin \hi model ($0.150\pm0.124$\,M$_\odot$/pc$^3$). The HI self-absorption dark matter halo has a longer scale length $R_c$ of $1.42\pm 3.48$\,kpc, versus $1.10\pm 1.81$\,kpc for the optically thin HI model. The median halo shape is spherical, at $q=1.0\pm0.6$ (self-absorbing \hi), while it is prolate at $q=1.5\pm0.6$ for the optically thin. Our best results were obtained for ESO\,274-G001 and UGC\,7321, for which we were able to measure the velocity dispersion in Paper III. These two galaxies have drastically different halo shapes, with one oblate and one strongly prolate. Overall, we find that the many assumptions required make this type of analysis susceptible to errors., Comment: Accepted for publication by MNRAS. Earlier papers in the series are in the refereeing process. High-res. version of this paper is available at www.astro.rug.nl/~vdkruit/Petersetal-V.pdf

We present deep near-infrared photometry and spectroscopy of the globular cluster 2MASS-GC03 projected in the Galactic disk using MMIRS on the Clay telescope (Las Campanas Observatory) and VISTA Variables in the Via Lactea survey (VVV) data. Most probable cluster member candidates were identified from near-infrared photometry. Out of ten candidates that were followed-up spectroscopically, five have properties of cluster members, from which we calculate = -0.9 +- 0.2 and a radial velocity of v_r > = -78 +- 12km/s. A distance of 10.8kpc is estimated from 3 likely RRLyrae members. Given that the cluster is currently at a distance of 4.2kpc from the Galactic center, the cluster's long survival time of an estimated 11.3 +- 1.2Gyr strengthens the case for its globular-cluster nature. The cluster has a hint of elongation in the direction of the Galactic center., Comment: Accepted for publication, MNRAS

We present SOFIA/FIFI-LS observations of the [CII] 158${\mu}$m cooling line across the nearby spiral galaxy NGC 6946. We combine these with UV, IR, CO, and H I data to compare [CII] emission to dust properties, star formation rate (SFR), H$_2$, and HI at 560pc scales via stacking by environment (spiral arms, interarm, and center), radial profiles, and individual, beam-sized measurements. We attribute $73\%$ of the [CII] luminosity to arms, and $19\%$ and $8\%$ to the center and interarm region, respectively. [CII]/TIR, [CII]/CO, and [CII]/PAH radial profiles are largely constant, but rise at large radii ($\gtrsim$8kpc) and drop in the center ("[CII] deficit"). This increase at large radii and the observed decline with the 70${\mu}$m/100${\mu}$m dust color are likely driven by radiation field hardness. We find a near proportional [CII]-SFR scaling relation for beam-sized regions, though the exact scaling depends on methodology. [CII] also becomes increasingly luminous relative to CO at low SFR (interarm or large radii), likely indicating more efficient photodissociation of CO and emphasizing the importance of [CII] as an H$_2$ and SFR tracer in such regimes. Finally, based on the observed [CII] and CO radial profiles and different models, we find ${\alpha}_{CO}$ to increase with radius, in line with the observed metallicity gradient. The low ${\alpha}_{CO}$ (galaxy average $\lesssim2\,M_{sun}\,pc^{-2}\,(K\,km\,s^{-1})^{-1}$) and low [CII]/CO ratios ($\sim$400 on average) imply little CO-dark gas across NGC 6946, in contrast to estimates in the Milky Way., Comment: 12 pages, 6 Figures, published in ApJ

The search for technosignatures from hypothetical galactic civilizations is going through a new phase of intense activity. For the first time, a significant fraction of the vast search space is expected to be sampled in the foreseeable future, potentially bringing informative data about the abundance of detectable extraterrestrial civilizations, or the lack thereof. Starting from the current state of ignorance about the galactic population of non-natural electromagnetic signals, we formulate a Bayesian statistical model to infer the mean number of radio signals crossing Earth, assuming either non-detection or the detection of signals in future surveys of the Galaxy. Under fairly noninformative priors, we find that not detecting signals within about $1$ kly from Earth, while suggesting the lack of galactic emitters or at best the scarcity thereof, is nonetheless still consistent with a probability exceeding $10$ \% that typically over $\sim 100$ signals could be crossing Earth, with radiated power analogous to that of the Arecibo radar, but coming from farther in the Milky Way. The existence in the Galaxy of potentially detectable Arecibo-like emitters can be reasonably ruled out only if all-sky surveys detect no such signals up to a radius of about $40$ kly, an endeavor requiring detector sensitivities thousands times higher than those of current telescopes. Conversely, finding even one Arecibo-like signal within $\sim 1000$ light years, a possibility within reach of current detectors, implies almost certainly that typically more than $\sim 100$ signals of comparable radiated power cross the Earth, yet to be discovered., Published in PNAS ahead of print October 1, 2018. Preprint has 13 pages, 7 figures + 7 pages of Supplementary Information with 5 figures

A number of Local Group dwarf galaxies are known to have two spatially segregated stellar metallicity populations, a centrally concentrated metal-rich population and a more extended metal-poor population. In this work we discuss mechanisms that lead to the formation of two spatially segregated metallicity populations. Using a set of high-resolution hydrodynamical simulations of Local Group-like environments, we select a sample of satellite and field galaxies, spanning the stellar mass range $10^6-10^9$M$_{\odot}$, that exhibit bimodality in their metallicity distributions. Among those, we identify a subsample with a strong spatial segregation in the two populations. We find three distinct mechanisms for their formation. In field dwarfs and in a small fraction of satellites, a merger causes the metal-poor stars to migrate to larger radii and encourages the available gas to sink to the centre of the dwarf. Most of the gas is subsequently blown out of the halo through star formation feedback, but the remaining gas is consumed in the formation of a metal-rich population. In the exclusive case of satellites that have retained some of their gas at infall, it is the compression of this gas by ram pressure near pericentre that triggers the formation of metal-rich stars, whilst simultaneously preventing star formation at larger radii through stripping. Additionally, in a small number of field and satellite dwarfs, interactions with gaseous filaments and other galaxies can result in the formation of a metal-rich population. Regardless of the formation mechanism, a history of mergers typically enhances the spatial segregation., Comment: 20 pages, 19 figures, accepted in MNRAS

The azimuthal variation of the HII region oxygen abundance in spiral galaxies is a key observable for understanding how quickly oxygen produced by massive stars can be dispersed within the surrounding interstellar medium. Observational constraints on the prevalence and magnitude of such azimuthal variations remain rare in the literature. Here, we report the discovery of pronounced azimuthal variations of HII region oxygen abundance in NGC 2997, a spiral galaxy at approximately 11.3 Mpc. Using 3D spectroscopic data from the TYPHOON Program, we study the HII region oxygen abundance at a physical resolution of 125 pc. Individual HII regions or complexes are identified in the 3D optical data and their strong emission line fluxes measured to constrain their oxygen abundances. We find 0.06 dex azimuthal variations in the oxygen abundance on top of a radial abundance gradient that is comparable to those seen in other star-forming disks. At a given radial distance, the oxygen abundances are highest in the spiral arms and lower in the inter-arm regions, similar to what has been reported in NGC 1365 using similar observations. We discuss whether the azimuthal variations could be recovered when the galaxy is observed at worse physical resolutions and lower signal-to-noise ratios., Comment: 13 pages, 17 figures, accepted to A&A

Context. The great majority of early-type dwarf galaxies, in the Local Group as well as in other galaxy groups, are found in the vicinity of much larger galaxies, making it hard to disentangle the role of internal versus external effects in driving their evolution., Aims. In order to minimize environmental effects and gain an insight into the internal mechanisms that shape the properties of these systems, we study one of the few dwarf spheroidal galaxies of the Local Group found in isolation: Cetus. This system is of particular interest since it does not follow the Local Group morphology-density relation., Methods. We obtained Very Large Telescope (VLT) FORS2 spectra (R similar to 2600) in the region of the nIR CaII triplet lines for 80 candidate red giant branch (RGB) stars. The analysis yielded line-of-sight velocities and metallicities ([Fe/H]) for 54 bona fide member stars. Results. The kinematic analysis shows that Cetus is a mainly pressure-supported (sigma(v) = 11.0(-1.3)(+1.6) km s(-1)), dark-matter-dominated system (M-1/2/L-V = 23.9(-8.9)(+9.)7 M-circle dot/L-circle dot) with no significant signs of internal rotation. We find Cetus to be a metal-poor system with a significant metallicity spread (median [Fe/H] = -1.71 dex, median-absolute-deviation = 0.49 dex), as expected for its stellar mass. We report the presence of a mild metallicity gradient compatible with those found in other dwarf spheroidals of the same luminosity; we trace the presence of a stellar population gradient also in the spatial distribution of stars in different evolutionary phases in ancillary SuprimeCam photometric data. There are tentative indications of two chemo-kinematically distinct sub-populations, with the more metal-poor stars showing a hotter kinematics than the metal-richer ones. Furthermore, the photometric dataset reveals the presence of a foreground population that most likely belongs to the Sagittarius stream., Conclusions. This study represents an important step forward in assessing the internal kinematics of the Cetus dwarf spheroidal galaxy as well as the first wide-area spectroscopic determination of its metallicity properties. With our analysis, Cetus adds to the growing scatter in stellar-dark matter halo properties in low-mass galactic systems. The presence of a metallicity gradient akin to those found in similarly luminous and pressure-supported systems inhabiting very different environments may hint at metallicity gradients in Local Group early-type dwarfs being driven by internal mechanisms.

The European physical journal / C 77(12), 831 (2017). doi:10.1140/epjc/s10052-017-5155-4, We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages.This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamLike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDCalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs,and the neutrino telescope likelihood package nulike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for resultsobtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application ofDarkBit’s advanced direct and indirect detection routines to a simple effective dark matter model., Published by Springer, Berlin

We examine the spatial distribution of the oldest and most metal poor stellar populations of Milky Way-sized galaxies using the APOSTLE cosmological hydrodynamical simulations of the Local Group. In agreement with earlier work, we find strong radial gradients in the fraction of the oldest (tform < 0.8 Gyr) and most metal poor ([Fe/H]< -2.5) stars, both of which increase outwards. The most metal poor stars form over an extended period of time; half of them form after z = 5.3, and the last 10% after z = 2.8. The age of the metal poor stellar population also shows significant variation with environment; a high fraction of them are old in the galaxy's central regions and an even higher fraction in some individual dwarf galaxies, with substantial scatter from dwarf to dwarf. Overall, over half of the stars that belong to both the oldest and most metal-poor population are found outside the solar circle. Somewhat counter-intuitively, we find that dwarf galaxies with a large fraction of metal poor stars that are very old are systems where metal poor stars are relatively rare, but where a substantial old population is present. Our results provide guidance for interpreting the results of surveys designed to hunt for the earliest and most pristine stellar component of our Milky Way., 14 pages, 9 figures. Submitted to MNRAS, revised after referee comments

We determine relative abundance ratios for the neutron-capture elements Zr, La, Ce, Nd, and Eu for a sample of 27 Galactic dwarf stars with -1.5 < [Fe/H], Comment: MNRAS in press

The forcing of interstellar turbulence, driven mainly by supernova explosions, is irrotational in nature, but the development of significant amounts of vorticity and helicity, accompanied by large-scale dynamo action, has been reported. Several earlier investigations examined vorticity production in simpler systems; here all the relevant processes can be considered simultaneously. We also investigate the mechanisms for the generation of net helicity and large-scale flow in the system. We use a three-dimensional, stratified, rotating and shearing local simulation domain of the size 1x1x2 kpc$^3$, forced with SN explosions occurring at the rate typical of the solar neighbourhood in the Milky Way. In addition to the nominal simulation run with realistic Milky Way parameters, we vary the rotation and shear rates, but keep the absolute value of their ratio fixed. Reversing the sign of shear vs. rotation allows us to separate the rotation- and shear-generated contributions. As in earlier studies, we find the generation of significant amounts of vorticity, with on average 65% of the kinetic energy being in the rotational modes. The vorticity production can be related to the baroclinicity of the flow, especially in the regions of hot, dilute clustered supernova bubbles. In these regions, the vortex stretching acts as a sink of vorticity. The net helicities produced by rotation and shear are of opposite signs for physically motivated rotation laws, with the solar neighbourhood parameters resulting in the near cancellation of the total net helicity. We also find the excitation of oscillatory mean flows, the strength and oscillation period of which depend on the Coriolis and shear parameters; we interpret these as signatures of the anisotropic kinetic (AKA) effect. We use the method of moments to fit for the turbulent transport coeffcients, and find $\alpha_{\rm AKA}$ values of the order 3-5 km/s., Comment: 19 pages, 12 figures, Astronomy and Astrophysics, accepted for publication

The origin of dark matter and the generation of neutrino masses could be related if neutrino masses arise from the spontaneous violation of ungauged lepton number. In this case the associated Nambu–Goldstone boson, the majoron, could acquire a mass from non-perturbative gravitational effects and play the role of DM. Here we report our cosmological and astrophysical constraints on majoron dark matter coming from Cosmic Microwave Background (CMB) and a variety of X- and γ -ray observations.

We present maps of the dust properties in the Small and Large Magellanic Clouds (SMC, LMC) from fitting Spitzer and Herschel observations with the \citet{DL07} dust model. We derive the abundance of the small carbonaceous grain (or polycyclic aromatic hydrocarbon; PAH) component. The global PAH fraction (q_pah, the fraction of the dust mass in the form of PAHs) is smaller in the SMC (1.0$^{+0.3}_{-0.3}$%) than in the LMC (3.3$^{+1.4}_{-1.3}$%). We measure the PAH fraction in different gas phases (H II regions, ionized gas outside of H II regions, molecular gas, and diffuse neutral gas). H II regions appear as distinctive holes in the spatial distribution of the PAH fraction. In both galaxies, the PAH fraction in the diffuse neutral medium is higher than in the ionized gas, but similar to the molecular gas. Even at equal radiation field intensity, the PAH fraction is lower in the ionized gas than in the diffuse neutral gas. We investigate the PAH life-cycle as a function of metallicity between the two galaxies. The PAH fraction in the diffuse neutral medium of the LMC is similar to that of the Milky Way ($\sim4.6$%), while it is significantly lower in the SMC. Plausible explanations for the higher PAH fraction in the diffuse neutral medium of the LMC compared to the SMC include: a more effective PAH production by fragmentation of large grains at higher metallicity, and/or the growth of PAHs in molecular gas., Comment: Accepted for publication in ApJ

We use our integrated SDSS photometry for 96 globular clusters in $g$ and $z$, as well as $r$ and $i$ photometry for a subset of 56 clusters, to derive the integrated colour-metallicity relation (CMR) for Galactic globular clusters. We compare this relation to previous work, including extragalactic clusters, and examine the influence of age, present-day mass function variations, structural parameters and the morphology of the horizontal branch on the relation. Moreover, we scrutinise the scatter introduced by foreground extinction (including differential reddening) and show that the scatter in the colour-metallicity relation can be significantly reduced combining two reddening laws from the literature. In all CMRs we find some low-reddening young GCs that are offset to the CMR. Most of these outliers are associated with the Sagittarius system. Simulations show that this is due less to age than to a different enrichment history. Finally, we introduce colour-metallicity relations based on the infrared Calcium triplet, which are clearly non-linear when compared to $(g^\prime-i^\prime)$ and $(g^\prime-z^\prime)$ colours., Accepted for publication in MNRAS. 17 pages, 18 figures

We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in comparison to the more commonly used Jeans modelling, and it allows us to derive the dark matter content in a robust way. Here we present our results for the Sculptor dwarf spheroidal galaxy, after testing our methods on mock data sets. We fit both the second and fourth velocity moment profile to break the mass-anisotropy degeneracy. For an NFW dark matter halo profile, we find that the mass of Sculptor within 1 kpc is M_1kpc = 1.03 \pm 0.07 x 10^8 Msol, and that its velocity anisotropy profile is tangentially biased and nearly constant with radius. The preferred concentration (c ~ 15) is low for its dark matter mass but consistent within the scatter found in N-body cosmological simulations. When we let the value of the central logarithmic slope \alpha{} vary, we find that the best-fit model has \alpha{} = 0, although an NFW cusp or shallower is consistent at 1 \sigma{} confidence level. On the other hand, very cuspy density profiles with logarithmic central slopes \alpha{} < -1.5 are strongly disfavoured for Sculptor., Comment: 18 pages, 12 figures, accepted by MNRAS http://mnras.oxfordjournals.org/content/433/4/3173

We analyze the distribution of G and K type stars towards the Galactic poles using RAVE and ELODIE radial velocities, 2MASS photometric star counts, and UCAC2 proper motions. The combination of photometric and 3D kinematic data allows us to disentangle and describe the vertical distribution of dwarfs, sub-giants and giants and their kinematics. We identify discontinuities within the kinematics and magnitude counts that separate the thin disk, thick disk and a hotter component. The respective scale heights of the thin disk and thick disk are 225$\pm$10 pc and 1048$\pm$36 pc. We also constrain the luminosity function and the kinematic distribution function. The existence of a kinematic gap between the thin and thick disks is incompatible with the thick disk having formed from the thin disk by a continuous process, such as scattering of stars by spiral arms or molecular clouds. Other mechanisms of formation of the thick disk such as `created on the spot' or smoothly `accreted' remain compatible with our findings., Comment: 15 pages, 13 EPS figures and 1 table

We present a new discrete chemo-dynamical axisymmetric modeling technique, which we apply to the dwarf spheroidal galaxy Sculptor. The major improvement over previous Jeans models is that realistic chemical distributions are included directly in the dynamical modelling of the discrete data. This avoids loss of information due to spatial binning and eliminates the need for hard cuts to remove contaminants and to separate stars based on their chemical properties. Using a combined likelihood in position, metallicity and kinematics, we find that our models naturally separate Sculptor stars into a metal-rich and a metal-poor population. Allowing for non-spherical symmetry, our approach provides a central slope of the dark matter density of $\gamma = 0.5 \pm 0.3$. The metal-rich population is nearly isotropic (with $\beta_r^{red} = 0.0\pm0.1$) while the metal-poor population is tangentially anisotropic (with $\beta_r^{blue} = -0.2\pm0.1$) around the half light radius of $0.26$ kpc. A weak internal rotation of the metal-rich population is revealed with $v_{max}/\sigma_0 = 0.15 \pm 0.15$. We run tests using mock data to show that a discrete dataset with $\sim 6000$ stars is required to distinguish between a core ($\gamma = 0$) and cusp ($\gamma = 1$), and to constrain the possible internal rotation to better than $1\,\sigma$ confidence with our model. We conclude that our discrete chemo-dynamical modelling technique provides a flexible and powerful tool to robustly constrain the internal dynamics of multiple populations, and the total mass distribution in a stellar system., Comment: Accepted by MNRAS

The stellar metallicity gradients of disc galaxies provide information on the disc assembly, star formation processes and chemical evolution. They also might store information on dynamical processes which could affect the distribution of chemical elements in the gas-phase and the stellar components. We studied the stellar metallicity gradients of stellar discs in a cosmological simulation. We explored the dependence of the stellar metallicity gradients on stellar age and the size and mass of the stellar discs. We used galaxies selected from a cosmological hydrodynamical simulation performed including a physically-motivated Supernova feedback and chemical evolution. The metallicity profiles were estimated for stars with different ages. We confront our numerical findings with results from the CALIFA Survey. The simulated stellar discs are found to have metallicity profiles with slopes in global agreement with observations. Low stellar-mass galaxies tend to have a larger variety of metallicity slopes. When normalized by the half-mass radius, the stellar metallicity gradients do not show any dependence and the dispersion increases significantly, regardless of galaxy mass. Galaxies with stellar masses around $10^{10}$M$_{\odot}$ show steeper negative metallicity gradients. The stellar metallicity gradients correlate with the half-mass radius. However, the correlation signal is not present when they are normalized by the half-mass radius. Stellar discs with positive age gradients are detected to have negative and positive metallicity gradients, depending on the relative importance of the recent star formation activity in the central regions. The large dispersions in the metallicity gradients as a function of stellar mass could be ascribed to the effects of dynamical processes such as mergers/interactions and/or migration as well as those regulating the conversion of gas into stars. [abridged], Comment: 9 pages, 9 figures, accepted for publication in A&A

We combine simultaneous constraints on stellar evolutionary status from asteroseismology, and on nitrogen abundances derived from large spectroscopic surveys, to follow nitrogen surface abundances all along the evolution of a low-mass star, comparing model expectations with data. After testing and calibrating the observed yields from the APOGEE survey, we first show that nitrogen surface abundances follow the expected trend after the first dredge-up occurred, i.e. that the more massive is the star the more nitrogen is enhanced. Moreover, the behaviour of nitrogen data along the evolution confirms the existence of non-canonical extra-mixing on the RGB for all low-mass stars in the field. But more surprisingly, the data indicate that nitrogen has been depleted between the RGB tip and the red clump. This may suggest that some nitrogen has been burnt near or at the He flash episode., Comment: 9 pages, 11 figures, accepted for publication in MNRAS