• ### Internal Structure of Giant and Icy Planets: Importance of Heavy Elements and Mixing(1705.09320)

June 7, 2018 astro-ph.EP
In this chapter we summarize current knowledge of the internal structure of giant planets. We concentrate on the importance of heavy elements and their role in determining the planetary composition and internal structure, in planet formation, and during the planetary long-term evolution. We briefly discuss how internal structure models are derived, present the possible structures of the outer planets in the Solar System, and summarise giant planet formation and evolution. Finally, we introduce giant exoplanets and discuss how they can be used to better understand giant planets as a class of planetary objects.
• ### An Analysis of Stochastic Jovian Oscillation Excitation by Moist Convection(1802.05406)

Feb. 15, 2018 astro-ph.EP
Recent observations of Jupiter have suggested the existence of global oscillatory modes at millihertz frequencies, yet the source mechanism responsible for driving these modes is still unknown. However, the energies necessary to produce observable surface oscillations have been predicted. Here we investigate if moist convection in Jupiter's upper atmosphere can be responsible for driving the global oscillations and what moist convective energy requirements are necessary to achieve these theoretical mode energies and surface amplitudes. We begin by creating a one-dimensional moist convective cloud model and find that the available kinetic energy of the rising cloud column falls below theoretical estimates of oscillations energies. That is, mode excitation cannot occur with a single storm eruption. We then explore stochastic excitation scenarios of the oscillations by moist convective storms. We find that mode energies and amplitudes can reach theoretical estimates if the storm energy available to the modes is more than just kinetic. In order for the modes to be excited, we find that they require $5 \times 10^{27}$ to 10$^{28}$ erg per day. However, even for a large storm eruption each day, the available kinetic energy from the storms falls two orders of magnitude short of the required driving energy. Although our models may oversimplify the true complexity of the coupling between Jovian storms and global oscillations, our findings reveal that enough thermal energy is associated with moist convection to drive the modes, should it be available to them.
• ### A semi-empirical model for magnetic braking of solar-type stars(1709.03054)

Sept. 10, 2017 astro-ph.SR
We develop new angular momentum evolution models for stars with masses of $0.5$ to $1.6~\rm M_\odot$ and from the pre-main-sequence (\rm PMS) through the end of their main-sequence (\rm MS) lifetime. The parametric models include magnetic braking based on numerical simulations of magnetised stellar winds, mass loss rate prescription, core-envelope decoupling as well as disk locking phenomena. We have also accounted for recent developments in modelling dramatically weakened magnetic braking in stars more evolved than the Sun. We fit the free parameters in our model by comparing model predictions to rotational distributions of a number of stellar clusters as well as individual field stars. Our model reasonably successfully reproduces the rotational behaviour of stars during the \rm PMS phase to the zero-age main-sequence (\rm ZAMS) spin up, sudden \rm ZAMS spin down, and convergence of the rotation rates afterwards. We find that including core-envelope decoupling improves our models especially for low-mass stars at younger ages. In addition, by accounting for the almost complete suppression of magnetic braking at slow spin periods, we provide better fits to observations of stellar rotations compared to previous models.
• ### Revisiting the pre-main-sequence evolution of stars I. Importance of accretion efficiency and deuterium abundance(1702.07901)

Feb. 25, 2017 astro-ph.SR, astro-ph.EP
Recent theoretical work has shown that the pre-main-sequence (PMS) evolution of stars is much more complex than previously envisioned. Instead of the traditional steady, one-dimensional solution, accretion may be episodic and not necessarily symmetrical, thereby affecting the energy deposited inside the star and its interior structure. Given this new framework, we want to understand what controls the evolution of accreting stars. We use the MESA stellar evolution code with various sets of conditions. In particular, we account for the (unknown) efficiency of accretion in burying gravitational energy into the protostar through a parameter, $\xi$, and we vary the amount of deuterium present. We confirm the findings of previous works that the evolution changes significantly with the amount of energy that is lost during accretion. We find that deuterium burning also regulates the PMS evolution. In the low-entropy accretion scenario, the evolutionary tracks in the H-R diagram are significantly different from the classical tracks and are sensitive to the deuterium content. A comparison of theoretical evolutionary tracks and observations allows us to exclude some cold accretion models ($\xi\sim 0$) with low deuterium abundances. We confirm that the luminosity spread seen in clusters can be explained by models with a somewhat inefficient injection of accretion heat. The resulting evolutionary tracks then become sensitive to the accretion heat efficiency, initial core entropy, and deuterium content. In this context, we predict that clusters with a higher D/H ratio should have less scatter in luminosity than clusters with a smaller D/H. Future work on this issue should include radiation-hydrodynamic simulations to determine the efficiency of accretion heating and further observations to investigate the deuterium content in star-forming regions. (abbrev.)
• ### Jupiter internal structure: the effect of different equations of state(1609.05460)

Nov. 13, 2016 astro-ph.EP
Heavy elements, even though its smaller constituent, are crucial to understand Jupiter formation history. Interior models are used to determine the amount of heavy elements in Jupiter interior, nevertheless this range is still subject to degeneracies due to uncertainties in the equations of state. Prior to Juno mission data arrival, we present Jupiter optimized calculations exploring the effect of different model parameters in the determination of Jupiter's core and heavy element's mass. We perform comparisons between equations of state published recently. The interior model of Jupiter is calculated from the equations of hydrostatic equilibrium, mass and energy conservation, and energy transport. The mass of the core and heavy elements is adjusted to match Jupiter's observational constrains radius and gravitational moments. We show that the determination of Jupiter interior structure is tied to the estimation of its gravitational moments and the accuracy of equations of state of hydrogen, helium and heavy elements. The location of the region where Helium rain occurs as well as its timescale are important to determine the distribution of heavy elements and helium in the interior of Jupiter. We show that differences find when modeling Jupiter's interior with recent EOS are more likely due to differences in the internal energy and entropy calculation. The consequent changes in the thermal profile lead to different estimations of the mass of the core and heavy elements, explaining differences in recently published Jupiter interior models. Our results help clarify differences find in Jupiter interior models and will help the interpretation of upcoming Juno data.
• ### Formation of dust-rich planetesimals from sublimated pebbles inside of the snow line(1610.09643)

Nov. 5, 2016 astro-ph.EP
Content: For up to a few millions of years, pebbles must provide a quasi-steady inflow of solids from the outer parts of protoplanetary disks to their inner regions. Aims: We wish to understand how a significant fraction of the pebbles grows into planetesimals instead of being lost to the host star. Methods:We examined analytically how the inward flow of pebbles is affected by the snow line and under which conditions dust-rich (rocky) planetesimals form. When calculating the inward drift of solids that is due to gas drag, we included the back-reaction of the gas to the motion of the solids. Results: We show that in low-viscosity protoplanetary disks (with a monotonous surface density similar to that of the minimum-mass solar nebula), the flow of pebbles does not usually reach the required surface density to form planetesimals by streaming instability. We show, however, that if the pebble-to-gas-mass flux exceeds a critical value, no steady solution can be found for the solid-to-gas ratio. This is particularly important for low-viscosity disks (alpha < 10^(-3)) where we show that inside of the snow line, silicate-dust grains ejected from sublimating pebbles can accumulate, eventually leading to the formation of dust-rich planetesimals directly by gravitational instability. Conclusions: This formation of dust-rich planetesimals may occur for extended periods of time, while the snow line sweeps from several au to inside of 1 au. The rock-to-ice ratio may thus be globally significantly higher in planetesimals and planets than in the central star.
• ### Condensation-inhibited convection in hydrogen-rich atmospheres: Stability against double-diffusive processes and thermal profiles for Jupiter, Saturn, Uranus, and Neptune(1610.05506)

Oct. 18, 2016 astro-ph.SR, astro-ph.EP
In an atmosphere, a cloud condensation region is characterized by a strong vertical gradient in the abundance of the related condensing species. On Earth, the ensuing gradient of mean molecular weight has relatively few dynamical consequences because N$_2$ is heavier than water vapor, so that only the release of latent heat significantly impacts convection. On the contrary, in an hydrogen dominated atmosphere (e.g. giant planets), all condensing species are significantly heavier than the background gas. This can stabilize the atmosphere against convection near a cloud deck if the enrichment in the given species exceeds a critical threshold. This raises two questions. What is transporting energy in such a stabilized layer, and how affected can the thermal profile of giant planets be? To answer these questions, we first carry out a linear analysis of the convective and double-diffusive instabilities in a condensable medium showing that an efficient condensation can suppress double-diffusive convection. This suggests that a stable radiative layer can form near a cloud condensation level, leading to an increase in the temperature of the deep adiabat. Then, we investigate the impact of the condensation of the most abundant species---water---with a steady-state atmosphere model. Compared to standard models, the temperature increase can reach several hundred degrees at the quenching depth of key chemical tracers. Overall, this effect could have many implications for our understanding of the dynamical and chemical state of the atmosphere of giant planets, for their future observations (with Juno for example), and for their internal evolution.
• ### Evolution of Protoplanetary Discs with Magnetically Driven Disc Winds(1609.00437)

Oct. 1, 2016 astro-ph.SR, astro-ph.EP
Aims: We investigate the evolution of protoplanetary discs (PPDs hereafter) with magnetically driven disc winds and viscous heating. Methods: We consider an initially massive disc with ~0.1 Msun to track the evolution from the early stage of PPDs. We solve the time evolution of surface density and temperature by taking into account viscous heating and the loss of the mass and the angular momentum by the disc winds within the framework of a standard alpha model for accretion discs. Our model parameters, turbulent viscosity, disc wind mass loss, and disc wind torque, which are adopted from local magnetohydrodynamical simulations and constrained by the global energetics of the gravitational accretion, largely depends on the physical condition of PPDs, particularly on the evolution of the vertical magnetic flux in weakly ionized PPDs. Results: Although there are still uncertainties concerning the evolution of the vertical magnetic flux remaining, surface densities show a large variety, depending on the combination of these three parameters, some of which are very different from the surface density expected from the standard accretion. When a PPD is in a "wind-driven accretion" state with the preserved vertical magnetic field, the radial dependence of the surface density can be positive in the inner region <1-10 au. The mass accretion rates are consistent with observations, even in the very low level of magnetohydrodynamical turbulence. Such a positive radial slope of the surface density gives a great impact on planet formation because (i)it inhibits the inward drift or even results in the outward drift of pebble/boulder-sized solid bodies, and (ii) it also makes the inward type-I migration of proto-planets slower or even reversed. Conclusions: The variety of our calculated PPDs should yield a wide variety of exoplanet systems.
• ### The Origin and Evolution of Saturn, with Exoplanet Perspective(1606.04510)

Sept. 8, 2016 astro-ph.EP
Saturn formed beyond the snow line in the primordial solar nebula that made it possible for it to accrete a large mass. Disk instability and core accretion models have been proposed for Saturn's formation, but core accretion is favored on the basis of its volatile abundances, internal structure, hydrodynamic models, chemical characteristics of protoplanetary disk, etc. The observed frequency, properties and models of exoplanets provide additional supporting evidence for core accretion. The heavy elements with mass greater than 4He make up the core of Saturn, but are presently poorly constrained, except for carbon. The C/H ratio is super-solar, and twice that in Jupiter. The enrichment of carbon and other heavy elements in Saturn and Jupiter requires special delivery mechanisms for volatiles to these planets. In this chapter we will review our current understanding of the origin and evolution of Saturn and its atmosphere, using a multi-faceted approach that combines diverse sets of observations on volatile composition and abundances, relevant properties of the moons and rings, comparison with the other gas giant planet, Jupiter, analogies to the extrasolar giant planets, as well as pertinent theoretical models.
• ### The radial dependence of pebble accretion rates: A source of diversity in planetary systems I. Analytical formulation(1604.01291)

April 14, 2016 astro-ph.EP
Context. The classical "planetesimal" accretion scenario for the formation of planets has recently evolved with the idea that "pebbles", centimeter- to meter-sized icy grains migrating in protoplanetary disks, can control planetesimal and/or planetary growth. Aims. We investigate how pebble accretion depends on disk properties and affects the formation of planetary systems Methods. We construct analytical models of pebble accretion onto planetary embryos that consistently account for the mass and orbital evolution of the pebble flow and reflect disk structure. Results. We derive simple formulas for pebble accretion rates in the so-called "settling" regime for planetary embryos with more than 100 km in size. For relatively smaller embryos or in outer disk regions, the accretion mode is 3D, meaning that the thickness of the pebble flow must be taken into account, and resulting in an accretion rate that is independent of the embryo mass. For larger embryos or in inner regions, the accretion is in a 2D mode, i.e., the pebble disk may be considered to be infinitely thin. We show that the radial dependence of the pebble accretion rata is different (even the sign of the power-law exponent changes) for different disk conditions such as the disk heating source (viscous heating or stellar irradiation), the drag law (Stokes or Epstein, and weak or strong coupling), and in the 2D or 3D accretion modes. We also discuss the effect of the sublimation/destruction of icy pebbles inside the snow line. Conclusions. Pebble accretion easily produces a large diversity of planetary systems. In other words, to correctly infer the results of planet formation through pebble accretion, detailed prescriptions of disk evolution and pebble growth, sublimation/destruction and migration are required.
• ### New constraints on Saturn's interior from Cassini astrometric data(1510.05870)

Nov. 26, 2015 astro-ph.EP
Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturn's tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturn's Love number, $k_2=0.390 \pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturn's interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planet's interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 \pm 0.74) \times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rhea's frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturn's moons might themselves experience large tidal dissipation.
• ### Suppression of type I migration by disk winds(1510.06010)

Oct. 20, 2015 astro-ph.EP
Planets less massive than Saturn tend to rapidly migrate inward in protoplanetary disks. This is the so-called type I migration. Simulations attempting to reproduce the observed properties of exoplanets show that type I migration needs to be significantly reduced over a wide region of the disk for a long time. However, the mechanism capable of suppressing type I migration over a wide region has remained elusive. The recently found turbulence-driven disk winds offer new possibilities. We investigate the effects of disk winds on the disk profile and type I migration for a range of parameters that describe the strength of disk winds. We also examine the in situ formation of close-in super-Earths in disks that evolve through disk winds. The disk profile, which is regulated by viscous diffusion and disk winds, was derived by solving the diffusion equation. We carried out a number of simulations and plot here migration maps that indicate the type I migration rate. We also performed N-body simulations of the formation of close-in super-Earths from a population of planetesimals and planetary embryos. We define a key parameter, Kw, which determines the ratio of strengths between the viscous diffusion and disk winds. For a wide range of Kw, the type I migration rate is presented in migration maps. These maps show that type I migration is suppressed over the whole close-in region when the effects of disk winds are relatively strong (Kw < 100). From the results of N-body simulations, we see that type I migration is significantly slowed down assuming Kw = 40. We also show that the results of N-body simulations match statistical orbital distributions of close-in super-Earths.
• ### Thermalizing a telescope in Antarctica: Analysis of ASTEP observations(1506.06009)

The installation and operation of a telescope in Antarctica represent particular challenges, in particular the requirement to operate at extremely cold temperatures, to cope with rapid temperature fluctuations and to prevent frosting. Heating of electronic subsystems is a necessity, but solutions must be found to avoid the turbulence induced by temperature fluctua- tions on the optical paths. ASTEP 400 is a 40 cm Newton telescope installed at the Concordia station, Dome C since 2010 for photometric observations of fields of stars and their exoplanets. While the telescope is designed to spread star light on several pixels to maximize photometric stability, we show that it is nonetheless sensitive to the extreme variations of the seeing at the ground level (between about 0.1 and 5 arcsec) and to temperature fluctuations between --30 degrees C and --80 degrees C. We analyze both day-time and night-time observations and obtain the magnitude of the seeing caused by the mirrors, dome and camera. The most important effect arises from the heating of the primary mirror which gives rise to a mirror seeing of 0.23 arcsec K--1 . We propose solutions to mitigate these effects.
• ### A reassessment of the in situ formation of close-in super-Earths(1504.03237)

April 13, 2015 astro-ph.EP
A large fraction of stars host one or multiple close-in super-Earth planets. There is an active debate about whether these planets formed in situ or at greater distances from the central star and migrated to their current position. It has been shown that part of their observed properties (e.g., eccentricity distribution) can be reproduced by N-body simulations of in situ formation starting with a population of protoplanets of high masses and neglecting the effects of the disk gas. We plan to reassess the in situ formation of close-in super-Earths through more complete simulations. We performed N-body simulations of a population of small planetary embryos and planetesimals that include the effects of disk-planet interactions (e.g., eccentricity damping, type I migration). In addition, we also consider the accretion of a primitive atmosphere from a protoplanetary disk. We find that planetary embryos grow very quickly well before the gas dispersal, and thus undergo rapid inward migration, which means that one cannot neglect the effects of a gas disk when considering the in-situ formation of close-in super-Earths. Owing to their rapid inward migration, super-Earths reach a compact configuration near the disk's inner edge whose distribution of orbital parameters matches the observed close-in super-Earths population poorly. On the other hand, simulations including eccentricity damping, but no type I migration, reproduce the observed distributions better. Including the accretion of an atmosphere does not help reproduce the bulk architecture of observations. Interestingly, we find that the massive embryos can migrate inside the disk edge while capturing only a moderately massive hydrogen/helium atmosphere. By this process they avoid becoming giant planets. The bulk of close-in super-Earths cannot form in situ, unless type I migration is suppressed in the entire disk inside 1 AU.
• The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.
• ### Giant Planets(1405.3752)

Jan. 20, 2015 cond-mat.other, astro-ph.EP
We review the interior structure and evolution of Jupiter, Saturn, Uranus and Neptune, and giant exoplanets with particular emphasis on constraining their global composition. Compared to the first edition of this review, we provide a new discussion of the atmospheric compositions of the solar system giant planets, we discuss the discovery of oscillations of Jupiter and Saturn, the significant improvements in our understanding of the behavior of material at high pressures and the consequences for interior and evolution models. We place the giant planets in our Solar System in context with the trends seen for exoplanets.
• ### Seismology of Giant Planets(1411.1740)

Nov. 11, 2014 astro-ph.EP
Seismology applied to giant planets could drastically change our understanding of their deep interiors, as it has happened with the Earth, the Sun, and many main-sequence and evolved stars. The study of giant planets' composition is important for understanding both the mechanisms enabling their formation and the origins of planetary systems, in particular our own. Unfortunately, its determination is complicated by the fact that their interior is thought not to be homogeneous, so that spectroscopic determinations of atmospheric abundances are probably not representative of the planet as a whole. Instead, the determination of their composition and structure must rely on indirect measurements and interior models. Giant planets are mostly fluid and convective, which makes their seismology much closer to that of solar-like stars than that of terrestrial planets. Hence, helioseismology techniques naturally transfer to giant planets. In addition, two alternative methods can be used: photometry of the solar light reflected by planetary atmospheres, and ring seismology in the specific case of Saturn. The current decade has been promising thanks to the detection of Jupiter's acoustic oscillations with the ground-based imaging-spectrometer SYMPA and indirect detection of Saturn's f-modes in its rings by the NASA Cassini orbiter. This has motivated new projects of ground-based and space-borne instruments that are under development. The purpose of this paper is to support the chapter Seismology of Giant Planets which will appear in May 2015 in the book Extraterrestrial Seismology, published in Cambridge University Press. Here, we review the key questions about jovian interiors as well as the most recent observational results, and refer the reader to the book chapter above-mentioned for further details.
• ### A non-grey analytical model for irradiated atmospheres. II: Analytical vs. numerical solutions(1311.6322)

Nov. 7, 2014 astro-ph.SR, astro-ph.EP
• ### Evolution of Exoplanets and their Parent Stars(1409.7477)

Sept. 26, 2014 astro-ph.SR, astro-ph.EP
Studying exoplanets with their parent stars is crucial to understand their population, formation and history. We review some of the key questions regarding their evolution with particular emphasis on giant gaseous exoplanets orbiting close to solar-type stars. For masses above that of Saturn, transiting exoplanets have large radii indicative of the presence of a massive hydrogen-helium envelope. Theoretical models show that this envelope progressively cools and contracts with a rate of energy loss inversely proportional to the planetary age. The combined measurement of planetary mass, radius and a constraint on the (stellar) age enables a global determination of the amount of heavy elements present in the planet interior. The comparison with stellar metallicity shows a correlation between the two, indicating that accretion played a crucial role in the formation of planets. The dynamical evolution of exoplanets also depends on the properties of the central star. We show that the lack of massive giant planets and brown dwarfs in close orbit around G-dwarfs and their presence around F-dwarfs are probably tied to the different properties of dissipation in the stellar interiors. Both the evolution and the composition of stars and planets are intimately linked.
• ### On the filtering and processing of dust by planetesimals 1. Derivation of collision probabilities for non-drifting planetesimals(1409.7328)

Sept. 26, 2014 astro-ph.SR, astro-ph.EP
Context. Circumstellar disks are known to contain a significant mass in dust ranging from micron to centimeter size. Meteorites are evidence that individual grains of those sizes were collected and assembled into planetesimals in the young solar system. Aims. We assess the efficiency of dust collection of a swarm of non-drifting planetesimals {\rev with radii ranging from 1 to $10^3$\,km and beyond. Methods. We calculate the collision probability of dust drifting in the disk due to gas drag by planetesimal accounting for several regimes depending on the size of the planetesimal, dust, and orbital distance: the geometric, Safronov, settling, and three-body regimes. We also include a hydrodynamical regime to account for the fact that small grains tend to be carried by the gas flow around planetesimals. Results. We provide expressions for the collision probability of dust by planetesimals and for the filtering efficiency by a swarm of planetesimals. For standard turbulence conditions (i.e., a turbulence parameter $\alpha=10^{-2}$), filtering is found to be inefficient, meaning that when crossing a minimum-mass solar nebula (MMSN) belt of planetesimals extending between 0.1 AU and 35 AU most dust particles are eventually accreted by the central star rather than colliding with planetesimals. However, if the disk is weakly turbulent ($\alpha=10^{-4}$) filtering becomes efficient in two regimes: (i) when planetesimals are all smaller than about 10 km in size, in which case collisions mostly take place in the geometric regime; and (ii) when planetary embryos larger than about 1000 km in size dominate the distribution, have a scale height smaller than one tenth of the gas scale height, and dust is of millimeter size or larger in which case most collisions take place in the settling regime. These two regimes have very different properties: we find that the local filtering efficiency $x_{filter,MMSN}$ scales with $r^{-7/4}$ (where $r$ is the orbital distance) in the geometric regime, but with $r^{-1/4}$ to $r^{1/4}$ in the settling regime. This implies that the filtering of dust by small planetesimals should occur close to the central star and with a short spread in orbital distances. On the other hand, the filtering by embryos in the settling regime is expected to be more gradual and determined by the extent of the disk of embryos. Dust particles much smaller than millimeter size tend only to be captured by the smallest planetesimals because they otherwise move on gas streamlines and their collisions take place in the hydrodynamical regime. Conclusions. Our results hint at an inside-out formation of planetesimals in the infant solar system because small planetesimals in the geometrical limit can filter dust much more efficiently close to the central star. However, even a fully-formed belt of planetesimals such as the MMSN only marginally captures inward-drifting dust and this seems to imply that dust in the protosolar disk has been filtered by planetesimals even smaller than 1 km (not included in this study) or that it has been assembled into planetesimals by other mechanisms (e.g., orderly growth, capture into vortexes). Further refinement of our work concerns, among other things: a quantitative description of the transition region between the hydro and settling regimes; an assessment of the role of disk turbulence for collisions, in particular in the hydro regime; and the coupling of our model to a planetesimal formation model.
• ### A non-grey analytical model for irradiated atmospheres. I: Derivation(1311.6597)

Nov. 26, 2013 astro-ph.SR, astro-ph.EP
Context. Semi-grey atmospheric models (with one opacity for the visible and one opacity for the infrared) are useful to understand the global structure of irradiated atmospheres, their dynamics and the interior structure and evolution of planets, brown dwarfs and stars. But when compared to direct numerical radiative transfer calculations for irradiated exoplanets, these models systematically overestimate the temperatures at low optical depth, independently of the opacity parameters. We wish to understand why semi-grey models fail at low optical depths, and provide a more accurate approximation to the atmospheric structure by accounting for the variable opacity in the infrared. Our analytical irradiated non-grey model is found to provide a range of temperatures that is consistent with that obtained by numerical calculations. We find that even for slightly non-grey thermal opacities the temperature structure differs significantly from previous semi-grey models. For small values of beta (expected when lines are dominant), we find that the non-grey effects are confined to low-optical depths. However, for beta larger than 0.5 (appropriate in the presence of bands with a wavelength-dependence smaller or comparable with the width of the Planck function), we find that the temperature structure is affected even down to infrared optical depths unity and deeper as a result of the so-called blanketing effect. The expressions that we derive may be used to provide a proper functional form for algorithms that invert the atmospheric properties from spectral information. Because a full atmospheric structure can be calculated directly, these expressions should be useful for simulations of the dynamics of these atmospheres and of the thermal evolution of the planets. Finally, they should be used to test full radiative transfer models and improve their convergence.
• ### Bulk Composition of GJ 1214b and other sub-Neptune exoplanets(1305.2629)

June 14, 2013 astro-ph.EP
GJ1214b stands out among the detected low-mass exoplanets, because it is, so far, the only one amenable to transmission spectroscopy. Up to date there is no consensus about the composition of its envelope although most studies suggest a high molecular weight atmosphere. In particular, it is unclear if hydrogen and helium are present or if the atmosphere is water dominated. Here, we present results on the composition of the envelope obtained by using an internal structure and evolutionary model to fit the mass and radius data. By examining all possible mixtures of water and H/He, with the corresponding opacities, we find that the bulk amount of H/He of GJ1214b is at most 7% by mass. In general, we find the radius of warm sub-Neptunes to be most sensitive to the amount of H/He. We note that all (Kepler-11b,c,d,f, Kepler-18b, Kepler-20b, 55Cnc-e, Kepler-36c and Kepler-68b) but two (Kepler-11e and Kepler-30b) of the discovered low-mass planets so far have less than 10% H/He. In fact, Kepler-11e and Kepler-30b have 10-18% and 5-15% bulk H/He. Conversely, little can be determined about the H2O or rocky content of sub-Neptune planets. We find that although a 100% water composition fits the data for GJ1214b, based on formation constraints the presence of heavier refractory material on this planet is expected, and hence, so is a component lighter than water required. A robust determination by transmission spectroscopy of the composition of the upper atmosphere of GJ1214b will help determine the extent of compositional segregation between the atmosphere and envelope.
• ### CoRoT: harvest of the exoplanet program(1306.0578)

June 3, 2013 astro-ph.EP
One of the objectives of the CoRoT mission is the search for transiting extrasolar planets using high-precision photometry, and the accurate characterization of their fundamental parameters. The CoRoT satellite consecutively observes crowded stellar fields since February 2007, in high-cadence precise photometry; periodic eclipses are detected and analysed in the stellar light curves. Then complementary observations using ground-based facilities allows establishing the nature of the transiting body and its mass. CoRoT has acquired more than 163,000 light curves and detected about 500 planet candidates. A fraction of them (5%) are confirmed planets whose masses are independently measured. Main highlights of the CoRoT discoveries are: i) the variety of internal structures in close-in giant planets, ii) the characterisation of the first known transiting rocky planet, CoRoT-7 b, iii) multiple constraints on the formation, evolution, role of tides in planetary systems.
• ### Interior Models of Saturn: Including the Uncertainties in Shape and Rotation(1302.6690)

Feb. 27, 2013 physics.geo-ph, astro-ph.EP
The accurate determination of Saturn's gravitational coefficients by Cassini could provide tighter constrains on Saturn's internal structure. Also, occultation measurements provide important information on the planetary shape which is often not considered in structure models. In this paper we explore how wind velocities and internal rotation affect the planetary shape and the constraints on Saturn's interior. We show that within the geodetic approach (Lindal et al., 1985, ApJ, 90, 1136) the derived physical shape is insensitive to the assumed deep rotation. Saturn's re-derived equatorial and polar radii at 100 mbar are found to be 54,445 $\pm$10 km and 60,365$\pm$10 km, respectively. To determine Saturn's interior we use {\it 1 D} three-layer hydrostatic structure models, and present two approaches to include the constraints on the shape. These approaches, however, result in only small differences in Saturn's derived composition. The uncertainty in Saturn's rotation period is more significant: with Voyager's 10h39mns period, the derived mass of heavy elements in the envelope is 0-7 M$_{\oplus}$. With a rotation period of 10h32mns, this value becomes $<4$ $M_{\oplus}$, below the minimum mass inferred from spectroscopic measurements. Saturn's core mass is found to depend strongly on the pressure at which helium phase separation occurs, and is estimated to be 5-20 M$_{\oplus}$. Lower core masses are possible if the separation occurs deeper than 4 Mbars. We suggest that the analysis of Cassini's radio occultation measurements is crucial to test shape models and could lead to constraints on Saturn's rotation profile and departures from hydrostatic equilibrium.
• ### Spitzer Infrared Observations and Independent Validation of the Transiting Super-Earth CoRoT-7b(1110.5336)

Oct. 24, 2011 astro-ph.EP
The detection and characterization of the first transiting super-Earth, CoRoT-7 b, has required an unprecedented effort in terms of telescope time and analysis. Although the star does display a radial velocity signal at the period of the planet, this has been difficult to disentangle from the intrinsic stellar variability, and pinning down the velocity amplitude has been very challenging. As a result, the precise value of the mass of the planet - and even the extent to which it can be considered to be confirmed - have been debated in the recent literature, with six mass measurements published so far based on the same spectroscopic observations, ranging from about 2 to 8 Earth masses. Here we report on an independent validation of the planet discovery, using one of the fundamental properties of a transit signal: its achromaticity. We observed four transits of CoRoT-7 b with Spitzer, in order to determine whether the depth of the transit signal in the near-infrared is consistent with that observed in the CoRoT bandpass, as expected for a planet. We detected the transit and found an average depth of 0.426 {\pm} 0.115 mmag at 4.5 {\mu}m, which is in good agreement with the depth of 0.350 {\pm} 0.011 mmag found by CoRoT. These observations place important constraints on the kinds of astrophysical false positives that could mimic the signal. Combining this with additional constraints reported earlier, we performed an exhaustive exploration of possible blends scenarios for CoRoT-7 b using the BLENDER technique. We are able to rule out the vast majority of false positives, and the remaining ones are found to be much less likely than a true transiting planet. We thus validate CoRoT-7 b as a bona-fide planet with a very high degree of confidence, independently of any radial-velocity information. Our Spitzer observations have additionally allowed us to significantly improve the ephemeris of the planet.