Data-driven methods play an increasingly important role in the field of astrophysics. In the context of large spectroscopic surveys of stars, data-driven methods are key in deducing physical parameters for millions of spectra in a short time. Convolutional Neural Networks (CNNs) enable us to connect observables (e.g., spectra, stellar magnitudes) to physical properties (atmospheric parameters, chemical abundances, or labels in general). We trained a CNN, adopting stellar atmospheric parameters and chemical abundances from APOGEE DR16 (resolution R~22500) data as training set labels. As input we used parts of the intermediate-resolution RAVE DR6 spectra (R~7500) overlapping with the APOGEE DR16 data as well as broad-band ALL_WISE and 2MASS photometry together with Gaia DR2 photometry and parallaxes. We derived precise atmospheric parameters Teff, log(g), [M/H] and chemical abundances of [Fe/H], [alpha/M], [Mg/Fe], [Si/Fe], [Al/Fe], and [Ni/Fe] for 420165 RAVE spectra. The precision typically amounts to 60K in Teff, 0.06 in log(g), and 0.02-0.04 dex for individual chemical abundances. Incorporating photometry and astrometry as additional constraints substantially improves the results in terms of accuracy and precision of the derived labels, as long as we operate in parts of the parameter space that is well covered by the training sample. Scientific validation confirms the robustness of the CNN results. We provide a catalog of CNN-trained atmospheric parameters and abundances along with their uncertainties for 420165 stars in the RAVE survey. CNN-based methods provide a powerful way to combine spectroscopic, photometric, and astrometric data, without applying any priors in the form of stellar evolutionary models. The developed procedure can extend the scientific output of RAVE spectra beyond DR6 to on-going and planned surveys such as Gaia RVS and 4MOST.

We combine high-resolution spectroscopic data from APOGEE-2 Survey Data Release 16 (DR16) with broad-band photometric data from several sources, as well as parallaxes from {\it Gaia} Data Release 2 (DR2). Using the Bayesian isochrone-fitting code {\tt StarHorse}, we derive distances, extinctions and astrophysical parameters for around 388,815 APOGEE stars, achieving typical distance uncertainties of $\sim 6\%$ for APOGEE giants, $\sim 2\%$ for APOGEE dwarfs, as well as extinction uncertainties of $\sim 0.07$ mag when all photometric information is available, and $\sim 0.17$ mag if optical photometry is missing. {\tt StarHorse} uncertainties vary with the input spectroscopic catalogue, with the available photometry, and with the parallax uncertainties. To illustrate the impact of our results, we show that, thanks to {\it Gaia} DR2 and the now larger sky coverage of APOGEE-2 (including APOGEE-South), we obtain an extended map of the Galactic plane, providing unprecedented coverage of the disk close to the Galactic mid-plane ($|Z_{Gal}|<1$ kpc) from the Galactic Centre out to $R_{\rm Gal}\sim 20$ kpc. The improvements in statistics as well as distance and extinction uncertainties unveil the presence of the bar in stellar density, as well as the striking chemical duality in the innermost regions of the disk, now clearly extending to the inner bulge. We complement this paper with distances and extinctions for stars in other public released spectroscopic surveys: 324,999 in GALAH DR2, 4,928,715 in LAMOST DR5, 408,894 in RAVE DR6, and 6,095 in GES DR3

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr.

Combining the precise parallaxes and optical photometry delivered by Gaia's second data release (Gaia DR2) with the photometric catalogues of PanSTARRS-1, 2MASS, and AllWISE, we derive Bayesian stellar parameters, distances, and extinctions for 265 million stars brighter than G=18. Because of the wide wavelength range used, our results substantially improve the accuracy and precision of previous extinction and effective temperature estimates. After cleaning our results for both unreliable input and output data, we retain 137 million stars, for which we achieve a median precision of 5% in distance, 0.20 mag in V-band extinction, and 245 K in effective temperature for G<14, degrading towards fainter magnitudes (12%, 0.20 mag, and 245 K at G=16; 16%, 0.23 mag, and 260 K at G=17, respectively). We find a very good agreement with the asteroseismic surface gravities and distances of 7000 stars in the Kepler, the K2-C3, and the K2-C6 fields, with stellar parameters from the APOGEE survey, as well as with distances to star clusters. Our results are available through the ADQL query interface of the Gaia mirror at the Leibniz-Institut f\"{u}r Astrophysik Potsdam (gaia.aip.de) and as binary tables at data.aip.de. As a first application, in this paper we provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps as a function of distance, and extensive density maps, demonstrating the potential of our value-added dataset for mapping the three-dimensional structure of our Galaxy. In particular, we see a clear manifestation of the Galactic bar in the stellar density distributions, an observation that can almost be considered a direct imaging of the Galactic bar.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The solar neighborhood contains disc stars that have recently crossed spiral arms in the Galaxy. We propose that boundaries in local velocity distributions separate stars that have recently crossed or been perturbed by a particular arm from those that haven't. Ridges in the stellar velocity distributions constructed from the second Gaia data release trace orbits that could have touched nearby spiral arms at apocentre or pericentre. The multiple ridges and arcs seen in local velocity distributions are consistent with the presence of multiple spiral features and different pattern speeds and imply that the outer Galaxy is flocculent rather than grand design.

We present a cosmography analysis of the Local Universe based on the recently released Two-Micron All-Sky Redshift Survey (2MRS). Our method is based on a Bayesian Networks Machine Learning algorithm (the Kigen-code) which self-consistently samples the initial density fluctuations compatible with the observed galaxy distribution and a structure formation model given by second order Lagrangian perturbation theory (2LPT). From the initial conditions we obtain an ensemble of reconstructed density and peculiar velocity fields which characterize the local cosmic structure with high accuracy unveiling nonlinear structures like filaments and voids in detail. Coherent redshift space distortions are consistently corrected within 2LPT. From the ensemble of cross-correlations between the reconstructions and the galaxy field and the variance of the recovered density fields we find that our method is extremely accurate up to k ~ 1 h Mpc^-1 and still yields reliable results down to scales of about 3-4 h^-1 Mpc. The motion of the local group we obtain within ~ 80 h^-1 Mpc (v_LG=522+-86 km s^-1, l_LG=291^o +- 16^o, b_LG=34^o+-8^o) is in good agreement with measurements derived from the CMB and from direct observations of peculiar motions and is consistent with the predictions of LambdaCDM.

We discuss the rest-frame UV emission between 5< z < 7 from the MareNostrum High-z Universe, a SPH simulation done with more than 2 billion particles. Cosmological simulations of galaxy formation generally overpredict the UV restframe luminosity function at high redshift, both at the bright and faint ends. In this Letter we explore a dust attenuation model where a larger extinction is applied to star populations younger than a given age, mimicking the effect of a clumpy interstellar medium. We show that this scenario fits reasonably well both the UV luminosity functions and the UV-continuum slopes derived from observations. The model assumes a large obscuration for stars younger than 25 Myr from the gas clouds where they should be embedded at their formation time. We find that the optical depth in these clouds should be between 30 and 100 times larger than the mean optical depth for the homogeneous part of the interstellar medium. These values are one order of magnitude larger than those estimated in local galaxies. Therefore, we conclude that LambdaCDM predictions for the high-z UV emission can accommodate the current observations if we consider a dust extinction model based on the assumption of a clumpy environment at high redshift.

Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.

We have used GADGET2 to simulate the formation of an elliptical galaxy in a cosmological dark matter halo with mass 3x10^12M_Sun/h. Using a stellar population synthesis model has allowed us to compute magnitudes, colours and surface brightness profiles. We have included a model to follow the growth of a central black hole and we have compared the results of simulations with and without feedback from AGNs. We have studied the interplay between cold gas accretion and merging in the development of galactic morphologies, the link between colour and morphology evolution, the effect of AGN feedback on the photometry of early type galaxies, the redshift evolution in the properties of quasar hosts, and the impact of AGN winds on the chemical enrichment of the intergalactic medium (IGM). We have found that the early phases of galaxy formation are driven by the accretion of cold filamentary flows, which form a disc at the centre of the dark matter halo. When the dark matter halo is sufficiently massive to support the propagation of a stable shock, cold accretion is shut down, and the star formation rate begins to decline. Mergers transform the disc into an elliptical galaxy, but also bring gas into the galaxy. Without a mechanism that removes gas from the merger remnants, the galaxy ends up with blue colours, atypical for its elliptical morphology. AGN feedback can solve this problem even with a fairly low heating efficiency. We have also demonstrated that AGN winds are potentially important for the metal enrichment of the IGM a high redshift.(abridged)

We investigate the origin and evolution of fossil groups in a concordance LCDM cosmological simulation. We consider haloes with masses between $(1-5)\times10^{13} \hMsun$ and study the physical mechanisms that lead to the formation of the large gap in magnitude between the brightest and the second most bright group member, which is typical for these fossil systems. Fossil groups are found to have high dark matter concentrations, which we can relate to their early formation time. The large magnitude-gaps arise after the groups have build up half of their final mass, due to merging of massive group members. We show that the existence of fossil systems is primarily driven by the relatively early infall of massive satellites, and that we do not find a strong environmental dependence for these systems. In addition, we find tentative evidence for fossil group satellites falling in on orbits with typically lower angular momentum, which might lead to a more efficient merger onto the host. We find a population of groups at higher redshifts that go through a ``fossil phase'': a stage where they show a large magnitude-gap, which is terminated by renewed infall from their environment.

We report some results from one of the largest hydrodynamical cosmological simulations of large scale structures that has been done up to date. The MareNostrum Universe SPH simulation consists of 2 billion particles (2 times 1024^3) in a cubic box of 500 h^-1 Mpc on a side. This simulation has been done in the MareNostrum parallel supercomputer at the Barcelona SuperComputer Center. Due to the large simulated volume and good mass resolution, our simulated catalog of dark matter halos comprises more than half a million objects with masses larger than a typical Milky Way galaxy halo. From this dataset we have studied several statistical properties such as the evolution of the halo mass function, the void distribution, the shapes of dark and gas halos and the large scale distribution of baryons.

We show by means of a high-resolution N-body simulation how the mass assembly histories of galaxy-size cold dark matter (CDM) halos depend on environment. Halos in high density environments form earlier and a higher fraction of their mass is assembled in major mergers, compared to low density environments. The distribution of the present-day specific mass aggregation rate is strongly dependent on environment. While in low density environments only ~20% of the halos are not accreting mass at the present epoch, this fraction rises to ~80% at high densities. At z=1 the median of the specific aggregation rate is ~4 times larger than at z=0 and almost independent on environment. All the dependences on environment found here are critically enhanced by local processes associated to subhalos because the fraction of subhalos increases as the environment gets denser. The distribution of the halo specific mass aggregation rate as well as its dependence on environment resemble the relations for the specific star formation rate distribution of galaxies. An analogue of the morphology-density relation is also present at the level of CDM halos, being driven by the halo major merging history. Nevertheless, baryonic processes are necessary in order to explain further details and the evolution of the star formation rate-, color- and morphology-environment relations.

We present the thermodynamics of a nonlocal chiral quark model with separable 4-fermion interaction for the case of U(3) flavor symmetry within a functional integral approach. The four free parameters of the model are fixed by the chiral condensate, and by the pseudoscalar meson properties (pion mass, kaon mass, pion decay constant). We discuss the T=0 equation of state (EoS) which describes quark confinement (zero quark matter pressure) below the critical chemical potential mu_c=333 MeV. The new result of the present approach is that the strange quark deconfinement is separated from the light quark one and occurs only at a higher chemical potential of mu_{c,s}=492 MeV. We compare the resulting EoS to bag model ones for two and three quark flavors, which have the phase transition to the vacuum with zero pressure also at mu_c. We study quark matter stars in general relativity theory assuming beta-equilibrium with electrons and show that for configurations with masses close to the maximum of stability at M=1.62 ... 1.64 M_{solar} strange quark matter can occur.