sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Apai, Daniel

Normalized to: Apai, D.

151 article(s) in total. 708 co-authors, from 1 to 44 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2055347
Identifying Exo-Earth Candidates in Direct Imaging Data through Bayesian Classification
Comments: 23 pages, 10 figures; published in AJ; re-submitted on 26 Feb 2020 to correct a typo in the metadata (no changes to the manuscript)
Submitted: 2019-10-29, last modified: 2020-02-26
Future space telescopes may be able to directly image $\sim$10 - 100 planets with sizes and orbits consistent with habitable surface conditions ("exo-Earth candidates" or EECs), but observers will face difficulty in distinguishing these from the potentially hundreds of non-habitable "false positives" which will also be detected. To maximize the efficiency of follow-up observations, a prioritization scheme must be developed to determine which planets are most likely to be EECs. In this paper, we present a Bayesian method for estimating the likelihood that any directly imaged extrasolar planet is a true exo-Earth candidate by interpreting the planet's apparent magnitude and separation in light of existing exoplanet statistics. As a specific application of this general framework, we use published estimates of the discovery yield of future space-based direct imaging mission concepts to conduct "mock surveys" in which we compute the likelihood that each detected planet is an EEC. We find that it will be difficult to determine which planets are EECs with $>50\%$ confidence using single-band photometry immediately upon their detection. The best way to reduce this ambiguity would be to constrain the planet's orbit by revisiting the system multiple times or through a radial velocity precursor survey. Astrometric or radial velocity constraints on the planet's mass would offer a lesser benefit. Finally, we show that a Bayesian approach to prioritizing targets would improve the follow-up efficiency of a direct imaging survey versus a blind approach using the same data. For example, the prioritized approach could reduce the amount of integration time required for the spectral detection (or rejection) of water absorption in most EECs by a factor of two.
[2]  [pdf] - 2072713
EDEN: Sensitivity Analysis and Transiting Planet Detection Limits for Nearby Late Red Dwarfs
Comments: Accepted to AJ
Submitted: 2020-02-23
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.
[3]  [pdf] - 2057751
Helios-r2 -- A new Bayesian, open-source retrieval model for brown dwarfs and exoplanet atmospheres
Comments: Accepted for publication in ApJ. The code is available under
Submitted: 2019-10-02, last modified: 2020-02-10
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.
[4]  [pdf] - 2061684
Cloud Atlas: High-precision HST/WFC3/IR Time-Resolved Observations of Directly-Imaged Exoplanet HD106906b
Comments: Accepted for publication to AJ
Submitted: 2020-01-22
HD106906b is an ~11$M_{\mathrm{Jup}}$, ~15Myr old directly-imaged exoplanet orbiting at an extremely large distance from its host star. The wide separation (7.11 arcsec) between HD106906b and its host star greatly reduces the difficulty in direct-imaging observations, making it one of the most favorable directly-imaged exoplanets for detailed characterization. In this paper, we present HST/WFC3/IR time-resolved observations of HD106906b in the F127M, F139M, and F153M bands. We have achieved ~1% precision in the lightcurves in all three bands. The F127M lightcurve demonstrates marginally-detectable ($2.7\sigma$ significance) variability with a best-fitting period of 4 hr, while the lightcurves in the other two bands are consistent with flat lines. We construct primary-subtracted deep images and use these images to exclude additional companions to HD106906 that are more massive than 4$M_{\mathrm{Jup}}$ and locate at projected distances of more than ~500 au. We measure the astrometry of HD106906b in two HST/WFC3 epochs and achieve precisions better than 2.5 mas. The position angle and separation measurements do not deviate from those in the 2004 HST/ACS/HRC images for more than $1\sigma$ uncertainty. We provide the HST/WFC3 astrometric results for 25 background stars that can be used as reference sources in future precision astrometry studies. Our observations also provide the first 1.4-micron water band photometric measurement for HD106906b. HD106906b's spectral energy distribution and the best-fitting BT-Settl model have an inconsistency in the 1.4-micron water absorption band, which highlights the challenges in modeling atmospheres of young planetary-mass objects.
[5]  [pdf] - 2025515
The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler
Comments: Accepted in ApJ
Submitted: 2019-05-21, last modified: 2019-12-13
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
[6]  [pdf] - 2057810
Cloud Atlas: Weak color modulations due to rotation in the planetary-mass companion GU Psc b and 11 other brown dwarfs
Comments: 15 pages, 7 figures, accepted for publication in AJ
Submitted: 2019-12-05
Among the greatest challenges in understanding ultra-cool brown dwarf and exoplanet atmospheres is the evolution of cloud structure as a function of temperature and gravity. In this study, we present the rotational modulations of GU Psc b -- a rare mid-T spectral type planetary-mass companion at the end of the L/T spectral type transition. Based on the HST/WFC3 1.1-1.67$\rm\, \mu m$ time-series spectra, we observe a quasi-sinusoidal light curve with a peak-to-trough flux variation of 2.7 % and a minimum period of eight hours. The rotation-modulated spectral variations are weakly wavelength-dependent, or largely gray between 1.1-1.67$\rm\,\mu$m. The gray modulations indicate that heterogeneous clouds are present in the photosphere of this low-gravity mid-T dwarf. We place the color and brightness variations of GU Psc b in the context of rotational modulations reported for mid-L to late-T dwarfs. Based on these observations, we report a tentative trend: mid-to-late T dwarfs become slightly redder in $J-H$ color with increasing $J$-band brightness, while L dwarfs become slightly bluer with increasing brightness. If this trend is verified with more T-dwarf samples, it suggests that in addition to the mostly gray modulations, there is a second-order spectral-type dependence on the nature of rotational modulations.
[7]  [pdf] - 2026283
ACCESS: A Visual to Near-infrared Spectrum of the Hot Jupiter WASP-43b with Evidence of $\rm H_2O$, but no evidence of Na or K
Comments: 27 pages, 18 figures, 7 tables. Accepted for publication in AJ. Updated affiliations
Submitted: 2019-11-08, last modified: 2019-11-11
We present a new ground-based visual transmission spectrum of the hot Jupiter WASP-43b, obtained as part of the ACCESS Survey. The spectrum was derived from four transits observed between 2015 and 2018, with combined wavelength coverage between 5,300 \r{A}-9,000 \r{A} and an average photometric precision of 708 ppm in 230 \r{A} bins. We perform an atmospheric retrieval of our transmission spectrum combined with literature HST/WFC3 observations to search for the presence of clouds/hazes as well as Na, K, H$\alpha$, and $\rm H_2O$ planetary absorption and stellar spot contamination over a combined spectral range of 5,318 \r{A}-16,420 \r{A}. We do not detect a statistically significant presence of Na I or K I alkali lines, or H$\alpha$ in the atmosphere of WASP-43b. We find that the observed transmission spectrum can be best explained by a combination of heterogeneities on the photosphere of the host star and a clear planetary atmosphere with $\rm H_2O$. This model yields a log-evidence of $8.26\pm0.42$ higher than a flat (featureless) spectrum. In particular, the observations marginally favor the presence of large, low-contrast spots over the four ACCESS transit epochs with an average covering fraction $f_\text{het} = 0.27^{+0.42}_{-0.16}$ and temperature contrast $\Delta T = 132\text{ K} \pm 132\text{ K}$. Within the planet's atmosphere, we recover a log $\rm H_2O$ volume mixing ratio of $-2.78^{+1.38}_{-1.47}$, which is consistent with previous $\rm H_2O$ abundance determinations for this planet.
[8]  [pdf] - 1975411
Cloud Atlas: Variability in and out of the Water Band in the Planetary-mass HD 203030B Points to Cloud Sedimentation in Low-gravity L Dwarfs
Comments: Accepted for publication in The Astrophysical Journal
Submitted: 2019-08-25
We use the Wide Field Camera 3 on the {\sl Hubble Space Telescope} to spectrophotometrically monitor the young L7.5 companion HD~203030B. Our time series reveal photometric variability at 1.27\,$\mu$m and 1.39\,$\mu$m on time scales compatible with rotation. We find a rotation period of $7.5^{+0.6}_{-0.5}$ h: comparable to those observed in other brown dwarfs and planetary-mass companions younger than 300 Myr. We measure variability amplitudes of $1.1\pm0.3\%$ (1.27\,$\mu$m) and $1.7\pm0.4\%$ (1.39\,$\mu$m), and a phase lag of 56$^\circ\pm$28$^\circ$ between the two light curves. We attribute the difference in photometric amplitudes and phases to a patchy cloud layer that is sinking below the level where water vapor becomes opaque. HD 203030B and the few other known variable young late-L dwarfs are unlike warmer (earlier-type and/or older) L dwarfs, for which variability is much less wavelength-dependent across the 1.1--1.7$\mu$m region. We further suggest that a sinking of the top-most cloud deck below the level where water or carbon monoxide gas become opaque may also explain the often enhanced variability amplitudes of even earlier-type low-gravity L dwarfs. Because these condensate and gas opacity levels are already well-differentiated in T dwarfs, we do not expect the same variability amplitude enhancement in young vs.\ old T dwarfs.
[9]  [pdf] - 1966860
3.8um Imaging of 400-600K Brown Dwarfs and Orbital Constraints for WISEP J045853.90+643452.6AB
Comments: Accepted for publication in The Astrophysical Journal, July 17 2019. This revision includes changes to the Appendix only. Additional new W1 photometry is given and Figure 11 is updated; Figure 13 contained erroneous data and is corrected; Figure 14 is new and shows the new relationships derived for transforming from ground-based L or M magnitudes to Spitzer [3.6] and [4.5]
Submitted: 2019-07-17, last modified: 2019-08-21
Half of the energy emitted by late-T- and Y-type brown dwarfs emerges at 3.5 < lambda um < 5.5. We present new L' (3.43 < lambda um < 4.11) photometry obtained at the Gemini North telescope for nine late-T and Y dwarfs, and synthesize L' from spectra for an additional two dwarfs. The targets include two binary systems which were imaged at a resolution of 0.25". One of these, WISEP J045853.90+643452.6AB, shows significant motion, and we present an astrometric analysis of the binary using Hubble Space Telescope, Keck Adaptive Optics, and Gemini images. We compare lambda ~4um observations to models, and find that the model fluxes are too low for brown dwarfs cooler than ~700K. The discrepancy increases with decreasing temperature, and is a factor of ~2 at T_eff=500K and ~4 at T_eff=400K. Warming the upper layers of a model atmosphere generates a spectrum closer to what is observed. The thermal structure of cool brown dwarf atmospheres above the radiative-convective boundary may not be adequately modelled using pure radiative equilibrium; instead heat may be introduced by thermochemical instabilities (previously suggested for the L- to T-type transition) or by breaking gravity waves (previously suggested for the solar system giant planets). One-dimensional models may not capture these atmospheres, which likely have both horizontal and vertical pressure/temperature variations.
[10]  [pdf] - 1953656
Thermal Infrared Imaging of MWC 758 with the Large Binocular Telescope: Planetary Driven Spiral Arms?
Comments: 10 pages, 4 figures, accepted for publication in ApJ
Submitted: 2019-07-15, last modified: 2019-08-12
Theoretical studies suggest that a giant planet around the young star MWC 758 could be responsible for driving the spiral features in its circumstellar disk. Here, we present a deep imaging campaign with the Large Binocular Telescope with the primary goal of imaging the predicted planet. We present images of the disk in two epochs in the $L^{\prime}$ filter (3.8 $\mu m$) and a third epoch in the $M^{\prime}$ filter (4.8 $\mu m$). The two prominent spiral arms are detected in each observation, which constitute the first images of the disk at $M^\prime$, and the deepest yet in $L^\prime$ ($\Delta L^\prime=$12.1 exterior to the disk at 5$\sigma$ significance). We report the detection of a S/N$\sim$3.9 source near the end of the Sourthern arm, and, from the source's detection at a consistent position and brightness during multiple epochs, we establish a $\sim$90% confidence-level that the source is of astrophysical origin. We discuss the possibilities that this feature may be a) an unresolved disk feature, and b) a giant planet responsible for the spiral arms, with several arguments pointing in favor of the latter scenario. We present additional detection limits on companions exterior to the spiral arms, which suggest that a $\lesssim$4 M$_{Jup}$ planet exterior to the spiral arms could have escaped detection. Finally, we do not detect the companion candidate interior to the spiral arms reported recently by Reggiani et al. (2018), although forward modelling suggests that such a source would have likely been detected.
[11]  [pdf] - 1923523
On the Mass Function, Multiplicity, and Origins of Wide-Orbit Giant Planets
Comments: 19 pages, 10 figures, accepted to ApJ
Submitted: 2019-04-12, last modified: 2019-07-26
A major outstanding question regarding the formation of planetary systems is whether wide-orbit giant planets form differently than close-in giant planets. We aim to establish constraints on two key parameters that are relevant for understanding the formation of wide-orbit planets: 1) the relative mass function and 2) the fraction of systems hosting multiple companions. In this study, we focus on systems with directly imaged substellar companions, and the detection limits on lower-mass bodies within these systems. First, we uniformly derive the mass probability distributions of known companions. We then combine the information contained within the detections and detection limits into a survival analysis statistical framework to estimate the underlying mass function of the parent distribution. Finally, we calculate the probability that each system may host multiple substellar companions. We find that 1) the companion mass distribution is rising steeply toward smaller masses, with a functional form of $N\propto M^{-1.3\pm0.3}$, and consequently, 2) many of these systems likely host additional undetected sub-stellar companions. Combined, these results strongly support the notion that wide-orbit giant planets are formed predominantly via core accretion, similar to the better studied close-in giant planets. Finally, given the steep rise in the relative mass function with decreasing mass, these results suggest that future deep observations should unveil a greater number of directly imaged planets.
[12]  [pdf] - 1921030
The HST Large Programme on NGC 6752. III. Detection of the Peak of the White Dwarf Luminosity Function
Comments: 10 pages, 1 table, 5 figures (3 in low resolution). Accepted for publication in MNRAS on 2019, July 12. Associated RGB-cube-WCS-fits file at
Submitted: 2019-07-14
We report on the white dwarf cooling sequence of the old globular cluster NGC 6752, which is chemically complex and hosts a blue horizontal branch. This is one of the last globular cluster white dwarf (WD) cooling sequences accessible to imaging by the Hubble Space Telescope. Our photometry and completeness tests show that we have reached the peak of the luminosity function of the WD cooling sequence, at a magnitude m_F606W=29.4+/-0.1, which is consistent with a formal age of ~14Gyr. This age is also consistent with the age from fits to the main-sequence turnoff (13-14Gyr), reinforcing our conclusion that we observe the expected accumulation of white dwarfs along the cooling sequence.
[13]  [pdf] - 1966846
Combining high contrast imaging and radial velocities to constrain the planetary architecture of nearby stars
Comments: 17 pages, 12 figures, accepted for publication in A&A
Submitted: 2019-07-09
Nearby stars are prime targets for exoplanet searches and characterization using a variety of detection techniques. Combining constraints from the complementary detection methods of high contrast imaging (HCI) and radial velocity (RV) can further constrain the planetary architectures of these systems because these methods place limits at different regions of the companion mass and semi-major axis parameter space. We aim to constrain the planetary architectures from the combination of HCI and RV data for 6 nearby stars within 6 pc: $\tau$ Ceti, Kapteyn's star, AX Mic, 40 Eri, HD 36395, and HD 42581. We compiled the sample from stars with available archival VLT/NACO HCI data at L$^{\prime}$ band (3.8 $\mu$m). The NACO data were fully reanalyzed using the state-of-the-art direct imaging pipeline PynPoint and combined with RV data from HARPS, Keck/HIRES, and CORALIE. A Monte Carlo approach was used to assess the completeness in the companion mass/semi-major axis parameter space from the combination of the HCI and RV data sets. We find that the HCI data add significant information to the RV constraints, increasing the completeness for certain companions masses/semi-major axes by up to 68 - 99% for 4 of the 6 stars in our sample, and by up to 1 - 13% for the remaining stars. The improvements are strongest for intermediate semi-major axes (15 - 40 AU), corresponding to the semi-major axes of the ice giants in our own solar system. The HCI mass limits reach 5 - 20 $M_{\textrm{Jup}}$ in the background-limited regime, depending on the age of the star. Through the combination of HCI and RV data, we find that stringent constraints can be placed on the possible substellar companions in these systems. Applying these methods systematically to nearby stars will quantify our current knowledge of the planet population in the solar neighborhood and inform future observations.
[14]  [pdf] - 1929831
A Thousand Earths: A Very Large Aperture, Ultralight Space Telescope Array for Atmospheric Biosignature Surveys
Comments: Accepted in the Astronomical Journal. 28 pages, 12 figures. More info on the project website:
Submitted: 2019-06-12
An outstanding, multi-disciplinary goal of modern science is the study of the diversity of potentially Earth-like planets and the search for life in them. This goal requires a bold new generation of space telescopes, but even the most ambitious designs yet hope to characterize several dozen potentially habitable planets. Such a sample may be too small to truly understand the complexity of exo-earths. We describe here a notional concept for a novel space observatory designed to characterize 1,000 transiting exo-earth candidates. The Nautilus concept is based on an array of inflatable spacecraft carrying very large diameter (8.5m), very low-weight, multi-order diffractive optical elements (MODE lenses) as light-collecting elements. The mirrors typical to current space telescopes are replaced by MODE lenses with a 10 times lighter areal density that are 100 times less sensitive to misalignments, enabling light-weight structure. MODE lenses can be cost-effectively replicated through molding. The Nautilus mission concept has a potential to greatly reduce fabrication and launch costs, and mission risks compared to the current space telescope paradigm through replicated components and identical, light-weight unit telescopes. Nautilus is designed to survey transiting exo-earths for biosignatures up to a distance of 300 pc, enabling a rigorous statistical exploration of the frequency and properties of life-bearing planets and the diversity of exo-earths.
[15]  [pdf] - 1871553
Cloud Atlas: Rotational Spectral Modulations and potential Sulfide Clouds in the Planetary-mass, Late T-type Companion Ross 458C
Comments: Accepted in ApJL
Submitted: 2019-03-26
Measurements of photometric variability at different wavelengths provide insights into the vertical cloud structure of brown dwarfs and planetary-mass objects. In seven Hubble Space Telescope consecutive orbits, spanning $\sim$10 h of observing time}, we obtained time-resolved spectroscopy of the planetary-mass T8-dwarf Ross 458C using the near-infrared Wide Field Camera 3. We found spectrophotometric variability with a peak-to-peak signal of 2.62$\pm$0.02 % (in the 1.10-1.60~$\mu$m white light curve). Using three different methods, we estimated a rotational period of 6.75$\pm$1.58~h for the white light curve, and similar periods for narrow $J$- and $H$- band light curves. Sine wave fits to the narrow $J$- and $H$-band light curves suggest a tentative phase shift between the light curves with wavelength when we allow different periods between both light curves. If confirmed, this phase shift may be similar to the phase shift detected earlier for the T6.5 spectral type 2MASS J22282889-310262. We find that, in contrast with 2M2228, the variability of Ross~458C shows evidence for a {color trend} within the narrow $J$-band, but gray variations in the narrow $H$-band. The spectral time-resolved variability of Ross 458C might be potentially due to heterogeneous sulfide clouds in the atmosphere of the object. Our discovery extends the study of spectral modulations of condensate clouds to the coolest T dwarfs, planetary-mass companions.
[16]  [pdf] - 1850053
Astro2020 White Paper: The L/T Transition
Submitted: 2019-03-15
The L/T transition is an important evolutionary phase in brown dwarf atmospheres, providing us with a unique opportunity to explore the effects of clouds, convection, winds, gravity and metallicity across a very narrow temperature range. Understanding these physical processes is critical for understanding ultracool atmospheres. In the next decade, we will answer three key questions regarding L/T transition atmospheres: 1. What is the physical mechanism behind the L/T transition? 2. What is the spatial extent of atmospheric structures at the L/T transition? 3. How do gravity and metallicity affect the L/T transition? The theory and methods developed for brown dwarfs will be used in the 2030s and beyond for solar-system age giant exoplanets and eventually habitable zone earth analogues. Developing these techniques now are crucial.
[17]  [pdf] - 1849323
Constraining Stellar Photospheres as an Essential Step for Transmission Spectroscopy of Small Exoplanets
Comments: Science white paper submitted in response to the the U.S. National Academies of Science, Engineering, and Medicine's call for community input to the Astro2020 Decadal Survey; 9 pages, 3 figures
Submitted: 2019-03-14
Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.
[18]  [pdf] - 1848642
A Statistical Comparative Planetology Approach to Maximize the Scientific Return of Future Exoplanet Characterization Efforts
Comments: White paper submitted in response to the solicitation of feedback for the "2020 Decadal Survey" by the National Academy of Sciences
Submitted: 2019-03-12
Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.
[19]  [pdf] - 1847056
Astro2020 Science White Paper: High-Resolution Spectroscopic Surveys of Ultracool Dwarf Stars & Brown Dwarfs
Comments: 9 pages, submitted to Astro2020 Science White Paper Call
Submitted: 2019-03-11
High resolution spectroscopy of the lowest-mass stars and brown dwarfs reveals their origins, multiplicity, compositions and physical properties, with implications for the star formation and chemical evolution history of the Milky Way. We motivate the need for high-resolution, infrared spectroscopic surveys to reach these faint sources.
[20]  [pdf] - 1847061
Discovery of Cold Brown Dwarfs or Free-Floating Giant Planets Close to the Sun
Comments: White Paper submitted for Astro2020 Science
Submitted: 2019-03-11
This White Paper describes the opportunities for discovery of Jupiter-mass objects with 300K atmospheres. The discovery and characterization of such cold objects is vital for understanding the low-mass terminus of the initial mass function and for optimizing the study of exoplanets by the next generation of large telescopes, space probes and space missions.
[21]  [pdf] - 1847026
Planet formation: The case for large efforts on the computational side
Comments: White paper submitted to the Astro2020 decadal survey
Submitted: 2019-03-11
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.
[22]  [pdf] - 1842515
Cloud Atlas: High-Contrast Time-Resolved Observations of Planetary-Mass Companions
Comments: Accepted for publication in AAS Journal
Submitted: 2019-01-31, last modified: 2019-02-19
Directly-imaged planetary-mass companions offer unique opportunities in atmospheric studies of exoplanets. They share characteristics of both brown dwarfs and transiting exoplanets, therefore, are critical for connecting atmospheric characterizations for these objects. Rotational phase mapping is a powerful technique to constrain the condensate cloud properties in ultra-cool atmospheres. Applying this technique to directly-imaged planetary-mass companions will be extremely valuable for constraining cloud models in low mass and surface gravity atmospheres and for determining the rotation rate and angular momentum of substellar companions. Here, we present Hubble Space Telescope Wide Field Camera 3 near-infrared time-resolved photometry for three planetary-mass companions, AB Pic B, 2M0122B, and 2M1207b. Using two-roll differential imaging and hybrid point spread function modeling, we achieve sub-percent photometric precision for all three observations. We find tentative modulations ($<\!\!2\sigma$) for AB Pic B and 2M0122B but cannot reach conclusive results on 2M1207b due to strong systematics. The relatively low significance of the modulation measurements cannot rule out the hypothesis that these planetary-mass companions have the same vertical cloud structures as brown dwarfs. Our rotation rate measurements, combined with archival period measurements of planetary-mass companions and brown dwarfs do not support a universal mass-rotation relation. The high precision of our observations and the high occurrence rates of variable low-surface gravity objects encourage high-contrast time-resolved observations with the James Webb Space Telescope.
[23]  [pdf] - 1825567
The HST Large Programme on NGC 6752. I. Serendipitous discovery of a dwarf Galaxy in background
Comments: 5 pages, 1 table, 3 figures (2 in low resolution). Published on MNRAS Letters and freely available Related press releases: and , (and only in italian: )
Submitted: 2019-02-01
As part of a large Hubble Space Telescope investigation aiming at reaching the faintest stars in the Galactic globular cluster NGC 6752, an ACS/WFC field was the subject of deep optical observations reaching magnitudes as faint as V~30. In this field we report the discovery of Bedin I a dwarf spheroidal galaxy too faint and too close to the core of NGC 6752 for detection in earlier surveys. As it is of broad interest to complete the census of galaxies in the local Universe, in this Letter we provide the position of this new object along with preliminary assessments of its main param eters. Assuming the same reddening as for NGC 6752, we estimate a distance modulus of (m-M)_0 = 29.70+/-0.13 from the observed red giant branch, i.e., 8.7 (+0.5 -0.7) Mpc, and size of ~840x340pc, about 1/5 the size of the LMC. A comparison of the observed colour-magnitude diagram with synthetic counterparts that account for the galaxy distance modulus, reddening and photometric errors, suggests the presence of an old (~13Gyr) and metal poor ([Fe/H]~-1.3) population. This object is most likely a relatively isolated satellite dwarf spheroidal galaxy of the nearby great spiral NGC 6744, or potentially the most distant isolated dwarf spheroidal known with a secure distance.
[24]  [pdf] - 1823111
The HST Large Programme on NGC6752. II. Multiple populations at the bottom of the main sequence probed in NIR
Comments: 7 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2019-01-22
Historically, multiple populations in Globular Clusters (GCs) have been mostly studied from ultraviolet and optical filters down to stars that are more massive than ~0.6 solar masses. Here we exploit deep near-infrared (NIR) photometry from the Hubble Space Telescope to investigate multiple populations among M-dwarfs in the GC NGC6752. We discovered that the three main populations (A, B and C), previously observed in the brightest part of the color-magnitude diagram, define three distinct sequences that run from the main-sequence (MS) knee towards the bottom of the MS (~0.15 solar masses). These results, together with similar findings on NGC2808, M4, and omega Centauri, demonstrate that multiple sequences of M-dwarfs are common features of the color-magnitude diagrams of GCs. The three sequences of low-mass stars in NGC6752 are consistent with stellar populations with different oxygen abundances. The range of [O/Fe] needed to reproduce the NIR CMD of NGC6752 is similar to the oxygen spread inferred from high-resolution spectroscopy of red-giant branch (RGB) stars. The relative numbers of stars in the three populations of M-dwarfs are similar to those derived among RGB and MS stars more massive than ~0.6 solar masses. As a consequence, the evidence that the properties of multiple populations do not depend on stellar mass is a constraint for the formation scenarios.
[25]  [pdf] - 1823042
ACCESS: Ground-based Optical Transmission Spectroscopy of the Hot Jupiter WASP-4b
Comments: 31 pages, 19 figures, 8 tables. Accepted for publication in AJ
Submitted: 2018-12-18
We present an optical transmission spectrum of the atmosphere of WASP-4b obtained through observations of four transits with Magellan/IMACS. Using a Bayesian approach to atmospheric retrieval, we find no evidence for scattering or absorption features in our transit spectrum. Our models include a component to model the transit light source effect (spectral contamination from unocculted spots on the stellar photosphere), which we show can have a marked impact on the observed transmission spectrum for reasonable spot covering fractions (< 5%); this is the first such analysis for WASP-4b. We are also able to fit for the size and temperature contrast of spots observed during the second and third transits, finding evidence for both small, cool and large, warm spot-like features on the photosphere. Finally, we compare our results to those published by Huitson et al. (2017) using Gemini/GMOS and May et al. (2018) using IMACS, and find that our data are in agreement.
[26]  [pdf] - 1830571
The Transit Light Source Effect II: The Impact of Stellar Heterogeneity on Transmission Spectra of Planets Orbiting Broadly Sun-like Stars
Comments: 30 pages, 13 figures, 6 tables, accepted for publication in AJ
Submitted: 2018-12-14
Transmission spectra probe exoplanetary atmospheres, but they can also be strongly affected by heterogeneities in host star photospheres through the transit light source effect. Here we build upon our recent study of the effects of unocculted spots and faculae on M-dwarf transmission spectra, extending the analysis to FGK dwarfs. Using a suite of rotating model photospheres, we explore spot and facula covering fractions for varying activity levels and the associated stellar contamination spectra. Relative to M dwarfs, we find that the typical variabilities of FGK dwarfs imply lower spot covering fractions, though they generally increase with later spectral types, from $\sim 0.1\%$ for F dwarfs to 2-4$\%$ for late-K dwarfs. While the stellar contamination spectra are considerably weaker than those for typical M dwarfs, we find that typically active G and K dwarfs produce visual slopes that are detectable in high-precision transmission spectra. We examine line offsets at H$\alpha$ and the Na and K doublets and find that unocculted faculae in K dwarfs can appreciably alter transit depths around the Na D doublet. We find that band-averaged transit depth offsets at molecular bands for CH$_{4}$, CO, CO$_{2}$, H$_{2}$O, N$_{2}$O, O$_{2}$, and O$_{3}$ are not detectable for typically active FGK dwarfs, though stellar TiO/VO features are potentially detectable for typically active late-K dwarfs. Generally, this analysis shows that inactive FGK dwarfs do not produce detectable stellar contamination features in transmission spectra, though active FGK host stars can produce such features and care is warranted in interpreting transmission spectra from these systems.
[27]  [pdf] - 1830552
Cloud Atlas: Hubble Space Telescope Near-Infrared Spectral Library of Brown Dwarfs, Planetary-mass companions, and hot Jupiters
Submitted: 2018-12-10
Bayesian atmospheric retrieval tools can place constraints on the properties of brown dwarfs and hot Jupiters atmospheres. To fully exploit these methods, high signal-to-noise spectral libraries with well-understood uncertainties are essential. We present a high signal-to-noise spectral library (1.10-1.69 microns) of the thermal emission of 76 brown dwarfs and hot Jupiters. All our spectra have been acquired with the Hubble Space Telescope's Wide Field Camera 3 instrument and its G141 grism. The near-infrared spectral types of these objects range from L4 to Y1. Eight of our targets have estimated masses below the deuterium-burning limit. We analyze the database to identify peculiar objects and/or multiple systems, concluding that this sample includes two very-low-surface-gravity objects and five intermediate-surface-gravity objects. In addition, spectral indices designed to search for composite atmosphere brown dwarfs, indicate that eight objects in our sample are strong candidates to have such atmospheres. None of these objects are overluminous, thus their composite atmospheres are unlikely a companion-induced artifact. Five of the eight confirmed candidates have been reported as photometrically variable, suggesting that composite atmospheric indices are useful in identifying brown dwarfs with strongly heterogeneous cloud covers. We compare hot Jupiters and brown dwarfs in a near-infrared color-magnitude diagram. We confirm that the coldest hot Jupiters in our sample have spectra similar to mid-L dwarfs, and the hottest hot Jupiters have spectra similar to those of M-dwarfs. Our sample provides a uniform dataset of a broad range of ultracool atmospheres, allowing large-scale, comparative studies, and providing a HST legacy spectral library.
[28]  [pdf] - 1794119
The LEECH Exoplanet Imaging Survey: Limits on Planet Occurrence Rates Under Conservative Assumptions
Comments: 31 pages, 13 figures, accepted to AJ
Submitted: 2018-10-24, last modified: 2018-12-06
We present the results of the largest $L^{\prime}$ ($3.8~\mu$m) direct imaging survey for exoplanets to date, the Large Binocular Telescope Interferometer (LBTI) Exozodi Exoplanet Common Hunt (LEECH). We observed 98 stars with spectral types from B to M. Cool planets emit a larger share of their flux in $L^{\prime}$ compared to shorter wavelengths, affording LEECH an advantage in detecting low-mass, old, and cold-start giant planets. We emphasize proximity over youth in our target selection, probing physical separations smaller than other direct imaging surveys. For FGK stars, LEECH outperforms many previous studies, placing tighter constraints on the hot-start planet occurrence frequency interior to $\sim20$ au. For less luminous, cold-start planets, LEECH provides the best constraints on giant-planet frequency interior to $\sim20$ au around FGK stars. Direct imaging survey results depend sensitively on both the choice of evolutionary model (e.g., hot- or cold-start) and assumptions (explicit or implicit) about the shape of the underlying planet distribution, in particular its radial extent. Artificially low limits on the planet occurrence frequency can be derived when the shape of the planet distribution is assumed to extend to very large separations, well beyond typical protoplanetary dust-disk radii ($\lesssim50$ au), and when hot-start models are used exclusively. We place a conservative upper limit on the planet occurrence frequency using cold-start models and planetary population distributions that do not extend beyond typical protoplanetary dust-disk radii. We find that $\lesssim90\%$ of FGK systems can host a 7 to 10 $M_{\mathrm{Jup}}$ planet from 5 to 50 au. This limit leaves open the possibility that planets in this range are common.
[29]  [pdf] - 1815083
VLT/SPHERE Multi-Wavelength High-Contrast Imaging of the HD 115600 Debris Disk: New Constraints on the Dust Geometry and the Presence of Young Giant Planets
Comments: 7 figures
Submitted: 2018-11-21
Young and dynamically active planetary systems can form disks of debris that are easier to image than the planets themselves. The morphology and evolution of these disks can help to infer the properties of the putative planets responsible for generating and shaping the debris structures. We present integral field spectroscopy and dual-band imaging from VLT/SPHERE ($1.0 - 1.7 \mu m$) of the debris disk around the young F2V/F3V star HD 115600. We aim to 1) characterize the geometry and composition of the debris ring, 2) search for thermal emission of young giant planets, and 3) in the absence of detected planets, to refine the inferred properties of plausible planets around HD 115600 to prepare future attempts to detect them. Using a different dust scattering model (ZODIPIC) than in the discovery paper (Henyey-Greenstein) to model the disk geometry, we find $a_0 = 46 \pm 2$ au for the disk's central radius and offsets $\Delta \alpha$, $\Delta \delta$ = $ -1.0 \pm 0.5, 0.5 \pm 0.5$ au. This offset is smaller than previously found, suggesting that unseen planets of lower masses could be sculpting the disk. Spectroscopy of the disk in Y-J bands with SPHERE shows reddish color, which becomes neutral or slightly blue in H band seen with GPI, broadly consistent with a mixed bulk disk composition of processed organics and water ice. While our observed field contains numerous background objects at wide separations, no exoplanet has been directly observed to a mass sensitivity limit of $2-3(5-7) M_{\rm J}$ between a projected separation of 40 and 200 au for hot (cold)-start models.
[30]  [pdf] - 1767545
ACCESS: A featureless optical transmission spectrum for WASP-19b from Magellan/IMACS
Comments: 20 pages (plus 5 for the Appendix), 17 figures, 5 tables. MNRAS, in press
Submitted: 2018-07-27, last modified: 2018-10-02
The short period ($0.94$-day) transiting exoplanet WASP-19b is an exceptional target for transmission spectroscopy studies, due to its relatively large atmospheric scale-height ($\sim 500$ km) and equilibrium temperature ($\sim 2100$ K). Here we report on six precise spectroscopic Magellan/IMACS observations, five of which target the full optical window from $0.45-0.9\mu$m and one targeting the $0.4-0.55\mu$m blue-optical range. Five of these datasets are consistent with a transmission spectrum without any significant spectral features, while one shows a significant slope as a function of wavelength, which we interpret as arising from photospheric heterogeneities in the star. Coupled with HST/WFC3 infrared observations, our optical/near-infrared measurements point to the presence of high altitude clouds in WASP-19b's atmosphere in agreement with previous studies. Using a semi-analytical retrieval approach, considering both planetary and stellar spectral features, we find a water abundance consistent with solar for WASP-19b and strong evidence for sub-solar abundances for optical absorbers such as TiO and Na; no strong optical slope is detected, which suggests that if hazes are present, they are much weaker than previously suggested. In addition, two spot-crossing events are observed in our datasets and analyzed, including one of the first unambiguously detected bright spot-crossing events on an exoplanet host star.
[31]  [pdf] - 1767366
The Near-Infrared Transmission Spectra of TRAPPIST-1 Planets b, c, d, e, f, and g and Stellar Contamination in Multi-Epoch Transit Spectra
Comments: 41 pages, 20 figures. Accepted for publication on AJ on 8/28/2018. Comments and suggestions are welcome
Submitted: 2018-02-06, last modified: 2018-08-30
The seven approximately Earth-sized transiting planets in the \object{TRAPPIST-1} system provide a unique opportunity to explore habitable zone and non-habitable zone small planets within the same system. Its habitable zone exoplanets -- due to their favorable transit depths -- are also worlds for which atmospheric transmission spectroscopy is within reach with the Hubble Space Telescope (HST) and with the James Webb Space Telescope (JWST). We present here an independent reduction and analysis of two \textit{HST} Wide Field Camera 3 (WFC3) near-infrared transit spectroscopy datasets for six planets (b through g). Utilizing our physically-motivated detector charge trap correction and a custom cosmic ray correction routine, we confirm the general shape of the transmission spectra presented by \textbf{\citet{deWit2016, deWit2018}}. Our data reduction approach leads to a 25\% increase in the usable data and reduces the risk of confusing astrophysical brightness variations (e.g., flares) with instrumental systematics. No prominent absorption features are detected in any individual planet's transmission spectra; by contrast, the combined spectrum of the planets shows a suggestive decrease around 1.4\,$\micron$ similar to an inverted water absorption feature. Including transit depths from \textit{K2}, the SPECULOOS-South Observatory, and \textit{Spitzer}, we find that the complete transmission spectrum is fully consistent with stellar contamination owing to the transit light source effect. These spectra demonstrate how stellar contamination can overwhelm planetary absorption features in low-resolution exoplanet transit spectra obtained by \textit{HST} and \textit{JWST} and also highlight the challenges in combining multi epoch observations for planets around rapidly rotating spotted stars.
[32]  [pdf] - 1743918
Retrieval of planetary and stellar properties in transmission spectroscopy with Aura
Comments: 22 pages, 8 figures, MNRAS, in press
Submitted: 2018-08-29
Transmission spectroscopy provides a powerful probe of the atmospheric properties of transiting exoplanets. To date, studies of exoplanets in transit have focused on inferring their atmospheric properties such as chemical compositions, cloud/haze properties, and temperature structures. However, surface inhomogeneities in the host stars of exoplanets in the form of cool spots and hot faculae can in principle imprint signatures on the observed planetary transit spectrum. Here we present Aura, a new retrieval paradigm for inferring both planetary and stellar properties from a transmission spectrum. We apply our retrieval framework to a sample of hot giant exoplanets to determine the significance of stellar heterogeneity and clouds/hazes in their spectra. The retrieval analyses distinguish four groups of planets. First, the spectra of WASP-6b and WASP-39b are best characterised by imprints of stellar heterogeneity and hazes and/or clouds. HD 209458b and HAT-P-12b comprise the second group for which there is weak evidence for stellar heterogeneity and a high significance of hazes and/or clouds. The third group constitutes HAT-P-1b and WASP-31b and shows weak evidence against stellar heterogeneity but weak to substantial indications of clouds/hazes. The fourth group -- WASP-19b, WASP-17b, and WASP-12b -- is fit best by molecular and alkali absorbers with H$_2$ scattering without evidence for stellar heterogeneity and weak to no evidence for clouds/hazes. Our retrieval methodology paves the way to simultaneous information on the star and planet from higher resolution spectra using future facilities such as the James Webb Space Telescope and large ground-based facilities.
[33]  [pdf] - 1732818
Magellan Adaptive Optics Imaging of PDS 70: Measuring the Mass Accretion Rate of a Young Giant Planet within a Gapped Disk
Comments: 6 pages, 5 figures, accepted to ApJL
Submitted: 2018-07-27
PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, M\"uller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$\alpha$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination of these filters allows for the accretion rate of the young planet to be inferred, as hot infalling hydrogen gas will emit strongly at H$\alpha$ over the optical continuum. We detected a source in H$\alpha$ at the position of PDS 70b on two sequential nights in May 2018, for which we establish a false positive probability of $<$0.1%. We conclude that PDS 70b is a young, actively accreting planet. We utilize the H$\alpha$ line luminosity to derive a mass accretion rate of $\dot M= 10^{-8\pm1}$ M$_{Jup}/yr$, where the large uncertainty is primarily due to the unknown amount of optical extinction from the circumstellar and circumplanetary disks. PDS 70b represents the second case of an accreting planet interior to a disk gap, and is among the early examples of a planet observed during its formation.
[34]  [pdf] - 1732722
An ALMA Survey of faint disks in the Chamaeleon I star-forming region: Why are some Class II disks so faint?
Comments: ApJ accepted
Submitted: 2018-06-12
ALMA surveys of nearby star-forming regions have shown that the dust mass in the disk is correlated with the stellar mass, but with a large scatter. This scatter could indicate either different evolutionary paths of disks or different initial conditions within a single cluster. We present ALMA Cycle 3 follow-up observations for 14 Class II disks that were low S/N detections or non-detections in our Cycle 2 survey of the $\sim 2$ Myr-old Chamaeleon I star-forming region. With 5 times better sensitivity, we detect millimeter dust continuum emission from six more sources and increase the detection rate to 94\% (51/54) for Chamaeleon I disks around stars earlier than M3. The stellar-disk mass scaling relation reported in \citet{pascucci2016} is confirmed with these updated measurements. Faint outliers in the $F_{mm}$--$M_*$ plane include three non-detections (CHXR71, CHXR30A, and T54) with dust mass upper limits of 0.2 M$_\oplus$ and three very faint disks (CHXR20, ISO91, and T51) with dust masses $\sim 0.5$ M$_\oplus$. By investigating the SED morphology, accretion property and stellar multiplicity, we suggest for the three millimeter non-detections that tidal interaction by a close companion ($<$100 AU) and internal photoevaporation may play a role in hastening the overall disk evolution. The presence of a disk around only the secondary star in a binary system may explain the observed stellar SEDs and low disk masses for some systems.
[35]  [pdf] - 1712979
The Exoplanet Population Observation Simulator. I - The Inner Edges of Planetary Systems
Comments: Accepted in AAS journals, code available on github
Submitted: 2018-05-21
The Kepler survey provides a statistical census of planetary systems out to the habitable zone. Because most planets are non-transiting, orbital architectures are best estimated using simulated observations of ensemble populations. Here, we introduce EPOS, the Exoplanet Population Observation Simulator, to estimate the prevalence and orbital architectures of multi-planet systems based on the latest Kepler data release, DR25. We estimate that at least 42% of sun-like stars have nearly coplanar planetary systems with 7 or more exoplanets. The fraction of stars with at least one planet within 1 au could be as high as 100% depending on assumptions about the distribution of single transiting planets. We estimate an occurrence rate of planets in the habitable zone around sun-like stars of eta_earth=36+-14%. The innermost planets in multi-planet systems are clustered around an orbital period of 10 days (0.1 au), reminiscent of the protoplanetary disk inner edge or could be explained by a planet trap at that location. Only a small fraction of planetary systems have the innermost planet at long orbital periods, with fewer than ~8% and ~3% having no planet interior to the orbit of Mercury and Venus, respectively. These results reinforce the view that the solar system is not a typical planetary system, but an outlier among the distribution of known exoplanetary systems. We predict that at least half of the habitable zone exoplanets are accompanied by (non-transiting) planets at shorter orbital periods, hence knowledge of a close-in exoplanet could be used as a way to optimize the search for Earth-size planets in the Habitable Zone with future direct imaging missions.
[36]  [pdf] - 1654591
Understanding Stellar Contamination in Exoplanet Transmission Spectra as an Essential Step in Small Planet Characterization
Comments: White paper submitted to the NAS Committee on Exoplanet Science Strategy
Submitted: 2018-03-23
Transmission spectroscopy during planetary transits is expected to be a major source of information on the atmospheres of small (approximately Earth-sized) exoplanets in the next two decades. This technique, however, is intrinsically affected by stellar spectral contamination caused by the fact that stellar photo- and chromospheres are not perfectly homogeneous. Such stellar contamination will often reach or exceed the signal introduced by the planetary spectral features. Finding effective methods to correct stellar contamination -- or at least to quantify its possible range -- for the most important exoplanets is a necessary step for our understanding of exoplanet atmospheres. This will require significantly deepening our understanding of stellar heterogeneity, which is currently limited by the available data.
[37]  [pdf] - 1654584
A comprehensive understanding of planet formation is required for assessing planetary habitability and for the search for life
Comments: White paper submitted to the NAS Committee on Exoplanet Science Strategy
Submitted: 2018-03-23
Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.
[38]  [pdf] - 1648852
Cloud Atlas: Rotational Modulations in the L/T Transition Brown Dwarf Companion HN Peg B
Comments: Accepted for publication in AAS Journal
Submitted: 2018-01-29, last modified: 2018-02-16
Time-resolved observations of brown dwarfs' rotational modulations provide powerful insights into the properties of condensate clouds in ultra-cool atmospheres. Multi-wavelength light curves reveal cloud vertical structures, condensate particle sizes, and cloud morphology, which directly constrain condensate cloud and atmospheric circulation models. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared G141 taken in six consecutive orbits observations of HN Peg B, an L/T transition brown dwarf companion to a G0V type star. The best-fit sine wave to the $1.1-1.7\mu$m broadband light curve has the amplitude of $1.206\pm0.025\%$ and period of $15.4\pm0.5$ hr. The modulation amplitude has no detectable wavelength dependence except in the 1.4 $\mu$m water absorption band, indicating that the characteristic condensate particle sizes are large ($>1\mu$m). We detect significantly ($4.4\sigma$) lower modulation amplitude in the 1.4$\mu$m water absorption band, and find that HN Peg B's spectral modulation resembles those of early T type brown dwarfs. We also describe a new empirical interpolation method to remove spectral contamination from the bright host star. This method may be applied in other high-contrast time-resolved observations with WFC3.
[39]  [pdf] - 1634196
The Transit Light Source Effect: False Spectral Features and Incorrect Densities for M-dwarf Transiting Planets
Comments: accepted for publication in ApJ
Submitted: 2017-11-15, last modified: 2018-01-24
Transmission spectra are differential measurements that utilize stellar illumination to probe transiting exoplanet atmospheres. Any spectral difference between the illuminating light source and the disk-integrated stellar spectrum due to starspots and faculae will be imprinted in the observed transmission spectrum. However, few constraints exist for the extent of photospheric heterogeneities in M dwarfs. Here, we model spot and faculae covering fractions consistent with observed photometric variabilities for M dwarfs and the associated 0.3-5.5 $\mu$m stellar contamination spectra. We find that large ranges of spot and faculae covering fractions are consistent with observations and corrections assuming a linear relation between variability amplitude and covering fractions generally underestimate the stellar contamination. Using realistic estimates for spot and faculae covering fractions, we find stellar contamination can be more than $10 \times$ larger than transit depth changes expected for atmospheric features in rocky exoplanets. We also find that stellar spectral contamination can lead to systematic errors in radius and therefore the derived density of small planets. In the case of the TRAPPIST-1 system, we show that TRAPPIST-1's rotational variability is consistent with spot covering fractions $f_{spot} = 8^{+18}_{-7}\%$ and faculae covering fractions $f_{fac} = 54^{+16}_{-46}\%$. The associated stellar contamination signals alter transit depths of the TRAPPIST-1 planets at wavelengths of interest for planetary atmospheric species by roughly 1-15 $\times$ the strength of planetary features, significantly complicating $JWST$ follow-up observations of this system. Similarly, we find stellar contamination can lead to underestimates of bulk densities of the TRAPPIST-1 planets of $\Delta(\rho) = -3^{+3}_{-8} \%$, thus leading to overestimates of their volatile contents.
[40]  [pdf] - 1623900
Life Beyond the Solar System: Observation and Modeling of Exoplanet Environments
Comments: 8 pages including title page and references, 0 figures. This is a white paper submitted to the National Academies call on the Astrobiology Science Strategy for the Search for Life in the Universe (
Submitted: 2018-01-23
The search for life on planets outside our solar system has largely been the province of the astrophysics community until recently. A major development since the NASA Astrobiology Strategy 2015 document (AS15) has been the integration of other NASA science disciplines (planetary science, heliophysics, Earth science) with ongoing exoplanet research in astrophysics. The NASA Nexus for Exoplanet System Science (NExSS) provides a forum for scientists to collaborate across disciplines to accelerate progress in the search for life elsewhere. Here we describe recent developments in these other disciplines, with a focus on exoplanet properties and environments, and the prospects for future progress that will be achieved by integrating emerging knowledge from astrophysics with insights from these fields.
[41]  [pdf] - 1641416
The Orbit of the Companion to HD 100453A: Binary-Driven Spiral Arms in a Protoplanetary Disk
Comments: 28 pages, 11 figures, Accepted to ApJ
Submitted: 2018-01-11
HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD 100453AB between 2015-2017 utilizing extreme adaptive optics systems on the Very Large Telescope and Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a=1.06"+/-0.09", e=0.17+/-0.07, and i=32.5+/- 6.5 degrees. We utilized publicly available ALMA CO data to constrain the inclination of the disk to i~28 degrees, which is relatively co-planar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamical and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-dirven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material.
[42]  [pdf] - 1630335
The HST large programme on $\omega$ Centauri - II. internal kinematics
Comments: 22 pages (2 in appendix), 8 figures, 6 tables, accepted for publication in ApJ
Submitted: 2018-01-04, last modified: 2018-01-08
In this second installment of the series, we look at the internal kinematics of the multiple stellar populations of the globular cluster $\omega$ Centauri in one of the parallel Hubble Space Telescope (HST) fields, located at about 3.5 half-light radii from the center of the cluster. Thanks to the over 15-year-long baseline and the exquisite astrometric precision of the HST cameras, well-measured stars in our proper-motion catalog have errors as low as $\sim 10\ \mu$as yr$^{-1}$, and the catalog itself extends to near the hydrogen-burning limit of the cluster. We show that second-generation (2G) stars are significantly more radially anisotropic than first-generation (1G) stars. The latter are instead consistent with an isotropic velocity distribution. In addition, 1G have excess systemic rotation in the plane of the sky with respect to 2G stars. We show that the six populations below the main-sequence (MS) knee identified in our first paper are associated to the five main population groups recently isolated on the upper MS in the core of cluster. Furthermore, we find both 1G and 2G stars in the field to be far from being in energy equipartition, with $\eta_{\rm 1G}=-0.007\pm0.026$ for the former, and $\eta_{\rm 2G}=0.074\pm0.029$ for the latter, where $\eta$ is defined so that the velocity dispersion $\sigma_\mu$ scales with stellar mass as $\sigma_\mu \propto m^{-\eta}$. The kinematical differences reported here can help constrain the formation mechanisms for the multiple stellar populations in $\omega$ Centauri and other globular clusters. We make our astro-photometric catalog publicly available.
[43]  [pdf] - 1637637
The HST large programme on $\omega$ Centauri -- III. Absolute proper motion
Comments: 17 pages, 8 figures, 5 tables, accepted for publication in ApJ
Submitted: 2018-01-04, last modified: 2018-01-08
In this paper we report a new estimate of the absolute proper motion (PM) of the globular cluster NGC 5139 ($\omega$ Cen) as part of the HST large program GO-14118+14662. We analyzed a field 17 arcmin South-West of the center of $\omega$ Cen and computed PMs with an epoch span of $\sim$15.1 years. We employed 45 background galaxies to link our relative PMs to an absolute reference-frame system. The absolute PM of the cluster in our field is: $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.341 \pm 0.028 , -6.557 \pm 0.043)$ mas yr$^{-1}$. Upon correction for the effects of viewing perspective and the known cluster rotation, this implies that for the cluster center of mass $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.238 \pm 0.028, -6.716 \pm 0.043)$ mas yr$^{-1}$. This measurement is direct and independent, has the highest random and systematic accuracy to date, and will provide an external verification for the upcoming Gaia Data Release 2. It also differs from most reported PMs for $\omega$ Cen in the literature by more than 5$\sigma$, but consistency checks compared to other recent catalogs yield excellent agreement. We computed the corresponding Galactocentric velocity, calculated the implied orbit of $\omega$ Cen in two different Galactic potentials, and compared these orbits to the orbits implied by one of the PM measurements available in the literature. We find a larger (by about 500 pc) perigalactic distance for $\omega$ Cen with our new PM measurement, suggesting a larger survival expectancy for the cluster in the Galaxy.
[44]  [pdf] - 1608531
Cloud Atlas: Discovery of Rotational Spectral Modulations in a Low-mass, L-type Brown Dwarf Companion to a Star
Submitted: 2017-10-23
Observations of rotational modulations of brown dwarfs and giant exoplanets allow the characterization of condensate cloud properties. As of now rotational spectral modulations have only been seen in three L-type brown dwarfs. We report here the discovery of rotational spectral modulations in LP261-75B, an L6-type intermediate surface gravity companion to an M4.5 star. As a part of the Cloud Atlas Treasury program we acquired time-resolved Wide Field Camera 3 grism spectroscopy (1.1--1.69~$\mu$m) of LP261-75B. We find gray spectral variations with the relative amplitude displaying only a weak wavelength dependence and no evidence for lower-amplitude modulations in the 1.4~$\mu$m water band than in the adjacent continuum. The likely rotational modulation period is 4.78$\pm$0.95 h, although the rotational phase is not well sampled. The minimum relative amplitude in the white light curve measured over the whole wavelength range is 2.41$\pm$0.14\%. We report an unusual light curve with seemingly three peaks approximately evenly distributed in rotational phase. The spectral modulations suggests that the upper atmosphere cloud properties in {LP261-75B} are similar to two other mid-L dwarfs of typical infrared colors, but differ from that of the extremely red L-dwarf WISE0047.
[45]  [pdf] - 1589591
The Science Case for an Extended Spitzer Mission
Comments: 75 pages. See page 3 for Table of Contents and page 4 for Executive Summary
Submitted: 2017-10-11
Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.
[46]  [pdf] - 1586917
Exploring Other Worlds: Science Questions for Future Direct Imaging Missions (EXOPAG SAG15 Report)
Comments: For the report with higher-resolution images and for updates and supporting documentations please visit
Submitted: 2017-08-09
The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.
[47]  [pdf] - 1586903
Enabling Narrow(est) IWA Coronagraphy with STIS BAR5 and BAR10 Occulters
Comments: 33 pages, 27 figures, under review
Submitted: 2017-07-26
The Space Telescope Imaging Spectrograph's (STIS) BAR5 coronagraphic occulter was designed to provide high-contrast, visible-light, imaging in close (>= 0.15") angular proximity to bright point-sources. This is the smallest inner working angle (IWA) possible with HST's suite of coronagraphically augmented instruments through its mission lifetime. The STIS BAR5 image plane occulter, however, was damaged (bent and deformed) pre-launch and had not been enabled for GO science use following the installation of the instrument in 1997, during HST servicing mission SM2. With the success of the HST GO 12923 program, discussed herein, we explored and verified the functionality and utility of the BAR5 occulter. Thus, despite its physical damage, with updates to the knowledge of the aperture mask metrology and target pointing requirements, a robust determination of achievable raw and PSF-subtracted stellocentric image contrasts and fidelity was conducted. We also investigated, and herein report on, the use of the BAR10 rounded corners as narrow-angle occulters and compare IWA vs. contrast performance for the BAR5, BAR10, and Wedge occulters. With that, we provide recommendations for the most efficacious BAR5 and BAR10 use on-orbit in support of GO science. With color-matched PSF-template subtracted coronagraphy, inclusive of a small (+/- 1/4 pixel) 3-point cross-bar dithering strategy we recommend, we find BAR5 can deliver effective ~ 0.2" IWA image contrast of ~ 4 x 10^-5 pixel^-1 to ~ 1 x 10^-8 pixel^-1 at 2". With the pointing updates (to the PDB SIAF.dat file and/or implemented through APT) that we identified, and with observing strategies we explored, we recommend the use of STIS BAR5 coronagraphy as a fully "supported" capability for unique GO science.
[48]  [pdf] - 1584555
ALMA Observations of the Young Substellar Binary System 2M1207
Comments: 10 pages, 6 figures, accepted for publication in AJ
Submitted: 2017-06-12
We present ALMA observations of the 2M1207 system, a young binary made of a brown dwarf with a planetary-mass companion at a projected separation of about 40 au. We detect emission from dust continuum at 0.89 mm and from the $J = 3 - 2$ rotational transition of CO from a very compact disk around the young brown dwarf. The small radius found for this brown dwarf disk may be due to truncation from the tidal interaction with the planetary-mass companion. Under the assumption of optically thin dust emission, we estimated a dust mass of 0.1 $M_{\oplus}$ for the 2M1207A disk, and a 3$\sigma$ upper limit of $\sim 1~M_{\rm{Moon}}$ for dust surrounding 2M1207b, which is the tightest upper limit obtained so far for the mass of dust particles surrounding a young planetary-mass companion. We discuss the impact of this and other non-detections of young planetary-mass companions for models of planet formation, which predict the presence of circum-planetary material surrounding these objects.
[49]  [pdf] - 1584514
An ALMA Survey of CO isotopologue emission from Protoplanetary Disks in Chamaeleon I
Comments: accepted for publication in ApJ
Submitted: 2017-06-11
The mass of a protoplanetary disk limits the formation and future growth of any planet. Masses of protoplanetary disks are usually calculated from measurements of the dust continuum emission by assuming an interstellar gas-to-dust ratio. To investigate the utility of CO as an alternate probe of disk mass, we use ALMA to survey $^{13}$CO and C$^{18}$O J = $3-2$ line emission from a sample of 93 protoplanetary disks around stars and brown dwarfs with masses from 0.03 -- 2 M$_{\odot}$ in the nearby Chamaeleon I star-forming region. We detect $^{13}$CO emission from 17 sources and C$^{18}$O from only one source. Gas masses for disks are then estimated by comparing the CO line luminosities to results from published disk models that include CO freeze-out and isotope-selective photodissociation. Under the assumption of a typical ISM CO-to-H$_2$ ratios of $10^{-4}$, the resulting gas masses are implausibly low, with an average gas mass of $\sim$ 0.05 M$_{Jup}$ as inferred from the average flux of stacked $^{13}$CO lines. The low gas masses and gas-to-dust ratios for Cha I disks are both consistent with similar results from disks in the Lupus star-forming region. The faint CO line emission may instead be explained if disks have much higher gas masses, but freeze-out of CO or complex C-bearing molecules is underestimated in disk models. The conversion of CO flux to CO gas mass also suffers from uncertainties in disk structures, which could affect gas temperatures. CO emission lines will only be a good tracer of the disk mass when models for C and CO depletion are confirmed to be accurate.
[50]  [pdf] - 1584184
Hubble Space Telescope astrometry of the closest brown dwarf binary system -- I. Overview and improved orbit
Comments: 19 pages, 9 figures, 9 tables. Accepted for publication in MNRAS on 2017 May 9
Submitted: 2017-06-01
Located at ~2pc, the L7.5+T0.5 dwarfs system WISE J104915.57-531906.1 (Luhman16AB) is the third closest system known to Earth, making it a key benchmark for detailed investigation of brown dwarf atmospheric properties, thermal evolution, multiplicity, and planet-hosting frequency. In the first study of this series -- based on a multi-cycle Hubble Space Telescope (HST) program -- we provide an overview of the project and present improved estimates of positions, proper motions, annual parallax, mass ratio, and the current best assessment of the orbital parameters of the A-B pair. Our HST observations encompass the apparent periastron of the binary at 220.5+/-0.2 mas at epoch 2016.402. Although our data seem to be inconsistent with recent ground-based astrometric measurements, we also exclude the presence of third bodies down to Neptune masses and periods longer than a year.
[51]  [pdf] - 1582785
Weather on Other Worlds. IV. H$\alpha$ emission and photometric variability are not correlated in L0$-$T8 dwarfs
Comments: 12 pages, 4 figures, accepted for publication in ApJ
Submitted: 2017-04-23
Recent photometric studies have revealed that surface spots that produce flux variations are present on virtually all L and T dwarfs. Their likely magnetic or dusty nature has been a much-debated problem, the resolution to which has been hindered by paucity of diagnostic multi-wavelength observations. To test for a correlation between magnetic activity and photometric variability, we searched for H$\alpha$ emission among eight L3$-$T2 ultra-cool dwarfs with extensive previous photometric monitoring, some of which are known to be variable at 3.6 $\mu$m or 4.5 $\mu$m. We detected H$\alpha$ only in the non-variable T2 dwarf 2MASS J12545393$-$0122474. The remaining seven objects do not show H$\alpha$ emission, even though six of them are known to vary photometrically. Combining our results with those for 86 other L and T dwarfs from the literature show that the detection rate of H$\alpha$ emission is very high (94$\%$) for spectral types between L0 and L3.5 and much smaller (20$\%$) for spectral types $\ge$L4, while the detection rate of photometric variability is approximately constant (30$\%-$55$\%$) from L0 to T8 dwarfs. We conclude that chromospheric activity, as evidenced by H$\alpha$ emission, and large-amplitude photometric variability are not correlated. Consequently, dust clouds are the dominant driver of the observed variability of ultra-cool dwarfs at spectral types at least as early as L0.
[52]  [pdf] - 1582302
The HST Large Program on Omega Centauri. I. Multiple stellar populations at the bottom of the main sequence probed in NIR-Optical
Comments: 13 pages, 9 figures, accepted for publication in MNRAS
Submitted: 2017-04-03
As part of a large investigation with Hubble Space Telescope to study the faintest stars within the globular cluster Omega Centauri, in this work we present early results on the multiplicity of its main sequence (MS) stars, based on deep optical and near-infrared observations. By using appropriate color-magnitude diagrams we have identified, for the first time, the two main stellar populations I, and II along the entire MS, from the turn-off towards the hydrogen-burning limit. We have compared the observations with suitable synthetic spectra of MS stars and conclude that the two MSs are consistent with stellar populations with different metallicity, helium, and light-element abundance. Specifically, MS-I corresponds to a metal-poor stellar population ([Fe/H]~-1.7) with Y~ 0.25 and [O/Fe]~0.30. The MS-II hosts helium-rich (Y~0.37-0.40) stars with metallicity ranging from [Fe/H]~-1.7 to -1.4. Below the MS knee (mF160W~19.5, our photometry reveals that each of the two main MSs hosts stellar subpopulations with different oxygen abundances, with very O-poor stars ([O/Fe]~-0.5) populating the MS-II. Such a complexity has never been observed in previous studies of M-dwarfs in globular clusters. A few months before the lunch of the James Webb Space Telescope, these results demonstrate the power of optical and near-infrared photometry in the study of multiple stellar populations in globular clusters.
[53]  [pdf] - 1574709
Parallaxes and Infrared Photometry of three Y0 dwarfs
Comments: 12 Pages, 10 figures, Accepted by MNRAS
Submitted: 2017-03-30
We have followed up the three Y0 dwarfs WISEPA J041022.71+150248.5, WISEPA J173835.53+273258.9 and WISEPC J205628.90+145953.3 using the UKIRT/WFCAM telescope/instruments. We find parallaxes that are more consistent and accurate than previously published values. We estimate absolute magnitudes in photometric pass-bands from $Y$ to $W3$ and find them to be consistent between the three Y0 dwarfs indicating the inherent cosmic absolute magnitude spread of these objects is small. We examine the MKO $J$ magnitudes over the four year time line and find small but significant monotonic variations. Finally we estimate physical parameters from a comparison of spectra and parallax to equilibrium and non-equilibrium models finding values consistent with solar metallicity, an effective temperature of 450-475\,K and log~g of 4.0-4.5.
[54]  [pdf] - 1557957
NGC 1980 is not a foreground population of Orion: Spectroscopic survey of young stars with low extinction in Orion A
Comments: 34 pages, 23 figures, and 6 tables, accepted by AJ
Submitted: 2017-03-20
We perform a spectroscopic survey of the foreground population in Orion A with MMT/Hectospec. We use these data, along with archival spectroscopic data and photometric data, to derive spectral types, extinction values, and masses for 691 stars. Using the Spitzer Space Telescope data, we characterize the disk properties of these sources. We identify 37 new transition disk (TD) objects, one globally depleted disk candidate, and 7 probable young debris disks. We discover an object with a mass less than 0.018-0.030 M$_{\odot}$, which harbors a flaring disk. Using the H$\alpha$ emission line, we characterize the accretion activity of the sources with disks, and confirm that fraction of accreting TDs is lower than that of optically thick disks (46$\pm$7% versus 73$\pm$9%, respectively). Using kinematic data from the Sloan Digital Sky Survey and APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC), we confirm that the foreground population shows similar kinematics to their local molecular clouds and other young stars in the same regions. Using the isochronal ages, we find that the foreground population has a median age around 1-2 Myr, which is similar to the one of other young stars in Orion A. Therefore, our results argue against the presence of a large and old foreground cluster in front of Orion A.
[55]  [pdf] - 1574658
A Physical Model-based Correction for Charge Traps in the Hubble Space Telescope's Wide Field Camera 3 Near-IR Detector and Applications to Transiting Exoplanets and Brown Dwarfs
Comments: 16 pages, 13 figures, accepted to Astronomical Journal
Submitted: 2017-03-03, last modified: 2017-03-16
The Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) near-IR channel is extensively used in time-resolved observations, especially for transiting exoplanet spectroscopy and brown dwarf and directly imaged exoplanet rotational phase mapping. The ramp effect is the dominant source of systematics in the WFC3 for time-resolved observations, which limits its photometric precision. Current mitigation strategies are based on empirical fits and require additional orbits "to help the telescope reach a thermal equilibrium". We show that the ramp effect profiles can be explained and corrected with high fidelity using charge trapping theories. We also present a model for this process that can be used to predict and to correct charge trap systematics. Our model is based on a very small number of parameters that are intrinsic to the detector. We find that these parameters are very stable between the different datasets, and we provide best-fit values. Our model is tested with more than 120 orbits ($\sim40$ visits) of WFC3 observations and is proved to be able to provide near photon noise limited corrections for observations made with both staring and scanning modes of transiting exoplanets as well as for starting-mode observations of brown dwarfs. After our model correction, the light curve of the first orbit in each visit has the same photometric precision as subsequent orbits, so data from the first orbit need no longer be discarded. Near IR arrays with the same physical characteristics (e.g., JWST/NIRCam) may also benefit from the extension of this model, if similar systematic profiles are observed.
[56]  [pdf] - 1542974
Probabilistic Constraints on the Mass and Composition of Proxima b
Comments: 6 pages, 3 figures, 1 table. Accepted for publication in ApJL
Submitted: 2017-02-08
Recent studies regarding the habitability, observability, and possible orbital evolution of the indirectly detected exoplanet Proxima b have mostly assumed a planet with $M \sim 1.3$ $M_\oplus$, a rocky composition, and an Earth-like atmosphere or none at all. In order to assess these assumptions, we use previous studies of the radii, masses, and compositions of super-Earth exoplanets to probabilistically constrain the mass and radius of Proxima b, assuming an isotropic inclination probability distribution. We find it is ~90% likely that the planet's density is consistent with a rocky composition; conversely, it is at least 10% likely that the planet has a significant amount of ice or an H/He envelope. If the planet does have a rocky composition, then we find expectation values and 95% confidence intervals of $\left<M\right>_\text{rocky} = 1.63_{-0.72}^{+1.66}$ $M_\oplus$ for its mass and $\left<R\right>_\text{rocky} = 1.07_{-0.31}^{+0.38}$ $R_\oplus$ for its radius.
[57]  [pdf] - 1533211
Atmospheric Compositions of Three Brown Dwarfs and Implications for their Formation Conditions
Comments: Submitted to ApJ on Oct 5, 2016. Comments welcome
Submitted: 2016-12-09
The formation mechanisms and chemical compositions of brown dwarfs likely span a diverse range. If they are formed predominantly in isolation by gravitational collapse like dwarf stars then their compositions might follow those of dwarf stars, spanning a relatively narrow range in metallicities (~ -0.4 to +0.4 dex) and predominantly oxygen-rich (C/O <~ 0.7). On the other hand, giant planets in the solar system are all super-solar in metallicity ([C/H] > 0.5~dex), which is thought to be a consequence of formation by core-accretion in a circumstellar disk. In this study we determine the atmospheric chemical compositions of three brown dwarfs and find them to be neither characteristic of dwarf stars nor giant planets. We derive high-precision atmospheric chemical abundances using high-S/N HST near-infrared spectra of three mid T Dwarfs which, together with two previously studied objects, display distinctly sub-solar metallicities and span C/O ratios of 0.4-1.0. These abundance patterns indicate either an old sub-stellar population and a multitude of formation environments or objects formed in circumstellar disks followed by ejection later. We show that for a simple synthetic model population of brown dwarfs with the same spectral-type as our sample the predicted age distribution is predominantly older than the Sun and, hence, more metal poor, consistent with our derived abundances. The diverse C/O ratios we find are indicative of either different chemical reservoirs in their formation environments or different formation pathways. Our results open the possibility that cool brown dwarfs may provide important probes of early galactic chemical evolution and inhomogeneity.
[58]  [pdf] - 1532938
ACCESS I: An Optical Transmission Spectrum of GJ 1214b Reveals a Heterogeneous Stellar Photosphere
Comments: Accepted for publication in ApJ
Submitted: 2016-12-01
GJ 1214b is the most studied sub-Neptune exoplanet to date. Recent measurements have shown its near-infrared transmission spectrum to be flat, pointing to a high-altitude opacity source in the exoplanet's atmosphere, either equilibrium condensate clouds or photochemical hazes. Many photometric observations have been reported in the optical by different groups, though simultaneous measurements spanning the entire optical regime are lacking. We present an optical transmission spectrum ($4,500-9,260$\AA) of GJ 1214b in 14 bins measured with Magellan/IMACS repeatedly over three transits. We measure a mean planet-to-star radius ratio of ${R_{p}/R_{s} = 0.1146\pm{2\times10^{-4}}}$ and mean uncertainty of $\sigma(R_{p}/R_{s})=8.7\times10^{-4}$ in the spectral bins. The optical transit depths are shallower on average than observed in the near-infrared. We present a model for jointly incorporating the effects of a composite photosphere and atmospheric transmission (CPAT) through the exoplanet's limb, and use it to examine the cases of absorber and temperature heterogeneities in the stellar photosphere. We find the optical and near-infrared measurements are best explained by the combination of (1) photochemical haze in the exoplanetary atmosphere with a mode particle size $r=0.1~\mu$m and haze-forming efficiency $f_{haze}=10 \%$ and (2) faculae in the unocculted stellar disk with a temperature contrast $\Delta T=354^{+46}_{-46}$ K, assuming 3.2% surface coverage. The CPAT model can be used to assess potential contributions of heterogeneous stellar photospheres to observations of exoplanet transmission spectra, which will be important for searches for spectral features in the optical.
[59]  [pdf] - 1532880
A candidate planetary-mass object with a photoevaporating disk in Orion
Comments: 6 pages, 4 figures. Accepted for publication in ApJL
Submitted: 2016-11-29
In this work, we report the discovery of a candidate planetary-mass object with a photoevaporating protoplanetary disk, Proplyd 133-353, which is near the massive star $\theta^{1}$ Ori C at the center of the Orion Nebula Cluster (ONC). The object was known to have extended emission pointing away from $\theta^{1}$ Ori C, indicating ongoing external photoevaporation. Our near-infrared spectroscopic data suggests that the central source of Proplyd 133-353 is substellar ($\sim$M9.5), might have a mass probably less than 13 Jupiter mass and an age younger than 0.5 Myr. Proplyd 133-353 shows a similar ratio of X-ray luminosity to stellar luminosity to other young stars in the ONC with a similar stellar luminosity, and has a similar proper motion to the mean one of confirmed ONC members. We propose that Proplyd 133-353 was formed in a very low-mass dusty cloud near $\theta^{1}$ Ori C as a second-generation of star formation, which can explain both its young age and the presence of its disk.
[60]  [pdf] - 1532346
Shadows and spirals in the protoplanetary disk HD 100453
Comments: A&A, accepted
Submitted: 2016-10-31, last modified: 2016-11-03
Understanding the diversity of planets requires to study the morphology and the physical conditions in the protoplanetary disks in which they form. We observed and spatially resolved the disk around the ~10 Myr old protoplanetary disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and near-infrared wavelengths, reaching an angular resolution of ~0.02", and an inner working angle of ~0.09". We detect polarized scattered light up to ~0.42" (~48 au) and detect a cavity, a rim with azimuthal brightness variations at an inclination of 38 degrees, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint spiral-like feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of ~119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by ~72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim, induces an azimuthal variation of the scale height that can contribute to the brightness variations. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. The origin of such a misalignment in HD 100453, and of the spirals, is unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit.
[61]  [pdf] - 1489973
Cloud Atlas: Discovery of Patchy Clouds and High-amplitude Rotational Modulations In a Young, Extremely Red L-type Brown Dwarf
Comments: 9 pages, 4 figures, Accepted for publication in ApJL on September 15th
Submitted: 2016-09-15, last modified: 2016-09-19
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 {\mu}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.
[62]  [pdf] - 1528217
A Super-Solar Metallicity For Stars With Hot Rocky Exoplanets
Comments: 10 pages, 7 figures, 1 table. Accepted for publication in AJ
Submitted: 2016-09-19
The host star metallicity provide a measure of the conditions in protoplanetary disks at the time of planet formation. Using a sample of over 20,000 Kepler stars with spectroscopic metallicities from the LAMOST survey, we explore how the exoplanet population depends on host star metallicity as a function of orbital period and planet size. We find that exoplanets with orbital periods less than 10 days are preferentially found around metal-rich stars ([Fe/H]~ 0.15 +- 0.05 dex). The occurrence rates of these hot exoplanets increases to ~30% for super-solar metallicity stars from ~10% for stars with a sub-solar metallicity. Cooler exoplanets, that resides at longer orbital periods and constitute the bulk of the exoplanet population with an occurrence rate of >~ 90%, have host-star metallicities consistent with solar. At short orbital periods, P<10 days, the difference in host star metallicity is largest for hot rocky planets (<1.7 R_Earth), where the metallicity difference is [Fe/H] =~ 0.25 +- 0.07 dex. The excess of hot rocky planets around metal-rich stars implies they either share a formation mechanism with hot Jupiters, or trace a planet trap at the protoplanetary disk inner edge which is metallicity-dependent. We do not find statistically significant evidence for a previously identified trend that small planets toward the habitable zone are preferentially found around low-metallicity stars. Refuting or confirming this trend requires a larger sample of spectroscopic metallicities.
[63]  [pdf] - 1510323
A Steeper than Linear Disk Mass-Stellar Mass Scaling Relation
Comments: ApJ, in press
Submitted: 2016-08-11
The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.
[64]  [pdf] - 1490786
Panchromatic Imaging of a Transitional Disk: The Disk of GM Aur in Optical and FUV Scattered Light
Comments: 40 pages, 11 figures, 4 tables, accepted to ApJ
Submitted: 2016-07-12
We have imaged GM Aur with HST, detected its disk in scattered light at 1400A and 1650A, and compared these with observations at 3300A, 5550A, 1.1 microns, and 1.6 microns. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300A shows no evidence of the 24AU radius cavity that has been previously observed in sub-mm observations. Comparison with dust grain opacity models indicates the surface of the entire disk is populated with sub-micron grains. We have compiled an SED from 0.1 microns to 1 mm, and used it to constrain a model of the star+disk system that includes the sub-mm cavity using the Monte Carlo Radiative Transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of sub-microns grains interior to the sub-mm cavity wall. We suggest one explanation for this which could be due to a planet of mass <9 Jupiter masses interior to 24 AU. A unique cylindrical structure is detected in the FUV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semi-minor axis, but does not resemble an accretion-driven jet. The structure is limb-brightened and extends 190 +/- 35 AU above the disk midplane. The inner radius of the limb-brightening is 40 +/- 10 AU, just beyond the sub-millimeter cavity wall.
[65]  [pdf] - 1457472
Direct Imaging Discovery of a Jovian Exoplanet Within a Triple Star System
Submitted: 2016-07-08
Direct imaging allows for the detection and characterization of exoplanets via their thermal emission. We report the discovery via imaging of a young Jovian planet in a triple star system and characterize its atmospheric properties through near-infrared spectroscopy. The semi-major axis of the planet is closer relative to that of its hierarchical triple star system than for any known exoplanet within a stellar binary or triple, making HD 131399 dynamically unlike any other known system. The location of HD 131399Ab on a wide orbit in a triple system demonstrates that massive planets may be found on long and possibly unstable orbits in multi-star systems. HD 131399Ab is one of the lowest mass (4+/-1 MJup) and coldest (850+/-50 K) exoplanets to have been directly imaged.
[66]  [pdf] - 1439844
$Extrasolar~Storms$: Pressure-dependent Changes In Light Curve Phase In Brown Dwarfs From Simultaneous $Hubble$ and $Spitzer$ Observations
Comments: 34 pages, 22 figures, accepted for publication in ApJ
Submitted: 2016-05-09
We present $Spitzer$/IRAC Ch1 and Ch2 monitoring of six brown dwarfs during 8 different epochs over the course of 20 months. For four brown dwarfs, we also obtained simulataneous $HST$/WFC3 G141 Grism spectra during two epochs and derived light curves in five narrow-band filters. Probing different pressure levels in the atmospheres, the multi-wavelength light curves of our six targets all exhibit variations, and the shape of the light curves evolves over the timescale of a rotation period, ranging from 1.4 h to 13 h. We compare the shapes of the light curves and estimate the phase shifts between the light curves observed at different wavelengths by comparing the phase of the primary Fourier components. We use state-of-the-art atmosphere models to determine the flux contribution of different pressure layers to the observed flux in each filter. We find that the light curves that probe higher pressures are similar and in phase, but are offset and often different from the light curves that probe lower pressures. The phase differences between the two groups of light curves suggest that the modulations seen at lower and higher pressures may be introduced by different cloud layers.
[67]  [pdf] - 1436374
Maps of Evolving Cloud Structures in Luhman 16AB from HST Time-Resolved Spectroscopy
Comments: 11 pages, 13 figures, accepted for publication in ApJ, videos of the Luhman 16A & B maps can be found in:
Submitted: 2016-05-05
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.
[68]  [pdf] - 1420072
The young stellar population of Lynds 1340. An infrared view
Comments: 73 pages, 33 figures, 15 tables. Accepted for publication in ApJS
Submitted: 2016-02-25
We present results of an infrared study of the molecular cloud Lynds 1340, forming three groups of low and intermediate-mass stars. Our goals are to identify and characterise the young stellar population of the cloud, study the relationships between the properties of the cloud and the emergent stellar groups, and integrate L1340 into the picture of the star-forming activity of our Galactic environment. We selected candidate young stellar objects from the Spitzer and WISE data bases using various published color criteria, and classified them based on the slope of the spectral energy distribution. We identified 170 Class II, 27 Flat SED, and Class 0/I sources. High angular resolution near-infrared observations of the RNO 7 cluster, embedded in L1340, revealed eight new young stars of near-infrared excess. The surface density distribution of young stellar objects shows three groups, associated with the three major molecular clumps of L1340, each consisting of less than 100 members, including both pre-main sequence stars and embedded protostars. New Herbig--Haro objects were identified in the Spitzer images. Our results demonstrate that L1340 is a prolific star-forming region of our Galactic environment in which several specific properties of the intermediate-mass mode of star formation can be studied in detail.
[69]  [pdf] - 1378914
High-Cadence, High-Contrast Imaging for Exoplanet Mapping: Observations of the HR 8799 Planets with VLT/SPHERE Satellite Spot-Corrected Relative Photometry
Comments: Astrophysical Journal, in press
Submitted: 2016-02-08, last modified: 2016-02-09
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%.
[70]  [pdf] - 1938164
Exoplanet Exploration Program Analysis Group (ExoPAG) Report to Paul Hertz Regarding Large Mission Concepts to Study for the 2020 Decadal Survey
Comments: 22 pages, 2 tables, no figures
Submitted: 2015-12-31
This is a joint summary of the reports from the three Astrophysics Program Analysis Groups (PAGs) in response to the "Planning for the 2020 Decadal Survey" charge given by the Astrophysics Division Director Paul Hertz. This joint executive summary contains points of consensus across all three PAGs. Additional findings specific to the individual PAGs are reported separately in the individual reports. The PAGs concur that all four large mission concepts identified in the white paper as candidates for maturation prior to the 2020 Decadal Survey should be studied in detail. These include the Far-IR Surveyor, the Habitable-Exoplanet Imaging Mission, the UV/Optical/IR Surveyor, and the X-ray Surveyor. This finding is predicated upon assumptions outlined in the white paper and subsequent charge, namely that 1) major development of future large flagship missions under consideration are to follow the implementation phases of JWST and WFIRST; 2) NASA will partner with the European Space Agency on its L3 Gravitational Wave Surveyor; 3) the Inflation Probe be classified as a probe-class mission to be developed according to the 2010 Decadal Survey report. If these key assumptions were to change, this PAG finding would need to be re-evaluated. The PAGs find that there is strong community support for the second phase of this activity - maturation of the four proposed mission concepts via Science and Technology Definition Teams (STDTs). The PAGs find that there is strong consensus that all of the STDTs contain broad and interdisciplinary representation of the science community. Finally, the PAGs find that there is community support for a line of Probe-class missions within the Astrophysics mission portfolio (condensed).
[71]  [pdf] - 1370286
Discovery of Rotational Modulations in the Planetary-Mass Companion 2M1207b: Intermediate Rotation Period and Heterogeneous Clouds in a Low Gravity Atmosphere
Comments: Accepted for publication in ApJ. 13 pages, 5 figures, 2 tables
Submitted: 2015-12-08
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.
[72]  [pdf] - 1359005
The LEECH Exoplanet Imaging Survey: Characterization of the Coldest Directly Imaged Exoplanet, GJ 504 b, and Evidence for Super-Stellar Metallicity
Comments: resubmitted to ApJ with referee's comments
Submitted: 2015-11-30
As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.
[73]  [pdf] - 1362624
First light of the VLT planet finder SPHERE. III. New spectrophotometry and astrometry of the HR8799 exoplanetary system
Comments: 13 pages, 12 figures, accepted for publication in A&A
Submitted: 2015-11-12
The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $\mu$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $\mu$m, 1.667 $\mu$m), $K1K2$ (2.110 $\mu$m, 2.251 $\mu$m), and broadband $J$ (1.245 $\mu$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $\mu$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Lo\`eve image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.
[74]  [pdf] - 1319667
An Increase in the Mass of Planetary Systems around Lower-Mass Stars
Comments: Accepted in ApJ. 11 pages, 9 figures, 2 tables (1 machine readable table included in source)
Submitted: 2015-10-08, last modified: 2015-10-20
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.
[75]  [pdf] - 1293696
Aeolus: A Markov--Chain Monte Carlo code for mapping ultracool atmospheres. An application on Jupiter and brown dwarf HST light curves
Comments: Accepted by ApJ; A high-quality version of this paper can be found at
Submitted: 2015-10-14
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.
[76]  [pdf] - 1300388
Discovery of a Two-Armed Spiral Structure in the Gapped Disk in HD 100453
Comments: 12 pages, 5 figures, accepted for publication in ApJL
Submitted: 2015-10-08
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.
[77]  [pdf] - 1300387
Discovery of an Edge-on Debris Disk with a Dust Ring and an Outer Disk Wing-tilt Asymmetry
Comments: 13 pages, 3 figures, accepted for publication in ApJL
Submitted: 2015-10-08
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.
[78]  [pdf] - 1300332
Cloud Structure of the Nearest Brown Dwarfs II: High-amplitude variability for Luhman 16 A and B in and out of the 0.99 micron FeH feature
Comments: 11 pages, 7 figures. Accepted for publication in the Astrophysical Journal
Submitted: 2015-09-21
The re-emergence of the 0.99 $\mu$m FeH feature in brown dwarfs of early- to mid-T spectral type has been suggested as evidence for cloud disruption where flux from deep, hot regions below the Fe cloud deck can emerge. The same mechanism could account for color changes at the L/T transition and photometric variability. We present the first observations of spectroscopic variability of brown dwarfs covering the 0.99 $\mu$m FeH feature. We observed the spatially resolved very nearby brown dwarf binary WISE J104915.57-531906.1 (Luhman 16AB), a late-L and early-T dwarf, with HST/WFC3 in the G102 grism at 0.8-1.15 $\mu$m. We find significant variability at all wavelengths for both brown dwarfs, with peak-to-valley amplitudes of 9.3% for Luhman 16B and 4.5% for Luhman 16A. This represents the first unambiguous detection of variability in Luhman 16A. We estimate a rotational period between 4.5 and 5.5 h, very similar to Luhman 16B. Variability in both components complicates the interpretation of spatially unresolved observations. The probability for finding large amplitude variability in any two brown dwarfs is less than 10%. Our finding may suggest that a common but yet unknown feature of the binary is important for the occurrence of variability. For both objects, the amplitude is nearly constant at all wavelengths except in the deep K I feature below 0.84 $\mu$m. No variations are seen across the 0.99 $\mu$m FeH feature. The observations lend strong further support to cloud height variations rather than holes in the silicate clouds, but cannot fully rule out holes in the iron clouds. We re-evaluate the diagnostic potential of the FeH feature as a tracer of cloud patchiness.
[79]  [pdf] - 1223763
The LEECH Exoplanet Imaging Survey. Further constraints on the planet architecture of the HR 8799 system
Comments: 10 pages, 3 figures, accepted for publication in A&A
Submitted: 2014-12-22, last modified: 2015-03-05
Context. Astrometric monitoring of directly-imaged exoplanets allows the study of their orbital parameters and system architectures. Because most directly-imaged planets have long orbital periods (>20 AU), accurate astrometry is challenging when based on data acquired on timescales of a few years and usually with different instruments. The LMIRCam camera on the LBT is being used for the LEECH survey to search for and characterize young and adolescent exoplanets in L' band, including their system architectures. Aims. We first aim to provide a good astrometric calibration of LMIRCam. Then, we derive new astrometry, test the predictions of the orbital model of 8:4:2:1 mean motion resonance proposed by Go\'zdziewski & Migaszewski, and perform new orbital fitting of the HR 8799 bcde planets. We also present deep limits on a putative fifth planet interior to the known planets. Methods. We use observations of HR 8799 and the Theta1 Ori C field obtained during the same run in October 2013. Results. We first characterize the distortion of LMIRCam. We determine a platescale and a true north orientation for the images of 10.707 +/- 0.012 mas/pix and -0.430 +/- 0.076 deg, respectively. The errors on the platescale and true north orientation translate into astrometric accuracies at a separation of 1 of 1.1 mas and 1.3 mas, respectively. The measurements for all planets are usually in agreement within 3 sigma with the ephemeris predicted by Go\'zdziewski & Migaszewski. The orbital fitting based on the new astrometric measurements favors an architecture for the planetary system based on 8:4:2:1 mean motion resonance. The detection limits allow us to exclude a fifth planet slightly brighter/more massive than HR 8799 b at the location of the 2:1 resonance with HR 8799 e (~9.5 AU) and about twice as bright as HR 8799 cde at the location of the 3:1 resonance with HR 8799 e (~7.5 AU).
[80]  [pdf] - 1231879
Volatile Delivery to Planets from Water-rich Planetesimals around Low Mass Stars
Comments: Accepted to ApJ
Submitted: 2015-02-25
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.
[81]  [pdf] - 1224007
The Inner Disk Structure, Disk-Planet Interactions, and Temporal Evolution in the Beta Pictoris System: A Two-Epoch HST/STIS Coronagraphic Study
Comments: Astrophysical Journal, in press; 70 pages; 21 figures
Submitted: 2015-01-13
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.
[82]  [pdf] - 919849
A stellar-mass-dependent drop in planet occurrence rates
Comments: 18 pages, 8 figures. Published in ApJ
Submitted: 2014-06-28, last modified: 2015-01-10
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.
[83]  [pdf] - 1223289
Stirring in massive, young debris discs from spatially resolved Herschel images
Comments: Accepted for publication in MNRAS, 22 pages, 7 figures, 6 tables (abstract abridged due to arXiv requirements)
Submitted: 2014-11-21
A significant fraction of main-sequence stars are encircled by dusty debris discs, where the short-lived dust particles are replenished through collisions between planetesimals. Most destructive collisions occur when the orbits of smaller bodies are dynamically stirred up, either by the gravitational effect of locally formed Pluto-sized planetesimals (self-stirring scenario), or via secular perturbation caused by an inner giant planet (planetary stirring). The relative importance of these scenarios in debris systems is unknown. Here we present new Herschel Space Observatory imagery of 11 discs selected from the most massive and extended known debris systems. All discs were found to be extended at far-infrared wavelengths, five of them being resolved for the first time. We evaluated the feasibility of the self-stirring scenario by comparing the measured disc sizes with the predictions of the model calculated for the ages of our targets. We concluded that the self-stirring explanation works for seven discs. However, in four cases, the predicted pace of outward propagation of the stirring front, assuming reasonable initial disc masses, was far too low to explain the radial extent of the cold dust. Therefore, for HD 9672, HD 16743, HD 21997, and HD 95086, another explanation is needed. We performed a similar analysis for {\ss} Pic and HR 8799, reaching the same conclusion. We argue that planetary stirring is a promising possibility to explain the disk properties in these systems. In HR 8799 and HD 95086 we may already know the potential perturber, since their known outer giant planets could be responsible for the stirring process. Our study demonstrates that among the largest and most massive debris discs self-stirring may not be the only active scenario, and potentially planetary stirring is responsible for destructive collisions and debris dust production in a number of systems.
[84]  [pdf] - 1451816
Weather on Other Worlds. II. Survey Results: Spots Are Ubiquitous on L and T Dwarfs
Comments: 42 pages, 11 figures, 2 tables, accepted by ApJ
Submitted: 2014-11-11
We present results from the "Weather on Other Worlds" Spitzer Exploration Science program to investigate photometric variability in L and T dwarfs, usually attributed to patchy clouds. We surveyed 44 L3-T8 dwarfs, spanning a range of $J-K_s$ colors and surface gravities. We find that 14/23 (61%; 95% confidence interval: 41%-78%) of our single L3-L9.5 dwarfs are variable with peak-to-peak amplitudes between 0.2% and 1.5%, and 5/16 (31%; 95% confidence interval: 14%-56%) of our single T0-T8 dwarfs are variable with amplitudes between 0.8% and 4.6%. After correcting for sensitivity, we find that 80% (95% confidence interval: 53%-100%) of L dwarfs vary by >0.2%, and 36% (95% confidence interval: 19%-52%) of T dwarfs vary by >0.4%. Given viewing geometry considerations, we conclude that photospheric heterogeneities causing >0.2% 3-5-micron flux variations are present on virtually all L dwarfs, and probably on most T dwarfs. A third of L dwarf variables show irregular light curves, indicating that L dwarfs may have multiple spots that evolve over a single rotation. Also, approximately a third of the periodicities are on time scales >10 h, suggesting that slowly-rotating brown dwarfs may be common. We observe an increase in the maximum amplitudes over the entire spectral type range, revealing a potential for greater temperature contrasts in T dwarfs than in L dwarfs. We find a tentative association (92% confidence) between low surface gravity and high-amplitude variability among L3-L5.5 dwarfs. Although we can not confirm whether lower gravity is also correlated with a higher incidence of variables, the result is promising for the characterization of directly imaged young extrasolar planets through variability.
[85]  [pdf] - 1223119
HST Rotational Spectral Mapping of Two L-Type Brown Dwarfs: Variability In and Out of Water Bands Indicates High-Altitude Haze Layers
Comments: 23 pages, 4 figures, accepted for publication in ApJ Letters
Submitted: 2014-11-11
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 $\mu$m and 1.7 $\mu$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $\mu$m and 1.4 $\mu$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 $\mu$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.
[86]  [pdf] - 1222922
Cloud structure of the nearest brown dwarfs: Spectroscopic variability of Luhman 16AB from the Hubble Space Telescope
Comments: 14 pages, 12 figures. Accepted for publication in ApJ
Submitted: 2014-10-31
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 $\mu$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 $\mu$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.
[87]  [pdf] - 886449
New Candidate Eruptive Young Stars in Lynds 1340
Comments: 14 pages, 3 figures; some typos are corrected and journal reference is added
Submitted: 2014-10-03, last modified: 2014-10-22
We report on the discovery of three candidate eruptive young stars, found during our comprehensive multi-wavelength study of the young stellar population of the dark cloud L1340. These stars are as follows. (1) IRAS 02224+7227 (2MASS 02270555+7241167, HH 487 S) exhibited FUor-like spectrum in our low-resolution optical spectra. The available photometric data restrict its luminosity to 23 L_sun < L_bol < 59 L_sun. (2) 2MASS 02263797+7304575, identified as a classical T Tauri star during our H alpha survey, exhibited an EXor type brightening in 2005 November, at the time of the SDSS observations of the region. (3) 2MASS 02325605+7246055, a low-mass embedded young star, associated with a fan-shaped infrared nebula, underwent an outburst between the DSS1 and DSS2 surveys, leading to the appearance of a faint optical nebula. Our [SII] and H alpha images, as well as the Spitzer IRAC 4.5 micron images revealed Herbig-Haro objects associated with this star. Our results suggest that amplitudes and time scales of outbursts do not necessarily correlate with the evolutionary stage of the stars.
[88]  [pdf] - 1215588
High contrast imaging at the LBT: the LEECH exoplanet imaging survey
Comments: 12 pages, 5 figures. Proceedings of the SPIE, 9148-20
Submitted: 2014-07-10
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 $\mu$m), as opposed to the shorter wavelength near-infrared bands (1-2.4 $\mu$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$\mu$m in preparation for JWST.
[89]  [pdf] - 830033
Cloud structure of brown dwarfs from spectroscopic variability observations
Comments: 4 pages, 1 figure, submitted to the proceedings of "Gaia and the unseen - the brown dwarf question", Torino, 24-26 March 2014, to be published in Memorie della Societa' Astronomica Italiana (SAIt), eds Ricky Smart, David Barrado, Jackie Faherty
Submitted: 2014-06-01
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.
[90]  [pdf] - 1208738
A Search for Companions to Brown Dwarfs in the Taurus and Chamaeleon Star Forming Regions
Comments: Astrophysical Journal, in press
Submitted: 2014-04-01
We present the results of a search for companions to young brown dwarfs in the Taurus and Chamaeleon I star forming regions (1/2-3 Myr). We have used WFPC2 on board HST to obtain F791W and F850LP images of 47 members of these regions that have spectral types of M6-L0 (0.01-0.1 Msun). An additional late-type member of Taurus, FU Tau (M7.25+M9.25), was also observed with adaptive optics at Keck Observatory. We have applied PSF subtraction to the primaries and have searched the resulting images for objects that have colors and magnitudes that are indicative of young low-mass objects. Through this process, we have identified promising candidate companions to 2MASS J04414489+2301513 (rho=0.105"/15 AU), 2MASS J04221332+1934392 (rho=0.05"/7 AU), and ISO 217 (rho=0.03"/5 AU). We reported the discovery of the first candidate in a previous study, showing that it has a similar proper motion as the primary through a comparison of astrometry measured with WFPC2 and Gemini adaptive optics. We have collected an additional epoch of data with Gemini that further supports that result. By combining our survey with previous high-resolution imaging in Taurus, Chamaeleon, and Upper Sco (10 Myr), we measure binary fractions of 14/93 = 0.15+0.05/-0.03 for M4-M6 (0.1-0.3 Msun) and 4/108 = 0.04+0.03/-0.01 for >M6 (<0.1 Msun) at separations of >10 AU. Given the youth and low density of these three regions, the lower binary fraction at later types is probably primordial rather than due to dynamical interactions among association members. The widest low-mass binaries (>100 AU) also appear to be more common in Taurus and Chamaeleon than in the field, which suggests that the widest low-mass binaries are disrupted by dynamical interactions at >10 Myr, or that field brown dwarfs have been born predominantly in denser clusters where wide systems are disrupted or inhibited from forming.
[91]  [pdf] - 1202590
Discovery of H-alpha Emission from the Close Companion Inside the Gap of Transitional Disk HD142527
Comments: 18 pages, 5 figures, accepted by The Astrophysical Journal Letters
Submitted: 2014-01-07
We utilized the new high-order 585 actuator Magellan Adaptive Optics system (MagAO) to obtain very high-resolution visible light images of HD142527 with MagAO's VisAO science camera. In the median seeing conditions of the 6.5m Magellan telescope (0.5-0.7"), we find MagAO delivers 24-19% Strehl at H-alpha (0.656 microns). We detect a faint companion (HD142527B) embedded in this young transitional disk system at just 86.3+/-1.9 mas (~12 AU) from the star. The companion is detected in both H-alpha and a continuum filter (Delta_mag=6.33+/-0.20 mag at H-alpha and 7.50+/-0.25 mag in the continuum filter). This provides confirmation of the tentative companion discovered by Biller and co-workers with sparse aperture masking at the 8m VLT. The H-alpha emission from the ~0.25 solar mass companion (EW=180 Angstroms) implies a mass accretion rate of ~5.9x10^-10 Msun/yr, and a total accretion luminosity of 1.2% Lsun. Assuming a similar accretion rate, we estimate that a 1 Jupiter mass gas giant could have considerably better (50-1,000x) planet/star contrasts at H-alpha than at H band (COND models) for a range of optical extinctions (3.4-0 mag). We suggest that ~0.5-5 Mjup extrasolar planets in their gas accretion phase could be much more luminous at H-alpha than in the NIR. This is the motivation for our new MagAO GAPplanetS survey for extrasolar planets.
[92]  [pdf] - 1202511
LEECH: A 100 Night Exoplanet Imaging Survey at the LBT
Comments: IAUS 299 Proceedings
Submitted: 2014-01-02
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.
[93]  [pdf] - 1202351
Brown Dwarf Photospheres are Patchy: A Hubble Space Telescope Near-infrared Spectroscopic Survey Finds Frequent Low-level Variability
Comments: Accepted for publication in ApJ, 16 pages, plus 25 pages appendix
Submitted: 2013-12-18
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 $\mu$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.
[94]  [pdf] - 1180055
ALMA observations of the molecular gas in the debris disk of the 30 Myr old star HD 21997
Comments: 8 pages, 4 figures, accepted for publication in ApJ
Submitted: 2013-10-18
The 30 Myr old A3-type star HD 21997 is one of the two known debris dust disks having a measurable amount of cold molecular gas. With the goal of understanding the physical state, origin, and evolution of the gas in young debris disks, we obtained CO line observations with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we report on the detection of 12CO and 13CO in the J=2-1 and J=3-2 transitions and C18O in the J=2-1 line. The gas exhibits a Keplerian velocity curve, one of the few direct measurements of Keplerian rotation in young debris disks. The measured CO brightness distribution could be reproduced by a simple star+disk system, whose parameters are r_in < 26 AU, r_out = 138 +/- 20 AU, M_*=1.8 +0.5 -0.2 M_Sun, and i = 32.6 +/- 3.1 degrees. The total CO mass, as calculated from the optically thin C18O line, is about (4-8) x 10^-2 M_Earth, while the CO line ratios suggest a radiation temperature on the order of 6-9 K. Comparing our results with those obtained for the dust component of the HD 21997 disk from the ALMA continuum observations by Mo\'or et al., we conclude that comparable amounts of CO gas and dust are present in the disk. Interestingly, the gas and dust in the HD 21997 system are not co-located, indicating a dust-free inner gas disk within 55 AU of the star. We explore two possible scenarios for the origin of the gas. A secondary origin, which involves gas production from colliding or active planetesimals, would require unreasonably high gas production rates and would not explain why the gas and dust are not co-located. We propose that HD 21997 is a hybrid system where secondary debris dust and primordial gas coexist. HD 21997, whose age exceeds both the model predictions for disk clearing and the ages of the oldest T Tauri-like or transitional gas disks in the literature, may be a key object linking the primordial and the debris phases of disk evolution.
[95]  [pdf] - 1180056
ALMA continuum observations of a 30 Myr old gaseous debris disk around HD 21997
Comments: Accepted for publication in ApJ Letters, 6 pages, 3 figures, 1 tables (emulateapj style)
Submitted: 2013-10-18
Circumstellar disks around stars older than 10 Myr are expected to be gas-poor. There are, however, two examples of old (30-40 Myr) debris-like disks containing a detectable amount of cold CO gas. Here we present ALMA and Herschel Space Observatory observations of one of these disks, around HD 21997, and study the distribution and origin of the dust and its connection to the gas. Our ALMA continuum images at 886um clearly resolve a broad ring of emission within a diameter of ~4.5 arcsec, adding HD 21997 to the dozen debris disks resolved at (sub)millimeter wavelengths. Modeling the morphology of the ALMA image with a radiative transfer code suggests inner and outer radii of ~55 and ~150 AU, and a dust mass of 0.09 M_Earth. Our data and modeling hints at an extended cold outskirt of the ring. Comparison with the morphology of the CO gas in the disk reveals an inner dust-free hole where gas nevertheless can be detected. Based on dust grain lifetimes, we propose that the dust content of this gaseous disk is of secondary origin produced by planetesimals. Since the gas component is probably primordial, HD 21997 is one of the first known examples of a hybrid circumstellar disk, a so-far little studied late phase of circumstellar disk evolution.
[96]  [pdf] - 1179067
A resolved debris disk around the candidate planet-hosting star HD95086
Comments: Accepted for publication in ApJ Letters, 6 pages, 2 figures, 3 tables (emulateapj style)
Submitted: 2013-09-06
Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star HD95086. The strong infrared excess of the system indicates that, similarly to HR8799, {\ss} Pic, and Fomalhaut, the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD95086 disk, and to investigate its possible interaction with the newly discovered planet, here we present new optical, infrared and millimeter observations. We detected no CO emission, excluding the possibility of an evolved gaseous primordial disk. Simple blackbody modeling of the spectral energy distribution suggests the presence of two spatially separate dust belts at radial distances of 6 and 64 AU. Our resolved images obtained with the Herschel Space Observatory reveal a characteristic disk size of ~6.0x5.4 arcsec (540x490 AU) and disk inclination of ~25 degree. Assuming the same inclination for the planet candidate's orbit, its re-projected radial distance from the star is 62 AU, very close to the blackbody radius of the outer cold dust ring. The structure of the planetary system at HD95086 resembles the one around HR8799. Both systems harbor a warm inner dust belt and a broad colder outer disk and giant planet(s) between the two dusty regions. Modelling implies that the candidate planet can dynamically excite the motion of planetesimals even out to 270 AU via their secular perturbation if its orbital eccentricity is larger than about 0.4. Our analysis adds a new example to the three known systems where directly imaged planet(s) and debris disks co-exist.
[97]  [pdf] - 1171661
Protoplanetary Disk Masses from Stars to Brown Dwarfs
Comments: ApJ, accepted (in press)
Submitted: 2013-05-29
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.
[98]  [pdf] - 1165320
HST Spectral Mapping of L/T Transition Brown Dwarfs Reveals Cloud Thickness Variations
Comments: Astrophysical Journal, in press
Submitted: 2013-03-17
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.
[99]  [pdf] - 1165207
Weather on Other Worlds I: Detection of Periodic Variability in the L3 Dwarf DENIS-P J1058.7-1548 with Precise Multi-Wavelength Photometry
Comments: 48 pages, 15 figures, accepted to ApJ
Submitted: 2013-03-12, last modified: 2013-03-14
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.
[100]  [pdf] - 1164874
Benchmark Tests for Markov Chain Monte Carlo Fitting of Exoplanet Eclipse Observations
Comments: 33 pages, 11 tables, 9 figures, accepted to ApJ
Submitted: 2013-02-26
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 $\phi_{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.
[101]  [pdf] - 1159521
The near-infrared spectral energy distribution of {\beta} Pictoris b
Comments: 19 pages, accepted in Astronomy and Astrophysics
Submitted: 2013-02-05
A gas giant planet has previously been directly seen orbiting at 8-10 AU within the debris disk of the ~12 Myr old star {\beta} Pictoris. The {\beta} Pictoris system offers the rare opportunity to study the physical and atmospheric properties of an exoplanet placed on a wide orbit and to establish its formation scenario. We obtained J (1.265 {\mu}m), H (1.66 {\mu}m), and M' (4.78 {\mu}m) band angular differential imaging of the system between 2011 and 2012. We detect the planetary companion in our four-epoch observations. We estimate J = 14.0 +- 0.3, H = 13.5 +- 0.2, and M' = 11.0 +- 0.3 mag. Our new astrometry consolidates previous semi-major axis (sma=8-10 AU) and excentricity (e <= 0.15) estimates of the planet. These constraints, and those derived from radial velocities of the star provides independent upper limits on the mass of {\beta} Pictoris b of 12 and 15.5 MJup for semi-major axis of 9 and 10 AU. The location of {\beta} Pictoris b in color-magnitude diagrams suggests it has spectroscopic properties similar to L0-L4 dwarfs. This enables to derive Log10(L/Lsun) = -3.87 +- 0.08 for the companion. The analysis with 7 PHOENIX-based atmospheric models reveals the planet has a dusty atmosphere with Teff = 1700 +- 100 K and log g = 4.0+- 0.5. "Hot-start" evolutionary models give a new mass of 10+3-2 MJup from Teff and 9+3-2 MJup from luminosity. Predictions of "cold-start" models are inconsistent with independent constraints on the planet mass. "Warm-start" models constrain the mass to M >= 6MJup and the initial entropies to values (Sinit >= 9.3Kb/baryon), intermediate between those considered for cold/hot-start models, but likely closer to those of hot-start models.
[102]  [pdf] - 588245
Vertical Atmospheric Structure in a Variable Brown Dwarf: Pressure-dependent Phase Shifts in Simultaneous Hubble Space Telescope-Spitzer Light Curves
Comments: Accepted to ApJL, 6 pages, 3 figures. Minor language updates from v1 to match published version
Submitted: 2012-10-24, last modified: 2012-11-08
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.
[103]  [pdf] - 1152459
Mapping Directly Imaged Giant Exoplanets
Comments: 51 pages, 14 figures
Submitted: 2012-10-25
With the increasing number of directly imaged giant exoplanets the current atmosphere models are often not capable of fully explaining the spectra and luminosity of the sources. A particularly challenging component of the atmosphere models is the formation and properties of condensate cloud layers, which fundamentally impact the energetics, opacity, and evolution of the planets. Here we present a suite of techniques that can be used to estimate the level of rotational modulations these planets may show. We propose that the time--resolved observations of such periodic photometric and spectroscopic variations of extrasolar planets due to their rotation can be used as a powerful tool to probe the heterogeneity of their optical surfaces. We address and discuss the following questions: a) what planet properties can be deduced from the light curve and/or spectra, and in particular can we determine rotation periods, spot--coverage, spot colors, spot spectra; b) what is the optimal configuration of instrument/wavelength/temporal sampling required for these measurements; and, c) can principal component analysis be used to invert the light curve and deduce the surface map of the planet. Our simulations describe the expected spectral differences between homogeneous (clear or cloudy) and patchy atmospheres, outline the significance of the dominant absorption features of water, methane, and CO and provide a method to distinguish these two types of atmospheres. Simulated photometry from current and future instruments is used to estimate the level of detectable photometric variations. We conclude that future instruments will be able to recover not only the rotation periods, cloud cover, cloud colors and spectra but even cloud evolution. We also show that a longitudinal map of the planet's atmosphere can be deduced from its disk--integrated light curves.
[104]  [pdf] - 1151463
Observational Constraints on the Stellar Radiation Field Impinging on Transitional Disk Atmospheres
Comments: Accepted for publication in ApJ
Submitted: 2012-09-18
Mid-infrared atomic and ionic line ratios measured in spectra of pre-main sequence stars are sensitive indicators of the hardness of the radiation field impinging on the disk surface. We present a low-resolution Spitzer IRS search for [Ar II] at 6.98 $\mu$m, [Ne II] at 12.81 $\mu$m, and [Ne III] 15.55 $\mu$m lines in 56 transitional disks. These objects, characterized by reduced near-infrared but strong far-infrared excess emission, are ideal targets to set constraints on the stellar radiation field onto the disk because their spectra are not contaminated by shock emission from jets/outflows or by molecular emission lines. After demonstrating that we can detect [Ne II] lines and recover their fluxes from the low-resolution spectra, here we report the first detections of [Ar II] lines towards protoplanetary disks. We did not detect [Ne III] emission in any of our sources. Our [Ne II]/[Ne III] line flux ratios combined with literature data suggest that a soft-EUV or X-ray spectrum produces these gas lines. Furthermore, the [Ar II]/[Ne II] line flux ratios point to a soft X-ray and/or soft-EUV stellar spectrum as the ionization source of the [Ar II] and [Ne II] emitting layer of the disk. If the soft X-ray component dominates over the EUV than we would expect larger photoevaporation rates hence a reduction of the time available to form planets.
[105]  [pdf] - 1117240
The Gray Needle: Large Grains in the HD 15115 Debris Disk from LBT/PISCES/Ks and LBTI/LMIRcam/L' Adaptive Optics Imaging
Comments: Apj-accepted March 27 2012; minor corrections
Submitted: 2012-03-12, last modified: 2012-05-28
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.
[106]  [pdf] - 1123342
An HST Imaging Survey of Low-Mass Stars in the Chamaeleon I Star Forming region
Comments: to appear on Astronomical Journal, accepted April 16, 2012 (AJ-10740)
Submitted: 2012-05-11
We present new HST/WFPC2 observations of 20 fields centered around T Tauri stars in the Chamaeleon I star forming region. Images have been obtained in the F631N ([OI]6300A), F656N (Ha) and F673N ([SII]6716A+6731A) narrow-band filters, plus the Johnson V-band equivalent F547M filter. We detect 31 T Tauri stars falling within our fields. We discuss the optical morphology of 10 sources showing evidence of either binarity, circumstellar material, or mass loss. We supplement our photometry with a compilation of optical, infrared and sub-millimeter data from the literature, together with new sub-mm data for three objects, to build the Spectral Energy Distributions (SED) of 19 single sources. Using an SED model fitting tool, we self-consistently estimate a number of stellar and disk parameters, while mass accretion rates are directly derived from our Ha photometry. We find that bolometric luminosities derived from dereddened optical data tend to be underestimated in systems with high alpha(2-24} IR spectral index, suggesting that disks seen nearly edge-on may occasionally be interpreted as low luminosity (and therefore more evolved) sources. On the other hand, the same alpha(2-24) spectral index, a tracer of the amount of dust in the warmer layers of the circumstellar disks, and the mass accretion rate appear to decay with the isocronal stellar age, suggesting that the observed age spread (~0.5-5 Myr) within the cluster is real. Our sample contains a few outliers that may have dissipated their circumstellar disks on shorter time-scale.
[107]  [pdf] - 1117237
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
Comments: Accepted to ApJ
Submitted: 2012-03-12, last modified: 2012-04-27
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.
[108]  [pdf] - 1116586
Orbital characterization of the \beta Pictoris b giant planet
Comments: 10 pages, 12 figures, accepted to A&A
Submitted: 2012-02-13
In June 2010, we confirmed the existence of a giant planet in the disk of the young star Beta Pictoris, located between 8 AU and 15 AU from the star. This young planet offers the rare opportunity to monitor a large fraction of the orbit using the imaging technique over a reasonably short timescale. Using the NAOS-CONICA adaptive-optics instrument (NACO) at the Very Large Telescope (VLT), we obtained repeated follow-up images of the Bpic system in the Ks and L' filters at four new epochs in 2010 and 2011. Complementing these data with previous measurements, we conduct a homogeneous analysis, which covers more than eight yrs, to accurately monitor the Bpic b position relative to the star. On the basis of the evolution of the planet's relative position with time, we derive the best-fit orbital solutions for our measurements. More reliable results are found with a Markov-chain Monte Carlo approach. The solutions favor a low-eccentricity orbit e < 0.17, with semi-major axis in the range 8--9 AU corresponding to orbital periods of 17--21 yrs. Our solutions favor a highly inclined solution with a peak around i=88.5+-1.7 deg, and a longitude of ascending node tightly constrained at Omega = -147.5+-1.5 deg. These results indicate that the orbital plane of the planet is likely to be above the midplane of the main disk, and compatible with the warp component of the disk being tilted between 3.5 deg and 4.0 deg. This suggests that the planet plays a key role in the origin of the inner warped-disk morphology of the Bpic disk. Finally, these orbital parameters are consistent with the hypothesis that the planet is responsible for the transit-like event observed in November 1981, and also linked to the cometary activity observed in the Bpic system.
[109]  [pdf] - 1116576
beta Pic b position relative to the Debris Disk
Comments: 13 pages, 6 figures, to appear in Astronomy and Astrophysics
Submitted: 2012-02-12
Context. We detected in 2009 a giant, close-by planet orbiting {\beta} Pic, a young star surrounded with a disk, extensively studied for more than 20 years. We showed that if located on an inclined orbit, the planet could explain several peculiarities of {\beta} Pictoris system. However, the available data did not permit to measure the inclination of {\beta} Pic b with respect to the disk, and in particular to establish in which component of the disk - the main, extended disk or the inner inclined component/disk-, the planet was located. Comparison between the observed planet position and the disk orientation measured on previous imaging data was not an option because of potential biases in the measurements. Aims. Our aim is to measure precisely the planet location with respect to the dust disk using a single high resolution image, and correcting for systematics or errors that degrades the precision of the disk and planet relative position measurements. Methods. We gathered new NaCo data at Ks band, with a set-up optimized to derive simultaneously the orientation(s) of the disk(s) and that of the planet. Results. We show that the projected position of {\beta} Pic b is above the midplane of the main disk. With the current data and knowledge on the system, this implies that {\beta} Pic b cannot be located in the main disk. The data rather suggest the planet being located in the inclined component.
[110]  [pdf] - 1092259
Resolving the inner regions of the HD97048 circumstellar disk with VLT/NACO polarimetric differential imaging
Comments: Accepted for publication in Astronomy & Astrophysics; 6 pages incl. 5 figures and 2 tables
Submitted: 2011-12-08
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.
[111]  [pdf] - 1083986
Molecular gas in young debris disks
Comments: Accepted for publication in ApJ Letters, 8 pages, 3 figures, 1 tables (emulateapj style)
Submitted: 2011-09-11
Gas-rich primordial disks and tenuous gas-poor debris disks are usually considered as two distinct evolutionary phases of the circumstellar matter. Interestingly, the debris disk around the young main-sequence star 49 Ceti possesses a substantial amount of molecular gas, and possibly represents the missing link between the two phases. Motivated to understand the evolution of the gas component in circumstellar disks via finding more 49 Ceti-like systems, we carried out a CO J=3-2 survey with Atacama Pathfinder EXperiment, targeting 20 infrared-luminous debris disks. These systems fill the gap between primordial and old tenuous debris disks in terms of fractional luminosity. Here we report on the discovery of a second 49 Ceti-like disk around the 30 Myr old A3-type star HD21997, a member of the Columba Association. This system was also detected in the CO(2-1) transition, and the reliable age determination makes it an even clearer example of an old gas-bearing disk than 49 Ceti. While the fractional luminosities of HD21997 and 49 Ceti are not particularly high, these objects seem to harbor the most extended disks within our sample. The double-peaked profiles of HD21997 were reproduced by a Keplerian disk model combined with the LIME radiative transfer code. Based on their similarities, 49 Ceti and HD21997 may be the first representatives of a so far undefined new class of relatively old (>~8 Myr), gaseous dust disks. From our results, neither primordial origin nor steady secondary production from icy planetesimals can unequivocally explain the presence of CO gas in the disk of HD21997.
[112]  [pdf] - 1053331
A peculiar young eruptive star in the dark cloud Lynds 1340
Comments: 15 pages, 4 figures, 1 table; accepted by ApJL
Submitted: 2011-04-13
We conducted a long-term optical photometric and spectroscopic monitoring of the strongly variable, accreting young sun-like star [KOS94] HA11, associated with the dark cloud Lynds 1340, that exhibited large amplitude (5-6 magnitudes in the I_C band) brightness variations on 2-3 years timescales, flat spectral energy distribution (SED), and extremely strong (300 < EW/Angstrom < 900) H alpha emission. In this Letter we describe the basic properties of the star, derived from our observations between 1999 and 2011, and put into context the observed phenomena. The observed variations in the emission spectra, near-infrared colors, and SED suggest that [KOS94] HA11 (spectral type: K7--M0) is an eruptive young star, possibly similar in nature to V1647 Ori: its large-scale photometric variations are governed by variable accretion rate, associated with variations in the inner disk structure. The star recently has undergone strong and rapid brightness variations, thus its further observations may offer a rare opportunity for studying structural and chemical rearrangements of the inner disk, induced by variable central luminosity.
[113]  [pdf] - 1053110
Measuring Organic Molecular Emission in Disks with Low Resolution Spitzer Spectroscopy
Submitted: 2011-04-01
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.
[114]  [pdf] - 1453874
Inner disc rearrangement revealed by dramatic brightness variations in the young star PV Cep
Comments: 11 pages, 4 figures, accepted for publication by MNRAS. 3 online tables added
Submitted: 2011-01-12, last modified: 2011-03-12
Young Sun-like stars at the beginning of the pre-main sequence (PMS) evolution are surrounded by accretion discs and remnant protostellar envelopes. Photometric and spectroscopic variations of these stars are driven by interactions of the star with the disc. Time scales and wavelength dependence of the variability carry information on the physical mechanisms behind these interactions. We conducted multi-epoch, multi-wavelength study of PV Cep, a strongly variable, accreting PMS star. By combining our own observations from 2004-2010 with archival and literature data, we show that PV Cep started a spectacular fading in 2005, reaching an I_C-band amplitude of 4 mag. Analysis of variation of the optical and infrared fluxes, colour indices, and emission line fluxes suggests that the photometric decline in 2005-2009 resulted from an interplay between variable accretion and circumstellar extinction: since the central luminosity of the system is dominated by accretion, a modest drop in the accretion rate could induce the drastic restructuring of the inner disc. Dust condensation in the inner disc region might have resulted in the enhancement of the circumstellar extinction.
[115]  [pdf] - 1453882
Evidence Against an Edge-On Disk Around the Extrasolar Planet 2MASS 1207 b and a New Thick Cloud Explanation for its Under-Luminosity
Comments: 35 pages, 9 figures, accepted to ApJ
Submitted: 2011-03-09
(Abridged) Since the discovery of the first directly-imaged, planetary-mass object, 2MASS 1207 b, several works have sought to explain a disparity between its observed temperature and luminosity. Given its known age, distance, and spectral type, 2MASS 1207 b is under-luminous by a factor of ~10 (~2.5 mags) when compared to standard models of brown-dwarf/giant-planet evolution. In this paper, we study three possible sources of 2MASS 1207 b's under-luminosity. First, we investigate Mohanty et al. (2007)'s hypothesis that a near edge-on disk might be responsible for 2MASS 1207 b's under-luminosity. We conclude that the hypothesis is unlikely due to the lack of variability seen in multi-epoch photometry and unnecessary due to the increasing sample of under-luminous brown-dwarfs/giant-exoplanets that cannot be explained by an edge-on disk. Next, we test the analogous possibility that a spherical shell of dust, could explain 2MASS 1207 b's under-luminosity. Models containing enough dust to create ~2.5 mags of extinction, placed at reasonable radii, are ruled out by our new Gemini/T-ReCS 8.7 micron photometric upper-limit for 2MASS 1207 b. Finally, we investigate the possibility that 2MASS 1207 b is intrinsically cooler than the commonly used AMES-DUSTY fits to its spectrum, and thus it is not, in fact, under-luminous. New, thick cloud model grids by Madhusudhan et al. (2011) fit 2MASS 1207 b's 1-10 micron SED well, but they do not quite fit its near-infrared spectrum. However, we suggest that with some "tuning", they might be capable of simultaneously reproducing 2MASS 1207 b's spectral shape and luminosity. In this case, the whole class of young, under-luminous brown-dwarfs/giant-exoplanets might be explained by atmospheres that are able to suspend thick, dusty clouds in their photospheres at lower temperatures than field brown-dwarfs.
[116]  [pdf] - 1051309
A Combined Subaru/VLT/MMT 1--5 Micron Study of Planets Orbiting HR 8799: Implications for Atmospheric Properties, Masses, and Formation
Comments: 27 pages, 15 figures, Accepted for publication in The Astrophysical Journal
Submitted: 2011-01-10, last modified: 2011-01-18
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.
[117]  [pdf] - 1042614
Structure and evolution of debris disks around F-type stars: I. Observations, database and basic evolutionary aspects
Comments: Accepted for publication in ApJS, 42 pages, 18 figures, 8 tables (emulateapj style)
Submitted: 2010-12-16
Although photometric and spectroscopic surveys with the Spitzer Space Telescope increased remarkably the number of well studied debris disks around A-type and Sun-like stars, detailed analyzes of debris disks around F-type stars remained less frequent. Using the MIPS camera and the IRS spectrograph we searched for debris dust around 82 F-type stars with Spitzer. We found 27 stars that harbor debris disks, nine of which are new discoveries. The dust distribution around two of our stars, HD 50571 and HD 170773, was found to be marginally extended on the 70um MIPS images. Combining the MIPS and IRS measurements with additional infrared and submillimeter data, we achieved excellent spectral coverage for most of our debris systems. We have modeled the excess emission of 22 debris disks using a single temperature dust ring model and of 5 debris systems with two-temperature models. The latter systems may contain two dust rings around the star. In accordance with the expected trends, the fractional luminosity of the disks declines with time, exhibiting a decay rate consistent with the range of model predictions. We found the distribution of radial dust distances as a function of age to be consistent with the predictions of both the self stirred and the planetary stirred disk evolution models. A more comprehensive investigation of the evolution of debris disks around F-type stars, partly based on the presented data set, will be the subject of an upcoming paper.
[118]  [pdf] - 1034665
First Results From VLT NACO Apodizing Phase Plate: 4-micron Images of the Exoplanet beta Pictoris b
Comments: Accepted for publication by ApJL
Submitted: 2010-09-02
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.
[119]  [pdf] - 1394950
Stellar-Mass-Dependent Disk Structure in Coeval Planet-Forming Disks
Comments: Astrophysical Journal, in press
Submitted: 2010-07-21
Previous studies suggest that the planet-forming disks around very-low-mass stars/brown dwarfs may be flatter than those around more massive stars, in contrast to model predictions of larger scale heights for gas-disks around lower-mass stars. We conducted a statistically robust study to determine whether there is evidence for stellar-mass-dependent disk structure in planet-forming disks. We find a statistically significant difference in the Spitzer/IRAC color distributions of disks around very-low-mass and low-mass stars all belonging to the same star-forming region, the Chamaeleon I star-forming region. We show that self consistently calculated disk models cannot fit the median spectral energy distributions (SEDs) of the two groups. These SEDs can be only explained by flatter disk models, consistent with the effect of dust settling in disks. We find that relative to the disk structure predicted for flared disks the required reduction in disk scale height is anti-correlated with the stellar mass, i.e. disks around lower-mass stars are flatter. Our results show that the initial and boundary conditions of planet formation are stellar-mass-dependent, an important finding that must be considered in planet formation models.
[120]  [pdf] - 1033148
A giant planet imaged in the disk of the young star Beta Pictoris
Comments: 4 pages, 2 figures. Published online 10 June 2010; 10.1126/science.1187187. To appear in Science
Submitted: 2010-06-16
Here we show that the ~10 Myr Beta Pictoris system hosts a massive giant planet, Beta Pictoris b, located 8 to 15 AU from the star. This result confirms that gas giant planets form rapidly within disks and validates the use of disk structures as fingerprints of embedded planets. Among the few planets already imaged, Beta Pictoris b is the closest to its parent star. Its short period could allow recording the full orbit within 17 years.
[121]  [pdf] - 1018739
Day-side z'-band emission and eccentricity of Wasp-12b
Comments: One more author added. Version accepted for publication on ApJL
Submitted: 2009-12-14, last modified: 2010-05-20
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.
[122]  [pdf] - 315997
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
Comments: accepted for publication on ApJ
Submitted: 2009-10-06
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{\omega} = 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.
[123]  [pdf] - 1002747
Constraining the orbit of the possible companion to Beta Pictoris: New deep imaging observations
Comments: 8 pages, 8 figures, 1 table, accepted for publication in Astronomy and Astrophysics
Submitted: 2009-06-30, last modified: 2009-07-01
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.
[124]  [pdf] - 1002653
Imaging search for the unseen companion to Eps Ind A -- Improving the detection limits with 4 micron observations
Comments: 10 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2009-06-22
Eps Ind A is one of the nearest sun-like stars, located only 3.6 pc away. It is known to host a binary brown dwarf companion, Eps Ind Ba/Bb, at a large projected separation of 6.7", but radial velocity measurements imply that an additional, yet unseen component is present in the system, much closer to Eps Ind A. Previous direct imaging has excluded the presence of any stellar or high-mass brown dwarf companion at small separations, indicating that the unseen companion may be a low-mass brown dwarf or high-mass planet. We present the results of a deep high-contrast imaging search for the companion, using active angular differential imaging (aADI) at 4 micron, a particularly powerful technique for planet searches around nearby and relatively old stars. We also develop an additional PSF reference subtraction scheme based on locally optimized combination of images (LOCI) to further enhance the detection limits. No companion is seen in the images, although we are sensitive to significantly lower masses than previously achieved. Combining the imaging data with the known radial velocity trend, we constrain the properties of the companion to within approximately 5-20 Mjup, 10-20 AU, and i > 20 deg, unless it is an exotic stellar remnant. The results also imply that the system is probably older than the frequently assumed age of ~1 Gyr.
[125]  [pdf] - 24998
The discovery of new warm debris disks around F-type stars
Comments: Accepted for publication in ApJ Letters, 13 pages, 2 figures, 2 tables
Submitted: 2009-06-05
We report the discovery of four rare debris disks with warm excesses around F-stars, significantly increasing the number of such systems known in the solar neighborhood. Three of the disks are consistent with the predictions of steady state planetesimal disk evolution models. The oldest source, HD 169666, displays a dust fractional luminosity too high to be in steady state and we suggest that this system recently underwent a transient event of dust production. In addition, two spectra of this star separated by ~three years show silicate emission features, indicative of submicron- to micron-sized grains. We argue that such small grains would be rapidly depleted and their presence in both spectra suggests that the production of small dust is continuous over at least on few years timescale. We predict that systems showing variable mid-infrared spectra, if exist, will provide valuable help in distinguishing the possible scenarios proposed for dust replenishment.
[126]  [pdf] - 22711
The low level of debris disk activity at the time of the Late Heavy Bombardment: a Spitzer study of Praesepe
Comments: 22 pages, 14 figures, 9 tables, emulateapj format; Accepted for publication in The Astrophysical Journal
Submitted: 2009-03-24
We present 24 micron photometry of the intermediate-age open cluster Praesepe. We assemble a catalog of 193 probable cluster members that are detected in optical databases, the Two Micron All Sky Survey (2MASS), and at 24 micron, within an area of ~ 2.47 square degrees. Mid-IR excesses indicating debris disks are found for one early-type and for three solar-type stars. Corrections for sampling statistics yield a 24 micron excess fraction (debris disk fraction) of 6.5 +- 4.1% for luminous and 1.9 +- 1.2% for solar-type stars. The incidence of excesses is in agreement with the decay trend of debris disks as a function of age observed for other cluster and field stars. The values also agree with those for older stars, indicating that debris generation in the zones that emit at 24 micron falls to the older 1-10 Gyr field star sample value by roughly 750 Myr. We discuss our results in the context of previous observations of excess fractions for early- and solar-type stars. We show that solar-type stars lose their debris disk 24 micron excesses on a shorter timescale than early-type stars. Simplistic Monte Carlo models suggest that, during the first Gyr of their evolution, up to 15-30% of solar-type stars might undergo an orbital realignment of giant planets such as the one thought to have led to the Late Heavy Bombardment, if the length of the bombardment episode is similar to the one thought to have happened in our Solar System. In the Appendix, we determine the cluster's parameters via boostrap Monte Carlo isochrone fitting, yielding an age of 757 Myr (+- 36 Myr at 1 sigma confidence) and a distance of 179 pc (+- 2 pc at 1 sigma confidence), not allowing for systematic errors.
[127]  [pdf] - 22377
Bridging the Gap Between Stars and Planets: The Formation and Early Evolution of Brown Dwarfs
Comments: 8 pages, 3 figures, Astro2010 Science White Paper (v2 has 2 references corrected/updated, and 3 repeated references in the v1 reference list removed)
Submitted: 2009-03-13, last modified: 2009-03-19
This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.
[128]  [pdf] - 900599
Do we really know the dust? Systematics and uncertainties of the mid-infrared spectral analysis methods
Comments: Accepted for publication in ApJ, 47 pages, 13 figures, 3 tables
Submitted: 2009-02-03
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.
[129]  [pdf] - 20249
MBM 12: young protoplanetary discs at high galactic latitude
Comments: 14 pages, accepted by Astronomy and Astrophysics
Submitted: 2009-01-12
(abridged) We present Spitzer infrared observations to constrain disc and dust evolution in young T Tauri stars in MBM 12, a star-forming cloud at high latitude with an age of 2 Myr and a distance of 275 pc. The region contains 12 T Tauri systems, with primary spectral types between K3 and M6; 5 are weak-line and the rest classical T Tauri stars. We first use MIPS and literature photometry to compile spectral energy distributions for each of the 12 members in MBM 12, and derive their IR excesses. The IRS spectra are analysed with the newly developed two-layer temperature distribution (TLTD) spectral decomposition method. For the 7 T Tauri stars with a detected IR excess, we analyse their solid-state features to derive dust properties such as mass-averaged grain size, composition and crystallinity. We find a spatial gradient in the forsterite to enstatite range, with more enstatite present in the warmer regions. The fact that we see a radial dependence of the dust properties indicates that radial mixing is not very efficient in the discs of these young T Tauri stars. The SED analysis shows that the discs in MBM 12, in general, undergo rapid inner disc clearing, while the binary sources have faster discevolution. The dust grains seem to evolve independently from the stellar properties, but are mildly related to disc properties such as flaring and accretion rates.
[130]  [pdf] - 315187
The Different Evolution of Gas and Dust in Disks around Sun-like and Cool Stars
Comments: submitted to the ApJ
Submitted: 2008-10-14, last modified: 2009-01-05
We present the first comprehensive, comparative Spitzer/IRS study of the dust and gas properties of disks around young sun-like stars (K1-M5) and cool stars/brown dwarfs (M5-M9). The comparison of these two large samples of over 60 sources reveal major differences in the evolution of both the dust and gas components.We report the first detections of organic molecules in disks around brown dwarfs. The detection rate statistics and the line flux ratios of HCN and C2H2 show a striking diference between the two samples, demonstrating a significant under-abundance of HCN relative to C2H2 in the disk surface of cool stars. We propose this to originate from the large difference in the UV-irradiation around the two types of sources. The statistical comparison of the 10 micron silicate emission features also reveals a significant difference between the two samples. Cool stars and brown dwarfs show weaker features arising from more processed silicate grains in the disk atmosphere. These findings complement previous indications of flatter disk structures and longer disk lifetimes around cool stars.
[131]  [pdf] - 6042
The Outer Cut-Off of the Giant Planet Population and the 6pc-Survey
Comments: To appear in the Conference Proceedings "Extreme Solar Systems"; 3 pages, 1 figure
Submitted: 2007-10-15
We present results from two high-contrast imaging surveys that exploit a novel technique, L-band angular differential imaging. Our first survey targeted 21 young stars in the Beta Pic and Tuc-Hor moving groups with VLT/NACO reaching typical sensitivities of <1 MJup at r > 20 AU. The statistical analysis of the null result demonstrates that the giant planet population is truncated at 30 AU or less (90% confidence level). Our second, on-going MMT/Clio survey utilizes the unique sensitivity achieved in the L-band for old planets to probe all M-dwarf stars within 6 pc. The proximity of these targets enables direct imaging of planets in orbits like Jupiter for the first time - a key step for directly imaging giant planets.
[132]  [pdf] - 5764
High-resolution polarimetry of Parsamian 21: revealing the structure of an edge-on FU Ori disc
Comments: Accepted for publication in the MNRAS. 16 pages, 18 figures and 5 tables
Submitted: 2007-10-07
We present the first high spatial resolution near-infrared direct and polarimetric observations of Parsamian 21, obtained with the VLT/NACO instrument. We complemented these measurements with archival infrared observations, such as HST/WFPC2 imaging, HST/NICMOS polarimetry, Spitzer IRAC and MIPS photometry, Spitzer IRS spectroscopy as well as ISO photometry. Our main conclusions are the following: (1) we argue that Parsamian 21 is probably an FU Orionis-type object; (2) Parsamian 21 is not associated with any rich cluster of young stars; (3) our measurements reveal a circumstellar envelope, a polar cavity and an edge-on disc; the disc seems to be geometrically flat and extends from approximately 48 to 360 AU from the star; (4) the SED can be reproduced with a simple model of a circumstellar disc and an envelope; (5) within the framework of an evolutionary sequence of FUors proposed by Green et al. (2006) and Quanz et al. (2007), Parsamian 21 can be classified as an intermediate-aged object.
[133]  [pdf] - 5515
Medium-separation binaries do not affect the first steps of planet formation
Comments: To appear in the Astrophysical Journal
Submitted: 2007-10-01, last modified: 2007-10-05
The first steps of planet formation are marked by the growth and crystallization of sub-micrometer-sized dust grains accompanied by dust settling toward the disk midplane. In this paper we explore whether the first steps of planet formation are affected by the presence of medium-separation stellar companions. We selected two large samples of disks around single and binary T Tauri stars in Taurus that are thought to have only a modest age spread of a few Myr. The companions of our binary sample are at projected separations between 10 and 450 AU with masses down to about 0.1 solar masses. We used the strength and shape of the 10 micron silicate emission feature as a proxy for grain growth and for crystallization respectively. The degree of dust settling was evaluated from the ratio of fluxes at two different mid-infrared wavelengths. We find no statistically significant difference between the distribution of 10 micron silicate emission features from single and binary systems. In addition, the distribution of disk flaring is indistinguishable between the single and binary system samples. These results show that the first steps of planet formation are not affected by the presence of a companion at tens of AU.
[134]  [pdf] - 5514
A Survey for Massive Giant Planets in Debris Disks with Evacuated Inner Cavities
Comments: Accepted for publication in ApJ
Submitted: 2007-10-01
The commonality of collisionally replenished debris around main sequence stars suggests that minor bodies are frequent around Sun-like stars. Whether or not debris disks in general are accompanied by planets is yet unknown, but debris disks with large inner cavities - perhaps dynamically cleared - are considered to be prime candidates for hosting large-separation massive giant planets. We present here a high-contrast VLT/NACO angular differential imaging survey for eight such cold debris disks. We investigated the presence of massive giant planets in the range of orbital radii where the inner edge of the dust debris is expected. Our observations are sensitive to planets and brown dwarfs with masses >3 to 7 Jupiter mass, depending on the age and distance of the target star. Our observations did not identify any planet candidates. We compare the derived planet mass upper limits to the minimum planet mass required to dynamically clear the inner disks. While we cannot exclude that single giant planets are responsible for clearing out the inner debris disks, our observations constrain the parameter space available for such planets. The non-detection of massive planets in these evacuated debris disks further reinforces the notion that the giant planet population is confined to the inner disk (<15 AU).
[135]  [pdf] - 1810
A novel L-band imaging search for giant planets in the Tucana and Beta Pictoris moving groups
Comments: 9 pages, 9 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2007-06-01
The paper presents results from a small survey of 22 young, nearby stars that was designed to detect substellar companions and ultimately giant extrasolar planets down to Jupiter masses. The targets are members of the Tucana and Beta Pictoris moving groups apart from the somewhat older star HIP 71395 that has a radial velocity trend suggesting a massive planet in large orbit. The survey was carried out in the L-band using adaptive optics assisted imaging with NAOS-CONICA (NACO) at the VLT. The chosen observation wavelength is well suited to search for close companions around young stars and delivers unprecedented detection limits. The presented technique reaches some of the best sensitivities as of today and is currently the most sensitive method for the contrast limited detection of substellar companions that are cooler than about 1000K. The companion to 51 Eri, GJ 3305, was found to be a very close binary on an eccentric orbit. No substellar companions were found around the target stars, although the method permitted to detect companions down to a few Jupiter masses at orbital distances of typically 5 astronomical units. A planet with a mass >1 M_Jup at distances >5 AU around AU Mic can be excluded at the time of our observations. The absence of detected planets sets constraints on the frequency distribution and maximum orbital distance of giant exoplanets. For example, a radial distribution power law index of 0.2 in combination with a maximum orbital radius exceeding 30 AU can be rejected at a 90% confidence level.
[136]  [pdf] - 89340
Disks around Brown Dwarfs and Cool Stars
Comments: Summary of the Cool Stars 14 Splinter Session "Disks around Brown Dwarfs and Cool Stars". 8 pages, 1 figure
Submitted: 2007-02-10
We review the current picture of disks around cool stars and brown dwarfs, including disk fractions, mass estimates, disk structure and dispersal, accretion, dust composition, and the debris disk phase. We discuss these in the framework of recent planet formation models.
[137]  [pdf] - 86210
Dust rings and filaments around the isolated young star V1331 Cygni
Comments: 15 pages ApJ preprint style including 3 figures, accepted for publication in ApJ (Feb. 2007)
Submitted: 2006-10-26
We characterize the small and large scale environment of the young star V1331 Cygni with high resolution HST/WFPC2 and Digitized Sky Survey images. In addition to a previously known outer dust ring (~30'' in diameter), the HST/WFPC2 scattered light image reveals an inner dust ring for the first time. This ring has a maximum radius of 6.5'' and is possibly related to a molecular envelope. Large-scale optical images show that V1331 Cyg is located at the tip of a long dust filament linking it to the dark cloud LDN 981. We discuss the origin of the observed dust morphology and analyze the object's relation to its parent dark cloud LDN 981. Finally, based on recent results from the literature, we investigate the properties of V1331 Cyg and conclude that in its current state the object does not show suffcient evidence to be characterized as an FU Ori object.
[138]  [pdf] - 85509
Massive binaries in high-mass star-forming regions: A multi-epoch radial velocity survey of embedded O-stars
Comments: To appear in the Astrophysical Journal (January 10, 2007, v654n 2 issue)
Submitted: 2006-10-03
We present the first multi-epoch radial velocity study of embedded young massive stars using near-infrared spectra obtained with ISAAC mounted at the ESO Very Large Telescope, with the aim to detect massive binaries. Our 16 targets are located in high-mass star-forming regions and many of them are associated to known ultracompact HI I regions, whose young age ensures that dynamic evolution of the clusters did not influence the intrinsic binarity rate. We identify two stars with about 90 km/s velocity differences between two epochs proving the presence of close massive binaries. The fact that 2 out of the 16 observed stars are binary systems suggests that at least 20% of the young massive stars are formed in close multiple systems, but may also be consistent with most if not all young massive stars being binaries. In addition, we show that the radial velocity dispersion of the full sample is about 35 km/s, significantly larger than our estimated uncertainty (25 km/s). This finding is consistent with similar measurements of the young massive cluster 30 Dor which might have a high intrinsic binary rate. Furthermore, we argue that virial cluster masses derived from the radial velocity dispersion of young massive stars may intrinsically overestimate the cluster mass due to the presence of binaries.
[139]  [pdf] - 82896
High Contrast L' Band Adaptive Optics Imaging to Detect Extrasolar Planets
Comments: 8 pages, 4 figures
Submitted: 2006-06-19
We are carrying out a survey to search for giant extrasolar planets around nearby, moderate-age stars in the mid-infrared L' and M bands (3.8 and 4.8 microns, respectively), using the Clio camera with the adaptive optics system on the MMT telescope. To date we have observed 7 stars, of a total 50 planned, including GJ 450 (distance about 8.55pc, age about 1 billion years, no real companions detected), which we use as our example here. We report the methods we use to obtain extremely high contrast imaging in L', and the performance we have obtained. We find that the rotation of a celestial object over time with respect to a telescope tracking it with an altazimuth mount can be a powerful tool for subtracting telescope-related stellar halo artifacts and detecting planets near bright stars. We have carried out a thorough Monte Carlo simulation demonstrating our ability to detect planets as small as 6 Jupiter masses around GJ 450. The division of a science data set into two independent parts, with companions required to be detected on both in order to be recognized as real, played a crucial role in detecting companions in this simulation. We mention also our discovery of a previously unknown faint stellar companion to another of our survey targets, HD 133002. Followup is needed to confirm this as a physical companion, and to determine its physical properties.
[140]  [pdf] - 80938
Nearby debris disk systems with high fractional luminosity reconsidered
Comments: 47 pages, 6 figures, AASTeX preprint (accepted for publication in the Astrophysical Journal); refined the text parts of table 6 and 7
Submitted: 2006-03-27, last modified: 2006-03-31
By searching the IRAS and ISO databases we compiled a list of 60 debris disks which exhibit the highest fractional luminosity values (fd>10^-4) in the vicinity of the Sun (d<120pc). Eleven out of these 60 systems are new discoveries. Special care was taken to exclude bogus disks from the sample. We computed the fractional luminosity values using available IRAS, ISO, and Spitzer data, and analysed the galactic space velocities of the objects. The results revealed that stars with disks of high fractional luminosity often belong to young stellar kinematic groups, providing an opportunity to obtain improved age estimates for these systems. We found that practically all disks with fd>5x10^-4 are younger than 100Myr. The distribution of the disks in the fractional luminosity versus age diagram indicates that (1) the number of old systems with high fd is lower than was claimed before; (2) there exist many relatively young disks of moderate fractional luminosity; and (3) comparing the observations with a current theoretical model of debris disk evolution a general good agreement could be found.
[141]  [pdf] - 79730
Crystalline silicates as a probe of disk formation history
Comments: Accepted for publication in ApJL
Submitted: 2006-02-07
We present a new perspective on the crystallinity of dust in protoplanetary disks. The dominant crystallization by thermal annealing happens in the very early phases of disk formation and evolution. Both the disk properties and the level of crystallinity are thereby directly linked to the properties of the molecular cloud core from which the star+disk system was formed. We show that, under the assumption of single star formation, rapidly rotating clouds produce disks which, after the main infall phase (i.e. in the optically revealed class II phase), are rather massive and have a high accretion rate but low crystallinity. Slowly rotating clouds, on the other hand, produce less massive disks with lower accretion rate, but high levels of crystallinity. Cloud fragmentation and the formation of multiple stars complicates the problem and necessitates further study. The underlying physics of the model is insufficiently understood to provide the precise relationship between crystallinity, disk mass and accretion rate. But the fact that with `standard' input physics the model produces disks which, in comparison to observations, appear to have either too high levels of crystallinity or too high disk masses, demonstrates that the comparison of these models to observations can place strong contraints on the disk physics. The question to ask is not why some sources are so crystalline, but why some other sources have such a low level of crystallinity.
[142]  [pdf] - 77778
The Onset of Planet Formation in Brown Dwarf Disks
Comments: Published in Science 2005, vol 310, 834; 3 pages in final format, 4 figures + 8 pages Supporting Online Material. For final typeset, see
Submitted: 2005-11-14
The onset of planet formation in protoplanetary disks is marked by the growth and crystallization of sub-micron-sized dust grains accompanied by dust settling toward the disk mid-plane. Here we present infrared spectra of disks around brown dwarfs and brown dwarf candidates. We show that all three processes occur in such cool disks in a way similar or identical to that in disks around low- and intermediate-mass stars. These results indicate that the onset of planet formation extends to disks around brown dwarfs, suggesting that planet formation is a robust process occurring in most young circumstellar disks.
[143]  [pdf] - 71275
HST/NICMOS observations of a proto-brown dwarf candidate
Comments: Astronomy & Astrophysics Letters, in press; 4 pages, 2 figures
Submitted: 2005-02-23
We present deep HST/NICMOS observations peering through the outflow cavity of the protostellar candidate IRAS 04381+2540 in the Taurus Molecular Cloud-1. A young stellar object as central source, a jet and a very faint and close (0.6") companion are identified. The primary and the companion have similar colours, consistent with strong reddening. We argue that the companion is neither a shock-excited knot nor a background star. The colour/magnitude information predicts a substellar upper mass limit for the companion, but the final confirmation will require spectroscopic information. Because of its geometry, young age and its rare low-mass companion, this system is likely to provide a unique insight into the formation of brown dwarfs.
[144]  [pdf] - 67396
Grain growth and dust settling in a brown dwarf disk: Gemini/T-ReCS observations of CFHT-BD-Tau 4
Comments: Accepted for publication in Astronomy & Astrophysics Letters, 4.5 pages with 1 figure
Submitted: 2004-09-13
We present accurate mid-infrared observations of the disk around the young, bona-fide brown dwarf CFHT-BD-Tau 4. We report GEMINI/T-ReCS measurements in the 7.9, 10.4 and 12.3 micron filters, from which we infer the presence of a prominent, broad silicate emission feature. The shape of the silicate feature is dominated by emission from 2 micron amorphous olivine grains. Such grains, being an order of magnitude larger than those in the interstellar medium, are a first proof of dust processing and grain growth in disks around brown dwarfs. The object's spectral energy distribution is below the prediction of the classical flared disk model but higher than that of the two-layer flat disk. A good match can be achieved by using an intermediate disk model with strongly reduced but non-zero flaring. Grain growth and dust settling processes provide a natural explanation for this disk geometry and we argue that such intermediate flaring might explain the observations of several other brown dwarf disks as well.
[145]  [pdf] - 65617
Evolution of Young Brown Dwarf Disks in the Mid-Infrared
Comments: 6 pages, 4 figures
Submitted: 2004-06-21, last modified: 2004-08-30
We have imaged two bona-fide brown dwarfs with TReCS/GEMINI-S and find mid-infrared excess emission that can be explained by optically thick dust disk models. In the case of the young ($\approx$2Myr) Cha H$\alpha$1 we measure fluxes at 10.4$\mu$m and 12.3$\mu$m that are fully consistent with a standard flared disk model and prominent silicate emission. For the $\approx$ 10Myr old brown dwarf 2MASS1207-3932 located in the TW Hydrae association we find excess emission at 8.7$\mu$m and 10.4$\mu$m with respect to its photosphere, and confirm disk accretion as likely cause of its strong activity. Disks around brown dwarfs likely last at least as long as their low-mass stellar counterparts in the T-Tauri phase. Grain growth, dust settling, and evolution of the geometry of brown dwarfs disks may appear on a timescale of 10Myr and can be witnessed by observations in the mid-infrared.
[146]  [pdf] - 65947
The hot core-ultracompact HII connection in G10.47+0.03
Comments: 12 pages, 5 figures, Accepted by Astronomy and Astrophysics
Submitted: 2004-07-06
We present infrared imaging and spectroscopic data of the complex massive star-forming region G10.47+0.03. The detection of seven mid-infrared (MIR) sources in our field combined with a sensitive Ks/ISAAC image allows to establish a very accurate astrometry, at the level of 0.3". Two MIR sources are found to be coincident with two ultracompact HII regions (UCHIIs) within our astrometric accuracy. Another MIR source lies very close to three other UCHII regions and to the hot molecular core (HMC) in G10.47+0.03. Spectroscopy of two of the most interesting MIR sources allows to identify the location and spectral type of the ionizing sources. We discuss in detail the relationship between the HMC, the UCHII regions and the nearby MIR source. The nature of the other MIR sources is also investigated.
[147]  [pdf] - 61186
FU Orionis: A Binary Star?
Comments: 10 pages, 2 figures, to be published in ApJ Letters
Submitted: 2003-11-27
By using the ALFA adaptive optics system at the 3.6m telescope of the Calar Alto Observatory we detected a faint red star in the apparent vicinity of FU Ori, the prototype of the FUor outburst stars. Independent confirmation of the detection is obtained from archival PUEO/CFHT images. The separation between the companion candidate and FU Ori is 0."50 and their brightness contrast is around 4 magnitudes. We discuss the possible nature of the newly detected star based on near-infrared photometry and its proper motion relative to FU Ori. The photometric data are consistent with a nearby late-type main sequence star, a background giant star, and a pre-main sequence star. On the basis of the proper motion and the stellar surface density in the direction towards FU Ori, we argue that the probabilities of the first two options are very low.
[148]  [pdf] - 60774
NACO Polarimetric Differential Imaging of TW Hya: A Sharp Look at the Closest T Tauri Disk
Comments: Accepted for publication in Astronomy & Astrophysics; 7 pages, 7 figures
Submitted: 2003-11-08
We present high-contrast imaging data on the disk of the classical T Tauri star TW Hya. The images were obtained through the polarimetric differential imaging technique with the adaptive optics system NACO. Our commissioning data show the presence of polarized disk emission between 0.1" and 1.4" from the star. We derive the first Ks-band radial polarized intensity distribution. We show that the polarized intensity compares well to shorter wavelengths surface brightness observations and confirm the previously reported gradual slope change around 0.8". These results show the potential of the new polarimetric differential imaging technique at 8m-class telescopes to map the inner regions of protoplanetary disks.
[149]  [pdf] - 57789
First Detection of Millimeter Dust Emission from Brown Dwarf Disks
Comments: 5 pages, accepted for ApJL
Submitted: 2003-07-03
We report results from the first deep millimeter continuum survey targeting Brown Dwarfs (BDs). The survey led to the first detection of cold dust in the disks around two young BDs (CFHT-BD-Tau 4 and IC348 613), with deep JCMT and IRAM observations reaching flux levels of a few mJy. The dust masses are estimated to be a few Earth masses assuming the same dust opacities as usually applied to TTauri stars.
[150]  [pdf] - 56668
The First Detailed Look at a Brown Dwarf Disk
Comments: 11 pages, 1 figure
Submitted: 2003-05-12
The combination of mid-infrared and recent submm/mm measurements allows us to set up the first comprehensive spectral energy distribution (SED) of the circumstellar material around a young Brown Dwarf. Simple arguments suggest that the dust is distributed in the form of a disk. We compare basic models to explore the disk parameters. The modeling shows that a flat disk geometry fits well the observations. A flared disk explains the SED only if it has a puffed-up inner rim and an inner gap much larger than the dust sublimation radius. Similarities and differences with disks around T Tauri stars are discussed.
[151]  [pdf] - 49841
Probing Dust around Brown Dwarfs: The naked LP 944-20 and the Disk of Cha Ha2
Submitted: 2002-06-13
We present the first mid-infrared (MIR) detection of a field brown dwarf (BD) and the first ground-based MIR measurements of a disk around a young BD candidate. We prove the absence of warm dust surrounding the field BD LP 944-20. In the case of the young BD candidate Cha Ha2, we find clear evidence for thermal dust emission from a disk. Surprisingly, the object does not exhibit any silicate feature as previously predicted. We show that the flat spectrum can be explained by an optically thick flat dust disk.