BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Eugene Chiang

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph.IM

Catalog footprint

What is connected

48works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

48 published item(s)

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) III: The vertical structure of debris disks

Debris disks -- collisionally sustained belts of dust and sometimes gas around main sequence stars -- are remnants of planet formation processes and are found in systems ${\gtrsim}10$ Myr old. Millimeter-wavelength observations are particularly important, as the grains probed by these observations are not strongly affected by radiation pressure and stellar winds, allowing them to probe the dynamics of large bodies producing dust. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is analyzing high-resolution observations of 24 debris disks to enable the characterization of debris disk substructures across a large sample for the first time. For the most highly inclined disks, it is possible to recover the vertical structure of the disk. We aim to model and analyze the most highly inclined systems in the ARKS sample in order to uniformly extract the vertical dust distributions for a sample of well-resolved debris disks. We employed both parametric and nonparametric methods to constrain the vertical dust distributions for the most highly inclined ARKS targets. We find a broad range of aspect ratios, revealing a wide diversity in vertical structure, with a range of best-fit parametric values of $0.0026 \leq h_{\rm HWHM} \leq 0.193$ and a median best-fit value of $h_{\rm HWHM}=0.021$. The results obtained by nonparametric modeling are generally consistent with the parametric modeling results. We find that five of the 13 disks are consistent with having total disk masses less than that of Neptune (17 $M_{\oplus}$), assuming stirring by internal processes (self-stirring and collisional and frictional damping). Furthermore, most systems show a significant preference for a Lorentzian vertical profile rather than a Gaussian.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IX: Gas-driven origin for the continuum arc in the debris disc of HD 121617

Debris discs were long considered to be largely gas-free environments governed by collisional fragmentation, gravitational stirring, and radiative forces. Recent CO detections show that gas is present, but its abundance and origin remain uncertain. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) revealed a narrow gas and dust ring in the disc HD 121617 with an asymmetric arc 40% brighter than the rest of the ring. We aim to constrain the total gas mass in HD 121617 assuming the dust arc is produced by hydrodynamical gas-dust interactions. We used the Dusty FARGO-ADSG code, modelling dust as Lagrangian particles, including radiation pressure and dust feedback, and varying the total gas mass. Simulations were compared to observations using radiative transfer. An unstable gas ring creates a size-dependent radial and azimuthal dust trap whose efficiency depends on gas mass. Two models, with 50 and 5 Earth masses of gas, reproduce both the ALMA band 7 arc and the outward offset of the VLT/SPHERE scattered-light ring via gas drag and radiation pressure. We infer a conservative gas-mass range of 2.5 to 250 Earth masses. If the ALMA asymmetry is caused by gas drag, the required gas mass compared with the observed CO implies substantial H2, consistent with primordial gas. HD 121617 would then be a hybrid disc between protoplanetary and debris stages. Since a planet could also create an arc, future observations are needed to distinguish these scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiwavelength Vertical Structure in the AU Mic Debris Disk: Characterizing the Collisional Cascade

Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $λ$ = 450 $μ$m dust continuum emission from the debris disk around the nearby M dwarf AU Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $λ= 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $μ$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AU Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $μ$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 °$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical Evidence of a Spiral Arm--Driving Planet in the MWC 758 Protoplanetary Disk

More than a dozen young stars host spiral arms in their surrounding protoplanetary disks. The excitation mechanisms of such arms are under debate. The two leading hypotheses -- companion-disk interaction and gravitational instability (GI) -- predict distinct motion for spirals. By imaging the MWC 758 spiral arm system at two epochs spanning ${\sim}5$ yr using the SPHERE instrument on the Very Large Telescope (VLT), we test the two hypotheses for the first time. We find that the pattern speeds of the spirals are not consistent with the GI origin. Our measurements further evince the existence of a faint "missing planet" driving the disk arms. The average spiral pattern speed is $0.\!^\circ22\pm0.\!^\circ03$ yr$^{-1}$, pointing to a driver at $172_{-14}^{+18}$ au around a $1.9$ $M_\odot$ central star if it is on a circular orbit. In addition, we witness time varying shadowing effects on a global scale that are likely originated from an inner disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Heavy-metal Jupiters by major mergers: metallicity vs. mass for giant planets

Some Jupiter-mass exoplanets contain $\sim$$100\, M_\oplus$ of metals, well above the $\sim$$10\, M_\oplus$ typically needed in a solid core to trigger giant planet formation by runaway gas accretion. We demonstrate that such `heavy-metal Jupiters' can result from planetary mergers near $\sim$10 au. Multiple cores accreting gas at runaway rates gravitationally perturb one another onto crossing orbits such that the average merger rate equals the gas accretion rate. Concurrent mergers and gas accretion implies the core mass scales with the total planet mass as $M_{\rm core} \propto M^{1/5}$ - heavier planets harbour heavier cores, in agreement with the observed mass-metallicity relation. While the average gas giant merges about once to double its core, others may merge multiple times, as merger trees grow chaotically. We show that the dispersion of outcomes inherent in mergers can reproduce the large scatter in observed planet metallicities, assuming $3-30\, M_\oplus$ pre-runaway cores. Mergers potentially correlate metallicity, eccentricity, and spin.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How Consumption and Repulsion Set Planetary Gap Depths and the Final Masses of Gas Giants

Planets open gaps in discs. Gap opening is typically modeled by considering the planetary Lindblad torque which repels disc gas away from the planet's orbit. But gaps also clear because the planet consumes local material. We present a simple, easy-to-use, analytic framework for calculating how gaps deplete and how the disc's structure as a whole changes by the combined action of Lindblad repulsion and planetary consumption. The final mass to which a gap-embedded gas giant grows is derived in tandem. The analytics are tested against 1D numerical experiments and calibrated using published multi-dimensional simulations. In viscous alpha discs, the planet, while clearing a gap, initially accretes practically all of the gas that tries to diffuse past, rapidly achieving super-Jupiter if not brown dwarf status. By contrast, in inviscid discs---that may still accrete onto their central stars by, say, magnetized winds---planets open deep, repulsion-dominated gaps. Then only a small fraction of the disc accretion flow is diverted onto the planet, which grows to a fraction of a Jupiter mass. Transitional disc cavities might be cleared by families of such low-mass objects opening inviscid, repulsion-dominated, overlapping gaps which allow most of the outer disc gas to flow unimpeded onto host stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Obliquity Constraints on an Extrasolar Planetary-Mass Companion

We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sub-Neptune Formation: The View from Resonant Planets

The orbital period ratios of neighbouring sub-Neptunes are distributed asymmetrically near first-order resonances. There are deficits of systems---"troughs" in the period ratio histogram---just short of commensurability, and excesses---"peaks"---just wide of it. We reproduce quantitatively the strongest peak-trough asymmetries, near the 3:2 and 2:1 resonances, using dissipative interactions between planets and their natal discs. Disc eccentricity damping captures bodies into resonance and clears the trough, and when combined with disc-driven convergent migration, draws planets initially wide of commensurability into the peak. The migration implied by the magnitude of the peak is modest; reductions in orbital period are $\sim$10\%, supporting the view that sub-Neptunes complete their formation more-or-less in situ. Once captured into resonance, sub-Neptunes of typical mass $\sim$$5$--$15 M_{\oplus}$ stay captured (contrary to an earlier claim), as they are immune to the overstability that afflicts lower mass planets. Driving the limited, short-scale migration is a gas disc depleted in mass relative to a solar-composition disc by 3--5 orders of magnitude. Such gas-poor but not gas-empty environments are quantitatively consistent with sub-Neptune core formation by giant impacts (and not, e.g., pebble accretion). While disc-planet interactions at the close of the planet formation era adequately explain the 3:2 and 2:1 asymmetries at periods $\gtrsim$ $5$--$15$ days, subsequent modification by stellar tides appears necessary at shorter periods, particularly for the 2:1.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The First Habitable Zone Earth-Sized Planet From TESS II: $Spitzer$ Confirms TOI-700 d

We present $Spitzer$ 4.5$μ$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Gemini Planet Imager view of the HD 32297 debris disk

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Circumplanetary Disk Dynamics in the Isothermal and Adiabatic Limits

Circumplanetary disks (CPDs) may be essential to the formation of planets, regulating their spin and accretion evolution. We perform a series of 3D hydrodynamics simulations in both the isothermal and adiabatic limits to systematically measure the rotation rates, sizes, and masses of CPDs as functions of $q_{\rm thermal}$, the ratio of the planet mass to the disk thermal mass. Our $q_{\rm thermal}$ ranges from 0.1 to 4; for our various disk temperatures, this corresponds to planet masses between 1 Earth mass and 4 Jupiter masses. Within this parameter space, we find that isothermal CPDs are disky and bound within $\sim$10\% of the planet's Bondi radius $r_{\rm B}$, with the innermost $\sim0.05\,r_{\rm B}$ in full rotational support. Adiabatic CPDs are spherical (therefore not actually "disks"), bound within $\sim0.2\,r_{\rm B}$, and mainly pressure-supported with rotation rates scaling linearly with $q_{\rm thermal}$; extrapolation suggests full rotational support of adiabatic envelopes at $\sim10\,q_{\rm thermal}$. Fast rotation and 3D supersonic flow render isothermal CPDs significantly different in structure from --- and orders of magnitude less massive than --- their 1D isothermal hydrostatic counterparts. Inside a minimum-mass solar nebula, even a maximally cooled, isothermal CPD around a 10 Earth-mass core may have less than 1 Earth mass, suggesting that gas giant formation may hinge on angular momentum transport processes in CPDs. Our CPD sizes and masses appear consistent with the regular satellites orbiting solar system giants.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Primer on Unifying Debris Disk Morphologies

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks. These interpretations suggest that the farthest reaches of planetary systems are perturbed by eccentric planets, possibly just a few Earth masses each.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bringing "The Moth" to Light: A Planet-Sculpting Scenario for the HD 61005 Debris Disk

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Correlations between compositions and orbits established by the giant impact era of planet formation

The giant impact phase of terrestrial planet formation establishes connections between super-Earths' orbital properties (semimajor axis spacings, eccentricities, mutual inclinations) and interior compositions (the presence or absence of gaseous envelopes). Using N-body simulations and analytic arguments, we show that spacings derive not only from eccentricities, but also from inclinations. Flatter systems attain tighter spacings, a consequence of an eccentricity equilibrium between gravitational scatterings, which increase eccentricities, and mergers, which damp them. Dynamical friction by residual disk gas plays a critical role in regulating mergers and in damping inclinations and eccentricities. Systems with moderate gas damping and high solid surface density spawn gas-enveloped super-Earths with tight spacings, small eccentricities, and small inclinations. Systems in which super-Earths coagulate without as much ambient gas, in disks with low solid surface density, produce rocky planets with wider spacings, larger eccentricities, and larger mutual inclinations. A combination of both populations can reproduce the observed distributions of spacings, period ratios, transiting planet multiplicities, and transit duration ratios exhibited by Kepler super-Earths. The two populations, both formed in situ, also help to explain observed trends of eccentricity vs. planet size, and bulk density vs. method of mass measurement (radial velocities vs. transit timing variations). Simplifications made in this study --- including the limited timespan of the simulations, and the approximate treatments of gas dynamical friction and gas depletion history --- should be improved upon in future work to enable a detailed quantitative comparison to the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dust dynamics in 2D gravito-turbulent disks

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} Ω> 1$, where $Ω$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} Ω\sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gap Opening in 3D: Single Planet Gaps

Giant planets can clear deep gaps when embedded in 2D (razor-thin) viscous circumstellar disks. We show by direct simulation that giant planets are just as capable of carving out gaps in 3D. Surface density maps are similar between 2D and 3D, even in detail. In particular, the scaling $Σ_{\rm gap} \propto q^{-2}$ of gap surface density with planet mass, derived from a global "zero-dimensional" balance of Lindblad and viscous torques, applies equally well to results obtained at higher dimensions. Our 3D simulations reveal extensive, near-sonic, meridional flows both inside and outside the gaps; these large-scale circulations might bear on disk compositional gradients, in dust or other chemical species. At high planet mass, gap edges are mildly Rayleigh unstable and intermittently shed streams of material into the gap - less so in 3D than in 2D.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How Spirals and Gaps Driven by Companions in Protoplanetary Disks Appear in Scattered Light at Arbitrary Viewing Angles

Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a disk by a companion varies with inclination and position angle in NIR scattered light. We compare the cases of a $3M_{\rm J}$ and a $0.1M_{\odot}$ companion, and make predictions suitable for testing with Gemini/GPI, VLT/NACO/SPHERE, and Subaru/HiCIAO/SCExAO. We find that the two trailing arms produced by an external perturber can have a variety of morphologies in inclined systems - they may appear as one trailing arm; two trailing arms on the same side of the disk; or two arms winding in opposite directions. The disk ring outside a planetary gap may also mimic spiral arms when viewed at high inclinations. We suggest potential explanations for the features observed in HH 30, HD 141569 A, AK Sco, HD 100546, and AB Aur. We emphasize that inclined views of companion-induced features cannot be converted into face-on views using simple and commonly practiced image deprojections.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Signatures of Gravitational Instability in Resolved Images of Protostellar Disks

Protostellar (class 0/I) disks, having masses comparable to those of their nascent host stars, and fed continuously from their natal infalling envelopes, are prone to gravitational instability (GI). Motivated by advances in near-infrared (NIR) adaptive optics imaging and mm-wave interferometry, we explore the observational signatures of GI in disks, using hydrodynamical and Monte Carlo radiative transfer simulations to synthesize NIR scattered light images and mm dust continuum maps. Spiral arms induced by GI, located at disk radii of hundreds of AUs, are local overdensities and have their photospheres displaced to higher altitudes above the disk midplane, arms therefore scatter more NIR light from their central stars than inter-arm regions, and are detectable at distances up to 1 kpc by Gemini/GPI, VLT/SPHERE, and Subaru/HiCIAO/SCExAO. By contrast, collapsed clumps formed by disk fragmentation have such strong local gravitational fields that their scattering photospheres are at lower altitudes, such fragments appear fainter than their surroundings in NIR scattered light. Spiral arms and streamers recently imaged in four FU Ori systems at NIR wavelengths resemble GI-induced structures and support the interpretation that FUors are gravitationally unstable protostellar disks. At mm wavelengths, both spirals and clumps appear brighter in thermal emission than the ambient disk and can be detected by ALMA at distances up to 0.4 kpc with one-hour integration times at ~0.1 arcsec resolution. Collapsed fragments having masses >=1 MJ can be detected by ALMA within ~10 minutes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two Transiting Low Density Sub-Saturns from K2

We report the discovery and confirmation of two sub-Saturn planets orbiting a bright (V = 11.3), metal-rich ([Fe/H] = 0.42 $\pm$ 0.04 dex) G3 dwarf in the K2 Campaign 2 field. The planets are 5.68 $\pm$ 0.56 Earth-radii and 7.82 $\pm$ 0.72 Earth-radii and have orbital periods of 20.8851 $\pm$ 0.0003 d and 42.3633$\pm$0.0006 d, near to the 2:1 mean-motion resonance. We obtained 32 radial velocities (RVs) with Keck/HIRES and detected the reflex motion due to EPIC-203771098b and c. These planets have masses of 21.0 $\pm$ 5.4 Earth-masses and 27.0 $\pm$ 6.9 Earth-masses, respectively. With low densities of 0.63 $\pm$ 0.25 g/cc and 0.31 $\pm$ 0.12 g/cc, respectively, the planets require thick envelopes of H/He to explain their large sizes and low masses. Interior structure models predict that the planets have fairly massive cores of 17.6 $\pm$ 4.3 Earth-masses and 16.1 $\pm$ 4.2 Earth-masses, respectively. They may have formed exterior to their present locations, accreted their H/He envelopes at large orbital distances, and migrated in as a resonant pair. The proximity to resonance, large transit depths, and host star brightness offer rich opportunities for TTV follow-up. Finally, the low surface gravities of the EPIC-203771098 planets make them favorable targets for transmission spectroscopy by HST, Spitzer, and JWST.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

$β$ Pictoris' inner disk in polarized light and new orbital parameters for $β$ Pictoris b

We present $H$-band observations of $β$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $β$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($ϕ=85.27°^{+0.26}_{-0.19}$) with a position angle $θ_{PA}=30.35°^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $θ_{PA}\approx29°$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $β$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $Ω=31.75°\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $β$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $β$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $β$ Pic b; 2) the mutual inclination of the inner disk and $β$ Pic $b$ is $4°$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $β$ Pic $b$ or self-stirring by the disk's parent bodies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Metallicity Recipe for Rocky Planets

Planets with sizes between those of Earth and Neptune divide into two populations: purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger gas-enveloped planets possessing voluminous and optically thick atmospheres. We show that whether a planet forms rocky or gas-enveloped depends on the solid surface density of its parent disk. Assembly times for rocky cores are sensitive to disk solid surface density. Lower surface densities spawn smaller planetary embryos; to assemble a core of given mass, smaller embryos require more mergers between bodies farther apart and therefore exponentially longer formation times. Gas accretion simulations yield a rule of thumb that a rocky core must be at least 2$M_\oplus$ before it can acquire a volumetrically significant atmosphere from its parent nebula. In disks of low solid surface density, cores of such mass appear only after the gas disk has dissipated, and so remain purely rocky. Higher surface density disks breed massive cores more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We test model predictions against observations, using planet radius as an observational proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface density. Theory can explain the observation that metal-rich stars host predominantly gas-enveloped planets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An M Dwarf Companion and Its Induced Spiral Arms in the HD 100453 Protoplanetary Disk

Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Breeding Super-Earths and Birthing Super-Puffs in Transitional Disks

The riddle posed by super-Earths (1-4$R_\oplus$, 2-20$M_\oplus$) is that they are not Jupiters: their core masses are large enough to trigger runaway gas accretion, yet somehow super-Earths accreted atmospheres that weigh only a few percent of their total mass. We show that this puzzle is solved if super-Earths formed late, as the last vestiges of their parent gas disks were about to clear. This scenario would seem to present fine-tuning problems, but we show that there are none. Ambient gas densities can span many (up to 9) orders of magnitude, and super-Earths can still robustly emerge after $\sim$0.1-1 Myr with percent-by-weight atmospheres. Super-Earth cores are naturally bred in gas-poor environments where gas dynamical friction has weakened sufficiently to allow constituent protocores to merge. So little gas is present at the time of core assembly that cores hardly migrate by disk torques: formation of super-Earths can be in situ. The picture --- that close-in super-Earths form in a gas-poor (but not gas-empty) inner disk, fed continuously by gas that bleeds inward from a more massive outer disk --- recalls the largely evacuated but still accreting inner cavities of transitional protoplanetary disks. We also address the inverse problem presented by super-puffs: an uncommon class of short-period planets seemingly too voluminous for their small masses (4-10$R_\oplus$, 2-6$M_\oplus$). Super-puffs easily acquire their thick atmospheres as dust-free, rapidly cooling worlds outside $\sim$1 AU where nebular gas is colder, less dense, and therefore less opaque. Unlike super-Earths which can form in situ, super-puffs migrated in to their current orbits; they are expected to form the outer links of mean-motion resonant chains, and to exhibit greater water content. We close by confronting observations and itemizing remaining questions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Eccentric Jupiters via Disk-Planet Interactions

Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude instability is consistent with that predicted analytically. Ellipticities can grow until they reach of order the disk's aspect ratio, beyond which the external Lindblad resonances that excite eccentricity are weakened by the planet's increasingly supersonic epicyclic motion. Forcing the planet to still larger eccentricities causes catastrophic eccentricity damping as the planet collides into gap walls. For standard parameters, the range of eccentricities for instability is modest; the threshold eccentricity for growth ($\sim$$0.04$) is not much smaller than the final eccentricity to which orbits grow ($\sim$$0.07$). If this threshold eccentricity can be lowered (perhaps by non-barotropic effects), and if the eccentricity driving documented here survives in 3D, it may robustly explain the low-to-moderate eccentricities $\lesssim 0.1$ exhibited by many giant planets (including Jupiter and Saturn), especially those without planetary or stellar companions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gemini Planet Imager Observations of the AU Microscopii Debris Disk: Asymmetries within One Arcsecond

We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband imaging polarimetry were obtained during the commissioning of GPI. In our broadband imaging polarimetry observations, we detect the disk only in total intensity and find asymmetries in the morphology of the disk between the southeast and northwest sides. The southeast side of the disk exhibits a bump at 1$''$ (10 AU projected separation) that is three times more vertically extended and three times fainter in peak surface brightness than the northwest side at similar separations. This part of the disk is also vertically offset by 69$\pm$30 mas to the northeast at 1$''$ when compared to the established disk mid-plane and consistent with prior ALMA and Hubble Space Telescope/STIS observations. We see hints that the southeast bump might be a result of detecting a horizontal sliver feature above the main disk that could be the disk backside. Alternatively when including the morphology of the northwest side, where the disk mid-plane is offset in the opposite direction $\sim$50 mas between 0$.''$4 and 1$.''$2, the asymmetries suggest a warp-like feature. Using our integral field spectroscopy data to search for planets, we are 50% complete for $\sim$4 $M_\mathrm{Jup}$ planets at 4 AU. We detect a source, resolved only along the disk plane, that could either be a candidate planetary mass companion or a compact clump in the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolved Millimeter-Wavelength Observations of Debris Disks around Solar-Type Stars

The presence of debris disks around young main sequence stars hints at the existence and structure of planetary systems. Millimeter-wavelength observations probe large grains that trace the location of planetesimal belts. The FEPS (Formation and Evolution of Planetary Systems) $Spitzer$ Legacy survey of nearby young solar analogues yielded a sample of five debris disk-hosting stars with millimeter flux suitable for interferometric follow-up. We present observations with the Submillimeter Array (SMA) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at ~2" resolution that spatially resolve the debris disks around these nearby ($d\sim$50 pc) stars. Two of the five disks (HD 377, HD 8907) are spatially resolved for the first time and one (HD 104860) is resolved at millimeter wavelengths for the first time. We combine our new observations with archival SMA and Atacama Large Millimeter/Submillimeter Array (ALMA) data to enable a uniform analysis of the full five-object sample. We simultaneously model the broad-band photometric data and resolved millimeter visibilities to constrain the dust temperatures and disk morphologies, and perform an MCMC analysis to fit for basic structural parameters. We find that the radii and widths of the cold outer belts exhibit properties consistent with scaled-up versions of the Solar System's Kuiper Belt. All the disks exhibit characteristic grain sizes comparable to the blowout size, and all the resolved observations of emission from large dust grains are consistent with an axisymmetric dust distribution to within the uncertainties. These results are consistent with comparable studies carried out at infrared wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spiral Arms in Gravitationally Unstable Protoplanetary Disks as Imaged in Scattered Light

Combining 3D smoothed-particle hydrodynamics and Monte Carlo radiative transfer calculations, we examine the morphology of spiral density waves induced by gravitational instability (GI) in protoplanetary disks, as they would appear in direct images at near-infrared (NIR) wavelengths. We find that systems with disk-to-star-mass ratios q=M_disk/M_star that are ~0.25 or more may produce prominent spiral arms in NIR imaging, remarkably resembling features observed in the MWC 758 and SAO 206462 systems. The contrast of GI-induced arms at NIR wavelengths can reach a factor of ~3, and their pitch angles are about 10-15 degree. The dominant azimuthal wavenumber of GI-induced spiral arms roughly obeys m~1/q in the range 2<~1/q<~8. In particular, a massive disk with q~0.5 can exhibit grand-design m=2 spirals. GI-induced arms are in approximate corotation with the local disk, and may therefore trap dust particles by pressure drag. Although GI can produce NIR spiral arms with morphologies, contrasts, and pitch angles similar to those reported in recent observations, it also makes other demands that may or may not be satisfied in any given system. A GI origin requires that the spirals be relatively compact, on scales <~100 AU; that the disk be massive, q>~0.25; and that the accretion rate Mdot be high, on the order of 1e-6 solar mass per year.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

To Cool is to Accrete: Analytic Scalings for Nebular Accretion of Planetary Atmospheres

Planets acquire atmospheres from their parent circumstellar disks. We derive a general analytic expression for how the atmospheric mass grows with time $t$, as a function of the underlying core mass $M_{\rm core}$ and nebular conditions, including the gas metallicity $Z$. Planets accrete as much gas as can cool: an atmosphere's doubling time is given by its Kelvin-Helmholtz time. Dusty atmospheres behave differently from atmospheres made dust-free by grain growth and sedimentation. The gas-to-core mass ratio (GCR) of a dusty atmosphere scales as GCR $\propto t^{0.4} M_{\rm core}^{1.7} Z^{-0.4} μ_{\rm rcb}^{3.4}$, where $μ_{\rm rcb} \propto 1/(1-Z)$ (for $Z$ not too close to 1) is the mean molecular weight at the innermost radiative-convective boundary. This scaling applies across all orbital distances and nebular conditions for dusty atmospheres; their radiative-convective boundaries, which regulate cooling, are not set by the external environment, but rather by the internal microphysics of dust sublimation, H$_2$ dissociation, and the formation of H$^-$. By contrast, dust-free atmospheres have their radiative boundaries at temperatures $T_{\rm rcb}$ close to nebular temperatures $T_{\rm out}$, and grow faster at larger orbital distances where cooler temperatures, and by extension lower opacities, prevail. At 0.1 AU in a gas-poor nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.9} M_{\rm core}^{1.6} Z^{-0.4} μ_{\rm rcb}^{3.3}$, while beyond 1 AU in a gas-rich nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.5} M_{\rm core}^1 Z^{-0.4}μ_{\rm rcb}^{2.2}$. We confirm our analytic scalings against detailed numerical models for objects ranging in mass from Mars (0.1 $M_\oplus$) to the most extreme super-Earths (10-20 $M_\oplus$), and explain why heating from planetesimal accretion cannot prevent the latter from undergoing runaway gas accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weak Turbulence in the HD 163296 Protoplanetary Disk Revealed by ALMA CO Observations

Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the turbulence in the disk around HD 163296, a nearby young A star, determined from ALMA Science Verification observations of four CO emission lines (the CO(3-2), CO(2-1), 13CO(2-1), and C18O(2-1) transitions). The different optical depths for these lines permit probes of non-thermal line-widths at a range of physical conditions (temperature and density) and depths into the disk interior. We derive stringent limits on the non-thermal motions in the upper layers of the outer disk such that any contribution to the line-widths from turbulence is <3% of the local sound speed. These limits are approximately an order of magnitude lower than theoretical predictions for full-blown MHD turbulence driven by the magneto-rotational instability, potentially suggesting that this mechanism is less efficient in the outer (R>30AU) disk than has been previously considered.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Class of Warm Jupiters with Mutually Inclined, Apsidally Misaligned, Close Friends

The orbits of giant extrasolar planets often have surprisingly small semi-major axes, large eccentricities, or severe misalignments between their normals and their host stars' spin axes. In some formation scenarios invoking Kozai-Lidov oscillations, an external planetary companion drives a planet onto an orbit having these properties. The mutual inclinations for Kozai-Lidov oscillations can be large and have not been confirmed observationally. Here we deduce that observed eccentric warm Jupiters with eccentric giant companions have mutual inclinations that oscillate between 35-65 deg. Our inference is based on the pairs' observed apsidal separations, which cluster near 90 deg. The near-orthogonality of periapse directions is effected by the outer companion's quadrupolar and octupolar potentials. These systems may be undergoing a stalled version of tidal migration that produces warm Jupiters over hot Jupiters, and provide evidence for a population of multi-planet systems that are not flat and have been sculpted by Kozai-Lidov oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast Modes and Dusty Horseshoes in Transitional Disks

The brightest transitional protoplanetary disks are often azimuthally asymmetric: their mm-wave thermal emission peaks strongly on one side. Dust overdensities can exceed $\sim$100:1, while gas densities vary by factors less than a few. We propose that these remarkable ALMA observations---which may bear on how planetesimals form---reflect a gravitational global mode in the gas disk. The mode is (1) fast---its pattern speed equals the disk's mean Keplerian frequency; (2) of azimuthal wavenumber $m=1$, displacing the host star from the barycenter; and (3) Toomre-stable. We solve for gas streamlines including the indirect stellar potential in the frame rotating with the pattern speed, under the drastic simplification that gas does not feel its own gravity. Near co-rotation, the gas disk takes the form of a horseshoe-shaped annulus. Dust particles with aerodynamic stopping times much shorter or much longer than the orbital period are dragged by gas toward the horseshoe center. For intermediate stopping times, dust converges toward a $\sim$45$^\circ$-wide arc on the co-rotation circle. Particles that do not reach their final accumulation points within disk lifetimes, either because of gas turbulence or long particle drift times, conform to horseshoe-shaped gas streamlines. Our mode is not self-consistent because we neglect gas self-gravity; still, we expect that trends between accumulation location and particle size, similar to those we have found, are generically predicted by fast modes and are potentially observable. Unlike vortices, global modes are not restricted in radial width to the pressure scale height; their large radial and azimuthal extents may better match observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravito-Turbulent Disks in 3D: Turbulent Velocities vs. Depth

Characterizing turbulence in protoplanetary disks is crucial for understanding how they accrete and spawn planets. Recent measurements of spectral line broadening promise to diagnose turbulence, with different lines probing different depths. We use 3D local hydrodynamic simulations of cooling, self-gravitating disks to resolve how motions driven by "gravito-turbulence" vary with height. We find that gravito-turbulence is practically as vigorous at altitude as at depth: even though gas at altitude is much too rarefied to be itself self-gravitating, it is strongly forced by self-gravitating overdensities at the midplane. The long-range nature of gravity means that turbulent velocities are nearly uniform vertically, increasing by just a factor of 2 from midplane to surface, even as the density ranges over nearly three orders of magnitude. The insensitivity of gravito-turbulence to height contrasts with the behavior of disks afflicted by the magnetorotational instability (MRI); in the latter case, non-circular velocities increase by at least a factor of 15 from midplane to surface, with various non-ideal effects only magnifying this factor. The distinct vertical profiles of gravito-turbulence vs. MRI turbulence may be used in conjunction with measurements of non-thermal linewidths at various depths to identify the source of transport in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Make Super-Earths, Not Jupiters: Accreting Nebular Gas onto Solid Cores at 0.1 AU and Beyond

Close-in super-Earths having radii 1--4 $R_\oplus$ may possess hydrogen atmospheres comprising a few percent by mass of their rocky cores. We determine the conditions under which such atmospheres can be accreted by cores from their parent circumstellar disks. Accretion from the nebula is problematic because it is too efficient: we find that 10-$M_\oplus$ cores embedded in solar metallicity disks tend to undergo runaway gas accretion and explode into Jupiters, irrespective of orbital location. The threat of runaway is especially dire at $\sim$0.1 AU, where solids may coagulate on timescales orders of magnitude shorter than gas clearing times; thus nascent atmospheres on close-in orbits are unlikely to be supported against collapse by planetesimal accretion. The time to runaway accretion is well approximated by the cooling time of the atmosphere's innermost convective zone, whose extent is controlled by where H$_2$ dissociates. Insofar as the temperatures characterizing H$_2$ dissociation are universal, timescales for core instability tend not to vary with orbital distance --- and to be alarmingly short for 10-$M_\oplus$ cores. Nevertheless, in the thicket of parameter space, we identify two scenarios, not mutually exclusive, that can reproduce the preponderance of percent-by-mass atmospheres for super-Earths at $\sim$0.1 AU, while still ensuring the formation of Jupiters at $\gtrsim 1$ AU. Scenario (a): planets form in disks with dust-to-gas ratios that range from $\sim$20$\times$ solar at 0.1 AU to $\sim$2$\times$ solar at 5 AU. Scenario (b): the final assembly of super-Earth cores from mergers of proto-cores --- a process that completes quickly at $\sim$0.1 AU once begun --- is delayed by gas dynamical friction until just before disk gas dissipates completely.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiwavelength Observations of the Candidate Disintegrating sub-Mercury KIC 12557548b

We present multiwavelength photometry, high angular resolution imaging, and radial velocities, of the unique and confounding disintegrating low-mass planet candidate KIC 12557548b. Our high angular resolution imaging, which includes spacebased HST/WFC3 observations in the optical, and groundbased Keck/NIRC2 observations in K'-band, allow us to rule-out background and foreground candidates at angular separations greater than 0.2 arcsec that are bright enough to be responsible for the transits we associate with KIC 12557548. Our radial velocity limit from Keck/HIRES allows us to rule-out bound, low-mass stellar companions to KIC 12557548 on orbits less than 10 years, as well as placing an upper-limit on the mass of the candidate planet of 1.2 Jupiter masses; therefore, the combination of our radial velocities, high angular-resolution imaging, and photometry are able to rule-out most false positive interpretations of the transits. Our precise multiwavelength photometry includes two simultaneous detections of the transit of KIC 12557548b using CFHT/WIRCam at 2.15 microns and the Kepler space telescope at 0.6 microns, as well as simultaneous null-detections of the transit by Kepler and HST/WFC3 at 1.4 microns. Our simultaneous HST/WFC3 and Kepler null-detections, provide no evidence for radically different transit depths at these wavelengths. Our simultaneous CFHT/WIRCam detections in the near-infrared and with Kepler in the optical reveal very similar transit depths (the average ratio of the transit depths at ~2.15 microns compared to ~0.6 microns is: 1.02 +/- 0.20). This suggests that if the transits we observe are due to scattering from single-size particles streaming from the planet in a comet-like tail, then the particles must be ~0.5 microns in radius or larger, which would favor that KIC 12557548b is a sub-Mercury, rather than super-Mercury, mass planet.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Catastrophic Evaporation of Rocky Planets

Short-period exoplanets can have dayside surface temperatures surpassing 2000 K, hot enough to vaporize rock and drive a thermal wind. Small enough planets evaporate completely. We construct a radiative-hydrodynamic model of atmospheric escape from strongly irradiated, low-mass rocky planets, accounting for dust-gas energy exchange in the wind. Rocky planets with masses < 0.1 M_Earth (less than twice the mass of Mercury) and surface temperatures > 2000 K are found to disintegrate entirely in < 10 Gyr. When our model is applied to Kepler planet candidate KIC 12557548b --- which is believed to be a rocky body evaporating at a rate of dM/dt > 0.1 M_Earth/Gyr --- our model yields a present-day planet mass of < 0.02 M_Earth or less than about twice the mass of the Moon. Mass loss rates depend so strongly on planet mass that bodies can reside on close-in orbits for Gyrs with initial masses comparable to or less than that of Mercury, before entering a final short-lived phase of catastrophic mass loss (which KIC 12557548b has entered). Because this catastrophic stage lasts only up to a few percent of the planet's life, we estimate that for every object like KIC 12557548b, there should be 10--100 close-in quiescent progenitors with sub-day periods whose hard-surface transits may be detectable by Kepler --- if the progenitors are as large as their maximal, Mercury-like sizes (alternatively, the progenitors could be smaller and more numerous). According to our calculations, KIC 12557548b may have lost ~70% of its formation mass; today we may be observing its naked iron core.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast Radial Flows in Transition Disk Holes

Protoplanetary "transition" disks have large, mass-depleted central cavities, yet also deliver gas onto their host stars at rates comparable to disks without holes. The paradox of simultaneous transparency and accretion can be explained if gas flows inward at much higher radial speeds inside the cavity than outside the cavity, since surface density (and by extension optical depth) varies inversely with inflow velocity at fixed accretion rate. Radial speeds within the cavity might even have to approach free-fall values to explain the huge surface density contrasts inferred for transition disks. We identify observational diagnostics of fast radial inflow in channel maps made in optically thick spectral lines. Signatures include (1) twisted isophotes in maps made at low systemic velocities and (2) rotation of structures observed between maps made in high-velocity line wings. As a test case, we apply our new diagnostic tools to archival ALMA data on the transition disk HD 142527, and uncover evidence for free-fall radial velocities inside its cavity. Although the observed kinematics are also consistent with a disk warp, the radial inflow scenario is preferred because it predicts low surface densities that appear consistent with recent observations of optically thin CO isotopologues in this disk. How material in the disk cavity sheds its angular momentum wholesale to fall freely onto the star is an unsolved problem; gravitational torques exerted by giant planets or brown dwarfs are briefly discussed as a candidate mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How Empty are Disk Gaps Opened by Giant Planets?

Gap clearing by giant planets has been proposed to explain the optically thin cavities observed in many protoplanetary disks. How much material remains in the gap determines not only how detectable young planets are in their birth environments, but also how strong corotation torques are, which impacts how planets can survive fast orbital migration. We determine numerically how the average surface density inside the gap, sigma_gap, depends on planet-to-star mass ratio q, Shakura-Sunyaev viscosity parameter alpha, and disk height-to-radius aspect ratio h/r. Our results are derived from our new GPU-accelerated Lagrangian hydrodynamical code PEnGUIn, and are verified by independent simulations with ZEUS90. For Jupiter-like planets, we find sigma_gap \propto q^-2.2 alpha^1.4 (h/r)^6.6, and for near brown dwarf masses, sigma_gap \propto q^-1 alpha^1.3 (h/r)^6.1. Surface density contrasts inside and outside gaps can be as large as 10^4, even when the planet does not accrete. We derive a simple analytic scaling, sigma_gap \propto q^-2 alpha^1 (h/r)^5, that compares reasonably well to empirical results, especially at low Neptune-like masses, and use discrepancies to highlight areas for progress.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Care and feeding of frogs

"Propellers" are features in Saturn's A ring associated with moonlets that open partial gaps. They exhibit non-Keplerian motion (Tiscareno 2010); the longitude residuals of the best-observed propeller, "Blériot," appear consistent with a sinusoid of period ~4 years. Pan and Chiang (2010) proposed that propeller moonlets librate in "frog resonances" with co-orbiting ring material. By analogy with the restricted three-body problem, they treated the co-orbital material as stationary in the rotating frame and neglected non-co-orbital material. Here we use simple numerical experiments to extend the frog model, including feedback due to the gap's motion, and drag associated with the Lindblad disk torques that cause Type I migration. Because the moonlet creates the gap, we expect the gap centroid to track the moonlet, but only after a time delay t_diff, the time for a ring particle to travel from conjunction with the moonlet to the end of the gap. We find that frog librations can persist only if t_diff exceeds the frog libration period P_lib, and if damping from Lindblad torques balances driving from co-orbital torques. If t_diff << P_lib, then the libration amplitude damps to zero. In the case of Blériot, the frog resonance model can reproduce the observed libration period P_lib ~ 4 yr. However, our simple feedback prescription suggests that Blériot's t_diff ~ 0.01P_lib, which is inconsistent with the observed libration amplitude of 260 km. We urge more accurate treatments of feedback to test the assumptions of our toy models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Colliding Planetary and Stellar Winds: Charge Exchange and Transit Spectroscopy in Neutral Hydrogen

When transiting their host stars, hot Jupiters absorb about 10% of the light in the wings of the stellar Lyman-alpha emission line. The absorption occurs at wavelengths Doppler-shifted from line center by +/- 100 km/s - larger than the thermal speeds with which partially neutral, 10E4 K hydrogen escapes from hot Jupiter atmospheres. It has been proposed that the absorption arises from 10E6 K hydrogen from the host stellar wind, made momentarily neutral by charge exchange with planetary H I. The +/-100 km/s velocities would then be attributed to the typical velocity dispersions of protons in the stellar wind - as inferred from spacecraft measurements of the Solar wind. To test this proposal, we perform 2D hydrodynamic simulations of colliding hot Jupiter and stellar winds, augmented by a chemistry module to compute the amount of hot neutral hydrogen produced by charge exchange. The Kelvin-Helmholtz instability mixes the two winds; in the mixing layer, charge exchange reactions establish, within tens of seconds. Enough hot neutral hydrogen is generated to reproduce the transit observations, and the amount of absorption converges with both spatial resolution and time. Our calculations support the idea that charge transfer between colliding winds correctly explains Lyman-alpha transit observations - modulo the effects of magnetic fields, which we do not model but which may suppress mixing. Other neglected effects include, in order of decreasing importance, rotational forces related to orbital motion, gravity, and stellar radiation pressure; we discuss quantitatively the errors introduced by our approximations. How hot stellar hydrogen cools when it collides with cold planetary hydrogen is also considered; a more careful treatment of how the mixing layer thermally equilibrates might explain the recent detection of Balmer H-alpha absorption in transiting hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From Dust to Planetesimals: Criteria for Gravitational Instability of Small Particles in Gas

Dust particles sediment toward the midplanes of protoplanetary disks, forming dust-rich sublayers encased in gas. What densities must the particle sublayer attain before it can fragment by self-gravity? We describe various candidate threshold densities. One of these is the Roche density, which is that required for a strengthless satellite to resist tidal disruption by its primary. Another is the Toomre density, which is that required for de-stabilizing self-gravity to defeat the stabilizing influences of pressure and rotation. We show that for sublayers containing aerodynamically well-coupled dust, the Toomre density exceeds the Roche density by many (up to about 4) orders of magnitude. We present 3D shearing box simulations of self-gravitating, stratified, dust-gas mixtures to test which of the candidate thresholds is relevant for collapse. All our simulations indicate that the larger Toomre density is required for collapse. This result is sensible because sublayers are readily stabilized by pressure. Sound-crossing times for thin layers are easily shorter than free-fall times, and the effective sound speed in dust-gas suspensions decreases only weakly with the dust-to-gas ratio (as the inverse square root). Our findings assume that particles are small enough that their stopping times in gas are shorter than all other timescales. Relaxing this assumption may lower the threshold for gravitational collapse back down to the Roche criterion. In particular, if the particle stopping time becomes longer than the sound-crossing time, sublayers may lose pressure support and become gravitationally unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Millimeter Emission Structure in the first ALMA Image of the AU Mic Debris Disk

We present 1.3 millimeter ALMA Cycle 0 observations of the edge-on debris disk around the nearby, ~10 Myr-old, M-type star AU Mic. These observations obtain 0.6 arcsec (6 AU) resolution and reveal two distinct emission components: (1) the previously known dust belt that extends to a radius of 40 AU, and (2) a newly recognized central peak that remains unresolved. The cold dust belt of mass about 1 lunar mass is resolved in the radial direction with a rising emission profile that peaks sharply at the location of the outer edge of the "birth ring" of planetesimals hypothesized to explain the midplane scattered light gradients. No significant asymmetries are discerned in the structure or position of this dust belt. The central peak identified in the ALMA image is ~6 times brighter than the stellar photosphere, which indicates an additional emission process in the inner regions of the system. Emission from a stellar corona or activity may contribute, but the observations show no signs of temporal variations characteristic of radio-wave flares. We suggest that this central component may be dominated by dust emission from an inner planetesimal belt of mass about 0.01 lunar mass, consistent with a lack of emission shortward of 25 microns and a location <3 AU from the star. Future millimeter observations can test this assertion, as an inner dust belt should be readily separated from the central star at higher angular resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stochastic flights of propellers

Kilometer-sized moonlets in Saturn's A ring create S-shaped wakes called "propellers" in surrounding material. The Cassini spacecraft has tracked the motions of propellers for several years and finds that they deviate from Keplerian orbits having constant semimajor axes. The inferred orbital migration is known to switch sign. We show using a statistical test that the time series of orbital longitudes of the propeller Blériot is consistent with that of a time-integrated Gaussian random walk. That is, Blériot's observed migration pattern is consistent with being stochastic. We further show, using a combination of analytic estimates and collisional N-body simulations, that stochastic migration of the right magnitude to explain the Cassini observations can be driven by encounters with ring particles 10-20 m in radius. That the local ring mass is concentrated in decameter-sized particles is supported on independent grounds by occultation analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Surface Layer Accretion in Conventional and Transitional Disks Driven by Far-Ultraviolet Ionization

Whether protoplanetary disks accrete at observationally significant rates by the magnetorotational instability (MRI) depends on how well ionized they are. Disk surface layers ionized by stellar X-rays are susceptible to charge neutralization by small condensates, ranging from ~0.01-micron-sized grains to angstrom-sized polycyclic aromatic hydrocarbons (PAHs). Ion densities in X-ray-irradiated surfaces are so low that ambipolar diffusion weakens the MRI. Here we show that ionization by stellar far-ultraviolet (FUV) radiation enables full-blown MRI turbulence in disk surface layers. Far-UV ionization of atomic carbon and sulfur produces a plasma so dense that it is immune to ion recombination on grains and PAHs. The FUV-ionized layer, of thickness 0.01--0.1 g/cm^2, behaves in the ideal magnetohydrodynamic limit and can accrete at observationally significant rates at radii > 1--10 AU. Surface layer accretion driven by FUV ionization can reproduce the trend of increasing accretion rate with increasing hole size seen in transitional disks. At radii < 1--10 AU, FUV-ionized surface layers cannot sustain the accretion rates generated at larger distance, and unless turbulent mixing of plasma can thicken the MRI-active layer, an additional means of transport is needed. In the case of transitional disks, it could be provided by planets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Forming Planetesimals by Gravitational Instability: I. The Role of the Richardson Number in Triggering the Kelvin-Helmholtz Instability

Gravitational instability (GI) of a dust-rich layer at the midplane of a gaseous circumstellar disk is one proposed mechanism to form planetesimals, the building blocks of rocky planets and gas giant cores. Self-gravity competes against the Kelvin-Helmholtz instability (KHI): gradients in dust content drive a vertical shear which risks overturning the dusty subdisk and forestalling GI. To understand the conditions under which the disk can resist the KHI, we perform 3D simulations of stratified subdisks in the limit that dust particles are small and aerodynamically well coupled to gas. This limit screens out the streaming instability and isolates the KHI. Each subdisk is assumed to have a vertical density profile given by a spatially constant Richardson number Ri. We vary Ri and the midplane dust-to-gas ratio mu and find that the critical Richardson number dividing KH-unstable from KH-stable flows is not unique; rather Ri_crit grows nearly linearly with mu for mu=0.3-10. Only for disks of bulk solar metallicity is Ri_crit ~ 0.2, close to the classical value. Our results suggest that a dusty sublayer can gravitationally fragment and presumably spawn planetesimals if embedded within a solar metallicity gas disk ~4x more massive than the minimum-mass solar nebula; or a minimum-mass disk having ~3x solar metallicity; or some intermediate combination of these two possibilities. Gravitational instability seems possible without resorting to the streaming instability or to turbulent concentration of particles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Forming Planetesimals by Gravitational Instability: II. How Dust Settles to its Marginally Stable State

Dust at the midplane of a circumstellar disk can become gravitationally unstable and fragment into planetesimals if the local dust-to-gas density ratio mu is sufficiently high. We simulate how dust settles in passive disks and ask how high mu can become. We settle the dust using a 1D code and test for dynamical stability using a 3D shearing box code. This scheme allows us to explore the behavior of small particles having short but non-zero stopping times in gas: 0 < t_stop << the orbital period. The streaming instability is thereby filtered out. Dust settles until shearing instabilities in the edges of the dust layer threaten to overturn the entire layer. In this state of marginal stability, mu=2.9 for a disk whose bulk (height-integrated) metallicity is solar. For a disk whose bulk metallicity is 4x solar, mu reaches 26.4. These maximum values of mu, which depend on the background radial pressure gradient, are so large that gravitational instability of small particles is viable in disks whose bulk metallicities are just a few (<4) times solar. Earlier studies assumed that dust settles until the Richardson number Ri is spatially constant. Our simulations are free of this assumption but provide support for it within the dust layer's edges, with the proviso that Ri increases with bulk metallicity in the same way that we found in Paper I. Only modest enhancements in bulk metallicity are needed to spawn planetesimals directly from small particles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Surface Layer Accretion in Transitional and Conventional Disks: From Polycyclic Aromatic Hydrocarbons to Planets

Transitional T Tauri disks have optically thin holes with radii > 10 AU, yet accrete up to the median T Tauri rate. Multiple planets inside the hole can torque the gas to high radial speeds over large distances, reducing the local surface density while maintaining accretion. Thus multi-planet systems, together with reductions in disk opacity due to grain growth, can explain how holes can be simultaneously transparent and accreting. There remains the problem of how outer disk gas diffuses into the hole. Here it has been proposed that the magnetorotational instability (MRI) erodes disk surface layers ionized by stellar X-rays. In contrast to previous work, we find that the extent to which surface layers are MRI-active is limited not by ohmic dissipation but by ambipolar diffusion, the latter measured by Am: the number of times a neutral hydrogen molecule collides with ions in a dynamical time. Simulations by Hawley & Stone showed that Am ~ 100 is necessary for ions to drive MRI turbulence in neutral gas. We calculate that in X-ray-irradiated surface layers, Am typically varies from ~0.001 to 1, depending on the abundance of charge-adsorbing polycyclic aromatic hydrocarbons, whose properties we infer from Spitzer observations. We conclude that ionization of H_2 by X-rays and cosmic rays can sustain, at most, only weak MRI turbulence in surface layers 1 to 10 g/cm^2 thick, and that accretion rates in such layers are too small compared to observed accretion rates for the majority of disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The propeller and the frog

"Propellers" in planetary rings are disturbances in ring material excited by moonlets that open only partial gaps. We describe a new type of co-orbital resonance that can explain the observed non-Keplerian motions of propellers. The resonance is between the moonlet underlying the propeller, and co-orbiting ring particles downstream of the moonlet where the gap closes. The moonlet librates within the gap about an equilibrium point established by co-orbiting material and stabilized by the Coriolis force. In the limit of small libration amplitude, the libration period scales linearly with the gap azimuthal width and inversely as the square root of the co-orbital mass. The new resonance recalls but is distinct from conventional horseshoe and tadpole orbits; we call it the "frog" resonance, after the relevant term in equine hoof anatomy. For a ring surface density and gap geometry appropriate for the propeller Blériot in Saturn's A ring, our theory predicts a libration period of ~4 years, similar to the ~3.7 year period over which Blériot's orbital longitude is observed to vary. These librations should be subtracted from the longitude data before any inferences about moonlet migration are made.

0 citations0 reviewsNo code linkedOpen access