Source author record

Daniel Apai

Daniel Apai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

An Integrative Analysis of the Rich Planetary System of the Nearby Star e Eridani: Ideal Targets For Exoplanet Imaging and Biosignature Searches

e Eridani, the fifth-closest Sun-like star, hosts at least three planets and could possibly harbor more. However, the veracity of the planet candidates in the system and its full planetary architecture remain unknown. Here we analyze the planetary architecture of e Eridani via DYNAMITE, a method providing an integrative assessment of the system architecture (and possibly yet-undetected planets) by combining statistical, exoplanet-population level knowledge with incomplete but specific information available on the system. DYNAMITE predicts the most likely location of an additional planet in the system based on the Kepler population demographic information from more than 2000 planets. Additionally, we analyze the dynamical stability of e Eridani system via N-body simulations. Our DYNAMITE and dynamical stability analyses provide support for planet candidates g, c, and f, and also predict one additional planet candidate with an orbital period between 549 -- 733 days, in the habitable zone of the system. We find that planet candidate f, if it exists, would also lie in the habitable zone. Our dynamical stability analysis also shows that the e Eridani planetary eccentricities, as reported, do not allow for a stable system, suggesting that they are lower. We introduce a new statistical approach for estimating the equilibrium and surface temperatures of exoplanets, based on a prior on the planetary albedo distribution. e Eridani is a rich planetary system with a possibility of containing two potentially habitable planets, and its vicinity to our Solar System makes it an important target for future imaging studies and biosignature searches.

preprint2022arXiv

The Scorpion Planet Survey: Wide-Orbit Giant Planets Around Young A-type Stars

The first directly imaged exoplanets indicated that wide-orbit giant planets could be more common around A-type stars. However, the relatively small number of nearby A-stars has limited the precision of exoplanet demographics studies to $\gtrsim$10%. We aim to constrain the frequency of wide-orbit giant planets around A-stars using the VLT/SPHERE extreme adaptive optics system, which enables targeting $\gtrsim$100 A-stars between 100$-$200 pc. We present the results of a survey of 84 A-stars within the nearby $\sim$5$-$17 Myr-old Sco OB2 association. The survey detected three companions$-$one of which is a new discovery (HIP75056Ab), whereas the other two (HD 95086b and HIP65426b) are now-known planets that were included without a priori knowledge of their existence. We assessed the image sensitivity and observational biases with injection and recovery tests combined with Monte Carlo simulations to place constraints on the underlying demographics. We measure a decreasing frequency of giant planets with increasing separation, with measured values falling between 10$-$2% for separations of 30$-$100 au, and 95% confidence-level (CL) upper limits of $\lesssim$45$-$8% for planets on 30$-$100 au orbits, and $\lesssim$5% between 200$-$500 au. These values are in excellent agreement with recent surveys of A-stars in the solar neighborhood$-$supporting findings that giant planets at $\lesssim$100 au are more frequent around A-stars than around solar-type hosts. Finally, the relatively low occurrence rate of super-Jupiters on wide orbits, the positive correlation with stellar mass, and the inverse correlation with orbital separation are consistent with core accretion being their dominant formation mechanism.

preprint2021arXiv

TESS Observations of the Luhman 16AB Brown Dwarf System: Rotational Periods, Lightcurve Evolution, and Zonal Circulation

Brown dwarfs were recently found to display rotational modulations, commonly attributed to cloud cover of varying thickness, possibly modulated by planetary-scale waves. However, the long-term, continuous, high-precision monitoring data to test this hypothesis for more objects is lacking. By applying our novel photometric approach to TESS data, we extract a high-precision lightcurve of the closest brown dwarfs, which form the binary system Luhman 16AB. Our observations, that cover about 100 rotations of Luhman 16B, display continuous lightcurve evolution. The periodogram analysis shows that the rotational period of the component that dominates the lightcurve is 5.28 h. We also find evidence for periods of 2.5 h, 6.94 h, and 90.8 h. We show that the 2.5 h and 5.28 h periods emerge from Luhman 16B and that they consist of multiple, slightly shifted peaks, revealing the presence of high-speed jets and zonal circulation in this object. We find that the lightcurve evolution is well fit by the planetary-scale waves model, further supporting this interpretation. We argue that the 6.94 h peak is likely the rotation period of Luhman 16A. By comparing the rotational periods to observed v sin(i) measurements, we show that the two brown dwarfs are viewed at angles close to their equatorial planes. We also describe a long-period (P~91 h) evolution in the lightcurve, which we propose emerges from the vortex-dominated polar regions. Our study paves the way toward direct comparisons of the predictions of global circulation models to observations via periodogram analysis.

preprint2020arXiv

Earths in Other Solar Systems N-body simulations: the Role of Orbital Damping in Reproducing the Kepler Planetary Systems

The population of exoplanetary systems detected by Kepler provides opportunities to refine our understanding of planet formation. Unraveling the conditions needed to produce the observed exoplanets will sallow us to make informed predictions as to where habitable worlds exist within the galaxy. In this paper, we examine using N-body simulations how the properties of planetary systems are determined during the final stages of assembly. While accretion is a chaotic process, trends in the ensemble properties of planetary systems provide a memory of the initial distribution of solid mass around a star prior to accretion. We also use EPOS, the Exoplanet Population Observation Simulator, to account for detection biases and show that different accretion scenarios can be distinguished from observations of the Kepler systems. We show that the period of the innermost planet, the ratio of orbital periods of adjacent planets, and masses of the planets are determined by the total mass and radial distribution of embryos and planetesimals at the beginning of accretion. In general, some amount of orbital damping, either via planetesimals or gas, during accretion is needed to match the whole population of exoplanets. Surprisingly, all simulated planetary systems have planets that are similar in size, showing that the "peas in a pod" pattern can be consistent with both a giant impact scenario and a planet migration scenario. The inclusion of material at distances larger than what Kepler observes has a profound impact on the observed planetary architectures, and thus on the formation and delivery of volatiles to possible habitable worlds.

preprint2020arXiv

EDEN: Sensitivity Analysis and Transiting Planet Detection Limits for Nearby Late Red Dwarfs

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.

preprint2020arXiv

First Images of the Protoplanetary Disk Around PDS 201

Scattered light imaging has revealed nearly a dozen circumstellar disks around young Herbig Ae/Be stars$-$enabling studies of structures in the upper disk layers as potential signs of on-going planet formation. We present the first images of the disk around the variable Herbig Ae star PDS 201 (V* V351 Ori), and an analysis of the images and spectral energy distribution through 3D Monte-Carlo radiative transfer simulations and forward modelling. The disk is detected in three datasets with LBTI/LMIRCam at the LBT, including direct observations in the $Ks$ and $L'$ filters, and an $L'$ observation with the 360$^\circ$ vector apodizing phase plate coronagraph. The scattered light disk extends to a very large radius of $\sim$250 au, which places it among the largest of such disks. Exterior to the disk, we establish detection limits on substellar companions down to $\sim$5 M$_{Jup}$ at $\gtrsim$1.5" ($\gtrsim$500 au), assuming the Baraffe et al. (2015) models. The images show a radial gap extending to $\sim$0.4" ($\sim$140 au at a distance of 340 pc) that is also evident in the spectral energy distribution. The large gap is a possible signpost of multiple high-mass giant planets at orbital distances ($\sim$60-100 au) that are unusually massive and widely-separated compared to those of planet populations previously inferred from protoplanetary disk substructures.

preprint2020arXiv

Helios-r2 -- A new Bayesian, open-source retrieval model for brown dwarfs and exoplanet atmospheres

We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.

preprint2020arXiv

Hidden Worlds: Dynamical Architecture Predictions of Undetected Planets in Multi-planet Systems and Applications to TESS Systems

Multi-planet systems produce a wealth of information for exoplanet science, but our understanding of planetary architectures is incomplete. Probing these systems further will provide insight into orbital architectures and formation pathways. Here we present a model to predict previously undetected planets in these systems via population statistics. The model considers both transiting and non-transiting planets, and can test the addition of more than one planet. Our tests show the model's orbital period predictions are robust to perturbations in system architectures on the order of a few percent, much larger than current uncertainties. Applying it to the multi-planet systems from TESS provides a prioritized list of targets, based on predicted transit depth and probability, for archival searches and for guiding ground-based follow-up observations hunting for hidden planets.

preprint2019arXiv

The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler

The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multi-planet systems, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100 seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/epos

preprint2016arXiv

Cloud Atlas: Discovery of Patchy Clouds and High-amplitude Rotational Modulations In a Young, Extremely Red L-type Brown Dwarf

Condensate clouds fundamentally impact the atmospheric structure and spectra of exoplanets and brown dwarfs but the connections between surface gravity, cloud structure, dust in the upper atmosphere, and the red colors of some brown dwarfs remain poorly understood. Rotational modulations enable the study of different clouds in the same atmosphere, thereby providing a method to isolate the effects of clouds. Here we present the discovery of high peak-to-peak amplitude (8%) rotational modulations in a low-gravity, extremely red (J-Ks=2.55) L6 dwarf WISEP J004701.06+680352.1 (W0047). Using the Hubble Space Telescope (HST) time-resolved grism spectroscopy we find a best-fit rotational period (13.20$\pm$0.14 hours) with a larger amplitude at 1.1 micron than at 1.7 micron. This is the third largest near-infrared variability amplitude measured in a brown dwarf, demonstrating that large-amplitude variations are not limited to the L/T transition but are present in some extremely red L-type dwarfs. We report a tentative trend between the wavelength dependence of relative amplitude, possibly proxy for small dust grains lofted in the upper atmosphere, and the likelihood of large-amplitude variability. By assuming forsterite as haze particle, we successfully explain the wavelength dependent amplitude with submicron-sized haze particles sizes of around 0.4 μm. W0047 links the earlier spectral and later spectral type brown dwarfs in which rotational modulations have been observed, the large amplitude variations in this object make this a benchmark brown dwarf for the study of cloud properties close to the L/T transition.

preprint2016arXiv

High-Cadence, High-Contrast Imaging for Exoplanet Mapping: Observations of the HR 8799 Planets with VLT/SPHERE Satellite Spot-Corrected Relative Photometry

Time-resolved photometry is an important new probe of the physics of condensate clouds in extrasolar planets and brown dwarfs. Extreme adaptive optics systems can directly image planets, but precise brightness measurements are challenging. We present VLT/SPHERE high-contrast, time-resolved broad H-band near-infrared photometry for four exoplanets in the HR 8799 system, sampling changes from night to night over five nights with relatively short integrations. The photospheres of these four planets are often modeled by patchy clouds and may show large-amplitude rotational brightness modulations. Our observations provide high-quality images of the system. We present a detailed performance analysis of different data analysis approaches to accurately measure the relative brightnesses of the four exoplanets. We explore the information in satellite spots and demonstrate their use as a proxy for image quality. While the brightness variations of the satellite spots are strongly correlated, we also identify a second-order anti-correlation pattern between the different spots. Our study finds that PCA-based KLIP reduction with satellite spot-modulated artificial planet-injection based photometry (SMAP) leads to a significant (~3x) gain in photometric accuracy over standard aperture-based photometry and reaches 0.1 mag per point accuracy for our dataset, the signal-to-noise of which is limited by small field rotation. Relative planet-to-planet photometry can be compared be- tween nights, enabling observations spanning multiple nights to probe variability. Recent high-quality relative H-band photometry of the b-c planet pair agree to about 1%.

preprint2016arXiv

Maps of Evolving Cloud Structures in Luhman 16AB from HST Time-Resolved Spectroscopy

WISE J104915.57-531906.1 is the nearest brown dwarf binary to our Solar system, consisting of two brown dwarfs in the L/T transition: Luhman 16A & B. In this paper we present the first map of Luhman 16A, and maps of Luhman 16B for two epochs. Our maps were created by applying Aeolus, a Markov-Chain Monte Carlo code that maps the top-of-the-atmosphere structure of brown dwarf and other ultracool atmospheres, to light curves of Luhman 16A & B using the Hubble Space Telescope's G141 and G102 grisms. Aeolus retrieved three or four spots in the top-of-the-atmosphere of Luhman 16A & B, with a surface coverage of 19%-32% (depending on an assumed rotational period of 5 hr or 8 hr) or 21%-38.5% (depending on the observational epoch) respectively. The brightness temperature of the spots of the best-fit models was ~200 K hotter than the background top-of-the-atmosphere. We compared our Luhman 16B map with the only previously published map. Interestingly, our map contained a large, cooler (DT~51 K) than the background top-of-the-atmosphere spot that lay at low latitudes, in agreement with the previous Luhman 16B map. Finally, we report the detection of a feature reappearing in Luhman 16B light curves that are separated by tens of hundreds of rotations from each other. We speculate this feature is related to top-of-the-atmosphere structures of Luhman 16B.

preprint2015arXiv

A stellar-mass-dependent drop in planet occurrence rates

The Kepler Spacecraft has discovered a large number of planets up to one-year periods and down to terrestrial sizes. While the majority of the target stars are main-sequence dwarfs of spectral type F, G, and K, Kepler covers stars with effective temperature as low as 2500 K, which corresponds to M stars. These cooler stars allow characterization of small planets near the habitable zone, yet it is not clear if this population is representative of that around FGK stars. In this paper, we calculate the occurrence of planets around stars of different spectral types as a function of planet radius and distance from the star, and show that they are significantly different from each other. We further identify two trends: First, the occurrence of Earth to Neptune-sized planets is successively higher toward later spectral types at all orbital periods probed by Kepler; Planets around M stars occur twice as frequently as around G stars, and thrice as frequently as around F stars. Second, a drop in planet occurrence is evident at all spectral types inward of a 10 day orbital period, with a plateau further out. By assigning to each spectral type a median stellar mass, we show that the distance from the star where this drop occurs is stellar mass dependent, and scales with semi-major axis as the cube root of stellar mass. By comparing different mechanisms of planet formation, trapping and destruction, we find that this scaling best matches the location of the pre-main-sequence co-rotation radius, indicating efficient trapping of migrating planets or planetary building blocks. These results demonstrate the stellar-mass dependence of the planet population, both in terms of occurrence rate and of orbital distribution. The prominent stellar-mass dependence of the inner boundary of the planet population shows that the formation or migration of planets is sensitive to the stellar parameters.

preprint2015arXiv

Aeolus: A Markov--Chain Monte Carlo code for mapping ultracool atmospheres. An application on Jupiter and brown dwarf HST light curves

Deducing the cloud cover and its temporal evolution from the observed planetary spectra and phase curves can give us major insight into the atmospheric dynamics. In this paper, we present Aeolus, a Markov-Chain Monte Carlo code that maps the structure of brown dwarf and other ultracool atmospheres. We validated Aeolus on a set of unique Jupiter Hubble Space Telescope (HST) light curves. Aeolus accurately retrieves the properties of the major features of the jovian atmosphere such as the Great Red Spot and a major 5um hot spot. Aeolus is the first mapping code validated on actual observations of a giant planet over a full rotational period. For this study, we applied Aeolus to J and H-bands HST light curves of 2MASSJ21392676+0220226 and 2MASSJ0136565+093347. Aeolus retrieves three spots at the top-of-the-atmosphere (per observational wavelength) of these two brown dwarfs, with a surface coverage of 21+-3% and 20.3+-1.5% respectively. The Jupiter HST light curves will be publicly available via ADS/VIZIR.

preprint2015arXiv

An Increase in the Mass of Planetary Systems around Lower-Mass Stars

Trends in the planet population with host star mass provide an avenue to constrain planet formation theories. We derive the planet radius distribution function for Kepler stars of different spectral types, sampling a range in host star masses. We find that M dwarf stars have 3.5 times more small planets (1.0-2.8 R_Earth) than main-sequence FGK stars, but two times fewer Neptune-sized and larger planets (>2.8 R_Earth). We find no systematic trend in the planet size distribution between spectral types F, G, and K to explain the increasing occurrence rates. Taking into account the mass-radius relationship and heavy-element mass of observed exoplanets, and assuming those are independent of spectral type, we derive the inventory of the heavy-element mass locked up in exoplanets at short orbits. The overall higher planet occurrence rates around M stars are not consistent with the redistribution of the same mass into more, smaller planets. At the orbital periods and planet radii where Kepler observations are complete for all spectral types, the average heavy-element mass locked up in exoplanets increases roughly inversely with stellar mass from 4 M_Earth in F stars to 5 M_Earth in G and K stars to 7 M_Earth in M stars. This trend stands in stark contrast with observed protoplanetary disk masses that decrease towards lower mass stars, and provides a challenge for current planet formation models. Neither models of in situ formation nor migration of fully-formed planets are consistent with these results. Instead, these results are indicative of large-scale inward migration of planetary building blocks --- either through type-I migration or radial drift of dust grains --- that is more efficient for lower mass stars, but does not result in significantly larger or smaller planets.

preprint2015arXiv

Discovery of a Two-Armed Spiral Structure in the Gapped Disk in HD 100453

We present VLT/SPHERE adaptive optics imaging in Y$-$, J$-$, H$-$, and K-bands of the HD 100453 system and the discovery of a two-armed spiral structure in a disk extending to 0.37" ($\sim$42 AU) from the star, with highly symmetric arms to the Northeast and Southwest. Inside of the spiral arms, we resolve a ring of emission from 0.18"-0.25" ($\sim$21-29 AU). By assuming that the ring is intrinsically circular we estimate an inclination of $\sim$34$^{o}$ from face-on. We detect dark crescents on opposite sides (NW and SE) which begin at 0.18" and continue to radii smaller than our inner working angle of 0.15", which we interpret as the signature of a gap at $\lesssim$21 AU that has likely been cleared by forming planets. We also detect the $\sim$120 AU companion HD 100453 B, and by comparing our data to 2003 HST/ACS and VLT/NACO images we estimate an orbital period of $\sim$850 yr. We discuss what implications the discovery of the spiral arms and finer structures of the disk may have on our understanding of the possible planetary system in HD 100453, and how the morphology of this disk compares to other related objects.

preprint2015arXiv

Discovery of an Edge-on Debris Disk with a Dust Ring and an Outer Disk Wing-tilt Asymmetry

Using VLT/SPHERE near-infrared dual-band imaging and integral field spectroscopy we discovered an edge-on debris disk around the 17\,Myr old A-type member of the Scorpius-Centaurus OB association HD 110058. The edge-on disk can be traced to about 0.6" or 65 AU projected separation. In its northern and southern wings, the disk shows at all wavelengths two prominent, bright and symmetrically placed knots at 0.3" or 32 AU from the star. We interpret these knots as a ring of planetesimals whose collisions may produce most of the dust observed in the disk. We find no evidence for a bow in the disk, but we identify a pair of symmetric, hook-like features in both wings. Based on similar features in the Beta Pictoris disk we propose that this wing-tilt asymmetry traces either an outer planetesimal belt that is inclined with respect to the disk midplane or radiation-pressure-driven dust blown out from a yet unseen, inner belt which is inclined with respect to the disk midplane. The misaligned inner or outer disk may be a result of interaction with a yet unseen planet. Overall, the disk geometry resembles the nearby disk around Beta Pictoris, albeit seen at smaller radial scales.

preprint2015arXiv

Discovery of Rotational Modulations in the Planetary-Mass Companion 2M1207b: Intermediate Rotation Period and Heterogeneous Clouds in a Low Gravity Atmosphere

Rotational modulations of brown dwarfs have recently provided powerful constraints on the properties of ultra-cool atmospheres, including longitudinal and vertical cloud structures and cloud evolution. Furthermore, periodic light curves directly probe the rotational periods of ultra-cool objects. We present here, for the first time, time-resolved high-precision photometric measurements of a planetary-mass companion, 2M1207b. We observed the binary system with HST/WFC3 in two bands and with two spacecraft roll angles. Using point spread function-based photometry, we reach a nearly photon-noise limited accuracy for both the primary and the secondary. While the primary is consistent with a flat light curve, the secondary shows modulations that are clearly detected in the combined light curve as well as in different subsets of the data. The amplitudes are 1.36% in the F125W and 0.78% in the F160W filters, respectively. By fitting sine waves to the light curves, we find a consistent period of $10.7^{+1.2}_{-0.6}$ hours and similar phases in both bands. The J- and H-band amplitude ratio of 2M1207b is very similar to a field brown dwarf that has identical spectral type but different J-H color. Importantly, our study also measures, for the first time, the rotation period for a directly imaged extra-solar planetary-mass companion.

preprint2015arXiv

The Inner Disk Structure, Disk-Planet Interactions, and Temporal Evolution in the Beta Pictoris System: A Two-Epoch HST/STIS Coronagraphic Study

We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.

preprint2015arXiv

Volatile Delivery to Planets from Water-rich Planetesimals around Low Mass Stars

Most models of volatile delivery to accreting terrestrial planets assume that the carriers for water are similar in water content to the carbonaceous chondrites in our Solar System. Here we suggest that the water content of primitive bodies in many planetary systems may actually be much higher, as carbonaceous chondrites have lost some of their original water due to heating from short-lived radioisotopes that drove parent body alteration. Using N-body simulations, we explore how planetary accretion would be different if bodies beyond the water line contained a water mass fraction consistent with chemical equilibrium calculations, and more similar to comets, as opposed to the more traditional water-depleted values. We apply this model to consider planet formation around stars of different masses and identify trends in the properties of Habitable Zone planets and planetary system architecture which could be tested by ongoing exoplanet census data collection. Comparison of such data with the model predicted trends will serve to evaluate how well the N-body simulations and the initial conditions used in studies of planetary accretion can be used to understand this stage of planet formation.

preprint2014arXiv

Cloud structure of brown dwarfs from spectroscopic variability observations

Recent discoveries of variable brown dwarfs have provided us with a new window into their three-dimensional cloud structure. The highest variables are found at the L/T transition, where the cloud cover is thought to break up, but variability has been found to occur also for both cloudy L dwarfs and (mostly) cloud-free mid T dwarfs. We summarize results from recent HST programs measuring the spectral variability of brown dwarfs in the near-infrared and compare to results from ground-based programs. We discuss the patchy cloud structure of L/T transition objects, for which it is becoming increasingly certain that the variability does not arise from cloud holes into the deep hot regions but from varying cloud thickness. We present a new patchy cloud model to explain the spectral variability of 2MASSJ21392676+0220226. We also discuss the curious multi-wavelength variability behavior of the recently discovered very nearby early T dwarf WISE J104915.57-531906.1B (Luhman 16B) and the mid T dwarf 2MASS J22282889-431026.

preprint2014arXiv

Cloud structure of the nearest brown dwarfs: Spectroscopic variability of Luhman 16AB from the Hubble Space Telescope

The binary brown dwarf WISE J104915.57$-$531906.1 (also Luhman 16AB), composed of a late L and early T dwarf, is a prototypical L/T transition flux reversal binary located at only 2 pc distance. Luhman 16B is a known variable whose light curves evolve rapidly. We present spatially resolved spectroscopic time-series of Luhman 16A and B covering 6.5 h using HST/WFC3 at 1.1 to 1.66 $μ$m. The small, count-dependent variability of Luhman 16A at the beginning of the observations likely stems from instrumental systematics; Luhman 16A appears non-variable above $\approx$0.4%. Its spectrum is well fit by a single cloud layer with intermediate cloud thickness (f_sed=2, Teff=1200 K). Luhman 16B varies at all wavelengths with peak-to-valley amplitudes of 7-11%. The amplitude and light curve shape changes over only one rotation period. The lowest relative amplitude is found in the deep water absorption band at 1.4 $μ$m, otherwise it mostly decreases gradually from the blue to the red edge of the spectrum. This is very similar to the other two known highly variable early T dwarfs. A two-component cloud model accounts for most of the variability, although small deviations are seen in the water absorption band. We fit the mean spectrum and relative amplitudes with a linear combination of two models of a warm, thinner cloud (Teff=1300 K, fsed=3) and a cooler, thicker cloud (Teff=1000-1100 K, f_sed=1), assuming out-of-equilibrium atmospheric chemistry. A cloud as for Luhman 16A but with holes cannot reproduce the variability of Luhman 16B, indicating more complex cloud evolution through the L/T transition. The projected separation of the binary has decreased by $\approx$0.3'' in 8 months.

preprint2014arXiv

High contrast imaging at the LBT: the LEECH exoplanet imaging survey

In Spring 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its $\sim$130-night campaign from the Large Binocular Telescope (LBT) atop Mt Graham, Arizona. This survey benefits from the many technological achievements of the LBT, including two 8.4-meter mirrors on a single fixed mount, dual adaptive secondary mirrors for high Strehl performance, and a cold beam combiner to dramatically reduce the telescope's overall background emissivity. LEECH neatly complements other high-contrast planet imaging efforts by observing stars at L' (3.8 $μ$m), as opposed to the shorter wavelength near-infrared bands (1-2.4 $μ$m) of other surveys. This portion of the spectrum offers deep mass sensitivity, especially around nearby adolescent ($\sim$0.1-1 Gyr) stars. LEECH's contrast is competitive with other extreme adaptive optics systems, while providing an alternative survey strategy. Additionally, LEECH is characterizing known exoplanetary systems with observations from 3-5$μ$m in preparation for JWST.

preprint2014arXiv

HST Rotational Spectral Mapping of Two L-Type Brown Dwarfs: Variability In and Out of Water Bands Indicates High-Altitude Haze Layers

We present time-resolved near-infrared spectroscopy of two L5 dwarfs, 2MASS J18212815+1414010 and 2MASS J15074759-1627386, observed with the Wide Field Camera 3 instrument on the Hubble Space Telescope (HST). We study the wavelength dependence of rotation-modulated flux variations between 1.1 $μ$m and 1.7 $μ$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $μ$m and 1.4 $μ$m vary at similar amplitudes as the adjacent continuum. This differs from the results of previous HST observations of L/T transition dwarfs, in which the water absorption at 1.4 $μ$m displays variations of about half of the amplitude at other wavelengths. We find that the relative amplitude of flux variability out of the water band with respect to that in the water band shows a increasing trend from the L5 dwarfs toward the early T dwarfs. We utilize the models of Saumon & Marley (2008) and find that the observed variability of the L5 dwarfs can be explained by the presence of spatially varying high-altitude haze layers above the condensate clouds. Therefore, our observations show that the heterogeneity of haze layers - the driver of the variability - must be located at very low pressures, where even the water opacity is negligible. In the near future, the rotational spectral mapping technique could be utilized for other atomic and molecular species to probe different pressure levels in the atmospheres of brown dwarfs and exoplanets and uncover both horizontal and vertical cloud structures.

preprint2014arXiv

LEECH: A 100 Night Exoplanet Imaging Survey at the LBT

In February 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its 100-night campaign from the Large Binocular Telescope atop Mount Graham in Arizona. LEECH nearly complements other high-contrast planet imaging efforts by observing stars in L' band (3.8 microns) as opposed to the shorter wavelength near-infrared bands (1-2.3 microns). This part of the spectrum offers deeper mass sensitivity for intermediate age (several hundred Myr-old) systems, since their Jovian-mass planets radiate predominantly in the mid-infrared. In this proceedings, we present the science goals for LEECH and a preliminary contrast curve from some early data.

preprint2013arXiv

Benchmark Tests for Markov Chain Monte Carlo Fitting of Exoplanet Eclipse Observations

Ground-based observations of exoplanet eclipses provide important clues to the planets' atmospheric physics, yet systematics in light curve analyses are not fully understood. It is unknown if measurements suggesting near-infrared flux densities brighter than models predict are real, or artifacts of the analysis processes. We created a large suite of model light curves, using both synthetic and real noise, and tested the common process of light curve modeling and parameter optimization with a Markov Chain Monte Carlo (MCMC) algorithm. With synthetic white-noise models, we find that input eclipse signals are generally recovered within 10% accuracy for eclipse depths greater than the noise amplitude, and to smaller depths for higher sampling rates and longer baselines. Red-noise models see greater discrepancies between input and measured eclipse signals, often biased in one direction. Finally, we find that in real data, systematic biases result even with a complex model to account for trends, and significant false eclipse signals may appear in a non-Gaussian distribution. To quantify the bias and validate an eclipse measurement, we compare both the planet-hosting star and several of its neighbors to a separately-chosen control sample of field stars. Re-examining the Rogers et al. (2009) Ks-band measurement of CoRoT-1b finds an eclipse $3190^{+370}_{-440}$ ppm deep centered at $ϕ_{me}$=$0.50418^{+0.00197}_{-0.00203}$. Finally, we provide and recommend the use of selected datasets we generated as a benchmark test for eclipse modeling and analysis routines, and propose criteria to verify eclipse detections.

preprint2013arXiv

Brown Dwarf Photospheres are Patchy: A Hubble Space Telescope Near-infrared Spectroscopic Survey Finds Frequent Low-level Variability

Condensate clouds strongly impact the spectra of brown dwarfs and exoplanets. Recent discoveries of variable L/T transition dwarfs argued for patchy clouds in at least some ultracool atmospheres. This study aims to measure the frequency and level of spectral variability in brown dwarfs and to search for correlations with spectral type. We used HST/WFC3 to obtain spectroscopic time series for 22 brown dwarfs of spectral types ranging from L5 to T6 at 1.1-1.7 $μ$m for $\approx$40 min per object. Using Bayesian analysis, we find 6 brown dwarfs with confident $(p>95\%)$ variability in the relative flux in at least one wavelength region at sub-percent precision, and 5 brown dwarfs with tentative $(p>68\%)$ variability. We derive a minimum variability fraction $f_{min}=27^{+11}_{-7}\%$ over all covered spectral types. The fraction of variables is equal within errors for mid L, late L and mid T spectral types; for early T dwarfs we do not find any confident variable but the sample is too small to derive meaningful limits. For some objects, the variability occurs primarily in the flux peak in the J or H band, others are variable throughout the spectrum or only in specific absorption regions. Four sources may have broad-band peak-to-peak amplitudes exceeding 1%. Our measurements are not sensitive to very long periods, inclinations near pole-on and rotationally symmetric heterogeneity. The detection statistics are consistent with most brown dwarf photospheres being patchy. While multiple-percent near-infrared variability may be rare and confined to the L/T transition, low-level heterogeneities are a frequent characteristic of brown dwarf atmospheres.

preprint2013arXiv

HST Spectral Mapping of L/T Transition Brown Dwarfs Reveals Cloud Thickness Variations

Most directly imaged giant exoplanets are fainter than brown dwarfs with similar spectra. To explain their relative underluminosity unusually cloudy atmospheres have been proposed. However, with multiple parameters varying between any two objects, it remained difficult to observationally test this idea. We present a new method, sensitive time-resolved Hubble Space Telescope near-infrared spectroscopy, to study two rotating L/T transition brown dwarfs (2M2139 and SIMP0136). The observations provide spatially and spectrally resolved mapping of the cloud decks of the brown dwarfs. The data allow the study of cloud structure variations while other parameters are unchanged. We find that both brown dwarfs display variations of identical nature: J- and H-band brightness variations with minimal color and spectral changes. Our light curve models show that even the simplest surface brightness distributions require at least three elliptical spots. We show that for each source the spectral changes can be reproduced with a linear combination of only two different spectra, i.e. the entire surface is covered by two distinct types of regions. Modeling the color changes and spectral variations together reveal patchy cloud covers consisting of a spatially heterogenous mix of low-brightness, low-temperature thick clouds and brighter, thin and warm clouds. We show that the same thick cloud patches seen in our varying brown dwarf targets, if extended to the entire photosphere, predict near-infrared colors/magnitudes matching the range occupied by the directly imaged exoplanets that are cooler and less luminous than brown dwarfs with similar spectral types. This supports the models in which thick clouds are responsible for the near infrared properties of these underluminous exoplanets.

preprint2013arXiv

Protoplanetary Disk Masses from Stars to Brown Dwarfs

We present SCUBA-2 850um observations for 7 very low mass stars (VLMS) and brown dwarfs (BDs): 3 in Taurus, 4 in the TWA, and all classical T Tauri (cTT) analogs. We detect 2 of the 3 Taurus disks, but none of the TWA ones. Our 3sigma limits correspond to a dust mass of 1.2 MEarth in Taurus and a mere 0.2 MEarth in the TWA (3--10x deeper than previous work). We combine our data with other sub-mm/mm surveys of Taurus, rho Oph and the TWA to investigate trends in disk mass and grain growth during the cTT phase. We find : (1) The minimum disk outer radius required to explain the upper envelope of sub-mm/mm fluxes is 100 AU for intermediate-mass stars, solar-types and VLMS, and 20 AU for BDs. (2) While the upper envelope of disk masses increases with Mstar from BDs to VLMS to solar-types, no increase is seen from solar-type to intermediate-mass stars. We propose this is due to enhanced photoevaporation around intermediate masses. (3) Many disks around Taurus and rho Oph intermediate-mass and solar-type stars evince an opacity index beta of 0--1, indicating large grains. Of the only four VLMS/BDs in these regions with multi-wavelength data, three are consistent with large grains, though optically thick disks are not ruled out. (4) For the TWA VLMS (TWA 30A, B), combining our fluxes with Mdot and ages suggests substantial grain growth by 10 Myr. The degree of grain growth in the TWA BDs (2M1207A, SSPM1102) remains largely unknown. (5) A Bayesian analysis shows that mean(log[Mdisk/Mstar]) = -2.4, roughly constant all the way from intermediate-mass stars to VLMS/BDs, and (6) the disk mass in close solar-type Taurus binaries is significantly lower than in singles (by a factor of 10), while that in wide solar-type Taurus binaries is closer to that in singles (lower by a factor of 3). (7) We discuss the implications for planet formation, and for the dependence of Mdot on Mstar.

preprint2013arXiv

Weather on Other Worlds I: Detection of Periodic Variability in the L3 Dwarf DENIS-P J1058.7-1548 with Precise Multi-Wavelength Photometry

Photometric monitoring from warm Spitzer reveals that the L3 dwarf DENIS-P J1058.7-1548 varies sinusoidally in brightness with a period of 4.25 +0.26/-0.16 hours and an amplitude of 0.388 +/- 0.043% (peak-to-valley) in the 3.6 micron band, confirming the reality of a 4.31 +/- 0.31 hour periodicity detected in J-band photometry from the SOAR telescope. The J-band variations are a factor of 2.17 +/- 0.35 larger in amplitude than those at 3.6 microns, while 4.5 micron Spitzer observations yield a 4.5 micron/3.6 micron amplitude ratio of only 0.23 +/- 0.15, consistent with zero 4.5 micron variability. This wide range in amplitudes indicates rotationally modulated variability due to magnetic phenomena and/or inhomogeneous cloud cover. Weak H-alpha emission indicates some magnetic activity, but it is difficult to explain the observed amplitudes by magnetic phenomena unless they are combined with cloud inhomogeneities (which might have a magnetic cause). However, inhomogenous cloudcover alone can explain all our observations, and our data align with theory in requiring that the regions with the thickest clouds also have the lowest effective temperature. Combined with published v sin(i) results, our rotation period yields a 95% confidence lower limit of 0.111 solar radii on the radius of DENIS-P J1058.7-1548, suggesting upper limits of 320 Myr and 0.055 solar masses on the age and mass. These limits should be regarded cautiously because of 3-sigma inconsistencies with other data; however, a lower limit of 45 degrees on the inclination is more secure. DENIS-P J1058.7-1548 is only the first of nearly two dozen low-amplitude variables discovered and analyzed by the Weather on Other Worlds project.

preprint2012arXiv

First Light LBT AO Images of HR 8799 bcde at 1.65 and 3.3 Microns: New Discrepancies between Young Planets and Old Brown Dwarfs

As the only directly imaged multiple planet system, HR 8799 provides a unique opportunity to study the physical properties of several planets in parallel. In this paper, we image all four of the HR 8799 planets at H-band and 3.3 microns with the new LBT adaptive optics system, PISCES, and LBTI/LMIRCam. Our images offer an unprecedented view of the system, allowing us to obtain H and 3.3$ micron photometry of the innermost planet (for the first time) and put strong upper-limits on the presence of a hypothetical fifth companion. We find that all four planets are unexpectedly bright at 3.3 microns compared to the equilibrium chemistry models used for field brown dwarfs, which predict that planets should be faint at 3.3 microns due to CH4 opacity. We attempt to model the planets with thick-cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to fit the 3.3 micron photometry increases the predicted L' (3.8 microns) flux enough that it is inconsistent with observations. In an effort to fit the SED of the HR 8799 planets, we construct mixtures of cloudy atmospheres, which are intended to represent planets covered by clouds of varying opacity. In this scenario, regions with low opacity look hot and bright, while regions with high opacity look faint, similar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. Our mixed cloud models reproduce all of the available data, but self-consistent models are still necessary to demonstrate their viability.

preprint2012arXiv

The Gray Needle: Large Grains in the HD 15115 Debris Disk from LBT/PISCES/Ks and LBTI/LMIRcam/L' Adaptive Optics Imaging

We present diffraction-limited \ks band and \lprime adaptive optics images of the edge-on debris disk around the nearby F2 star HD 15115, obtained with a single 8.4 m primary mirror at the Large Binocular Telescope. At \ks band the disk is detected at signal-to-noise per resolution element (SNRE) \about 3-8 from \about 1-2\fasec 5 (45-113 AU) on the western side, and from \about 1.2-2\fasec 1 (63-90 AU) on the east. At \lprime the disk is detected at SNRE \about 2.5 from \about 1-1\fasec 45 (45-90 AU) on both sides, implying more symmetric disk structure at 3.8 \microns . At both wavelengths the disk has a bow-like shape and is offset from the star to the north by a few AU. A surface brightness asymmetry exists between the two sides of the disk at \ks band, but not at \lprime . The surface brightness at \ks band declines inside 1\asec (\about 45 AU), which may be indicative of a gap in the disk near 1\asec. The \ks - \lprime disk color, after removal of the stellar color, is mostly grey for both sides of the disk. This suggests that scattered light is coming from large dust grains, with 3-10 \microns -sized grains on the east side and 1-10 \microns dust grains on the west. This may suggest that the west side is composed of smaller dust grains than the east side, which would support the interpretation that the disk is being dynamically affected by interactions with the local interstellar medium.

preprint2012arXiv

Vertical Atmospheric Structure in a Variable Brown Dwarf: Pressure-dependent Phase Shifts in Simultaneous Hubble Space Telescope-Spitzer Light Curves

Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1-1.7 microns and broadband 4.5 microns. Spectroscopic observations were taken with Wide Field Camera 3 on board the Hubble Space Telescope and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability with a period of 1.4 hr at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ from wavelength to wavelength with a maximum difference of more than half of a rotational period. We compare the spectra with atmospheric models of different cloud prescriptions, from which we determine the pressure levels probed at different wavelengths. We find that the phase lag increases with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that could cause this trend of light curve phase with probed pressure level. These observations are the first to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus are an ideal test case for realistic three dimensional simulations of the atmospheric structure with clouds in brown dwarfs and extrasolar planets.

preprint2011arXiv

A Combined Subaru/VLT/MMT 1--5 Micron Study of Planets Orbiting HR 8799: Implications for Atmospheric Properties, Masses, and Formation

We present new 1--1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3--5 micron MMT/Clio data first presented by Hinz et al. (2010). Our VLT/NaCo data yields a detection of a fourth planet at a projected separation of ~ 15 AU -- "HR 8799e". We also report new, albeit weak detections of HR 8799b at 1.03 microns and 3.3 microns. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR8799c, and possibly HR 8799d, have near-to-mid IR colors/magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets' SEDs far more accurately and without the need for rescaling the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4--4.5, Teff = 900K, while HR 8799c, d, and (by inference) e have log(g) = 4--4.5, Teff = 1000--1200K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6--7 Mj, 7--10 Mj, 7--10 Mj and 7--10 Mj. 'Patchy" cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.

preprint2011arXiv

Measuring Organic Molecular Emission in Disks with Low Resolution Spitzer Spectroscopy

We explore the extent to which Spitzer IRS spectra taken at low spectral resolution can be used in quantitative studies of organic molecular emission from disks surrounding low mass young stars. We use Spitzer IRS spectra taken in both the high and low resolution modules for the same sources to investigate whether it is possible to define line indices that can measure trends in the strength of the molecular features in low resolution data. We find that trends in HCN emission strength seen in the high resolution data can be recovered in low resolution data. In examining the factors that influence the HCN emission strength, we find that the low resolution HCN flux is modestly correlated with stellar accretion rate and X-ray luminosity. Correlations of this kind are perhaps expected based on recent observational and theoretical studies of inner disk atmospheres. Our results demonstrate the potential of using the large number of low resolution disk spectra that reside in the Spitzer archive to study the factors that influence the strength of molecular emission from disks. Such studies would complement results for the much smaller number of circumstellar disks that have been observed at high resolution with IRS.

preprint2011arXiv

Resolving the inner regions of the HD97048 circumstellar disk with VLT/NACO polarimetric differential imaging

Circumstellar disks are the cradles of planetary systems and their physical and chemical properties directly influence the planet formation process. As most planets supposedly form in the inner disk regions, i.e., within a few tens of AU, it is crucial to study circumstellar disk on these scales to constrain the conditions for planet formation. Our aims are to characterize the inner regions of the circumstellar disk around the young Herbig Ae/Be star HD97048 in polarized light. We use VLT/NACO to observe HD97048 in polarimetric differential imaging (PDI) mode in the H and Ks band. We spatially resolve the disk around HD97048 in polarized flux in both filters on scales between ~0.1"-1.0" corresponding to the inner ~16-160 AU. Fitting isophots to the flux calibrated H-band image between 13 - 14 mag/arcsec^2 and 14 - 15 mag/arcsec^2 we derive a apparent disk inclination angle of 34+-5 deg and 47+-2 deg, respectively. The disk position angle in both brightness regimes is almost identical and roughly 80 deg. Along the disk major axis the surface brightness of the polarized flux drops from ~11 mag/arcsec^2 at ~0.1" (~16 AU) to ~15.3 mag/arcsec^2 at ~1.0" (~160 AU). The brightness profiles along the major axis are fitted with power-laws falling off as ~r^(-1.78+-0.02) in H and ~r^(-2.34+-0.04) in Ks. As the surface brightness drops off more rapidly in Ks compared to H the disks becomes relatively bluer at larger separations possibly indicating changing dust grain properties as a function of radius. For the first time the inner ~0.1"-1.0" (~16-160 AU) of the surface layer of the HD97048 circumstellar disk have been imaged in scattered light demonstrating the power of ground-based imaging polarimetry. Our data fill an important gap in a large collection of existing data including resolved thermal dust and PAH emission images as well as resolved gas emission lines.

preprint2010arXiv

Day-side z'-band emission and eccentricity of Wasp-12b

We report the detection of the eclipse of the very-hot Jupiter WASP-12b via z'-band time-series photometry obtained with the 3.5-meter ARC telescope at Apache Point Observatory. We measure a decrease in flux of 0.082+/-0.015% during the passage of the planet behind the star. That planetary flux is equally well reproduced by atmospheric models with and without extra absorbers, and blackbody models with f > 0.585+/-0.080. It is therefore necessary to measure the planet at other wavelengths to further constrain its atmospheric properties. The eclipse appears centered at phase = 0.5100 (+0.0072,-0.0061), consistent with an orbital eccentricity of |e cos w| = 0.016 (+0.011,-0.009) (see note at end of Section 4). If the orbit of the planet is indeed eccentric, the large radius of WASP-12b can be explained by tidal heating.

preprint2010arXiv

First Results From VLT NACO Apodizing Phase Plate: 4-micron Images of the Exoplanet beta Pictoris b

Direct imaging of exoplanets requires both high contrast and high spatial resolution. Here, we present the first scientific results obtained with the newly commissioned Apodizing Phase Plate coronagraph (APP) on VLT/NACO. We detected the exoplanet beta Pictoris b in the narrow band filter centered at 4.05 micron (NB4.05). The position angle (209.13 +- 2.12 deg) and the projected separation to its host star (0."354 +- 0".012, i.e., 6.8 +- 0.2 AU at a distance of 19.3 pc) are in good agreement with the recently presented data from Lagrange et al. (2010). Comparing the observed NB4.05 magnitude of 11.20 +- 0.23 mag to theoretical atmospheric models we find a best fit with a 7-10 M_Jupiter object for an age of 12 Myr, again in agreement with previous estimates. Combining our results with published L' photometry we can compare the planet's [L' - NB4.05] color to that of cool field dwarfs of higher surface gravity suggesting an effective temperature of ~1700 K. The best fit theoretical model predicts an effective temperature of ~1470 K, but this difference is not significant given our photometric uncertainties. Our results demonstrate the potential of NACO/APP for future planet searches and provides independent confirmation as well as complementary data for beta Pictoris b.

preprint2009arXiv

Constraining the orbit of the possible companion to Beta Pictoris: New deep imaging observations

We recently reported on the detection of a possible planetary-mass companion to Beta Pictoris at a projected separation of 8 AU from the star, using data taken in November 2003 with NaCo, the adaptive-optics system installed on the Very Large Telescope UT4. Eventhough no second epoch detection was available, there are strong arguments to favor a gravitationally bound companion rather than a background object. If confirmed and located at a physical separation of 8 AU, this young, hot (~1500 K), massive Jovian companion (~8 Mjup) would be the closest planet to its star ever imaged, could be formed via core-accretion, and could explain the main morphological and dynamical properties of the dust disk. Our goal was to return to Beta Pic five years later to obtain a second-epoch observation of the companion or, in case of a non-detection, constrain its orbit. Deep adaptive-optics L'-band direct images of Beta Pic and Ks-band Four-Quadrant-Phase-Mask (4QPM) coronagraphic images were recorded with NaCo in January and February 2009. We also use 4QPM data taken in November 2004. No point-like signal with the brightness of the companion candidate (apparent magnitudes L'=11.2 or Ks ~ 12.5) is detected at projected distances down to 6.5 AU from the star in the 2009 data. As expected, the non-detection does not allow to rule out a background object; however, we show that it is consistent with the orbital motion of a bound companion that got closer to the star since first observed in 2003 and that is just emerging from behind the star at the present epoch. We place strong constraints on the possible orbits of the companion and discuss future observing prospects.

preprint2009arXiv

Do we really know the dust? Systematics and uncertainties of the mid-infrared spectral analysis methods

The spectral region around 10 micrometer, showing prominent dust emission bands, is commonly used to derive the chemical composition of protoplanetary dust. Different analysis methods have been proposed for this purpose, but so far, no comparative test has been performed to test their validity. We calculated model spectra of disk models with different geometries and central sources, using a 2D radiative transfer code. These spectra were then fitted in a blind test using four different spectral decomposition methods. We studied the effect of disk structure (flared vs. flat), inclination angle, size of the inner disk hole and stellar luminosity on the fitted chemical composition. Our results show that the derived dust compositions by all methods deviate systematically from the real chemical composition. Out of the four tested spectral decomposition methods, our new two-layer temperature distribution method, differs the least from the input dust composition and the results show the weakest systematic effects. The reason for the deviations of the results given by other methods lies in their simplifying assumptions (e.g. single average grain temperature or one component continuum). We also tested the influence of different noise levels on the results of the spectral decomposition methods. We find that, for ground-based observations (8-13 micrometer) the expected uncertainty in the value of the crystallinity is about 11% for a signal-to-noise ratio of 100, while for space-based observations (7-17 micrometer) the uncertainty is about 5%. On the basis of our results, we propose a recipe for the analysis and interpretation of dust spectroscopy data in the mid-infrared which should be especially valuable for analysing Spitzer data and ground-based infrared spectroscopy data in the 10 micrometer window.

preprint2009arXiv

Ks-band detection of thermal emission and color constraints to CoRoT-1b: A low-albedo planet with inefficient atmospheric energy redistribution and a temperature inversion

We report the detection in Ks-band of the secondary eclipse of the hot Jupiter CoRoT-1b, from time series photometry with the ARC 3.5-m telescope at Apache Point Observatory. The eclipse shows a depth of 0.336+/-0.042 percent and is centered at phase 0.5022 (+0.0023,-0.0027), consistent with a zero eccentricity orbit ecosω = 0.0035 (+0.0036,-0.0042). We perform the first optical to near-infrared multi-band photometric analysis of an exoplanet's atmosphere and constrain the reflected and thermal emissions by combining our result with the recent 0.6, 0.71, and 2.09 micron secondary eclipse detections by Snellen et al. (2009), Gillon et al. (2009), and Alonso et al. (2009a). Comparing the multi-wavelength detections to state-of-the-art radiative-convective chemical-equilibrium atmosphere models, we find the near-infrared fluxes difficult to reproduce. The closest blackbody-based and physical models provide the following atmosphere parameters: a temperature T = 2454 (+84,-170) K, a very low Bond albedo A_B = 0.000 (+0.087,-0.000), and an energy redistribution parameter P_n = 0.1, indicating a small but nonzero amount of heat transfer from the day- to night-side. The best physical model suggests a thermal inversion layer with an extra optical absorber of opacity kappa_e =0.05cm^2g^-1, placed near the 0.1-bar atmospheric pressure level. This inversion layer is located ten times deeper in the atmosphere than the absorbers used in models to fit mid-infrared Spitzer detections of other irradiated hot Jupiters.