We develop a data-driven model to map stellar parameters (effective temperature, surface gravity and metallicity) accurately and precisely to broad-band stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a dataset with spectroscopic labels from LAMOST, APOGEE and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS 1, 2MASS and WISE. Testing the model on these datasets leads to an excellent fit and a precise -- and by construction accurate -- prediction of the color-magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, temperature-dependent reddening vector. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model may be obtained at https://doi.org/10.5281/zenodo.3902382.

Stars born at the same time in the same place should have formed from gas of the same element composition. But most stars subsequently disperse from their birth siblings, in orbit and orbital phase, becoming 'field stars'. Here we explore and provide direct observational evidence for this process in the Milky Way disc, by quantifying the probability that orbit-similarity among stars implies indistinguishable metallicity. We define the orbit similarity among stars through their distance in action-angle space, $\Delta (J,\theta)$, and their abundance similarity simply by $\Delta$[Fe/H]. Analyzing a sample of main sequence stars from Gaia DR2 and LAMOST, we find an excess of pairs with the same metallicity ($\Delta\mathrm{[Fe/H]}<0.1$) that extends to remarkably large separations in $\Delta (J,\theta)$ that correspond to nearly 1 kpc distances. We assess the significance of this effect through a mock sample, drawn from a smooth and phase-mixed orbit distribution. Through grouping such star pairs into associations with a friend-of-friends algorithm linked by $\Delta (J,\theta)$, we find 100s of mono-abundance groups with $\ge 3$ (to $\gtrsim 20$) members; these groups -- some clusters, some spread across the sky -- are over an order-of-magnitude more abundant than expected for a smooth phase-space distribution, suggesting that we are witnessing the 'dissolution' of stellar birth associations into the field.

We present a mock stellar catalog, matching in volume, depth and data model the content of the planned Gaia early data release 3 (Gaia EDR3). We have generated our catalog (GeDR3mock) using galaxia, a tool to sample stars from an underlying Milky Way (MW) model or from N-body data. We used an updated Besan\c{c}on Galactic model together with the latest PARSEC stellar evolutionary tracks, now also including white dwarfs. We added the Magellanic clouds and realistic open clusters with internal rotation. We empirically modelled uncertainties based on Gaia DR2 (GDR2) and scaled them according to the longer baseline in Gaia EDR3. The apparent magnitudes were reddened according to a new selection of 3D extinction maps. To help with the Gaia selection function we provide all-sky magnitude limit maps in G and BP for a few relevant GDR2 subsets together with the routines to produce these maps for user-defined subsets. We supplement the catalog with photometry and extinctions in non-Gaia bands. The catalog is available in the Virtual Observatory and can be queried just like the actual Gaia EDR3 will be. We highlight a few capabilities of the Astronomy Data Query Language (ADQL) with educative catalog queries. We use the data extracted from those queries to compare GeDR3mock to GDR2, which emphasises the importance of adding observational noise to the mock data. Since the underlying truth, e.g. stellar parameters, is know in GeDR3mock, it can be used to construct priors as well as mock data tests for parameter estimation. All code, models and data used to produce GeDR3mock are linked and contained in galaxia_wrap, a python package, representing a fast galactic forward model, able to project MW models and N-body data into realistic Gaia observables.

Constraining parameters such as the initial mass function high-mass slope and the frequency of type Ia supernovae is of critical importance in the ongoing quest to understand galactic physics and create realistic hydrodynamical simulations. In this paper, we demonstrate a method to precisely determine these using individual chemical abundances from a large set of stars, coupled with some estimate of their ages. Inference is performed via the simple chemical evolution model Chempy in a Bayesian framework, marginalizing over each star's specific interstellar medium parameters, including an element-specific `model error' parameter to account for inadequacies in our model. Hamiltonian Monte Carlo (HMC) methods are used to sample the posterior function, made possible by replacing Chempy with a trained neural network at negligible error. The approach is tested using data from both Chempy and the IllustrisTNG simulation, showing sub-percent agreement between inferred and true parameters using data from up to 1600 individual stellar abundances. For IllustrisTNG, strongest constraints are obtained from metal ratios, competitive with those from other methods including star counts. Analysis using a different set of nucleosynthetic yields shows that incorrectly assumed yield models can give non-negligible bias in the derived parameters; this is reduced by our model errors, which further show how well the yield tables match data. We also find a significant bias from analyzing only a small set of stars, as is often done in current analyses. The method can be easily applied to observational data, giving tight bounds on key galactic parameters from chemical abundances alone.

It is well established that the chemical structure of the Milky Way exhibits a bimodality with respect to the $\alpha$-enhancement of stars at a given [Fe/H]. This has been studied largely based on a bulk $\alpha$ abundance, computed as a summary of several individual $\alpha$-elements. Inspired by the expected subtle differences in their nucleosynthetic origins, here we probe the higher level of granularity encoded in the inter-family [Mg/Si] abundance ratio. Using a large sample of stars with APOGEE abundance measurements, we first demonstrate that there is additional information in this ratio beyond what is already apparent in [$\alpha$/Fe] and [Fe/H] alone. We then consider Gaia astrometry and stellar age estimates to empirically characterize the relationships between [Mg/Si] and various stellar properties. We find small but significant trends between this ratio and $\alpha$-enhancement, age, [Fe/H], location in the Galaxy, and orbital actions. To connect these observed [Mg/Si] variations to a physical origin, we attempt to predict the Mg and Si abundances of stars with the galactic chemical evolution model Chempy. We find that we are unable to reproduce abundances for the stars that we fit, which highlights tensions between the yield tables, the chemical evolution model, and the data. We conclude that a more data-driven approach to nucleosynthetic yield tables and chemical evolution modeling is necessary to maximize insights from large spectroscopic surveys.

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.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way.

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.

We will study the formation history of the Milky Way, and the earliest phases of its chemical enrichment, with a sample of more than 1.5 million stars at high galactic latitude. Elemental abundances of up to 20 elements with a precision of better than 0.2 dex will be derived for these stars. The sample will include members of kinematically coherent substructures, which we will associate with their possible birthplaces by means of their abundance signatures and kinematics, allowing us to test models of galaxy formation. Our target catalogue is also expected to contain 30,000 stars at a metallicity of less than one hundredth that of the Sun. This sample will therefore be almost a factor of 100 larger than currently existing samples of metal-poor stars for which precise elemental abundances are available (determined from high-resolution spectroscopy), enabling us to study the early chemical evolution of the Milky Way in unprecedented detail.

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.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

We aim at the discovery of new planetary systems by exploiting the transit light curve results from TESS orbital observatory's Sector 1 and 2 observations and validating them with precise Doppler measurements obtained from archival HARPS data. Taking advantage of the reported TESS transit events around GJ143 (TOI 186) and HD23472 (TOI 174) we model their HARPS precise Doppler measurements and derive orbital parameters for these two systems. For the GJ143 system TESS has reported only a single transit, and thus its period is unconstrained from photometry. Our RV analysis of GJ143 reveal the full Keplerian solution of the system, which is consistent with an eccentric planet with a mass almost twice that of Neptune and a period of $P_{\rm b}$ = $35.59_{-0.01}^{+0.01}$ days. Our estimates of the GJ143 b planet are fully consistent with the transit timing from TESS. We confirm the two-planet system around HD23472, which according to our analysis is composed of two Neptune mass planets in a possible 5:3 MMR.

For the vast majority of stars in the second Gaia data release, reliable distances cannot be obtained by inverting the parallax. A correct inference procedure must instead be used to account for the nonlinearity of the transformation and the asymmetry of the resulting probability distribution. Here we infer distances to essentially all 1.33 billion stars with parallaxes published in the second \gaia\ data release. This is done using a weak distance prior that varies smoothly as a function of Galactic longitude and latitude according to a Galaxy model. The irreducible uncertainty in the distance estimate is characterized by the lower and upper bounds of an asymmetric confidence interval. Although more precise distances can be estimated for a subset of the stars using additional data (such as photometry), our goal is to provide purely geometric distance estimates, independent of assumptions about the physical properties of, or interstellar extinction towards, individual stars. We analyse the characteristics of the catalogue and validate it using clusters. The catalogue can be queried on the Gaia archive using ADQL at http://gea.esac.esa.int/archive/ and downloaded from http://www.mpia.de/~calj/gdr2_distances.html .

To fully harvest the rich library of stellar elemental abundance data available, we require reliable models that facilitate our interpretation of them. Galactic chemical evolution (GCE) models are one such set, and a key part of which are the selection of chemical yields from different nucleosynthetic enrichment channels, predominantly asymptotic giant branch (AGB) stars, Type Ia supernovae (SNe Ia), and core-collapse supernovae (CC-SNe). Here, we present a scoring system for yield tables based on their ability to reproduce proto-solar abundances within a simple parametrisation of the GCE modelling software Chempy, which marginalises over galactic parameters describing simple stellar populations (SSPs) and interstellar medium physics. Two statistical scoring methods are presented, based on Bayesian evidence and leave-one-out cross-validation and are applied to five CC-SN tables, (a) for all mutually available elements and (b) for a subset of the 9 most abundant elements. We find that the yields of Prantzos et al. (P18, including stellar rotation) and Chieffi & Limongi (C04) best reproduce proto-solar abundances for the two cases, respectively. The inferred best-fit SSP parameters for (b) are $\alpha_\mathrm{IMF}=-2.45^{+0.15}_{-0.11}$ for the initial mass function high-mass slope and $\mathrm{N}_\mathrm{Ia}=1.29^{+0.45}_{-0.31}\times10^{-3}$ M$_\odot^{-1}$ for the SN Ia normalisation, which are broadly consistent across tested yield tables. Additionally, we demonstrate how Chempy can be used to dramatically improve elemental abundance predictions of hydrodynamical simulations by plugging tailored best-fit SSP parameters into a Milky Way analogue from Gutcke & Springel. Our code, including a comprehensive tutorial, is freely available and can additionally provide SSP enrichment tables for any combination of parameters and yield tables.

Passing stars may play an important role in the evolution of our solar system. We search for close stellar encounters to the Sun among all 7.2 million stars in Gaia-DR2 that have six-dimensional phase space data. We characterize encounters by integrating their orbits through a Galactic potential and propagating the correlated uncertainties via a Monte Carlo resampling. After filtering to remove spurious data, we find 694 stars that have median (over uncertainties) closest encounter distances within 5 pc, all occurring within 15 Myr from now. 26 of these have at least a 50% chance of coming closer than 1 pc (and 7 within 0.5 pc), all but one of which are newly discovered here. We further confirm some and refute several other previously-identified encounters, confirming suspicions about their data. The closest encounter in the sample is Gl 710, which has a 95% probability of coming closer than 0.08 pc (17 000 AU). Taking mass estimates from Gaia astrometry and multiband photometry for essentially all encounters, we find that Gl 710 also has the largest impulse on the Oort cloud. Using a Galaxy model, we compute the completeness of the Gaia-DR2 encountering sample as a function of perihelion time and distance. Only 15% of encounters within 5 pc occurring within +/- 5 Myr of now have been identified, mostly due to the lack of radial velocities for faint and/or cool stars. Accounting for the incompleteness, we infer the present rate of encounters within 1 pc to be 19.7 +/- 2.2 per Myr, a quantity expected to scale quadratically with the encounter distance out to at least several pc. Spuriously large parallaxes in our sample from imperfect filtering would tend to inflate both the number of encounters found and this inferred rate. The magnitude of this effect is hard to quantify.

We present a mock catalog of Milky Way stars, matching in volume and depth the content of the Gaia data release 2 (GDR2). We generated our catalog using Galaxia, a tool to sample stars from a Besancon Galactic model, together with a realistic 3D dust extinction map. The catalog mimicks the complete GDR2 data model and contains most of the entries in the Gaia source catalog: 5-parameter astrometry, 3-band photometry, radial velocities, stellar parameters, and associated scaled nominal uncertainty estimates. In addition, we supplemented the catalog with extinctions and photometry for non-Gaia bands. This catalog can be used to prepare GDR2 queries in a realistic runtime environment, and it can serve as a Galactic model against which to compare the actual GDR2 data in the space of observables. The catalog is hosted through the virtual observatory GAVO's Heidelberg data center service and thus can be queried using ADQL as for GDR2 data.

In the age of high-resolution spectroscopic stellar surveys of the Milky Way, the number of stars with detailed abundances of multiple elements is rapidly increasing. These elemental abundances are directly influenced by the evolutionary history of the Galaxy, but this can be difficult to interpret without an absolute timeline of the abundance enrichment. We present age-abundance trends for [M/H], [{\alpha}/M], and 17 individual elements using a sample of 721 solar neighbourhood Hipparcos red giant stars observed by APOGEE. These age trends are determined through a Bayesian hierarchical modelling method presented by Feuillet et al. (2016). We confirm that the [{\alpha}/M]- age relation in the solar neighbourhood is steep and relatively narrow (0.20 dex age dispersion), as are the [O/M]- and [Mg/M]-age relations. The age trend of [C/N] is steep and smooth, consistent with stellar evolution. The [M/H]-age relation has a mean age dispersion of 0.28 dex and a complex overall structure. The oldest stars in our sample are those with the lowest and highest metallicities, while the youngest stars are those with solar metallicity. These results provide strong constraints on theoretical models of Galactic chemical evolution (GCE). We compare them to the predictions of one-zone GCE mod- els and multi-zone mixtures, both analytic and numerical. These comparisons support the hypothesis that the solar neighbourhood is composed of stars born at a range of Galactocentric radii, and that the most metal-rich stars likely migrated from a region with earlier and more rapid star formation such as the inner Galaxy.

In this thesis, the field star Initial Mass Function (IMF) and chemical evolution parameters for the Milky Way (MW) are derived using a forward modelling technique in combination with Bayesian statistics. Starting from a local MM disc model, observations of stellar samples in the Solar Neighbourhood are synthesised and compared to the corresponding volume-complete observational samples of Hipparcos stars. The resulting IMF, derived from observations in the range from 0.5 to 8Msun, is a two-slope broken power law with powers of -1.49 +- 0.08 and -3.02 +- 0.06 for the low-mass slope and the high-mass slope, respectively, with a break at 1.39 +- 0.05Msun. In order to constrain the IMF for stars more massive than 8Msun, a fast and flexible chemical enrichment code, Chempy, was developed, which is also able to reproduce spatial and stellar population selections of observational samples. The inferred high-mass slope for stellar masses above 6Msun is -2.28 +- 0.09, accounting for the systematic effects of different yield sets from the literature. This shows that constraints from chemical modelling, similarly to hydrodynamical simulations of the Galaxy, demand a Salpeter high-mass index. This is hard to recover from star count analysis given the rareness of high-mass stars.

Elemental abundances of stars are the result of the complex enrichment history of their galaxy. Interpretation of observed abundances requires flexible modeling tools to explore and quantify the information about Galactic chemical evolution (GCE) stored in such data. Here we present Chempy, a newly developed code for GCE modeling, representing a parametrized open one-zone model within a Bayesian framework. A Chempy model is specified by a set of 5-10 parameters that describe the effective galaxy evolution along with the stellar and star-formation physics: e.g. the star-formation history, the feedback efficiency, the stellar initial mass function (IMF) and the incidence of supernova type Ia (SN Ia). Unlike established approaches, Chempy can sample the posterior probability distribution in the full model parameter space and test data-model matches for different nucleosynthetic yield sets. We extend Chempy to a multi-zone scheme. As an illustrative application, we show that interesting parameter constraints result from only the ages and elemental abundances of Sun, Arcturus and the present-day interstellar medium (ISM). For the first time, we use such information to infer IMF parameter via GCE modeling, where we properly marginalize over nuisance parameters and account for different yield sets. We find that of the IMF $11.6_{-1.6}^{+2.1}$ % explodes as core-collapse SN, compatible with Salpeter 1955. We also constrain the incidence of SN Ia per 10^3 Msun to 0.5-1.4. At the same time, this Chempy application shows persistent discrepancies between predicted and observed abundances for some elements, irrespective of the chosen yield set. These cannot be remedied by any variations of Chempy's parameters and could be an indication for missing nucleosynthetic channels. Chempy should be a powerful tool to confront predictions from stellar nucleosynthesis with far more complex abundance data sets.

Our detailed analytic local disc model (JJ-model) quantifies the interrelation between kinematic properties (e.g. velocity dispersions and asymmetric drift), spatial parameters (scale-lengths and vertical density profiles), and properties of stellar sub-populations (age and abundance distributions). Any consistent radial extension of the disc evolution model should predict specific features in the different distribution functions and in their correlations. Large spectroscopic surveys (SEGUE, RAVE, APOGEE, Gaia-ESO) allow significant constraints on the long-term evolution of the thin disc. We discuss the qualitative difference of correlations (like the alpha-enhancement as function of metallicity) and distribution functions (e.g. in [Mg/H] or [Fe/H]) for the construction of a disc model. In the framework of the JJ-model we build a local chemical enrichment model and show that significant vertical gradients for main sequence and red clump stars are expected in the thin disc. A Jeans analysis of the asymmetric drift provides a link to the radial structure of the disc. The derived metallicity-dependent radial scale-lengths can be combined in the future with the abundance distributions at different Galactocentric distances to construct full disc models. We expect to be able to constrain possible scenarios of inside-out growth of the thin disc and to characterise those populations, which require significant radial migration.

A new sample of stars, representative of the solar neighbourhood luminosity function, is constructed from the Hipparcos catalogue and the Fifth Catalogue of Nearby Stars. We have cross-matched to sources in the 2MASS catalogue so that for all stars individually determined Near Infrared photometry (NIR) is available on a homogeneous system (typically K_s). The spatial completeness of the sample has been carefully determined by statistical methods, and the NIR luminosity function of the stars has been derived by direct star counts. We find a local volume luminosity of 0.121 +/- 0.004 L_K_sun/(pc**3), corresponding to a volumetric mass-to-light ratio of M/L_K = 0.31 +/- 0.02 M_sun/L_K_sun, where giants contribute 80 per cent to the light but less than 2 per cent to the stellar mass. We derive the surface brightness of the solar cylinder with the help of a vertical disc model. We find a surface brightness of 99 L_K_sun/(pc**2) with an uncertainty of approximately 10 %. This corresponds to a mass-to-light ratio for the solar cylinder of M/L_K = 0.34 M_sun/L_K_sun. The mass-to-light ratio for the solar cylinder is only 10% larger than the local value despite the fact that the local population has a much larger contribution of young stars. It turns out that the effective scale heights of the lower main sequence carrying most of the mass is similar to that of the giants, which are dominating the NIR light. The corresponding colour for the solar cylinder is V-K=2.89 mag compared to the local value of V-K = 2.46 mag. An extrapolation of the local surface brightness to the whole Milky Way yields a total luminosity of M_K = -24.2 mag. The Milky Way falls in the range of K band Tully-Fisher (TF) relations from the literature.

We use our vertical Milky Way disc model together with Galaxia to create mock observations of stellar samples in the solar neighbourhood. We compare these to the corresponding volume complete observational samples of dereddened and binary accounted data from Hipparcos and the Catalogue of Nearby Stars. Sampling the likelihood in the parameter space we determine a new fiducial IMF considering constraints from dwarf and giant stars. The resulting IMF observationally backed in the range from 0.5 to 10 Msun is a two slope broken power law with -1.49 +- 0.08 for the low mass slope, a break at 1.39 +- 0.05 Msun and a high mass slope of -3.02 +- 0.06. The Besancon group also converging to a similar IMF even though their observational sample being quite different to ours shows that the forward modelling technique is a powerful diagnostic to test theoretical concepts like the local field star IMF.