Search

We examine the distribution of radio emission from ~42,000 quasars from the Sloan Digital Sky Survey, as measured in the LOFAR Two-Metre Sky Survey (LoTSS). We present a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: active galactic nuclei (jets) and star-formation. Our two-component model provides an excellent match to the observed radio flux density distributions across a wide range of redshifts and quasar optical luminosities; this suggests that the jet-launching mechanism operates in all quasars but with different powering efficiency. The wide distribution of jet powers allows for a smooth transition between the 'radio-quiet' and 'radio-loud' quasar regimes, without need for any explicit bimodality. The best-fit model parameters indicate that the star-formation rate of quasar host galaxies correlates strongly with quasar luminosity and also increases with redshift at least out to z~2. For a model where star-formation rate scales as $SFR \propto L_{bol}^\alpha (1+z)^\beta$, we find $\alpha = 0.47 \pm 0.01$ and $\beta = 1.61 \pm 0.05$, in agreement with far-infrared studies. Quasars contribute ~0.15 per cent of the cosmic star-formation rate density at z=0.5, rising to 0.4 per cent by z=2. The typical radio jet power is seen to increase with both increasing optical luminosity and black hole mass independently, but does not vary with redshift, suggesting intrinsic properties govern the production of the radio jets. We discuss the implications of these results for the triggering of quasar activity and the launching of jets., Comment: Accepted for publication in MNRAS. 21 pages

We study the radio spectral properties of 2,094 star-forming galaxies (SFGs) by combining our early science data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey with VLA, GMRT radio data, and rich ancillary data in the COSMOS field. These SFGs are selected at VLA 3GHz, and their flux densities from MeerKAT 1.3GHz and GMRT 325MHz imaging data are extracted using the "super-deblending" technique. The median radio spectral index is $\alpha_{\rm 1.3GHz}^{\rm 3GHz}=-0.80\pm0.01$ without significant variation across the rest-frame frequencies ~1.3-10GHz, indicating radio spectra dominated by synchrotron radiation. On average, the radio spectrum at observer-frame 1.3-3GHz slightly steepens with increasing stellar mass with a linear fitted slope of $\beta=-0.08\pm0.01$, which could be explained by age-related synchrotron losses. Due to the sensitivity of GMRT 325MHz data, we apply a further flux density cut at 3GHz ($S_{\rm 3GHz}\ge50\,\mu$Jy) and obtain a sample of 166 SFGs with measured flux densities at 325MHz, 1.3GHz, and 3GHz. On average, the radio spectrum of SFGs flattens at low frequency with the median spectral indices of $\alpha^{\rm 1.3GHz}_{\rm 325MHz}=-0.59^{+0.02}_{-0.03}$ and $\alpha^{\rm 3.0GHz}_{\rm 1.3GHz}=-0.74^{+0.01}_{-0.02}$. At low frequency, our stacking analyses show that the radio spectrum also slightly steepens with increasing stellar mass. By comparing the far-infrared-radio correlations of SFGs based on different radio spectral indices, we find that adopting $\alpha_{\rm 1.3GHz}^{\rm 3GHz}$ for $k$-corrections will significantly underestimate the infrared-to-radio luminosity ratio ($q_{\rm IR}$) for >17% of the SFGs with measured flux density at the three radio frequencies in our sample, because their radio spectra are significantly flatter at low frequency (0.33-1.3GHz)., Comment: 15 pages, 6 figures, 1 table; accepted for publication in MNRAS on 2021 Aug 4

We present results from the NIRVANDELS survey investigating the gas-phase metallicity ($\mathrm{Z}_{\mathrm{gas}}$, tracing O/H) and stellar metallicity ($Z_{\star}$, tracing Fe/H) of 33 star-forming galaxies at redshifts $2.95 < z < 3.80$. Based on a combined analysis of deep optical and near-IR spectra, tracing the rest-frame far ultraviolet and rest-frame optical respectively, we present the first simultaneous determination of the stellar and gas-phase mass-metallicity relationships (MZRs) at $z\simeq3.4$. In both cases, we find that metallicity increases with increasing stellar mass ($M_{\star}$), and that the power-law slope at $M_{\star} \lesssim 10^{10} \mathrm{M}_{\odot}$ of both MZRs scales as $Z \propto M_{\star}^{0.3}$. Comparing the stellar and gas-phase MZRs, we present direct evidence for super-solar O/Fe ratios (i.e., $\alpha$-enhancement) at $z>3$, finding $\mathrm{(O/Fe)}\simeq (2.54 \pm 0.38) \times \mathrm{(O/Fe)}_{\odot}$, with no clear dependence on $M_{\star}$., Comment: 18 pages, 11 figures, MNRAS accepted

The search for extrasolar planets has driven rapid advances in instrumentation, resulting in cameras such as SPHERE at the VLT, GPI at Gemini South and SCExAO at Subaru, capable of achieving very high contrast ($\sim10^{6}$) around bright stars with small inner working angles ($\sim 0\farcs{1}$). The optimal exploitation of data from these instruments depends on the availability of easy-to-use software to process and analyse their data products. We present a pure-python pipeline, {\sc precision}, which provides fast, memory-efficient reduction of data from the SPHERE/IRDIS near-infrared imager, and can be readily extended to other instruments. We apply {\sc precision} to observations of the extreme red supergiant VX~Sgr, the inner outflow of which is revealed to host complex, asymmetric structure in the near-IR. In addition, optical polarimetric imaging reveals clear extended polarised emission on $\sim0.5^{\prime\prime}$ scales which varies significantly with azimuth, confirming the asymmetry. While not conclusive, this could suggest that the ejecta are confined to a disc or torus, which we are viewing nearly face on, although other non-spherical or clumpy configurations remain possible. VX~Sgr has no known companions, making such a geometry difficult to explain, as there is no obvious source of angular momentum in the system., 13 pages, 5 figures. Accepted by MNRAS

We describe the observation and confirmation of nine new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on (i) galaxy group and cluster samples, and (ii) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope to acquire a spectroscopic follow-up of 21 of these candidates. Through an analysis of these spectroscopic follow-up data, we confirmed nine new lensing systems and rejected two candidates, and the analysis was inconclusive on 10 candidates. For each of the confirmed systems, our report measured spectroscopic properties, estimated source image-lens separations, and estimated enclosed masses as well. The sources that we targeted have an i-band surface brightness range of i ∼ 22-24 mag arcsec and a spectroscopic redshift range of z ∼ 0.8−2.6. The lens galaxies have a photometric redshift range of z ∼ 0.3−0.7. The lensing systems range in source image-lens separation from 2 to 9 arcsec and in enclosed mass from 10 to 10 M, Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundac¸ão Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência e Tecnologia, the Deutsche Forschungsge-meinschaft, and the Collaborating Institutions in the DES. The DES data management system is supported by the National Science Foundation under Grant Number AST-1138766. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2012-39559, ESP2013-48274, FPA2013-47986, and Centro de Excelencia Severo Ochoa SEV-2012-0234, some of which include ERDF funds from the European Union. This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics.

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall)., Comment: 16 pages, 15 figures, submitted to MNRAS

We investigate the effect of microlensing on the standardisation of strongly lensed Type Ia supernovae (GLSNe Ia). We present predictions for the amount of scatter induced by microlensing across a range of plausible strong lens macromodels. We find that lensed images in regions of low convergence, shear and stellar density are standardisable, where the microlensing scatter is < 0.15 magnitudes, comparable to the intrinsic dispersion of for a typical SN Ia. These standardisable configurations correspond to asymmetric lenses with an image located far outside the Einstein radius of the lens. Symmetric and small Einstein radius lenses (< 0.5 arcsec) are not standardisable. We apply our model to the recently discovered GLSN Ia iPTF16geu and find that the large discrepancy between the observed flux and the macromodel predictions from More et al. (2017) cannot be explained by microlensing alone. Using the mock GLSNe Ia catalogue of Goldstein et al. (2017), we predict that ~ 22% of GLSNe Ia discovered by LSST will be standardisable, with a median Einstein radius of 0.9 arcseconds and a median time-delay of 41 days. By breaking the mass-sheet degeneracy the full LSST GLSNe Ia sample will be able to detect systematics in H0 at the 0.5% level., 11 pages, 8 Figures. Accepted by MNRAS May 17 2018

We present a new method that incorporates the horizontal branch morphology into synthetic colour-magnitude diagram based star formation history determinations. This method, we call MORGOTH, self-consistently takes into account all the stellar evolution phases up to the early asymptothic giant branch, flexibly modelling red giant branch mass loss. We test MORGOTH on a range of synthetic populations, and find that the inclusion of the horizontal branch significantly increases the precision of the resulting star formation histories. When the main sequence turn-off is detected, MORGOTH can fit the star formation history and the red giant branch mass loss at the same time, efficiently breaking this degeneracy. As part of testing MORGOTH, we also model the observed colour-magnitude diagram of the well studied Sculptor dwarf spheroidal galaxy. We recover a new more detailed star formation history for this galaxy. Both the new star formation history and the red giant branch mass loss we determined for Sculptor with MORGOTH are in good agreement with previous analyses, thus demonstrating the power of this new approach., 12 pages, 11 figures, accepted for publication in MNRAS

We study the environments of low- and high- excitation radio galaxies (LERGs and HERGs respectively) in the redshift range $0.01 < z < 0.4$, using a sample of 399 radio galaxies and environmental measurements from the Galaxy And Mass Assembly (GAMA) survey. In our analysis we use the fifth nearest neighbour density ($\Sigma_{5}$) and the GAMA galaxy groups catalogue (G3Cv6) and construct control samples of galaxies matched in {\update stellar mass and colour} to the radio-detected sample. We find that LERGs and HERGs exist in different environments and that this difference is dependent on radio luminosity. High-luminosity LERGs ($L_{\rm NVSS} \gtrsim 10^{24}$ W Hz$^{-1}$) lie in much denser environments than a matched radio-quiet control sample (about three times as dense, as measured by $\Sigma_{5}$), and are more likely to be members of galaxy groups ($82^{+5}_{-7}$ percent of LERGs are in GAMA groups, compared to $58^{+3}_{-3}$ percent of the control sample). In contrast, the environments of the HERGs and lower luminosity LERGs are indistinguishable from that of a matched control sample. Our results imply that high-luminosity LERGs lie in more massive haloes than non-radio galaxies of similar stellar mass and colour, in agreement with earlier studies (Wake et al. 2008; Donoso et al. 2010). When we control for the preference of LERGs to be found in groups, both high- and low- luminosity LERGs are found in higher-mass haloes ($\sim 0.2$ dex; at least 97 percent significant) than the non-radio control sample., Comment: 18 pages, 11 figures. Accepted for publication by MNRAS

The relationships between stellar mass, gas-phase metallicity and star formation rate (i.e. the Mass-Metallicity, MZR, and the Fundamental Metallcity Relation, FMR) in the local Universe are revisited by fully anchoring the metallicity determination for SDSS galaxies on the Te abundance scale defined exploiting the strong-line metallicity calibrations presented in Curti et al. (2017). Self-consistent metallicity measurements allow a more unbiased assessment of the scaling relations involving M, Z and SFR, which provide powerful constraints for the chemical evolution models. We parametrise the MZR with a new functional form which allows us to better characterise the turnover mass. The slope and saturation metallicity are in good agreement with previous determinations of the MZR based on the Te method, while showing significantly lower normalisation compared to those based on photoionisation models. The Z- SFR dependence at fixed stellar mass is also investigated, being particularly evident for highly star forming galaxies, where the scatter in metallicity is reduced up to a factor of ~30%. A new parametrisation of the FMR is given by explicitly introducing the SFR-dependence of the turnover mass into the MZR. The residual scatter in metallicity for the global galaxy population around the new FMR is 0.054 dex. The new FMR presented in this work represents a useful local benchmark to compare theoretical predictions and observational studies (of both local and high-redshift galaxies) whose metallicity measurements are tied to the abundance scale defined by the Te method, hence allowing to properly assess its evolution with cosmic time., Comment: 21 pages, 13 figures. Accepted for publication on MNRAS. V2: Equation (3) corrected

A new family of self-consistent DF-based models of stellar systems is explored. The stellar component of the models is described by a distribution function (DF) depending on the action integrals, previously used to model the Fornax dwarf spheroidal galaxy (dSph). The stellar component may cohabit with either a dark halo, also described by a DF, or with a massive central black hole. In all cases we solve for the model's self-consistent potential. Focussing on spherically symmetric models, we show how the stellar observables vary with the anisotropy prescribed by the DF, with the dominance and nature of the dark halo, and with the mass of the black hole. We show that precise fits to the observed surface brightness profiles of four globular clusters can be obtained for a wide range of prescribed velocity anisotropies. We also obtain precise fits to the observed projected densities of four dSphs. Finally, we present a three-component model of the Scupltor dSph with distinct DFs for the red and blue horizontal branch stars and the dark matter halo., submitted to MNRAS

We use classifications provided by citizen scientists though Galaxy Zoo to investigate the correlation between bulge size and arm winding in spiral galaxies. Whilst the traditional spiral sequence is based on a combination of both measures, and is supposed to favour arm winding where disagreement exists, we demonstrate that, in modern usage, the spiral classifications Sa-Sd are predominantly based on bulge size, with no reference to spiral arms. Furthermore, in a volume limited sample of galaxies with both automated and visual measures of bulge prominence and spiral arm tightness, there is at best a weak correlation between the two. Galaxies with small bulges have a wide range of arm winding, while those with larger bulges favour tighter arms. This observation, interpreted as revealing a variable winding speed as a function of bulge size, may be providing evidence that the majority of spiral arms are not static density waves, but rather wind-up over time. This suggests the "winding problem" could be solved by the constant reforming of spiral arms, rather than needing a static density wave. We further observe that galaxies exhibiting strong bars tend have more loosely wound arms at a given bulge size than unbarred spirals. This observations suggests that the presence of a bar may slow the winding speed of spirals, and may also drive other processes (such as density waves) which generate spiral arms. It is remarkable that after over 170 years of observations of spiral arms in galaxies our understanding of them remains incomplete., Comment: 16 pages, 8 figures, accepted for publication in MNRAS

Using deep photometric data from CFHT/Megacam, we study the morphology and density of the GD-1 stream, one of the longest and coldest stellar streams in the Milky Way. Our deep data recovers the lower main sequence of the stream with unprecedented quality, clearly separating it from Milky Way foreground and background stars. An analysis of the distance to different parts of the stream shows that GD-1 lies at a heliocentric distance between 8 and 10 kpc, with only a shallow gradient across 45 deg on the sky. Matched filter maps of the stream density show clear density variations, such as deviations from a single orbital track and tentative evidence for stream fanning. We also detect a clear under-density in the middle of the stream track at $\varphi_{1}$=-45 deg surrounded by overdense stream segments on either side. This location is a promising candidate for the elusive missing progenitor of the GD-1 stream. We conclude that the GD-1 stream has clearly been disturbed by interactions with the Milky Way disk or other sub-halos., Comment: submitted to MNRAS, comments welcome

Faint tidal features around galaxies record their merger and interaction histories over cosmic time. Due to their low surface brightnesses and complex morphologies, existing automated methods struggle to detect such features and most work to date has heavily relied on visual inspection. This presents a major obstacle to quantitative study of tidal debris features in large statistical samples, and hence the ability to be able to use these features to advance understanding of the galaxy population as a whole. This paper uses convolutional neural networks (CNNs) with dropout and augmentation to identify galaxies in the CFHTLS-Wide Survey that have faint tidal features. Evaluating the performance of the CNNs against previously-published expert visual classifications, we find that our method achieves high (76%) completeness and low (20%) contamination, and also performs considerably better than other automated methods recently applied in the literature. We argue that CNNs offer a promising approach to effective automatic identification of low surface brightness tidal debris features in and around galaxies. When applied to forthcoming deep wide-field imaging surveys (e.g. LSST, Euclid), CNNs have the potential to provide a several order-of-magnitude increase in the sample size of morphologically-perturbed galaxies and thereby facilitate a much-anticipated revolution in terms of quantitative low surface brightness science., 16 pages, 18 figures, accepted for publication in MNRAS

We present a spatially resolved study of the relation between dust and metallicity in the nearby spiral galaxies M 101 (NGC 5457) and NGC 628 (M 74). We explore the relation between the chemical abundances of their gas and stars with their dust content and their chemical evolution. The empirical spatially resolved oxygen effective yield and the gas-to-dust mass ratio (GDR) across both disc galaxies are derived, sampling 1 dex in oxygen abundance. We find that the metal budget of the NGC 628 disc and most of the M 101 disc appears consistent with the predictions of the simple model of chemical evolution for an oxygen yield between half and one solar, whereas the outermost region (R ≥ 0.8 R 25 ) of M 101 presents deviations suggesting the presence of gas flows. The GDR-metallicity relation shows a two slopes behaviour, with a break at 12 + log(O/H) ≈ 8.4, a critical metallicity predicted by theoretical dust models when stardust production equals grain growth. A relation between GDR and the fraction of molecular to total gas, ΣH 2 }/Σ gas is also found. We suggest an empirical relationship between GDR and the combination of 12 + log(O/H), for metallicity, and ΣH 2 /Σ gas , a proxy for the molecular clouds fraction. The GDR is closely related with metallicity at low abundance and with ΣH 2 /Σ gas for higher metallicities suggesting interstellar medium dust growth. The ratio Σ dust /Σ star correlates well with 12+ log(O/H) and strongly with log(N/O) in both galaxies. For abundances below the critical one, the 'stardust' production gives us a constant value suggesting a stellar dust yield similar to the oxygen yield. © 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., This work was partially supported by the Spanish Ministerio de Economia y Competitividad under grants AYA2016-79724-C4-4-P and AYA2016-79724-C4-3-P, and excellence project PEX2011-FQM-7058 of Junta de Andalucia (Spain). MR acknowledges support by the research projects AYA2014-53506-P and AYA2017-84897-P from the Spanish Ministerio de Economia y Competitividad and Junta de Andalucia grant FQM108; support from the European Regional Development Funds (FEDER) is acknowledged.

We use the Sloan Digital Sky Survey (SDSS)–Gaia catalogue to search for substructure in the stellar halo. The sample comprises 62 133 halo stars with full phase space coordinates and extends out to heliocentric distances of ∼10 kpc. As actions are conserved under slow changes of the potential, they permit identification of groups of stars with a common accretion history. We devise a method to identify halo substructures based on their clustering in action space, using metallicity as a secondary check. This is validated against smooth models and numerical constructed stellar haloes from the Aquarius simulations. We identify 21 substructures in the SDSS–Gaia catalogue, including seven high-significance, high-energy and retrograde ones. We investigate whether the retrograde substructures may be material stripped off the atypical globular cluster ω Centauri. Using a simple model of the accretion of the progenitor of the ω Centauri, we tentatively argue for the possible association of up to five of our new substructures (labelled Rg1, Rg3, Rg4, Rg6 and Rg7) with this event. This sets a minimum mass of 5 × 108 M⊙ for the progenitor, so as to bring ω Centauri to its current location in action–energy space. Our proposal can be tested by high-resolution spectroscopy of the candidates to look for the unusual abundance patterns possessed by ω Centauri stars.

Investigating the chemical homogeneity of stars born from the same molecular cloud at virtually the same time is very important for our understanding of the chemical enrichment of the interstellar medium and with it the chemical evolution of the Galaxy. One major cause of inhomogeneities in the abundances of open clusters is stellar evolution of the cluster members. In this work, we investigate variations in the surface chemical composition of member stars of the old open clusterM67 as a possible consequence of atomic diffusion effects taking place during the main-sequence phase. The abundances used are obtained from high-resolution UVES/FLAMES spectra within the framework of the Gaia-ESO Survey. We find that the surface abundances of stars on the main sequence decrease with increasing mass reaching a minimum at the turn-off. After deepening of the convective envelope in subgiant branch stars, the initial surface abundances are restored.We found themeasured abundances to be consistent with the predictions of stellar evolutionary models for a cluster with the age and metallicity of M67. Our findings indicate that atomic diffusion poses a non-negligible constraint on the achievable precision of chemical tagging methods.© 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., This work was supported by Sonderforschungsbereich SFB 881 'The Milky Way System' (subproject B5) of the German Research Foundation (DFG). AB acknowledges support from the SFB 881 visitor program. CBM thanks Elisabetta Caffau for fruitful discussions and Don VandenBerg for providing the isochrones shown in Fig 1. The authors would like to thank the anonymous referee for the constructive comments that helped in improving the present work.r MTC acknowledges the financial support from the Spanish Ministerio de Economia y Competitividad, through grant AYA2016-75931. AD and GT acknowledge support by the Research Council of Lithuania (MIP-082/2015). UH acknowledges support from the Swedish National Space Board (SNSB/Rymdstyrelsen). TB was supported by the project grant 'The New Milky' from the Knut and Alice Wallenberg foundation. RS acknowledges support from the Polish Ministry of Science and Higher Education. SGS and VA acknowledge the support by Fundacao para a Ciencia e Tecnologia (FCT) through national funds and a research grant (project ref. UID/FIS/04434/2013, and PTDC/FISAST/7073/2014). SGS also acknowledges the support from FCT through Investigador FCT contract of reference IF/00028/2014 and POPH/FSE (EC) by FEDER funding through the program 'Programa Operacional de Factores de Competitividade - COMPETE'. VA acknowledges the support from FCT through Investigador FCT contract IF/00650/2015/CP1273/CT0001. EJA acknowledges partial support from the Spanish Ministry for Economy and Competitiveness and FEDER funds through grant AYA2013-40611-P. ARC is supported through an Australian Research Council Discovery Project under grant DP160100637. CL thanks the Swiss National Science Foundation for supporting this research through the Ambizione grant number PZ00P2 168065. AK acknowledges support from the Swedish National Space Board (SNSB). XF acknowledges financial support from Premiale 2015 MITiC (PI B/Garilli). SLM acknowledges support from the Australian Research Council through grant DE140100598. AB acknowledges support from the Millennium Science Initiative (Chilean Ministry of Economy).r Based on data products from observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 188.B-3002. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF/Osservatorio Astrofisico di Arcetri. These data have been obtained from the Gaia-ESO Survey Data Archive, prepared and hosted by the Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Science and Technology Facilities Council. This work was partly supported by the European Union FP7 programme through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge the support from INAF and Ministero dell' Istruzione, dell' Universita' e della Ricerca (MIUR) in the form of the grant 'Premiale VLT 2012', and through PRIN-INAF 2014 'The Gaia-ESO Survey'. The results presented here benefit from discussions held during the Gaia-ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme. This research has made use of the SIMBAD data base, operated at CDS, Strasbourg, France and NASA's Astrophysics Data System.

We report the discovery of a thin stellar stream - which we name the Jet stream - crossing the constellations of Hydra and Pyxis. The discovery was made in data from the SLAMS survey, which comprises deep $g$ and $r$ imaging for a $650$ square degree region above the Galactic disc performed by the CTIO Blanco + DECam. SLAMS photometric catalogues will be made publicly available. The stream is approximately 0.18 degrees wide and 10 degrees long, though it is truncated by the survey footprint. Its colour-magnitude diagram is consistent with an old, metal-poor stellar population at a heliocentric distance of approximately 29 kpc. We corroborate this measurement by identifying a spatially coincident overdensity of likely blue horizontal branch stars at the same distance. There is no obvious candidate for a surviving stream progenitor., MNRAS accepted version

We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies., 19 pages, 14 figures, Accepted MNRAS

We present a comparison of nine galaxy formation models, eight semi-analytical, and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of comoving width 125h−1 Mpc, with a dark-matter particle mass of 1.24 × 109h−1M) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some ‘memory’ of any previous calibration that served as the starting point (especially for the manually tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3σ. The second calibration extended the observational data to include the z = 2 SMF alongside the z ∼ 0 star formation rate function, cold gas mass, and the black hole–bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models, we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

We present the results of two-component (disc+bar) and three-component (disc+bar+bulge) multiwavelength 2D photometric decompositions of barred galaxies in five SDSS bands ($ugriz$). This sample of $\sim$3,500 nearby ($z, 26 pages, 21 figures; accepted for publication in MNRAS

We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess the significance of halo substructures based on their clustering in velocity space. The two most substantial structures are multiple wraps of a stream which has undergone considerable phase mixing (S1, with 94 members) and a kinematically cold stream (S2, with 61 members). The member stars of S1 have a median position of ($X,Y,Z$) = ($8.12, -0.22, 2.75$) kpc and a median metallicity of [Fe/H] $= -1.78$. The stars of S2 have median coordinates ($X,Y,Z$) = ($8.66, 0.30, 0.77$) kpc and a median metallicity of [Fe/H] $= -1.91$. They lie in velocity space close to some of the stars in the stream reported by Helmi et al. (1999). By modelling, we estimate that both structures had progenitors with virial masses $\approx 10^{10} M_\odot$ and infall times $\gtrsim 9$ Gyr ago. Using abundance matching, these correspond to stellar masses between $10^6$ and $10^7 M_\odot$. These are somewhat larger than the masses inferred through the mass-metallicity relation by factors of 5 to 15. Additionally, we identify two further substructures (S3 and S4 with 55 and 40 members) and two clusters or moving groups (C1 and C2 with 24 and 12) members. In all 6 cases, clustering in kinematics is found to correspond to clustering in both configuration space and metallicity, adding credence to the reliability of our detections., MNRAS, in press

We report the discovery of a nearby dwarf galaxy in the constellation of Hydrus, between the Large and the Small Magellanic Clouds. Hydrus 1 is a mildy elliptical ultra-faint system with luminosity $M_V\sim$ -4.7 and size $\sim$ 50 pc, located 28 kpc from the Sun and 24 kpc from the LMC. From spectroscopy of $\sim$ 30 member stars, we measure a velocity dispersion of 2.7 km/s and find tentative evidence for a radial velocity gradient consistent with 3 km/s rotation. Hydrus 1's velocity dispersion indicates that the system is dark matter dominated, but its dynamical mass-to-light ratio M/L $\sim$ 66 is significantly smaller than typical for ultra-faint dwarfs at similar luminosity. The kinematics and spatial position of Hydrus~1 make it a very plausible member of the family of satellites brought into the Milky Way by the Magellanic Clouds. While Hydrus 1's proximity and well-measured kinematics make it a promising target for dark matter annihilation searches, we find no evidence for significant gamma-ray emission from Hydrus 1. The new dwarf is a metal-poor galaxy with a mean metallicity [Fe/H]=-2.5 and [Fe/H] spread of 0.4 dex, similar to other systems of similar luminosity. Alpha-abundances of Hyi 1 members indicate that star-formation was extended, lasting between 0.1 and 1 Gyr, with self-enrichment dominated by SN Ia. The dwarf also hosts a highly carbon-enhanced extremely metal-poor star with [Fe/H] $\sim$ -3.2 and [C/Fe] $\sim$ +3.0., submitted to MNRAS; 21 pages

We combine Gaia data release 1 astrometry with Sloan Digital Sky Survey (SDSS) images taken some ~10-15 years earlier, to measure proper motions of stars in the halo of our Galaxy. The SDSS-Gaia proper motions have typical statistical errors of 2 mas/yr down to r ~ 20 mag, and are robust to variations with magnitude and colour. Armed with this exquisite set of halo proper motions, we identify RR Lyrae, blue horizontal branch (BHB), and K giant stars in the halo, and measure their net rotation with respect to the Galactic disc. We find evidence for a gently rotating prograde signal (< $V_��$ > ~ 5-25 km/s) in the halo stars, which shows little variation with Galactocentric radius out to 50 kpc. The average rotation signal for the three populations is < $V_��$ > = 14 +/- 2 +/- 10 (syst.) km/s. There is also tentative evidence for a kinematic correlation with metallicity, whereby the metal richer BHB and K giant stars have slightly stronger prograde rotation than the metal poorer stars. Using the Auriga simulation suite we find that the old (T >10 Gyr) stars in the simulated halos exhibit mild prograde rotation, with little dependence on radius or metallicity, in general agreement with the observations. The weak halo rotation suggests that the Milky Way has a minor in situ halo component, and has undergone a relatively quiet accretion history., 17 pages, 13 figures. MNRAS accepted

We identify a halo substructure in the Tycho Gaia Astrometric Solution (TGAS) data set, cross-matched with the Radial Velocity Experiment (RAVE-on) data release. After quality cuts, the stars with large radial action (JR > 800 km s−1 kpc) are extracted. A subset of these stars is clustered in longitude and velocity and can be selected with further cuts. The 14 stars are centred on (X, Y, Z) ≈ (9.0, −1.0, −0.6) kpc and form a coherently moving structure in the halo with median (vR, vϕ, vz) = (167.33, 0.86, −94.85) km s−1. They are all metal-poor giants with median [Fe/H] = −0.83. To guard against the effects of distance errors, we compute spectrophotometric distances for 8 out of the 14 stars where this is possible. We find that six of the stars are still comoving. These six stars also have a much tighter [Fe/H] distribution ∼−0.7 with one exception ([Fe/H] = −2.12). We conclude that the existence of the comoving cluster is stable against changes in distance estimation and conjecture that this is the dissolving remnant of a globular cluster accreted long ago., GCM thanks the Boustany Foundation, Cambridge Commonwealth, European & International Trust, and Issac Newton Studentship for their support on his work. SEK thanks the United Kingdom Science and Technology Council (STFC) for the award of an Ernest Rutherford fellowship (grant number ST/N004493/1).

We present a clustering comparison of 12 galaxy formation models [including semi-analytic models (SAMs) and halo occupation distribution (HOD) models] all run on halo catalogues and merger trees extracted from a single Λ cold dark matter N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the two-point correlation functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 h(−1) Mpc. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs., Monthly Notices of the Royal Astronomical Society, 469 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Monthly Notices of the Royal Astronomical Society, 464 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged], Comment: 13 pages, 12 figures, MNRAS submitted

Observational searches for faint active nuclei at $z > 6$ have been extremely elusive, with a few candidates whose high-$z$ nature is still to be confirmed. Interpreting this lack of detections is crucial to improve our understanding of high-$z$ supermassive black holes (SMBHs) formation and growth. In this work, we present a model for the emission of accreting BHs in the X-ray band, taking into account super-Eddington accretion, which can be very common in gas-rich systems at high-$z$. We compute the spectral energy distribution for a sample of active galaxies simulated in a cosmological context, which represent the progenitors of a $z \sim 6$ SMBH with $M_{\rm BH} \sim 10^9 \, M_\odot$. We find an average Compton thick fraction of $\sim 45\%$ and large typical column densities ($N_H \gtrsim 10^{23} \rm \, cm^2$). However, faint progenitors are still luminous enough to be detected in the X-ray band of current surveys. Even accounting for a maximum obscuration effect, the number of detectable BHs is reduced at most by a factor 2. In our simulated sample, observations of faint quasars are mainly limited by their very low active fraction ($f_{\rm act} \sim 1 \%$), which is the result of short, super-critical growth episodes. We suggest that to detect high-$z$ SMBHs progenitors, large area surveys with shallower sensitivities, such as Cosmos Legacy and XMM-LSS+XXL, are to be preferred with respect to deep surveys probing smaller fields, such as CDF-S.

We present the discovery of stellar tidal tails around the Large and the Small Magellanic Clouds in the Gaia DR1 data. In between the Clouds, their tidal arms are stretched towards each other to form an almost continuous stellar bridge. Our analysis relies on the exquisite quality of the Gaia's photometric catalogue to build detailed star-count maps of the Clouds. We demonstrate that the Gaia DR1 data can be used to detect variable stars across the whole sky, and in particular, RR Lyrae stars in and around the LMC and the SMC. Additionally, we use a combination of Gaia and Gale to follow the distribution of Young Main Sequence stars in the Magellanic System. Viewed by Gaia, the Clouds show unmistakable signs of interaction. Around the LMC, clumps of RR Lyrae are observable as far as ~20 degrees, in agreement with the most recent map of Mira-like stars reported in Deason et al (2016). The SMC's outer stellar density contours show a characteristic S-shape, symptomatic of the on-set of tidal stripping. Beyond several degrees from the center of the dwarf, the Gaia RR Lyrae stars trace the Cloud's trailing arm, extending towards the LMC. This stellar tidal tail mapped with RR Lyrae is not aligned with the gaseous Magellanic Bridge, and is shifted by some ~5 degrees from the Young Main Sequence bridge. We use the offset between the bridges to put constraints on the density of the hot gaseous corona of the Milky Way., Comment: Submitted to MNRAS

We present a population study of the star formation history of 1244 Type 2 AGN host galaxies, compared to 6107 inactive galaxies. A Bayesian method is used to determine individual galaxy star formation histories, which are then collated to visualise the distribution for quenching and quenched galaxies within each population. We find evidence for some of the Type 2 AGN host galaxies having undergone a rapid drop in their star formation rate within the last 2 Gyr. AGN feedback is therefore important at least for this population of galaxies. This result is not seen for the quenching and quenched inactive galaxies whose star formation histories are dominated by the effects of downsizing at earlier epochs, a secondary effect for the AGN host galaxies. We show that histories of rapid quenching cannot account fully for the quenching of all the star formation in a galaxy's lifetime across the population of quenched AGN host galaxies, and that histories of slower quenching, attributed to secular (non-violent) evolution, are also key in their evolution. This is in agreement with recent results showing both merger-driven and non-merger processes are contributing to the co-evolution of galaxies and supermassive black holes. The availability of gas in the reservoirs of a galaxy, and its ability to be replenished, appear to be the key drivers behind this co-evolution., 12 pages, 11 figures

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed)., 32 pages, 15 figures; a revised Figure 1 and minor changes, to match the published MNRAS version

Approximately 20% of early-type galaxies host small nuclear stellar discs that are tens to a few hundred parsecs in size. Such discs are expected to be easily disrupted during major galactic encounters, hence their age serve to constrain their assembly history. We use VIMOS integral-field spectroscopic observations for the intermediate-mass E0 galaxy NGC4458 and age-date its nuclear disc via high-resolution fitting of various model spectra. We find that the nuclear disc is at least 6 Gyr old. A clue to gain narrow limits to the stellar age is our knowledge of the nuclear disc contribution to the central surface brightness. The presence of an old nuclear disk, or the absence of disruptive encounters since z~0.6, for a small galaxy such as NGC4458 which belongs to the Virgo cluster, may be consistent with a hierarchical picture for galaxy formation where the smallest galaxies assembles earlier and the crowded galactic environments reduce the incidence of galaxy mergers. On the other hand, NGC4458 displays little or no bulk rotation except for a central kpc-scale kinematically-decoupled core. Slow rotation and decoupled core are usually explained in terms of mergers. The presence and age of the nuclear disc constraint these mergers to have happened at high redshift., Comment: 10 pages, 8 figures. Accepted for publication on MNRAS

We explore the possibility of detecting hydrogen radio recombination lines from 0 < z < 10 quasars. We compute the expected Hnalpha flux densities as a function of absolute magnitude and redshift by considering (i) the range of observed AGN spectral indices from UV to X-ray bands, (ii) secondary ionizations from X-ray photons, and (iii) stimulated emission due to nonthermal radiation. All these effects are important to determine the line fluxes. We find that the combination of slopes: alpha_X,hard = -1.11, alpha_X,soft = -0.7, alpha_EUV = -1.3, alpha_UV = -1.7, maximizes the expected flux, f_Hnalpha = 10 microJy for z = 7 quasars with M_AB = -27 in the n = 50 lines; allowed SED variations produce variations by a factor of 3 around this value. Secondaries boost the line intensity by a factor of 2 to 4, while stimulated emission in high-z quasars with M_AB = -26 provides an extra boost to RRL flux observed at nu = 1 GHz if recombinations arise in HII regions with T_e = 10^3-5 K, n_e = 10^3-5 cm^-3. We compute the sensitivity required for a 5sigma detection of Hnalpha lines using the SKA, finding that the SKA-MID could detect sources with M_AB < -27 (M_AB < -26) at z < 8 (z < 3) in less than 100 hrs of observing time. These observations could open new paths to searches for obscured SMBH progenitors, complementing X-ray, optical/IR and sub-mm surveys., 11 pages, 9 figures; to be published in Monthly Notices of the Royal Astronomical Society Main Journal

We present the first results of a survey for high redshift, z $\ge$ 6, quasars using izY multi-colour photometric observations from the Dark Energy Survey (DES). Here we report the discovery and spectroscopic confirmation of the $\rm z_{AB}, Y_{AB}$ = 20.2, 20.2 (M$_{1450}$ = $-$26.5) quasar DES J0454$-$4448 with an emission line redshift of z = 6.10$\pm$0.03 and a HI near zone size of 4.6 $\pm$ 1.7 Mpc.The quasar was selected as an i-band drop out with i$-$z = 2.46 and z$_{AB} < 21.5$ from an area of $\rm \sim$300 deg$^2$. It is the brightest of our 43 candidates and was identified for follow-up spectroscopically solely based on the DES i$-$z and z$-$Y colours. The quasar is detected by WISE and has $W1_{AB} = 19.68$. The discovery of one spectroscopically confirmed quasar with 5.7 $$ 50-100 new quasars with z $>$ 6 including 3-10 with z $>$ 7 dramatically increasing the numbers of quasars currently known that are suitable for detailed studies including determination of the neutral HI fraction of the intergalactic medium (IGM) during the epoch of Hydrogen reionization., Comment: 9 pages, 13 figures, this is a pre-copyedited, author-produced version of an article accepted for publication in MNRAS following peer review

Monthly Notices of the Royal Astronomical Society, 448 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $\mu_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $\alpha_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $\alpha_{\rm CO}(z=2)\sim 1$ versus $\alpha_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$\mu_*$ systems., Comment: 22 pages, MNRAS accepted