scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: Earth and Planetary Astrophysics in 2016"


Journal ArticleDOI
TL;DR: The first high-resolution sub-mm survey of both dust and gas for a large population of protoplanetary disks is presented in this article, where the authors use ALMA to survey 89 disks around stars with a few Martian masses (0.2--0.4~$M_{\oplus}$).
Abstract: We present the first high-resolution sub-mm survey of both dust and gas for a large population of protoplanetary disks. Characterizing fundamental properties of protoplanetary disks on a statistical level is critical to understanding how disks evolve into the diverse exoplanet population. We use ALMA to survey 89 protoplanetary disks around stars with $M_{\ast}>0.1~M_{\odot}$ in the young (1--3~Myr), nearby (150--200~pc) Lupus complex. Our observations cover the 890~$\mu$m continuum and the $^{13}$CO and C$^{18}$O 3--2 lines. We use the sub-mm continuum to constrain $M_{\rm dust}$ to a few Martian masses (0.2--0.4~$M_{\oplus}$) and the CO isotopologue lines to constrain $M_{\rm gas}$ to roughly a Jupiter mass (assuming ISM-like $\rm {[CO]/[H_2]}$ abundance). Of 89 sources, we detect 62 in continuum, 36 in $^{13}$CO, and 11 in C$^{18}$O at $>3\sigma$ significance. Stacking individually undetected sources limits their average dust mass to $\lesssim6$ Lunar masses (0.03~$M_{\oplus}$), indicating rapid evolution once disk clearing begins. We find a positive correlation between $M_{\rm dust}$ and $M_{\ast}$, and present the first evidence for a positive correlation between $M_{\rm gas}$ and $M_{\ast}$, which may explain the dependence of giant planet frequency on host star mass. The mean dust mass in Lupus is 3$\times$ higher than in Upper Sco, while the dust mass distributions in Lupus and Taurus are statistically indistinguishable. Most detected disks have $M_{\rm gas}\lesssim1~M_{\rm Jup}$ and gas-to-dust ratios $<100$, assuming ISM-like $\rm {[CO]/[H_2]}$ abundance; unless CO is very depleted, the inferred gas depletion indicates that planet formation is well underway by a few Myr and may explain the unexpected prevalence of super-Earths in the exoplanet population.

508 citations


Journal ArticleDOI
TL;DR: In this paper, the authors present astrophysical false positive probability calculations for every Kepler Object of Interest (KOI) and demonstrate the first large-scale demonstration of a fully automated transiting planet validation procedure.
Abstract: We present astrophysical false positive probability calculations for every Kepler Object of Interest (KOI)---the first large-scale demonstration of a fully automated transiting planet validation procedure. Out of 7056 KOIs, we determine that 1935 have probabilities <1% to be astrophysical false positives, and thus may be considered validated planets. 1284 of these have not yet been validated or confirmed by other methods. In addition, we identify 428 KOIs likely to be false positives that have not yet been identified as such, though some of these may be a result of unidentified transit timing variations. A side product of these calculations is full stellar property posterior samplings for every host star, modeled as single, binary, and triple systems. These calculations use 'vespa', a publicly available Python package able to be easily applied to any transiting exoplanet candidate.

413 citations


Journal ArticleDOI
TL;DR: In this paper, an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region is presented.
Abstract: The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

350 citations


Journal ArticleDOI
TL;DR: In this paper, the authors show that the orbits of distant Kuiper belt objects cluster not only in argument of perihelion, but also in physical space, and they demonstrate that the positions and orbital planes of the objects are tightly confined and that such a clustering has only a probability of 0.007% to be due to chance.
Abstract: Recent analyses have shown that distant orbits within the scattered disk population of the Kuiper belt exhibit an unexpected clustering in their respective arguments of perihelion. While several hypotheses have been put forward to explain this alignment, to date, a theoretical model that can successfully account for the observations remains elusive. In this work we show that the orbits of distant Kuiper belt objects cluster not only in argument of perihelion, but also in physical space. We demonstrate that the perihelion positions and orbital planes of the objects are tightly confined and that such a clustering has only a probability of 0.007% to be due to chance, thus requiring a dynamical origin. We find that the observed orbital alignment can be maintained by a distant eccentric planet with mass greater than ~10 Earth masses, whose orbit lies in approximately the same plane as those of the distant Kuiper belt objects, but whose perihelion is 180 degrees away from the perihelia of the minor bodies. In addition to accounting for the observed orbital alignment, the existence of such a planet naturally explains the presence of high perihelion Sedna-like objects, as well as the known collection of high semimajor axis objects with inclinations between 60 and 150 degrees whose origin was previously unclear. Continued analysis of both distant and highly inclined outer solar system objects provides the opportunity for testing our hypothesis as well as further constraining the orbital elements and mass of the distant planet.

268 citations


Journal ArticleDOI
TL;DR: In this article, a review of recent progress in high-contrast imaging with particular emphasis on observational results, discoveries near and below the deuterium-burning limit, and a practical overview of large-scale surveys and dedicated instruments is presented.
Abstract: High-contrast adaptive optics imaging is a powerful technique to probe the architectures of planetary systems from the outside-in and survey the atmospheres of self-luminous giant planets. Direct imaging has rapidly matured over the past decade and especially the last few years with the advent of high-order adaptive optics systems, dedicated planet-finding instruments with specialized coronagraphs, and innovative observing and post-processing strategies to suppress speckle noise. This review summarizes recent progress in high-contrast imaging with particular emphasis on observational results, discoveries near and below the deuterium-burning limit, and a practical overview of large-scale surveys and dedicated instruments. I conclude with a statistical meta-analysis of deep imaging surveys in the literature. Based on observations of 384 unique and single young ($\approx$5--300~Myr) stars spanning stellar masses between 0.1--3.0~\Msun, the overall occurrence rate of 5--13~\Mjup \ companions at orbital distances of 30--300~AU is 0.6$^{+0.7}_{-0.5}$\% assuming hot-start evolutionary models. The most massive giant planets regularly accessible to direct imaging are about as rare as hot Jupiters are around Sun-like stars. Dividing this sample into individual stellar mass bins does not reveal any statistically-significant trend in planet frequency with host mass: giant planets are found around 2.8$^{+3.7}_{-2.3}$\% of BA stars, $<$4.1\% of FGK stars, and $<$3.9\% of M dwarfs. Looking forward, extreme adaptive optics systems and the next generation of ground- and space-based telescopes with smaller inner working angles and deeper detection limits will increase the pace of discovery to ultimately map the demographics, composition, evolution, and origin of planets spanning a broad range of masses and ages.

261 citations


Journal ArticleDOI
Xue-Ning Bai1
TL;DR: In this paper, the role of wind-driven accretion and wind mass loss in the formation of PPDs is discussed. But the authors focus on the global evolution of the PPD.
Abstract: A global evolution picture of protoplanetary disks (PPDs) is key to understanding almost every aspect of planet formation, where standard alpha-disk models have been constantly employed for its simplicity. In the mean time, disk mass loss has been conventionally attributed to photoevaporation, which controls disk dispersal. However, a paradigm shift towards accretion driven by magnetized disk winds has been realized in the recent years, thanks to studies of non-ideal magneto-hydrodynamic effects in PPDs. I present a framework of global PPD evolution aiming to incorporate these advances, highlighting the role of wind-driven accretion and wind mass loss. Disk evolution is found to be largely dominated by wind-driven processes, and viscous spreading is suppressed. The timescale of disk evolution is controlled primarily by the amount of external magnetic flux threading the disks, and how rapidly the disk loses the flux. Rapid disk dispersal can be achieved if the disk is able to hold most of its magnetic flux during the evolution. In addition, because wind launching requires sufficient level of ionization at disk surface (mainly via external far-UV radiation), wind kinematics is also affected by far-UV penetration depth and disk geometry. For typical disk lifetime of a few Myrs, the disk loses approximately the same amount of mass through the wind as through accretion onto the protostar, and most of the wind mass loss proceeds from the outer disk via a slow wind. Fractional wind mass loss increases with increasing disk lifetime. Significant wind mass loss likely substantially enhances the dust to gas mass ratio, and promotes planet formation.

243 citations


Journal ArticleDOI
TL;DR: In this article, a chain of models, linking the formation of a planet to its observable present-day spectrum, is presented, including the planet's formation and migration, its long-term thermodynamic evolution, a variety of disk chemistry models, a non-gray atmospheric model, and a radiometric model to obtain simulated spectroscopic observations.
Abstract: The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a "chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with JWST and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) Envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O 1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed.

213 citations


Journal ArticleDOI
TL;DR: In this paper, the authors connect models of dust evolution and planetesimal formation to identify regions of protoplanetary discs that are favorable for the formation of kilometre-sized bodies and the first planetary embryos.
Abstract: The consistency of planet formation models suffers from the disconnection between the regime of small and large bodies. This is primarily caused by so-called growth barriers: the direct growth of larger bodies is halted at centimetre-sized objects and particular conditions are required for the formation of larger, gravitationally bound planetesimals. We aim to connect models of dust evolution and planetesimal formation to identify regions of protoplanetary discs that are favourable for the formation of kilometre-sized bodies and the first planetary embryos. We combine semi-analytical models of viscous protoplanetary disc evolution, dust growth and drift including backreaction of the dust particles on the gas, and planetesimal formation via the streaming instability into one numerical code. We investigate how planetesimal formation is affected by the mass of the protoplanetary disc, its initial dust content, and the stickiness of dust aggregates. We find that the dust growth and drift leads to a global redistribution of solids. The pile-up of pebbles in the inner disc provides local conditions where the streaming instability is effective. Planetesimals form in an annulus with its inner edge lying between 0.3 AU and 1 AU and its width ranging from 0.3 AU to 3 AU. The resulting surface density of planetesimals follows a radial profile that is much steeper than the initial disc profile. These results support formation of terrestrial planets in the solar system from a narrow annulus of planetesimals, which reproduces their peculiar mass ratios.

213 citations


Journal ArticleDOI
TL;DR: In this article, scattered light images of the TW Hya disk performed with SPHERE in PDI mode at 063, 079, 124 and 162 micron were presented.
Abstract: We present scattered light images of the TW Hya disk performed with SPHERE in PDI mode at 063, 079, 124 and 162 micron We also present H2/H3-band ADI observations Three distinct radial depressions in the polarized intensity distribution are seen, around 85, 21, and 6~au The overall intensity distribution has a high degree of azimuthal symmetry; the disk is somewhat brighter than average towards the South and darker towards the North-West The ADI observations yielded no signifiant detection of point sources in the disk Our observations have a linear spatial resolution of 1 to 2au, similar to that of recent ALMA dust continuum observations The sub-micron sized dust grains that dominate the light scattering in the disk surface are strongly coupled to the gas We created a radiative transfer disk model with self-consistent temperature and vertical structure iteration and including grain size-dependent dust settling This method may provide independent constraints on the gas distribution at higher spatial resolution than is feasible with ALMA gas line observations We find that the gas surface density in the "gaps" is reduced by 50% to 80% relative to an unperturbed model Should embedded planets be responsible for carving the gaps then their masses are at most a few 10 Mearth The observed gaps are wider, with shallower flanks, than expected for planet-disk interaction with such low-mass planets If forming planetary bodies have undergone collapse and are in the "detachted phase" then they may be directly observable with future facilities such as METIS at the E-ELT

202 citations


Journal ArticleDOI
TL;DR: In this article, the authors report the results of the statistical analysis of planetary signals discovered in MOA-II microlensing survey alert system events from 2007 to 2012, finding significant planetary signals in 23 of the 1474 alert events that are well characterized by the MOA II survey data alone, and combine this analysis with the previous analyses of Gould et al. and Cassan et al., bringing the total sample to 30 planets.
Abstract: We report the results of the statistical analysis of planetary signals discovered in MOA-II microlensing survey alert system events from 2007 to 2012. We determine the survey sensitivity as a function of planet-star mass ratio, $q$, and projected planet-star separation, $s$, in Einstein radius units. We find that the mass ratio function is not a single power-law, but has a change in slope at $q \sim 10^{-4}$, corresponding to $\sim 20 M_{\oplus}$ for the median host star mass of $\sim 0.6 M_{\odot}$. We find significant planetary signals in 23 of the 1474 alert events that are well characterized by the MOA-II survey data alone. Data from other groups are used only to characterize planetary signals that have been identified in the MOA data alone. The distribution of mass ratios and separations of the planets found in our sample are well fit by a broken power-law model of the form $dN_{\rm pl}/(d{\rm log} q\ d{\rm log} s) = A (q/q_{\rm br})^n s^m \, {\rm dex}^{-2}$ for $q > q_{\rm br}$ and $dN_{\rm pl}/(d{\rm log} q\ d{\rm log} s) = A (q/q_{\rm br})^p s^m \, {\rm dex}^{-2}$ for $q < q_{\rm br}$, where $q_{\rm br}$ is the mass ratio of the break. We also combine this analysis with the previous analyses of Gould et al. and Cassan et al., bringing the total sample to 30 planets. This combined analysis yields $A = 0.61^{+0.21}_{-0.16}$, $n =-0.93\pm 0.13$, $m = 0.49_{-0.49}^{+0.47}$ and $p = 0.6^{+0.5}_{-0.4}$ for $q_{\rm br}\equiv 1.7\times 10^{-4}$. The unbroken power law model is disfavored with a $p$-value of 0.0022, which corresponds to a Bayes factor of 27 favoring the broken power-law model. These results imply that cold Neptunes are likely to be the most common type of planets beyond the snow line.

196 citations


Journal ArticleDOI
TL;DR: The nature of the early Martian climate is one of the major unanswered questions of planetary science as discussed by the authors, but a new wave of orbital and in situ observations and improvements in climate modeling have led to significant advances over the last decade.
Abstract: The nature of the early Martian climate is one of the major unanswered questions of planetary science. Key challenges remain, but a new wave of orbital and in situ observations and improvements in climate modeling have led to significant advances over the last decade. Multiple lines of geologic evidence now point to an episodically warm surface during the late Noachian and early Hesperian periods 3-4 Ga. The low solar flux received by Mars in its first billion years and inefficiency of plausible greenhouse gases such as CO2 means that the steady-state early Martian climate was likely cold. A denser CO2 atmosphere would have caused adiabatic cooling of the surface and hence migration of water ice to the higher altitude equatorial and southern regions of the planet. Transient warming caused melting of snow and ice deposits and a temporarily active hydrological cycle, leading to erosion of the valley networks and other fluvial features. Precise details of the warming mechanisms remain unclear, but impacts, volcanism and orbital forcing all likely played an important role. The lack of evidence for glaciation across much of Mars' ancient terrain suggests the late Noachian surface water inventory was not sufficient to sustain a northern ocean. While mainly inhospitable on the surface, early Mars may nonetheless have presented significant opportunities for the development of microbial life.

Journal ArticleDOI
TL;DR: In this paper, high precision HST WFC3 observations of the dayside emission spectrum of the hot Jupiter HD209458b were used to detect a water band in absorption at 6.2 sigma confidence.
Abstract: The nature of the vertical thermal structure of hot Jupiter atmospheres is one of the key questions raised by the characterization of transiting exoplanets over the last decade. There have been claims that many hot Jupiter's exhibit vertical profiles with increasing temperature with decreasing pressure in the infrared photosphere that leads to the reversal of molecular absorption bands into emission features (an inversion). However, these claims have been based on broadband photometry rather than the unambiguous identification of emission features with spectroscopy, and the chemical species that could cause the thermal inversions by absorbing stellar irradiation at high altitudes have not been identified despite extensive theoretical and observational effort. Here we present high precision HST WFC3 observations of the dayside emission spectrum of the hot Jupiter HD209458b; the first exoplanet suggested to have a thermal inversion. Our observations resolve a water band in absorption at 6.2 sigma confidence. When combined with Spitzer photometry the data are indicative of a monotonically decreasing temperature with pressure over the range 1-0.001 bar at 7.7 sigma confidence. We test the robustness of our results by exploring a variety of model assumptions including the temperature profile parameterization, presence of a cloud, and choice of Spitzer data reduction. We also introduce a new analysis method, "chemical retrieval-on-retrieval", to determine the elemental abundances from the spectrally retrieved mixing ratios with thermochemical self-consistency and find plausible abundances consistent with solar metallicity (0.06 - 10 x solar) and carbon-to-oxygen ratios less than unity. This work suggests that high-precision spectrophotometric results are required to robustly infer thermal structures and compositions of extra-solar planet atmospheres.

Journal ArticleDOI
TL;DR: In this article, the authors used the Spitzer Space Telescope to obtain a total of three phase curve observations of WASP-43b (P=0.813 days) at 3.6 and 4.5 microns.
Abstract: Previous measurements of heat redistribution efficiency (the ability to transport energy from a planet's highly irradiated dayside to its eternally dark nightside) show considerable variation between exoplanets. Theoretical models predict a positive correlation between heat redistribution efficiency and temperature for tidally locked planets; however, recent HST WASP-43b spectroscopic phase curve results are inconsistent with current predictions. Using the Spitzer Space Telescope, we obtained a total of three phase curve observations of WASP-43b (P=0.813 days) at 3.6 and 4.5 microns. The first 3.6 micron visit exhibits spurious nightside emission that requires invoking unphysical conditions in our cloud-free atmospheric retrievals. The two other visits exhibit strong day-night contrasts that are consistent with the HST data. To reconcile the departure from theoretical predictions, WASP-43b would need to have a high-altitude, nightside cloud/haze layer blocking its thermal emission. Clouds/hazes could be produced within the planet's cool, nearly retrograde mid-latitude flows before dispersing across its nightside at high altitudes. Since mid-latitude flows only materialize in fast-rotating ($\lesssim1$ day) planets, this may explain an observed trend connecting measured day-night contrast with planet rotation rate that matches all current Spitzer phase curve results. Combining independent planetary emission measurements from multiple phases, we obtain a precise dayside hemisphere H2O abundance ($2.5\times 10^{-5} - 1.1\times 10^{-4}$ at 1$\sigma$ confidence) and, assuming chemical equilibrium and a scaled solar abundance pattern, we derive a corresponding metallicity estimate that is consistent with being solar (0.4 -- 1.7). Using the retrieved global CO+CO2 abundance under the same assumptions, we estimate a comparable metallicity of 0.3 - 1.7$\times$ solar. (Abridged)

Journal ArticleDOI
TL;DR: In this paper, a 3D Global Climate Model (GCM) was used to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect.
Abstract: Radial velocity monitoring has found the signature of a $M \sin i = 1.3$~M$_\oplus$ planet located within the Habitable Zone (HZ) of Proxima Centauri \citep{Anglada16}. Despite a hotter past and an active host star the planet Proxima~b could have retained enough volatiles to sustain surface habitability \citep{Ribas2016}. Here we use a 3D Global Climate Model (GCM) to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect. We find that a broad range of atmospheric compositions allow surface liquid water. On a tidally-locked planet with sufficient surface water inventory, liquid water is always present, at least in the substellar region. With a non-synchronous rotation, this requires a minimum greenhouse warming ($\sim$10~mbar of CO$_2$ and 1~bar of N$_2$). If the planet is dryer, $\sim$0.5~bar/1.5~bars of CO$_2$ (respectively for asynchronous/synchronous rotation) suffice to prevent the trapping of any arbitrary small water inventory into polar/nightside ice caps. We produce reflection/emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible by direct imaging with forthcoming large telescopes. The angular separation of $7 \lambda/D$ at 1~$\mu$m (with the E-ELT) and a contrast of $\sim$10$^{-7}$ will enable high-resolution spectroscopy and the search for molecular signatures, including H$_2$O, O$_2$, and CO$_2$. The observation of thermal phase curves can be attempted with JWST, thanks to a contrast of $2\times10^{-5}$ at 10~$\mu$m. Proxima~b will also be an exceptional target for future IR interferometers. Within a decade it will be possible to image Proxima~b and possibly determine whether this exoplanet's surface is habitable.

Journal ArticleDOI
TL;DR: In this paper, the authors examined the location and shape of the brown-dwarf desert in both the period-mass and period-radius planes, using the currently available large samples of planets.
Abstract: A few studies have reported a significant dearth of exoplanets with Neptune mass and radius with orbital periods below $2$--$4$ d. This cannot be explained by observational biases because many Neptunian planets with longer orbital periods have been detected. The existence of this desert is similar to the appearance of the so-called brown-dwarf desert that suggests different formation mechanisms of planets and stellar companions with short orbital periods. Similarly, the Neptunian desert might indicate different mechanisms of formation and evolution for hot Jupiters and short-period super-Earths. We here follow a previous study and examine the location and shape of the desert in both the period-mass and period-radius planes, using the currently available large samples of planets. The desert in the period-mass plane has a relatively sharp upper edge, with a planetary mass that is inversely proportional to the planetary orbital period, while the lower, somewhat blurred, boundary is located along masses that are apparently linearly proportional to the period. The desert in the period-radius plane of the transiting planets is less clear. It seems as if the radius along the upper boundary is inversely proportional to the period to the power of one-third, while the lower boundary shows a radius that is proportional to the period to the power of two-thirds. The combination of the two upper bounds of the desert, in the period-mass and period-radius planes, yields a planetary mass-radius relation of $R_{\rm p}/R_{\rm Jup}\simeq (1.2\pm0.3)(M_{\rm p}/M_{\rm Jup})^{0.27\pm0.11}\,$ for $ 0.1\lesssim M_{\rm p}/M_{\rm Jup}\lesssim 1$. The derived shape of the desert, which might extend up to periods of $5$--$10$ d, could shed some light on the formation and evolution of close-in planets.

Journal ArticleDOI
TL;DR: In this paper, the authors investigate dust structures in the (sub)mm continuum emission of four disks (TW Hya, HL Tau, HD 163296 and DM Tau) with the highest spatial resolution ALMA observations to date.
Abstract: An unsolved problem in step-wise core-accretion planet formation is that rapid radial drift in gas-rich protoplanetary disks should drive mm-/meter-sized particles inward to the central star before large bodies can form. One promising solution is to confine solids within small scale structures. Here we investigate dust structures in the (sub)mm continuum emission of four disks (TW Hya, HL Tau, HD 163296 and DM Tau), a sample of disks with the highest spatial resolution ALMA observations to date. We retrieve the surface brightness distributions using synthesized images and fitting visibilities with analytical functions. We find that the continuum emission of the four disks is ~axi-symmetric but rich in 10-30AU-sized radial structures, possibly due to physical gaps, surface density enhancements or localized dust opacity variations within the disks. These results suggest that small scale axi-symmetric dust structures are likely to be common, as a result of ubiquitous processes in disk evolution and planet formation. Compared with recent spatially resolved observations of CO snowlines in these same disks, all four systems show enhanced continuum emission from regions just beyond the CO condensation fronts, potentially suggesting a causal relationship between dust growth/trapping and snowlines.

Journal ArticleDOI
TL;DR: In this article, the authors provide self-consistent atmospheric models and synthetic JWST observations for prime exoplanet targets in order to identify spectral regions of interest and estimate the number of transits needed to distinguish between model setups.
Abstract: The James Webb Space Telescope will enable astronomers to obtain exoplanet spectra of unprecedented precision. Especially the MIRI instrument may shed light on the nature of the cloud particles obscuring planetary transmission spectra in the optical and near-infrared. We provide self-consistent atmospheric models and synthetic JWST observations for prime exoplanet targets in order to identify spectral regions of interest and estimate the number of transits needed to distinguish between model setups. We select targets which span a wide range in planetary temperature and surface gravity, ranging from super-Earths to giant planets, and have a high expected SNR. For all targets we vary the enrichment, C/O ratio, presence of optical absorbers (TiO/VO) and cloud treatment. We calculate atmospheric structures and emission and transmission spectra for all targets and use a radiometric model to obtain simulated observations. We analyze JWST's ability to distinguish between various scenarios. We find that in very cloudy planets such as GJ 1214b less than 10 transits with NIRSpec may be enough to reveal molecular features. Further, the presence of small silicate grains in atmospheres of hot Jupiters may be detectable with a single JWST MIRI transit. For a more detailed characterization of such particles less than 10 transits are necessary. Finally, we find that some of the hottest hot Jupiters are well fitted by models which neglect the redistribution of the insolation and harbor inversions, and that 1-4 eclipse measurements with NIRSpec are needed to distinguish between the inversion models. Wet thus demonstrate the capabilities of JWST for solving some of the most intriguing puzzles in current exoplanet atmospheric research. Further, by publishing all models calculated for this study we enable the community to carry out similar or retrieval analyses for all planets included in our target list.

Journal ArticleDOI
TL;DR: In this paper, the authors investigate the combined role of radiative braking and stellar wind interactions using numerical simulations with the EVaporating Exoplanet code (EVE) and derive atmospheric and stellar properties through a direct comparison of simulated and observed spectra.
Abstract: The warm Neptune GJ436b was observed with HST/STIS at three different epochs in the stellar Ly-alpha line, showing deep, repeated transits caused by a giant exosphere of neutral hydrogen. The low radiation pressure from the M-dwarf host star was shown to play a major role in the dynamics of the escaping gas. Yet by itself it cannot explain the time-variable spectral features detected in each transit. Here we investigate the combined role of radiative braking and stellar wind interactions using numerical simulations with the EVaporating Exoplanet code (EVE) and we derive atmospheric and stellar properties through the direct comparison of simulated and observed spectra. Our simulations match the last two epochs well. The observed sharp early ingresses come from the abrasion of the planetary coma by the stellar wind. Spectra observed during the transit can be produced by a dual exosphere of planetary neutrals (escaped from the upper atmosphere of the planet) and neutralized protons (created by charge-exchange with the stellar wind). We find similar properties at both epochs for the planetary escape rate (2.5x10$^{8}$ g/s), the stellar photoionization rate (2x10$^{-5}$ /s), the stellar wind bulk velocity (85 km/s), and its kinetic dispersion velocity (10 km/s). We find high velocities for the escaping gas (50-60 km/s) that may indicate MHD waves that dissipate in the upper atmosphere and drive the planetary outflow. In the last epoch the high density of the stellar wind (3x10$^{3}$ /cm3) led to the formation of an exospheric tail mainly composed of neutralized protons. The observations of GJ436 b allow for the first time to clearly separate the contributions of radiation pressure and stellar wind and to probe the regions of the exosphere shaped by each mechanism.

Journal ArticleDOI
TL;DR: In this paper, a numerical algorithm for stiff mutual drag force was used to simulate small particles with significantly higher resolutions and longer simulation times than in previous investigations, and it was shown that particles of dimensionless stopping time (i.e., 10-2 and 10-3 ) can concentrate themselves via the streaming instability at a solid abundance of a few percent.
Abstract: Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform simulations of small particles with significantly higher resolutions and longer simulation times than in previous investigations. We find that particles of dimensionless stopping time $\tau_\mathrm{s} = 10^{-2}$ and $10^{-3}$ -- representing mm- and cm-sized particles interior of the water ice line -- concentrate themselves via the streaming instability at a solid abundance of a few percent. We thus revise a previously published critical solid abundance curve for the regime of $\tau_\mathrm{s} \ll 1$. The solid density in the concentrated regions reaches values higher than the Roche density, indicating that direct collapse of particles down to mm sizes into planetesimals is possible. Our results hence bridge the gap in particle size between direct dust growth limited by bouncing and the streaming instability.

Journal ArticleDOI
TL;DR: 1-D coupled climate-photochemical models are used to generate self-consistent atmospheres for several evolutionary scenarios, including high-O2, high-CO2, and more Earth-like atmospheres, with both oxic and anoxic compositions, and show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone.
Abstract: Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here we use 1D coupled climate-photochemical models to generate self-consistent atmospheres for evolutionary scenarios predicted in our companion paper (Barnes et al., 2016). These include high-O2, high-CO2, and more Earth-like atmospheres, with either oxidizing or reducing compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b, but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere, and JWST observations longward of 7 microns could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST, but may be within the reach of larger aperture telescopes. Direct imaging spectra may detect O4, which is diagnostic of massive water loss and O2 retention, rather than a photosynthesis. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 {\mu}m would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4 and either photosynthetically produced O2 or a hydrocarbon haze layer.

Journal ArticleDOI
TL;DR: In this article, the authors used a 3D Earth-analog GCM with self-consistent relationships between stellar metallicity, stellar effective temperature, and the planetary orbital/rotational period to estimate the inner edge of the habitable zone for synchronously rotating terrestrial planets around late-K and M-dwarf stars.
Abstract: Terrestrial planets at the inner edge of the habitable zone of late-K and M-dwarf stars are expected to be in synchronous rotation, as a consequence of strong tidal interactions with their host stars. Previous global climate model (GCM) studies have shown that, for slowly-rotating planets, strong convection at the substellar point can create optically thick water clouds, increasing the planetary albedo, and thus stabilizing the climate against a thermal runaway. However these studies did not use self-consistent orbital/rotational periods for synchronously rotating planets placed at different distances from the host star. Here we provide new estimates of the inner edge of the habitable zone for synchronously rotating terrestrial planets around late-K and M-dwarf stars using a 3-D Earth-analog GCM with self-consistent relationships between stellar metallicity, stellar effective temperature, and the planetary orbital/rotational period. We find that both atmospheric dynamics and the efficacy of the substellar cloud deck are sensitive to the precise rotation rate of the planet. Around mid-to-late M-dwarf stars with low metallicity, planetary rotation rates at the inner edge of the HZ become faster, and the inner edge of the habitable zone is farther away from the host stars than in previous GCM studies. For an Earth-sized planet, the dynamical regime of the substellar clouds begins to transition as the rotation rate approaches ~10 days. These faster rotation rates produce stronger zonal winds that encircle the planet and smear the substellar clouds around it, lowering the planetary albedo, and causing the onset of the water-vapor greenhouse climatic instability to occur at up to ~25% lower incident stellar fluxes than found in previous GCM studies. For mid-to-late M-dwarf stars with high metallicity and for mid-K to early-M stars, we agree with previous studies.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface, and concluded that ProExima b is a viable candidate habitable planet.
Abstract: Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

Journal ArticleDOI
TL;DR: In this article, the authors present the first global simulations of a weakly ionized disk that exhibits large-scale magnetized winds and observe self-organization, resulting in axisymmetric rings of density and associated pressurebumps.
Abstract: The planet-forming region of protoplanetary disks is cold, dense, and therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence is thought to be mostly absent, and another mechanism has to be found to explain gas accretion. It has been proposed that magnetized winds, launched from the ionized disk surface, could drive accretion in the presence of a large-scale magnetic field. The efficiency and the impact of these surface winds on the disk structure is still highly uncertain. We present the first global simulations of a weakly ionized disk that exhibits large-scale magnetized winds. We also study the impact of self-organization, which was previously demonstrated only in non-stratified models. We perform numerical simulations of stratified disks with the PLUTO code. We compute the ionization fraction dynamically, and account for all three non-ideal MHD effects: ohmic and ambipolar diffusions, and the Hall drift. Simplified heating and cooling due to non-thermal radiation is also taken into account in the disk atmosphere. We find that disks can be accreting or not, depending on the configuration of the large-scale magnetic field. Magnetothermal winds, driven both by magnetic acceleration and heating of the atmosphere, are obtained in the accreting case. In some cases, these winds are asymmetric, ejecting predominantly on one side of the disk. The wind mass loss rate depends primarily on the average ratio of magnetic to thermal pressure in the disk midplane. The non-accreting case is characterized by a meridional circulation, with accretion layers at the disk surface and decretion in the midplane. Finally, we observe self-organization, resulting in axisymmetric rings of density and associated pressure "bumps". The underlying mechanism and its impact on observable structures are discussed.

Journal ArticleDOI
TL;DR: In this paper, the authors proposed a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds, and showed that due to the slowness of the CO/CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary.
Abstract: The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BD) since their first detections twenty years ago, has always been the key role played by micron-size condensates, called "dust" or "clouds", in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transition, and a physical first-principle understanding of this transition is lacking. In this paper, we propose a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds. We show that, due to the slowness of the CO/CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary. The induced small-scale turbulent energy transport reduces the temperature gradient in the atmosphere, explaining the observed increase in near infrared J-H and J-K colors of L dwarfs and hot EGPs, while a warming up of the deep atmosphere along the L to T transition, as the CO/CH4 instability vanishes, naturally solves the two aforementioned puzzles, and provides a physical explanation of the L to T transition. This new picture leads to a drastic revision of our understanding of BD and EGP atmospheres and their evolution.

Journal ArticleDOI
TL;DR: In this paper, the authors uncover a dramatic distinction in the prevalence of sub-Jovian companions, between systems that contain hot Jupiters (periods inward of 10 days) and those that host warm Jupsiters(periods between 10 and 200 days).
Abstract: Exploiting the Kepler transit data, we uncover a dramatic distinction in the prevalence of sub-Jovian companions, between systems that contain hot Jupiters (periods inward of 10 days) and those that host warm Jupiters (periods between 10 and 200 days). Hot Jupiters, with the singular exception of WASP-47b, do not have any detectable inner or outer planetary companions (with periods inward of 50 days and sizes down to $2 R_{\rm Earth}$). Restricting ourselves to inner companions, our limits reach down to $1 R_{\rm Earth}$. In stark contrast, half of the warm Jupiters are closely flanked by small companions. Statistically, the companion fractions for hot and warm Jupiters are mutually exclusive, particularly in regard to inner companions. The high companion fraction of warm Jupiters also yields clues to their formation. The warm Jupiters that have close-by siblings should have low orbital eccentricities and low mutual inclinations. The orbital configurations of these systems are reminiscent of those of the low-mass, close-in planetary systems abundantly discovered by the Kepler mission. This, and other arguments, lead us to propose that these warm Jupiters are formed in-situ. There are indications that there may be a second population of warm Jupiters with different characteristics. In this picture, WASP-47b could be regarded as the extending tail of the in-situ warm Jupiters into the hot Jupiter region, and does not represent the generic formation route for hot Jupiters.

Journal ArticleDOI
TL;DR: In this article, the first year of the NASA K2 mission (Campaigns 0-4) was used to discover 197 candidates for Earth-like terrestrial worlds, with median values of R_P = 2.3 R_E, P=8.6 d, Tef = 5300 K, and Kp=12.7 mag.
Abstract: We present 197 planet candidates discovered using data from the first year of the NASA K2 mission (Campaigns 0-4), along with the results of an intensive program of photometric analyses, stellar spectroscopy, high-resolution imaging, and statistical validation. We distill these candidates into sets of 104 validated planets (57 in multi-planet systems), 30 false positives, and 63 remaining candidates. Our validated systems span a range of properties, with median values of R_P = 2.3 R_E, P=8.6 d, Tef = 5300 K, and Kp=12.7 mag. Stellar spectroscopy provides precise stellar and planetary parameters for most of these systems. We show that K2 has increased by 30% the number of small planets known to orbit moderately bright stars (1-4 R_E, Kp=9-13 mag). Of particular interest are 37 planets smaller than 2 R_E, 15 orbiting stars brighter than Kp=11.5, five receiving Earth-like irradiation levels, and several multi-planet systems -- including four planets orbiting the M dwarf K2-72 near mean-motion resonances. By quantifying the likelihood that each candidate is a planet we demonstrate that our candidate sample has an overall false positive rate of 15-30%, with rates substantially lower for small candidates ( 8 R_E and/or with P < 3 d. Extrapolation of the current planetary yield suggests that K2 will discover between 500-1000 planets in its planned four-year mission -- assuming sufficient follow-up resources are available. Efficient observing and analysis, together with an organized and coherent follow-up strategy, is essential to maximize the efficacy of planet-validation efforts for K2, TESS, and future large-scale surveys.

Journal ArticleDOI
TL;DR: In this paper, the authors show that the non-resonant-to-to--resonance ratio obtained with the grainy migration is higher, up to ~10 times higher for the range of parameters investigated here, than in a model with smooth migration, and they find that the probability for an outer disk object to end up on a stable orbit in the Kuiper belt is ~10^{-3}.
Abstract: The Kuiper belt is a population of icy bodies beyond the orbit of Neptune. The complex orbital structure of the Kuiper belt, including several categories of objects inside and outside of resonances with Neptune, emerged as a result of Neptune's migration into an outer planetesimal disk. An outstanding problem with the existing migration models is that they invariably predict excessively large resonant populations, while observations show that the non-resonant orbits are in fact common (e.g., the main belt population is 2-4 times larger than Plutinos in the 3:2 resonance). Here we show that this problem can resolved if it is assumed that Neptune's migration was grainy, as expected from scattering encounters of Neptune with massive planetesimals. The grainy migration acts to destabilize resonant bodies with large libration amplitudes, a fraction of which ends up on stable non-resonant orbits. Thus, the non-resonant--to--resonant ratio obtained with the grainy migration is higher, up to ~10 times higher for the range of parameters investigated here, than in a model with smooth migration. The grainy migration leads to a narrower distribution of the libration amplitudes in the 3:2 resonance. The best fit to observations is obtained when it is assumed that the outer planetesimal disk below 30 AU contained 1000-4000 Plutos, or ~1000 bodies twice as massive as Pluto. We find that the probability for an outer disk object to end up on a stable orbit in the Kuiper belt is ~10^{-3}. Together, these results are consistent with having only two (known) Pluto-mass bodies in the Kuiper belt today (Pluto and Eris). The combined mass of Pluto-class objects in the original disk was ~2-8 Earth masses (M_Earth), which represents 10-40% of the estimated disk mass (M_disk ~ 20 M_Earth).

Journal ArticleDOI
TL;DR: In this article, the decorrelated light curves are searched for transits with an adapted version of the CoRoT alarm pipeline, which is optimized for bright stars for which spectroscopic follow-up is possible.
Abstract: With only two functional reaction wheels, Kepler cannot maintain stable pointing at its original target field and entered a new mode of observation called K2. Our method is based on many years of experience in planet hunting for the CoRoT mission. Due to the unstable pointing, K2 light curves present systematics that are correlated with the target position in the CCD. Therefore, our pipeline also includes a decorrelation of this systematic noise. Our pipeline is optimised for bright stars for which spectroscopic follow-up is possible. We achieve a maximum precision on 6 hours of 6 ppm. The decorrelated light curves are searched for transits with an adapted version of the CoRoT alarm pipeline. We present 172 planetary candidates and 327 eclipsing binary candidates from campaigns 1, 2, 3, 4, 5 and 6 of K2. Both the planetary candidates and eclipsing binary candidates lists are made public to promote follow-up studies. The light curves will also be available to the community.

Journal ArticleDOI
TL;DR: In this paper, the authors review the theory of atmospheric chemistry in all classes of exoplanets discovered to date, from highly irradiated gas giants, ice giants, and super-Earths, to directly imaged giant planets at large orbital separations.
Abstract: Characterizing the atmospheres of extrasolar planets is the new frontier in exoplanetary science. The last two decades of exoplanet discoveries have revealed that exoplanets are very common and extremely diverse in their orbital and bulk properties. We now enter a new era as we begin to investigate the chemical diversity of exoplanets, their atmospheric and interior processes, and their formation conditions. Recent developments in the field have led to unprecedented advancements in our understanding of atmospheric chemistry of exoplanets and the implications for their formation conditions. We review these developments in the present work. We review in detail the theory of atmospheric chemistry in all classes of exoplanets discovered to date, from highly irradiated gas giants, ice giants, and super-Earths, to directly imaged giant planets at large orbital separations. We then review the observational detections of chemical species in exoplanetary atmospheres of these various types using different methods, including transit spectroscopy, doppler spectroscopy, and direct imaging. In addition to chemical detections, we discuss the advances in determining chemical abundances in these atmospheres and how such abundances are being used to constrain exoplanetary formation conditions and migration mechanisms. Finally, we review recent theoretical work on the atmospheres of habitable exoplanets, followed by a discussion of future outlook of the field.

Journal ArticleDOI
TL;DR: In this article, the authors presented results from the analysis of 82 new radial velocity observations of this system obtained with HARPS-N, together with the existing 14 HIRES data points, and concluded that Kepler-21b belongs to the population of terrestrial planets with iron, magnesium silicate interiors, which have lost the majority of their envelope volatiles via stellar winds or gravitational escape.
Abstract: HD 179070, aka Kepler-21, is a V = 8.25 F6IV star and the brightest exoplanet host discovered by Kepler. An early detailed analysis by Howell et al. (2012) of the first thirteen months (Q0 - Q5) of Kepler light curves revealed transits of a planetary companion, Kepler-21b, with a radius of about 1.60 +/- 0.04 R_earth and an orbital period of about 2.7857 days. However, they could not determine the mass of the planet from the initial radial velocity observations with Keck-HIRES, and were only able to impose a 2-sigma upper limit of 10 M_earth. Here we present results from the analysis of 82 new radial velocity observations of this system obtained with HARPS-N, together with the existing 14 HIRES data points. We detect the Doppler signal of Kepler-21b with a radial velocity semi-amplitude K = 2.00 +/- 0.65 m/s, which corresponds to a planetary mass of 5.1 +/- 1.7 M_earth. We also measure an improved radius for the planet of 1.639 (+0.019, -0.015) R_earth, in agreement with the radius reported by Howell et al. (2012). We conclude that Kepler-21b, with a density of 6.4 +/- 2.1 g/cm^3, belongs to the population of terrestrial planets with iron, magnesium silicate interiors, which have lost the majority of their envelope volatiles via stellar winds or gravitational escape. The radial velocity analysis presented in this paper serves as example of the type of analysis that will be necessary to confirm the masses of TESS small planet candidates.