Using atomistic simulations, we characterize the adsorption process of organic molecules on carbon nanoparticles, both of which have been reported to be abundant in the interstellar medium (ISM). It is found that the aromatic organics are adsorbed more readily than the aliphatic ones. This selectivity would favor the formation of polycyclic aromatic hydrocarbons (PAHs) or fullerene-like structures in the ISM due to structural similarity. It is also observed in our simulations that the molecules form a monolayer over the nanoparticle surface before stacking up in aggregates. This suggests a possible layer-by-layer formation process of onion-like nanostructures in the ISM. These findings reveal the possible role of carbon nanoparticles as selective catalysts that could provide reaction substrates for the formation of interstellar PAHs, high-fullerenes and soots from gas-phase molecules.

The origin of the martian methane is still poorly understood. A plausible explanation is that methane could have been produced either by hydrothermal alteration of basaltic crust or by serpentinization of ultramafic rocks producing hydrogen and reducing crustal carbon into methane. Once formed, methane storage on Mars is commonly associated with the presence of hidden clathrate reservoirs. Here, we alternatively suggest that chabazite and clinoptilolite, which belong to the family of zeolites, may form a plausible storage reservoir of methane in the martian subsurface. Because of the existence of many volcanic terrains, zeolites are expected to be widespread on Mars and their Global Equivalent Layer may range up to more than $\sim$1 km, according to the most optimistic estimates. If the martian methane present in chabazite and clinoptilolite is directly sourced from an abiotic source in the subsurface, the destabilization of a localized layer of a few millimeters per year may be sufficient to explain the current observations. The sporadic release of methane from these zeolites requires that they also remained isolated from the atmosphere during its evolution. The methane release over the ages could be due to several mechanisms such as impacts, seismic activity or erosion. If the methane outgassing from excavated chabazite and/or clinoptilolite prevails on Mars, then the presence of these zeolites around Gale Crater could explain the variation of methane level observed by Mars Science Laboratory.

We review the reservoirs of methane clathrates that may exist in the different bodies of the Solar System. Methane was formed in the interstellar medium prior to having been embedded in the protosolar nebula gas phase. This molecule was subsequently trapped in clathrates that formed from crystalline water ice during the cooling of the disk and incorporated in this form in the building blocks of comets, icy bodies, and giant planets. Methane clathrates may play an important role in the evolution of planetary atmospheres. On Earth, the production of methane in clathrates is essentially biological, and these compounds are mostly found in permafrost regions or in the sediments of continental shelves. On Mars, methane would more likely derive from hydrothermal reactions with olivine-rich material. If they do exist, martian methane clathrates would be stable only at depth in the cryosphere and sporadically release some methane into the atmosphere via mechanisms that remain to be determined.

Grand Canonical Monte Carlo simulations are used to reproduce the N$_2$/CO ratio ranging between 1.7 $\times$ 10$^{-3}$ and 1.6 $\times$ 10$^{-2}$ observed {\it in situ} in the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA mass spectrometer aboard the Rosetta spacecraft, assuming that this body has been agglomerated from clathrates in the protosolar nebula. Simulations are done using an elaborated interatomic potentials for investigating the temperature dependence of the trapping within a multiple guest clathrate formed from a gas mixture of CO and N$_2$ in proportions corresponding to those expected for the protosolar nebula. By assuming that 67P/Churyumov-Gerasimenko agglomerated from clathrates, our calculations suggest the cometary grains must have been formed at temperatures ranging between $\sim$31.8 and 69.9 K in the protosolar nebula to match the N$_2$/CO ratio measured by the ROSINA mass spectrometer. The presence of clathrates in Jupiter family comets could then explain the potential N$_2$ depletion (factor up to $\sim$87 compared to the protosolar value) measured in 67P/Churyumov-Gerasimenko.

Motivated by the possibility of comparing theoretical predictions of Lake Vostok's composition with future in situ measurements, we investigate the composition of clathrates that are expected to form in this environment from the air supplied to the lake by melting ice. In order to establish the best possible correlation between the lake water composition with that of air clathrates formed in situ, we use a statistical thermodynamic model based on the description of the guest-clathrate interaction by a spherically averaged Kihara potential with a nominal set of potential parameters. We determine the fugacities of the different volatiles present in the lake by defining a "pseudo" pure substance dissolved in water owning the average properties of the mixture and by using the Redlich-Kwong equation of state to mimic its thermodynamic behavior. Irrespective of the clathrate structure considered in our model, we find that xenon and krypton are strongly impoverished in the lake water (a ratio in the 0.04--0.1 range for xenon and a ratio in the $\sim$0.15--0.3 range for krypton), compared to their atmospheric abundances. Argon and methane are also found depleted in the Lake Vostok water by factors in the 0.5--0.95 and 0.3--0.5 ranges respectively, compared to their atmospheric abundances. On the other hand, the carbone dioxide abundance is found substantially enriched in the lake water compared to its atmospheric abundance (by a factor in the 1.6--5 range at 200 residence times). The comparison of our predictions of the CO$_2$ and CH$_4$ mole fractions in Lake Vostok with future in situ measurements will allow disentangling between the possible supply sources.

A striking feature of the atmosphere of Titan is that no heavy noble gases other than argon were detected by the Gas Chromatograph Mass Spectrometer (GCMS) aboard the Huygens probe during its descent to Titan's surface in January 2005. Here we provide an explanation of the mysterious absence or rarity of these noble gases in Titan's atmosphere: the thermodynamic conditions prevailing at the surface-atmosphere interface of the satellite allow the formation of multiple guest clathrates that preferentially store some species, including all heavy noble gases, over others. The clean water ice needed for formation of these clathrates could be delivered by successive episodes of cryovolcanic lavas that have been hypothesized to regularly cover the surface of Titan. The formation of clathrates in the porous lavas and their propensity for trapping Ar, Kr and Xe would progressively remove these species from the atmosphere of Titan over its history. In some circumstances, a global clathrate crust with an average thickness not exceeding a few meters could be sufficient on Titan for a complete removal of the heavy noble gases from the atmosphere.

Favored theories of giant planet formation center around two main paradigms, namely the core accretion model and the gravitational instability model. These two formation scenarios support the hypothesis that the giant planet metallicities should be higher or equal to that of the parent star. Meanwhile, spectra of the transiting hot Jupiter HD189733b suggest that carbon and oxygen abundances range from depleted to enriched with respect to the star. Here, using a model describing the formation sequence and composition of planetesimals in the protoplanetary disk, we determine the range of volatile abundances in the envelope of HD189733b that is consistent with the 20--80 Earth-masses of heavy elements estimated to be present in the planet's envelope. We then compare the inferred carbon and oxygen abundances to those retrieved from spectroscopy and we find a range of supersolar values that directly fit both spectra and internal structure models. In some cases, we find that the apparent contradiction between the subsolar elemental abundances and the {mass of heavy elements predicted in HD189733b by internal structure models} can be explained by the presence of large amounts of carbon molecules in the form of polycyclic aromatic hydrocarbons and soots in the upper layers of the envelope, as suggested by recent photochemical models. A diagnostic test that would confirm the presence of these compounds in the envelope is the detection of acetylene.

Examination of ambient thermodynamic conditions suggest that clathrate hydrates could exist in the martian permafrost, on the surface and in the interior of Titan, as well as in other icy satellites. Clathrate hydrates probably formed in a significant fraction of planetesimals in the solar system. Thus, these crystalline solids may have been accreted in comets, in the forming giant planets and in their surrounding satellite systems. In this work, we use a statistical thermodynamic model to investigate the composition of clathrate hydrates that may have formed in the primordial nebula. In our approach, we consider the formation sequence of the different ices occurring during the cooling of the nebula, a reasonable idealization of the process by which volatiles are trapped in planetesimals. We then determine the fractional occupancies of guests in each clathrate hydrate formed at given temperature. The major ingredient of our model is the description of the guest-clathrate hydrate interaction by a spherically averaged Kihara potential with a nominal set of parameters, most of which being fitted on experimental equilibrium data. Our model allows us to find that Kr, Ar and N$_2$ can be efficiently encaged in clathrate hydrates formed at temperatures higher than $\sim$ 48.5 K in the primitive nebula, instead of forming pure condensates below 30 K. However, we find at the same time that the determination of the relative abundances of guest species incorporated in these clathrate hydrates strongly depends on the choice of the parameters of the Kihara potential and also on the adopted size of cages. Indeed, testing different potential parameters, we have noted that even minor dispersions between the different existing sets can lead to non-negligible variations in the determination of the volatiles trapped in clathrate hydrates formed in the primordial nebula.

The difference between the measured atmospheric abundances of neon, argon, krypton and xenon for Venus, the Earth and Mars is striking. Because these abundances drop by at least two orders of magnitude as one moves outward from Venus to Mars, the study of the origin of this discrepancy is a key issue that must be explained if we are to fully understand the different delivery mechanisms of the volatiles accreted by the terrestrial planets. In this work, we aim to investigate whether it is possible to quantitatively explain the variation of the heavy noble gas abundances measured on Venus, the Earth and Mars, assuming that cometary bombardment was the main delivery mechanism of these noble gases to the terrestrial planets. To do so, we use recent dynamical simulations that allow the study of the impact fluxes of comets upon the terrestrial planets during the course of their formation and evolution. Assuming that the mass of noble gases delivered by comets is proportional to rate at which they collide with the terrestrial planets, we show that the krypton and xenon abundances in Venus and the Earth can be explained in a manner consistent with the hypothesis of cometary bombardment. In order to explain the krypton and xenon abundance differences between the Earth and Mars, we need to invoke the presence of large amounts of CO2-dominated clathrates in the Martian soil that would have efficiently sequestered these noble gases.

Recent observations have evidenced traces of methane CH4 heterogeneously distributed in the martian atmosphere. However, because the lifetime of CH4 in the atmosphere of Mars is estimated to be around 300-600 years on the basis of photochemistry, its release from a subsurface reservoir or an active primary source of methane have been invoked in the recent literature. Among the existing scenarios, it has been proposed that clathrate hydrates located in the near subsurface of Mars could be at the origin of the small quantities of the detected CH4. Here, we accurately determine the composition of these clathrate hydrates, as a function of temperature and gas phase composition, by using a hybrid statistical thermodynamic model based on experimental data. Compared to other recent works, our model allows us to calculate the composition of clathrate hydrates formed from a more plausible composition of the martian atmosphere by considering its main compounds, i.e. carbon dioxyde, nitrogen and argon, together with methane. Besides, because there is no low temperature restriction in our model, we are able to determine the composition of clathrate hydrates formed at temperatures corresponding to the extreme ones measured in the polar caps. Our results show that methane enriched clathrate hydrates could be stable in the subsurface of Mars only if a primitive CH4-rich atmosphere has existed or if a subsurface source of CH4 has been (or is still) present.

We describe a scenario of Titan's formation matching the constraints imposed by its current atmospheric composition. Assuming that the abundances of all elements, including oxygen, are solar in the outer nebula, we show that the icy planetesimals were agglomerated in the feeding zone of Saturn from a mixture of clathrates with multiple guest species, so-called stochiometric hydrates such as ammonia hydrate, and pure condensates. We also use a statistical thermodynamic approach to constrain the composition of multiple guest clathrates formed in the solar nebula. We then infer that krypton and xenon, that are expected to condense in the 20-30 K temperature range in the solar nebula, are trapped in clathrates at higher temperatures than 50 K. Once formed, these ices either were accreted by Saturn or remained embedded in its surrounding subnebula until they found their way into the regular satellites growing around Saturn. In order to explain the carbon monoxide and primordial argon deficiencies of Titan's atmosphere, we suggest that the satellite was formed from icy planetesimals initially produced in the solar nebula and that were partially devolatilized at a temperature not exceeding 50 K during their migration within Saturn's subnebula. The observed deficiencies of Titan's atmosphere in krypton and xenon could result from other processes that may have occurred both prior or after the completion of Titan. Thus, krypton and xenon may have been sequestrated in the form of XH3+ complexes in the solar nebula gas phase, causing the formation of noble gas-poor planetesimals ultimately accreted by Titan. Alternatively, krypton and xenon may have also been trapped efficiently in clathrates located on the satellite's surface or in its atmospheric haze.

In this paper, we use a statistical thermodynamic approach to quantify the efficiency with which clathrates on the surface of Titan trap noble gases. We consider different values of the Ar, Kr, Xe, CH4, C2H6 and N2 abundances in the gas phase that may be representative of Titan's early atmosphere. We discuss the effect of the various parameters that are chosen to represent the interactions between the guest species and the ice cage in our calculations. We also discuss the results of varying the size of the clathrate cages. We show that the trapping efficiency of clathrates is high enough to significantly decrease the atmospheric concentrations of Xe and, to a lesser extent, of Kr, irrespective of the initial gas phase composition, provided that these clathrates are abundant enough on the surface of Titan. In contrast, we find that Ar is poorly trapped in clathrates and, as a consequence, that the atmospheric abundance of argon should remain almost constant. We conclude that the mechanism of trapping noble gases via clathration can explain the deficiency in primordial Xe and Kr observed in Titan's atmosphere by Huygens, but that this mechanism is not sufficient to explain the deficiency in Ar.

We use a statistical thermodynamic approach to determine the composition of clathrate hydrates which may form from a multiple compound gas whose composition is similar to that of Titan's atmosphere. Assuming that noble gases are initially present in this gas phase, we calculate the ratios of xenon, krypton and argon to species trapped in clathrate hydrates. We find that these ratios calculated for xenon and krypton are several orders of magnitude higher than in the coexisting gas at temperature and pressure conditions close to those of Titan's present atmosphere at ground level. Furthermore we show that, by contrast, argon is poorly trapped in these ices. This trapping mechanism implies that the gas-phase is progressively depleted in xenon and krypton when the coexisting clathrate hydrates form whereas the initial abundance of argon remains almost constant. Our results are thus compatible with the deficiency of Titan's atmosphere in xenon and krypton measured by the {\it Huygens} probe during its descent on January 14, 2005. However, in order to interpret the subsolar abundance of primordial Ar also revealed by {\it Huygens}, other processes that occurred either during the formation of Titan or during its evolution must be also invoked.

The Two Micron All Sky Survey, along with the Stellar Population Synthesis Model of the Galaxy, developed in Besancon, is used to calculate the extinction distribution along different lines of sight. By combining many lines of sight, the large scale distribution of interstellar material can be deduced. The Galaxy model is used to provide the intrinsic colour of stars and their probable distances, so that the near infrared colour excess, and hence the extinction, may be calculated and its distance evaluated. Such a technique is dependent on the model used, however we are able to show that moderate changes in the model parameters result in insignificant changes in the predicted extinction. This technique has now been applied to over 64000 lines of sight, each separated by 15 arcmin, in the inner Galaxy (|l|<=100 deg, |b|<=10 deg). Using our extinction map, we have derived the main characteristics of the large scale structure of the dust distribution: scale height and warp of the ISM disc as well as the angle of the dust in the Galactic Bar. This resulting extinction map will be useful for studies of the inner Galaxy and its stellar populations.

The bulge is a region of the Galaxy which is of tremendous interest for understanding Galaxy formation. However, measuring photometry and kinematics in it raises several inherent issues, like high extinction in the visible and severe crowding. Here we attempt to estimate the problem of the visibility of the bulge at optical wavelengths, where large CCD mosaics allow to easily cover wide regions from the ground, and where future astrometric missions are planned. Assuming the Besancon Galaxy model and high resolution extinction maps, we estimate the stellar density as a function of longitude, latitude and apparent magnitude and we deduce the possibility of reaching and measuring bulge stars. The method is applied to three Gaia instruments, the BBP and MBP photometers, and the RVS spectrograph. We conclude that, while in the BBP most of the bulge will be accessible, in the MBP there will be a small but significant number of regions where bulge stars will be detected and accurately measured in crowded fields. Assuming that the RVS spectra may be extracted in moderately crowded fields, the bulge will be accessible in most regions apart from the strongly absorbed inner plane regions, because of high extinction, and in low extinction windows like the Baades's window where the crowding is too severe.

We attempt to study the characteristics of the different stellar populations present in the Galactic central region. A Monte Carlo method is used to simultaneously fit 11 thin disc and triaxial outer bulge density parameters on (Ks,J-Ks) star count data in almost 100 windows from the DENIS near infrared large scale survey at -8 deg < l< 12 deg and -4 deg < b < 4 deg. Various bulge density profiles and luminosity functions were tested uing a population synthesis scheme. The best models, selected by a maximum likelihood test, give the following description: the outer bulge is proxy, prolate, and oriented 10.6+/-3 deg with respect to the Sun-center direction. It seems that the main bulge population is not older than 10 Gyr, but this preliminary result needs further work to be confirmed. A significant central hole is found in the middle of the thin disc. We discuss these results in regard to previous findings and the scenario of bulge formation.

Since the Hipparcos mission and recent large scale surveys in the optical and the near-infrared, new constraints have been obtained on the structure and evolution history of the Milky Way. The population synthesis approach is a useful tool to interpret such data sets and to test scenarios of evolution of the Galaxy. We present here new constraints on evolution parameters obtained from the Besancon model of population synthesis and analysis of optical and near-infrared star counts. The Galactic potential is computed self-consistently, in agreement with Hipparcos results and the observed rotation curve. Constraints are posed on the outer bulge structure, the warped and flared disc, the thick disc and the spheroid populations. The model is tuned to produce reliable predictions in the visible and the near-infrared in wide photometric bands from U to K. Finally, we describe applications such as photometric and astrometric simulations and a new classification tool based on a Bayesian probability estimator, which could be used in the framework of Virtual Observatories. As examples, samples of simulated star counts at different wavelengths and directions are also given.

Garzon et al. (1997), Lopez-Corredoira et al. (1999) and Hammersley et al. (2000) have identified in TMGS and DENIS data a large excess of stars at l=27 deg and b=0 deg which might correspond to an in-plane bar. We compared near infrared CAIN star counts and simulations from the Besancon Model of Galaxy on 15 fields between 15 deg and 45 deg in longitude and -2 deg and 2 deg in latitude. Comparisons confirm the existence of an overdensity at longitudes lower than 27 deg which is inhomogeneous and decreases very steeply off the Galactic plane. The observed excess in the star distribution over the predicted density is even higher than 100%. Its distance from the sun is estimated to be lower than 6 kpc. If this overdensity corresponds to the stellar population of the bar, we estimate its half-length to 3.9 +/ -0.4 kpc and its angle from the Sun-center direction to 45 +/- 9 degrees.