A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

We use Gaia DR2 magnitudes, colours and parallaxes for stars with G<12 to explore a 15-dimensional space that includes simultaneously the initial mass function (IMF) and a non-parametric star formation history (SFH) for the Galactic disc. This inference is performed by combining the Besancon Galaxy Model fast approximate simulations (BGM FASt) and an approximate Bayesian computation algorithm. We find in Gaia DR2 data an imprint of a star formation burst 2-3 Gyr ago, in the Galactic thin disc domain, and a present star formation rate (SFR) of about 1 Msun. Our results show a decreasing trend of the SFR from 9-10 Gyr to 6-7 Gyr ago. This is consistent with the cosmological star formation quenching observed at redshifts z<1.8. This decreasing trend is followed by a SFR enhancement starting at about 5 Gyr ago and continuing until about 1 Gyr ago which is detected with high statistical significance by discarding the null hypothesis of an exponential SFH with a p-value=0.002. We estimate, from our best fit model, that about the 50% of the mass used to generate stars, along the thin disc life, was expended in the period from 5 to 1 Gyr ago. The timescale and the amount of stellar mass generated during the SFR enhancement event lead us to hypothesise that its origin, currently under investigation, is not intrinsic to the disc. Thus, an external perturbation is needed for its explanation. Additionally we find, for the thin disc, a slope of the IMF of $\alpha_3 \approx 2$ for masses M>1.53 Msun and $\alpha_2 \approx 1.3$ for the mass range between 0,5 and 1,53 Msun. This is the first time that we consider a non-parametric SFH for the thin disc in the Besancon Galaxy Model. This new step, together with the capabilities of the Gaia DR2 parallaxes to break degeneracies between different stellar populations, allow us to better constrain the SFH and the IMF.

SN 1572 (Tycho Brahe's supernova) clearly belongs to the Ia (thermonuclear) type. It was produced by the explosion of a white dwarf in a binary system. Its remnant has been the first of this type to be explored in search of a possible surviving companion, the mass donor that brought the white dwarf to the point of explosion. A high peculiar motion with respect to the stars at the same location in the Galaxy, mainly due to the orbital velocity at the time of the explosion, is a basic criterion for the detection of such companions. Radial velocities from the spectra of the stars close to the geometrical center of Tycho's supernova remnant, plus proper motions of the same stars, obtained by astrometry with the {\it Hubble Space Telescope}, have been used so far. In addition, a detailed chemical analysis of the atmospheres of a sample of candidate stars had been made. However, the distances to the stars, remained uncertain. Now, the Second {\it Gaia} Data Release (DR2) provides unprecedent accurate distances and new proper motions for the stars can be compared with those made from the {\it HST}. We consider the Galactic orbits that the candidate stars to SN companion would have in the future. We do this to explore kinematic peculiarity. We also locate a representative sample of candidate stars in the Toomre diagram. Using the new data, we reevaluate here the status of the candidates suggested thus far, as well as the larger sample of the stars seen in the central region of the remnant.

Context: The Initial Mass Function (IMF) plays a crucial role on galaxy evolution and its implications on star formation theory make it a milestone for the next decade. It is in the intermediate and high mass ranges where the uncertainties of the IMF are larger. This is a major subject of debate and analysis both for Galactic and extragalactic science. Aims: Our goal is to constrain the IMF of the Galactic thin disc population using both Galactic Classical Cepheids and Tycho-2 data. Methods: For the first time the Besan\c{c}on Galaxy Model (BGM) has been used to characterise the Galactic population of the Classical Cepheids. We have modified the age configuration in the youngest populations of the BGM thin disc model to avoid artificial discontinuities in the age distribution of the simulated Cepheids. Three statistical methods, optimized for different mass ranges, have been developed and applied to search for the best IMF that fits the observations. This strategy allows us to quantify variations in the Star Formation History (SFH), the stellar density at Sun position and the thin disc radial scale length. A rigorous treatment of unresolved multiple stellar systems has been undertaken adopting a spatial resolution according to the catalogues used. Results: For intermediate masses, our study favours a composite field-star IMF slope of $\alpha=3.2$ for the local thin disc, excluding flatter values such as the Salpeter IMF ($\alpha=2.35$). Moreover, a constant Star Formation History is definitively excluded, the three statistical methods considered here show that it is inconsistent with the observational data. Conclusions: Using field stars and Galactic Classical Cepheids, we have found, above $1M_\odot$, an IMF steeper than the canonical stellar IMF of associations and young clusters. This result is consistent with the predictions of the Integrated Galactic IMF.

We develop a new theoretical framework to generate Besan\c{c}on Galaxy Model fast approximate simulations (BGM FASt) to address fundamental questions of the Galactic structure and evolution performing multi-parameter inference. As a first application of our strategy we simultaneously infer the IMF, the star formation history and the stellar mass density in the Solar Neighbourhood. The BGM FASt strategy is based on a reweighing scheme, that uses a specific pre-sampled simulation, and on the assumption that the distribution function of the generated stars in the Galaxy can be described by an analytical expression. To validate BGM FASt we execute a set of tests. Finally, we use BGM FASt with an approximate Bayesian computation algorithm to obtain the posterior PDF of the inferred parameters, by comparing synthetic versus Tycho-2 colour-magnitude diagrams. Results: The validation shows a very good agreement between BGM FASt and the standard BGM, with BGM FASt being $\approx 10^4$ times faster. By analysing Tycho-2 data we obtain a thin disc star formation history decreasing in time and a present rate of $1.2 \pm 0.2 M_\odot/yr$. The resulting total stellar mass density in the Solar Neighbourhood is $0.051_{-0.005}^{+0.002} M_\odot/pc^3$ and the local dark matter density is $0.012 \pm 0.001 M_\odot/pc^3$. For the composite IMF we obtain a slope of $\alpha_2={2.1}_{-0.3}^{+0.1}$ in the mass range between $0.5 M_\odot$ and $1.53M_\odot$. The results of the slope at the high mass range are trustable up to $4M_\odot$ and highly depend on the choice of the extinction map (obtaining $\alpha_3={2.9}_{-0.2}^{+0.2}$ and $\alpha_3={3.7}_{-0.2}^{+0.2}$ respectively, for two different extinction maps). Systematic uncertainties are not included. Conclusions: The good performance of BGM FASt demonstrates that it is a very valuable tool to perform multi-parameter inference using Gaia data releases.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of 150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions.

We use the Gaia-TGAS data to compare the transverse velocities in Galactic longitude (coming from proper motions and parallaxes) in the Milky Way disk for negative and positive longitudes as a function of distance. The transverse velocities are strongly asymmetric and deviate significantly from the expectations for an axisymmetric Galaxy. The value and sign of the asymmetry changes at spatial scales of several tens of degrees in Galactic longitude and about 0.5 kpc in distance. The asymmetry is statistically significant at 95% confidence level for 57% of the region probed, which extends up to ~1.2 kpc. A percentage of 24% of the region studied shows absolute differences at this confidence level larger than 5 km/s and 7% larger than 10 km/s. The asymmetry pattern shows mild variations in the vertical direction and with stellar type. A first qualitative comparison with spiral arm models indicates that the arms are unlikely to be the main source of the asymmetry. We briefly discuss alternative origins. This is the first time that global all-sky asymmetries are detected in the Milky Way kinematics, beyond the local neighbourhood, and with a purely astrometric sample.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

We report Herschel Spitzer and Wise observations of 44 z=4.8 optically selected active galactic nuclei (AGNs). This sample contains the highest mass black holes (BHs) at this redshift. Ten of the objects were detected by Herschel and five show emission that is not associated with the AGNs. The star formation (SF) luminosity (LSF) obtained by fitting the spectral energy distribution (SED) with standard SF templates is in the range 10^(46.62-47.21) erg/s corresponding to SF rates of 1090-4240 Msun/yr. Fitting with luminous submillimeter galaxy SEDs gives SF rates that are smaller by 0.05 dex when using all bands and 0.1 dex when ignoring the 250-mic band. A 40K gray-body fits to only the 500-mic fluxes reduce LSF by about a factor two. Stacking analysis of 29 undetected sources give significant signals in all three bands. A SF template fit indicates LSF=10^(46.19-46.23) erg/s depending on the assume AGN contribution. A 40K fit to the stacked 500-mic flux gives LSF=10^(45.95) erg/s. The mean BH mass (mbh) and AGN luminosity (LAGN) of the detected sources are significantly higher than those of the undetected ones. The spectral differences are seen all the way from UV to far infrared wavelengths. The mean optical-UV spectra are similar to the ones predicted for thin accretion disks around BHs with the measured masses and accretion rates. We suggest two alternative explanations to the correlation of LSF, LAGN and mbh, one involving no AGN feedback and the second a moderate feedback that affects, but not totally quench SF in 3/4 of the sources. We compare our LSF and LAGN to lower redshift samples and show a new correlation between LSF and mbh. We also examine several rather speculative ideas about the host galaxy properties including the possibility that the detected sources are above the SF mass sequence (MS) at z=4.8, perhaps in mergers, and most of the undetected sources are on the MS.

We report new Herschel observations of 25 z=4.8 extremely luminous optically selected active galactic nuclei (AGNs). Five of the sources have extremely large star forming (SF) luminosities, L_SF, corresponding to SF rates (SFRs) of 2800-5600 M_sol/yr assuming a Salpeter IMF. The remaining sources have only upper limits on their SFRs but stacking their Herschel images results in a mean SFR of 700 +/- 150 M_sol/yr. The higher SFRs in our sample are comparable to the highest observed values so far, at any redshift. Our sample does not contain obscured AGNs, which enables us to investigate several evolutionary scenarios connecting super-massive black holes and SF activity in the early universe. The most probable scenario is that we are witnessing the peak of SF activity in some sources and the beginning of the post-starburst decline in others. We suggest that all 25 sources, which are at their peak AGN activity, are in large mergers. AGN feedback may be responsible for diminishing the SF activity in 20 of them but is not operating efficiently in 5 others.

We present a detailed investigation of the near-to-far infrared (IR) spectral energy distribution (SED) of a large sample of Spitzer-observed active galactic nuclei (AGN). We fitted the spectra of 51 narrow line Seyfert 1 galaxies (NLS1s) and 64 broad line Seyfert 1 galaxies (BLS1s) using a three component model: a clumpy torus, a dusty narrow line region (NLR) and hot pure-graphite dust located in the outer part of the broad line region (BLR). The fitting is performed on star formation (SF) subtracted SEDs using SF templates that take into account the entire range of possible host galaxy properties. We find that the mid-IR intrinsic emission of NLS1s and BLS1s are very similar, regardless of the AGN luminosity, with long wavelength downturn at around 20-25 um. We present a detailed model of the hot dust component that takes into account the distribution of dust temperature within the clouds and their emission line spectrum. The hot dust continuum provides a very good fit to the observed near-IR continuum spectrum. Most line emission in this component is dramatically suppressed, except MgII(2798A) and HeI lines that are still contributing significantly to the total BLR spectrum. We calculate the covering factors of all the AGN components and show that the covering factor of the hot-dust clouds is about 0.15-0.35, similar to the covering factor of the torus, and is anti-correlated with the source luminosity and the normalized accretion rate.

We fitted the optical to mid-infrared (MIR) spectral energy distributions (SEDs) of ~15000 type-I, 0.75<z<2, active galactic nuclei (AGNs) in an attempt to constrain the properties of the physical component responsible for the rest-frame near-infrared (NIR) emission. We combine optical spectra from the Sloan Digital Sky Survey (SDSS) and MIR photometry from the preliminary data release of the Wide Infrared Survey Explorer (WISE). The sample spans a large range of AGN properties: luminosity, black hole mass, and accretion rate. Our model has two components: a UV-optical continuum source and very hot, pure-graphite dust clouds. We present the luminosity of the hot-dust component and its covering factor, for all sources, and compare it with the intrinsic AGN properties. We find that the hot-dust component is essential to explain the (rest) NIR emission in almost all AGNs in our sample, and that it is consistent with clouds containing pure-graphite grains and located between the dust-free broad line region (BLR) and the "standard" torus. The covering factor of this component has a relatively narrow distribution around a peak value of ~0.13, and it correlates with the AGN bolometric luminosity. We suggest that there is no significant correlation with either black hole mass or normalized accretion rate. The fraction of hot-dust-poor AGNs in our sample is ~15-20%, consistent with previous studies. We do not find a dependence of this fraction on redshift or source luminosity.

We fitted Spitzer/IRS ~2-35 um spectra of 26 luminous QSOs in attempt to define the main emission components. Our model has three major components: a clumpy torus, dusty narrow line region (NLR) clouds and a blackbody-like dust. The models utilize the clumpy torus of Nenkova et al. (2008) and are the first to allow its consistent check in type-I AGNs. Single torus models and combined torus-NLR models fail to fit the spectra of most sources but three component models adequately fit the spectra of all sources. We present torus inclination, cloud distribution, covering factor and torus mass for all sources and compare them with bolometric luminosity, black hole mass and accretion rate. The torus covering factor and mass are found to be correlated with the bolometric luminosity of the sources. We find that a substantial amount of the ~2-7 um radiation originates from a hot dust component, which likely situated in the innermost part of the torus. The luminosity radiated by this component and its covering factor are comparable to those of the torus. We quantify the emission by the NLR clouds and estimate their distance from the center. The distances are ~700 times larger than the dust sublimation radius and the NLR covering factor is about 0.07. The total covering factor by all components is in good agreement with the known AGN type-I:type-II ratio.