The Effects on a Protoplanetary Disk by a Forming Giant Planet Induced Gap

Formation of a disc gap induced by a planet: effect of the deviation from Keplerian disc rotation

Authors:

Kanangawa et al

Abstract:

The gap formation induced by a giant planet is important in the evolution of the planet and the protoplanetary disc. We examine the gap formation by a planet with a new formulation of one-dimensional viscous discs which takes into account the deviation from Keplerian disc rotation due to the steep gradient of the surface density. This formulation enables us to naturally include the Rayleigh stable condition for the disc rotation. It is found that the derivation from Keplerian disc rotation promotes the radial angular momentum transfer and makes the gap shallower than in the Keplerian case. For deep gaps, this shallowing effect becomes significant due to the Rayleigh condition. In our model, we also take into account the propagation of the density waves excited by the planet, which widens the range of the angular momentum deposition to the disc. The effect of the wave propagation makes the gap wider and shallower than the case with instantaneous wave damping. With these shallowing effects, our one-dimensional gap model is consistent with the recent hydrodynamic simulations.

The Gaia Mission, Binary Stars and Exoplanets

The Gaia Mission, Binary Stars and Exoplanets

Authors:

Eyer et al

Abstract:

On the 19th of December 2013, the Gaia spacecraft was successfully launched by a Soyuz rocket from French Guiana and started its amazing journey to map and characterise one billion celestial objects with its one billion pixel camera. In this presentation, we briefly review the general aims of the mission and describe what has happened since launch, including the Ecliptic Pole scanning mode. We also focus especially on binary stars, starting with some basic observational aspects, and then turning to the remarkable harvest that Gaia is expected to yield for these objects.

Resonances in Retrograde Circumbinary Disks

Resonances in retrograde circumbinary discs

Authors:

Nixon et al

Abstract:

We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gravitational potential of the binary which rotate in a retrograde sense to the binary orbit and so rotate progradely with respect to this disc, allowing a resonant interaction to occur between the binary and the disc. The resulting resonant torques distinctly alter the disc response from the circular binary case. We describe results of three-dimensional hydrodynamic simulations to explore this effect and categorise the response of the disc in terms of modes whose strengths vary as a function of binary mass ratio and eccentricity. These mode strengths are weak compared to the largest mode strengths expected in the prograde case where the binary and disc rotate in the same sense. However, for sufficiently high binary eccentricity, resonant torques open a gap in a retrograde circumbinary disc, while permitting gas inflow on to the binary via gas streams. The inflow results in a time varying accretion rate on to the binary that is modulated over the binary orbital period, as was previously found to occur in the prograde case.

A High Deuterium Ratio is From Continued D Fractionation in Protoplanetary Disk

Chemistry in disks

Authors:

Teague et al

Abstract:

Aims.

We study the deuteration and ionization structure of the DM Tau disk via interferometric observations and modelling of the key molecular ions, HCO+ and DCO+.

Methods.

The Plateau de Bure Array is used to observe DM Tau in lines of HCO+ (1−0), (3−2) and DCO+ (3−2) with a ~ 1.5′′ angular and ~0.2 km s-1 spectral resolution. Using a power-law fitting approach the observed column densities profiles are derived and thus the isotopic ratio RD = DCO+/HCO+. Chemical modelling allowed an exploration of the sensitivity of HCO+ and DCO+ abundances to physical parameters out with temperature. A steady state approximation was employed to observationally constrain the ionization fraction x(e−).

Results.

Fitting of radiative transfer models suggests that there is a chemical hole in HCO+ and DCO+, extending up to 50 AU from the star. More work is required to discern the cause of this. The observed column densities of HCO+ and DCO+ at 100 AU were (9.8+0.3-0.7) × 1012 and (1.2 ± 0.7) × 1012 cm-2 respectively. Where both HCO+ and DCO+ were present, RD was found to increase radially from 0.1 at 50 AU to 0.2 at 450 AU. This behaviour was well reproduced by the chemical model. The X-ray luminosity of the central star, the interstellar UV and CO depletion were found to be the most important physical parameters controlling the abundances of HCO+ and DCO+. Differences in the vertical extent of HCO+ and DCO+ molecular layers resulted in different responses to changing physical parameters, manifesting as radial gradients in RD. The ionization fraction was found to be x(e−) ~ 10-7 in the molecular layer, comparable to the disk averaged value. Modelling shows that while HCO+ is the most dominant charged molecular ion in our disk model, atomic ions, such as C+, S+, H+, Na+ and Mg+, dominate the charge in both the molecular layer and disk atmosphere.

Conclusions.

A high value of RD is indicative of continued deuterium fractionation in a protoplanetary disk after pre/protostellar phases. Radial properties of RD can be employed to discern the importance of ionization from X-rays and UV, thus necessitating the need for more, high resolution observations of DCO+ and other deuterated species in disks. A steady-state approach commonly adopted for constraining ionization degree in prestellar cores is not applicable for disks where accurate determination of the ionization fraction in the molecular layer requires knowledge of the atomic ions present as molecular ions are relatively sparse.

Kepler-447b: A VERY Weird Hot Jupiter

Kepler-447b: a hot-Jupiter with an extremely grazing transit

Authors:

Lilo-Box et al

Abstract:

We present the radial velocity confirmation of the extrasolar planet Kepler-447b, initially detected as a candidate by the Kepler mission. In this work, we analyze its transit signal and the radial velocity data obtained with the Calar Alto Fiber-fed Echelle spectrograph (CAFE). By simultaneously modeling both datasets, we obtain the orbital and physical properties of the system. According to our results, Kepler-447b is a Jupiter-mass planet (Mp=1.37+0.48−0.46MJup), with an estimated radius of Rp=1.65+0.59−0.56RJup (uncertainties provided in this work are 3σ unless specified). This translates into a sub-Jupiter density. The planet revolves every ∼7.8 days around a G8V star with detected activity in the Kepler light curve. Kepler-447b transits its host with a large impact parameter (b=1.076+0.112−0.086), being one of the few planetary grazing transits confirmed so far and the first in the Kepler large crop of exoplanets. We estimate that only around 20% of the projected planet disk occults the stellar disk. The relatively large uncertainties in the planet radius are due to the large impact parameter and short duration of the transit. Within the transit time interval, we find the presence of large (somehow modulated) outliers during the transit. We propose and analyze different scenarios that could explain these brighter data points, including instrumental effects, additional perturbing bodies, stellar pulsations, rotation of a non-spherical planet, and spot-crossing events. However, short-cadence photometric data (at the 1 minute level) is still needed to unveil the nature of this observational effect.

The HD 130322 ExoPlanetary System Revisited

Refined Properties of the HD 130322 Planetary System

Authors:

Hinkel et al

Abstract:

Exoplanetary systems closest to the Sun, with the brightest host stars, provide the most favorable opportunities for characterization studies of the host star and their planet(s). The Transit Ephemeris Refinement and Monitoring Survey uses both new radial velocity measurements and photometry in order to greatly improve planetary orbit uncertainties and the fundamental properties of the star, in this case HD 130322. The only companion, HD 130322b, orbits in a relatively circular orbit, e = 0.029 every ~10.7 days. Radial velocity measurements from multiple sources, including 12 unpublished from the Keck I telescope, over the course of ~14 years have reduced the uncertainty in the transit midpoint to ~2 hours. The transit probability for the b-companion is 4.7%, where M_p sin i = 1.15 M_J and a = 0.0925 AU. In this paper, we compile photometric data from the T11 0.8m Automated Photoelectric Telescope at Fairborn Observatory taken over ~14 years, including the constrained transit window, which results in a dispositive null result for both full transit exclusion of HD 130322b to a depth of 0.017 mag and grazing transit exclusion to a depth of ~0.001 mag. Our analysis of the starspot activity via the photometric data reveals a highly accurate stellar rotation period: 26.53 +/-0.70 days. In addition, the brightness of the host with respect to the comparison stars is anti-correlated with the Ca II H and K indices, typical for a young solar-type star.

B Class Stars in the Quintuplet and Arches Galactic Center Clusters

Circumstellar discs in Galactic centre clusters: Disc-bearing B-type stars in the Quintuplet and Arches clusters

Authors:

Stolte et al

Abstract:

We investigate the circumstellar disc fraction as determined from L-band excess observations of the young, massive Arches and Quintuplet clusters residing in the central molecular zone of the Milky Way. The Quintuplet cluster was searched for L-band excess sources for the first time. We find a total of 26 excess sources in the Quintuplet cluster and 21 in the Arches cluster, of which 13 are new detections. With the aid of proper motion membership samples, the disc fraction of the Quintuplet cluster was derived for the first time to be 4.0 +/- 0.7%. There is no evidence for a radially varying disc fraction in this cluster. In the case of the Arches cluster, a disc fraction of 9.2 +/- 1.2% approximately out to the cluster's predicted tidal radius, r less than 1.5 pc, is observed. This excess fraction is consistent with our previously found disc fraction in the cluster in the radial range 0.3 less than r less than 0.8 pc. In both clusters, the host star mass range covers late A- to early B-type stars, 2 less than M less than 15 Msun, as derived from J-band photospheric magnitudes. We discuss the unexpected finding of dusty circumstellar discs in these UV intense environments in the context of primordial disc survival and formation scenarios of secondary discs. We consider the possibility that the L-band excess sources in the Arches and Quintuplet clusters could be the high-mass counterparts to T Tauri pre-transitional discs. As such a scenario requires a long pre-transitional disc lifetime in a UV intense environment, we suggest that mass transfer discs in binary systems are a likely formation mechanism for the B-star discs observed in these starburst clusters.

Direct Detection of Exoplanet KIC 12557548b's Dust Cloud

Direct evidence for an evolving dust cloud from the exoplanet KIC 12557548 b

Authors:

Bochinski et al

Abstract:

We present simultaneous multi-color optical photometry using ULTRACAM of the transiting exoplanet KIC 12557548 b (also known as KIC 1255 b). This reveals, for the first time, the color dependence of the transit depth. Our g and z transits are similar in shape to the average Kepler short-cadence profile, and constitute the highest-quality extant coverage of individual transits. Our Night 1 transit depths are 0.85 +/- 0.04% in z; 1.00 +/- 0.03% in g; and 1.1 +/- 0.3% in u. We employ a residual-permutation method to assess the impact of correlated noise on the depth difference between the z and g bands and calculate the significance of the color dependence at 3.2{\sigma}. The Night 1 depths are consistent with dust extinction as observed in the ISM, but require grain sizes comparable to the largest found in the ISM: 0.25-1{\mu}m. This provides direct evidence in favor of this object being a disrupting low-mass rocky planet, feeding a transiting dust cloud. On the remaining four nights of observations the object was in a rare shallow-transit phase. If the grain size in the transiting dust cloud changes as the transit depth changes, the extinction efficiency is expected to change in a wavelength- and composition-dependent way. Observing a change in the wavelength-dependent transit depth would offer an unprecedented opportunity to determine the composition of the disintegrating rocky body KIC 12557548 b. We detected four out-of-transit u band events consistent with stellar flares.

Convection and Mixing in Giant Planet Evolution

Convection and Mixing in Giant Planet Evolution

Authors:

Vazan et al

Abstract:

The primordial internal structures of gas giant planets are unknown. Often giant planets are modeled under the assumption that they are adiabatic, convective, and homogeneously mixed, but this is not necessarily correct. In this work, we present the first self-consistent calculation of convective transport of both heat and material as the planets evolve. We examine how planetary evolution depends on the initial composition and its distribution, whether the internal structure changes with time, and if so, how it affects the evolution. We consider various primordial distributions, different compositions, and different mixing efficiencies and follow the distribution of heavy elements in a Jupiter-mass planet as it evolves. We show that a heavy-element core cannot be eroded by convection if there is a sharp compositional change at the core-envelope boundary. If the heavy elements are initially distributed within the planet according to some compositional gradient, mixing occurs in the outer regions resulting in a compositionally homogeneous outer envelope. Mixing of heavy materials that are injected in a convective gaseous envelope are found to mix efficiently. Our work demonstrates that the primordial internal structure of a giant planet plays a substantial role in determining its long-term evolution and that giant planets can have non-adiabatic interiors. These results emphasize the importance of coupling formation, evolution, and internal structure models of giant planets self-consistently.

HD 131835 has an Unusually Warm Debris Disk

Discovery of Resolved Debris Disk Around HD 131835

Authors:

Hung et al

Abstract:

We report the discovery of the resolved disk around HD 131835 and present the analysis and modeling of its thermal emission. HD 131835 is a ~15 Myr A2 star in the Scorpius-Centaurus OB association at a distance of 122.7 +16.2 -12.8 parsec. The extended disk has been detected to ~1.5" (200 AU) at 11.7 {\mu}m and 18.3 {\mu}m with T-ReCS on Gemini South. The disk is inclined at an angle of ~75{\deg} with the position angle of ~61{\deg}. The flux of HD 131835 system is 49.3+-7.6 mJy and 84+-45 mJy at 11.7 {\mu}m and 18.3 {\mu}m respectively. A model with three grain populations gives a satisfactory fit to both the spectral energy distribution and the images simultaneously. This best-fit model is composed of a hot continuous power-law disk and two rings. We characterized the grain temperature profile and found that the grains in all three populations are emitting at temperatures higher than blackbodies. In particular, the grains in the continuous disk are unusually warm; even when considering small graphite particles as the composition.

A Pulsating Hot Subdwarf B Star With a Brown Dwarf Companion

An eclipsing post common-envelope system consisting of a pulsating hot subdwarf B star and a brown dwarf companion

Authors:

Schaffenroth et al

Abstract:

Hot subdwarf B stars (sdBs) are evolved, core helium-burning objects located on the extreme horizontal branch. Their formation history is still puzzling as the sdB progenitors must lose nearly all of their hydrogen envelope during the red-giant phase. About half of the known sdBs are in close binaries with periods from 1.2 h to a few days, a fact that implies they experienced a common-envelope phase. Eclipsing hot subdwarf binaries (also called HW Virginis systems) are rare but important objects for determining fundamental stellar parameters. Even more significant and uncommon are those binaries containing a pulsating sdB, as the mass can be determined independently by asteroseismology.

Here we present a first analysis of the eclipsing hot subdwarf binary V2008-1753. The light curve shows a total eclipse, a prominent reflection effect, and low--amplitude pulsations with periods from 150 to 180 s. An analysis of the light-- and radial velocity (RV) curves indicates a mass ratio close to q=0.146, an RV semi-amplitude of K=54.6kms−1, and an inclination of i=86.8∘. Combining these results with our spectroscopic determination of the surface gravity, logg=5.83, the best--fitting model yields an sdB mass of 0.47M⊙ and a companion mass of 69 MJup. As the latter mass is below the hydrogen-burning limit, V2008-1753 represents the first HW Vir system known consisting of a pulsating sdB and a brown dwarf companion. Consequently, it holds great potential for better constraining models of sdB binary evolution and asteroseismology.

Eclipsing Subdwarf B Binary HS0705+6700 Appears to Have Exoplanets

The eclipsing binary HS0705+6700 and the search for circumbinary objects

Authors:

Pulley et al

Abstract:

HS0705+6700 (also identified as V470 Cam) is a short period (2.3 h) post common envelope detached eclipsing sdB binary system which exhibits transit time variations (TTVs) of a cyclical nature. We report a further 25 timings of light minima and show that our new TTVs support and extend this cyclical pattern to 1.6 periods. We examine possible causes of the observed TTVs and confirm that the presence of a third, and possibly a fourth, body could provide an elegant explanation of these cyclical variations. However other non-circumbinary mechanisms, e.g. Applegate magnetic dynamo effects, will remain possible contenders until sufficient data has been accumulated to demonstrate that the periodicity of the TTVs is time independent.

Detecting the PaleoArchean Geodynamo & What it Means for Habitability

Detecting the oldest geodynamo and attendant shielding from the solar wind: Implications for habitability

Authors:

Taduno et al

Abstract:

The onset and nature of the earliest geomagnetic field is important for understanding the evolution of the core, atmosphere and life on Earth. A record of the early geodynamo is preserved in ancient silicate crystals containing minute magnetic inclusions. These data indicate the presence of a geodynamo during the Paleoarchean, between 3.4 and 3.45 billion years ago. While the magnetic field sheltered Earth's atmosphere from erosion at this time, standoff of the solar wind was greatly reduced, and similar to that during modern extreme solar storms. These conditions suggest that intense radiation from the young Sun may have modified the atmosphere of the young Earth by promoting loss of volatiles, including water. Such effects would have been more pronounced if the field were absent or very weak prior to 3.45 billion years ago, as suggested by some models of lower mantle evolution. The frontier is thus trying to obtain geomagnetic field records that are much greater than 3.45 billion-years-old, as well as constraining solar wind pressure for these times. In this review we suggest pathways for constraining these parameters and the attendant history of Earth's deep interior, hydrosphere and atmosphere. In particular, we discuss new estimates for solar wind pressure for the first 700 million years of Earth history, the competing effects of magnetic shielding versus solar ion collection, and bounds on the detection level of a geodynamo imposed by the presence of external fields. We also discuss the prospects for constraining Hadean-Paleoarchean magnetic field strength using paleointensity analyses of zircons.

A new way to Study Carbon-rich Exoplanet Atmospheres

A new chemical scheme to study carbon-rich exoplanet atmospheres

Authors:

Venot et al

Abstract:

Atmospheres with a high C/O ratio are expected to contain an important quantity of hydrocarbons, including heavy molecules (with more than 2 carbon atoms). To study correctly these C-rich atmospheres, a chemical scheme adapted to this composition is necessary. We have implemented a chemical scheme that can describe the kinetics of species with up to 6 carbon atoms. This chemical scheme has been developed with specialists of combustion and validated through experiments on a wide range of T and P. This chemical network is available on the online database KIDA. We have created a grid of 12 models to explore different thermal profiles and C/O ratios. For each of them, we have compared the chemical composition determined with a C0-C2 chemical scheme (species with up to 2 carbon atoms) and with the C0-C6 scheme. We found no difference in the results obtained with the two schemes when photolyses are not included in the model, whatever the temperature of the atmosphere. In contrast, when there is photochemistry, differences can appear in the upper atmosphere. These differences are found for all the tested PT profiles in the case that the C/O ratio is above 1. When the C/O ratio of the atmosphere is solar, differences are only found at temperatures lower than 1000K. The differences linked to the use of different chemical schemes do not have important influence on the synthetic spectra. However, we have confirmed that C2H2 and HCN as possible tracers of warm C-rich atmospheres. The use of this new chemical scheme is mandatory to model atmospheres with a high C/O ratio and, in particular, if one is interested in studying in details the photochemistry. If one is just interested in the synthetic spectra, the use of a smaller scheme may be sufficient.

Terrestrial Planet Formation

Earth and Terrestrial Planet Formation

Authors:

Jacobson et al

Abstract:

The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with the nebula gas lead to remarkable new accretion mechanisms for the formation of planetesimals and the growth of planetary embryos.

Many of the popular models for the later stages of planet formation are presented. The classical models with the giant planets on fixed orbits are not consistent with the known history of the Solar System, fail to create a high Earth/Mars mass ratio, and, in many cases, are also internally inconsistent. The successful Grand Tack model creates a small Mars, a wet Earth, a realistic asteroid belt and the mass-orbit structure of the terrestrial planets.

In the Grand Tack scenario, growth curves for Earth most closely match a Weibull model. The feeding zones, which determine the compositions of Earth and Venus follow a particular pattern determined by Jupiter, while the feeding zones of Mars and Theia, the last giant impactor on Earth, appear to randomly sample the terrestrial disk. The late accreted mass samples the disk nearly evenly.

Will Structure of HR 8799's Inner Debris Disk Reveal Water Delivery Mechanism for Potential Terrestrial Exoplanets?

The dynamical structure of HR 8799's inner debris disk

Authors:

Contro et al

Abstract:

The HR 8799 system, with its four giant planets and two debris belts, has an architecture closely mirroring that of our Solar system where the inner, warm asteroid belt and outer, cool Edgeworth-Kuiper belt bracket the giant planets. As such, it is a valuable laboratory for examining exoplanetary dynamics and debris disk-exoplanet interactions. Whilst the outer debris belt of HR 8799 has been well resolved by previous observations, the spatial extent of the inner disk remains unknown. This leaves a significant question mark over both the location of the planetesimals responsible for producing the belt's visible dust and the physical properties of those grains. We have performed the most extensive simulations to date of the inner, unresolved debris belt around HR 8799, using UNSW Australia's Katana supercomputing facility to follow the dynamical evolution of a model inner disk comprising 300,298 particles for a period of 60 million years. These simulations have enabled the characterisation of the extent and structure of the inner disk in detail, and will in future allow us to provide a first estimate of the small-body impact rate and water delivery prospects for possible (as-yet undetected) terrestrial planet(s) in the inner system.

The Implications of Colliding Stellar Winds in Low Mass Binary Star Systems for Habitable Worlds

Colliding Winds in Low-Mass Binary Star Systems: wind interactions and implications for habitable planets

Authors:

Johnstone et al

Abstract:

Context.

In binary star systems, the winds from the two components impact each other, leading to strong shocks and regions of enhanced density and temperature. Potentially habitable circumbinary planets must continually be exposed to these interactions regions.

Aims.

We study, for the first time, the interactions between winds from low-mass stars in a binary system, to show the wind conditions seen by potentially habitable circumbinary planets.

Methods.

We use the advanced 3D numerical hydrodynamic code Nurgush to model the wind interactions of two identical winds from two solar mass stars with circular orbits and a binary separation of 0.5 AU. As input into this model, we use a 1D hydrodynamic simulation of the solar wind, run using the Versatile Advection Code. We derive the locations of stable and habitable orbits in this system to explore what wind conditions potentially habitable planets will be exposed to during their orbits.

Results.

Our wind interaction simulations result in the formation of two strong shock waves separated by a region of enhanced density and temperature. The wind-wind interaction region has a spiral shape due to Coriolis forces generated by the orbital motions of the two stars. The stable and habitable zone in this system extends from approximately 1.4 AU to 2.4 AU. Habitable planets have to pass through strong shock waves several times per orbit and spend a significant amount of time embedded in the higher density matter between the shocks. The enhanced density in the wind-wind interaction region is likely to lead to a 20% decrease in the size of a planet’s magnetosphere.

Conclusions.

Our results indicate that wind-wind interactions are likely to influence the magnetospheres and upper atmospheres of circumbinary planets and could have moderate implications for the development of habitable planetary environments.

How Habitable Planets get Their Water

Planetary Systems and the Formation of Habitable Planets

Authors:

Dvorak et al

Abstract:

As part of a national scientific network 'Pathways to Habitability' the formation of planets and the delivery of water onto these planets is a key question as water is essential for the development of life. In the first part of the paper we summarize the state of the art of planet formation - which is still under debate in the astronomical community - before we show our results on this topic. The outcome of our numerical simulations depends a lot on the choice of the initial distribution of planetesimals and planetary embryos after gas disappeared in the protoplanetary disk. We also take into account that some of these planetesimals of sizes in the order of the mass of the Moon already contained water; the quantity depends on the distance from the Sun - close-by bodies are dry, but starting from a distance of about 2 AU they can contain substantial amounts of water. We assume that the gas giants and terrestrial planets are already formed when we check the collisions of the small bodies containing water (in the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the crust and on the surface of Earth. In the second part we discuss in more detail how the formation of larger bodies after a collision may happen as the outcome depends on parameters like collision velocity, impact angle, and the materials involved. We present results obtained by SPH (Smooth Particle Hydrodynamics) simulations. We briefly describe this method and show different scenarios with respect to the formed bodies, possible fragmentation and the water content before and after the collision. In an appendix we discuss detection methods for extrasolar planets (close to 2000 such objects have been discovered so far).

X Ray and Cosmic Ray Ionization of TW Hyadrae's Protoplanetary Disk

CONSTRAINING THE X-RAY AND COSMIC-RAY IONIZATION CHEMISTRY OF THE TW Hya PROTOPLANETARY DISK: EVIDENCE FOR A SUB-INTERSTELLAR COSMIC-RAY RATE

Authors:

Cleeves et al

Abstract:

We present an observational and theoretical study of the primary ionizing agents (cosmic rays (CRs) and X-rays) in the TW Hya protoplanetary disk. We use a set of resolved and unresolved observations of molecular ions and other molecular species, encompassing 11 lines total, in concert with a grid of disk chemistry models. The molecular ion constraints comprise new data from the Submillimeter Array on HCO+, acquired at unprecedented spatial resolution, and data from the literature, including ALMA observations of N2H+. We vary the model incident CR flux and stellar X-ray spectra and find that TW Hya's HCO+ and N2H+ emission are best-fit by a moderately hard X-ray spectra, as would be expected during the "flaring" state of the star, and a low CR ionization rate, ζCR lsim 10–19 s–1. This low CR rate is the first indication of the presence of CR exclusion by winds and/or magnetic fields in an actively accreting T Tauri disk system. With this new constraint, our best-fit ionization structure predicts a low turbulence "dead-zone" extending from the inner edge of the disk out to 50-65 AU. This region coincides with an observed concentration of millimeter grains, and we propose that the inner region of TW Hya is a dust (and possibly planet) growth factory as predicted by previous theoretical work.

Is There a Link Between Exoplanets & Circumstellar Debris Disks?

PURSUING THE PLANET–DEBRIS DISK CONNECTION: ANALYSIS OF UPPER LIMITS FROM THE ANGLO-AUSTRALIAN PLANET SEARCH

Authors:

Wittenmyer et al

Abstract:

Solid material in protoplanetary disks will suffer one of two fates after the epoch of planet formation; either being bound up into planetary bodies, or remaining in smaller planetesimals to be ground into dust. These end states are identified through detection of sub-stellar companions by periodic radial velocity (or transit) variations of the star, and excess emission at mid- and far-infrared wavelengths, respectively. Since the material that goes into producing the observable outcomes of planet formation is the same, we might expect these components to be related both to each other and their host star. Heretofore, our knowledge of planetary systems around other stars has been strongly limited by instrumental sensitivity. In this work, we combine observations at far-infrared wavelengths by IRAS, Spitzer, and Herschel with limits on planetary companions derived from non-detections in the 16 year Anglo-Australian Planet Search to clarify the architectures of these (potential) planetary systems and search for evidence of correlations between their constituent parts. We find no convincing evidence of such correlations, possibly owing to the dynamical history of the disk systems, or the greater distance of the planet-search targets. Our results place robust limits on the presence of Jupiter analogs which, in concert with the debris disk observations, provides insights on the small-body dynamics of these nearby systems.

Magnetic Turbulence and Thermodynamics in the Inner Region of Protoplanetary Disks

Magnetic Turbulence and Thermodynamics in the Inner Region of Protoplanetary Discs

Author:

Hirose

Abstract:

Using radiation magnetohydrodynamics simulations with realistic opacities and equation of state, and zero net magnetic flux, we have explored thermodynamics in the inner part of protoplanetary discs where magnetic turbulence is expected. The thermal equilibrium curve consists of the upper, lower, and middle branches. The upper (lower) branch corresponds to hot (cool) and optically very (moderately) thick discs, respectively, while the middle branch is characterized by convective energy transport near the midplane. Convection is also the major energy transport process near the low surface density end of the upper branch. There, convective motion is fast with Mach numbers reaching ≳0.01, and enhances both magnetic turbulence and cooling, raising the ratio of vertically-integrated shear stress to vertically-integrated pressure by a factor of several. This convectively enhanced ratio seems a robust feature in accretion discs having an ionization transition. We have also examined causes of the S-shaped thermal equilibrium curve, as well as the thermal stability of the equilibrium solutions. Finally, we compared our results with the disc instability models used to explain FU Ori outbursts. Although the thermal equilibrium curve in our results also exhibits bistability, the surface density contrast across the bistability is an order of magnitude smaller, and the stress-to-pressure ratios in both upper and lower branches are two orders of magnitude greater, than those favored in the disc instability models. It therefore appears likely that FU Ori outbursts are not due solely to a thermal-viscous limit cycle resulting from accretion driven by local magnetic turbulence.

Signatures of Magnetic Driven Turbulence in Protoplanetary Disks

Signatures of MRI-Driven Turbulence in Protoplanetary Disks: Predictions for ALMA Observations

Authors:

Simon et al

Abstract:

Spatially resolved observations of molecular line emission provide unique constraints on protoplanetary disk turbulence. Using local non-ideal MHD simulations and radiative transfer calculations, we assess the ability of ALMA observations to robustly detect and characterize disk turbulence. We specifically predict the outcome of the magnetorotational instability (MRI) in the disk around HD 163296, a promising observational target. We find that the MRI can support the observed level of accretion if the outer disk surface is ionized by far-UV photons and threaded by a weak net vertical magnetic field. We identify two classes of MRI solution - dynamo solutions in which the surface magnetic field reverses periodically, and non-dynamo solutions in which much of the Maxwell stress is steady and large scale. In both classes the small-scale turbulence increases in strength with height above the mid-plane, and can be represented as a microturbulent component. Using vertical profiles of the turbulent velocity from simulations at different radii, we use radiative transfer calculations to quantify the observational signatures. We show that the peak to line center flux ratio is a robust diagnostic of turbulence that is only mildly degenerate with uncertainties in disk temperature. For the CO(3-2) line variations in the predicted peak-to-trough ratio between our most and least turbulent models are ~15%. We develop predictions for other molecular lines and for channel maps whose morphology allows for independent constraints on turbulence.

Type I Planet Migration in a Magnetized Disk

Type I Planet Migration in a Magnetized Disk. II. Effect of Vertical Angular Momentum Transport

Authors:

Bans et al

Abstract:

We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a linear perturbation analysis in the ideal-MHD limit. We focus on wind-driving disks, in which a magnetic torque ∝B0z∂B0φ/∂z (where B0z and B0φ are the equilibrium vertical and azimuthal field components) induces vertical angular momentum transport. We derive the governing differential equation for the disk response and identify its resonances and turning points. For a disk containing a slightly subthermal, pure-B0z field, the total 3D torque is close to its value in the 2D limit but remains lower than the hydrodynamic torque. In contrast with the 2D pure-B0φ field model considered by Terquem (2003), inward migration is not reduced in this case when the field amplitude decreases with radius. The presence of a subdominant B0φ component whose amplitude increases from zero at z=0 has little effect on the torque when acting alone, but in conjunction with a B0z component it gives rise to a strong torque that speeds up the inward migration by a factor ≳200. This factor could, however, be reduced in a real disk by dissipation and magnetic diffusivity effects. Unlike all previously studied disk migration models, in the B0z+∂B0φ/∂z case the dominant contributions to the torque add with the same sign from the two sides of the planet. We attribute this behavior to a new mode of interaction wherein a planet moves inward by plugging into the disk's underlying angular momentum transport mechanism.

ExoPlanetary Migration in a Magnetized Disk Part I

Type I Planet Migration in a Magnetized Disk. I. Effect of Large-Scale Vertical and Azimuthal Field Components

Authors:

Uribe et al

Abstract:

We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a combination of numerical simulations in 2D and 3D and a linear perturbation analysis. Steady-state models of such disks require the inclusion of magnetic diffusivity. To make progress using ideal MHD, we focus on simplified field configurations, involving purely vertical (Bz) and azimuthal (Bφ) field components and a combination of the two. For each of the models we calculate the locations of the relevant resonances and of the turning points, which delineate the propagation regions of the MHD waves that transport angular momentum from the planet to the disk. We use both numerical and semianalytic methods to evaluate the cumulative back torque acting on the planet, and explore the effect of spatial gradients in the disk's physical variables on the results. We conclude that, under realistic (3D) circumstances, a large-scale magnetic field can slow down the inward migration that characterizes the underlying unmagnetized disk --- by up to a factor of ∼2 when the magnetic pressure approaches the thermal pressure --- but it cannot reverse it. A previous inference that a pure-Bϕ field whose amplitude decreases fast enough with radius leads to outward migration applies only in 2D. In fact, we find that, in 3D, a pure-Bϕ disk undergoes a rapid transition to turbulence on account of a magnetorotational instability that is triggered by the planet-induced appearance of a weak Bz component.

Best Picture Yet by Hubble Space Telescope of beta Pictoris Circumstellar Disk

link.

Royal Tyrrell Museum Lecture on Meteorite Clues to Planetary Formation

There is no Circumbinary Brown Dwarf Around V471 Tauri

The First Science Results from SPHERE: Disproving the Predicted Brown Dwarf around V471 Tau

Authors:

Hardy et al

Abstract:

Variations of eclipse arrival times have recently been detected in several post common envelope binaries consisting of a white dwarf and a main sequence companion star. The generally favoured explanation for these timing variations is the gravitational pull of one or more circumbinary substellar objects periodically moving the center of mass of the host binary. Using the new extreme-AO instrument SPHERE, we image the prototype eclipsing post-common envelope binary V471 Tau in search of the brown dwarf that is believed to be responsible for variations in its eclipse arrival times. We report that an unprecedented contrast of 12.1 magnitudes in the H band at a separation of 260 mas was achieved, but resulted in a non-detection. This implies that there is no brown dwarf present in the system unless it is three magnitudes fainter than predicted by evolutionary track models, and provides damaging evidence against the circumbinary interpretation of eclipse timing variations. In the case of V471 Tau, a more consistent explanation is offered with the Applegate mechanism, in which these variations are prescribed to changes in the quadrupole moment within the main-sequence star

28 T Dwarf Candidates in ρ Oph Dark Cloud L 1688

Searching for T dwarfs in the ρ Oph dark cloud L 1688

Authors:

Chiang et al

Abstract:

We present a list of T dwarf candidates in the dark cloud L 1688 in the ρ Oph star-forming region. These candidates are selected with infrared colours sensitive to T dwarf characteristics of methane absorptions and of cool atmospheres. The 1.6-μm methane feature is diagnosed by on–off imaging using an H-band and an intermediate-band methane filter, calibrated to a set of known brown dwarfs of M, L, and T types in the field. Another methane feature at 3.3 μm is traced with the Spitzer/Infrared Array Camera (IRAC) [3.6] − [4.5] colour. For cool atmospheres, the H − [4.5] and K − [4.5] colours are utilized. With an additional criterion of mid-infrared brightness to eliminate extragalactic interlopers, a total of 28 T dwarf candidates have been identified. A comprehensive assessment was conducted to estimate the level of contamination of our sample by young stellar variability, by extragalactic sources sharing the same colour behaviour, or by foreground T dwarfs. Though extragalactic sources may contribute up to about half of the false positives, our candidates show close spatial association with the dark cloud, rather than randomly distributed as a background population would have been. Furthermore, even though our candidates are not selected a priori by a colour–magnitude relation, they mostly follow the 1 Myr isochrones, ascertaining their youth. Our selection methodology provides guidance to search for T dwarfs in other star-forming regions. Our candidate list, when comparing with those in the literature, which often rely on a single criterion on cool temperature or methane, is more conservative but should be more secure for follow-up spectroscopic confirmation of a T dwarf sample at the early evolutionary stage.

Multi-epoch Monitoring of Extremely Cool Brown Dwarfs

The Brown-dwarf Atmosphere Monitoring (BAM) Project II: Multi-epoch monitoring of extremely cool brown dwarfs

Authors:

Rajan et al

Abstract:

With the discovery of Y dwarfs by the WISE mission, the population of field brown dwarfs now extends to objects with temperatures comparable to those of Solar System planets. To investigate the atmospheres of these newly identified brown dwarfs, we have conducted a pilot study monitoring an initial sample of three late T-dwarfs (T6.5, T8 and T8.5) and one Y-dwarf (Y0) for infrared photometric variability at multiple epochs. With J-band imaging, each target was observed for a period of 1.0h to 4.5h per epoch, which covers a significant fraction of the expected rotational period. These measurements represent the first photometric monitoring for these targets. For three of the four targets (2M1047, Ross 458C and WISE0458), multi-epoch monitoring was performed, with the time span between epochs ranging from a few hours to ~2 years. During the first epoch, the T8.5 target WISE0458 exhibited variations with a remarkable min-to-max amplitude of 13%, while the second epoch light curve taken ~2 years later did not note any variability to a 3% upper limit. With an effective temperature of ~600 K, WISE0458 is the coldest variable brown dwarf published to-date, and combined with its high and variable amplitude makes it a fascinating target for detailed follow-up. The three remaining targets showed no significant variations, with a photometric precision between 0.8% and 20.0%, depending on the target brightness. Combining the new results with previous multi-epoch observations of brown dwarfs with spectral types of T5 or later, the currently identified variables have locations on the colour-colour diagram better matched by theoretical models incorporating cloud opacities rather than cloud-free atmospheres. This preliminary result requires further study to determine if there is a definitive link between variability among late-T dwarfs and their location on the colour-colour diagram.

Assessing the Earth Analogues in the Later Stages of Their Habitable Lifetimes

In Search of Future Earths: Assessing the possibility of finding Earth analogues in the later stages of their habitable lifetimes

Authors:

O'Malley-James et al

Abstract:

Earth will become uninhabitable within 2-3 Gyr as a result of the moving boundaries of the habitable zone caused by the increasing luminosity of the Sun. Predictions about the future of habitable conditions on Earth include a decline in species diversity and habitat extent, ocean loss and changes in the magnitudes of geochemical cycles. However, testing these predictions on the present-day Earth is difficult. The discovery of a planet that is a near analogue to the far future Earth could provide a means to test these predictions. Such a planet would need to have an Earth-like biosphere history, requiring it to have been in its system's habitable zone (HZ) for Gyr-long periods during the system's past, and to be approaching the inner-edge of the HZ at present. Here we assess the possibility of finding this very specific type of exoplanet and discuss the benefits of analysing older Earths in terms of improving our understanding of long-term geological and bio-geological processes. As an illustrative example, G stars within 10 parsecs are assessed as potential old-Earth-analogue hosts. Surface temperature estimates for hypothetical inner-HZ Earth analogues are used to determine whether any such planets in these systems would be at the right stage in their late-habitable lifetimes to exhibit detectable biosignatures. Predictions from planet formation studies and biosphere evolution models suggest that only 0.36% of G stars in the solar neighbourhood could host an old-Earth-analogue. However, if the development of an Earth-like biosphere is assumed to be rare, requiring a sequence of low-probability events to occur, then such planets are unlikely to be found in the solar neighbourhood - although 1000s could be present in the galaxy as a whole.

Can we Figure out the Internal Structure of Rocky Exoplanets From Their Mass and Radius?

Can we constrain interior structure of rocky exoplanets from mass and radius measurements?

Authors:

Dorn et al

Abstract:

We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. We provide a general methodology of analyzing interior structures of exoplanets that may help to understand how interior models are distributed among star systems. The methodology we propose is sufficiently general to allow its future extension to more complex internal structure
including hydrogen- and water-rich exoplanets.

Delamnination of the Geosphere in SuperEarths

Delamination in super-Earths extrapolated from the Earth model

Authors:

Shoji et al

Abstract:

It is suggested that the delamination process, in which the mantle lithosphere is peeled into the asthenosphere, contributes to the topographies and magmatism of the Earth. We investigated the vigorousness of the delamination in super-Earths by applying the Earth model to planets of heavy mass. Delamination is induced in planets of mass 5M⊕5M⊕ by the negative buoyancy of the mantle lithosphere. However, assuming pressure dependent rheology, the thermal Rayleigh number decreases due to the high pressure in super-Earths and thus the magnitude of convection in the Moho decreases. Because reduced convection in the Moho weakens the peeling of the mantle lithosphere, the delaminated area is narrower. The magnitude of the heat flux caused by the delamination process is also reduced in planets large in size compared with Earth. Although further work is needed, our model indicates that delamination can transfer more heat than the conduction of the lithosphere if the planet's mass is less than 5M⊕.

V471 Tauri's Circumbinary Brown Dwarf is Missing

Some pairs of stars consist of two normal stars with slightly different masses. When the star of slightly higher mass ages and expands to become a red giant, material is transferred to other star and ends up surrounding both stars in a huge gaseous envelope. When this cloud disperses the two move closer together and form a very tight pair with one white dwarf , and one more normal star.

One such stellar pair is called V471 Tauri. It is a member of the Hyades star cluster in the constellation of Taurus and is estimated to be around 600 million years old and about 163 light-years from Earth. The two stars are very close and orbit each other every 12 hours. Twice per orbit one star passes in front of the other -- which leads to regular changes in the brightness of the pair observed from Earth as they eclipse each other.

A team of astronomers led by Adam Hardy (Universidad Valparaíso, Valparaíso, Chile) first used the ULTRACAM system on ESO's New Technology Telescope to measure these brightness changes very precisely. The times of the eclipses were measured with an accuracy of better than two seconds -- a big improvement on earlier measurements.

The eclipse timings were not regular, but could be explained well by assuming that there was a brown dwarf orbiting both stars whose gravitational pull was disturbing the orbits of the stars. They also found hints that there might be a second small companion object.

Up to now however, it has been impossible to actually image a faint brown dwarf so close to much brighter stars. But the power of the newly installed SPHERE instrument on ESO's Very Large Telescope allowed the team to look for the first time exactly where the brown dwarf companion was expected to be. But they saw nothing, even though the very high quality images from SPHERE should have easily revealed it.

"There are many papers suggesting the existence of such circumbinary objects, but the results here provide damaging evidence against this hypothesis," remarks Adam Hardy.

link.

HD 69830: A 7 Billion Year Old Star With Three Neptunes and an Asteroid Belt Closer Than the Snowline

STELLAR PARAMETERS FOR HD 69830, A NEARBY STAR WITH THREE NEPTUNE MASS PLANETS AND AN ASTEROID BELT

Authors:

Tanner et al

Abstract:

We used the CHARA Array to directly measure the angular diameter of HD 69830, home to three Neptune mass planets and an asteroid belt. Our measurement of 0.674 ± 0.014 mas for the limb-darkened angular diameter of this star leads to a physical radius of R * = 0.9058 ± 0.0190 R☉ and luminosity of L * = 0.622 ± 0.014 L☉ when combined with a fit to the spectral energy distribution of the star. Placing these observed values on an Hertzsprung-Russel diagram along with stellar evolution isochrones produces an age of 10.6 ± 4 Gyr and mass of 0.863 ± 0.043 M☉. We use archival optical echelle spectra of HD 69830 along with an iterative spectral fitting technique to measure the iron abundance ([Fe/H] = –0.04 ± 0.03), effective temperature (5385 ± 44 K), and surface gravity (log g = 4.49 ± 0.06). We use these new values for the temperature and luminosity to calculate a more precise age of 7.5 ± 3 Gyr. Applying the values of stellar luminosity and radius to recent models on the optimistic location of the habitable zone produces a range of 0.61-1.44 AU; partially outside the orbit of the furthest known planet (d) around HD 69830. Finally, we estimate the snow line at a distance of 1.95 ± 0.19 AU, which is outside the orbit of all three planets and its asteroid belt.

Orbital Instability of Close-in exomoons in Non-coplanar Systems

Orbital instability of close-in exomoons in non-coplanar systems

Authors:

Hong et al

Abstract:

This work shows the dynamical instability that can happen to close-in satellites when planet oblateness is not accounted for in non-coplanar multiplanet systems. Simulations include two secularly interacting Jupiter-mass planets mutually inclined by 10 degrees, with the host planet either oblate or spherical. With a spherical host planet, moons within a critical planetocentric distance experience high inclinations and in some cases high eccentricities, while more distant moons orbit stably with low inclinations and eccentricities, as expected. These counter-intuitive dynamical phenomena disappear with an oblate host planet, in which case the moons' Laplace plane transitions from the host planet's equatorial plane to the host planet's precessing orbital plane as their semi-major axes increase, and all moons are dynamically stable with very mild changes in orbits. Direct perturbation from the perturbing planet has been investigated and ruled out as an explanation for the behavior of the innermost satellites, therefore leaving the central star's perturbation as the cause. Instability occurs while the nodal precession of the satellite and the central star (as seen from the host planet's frame) approaches the 1:1 secular resonance. In non-coplanar systems, around a non-oblate planet, the nodal precession of the moon becomes slow and comparable to that of the planet, giving rise to resonant configurations. The above effect needs to be taken into account in setting up numerical simulations.

Detecting Exomoons

On the detection of Exomoons

Authors:

Hippke et al

Abstract:

Despite the discovery of thousands of exoplanets, no exomoons have been detected so far. We test a recently developed method for exomoon search, the "orbital sampling effect" (OSE), using the full exoplanet photometry from the Kepler Space Telescope. The OSE is applied to phase-folded transits, for which we present a framework to detect false positives, and discuss four candidates which pass several of our tests. Using numerical simulations, we inject exomoon signals into real Kepler data and retrieve them, showing that under favorable conditions, exomoons can be found with Kepler and the OSE method. In addition, we super-stack a large sample of Kepler planets to search for the average exomoon OSE and the accompanying increase in noise, the "scatter peak". We find significant exomoon presence for planets with 35d less than P less than 80d, with an average dip per planet of 6+-2ppm, about the radius of Ganymede.

Spectral Classification of ExoPlanet Candidates

The power of low-resolution spectroscopy: On the spectral classification of planet candidates in the ground-based CoRoT follow-up

Authors:

Ammler-von Eiff et al

Abstract:

Planetary transits detected by the CoRoT mission can be mimicked by a low-mass star in orbit around a giant star. Spectral classification helps to identify the giant stars and also early-type stars which are often excluded from further follow-up.

We study the potential and the limitations of low-resolution spectroscopy to improve the photometric spectral types of CoRoT candidates. In particular, we want to study the influence of the signal-to-noise ratio (SNR) of the target spectrum in a quantitative way. We built an own template library and investigate whether a template library from the literature is able to reproduce the classifications. Including previous photometric estimates, we show how the additional spectroscopic information improves the constraints on spectral type.

Low-resolution spectroscopy (R≈1000) of 42 CoRoT targets covering a wide range in SNR (1-437) and of 149 templates was obtained in 2012-2013 with the Nasmyth spectrograph at the Tautenburg 2m telescope. Spectral types have been derived automatically by comparing with the observed template spectra. The classification has been repeated with the external CFLIB library.

The spectral class obtained with the external library agrees within a few sub-classes when the target spectrum has a SNR of about 100 at least. While the photometric spectral type can deviate by an entire spectral class, the photometric luminosity classification is as close as a spectroscopic classification with the external library. A low SNR of the target spectrum limits the attainable accuracy of classification more strongly than the use of external templates or photometry. Furthermore we found that low-resolution reconnaissance spectroscopy ensures that good planet candidates are kept that would otherwise be discarded based on photometric spectral type alone.

MOA-2009-BLG-319 Does not Have More Exoplanets Neptune Size or Larger

Constraint on Additional Planets in Planetary Systems Discovered through the Channel of High-magnification Gravitational Microlensing Events

Authors:

Shin et al

Abstract:

High-magnification gravitational microlensing events provide an important channel of detecting planetary systems with multiple giants located at their birth places. In order to investigate the potential existence of additional planets, we reanalyze the light curves of the eight high-magnification microlensing events for each of which a single planet was previously detected. The analyzed events include OGLE-2005-BLG-071, OGLE-2005-BLG-169, MOA-2007-BLG-400, MOA-2008-BLG-310, MOA-2009-BLG-319, MOA-2009-BLG-387, MOA-2010-BLG-477, and MOA-2011-BLG-293. We find that including an additional planet improves fits with Δχ2 less than80 for seven out of eight analyzed events. For MOA-2009-BLG-319, the improvement is relatively big with Δχ2∼143. From inspection of the fits, we find that the improvement of the fits is attributed to systematics in data. Although no clear evidence of additional planets is found, it is still possible to constrain the existence of additional planets in the parameter space. For this purpose, we construct exclusion diagrams showing the confidence levels excluding the existence of an additional planet as a function of its separation and mass ratio. We also present the exclusion ranges of additional planets with 90\% confidence level for Jupiter, Saturn, and Uranus-mass planets.

Direct Imaging of Exoplanets With Extremely Large Telescopes

Direct Imaging of Exoplanets Without Background Subtraction: Implications for ELTs

Author:

Frazin

Abstract:

The ultra-high contrast capability required to form images of other solar systems is arguably the highest-profile challenge in astronomy today. The current high-contrast imaging efforts all require background subtraction to separate the planetary image from the image of the host star. Background estimation is difficult due to the presence of non-common path aberrations (NCPAs) that change with time. The only major source of information that is not being utilized by current efforts is the random encoding of the planetary image and the NCPAs by the atmosphere on millisecond time-scales. Here, a method that utilizes this information in order to avoid background subtraction altogether is proposed. This new paradigm will allow simultaneous estimation of the time-dependent NCPAs and the planetary image via rigorous statistical inference procedures. These procedures are fully compatible with other information sources, such as diurnal field rotation and spectral diversity. Given the open-ended nature of the background subtraction issues, the ideas explained herein may well the key to imaging habitable planets with Extremely Large Telescopes (ELTs). Fully exploiting the information content of millisecond exposures will require significant design modifications of the ELT wavefront sensors and science camera systems, if ultra-high contrast imaging is to be priority.

What Would Jupiter's Signature be as an Exoplanet?

Jupiter as an exoplanet: UV to NIR transmission spectrum reveals hazes, a Na layer and possibly stratospheric H2O-ice clouds

Authors:

Montañes-Rodriguez et al

Abstract:

Currently, the analysis of transmission spectra is the most successful technique to probe the chemical composition of exoplanet atmospheres. But the accuracy of these measurements is constrained by observational limitations and the diversity of possible atmospheric compositions. Here we show the UV-VIS-IR transmission spectrum of Jupiter, as if it were a transiting exoplanet, obtained by observing one of its satellites, Ganymede, while passing through Jupiter's shadow i.e., during a solar eclipse from Ganymede. The spectrum shows strong extinction due to the presence of clouds (aerosols) and haze in the atmosphere, and strong absorption features from CH4. More interestingly, the comparison with radiative transfer models reveals a spectral signature, which we attribute here to a Jupiter stratospheric layer of crystalline H2O ice. The atomic transitions of Na are also present. These results are relevant for the modeling and interpretation of giant transiting exoplanets. They also open a new technique to explore the atmospheric composition of the upper layers of Jupiter's atmosphere.

Determining the Orbits of Exoplanets and Brown Dwarfs, Including Fomalhaut b

Constraining the orbits of sub-stellar companions imaged over short orbital arcs

Authors:

Pearce et al

Abstract:

Imaging a star's companion at multiple epochs over a short orbital arc provides only four of the six coordinates required for a unique orbital solution. Probability distributions of possible solutions are commonly generated by Monte Carlo (MCMC) analysis, but these are biased by priors and may not probe the full parameter space. We suggest alternative methods to characterise possible orbits, which compliment the MCMC technique. Firstly the allowed ranges of orbital elements are prior-independent, and we provide means to calculate these ranges without numerical analyses. Hence several interesting constraints (including whether a companion even can be bound, its minimum possible semi-major axis and its minimum eccentricity) may be quickly computed using our relations as soon as orbital motion is detected. We also suggest an alternative to posterior probability distributions as a means to present possible orbital elements, namely contour plots of elements as functions of line of sight coordinates. These plots are prior-independent, readily show degeneracies between elements and allow readers to extract orbital solutions themselves. This approach is particularly useful when there are other constraints on the geometry, for example if a companion's orbit is assumed to be aligned with a disc. As examples we apply our methods to several imaged sub-stellar companions including Fomalhaut b, and for the latter object we show how different origin hypotheses affect its possible orbital solutions. We also examine visual companions of A- and G-type main sequence stars in the Washington Double Star Catalogue, and show that ≳50 per cent must be unbound.

Giant ExoPlanets may Only Exist if the Protoplanetary Disks Last Long Enough

Protoplanetary disk lifetimes vs stellar mass and possible implications for giant planet populations

Authors:

Ribas et al

Abstract:

We study the dependence of protoplanetary disk evolution on stellar mass using a large sample of young stellar objects in nearby young star-forming regions. We update the protoplanetary disk fractions presented in our recent work (paper I of this series) derived for 22 nearby (less than 500 pc) associations between 1 and 100 Myr. We use a subsample of 1 428 spectroscopically confirmed members to study the impact of stellar mass on protoplanetary disk evolution. We divide this sample into two stellar mass bins (2 M⊙ boundary) and two age bins (3 Myr boundary), and use infrared excesses over the photospheric emission to classify objects in three groups: protoplanetary disks, evolved disks, and diskless. The homogeneous analysis and bias corrections allow for a statistically significant inter-comparison of the obtained results. We find robust statistical evidence of disk evolution dependence with stellar mass. Our results, combined with previous studies on disk evolution, confirm that protoplanetary disks evolve faster and/or earlier around high-mass (greater than 2 M⊙) stars. We also find a roughly constant level of evolved disks throughout the whole age and stellar mass spectra. We conclude that protoplanetary disk evolution depends on stellar mass. Such a dependence could have important implications for gas giant planet formation and migration, and could contribute to explaining the apparent paucity of hot Jupiters around high-mass stars.

Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data

Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data

Authors:

Brewer et al

Abstract:

Inferring the number of planets N in an exoplanetary system from radial velocity (RV) data is a challenging task. Recently, it has become clear that RV data can contain periodic signals due to stellar activity, which can be difficult to distinguish from planetary signals. However, even doing the inference under a given set of simplifying assumptions (e.g. no stellar activity) can be difficult. It is common for the posterior distribution for the planet parameters, such as orbital periods, to be multimodal and to have other awkward features. In addition, when N is unknown, the marginal likelihood (or evidence) as a function of N is required. Rather than doing separate runs with different trial values of N, we propose an alternative approach using a trans-dimensional Markov Chain Monte Carlo method within Nested Sampling. The posterior distribution for N can be obtained with a single run. We apply the method to ν Oph and Gliese 581, finding moderate evidence for additional signals in ν Oph with periods of 36.11 ± 0.034 days, 75.58 ± 0.80 days, and 1709 ± 183 days; the posterior probability that at least one of these exists is 85%. The results also suggest Gliese 581 hosts many (7-15) "planets" (or other causes of other periodic signals), but only 4-6 have well determined periods. The analysis of both of these datasets shows phase transitions exist which are difficult to negotiate without Nested Sampling.

Using Venus' Solar Transit to Study Exoplanetary Transits

A new view on exoplanet transits: Transit of Venus described using three-dimensional solar atmosphere Stagger-grid simulations

Authors:

Chiavassa et al

Abstract:

Stellar activity and, in particular, convection-related surface structures, potentially cause fluctuations that can affect the transit light curves. Surface convection simulations can help the interpretation of ToV. We used realistic three-dimensional radiative hydrodynamical simulation of the Sun from the Stagger-grid and synthetic images computed with the radiative transfer code Optim3D to provide predictions for the transit of Venus in 2004 observed by the satellite ACRIMSAT. We computed intensity maps from RHD simulation of the Sun and produced synthetic stellar disk image. We computed the light curve and compared it to the ACRIMSAT observations and also to the light curves obtained with solar surface representations carried out using radial profiles with different limb-darkening laws. We also applied the same spherical tile imaging method to the observations of center-to-limb Sun granulation with HINODE. We managed to explain ACRIMSAT observations of 2004 ToV and showed that the granulation pattern causes fluctuations in the transit light curve. We evaluated the contribution of the granulation to the ToV. We showed that the granulation pattern can partially explain the observed discrepancies between models and data. This confirms that the limb-darkening and the granulation pattern simulated in 3D RHD Sun represent well what is imaged by HINODE. In the end, we found that the Venus's aureole contribution during ToV is less intense than the solar photosphere, and thus negligible. Being able to explain consistently the data of 2004 ToV is a new step forward for 3D RHD simulations that are becoming essential for the detection and characterization of exoplanets. They show that the granulation have to be considered as an intrinsic incertitude, due to the stellar variability, on precise measurements of exoplanet transits of, most likely, planets with small diameters.

How Good are Amateur Observatories at Transit Detection of Exoplanets?

Benchmarking the power of amateur observatories for the TTV exoplanets detection

Authors:

Baluev et al

Abstract:

We perform an analysis of ~80000 photometric measurements for the following 10 stars hosting transiting planets: WASP-2, -4, -5, -52, Kelt-1, CoRoT-2, XO-2, TrES-1, HD 189733, GJ 436. Our analysis includes mainly transit lightcurves from the Exoplanet Transit Database, public photometry from the literature, and some proprietary photometry privately supplied by other authors. Half of these lightcurves were obtained by amateurs. From this photometry we derive 306 transit timing measurements, as well as improved planetary transit parameters.

Additionally, for 5 of these 10 stars we present a set of radial velocity measurements obtained from the spectra stored in the HARPS and HARPS-N archives using the HARPS-TERRA pipeline.

Our analysis of these TTV and RV data did not reveal significant hints of additional orbiting bodies in almost all of the cases. In the WASP-4 case, we found hints of marginally significant TTV signals having amplitude 10-20 sec, although their parameters are model-dependent and uncertain, while radial velocities did not reveal statistically significant Doppler signals.

The Outer Disk Around T Cha

High resolution observations of the outer disk around T Cha: the view from ALMA

Authors:

Huelamo et al

Abstract:

T Cha is a young star surrounded by a transitional disk with signatures of planet formation. We have obtained high-resolution and high-sensitivity ALMA observations of T Cha in the CO(3--2), 13CO(3--2), and CS(7--6) emission lines to reveal the spatial distribution of the gaseous disk around the star. In order to study the dust within the disk we have also obtained continuum images at 850μm from the line-free channels. We have spatially resolved the outer disk around T Cha. Using the CO(3-2) emission we derive a radius of ∼230 AU. We also report the detection of the 13CO(3-2) and the CS(7-8) molecular emissions, which show smaller radii than the CO(3-2) detection. The continuum observations at 850μm allow the spatial resolution of the dusty disk, which shows two emission bumps separated by ∼40AU, consistent with the presence of a dust gap in the inner regions of the disk, and an outer radius of ∼80AU. Therefore, T Cha is surrounded by a compact dusty disk and a larger and more diffuse gaseous disk, as previously observed in other young stars. The continuum intensity profiles are different at both sides of the disk suggesting possible dust asymmetries. We derive an inclination of i(deg)=67±5, and a position angle of PA (deg)= 113±6, for both the gas and dust disks. The comparison of the ALMA data with radiative transfer models shows that the gas and dust components can only be simultaneously reproduced when we include a tapered edge prescription for the surface density profile. The best model suggests that most of the disk mass is placed within a radius of R< 50AU. Finally, we derive a dynamical mass for the central object of M∗=1.5±0.2M⊙, comparable to the one estimated with evolutionary models for an age of ∼10Myr.

Protoplanetary Disks in the Young NGC 2024 Cluster

Protoplanetary Disk Masses in the Young NGC 2024 Cluster

Authors:

Mann et al

Abstract:

We present the results from a Submillimeter Array survey of the 887 micron continuum emission from the protoplanetary disks around 95 young stars in the young cluster NGC 2024. Emission was detected from 22 infrared sources, with flux densities from ~5 to 330 mJy; upper limits (at 3sigma) for the other 73 sources range from 3 to 24 mJy. For standard assumptions, the corresponding disk masses range from ~0.003 to 0.2Msolar, with upper limits at 0.002--0.01Msolar. The NGC 2024 sample has a slightly more populated tail at the high end of its disk mass distribution compared to other clusters, but without more information on the nature of the sample hosts it remains unclear if this difference is statistically significant or a superficial selection effect. Unlike in the Orion Trapezium, there is no evidence for a disk mass dependence on the (projected) separation from the massive star IRS2b in the NGC 2024 cluster. We suggest that this is due to either the cluster youth or a comparatively weaker photoionizing radiation field.

The Effects of Shockwaves in Protoplanetary Disks

EVAPORATION OF GRAIN-SURFACE SPECIES BY SHOCK WAVES IN A PROTOPLANETARY DISK

Authors:

Aota et al

Abstract:

Recent Atacama Large Millimeter/submillimeter Array observations of young protostellar objects detected warm SO emission, which could be associated with a forming protostellar disk. In order to investigate if such warm gas can be produced by accretion shock onto the forming disk, we calculate the sputtering and thermal desorption of various grain-surface species in one-dimensional shock waves. We find that thermal desorption is much more efficient than the sputtering in the post-shock region. While H2O can be thermally desorbed, if the accretion velocity is larger than 8 km s–1 with the pre-shock gas number density of 109 cm–3, SO is desorbed if the accretion velocity gsim2 km s–1 and gsim4 km s–1, with the pre-shock density of 109 cm–3 and 108 cm–3, respectively. We also find that the column density of hydrogen nuclei in warm post-shock gas is N warm ~ 1021 cm–2.

Characterizing Transiting ExoPlanet Atmospheres through 2025

Characterizing Transiting Planet Atmospheres through 2025

Authors:

Cowan et al

Abstract:

We have only been able to comprehensively characterize the atmospheres of a handful of transiting planets, because most orbit faint stars. TESS will discover transiting planets orbiting the brightest stars, enabling, in principle, an atmospheric survey of 10^2 to 10^3 bright hot Jupiters and warm sub-Neptunes. Uniform observations of such a statistically significant sample would provide leverage to understand---and learn from---the diversity of short-period planets. We argue that the best way to maximize the scientific returns of TESS is with a follow-up space mission consisting of a ~1 m telescope with an optical--NIR spectrograph: it could measure molecular absorption for non-terrestrial planets, as well as eclipses and phase variations for the hottest jovians. Such a mission could observe up to 10^3 transits per year, thus enabling it to survey a large fraction of the bright (J less than 11) TESS planets. JWST could be used to perform detailed atmospheric characterization of the most interesting transiting targets (transit, eclipse, and---when possible---phase-resolved spectroscopy). TESS is also expected to discover a few temperate terrestrial planets transiting nearby M-Dwarfs. Characterizing these worlds will be time-intensive: JWST will need months to provide tantalizing constraints on the presence of an atmosphere, planetary rotational state, clouds, and greenhouse gases. Future flagship missions should be designed to provide better constraints on the habitability of M-Dwarf temperate terrestrial planets.

BD+49 828b, HD 95127b & HD 216536b: Three Giant ExoPlanets Orbting Red Giants

Three red giants with substellar-mass companions

Authors:

Niedzielski et al

Abstract:

We present three giant stars from the ongoing Penn State-Toru\'n Planet Search with the Hobby-Eberly Telescope, which exhibit radial velocity variations that point to a presence of planetary --mass companions around them. BD+49 828 is a M=1.52±0.22 M⊙ K0 giant with a msini=1.6+0.4−0.2 MJ minimum mass companion in a=4.2+0.32−0.2 AU (2590+300−180d), e=0.35+0.24−0.10 orbit. HD 95127, a logL/L⊙=2.28±0.38, R=20±9 R⊙, M=1.20±0.22 M⊙ K0 giant has a msini=5.01+0.61−0.44 MJ minimum mass companion in a=1.28+0.01−0.01 AU (482+5−5d), e=0.11+0.15−0.06 orbit. Finally, HD 216536, is a M=1.36±0.38 M⊙ K0 giant with a msini=1.47+0.20−0.12 MJ minimum mass companion in a=0.609+0.002−0.002 AU (148.6+0.7−0.7d), e=0.38+0.12−0.10 orbit. Both, HD 95127 b and HD 216536 b in their compact orbits, are very close to the engulfment zone and hence prone to ingestion in the near future. BD+49 828 b is among the longest period planets detected with the radial velocity technique until now and it will remain unaffected by stellar evolution up to a very late stage of its host. We discuss general properties of planetary systems around evolved stars and planet survivability using existing data on exoplanets in more detail.

Can Planetoids Impacting on White Dwarfs Cause Type Ia Supernova?

Death by Dynamics: Planetoid-Induced Explosions on White Dwarfs

Authors:

Di Stefano et al

Abstract:

At intervals as short as ten thousand years, each white dwarf (WD) passes within a solar radius of a planetoid, i.e., a comet, asteroid, or planet. Gravitational tidal forces tear the planetoid apart; its metal-rich debris falls onto the WD, enriching the atmosphere. A third of WDs exhibit atmospheric "pollution". For roughly every hundred planetoid disruptions, a planetoid collides with a WD. We simulate a small number of collisions, in which "death-by-dynamics" refers to the fate of the planetoid. We also compute the energies and likely durations of a broad sample of collision events, and identify detection strategies at optical and X-ray wavelengths. Collisions with the most massive planetoids can be detected in external galaxies. Some may trigger nuclear burning. If one in ∼107−108 of WD-planetoid collisions creates the conditions needed for a Type Ia supernova (SN~Ia), "death-by-dynamics" would also refer to the fate of the WD, and could provide a novel channel for the production of SN~Ia. We consider the circumstances under which the rate of SNe~Ia can be increased by interactions with planetoids.

AAAS Alternate Biochemistries Panel Tomorrow @ 3 PM in San Jose

Ideas about directing evolution of life forms on Earth and finding life on other planets are rapidly morphing from science-fiction fantasy into mainstream science, says David Lynn, a chemist at Emory University.

"These areas of science are rapidly coming of age because of our increasing knowledge and advancing technology. It's an exciting time. We're on the threshold of answering fundamental questions including: What is life? Are there forms of life that we haven't even yet imagined? Are we alone in the universe?"

A panel discussion, "Searching for Alternative Chemistries of Life on Earth and Throughout the Universe," is set for Friday, February 13, at 3 pm, during the annual meeting of the American Association for the Advancement of Science (AAAS) in San Jose. Lynn co-organized the panel with Jay Goodwin, a senior Emory research fellow and an AAAS Science and Technology Policy Fellow.

link.

How Fast do Hot Jupiters Evaporate?

High-energy irradiation and mass loss rates of hot Jupiters in the solar neighborhood

Authors:

Salz et al

Abstract:

Giant gas planets in close proximity to their host stars experience strong irradiation. In extreme cases photoevaporation causes a transonic, planetary wind and the persistent mass loss can possibly affect the planetary evolution. We have identified nine hot Jupiter systems in the vicinity of the Sun, in which expanded planetary atmospheres should be detectable through Lyman alpha transit spectroscopy according to predictions. We use X-ray observations with Chandra and XMM-Newton of seven of these targets to derive the high-energy irradiation level of the planetary atmospheres and the resulting mass loss rates. We further derive improved Lyman alpha luminosity estimates for the host stars including interstellar absorption. According to our estimates WASP-80 b, WASP-77 b, and WASP-43 b experience the strongest mass loss rates, exceeding the mass loss rate of HD 209458 b, where an expanded atmosphere has been confirmed. Furthermore, seven out of nine targets might be amenable to Lyman alpha transit spectroscopy. Finally, we check the possibility of angular momentum transfer from the hot Jupiters to the host stars in the three binary systems among our sample, but find only weak indications for increased stellar rotation periods of WASP-77 and HAT-P-20.

The Cause of Hot Jupiter Inflation

Hot-Jupiter Inflation due to Deep Energy Deposition

Authors:

Ginzburg et al

Abstract:

Some extrasolar giant planets in close orbits---"hot Jupiters"---exhibit larger radii than that of a passively cooling planet. The extreme irradiation Leq these hot Jupiters receive from their close in stars creates a thick isothermal layer in their envelopes, which slows down their convective cooling, allowing them to retain their inflated size for longer. This is yet insufficient to explain the observed sizes of the most inflated planets. Some models invoke an additional power source, deposited deep in the planet's envelope. Here we present an analytical model for the cooling of such irradiated, and internally heated gas giants. We show that a power source Ldep, deposited at an optical depth τdep, creates an exterior convective region, between optical depths Leq/Ldep and τdep, beyond which a thicker isothermal layer exists, which in extreme cases may extend to the center of the planet. This convective layer, which occurs only for Ldepτdep>Leq, further delays the cooling of the planet. Such a planet is equivalent to a planet irradiated with Leq(1+Ldepτdep/Leq)β, where β≈0.35 is an effective power-law index describing the radiative energy density as function of the optical depth for a convective planet U∝τβ. Our simple analytical model reproduces the main trends found in previous numerical works, and provides an intuitive understanding. We derive scaling laws for the cooling rate of the planet, its central temperature, and radius. These scaling laws can be used to estimate the effects of tidal or Ohmic dissipation, wind shocks, or any other mechanism involving energy deposition, on sizes of hot Jupiters.

KOI-614b, KOI-206b, and KOI-680b: A Warm Jupiter and two Highly Inflated Gas Giants

SOPHIE velocimetry of Kepler transit candidates. XV. KOI-614b, KOI-206b, and KOI-680b: a massive warm Jupiter orbiting a G0 metallic dwarf and two highly inflated planets with a distant companion around evolved F-type stars

Authors:

Almenara et al

Abstract:

We report the validation and characterization of three new transiting exoplanets using SOPHIE radial velocities: KOI-614b, KOI-206b, and KOI-680b. KOI-614b has a mass of 2.86±0.35 MJup and a radius of 1.13+0.26−0.18 RJup, and it orbits a G0, metallic ([Fe/H]=0.35±0.15) dwarf in 12.9 days. Its mass and radius are familiar and compatible with standard planetary evolution models, so it is one of the few known transiting planets in this mass range to have an orbital period over ten days. With an equilibrium temperature of Teq=1000±45 K, this places KOI-614b at the transition between what is usually referred to as "hot" and "warm" Jupiters. KOI-206b has a mass of 2.82±0.52 MJup and a radius of 1.45±0.16 RJup, and it orbits a slightly evolved F7-type star in a 5.3-day orbit. It is a massive inflated hot Jupiter that is particularly challenging for planetary models because it requires unusually large amounts of additional dissipated energy in the planet. On the other hand, KOI-680b has a much lower mass of 0.84±0.15 MJup and requires less extra-dissipation to explain its uncommonly large radius of 1.99±0.18 RJup. It is one of the biggest transiting planets characterized so far, and it orbits a subgiant F9-star well on its way to the red giant stage, with an orbital period of 8.6 days. With host stars of masses of 1.46±0.17 M⊙ and 1.54±0.09 M⊙, respectively, KOI-206b, and KOI-680b are interesting objects for theories of formation and survival of short-period planets around stars more massive than the Sun. For those two targets, we also find signs of a possible distant additional companion in the system.

Fomalhaut b a False Positive?

The companion candidate near Fomalhaut - a background neutron star?

Authors:

Neuhaeuser et al

Abstract:

The directly detected planetary mass companion candidate close to the young, nearby star Fomalhaut is a subject of intense discussion. While the detection of common proper motion led to the interpretation as Jovian-mass companion, later non-detections in the infrared raised doubts. Recent astrometric measurements indicate a belt crossing or highly eccentric orbit for the object, if a companion, making the planetary interpretation potentially even more problematic. In this study we discuss the possibility of Fomalhaut\,b being a background object with a high proper motion. By analysing the available photometric and astrometric data of the object, we show that they are fully consistent with a neutron star: Neutron stars are faint, hot (blue), and fast moving. Neutron stars with an effective temperature of the whole surface area being 112,000 K to 126,500 K (with small to negligible extinction) at a distance of roughly 11 pc (best fit) would be consistent with all observables, namely with the photometric detections in the optical, with the upper limits in the infrared and X-rays, as well as with the astrometry (consistent with a distances of 11 pc or more and high proper motion as typical for neutron stars) as well as with non-detection of pulsation (not beamed). We consider the probability of finding an unrelated object or even a neutron star nearby and mostly co-aligned in proper motion with Fomalhaut A and come to the conclusion that this is definitely well possible.

Doubts Cast on hot Neptune Gliese 436b's Multi Exoplanetary System

On the GJ 436 planetary system

Authors:

Maciejewski et al

Abstract:

The GJ 436 system contains a transiting planet GJ 436 b which is a hot analogue of Neptune on an eccentric orbit. Recently, two additional transiting sub-Earth planets have been postulated in the literature. We observed three transits of GJ 436 b over the course of 3 years using two-meter class telescopes, each with a photometric precision better than one millimagnitude. We studied system dynamics based on the existence of the additional planets. We redetermined system parameters, which were in agreement with those found in the literature. We refined the orbital period of GJ 436 b and found no evidence of transit timing variations. The orbital motion of the GJ 436 c planet candidate was found to be significantly affected by the planet b with variations in transit times at a level of 20 minutes. As the orbital period of the GJ 436 d planet candidate remains unknown, our numerical experiments rule out orbits in low-order resonances with GJ 436 b. The GJ 436 system with the hot Neptune and additional two Earth-like planets, if confirmed, would be an important laboratory for studies of formation and evolution of planetary systems.

Revisiting the GJ436 system

Revisiting Spitzer transit observations with Independent Component Analysis: new results for the GJ436 system

Authors:

Morello et al

Abstract:

We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 {\mu}m of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to extract conclusively an atmospheric signal; those analyses also reported discrepant results, hence the necessity of this reanalysis. The method we used has been proposed in Morello et al. (2014) to analyze 3.6 {\mu}m transit light-curves of the hot Jupiter HD189733b; it performes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series read by individual pixels, from the same photometric observation. Our method only assumes the independence of instrumental and astrophysical signals, and therefore guarantees a higher degree of objectivity compared to parametric detrending techniques published in the literature. The datasets we analyzed in this paper represent a more challenging test compared to the previous ones. Contrary to previous results reported in the literature, our results (1) do not support any detectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable at the level 10e-4 in the IR, and (3) the transit depth measurements at the two wavelengths are consistent within 1{\sigma}. We also (4) detect a possible transit duration variation (TDV) of 80 s (2 {\sigma} significance level), that has not been pointed out in the literature, and (5) confirm no transit timing variations (TTVs) greater than 30 s.