We study the shape of the thick disc using photometric data at high and intermediate latitudes from SDSS and 2MASS surveys. We use the population synthesis approach using an ABC-MCMC method to characterize the thick disc shape, scale height, scale length, local density and flare, and we investigate the extend of the thick disc formation period by simulating several formation episodes. We find that the vertical variation in density is not exponential, but much closer to an hyperbolic secant squared. Assuming a single formation epoch, the thick disc is better fitted with a sech2 scale height of 470 pc and a scale length of 2.3 kpc. However if one simulates two successive formation episodes, mimicking an extended formation period, the older episode has a higher scale height and a larger scale length than the younger episode, indicating a contraction during the collapse phase. The scale height decreases from 800 pc to 340 pc, and the scale length from 3.2 kpc to 2 kpc. The star formation increases from the old episode to the young one. During the fitting process, the halo parameters are also determined. The constraint on the halo shows that a transition between a inner and outer halo, if exists, cannot be at a distance of less than about 30 kpc, which is the limit of our investigation using turnoff halo stars. Finally, we show that extrapolating the thick disc towards the bulge region explains well the stellar populations observed there, that there is no longer need to invoke a classical bulge. To explain these results, the most probable scenario for the thick disc is that it formed while the Galaxy was gravitationally collapsing from well mixed gas-rich giant clumps sustained by high turbulence which awhile prevent a thin disc to form, as proposed by Bournaud et al. (2009). This scenario explains the observations in the thick disc region as well as in the bulge region. (abridged)

The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allows us to test hypotheses on the star formation history, star evolution, and chemical and dynamical evolution of the Galaxy. The new version of the Besancon Galaxy model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the initial mass function (IMF) and star formation rate (SFR) of the Galactic disc. We present a new strategy for the generation of thin disc stars, which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. The local dynamical self-consistency is maintained in this new scheme. We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. For the first time, a full sky comparison is performed between BGM and data. This strategy allows us to constrain the IMF slope at high masses, which is found to be close to 3.0 and excludes a shallower slope such as Salpeter's one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 Mo/pc^3 implying that there is no compelling evidence for the significant amount of dark matter in the disc. While the model is fitted to Tycho-2 data, which is a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

A numerical model of the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect for objects defined in terms of a triangular mesh is described. The algorithm requires that each surface triangle can be handled independently, which implies the use of a 1D thermal model. Insolation of each triangle is determined by an optimized ray-triangle intersection search. Surface temperature is modeled with a spectral approach; imposing a quasi-periodic solution we replace heat conduction equation by the Helmholtz equation. Nonlinear boundary conditions are handled by an iterative, FFT based solver. The results resolve the question of the YORP effect in rotation rate independence on conductivity within the nonlinear 1D thermal model regardless of the accuracy issues and homogeneity assumptions. A seasonal YORP effect in attitude is revealed for objects moving on elliptic orbits when a nonlinear thermal model is used.