In this white paper, we present a cross-section of important scientific questions that remain partially or completely unanswered, ranging from Titan exosphere to the deep interior, and we detail which instrumentation and mission scenarios should be used to answer them. Our intention is to formulate the science goals for the next generation of planetary missions to Titan in order to prepare the future exploration of the moon. The ESA L-class mission concept that we propose is composed of a Titan orbiter and at least an in situ element (lake lander and/or drone(s)).

In order to examine how the terrestrial life emerged, a number of laboratory experiments have been conducted since the 1950s. Methane has been one of the key molecules in these studies. In earlier studies, strongly reducing gas mixtures containing methane and ammonia were mainly used to simulate possible reactions in primitive Earth atmosphere, and amino acids and other organic compounds were detected. Since the primitive Earth atmosphere was estimated to be less reducing, contribution of extraterrestrial organics to the origin of life is considered quite important. Extraterrestrial organic chemistry has been experimentally and theoretically studied intensively, including laboratory experiments simulating interstellar molecular reactions. Endogenous and exogenous organics should have been supplied to the primitive ocean. Now submarine hydrothermal systems are considered one of the plausible sites of generation of life. Experiments simulating submarine hydrothermal systems where methane played an important role are now intensively being conducted. We have recognized the importance of such studies on possible reactions in other solar system bodies to understand the origins of life. Titan and other icy bodies, where methane plays significant roles, are especially good targets to be studied. In the case of Titan, not only methane-containing atmospheres but also liquidospheres composed of methane and other hydrocarbons have been used in simulation experiments. This paper summarizes experiments simulating various terrestrial and extraterrestrial environments, and possible roles of methane in chemical evolution are discussed.

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.

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

We have investigated dissociative recombination (DR) of NH$^+$ with electrons using a merged beams configuration at the TSR heavy-ion storage ring located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. We present our measured absolute merged beams recombination rate coefficient for collision energies from 0 to 12 eV. From these data we have extracted a cross section which we have transformed to a plasma rate coefficient for the collisional plasma temperature range from $T_{\rm pl} = 10$ to $18000$ K. We show that the NH$^+$ DR rate coefficient data in current astrochemical models are underestimated by up to a factor of $\sim 9$. Our new data will result in predicted NH$^+$ abundances lower than calculated by present models. This is in agreement with the sensitivity limits of all observations attempting to detect NH$^+$ in interstellar clouds.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

Predictions of astrochemical models depend strongly on the reaction rate coefficients used in the simulations. We reviewed a number of key reactions for the chemistry of nitrogen-bearing species in the dense interstellar medium and proposed new reaction rate coefficients for those reactions. The details of the reviews are given in the form of a datasheet associated with each reaction. The new recommended rate coefficients are given with an uncertainty and a temperature range of validity and will be included in KIDA (http://kida.obs.u-bordeaux1.fr).

We have measured dissociative recombination of HCl+ with electrons using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. We present the measured absolute merged beams recombination rate coefficient for collision energies from 0 to 4.5 eV. We have also developed a new method for deriving the cross section from the measurements. Our approach does not suffer from approximations made by previously used methods. The cross section was transformed to a plasma rate coefficient for the electron temperature range from T=10 to 5000 K. We show that the previously used HCl+ DR data underestimate the plasma rate coefficient by a factor of 1.5 at T=10 K and overestimate it by a factor of 3.0 at T=300 K. We also find that the new data may partly explain existing discrepancies between observed abundances of chlorine-bearing molecules and their astrochemical models.

Absolute photodetachment cross sections of two anions of astrophysical importance CN- and C3N- were measured to be (1.18 +- (0.03)_stat (0.17)_sys) * 10^-17 cm^2 and (1.43 +- (0.14)_stat (0.37)_sys) * 10^-17 cm^2 respectively at the ultraviolet wavelength of 266 nm (4.66 eV). These relatively large values of the cross sections imply that photodetachment can play a major role in the destruction mechanisms of these anions particularly in photon-dominated regions. We have therefore carried out model calculations using the newly measured cross sections to investigate the abundance of these molecular anions in the cirumstellar envelope of the carbon-rich star IRC+10216. The model predicts the relative importance of the various mechanisms of formation and destruction of these species in different regions of the envelope. UV photodetachment was found to be the major destruction mechanism for both CN- and C3N- anions in those regions of the envelope, where they occur in peak abundance. It was also found that photodetachment plays a crucial role in the degradation of these anions throughout the circumstellar envelope.

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

We observed the 9-cm Lambda-doubling lines of CH along the dense filament of TMC-1. The CH column densities were compared with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation of the observed OH/CH abundance ratios. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0e-8, to the southeast where it reaches a value of 2.0e-8. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios of ~ 16 - 20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase provides a natural explanation for the observation of "early-type" molecules in this region.

We survey the current situation regarding chemical modelling of the synthesis of molecules in the interstellar medium. The present state of knowledge concerning the rate coefficients and their uncertainties for the major gas-phase processes -- ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative recombination -- is reviewed. Emphasis is placed on those reactions that have been identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances of the species observed in the interstellar medium. These sensitivity analyses have been carried out for gas-phase models of three representative, molecule-rich, astronomical sources: the cold dense molecular clouds TMC-1 and L134N, and the expanding circumstellar envelope IRC +10216. Our review has led to the proposal of new values and uncertainties for the rate coefficients of many of the key reactions. The impact of these new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar dust particles also influence the observed abundances of molecules in the interstellar medium. Their role is included in gas-grain, as distinct from gas-phase only, models. We review the methods for incorporating both accretion onto, and reactions on, the surfaces of grains in such models, as well as describing some recent experimental efforts to simulate and examine relevant processes in the laboratory. These efforts include experiments on the surface-catalysed recombination of hydrogen atoms, on chemical processing on and in the ices that are known to exist on the surface of interstellar grains, and on desorption processes, which may enable species formed on grains to return to the gas-phase.

Spectral line surveys are useful since they allow identification of new molecules and new lines in uniformly calibrated data sets. Nonetheless, large portions of the sub-millimetre spectral regime remain unexplored due to severe absorptions by H2O and O2 in the terrestrial atmosphere. The purpose of the measurements presented here is to cover wavelength regions at and around 0.55 mm -- regions largely unobservable from the ground. Using the Odin astronomy/aeronomy satellite, we performed the first spectral survey of the Orion KL molecular cloud core in the bands 486--492 and 541--576 GHz with rather uniform sensitivity (22--25 mK baseline noise). Odin's 1.1 m size telescope, equipped with four cryo-cooled tuneable mixers connected to broad band spectrometers, was used in a satellite position-switching mode. Two mixers simultaneously observed different 1.1 GHz bands using frequency steps of 0.5 GHz (25 hours each). An on-source integration time of 20 hours was achieved for most bands. The entire campaign consumed ~1100 orbits, each containing one hour of serviceable astro-observation. We identified 280 spectral lines from 38 known interstellar molecules (including isotopologues) having intensities in the range 80 to 0.05 K. An additional 64 weak lines remain unidentified. Apart from the ground state rotational 1(1,0)--1(0,1) transitions of ortho-H2O, H218O and H217O, the high energy 6(2,4)--7(1,7) line of para-H2O and the HDO(2,0,2--1,1,1) line have been observed, as well as the 1,0--0,1 lines from NH3 and its rare isotopologue 15NH3. We suggest assignments for some unidentified features, notably the new interstellar molecules ND and SH-. Severe blends have been detected in the line wings of the H218O, H217O and 13CO lines changing the true linewidths of the outflow emission.

We investigate the physical and chemical conditions in a typical star forming region, including an unbiased search for new molecules in a spectral region previously unobserved. Due to its proximity, the Orion KL region offers a unique laboratory of molecular astrophysics in a chemically rich, massive star forming region. Several ground-based spectral line surveys have been made, but due to the absorption by water and oxygen, the terrestrial atmosphere is completely opaque at frequencies around 487 and 557 GHz. To cover these frequencies we used the Odin satellite to perform a spectral line survey in the frequency ranges 486-492 GHz and 541-577 GHz, filling the gaps between previous spectral scans. Odin's high main beam efficiency and observations performed outside the atmosphere make our intensity scale very well determined. We observed 280 spectral lines from 38 molecules including isotopologues, and, in addition, 64 unidentified lines. The beam-averaged emission is dominated by CO, H2O, SO2, SO, 13CO and CH3OH. Species with the largest number of lines are CH3OH, (CH33)2O, SO2, 13CH3OH, CH3CN and NO. Six water lines are detected including the ground state rotational transition o-H2O, its isotopologues o-H218O and o-H217O, the Hot Core tracing p-H2O transition 6(2,4)-7(1,7), and the 2(0, 2)-1(1,1) transition of HDO. Other lines of special interest are the 1_0-0_0 transition of NH3 and its isotopologue 15NH3. Isotopologue abundance ratios of D/H, 12C/13C, 32S/34S, 34S/33S, and 18O/17O are estimated. The temperatures, column densities and abundances in the various subregions are estimated, and we find very high gas-phase abundances of H2O, NH3, SO2, SO, NO, and CH3OH. A comparison with the ice inventory of ISO sheds new light on the origin of the abundant gas-phase molecules.