We present here the first release of the open-source python package ExoTETHyS, which aims to provide a stand-alone set of tools for modeling spectro-photometric observations of the transiting exoplanets. In particular, we describe: (1) a new calculator of stellar limb-darkening coefficients that outperforms the existing software by one order of magnitude in terms of light-curve model accuracy, i.e., down to <10 parts per million (ppm); (2) an exact transit light-curve generator based on the entire stellar intensity profile rather than limb-darkening coefficients. New tools will be added in later releases to model various effects in exoplanetary transits and eclipsing binaries. ExoTETHyS is a reference package for high-precision exoplanet atmospheric spectroscopy with the upcoming JWST and ARIEL missions.

"Hot super-Earths" (or "Mini-Neptunes") between 1 and 4 times Earth's size with period shorter than 100 days orbit 30-50\% of Sun-like type stars. Their orbital configuration -- measured as the period ratio distribution of adjacent planets in multi-planet systems -- is a strong constraint for formation models. Here we use N-body simulations with synthetic forces from an underlying evolving gaseous disk to model the formation and long-term dynamical evolution of super-Earth systems. While the gas disk is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disk. These resonant chains are far more compact than the observed super-Earth systems. Once the gas dissipates resonant chains may become dynamically unstable. They undergo a phase of giant impacts that spreads the systems out. Disk turbulence has no measurable effect on the outcome. Our simulations match observations if a small fraction of resonant chains remain stable, while most super-Earths undergo a late dynamical instability. Our statistical analysis restricts the contribution of stable systems to less than $25\%$. Our results also suggest that the large fraction of observed single planet systems does not necessarily imply any dichotomy in the architecture of planetary systems. Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in $\sim 5\%$ of the cases. This leads to a mystery: in our simulations only 50-60\% of resonant chains became unstable whereas at least 75\% (and probably 90-95\%) must be unstable to match observations.

Planetary embryos embedded in gaseous protoplanetary disks undergo Type I orbital migration. Migration can be inward or outward depending on the local disk properties but, in general, only planets more massive than several $M_\oplus$ can migrate outward. Here we propose that an embryo's migration history determines whether it becomes a hot super-Earth or the core of a giant planet. Systems of hot super-Earths (or mini-Neptunes) form when embryos migrate inward and pile up at the inner edge of the disk. Giant planet cores form when inward-migrating embryos become massive enough to switch direction and migrate outward. We present simulations of this process using a modified N-body code, starting from a swarm of planetary embryos. Systems of hot super-Earths form in resonant chains with the innermost planet at or interior to the disk inner edge. Resonant chains are disrupted by late dynamical instabilities triggered by the dispersal of the gaseous disk. Giant planet cores migrate outward toward zero-torque zones, which move inward and eventually disappear as the disk disperses. Giant planet cores migrate inward with these zones and are stranded at ~1-5 AU. Our model reproduces several properties of the observed extra-solar planet populations. The frequency of giant planet cores increases strongly when the mass in solids is increased, consistent with the observed giant exoplanet - stellar metallicity correlation. The frequency of hot super-Earths is not a function of stellar metallicity, also in agreement with observations. Our simulations can reproduce the broad characteristics of the observed super-Earth population.

Rate constants for the C(3P) + CH3OH reaction have been measured in a continuous supersonic flow reactor over the range 50 K to 296 K. C(3P) was created by the in-situ pulsed laser photolysis of CBr4, a multiphoton process which also produced some C(1D), allowing us to investigate simultaneously the low temperature kinetics of the C(1D) + CH3OH reaction. C(1D) atoms were followed by an indirect chemiluminescent tracer method in the presence of excess CH3OH. C(3P) atoms were detected by the same chemiluminescence technique and also by direct vacuum ultra-violet laser induced fluorescence (VUV LIF). Secondary measurements of product H(2S) atom formation have been undertaken allowing absolute H atom yields to be obtained by comparison with a suitable reference reaction. In parallel, statistical calculations have been performed based on ab-initio calculations of the complexes, adducts and transition states (TSs) relevant to the title reaction. By comparison with the experimental H atom yields, the preferred reaction pathways could be determined, placing important constraints on the statistical calculations. The experimental and theoretical work are in excellent agreement, predicting a negative temperature dependence of the rate constant increasing from 2.2 x 10-11 cm3 molecule-1 s-1 at 296 K to 20.0 x 10-11 cm3 molecule-1 s-1 at 50 K. CH3 and HCO are found to be the major products under our experimental conditions. As this reaction is not considered in current astrochemical networks, its influence on interstellar methanol abundances is tested using a model of dense interstellar clouds.

It has been proposed that the observed systems of hot super-Earths formed in situ from high-mass disks. By fitting a disk profile to the entire population of Kepler planet candidates, Chiang & Laughlin (2013) constructed a "minimum-mass extrasolar nebula" with surface density profile Sigma r^-1.6. Here we use multiple-planet systems to show that it is inconsistent to assume a universal disk profile. Systems with 3-6 low-mass planets (or planet candidates) produce a diversity of minimum-mass disks with surface density profiles ranging from Sigma r^-3.2 to Sigma r^0.5 (5th-95th percentile). By simulating the transit detection of populations of synthetic planetary systems designed to match the properties of observed super-Earth systems, we show that a universal disk profile is statistically excluded at high confidence. Rather, the underlying distribution of minimum-mass disks is characterized by a broad range of surface density slopes. Models of gaseous disks can only explain a narrow range of slopes (roughly between r^0 and r^-1.5). Yet accretion of terrestrial planets in a gas-free environment preserves the initial radial distribution of building blocks. The known systems of hot super-Earths must therefore not represent the structure of their parent gas disks and can not have predominantly formed in situ. We instead interpret the diversity of disk slopes as the imprint of a process that re-arranged the solids relative to the gas in the inner parts of protoplanetary disks. A plausible mechanism is inward type 1 migration of Mars- to Earth-mass planetary embryos, perhaps followed by a final assembly phase.

Earth-mass bodies are expected to undergo Type I migration directed either inward or outward depending on the thermodynamical state of the protoplanetary disc. Zones of convergent migration exist where the Type I torque cancels out. We study the evolution of multiple protoplanets of a few Earth masses embedded in a non-isothermal protoplanetary disc. The protoplanets are located in the vicinity of a convergence zone located at the transition between two different opacity regimes. Inside the convergence zone, Type I migration is directed outward and outside the zone migration is directed inward. We used a grid-based hydrodynamical code that includes radiative effects. We performed simulations varying the initial number of embryos and tested the effect of including stochastic forces to mimic the effects resulting from turbulence. We also performed N-body runs calibrated on hydrodynamical calculations to follow the evolution on Myr timescales. For a small number of initial embryos (N = 5-7) and in the absence of stochastic forcing, the population of protoplanets migrates convergently toward the zero-torque radius and forms a stable resonant chain that protects embryos from close encounters. In systems with a larger initial number of embryos, or in which stochastic forces were included, these resonant configurations are disrupted. This in turn leads to the growth of larger cores via a phase of giant impacts, after which the system settles to a new stable resonant configuration. Giant planets cores with masses of 10 Earth masses formed in about half of the simulations with initial protoplanet masses of m_p = 3 Earth masses but in only 15% of simulations with m_p = 1 Earth mass. This suggests that if ~2-3 Earth mass protoplanets can form in less than ~1 Myr, convergent migration and giant collisions can grow giant planet cores at Type I migration convergence zones.

Using N-body simulations with planet-disk interactions, we present a mechanism capable of forming compact systems of hot super Earths such as Kepler 11. Recent studies show that outward migration is common in the inner parts of radiative disks. However we show that two processes naturally tip the balance in favor of inward migration. First the corotation torque is too weak to generate outward migration for planetary embryos less massive than $4\mearth$. Second, system of multiple embryos generate sustained non-zero eccentricities that damp the corotation torque and again favor inward migration. Migration and accretion of planetary embryos in realistic disks naturally produce super Earths in resonant chains near the disk inner edge. Their compact configuration is similar to the observed systems.

Earth-mass planets embedded in gaseous protoplanetary disks undergo Type I orbital migration. In radiative disks an additional component of the corotation torque scaling with the entropy gradient across the horseshoe region can counteract the general inward migration, Type I migration can then be directed either inward or outward depending on the local disk properties. Thus, special locations exist in the disk toward which planets migrate in a convergent way. Here we present N-body simulations of the convergent migration of systems of low-mass (M=1-10 m_earth) planets. We show that planets do not actually converge in convergence zones. Rather, they become trapped in chains of mean motion resonances (MMRs). This causes the planets' eccentricities to increase to high enough values to affect the structure of the horseshoe region and weaken the positive corotation torque. The zero-torque equilibrium point of the resonant chain of planets is determined by the sum of the attenuated corotation torques and unattenuated differential Lindblad torques acting on each planet. The effective convergence zone is shifted inward. Systems with several planets can experience stochastic migration as a whole due to continuous perturbations from planets entering and leaving resonances.

Eclipse exoplanet spectroscopy has yielded detection of H_2O, CH_4, CO_2 and CO in the atmosphere of hot jupiters and neptunes. About 40 large terrestrial planets are announced or confirmed, two of which are transiting, and another deemed habitable. Hence the potential for eclipse spectroscopy of terrestrial planets with James Webb Space Telescope (JWST) has become an active field of study. We explore the parameter space (type of stars, planet orbital periods and types, and instruments/wavelengths) in terms of the signal-to-noise ratio (S/N) achievable on the detection of spectroscopic features. We use analytic formula and model data for both the astrophysical scene and the instrument, to plot S/N contour maps, while indicating how the S/N scales with the fixed parameters. We systematically compare stellar photon noise-only figures with ones including detailed instrumental and zodiacal noises. Likelihood of occurring targets is based both on model and catalog star population of the solar neighborhood. The 9.6 micron ozone band is detectable (S/N = 3) with JWST, for a warm super-earth 6.7 pc away, using ~2% of the 5-year nominal mission time (summing observations, M4V and lighter host star for primary eclipses, M5V for secondary). If every star up to this mass limit and distance were to host a habitable planet, there should be statistically ~1 eclipsing case. Investigation of systematic noises in the co-addition of 5 years worth-, tens of days separated-, hours-long observations is critical, complemented by dedicated characterisation of the instruments, currently in integration phase. The census of nearby transiting habitable planets must be complete before the beginning of science operations.