-
Cavitating bubbles in condensing gas as a means of forming clumps, chondrites, and planetesimals
Authors:
Eugene Chiang
Abstract:
Vaporized metal, silicates, and ices on the verge of re-condensing into solid or liquid particles appear in many contexts: behind shocks, in impact ejecta, and within the atmospheres and outflows of stars, disks, planets, and minor bodies. We speculate that a condensing gas might fragment, forming overdensities within relative voids, from a radiation-condensation instability. Seeded with small the…
▽ More
Vaporized metal, silicates, and ices on the verge of re-condensing into solid or liquid particles appear in many contexts: behind shocks, in impact ejecta, and within the atmospheres and outflows of stars, disks, planets, and minor bodies. We speculate that a condensing gas might fragment, forming overdensities within relative voids, from a radiation-condensation instability. Seeded with small thermal fluctuations, a condensible gas will exhibit spatial variations in the density of particle condensates. Regions of higher particle density may radiate more, cooling faster. Faster cooling leads to still more condensation, lowering the local pressure. Regions undergoing runaway condensation may collapse under the pressure of their less condensed surroundings. Particle condensates will compactify with collapsing regions, into overdense clumps or macroscopic solids (planetesimals). As a first step toward realizing this hypothetical instability, we calculate the evolution of a small volume of condensing silicate vapor -- a spherical test "bubble" embedded in a background medium whose pressure and radiation field are assumed fixed for simplicity. Such a bubble condenses and collapses upon radiating its latent heat to the background, assuming its energy loss is not stopped by background irradiation. Collapse speeds can range up to sonic, similar to cavitation in terrestrial settings. Adding a non-condensible gas like hydrogen to the bubble stalls the collapse. We discuss whether cavitation can provide a way for mm-sized chondrules and refractory solids to assemble into meteorite parent bodies, focusing on CB/CH chondrites whose constituents likely condensed from silicate/metal vapor released from the most energetic asteroid collisions.
△ Less
Submitted 4 September, 2024;
originally announced September 2024.
-
The Resonant Remains of Broken Chains from Major and Minor Mergers
Authors:
Rixin Li,
Eugene Chiang,
Nick Choksi,
Fei Dai
Abstract:
Observations with the Transiting Exoplanet Survey Satellite and Kepler have revealed that practically all close-in sub-Neptunes form in mean-motion resonant chains, most of which unravel on timescales of 100 Myr. Using a series of N-body integrations, we study how planetary collisions resulting from the destabilization of resonant chains produce the distribution of orbital periods observed among m…
▽ More
Observations with the Transiting Exoplanet Survey Satellite and Kepler have revealed that practically all close-in sub-Neptunes form in mean-motion resonant chains, most of which unravel on timescales of 100 Myr. Using a series of N-body integrations, we study how planetary collisions resulting from the destabilization of resonant chains produce the distribution of orbital periods observed among mature systems, focusing on the resonant fine structures that remain post-instability. In their natal chains, planets near first-order resonances have period ratios just wide of perfect commensurability, driven there by disk migration and eccentricity damping. Sufficiently large resonant libration amplitudes (of unknown origin) are needed to trigger instability. Ensuing collisions between planets ("major mergers") erode but do not completely eliminate resonant pairs; survivors which avoid mergers show up as narrow "peaks" just wide of commensurability in the histogram of neighboring-planet period ratios. Merger products exhibit a broad range of period ratios, with each resonant peak in the histogram spawning a continuum of ratios wide of the given resonance. These continua may fill in period ratios between relatively closely separated resonances such as the 5:4, 4:3, and 3:2, but may fail to bridge the relatively wide gap between the 3:2 and the 2:1. Thus a "trough" manifests just short of the 2:1 resonance (and only the 2:1 resonance), as observed. Major mergers are not perfect, and generate collisional debris which undergoes "minor mergers" with planets, in many cases further widening resonant pairs. With all this dynamical activity, free eccentricities of resonant pairs, and by extension the phases of their transit timing variations (TTVs), are readily excited. Because non-resonant planets are merger products, they are predicted to have higher masses than resonant planets.
△ Less
Submitted 19 August, 2024;
originally announced August 2024.
-
The Prevalence of Resonance Among Young, Close-in Planets
Authors:
Fei Dai,
Max Goldberg,
Konstantin Batygin,
Jennifer van Saders,
Eugene Chiang,
Nick Choksi,
Rixin Li,
Erik A. Petigura,
Gregory J. Gilbert,
Sarah C. Millholland,
Yuan-Zhe Dai,
Luke Bouma,
Lauren M. Weiss,
Joshua N. Winn
Abstract:
Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common…
▽ More
Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common in younger systems. This prediction can now be tested, thanks to recent discoveries of young planets, particularly those in stellar clusters, by NASA's {\it TESS} mission. We divided the known planetary systems into three age groups: young ($<$100-Myr-old), adolescent (0.1-1-Gyr-old), and mature ($>1$-Gyr-old). The fraction of neighboring planet pairs having period ratios within a few percent of a first-order commensurability (e.g.~4:3, 3:2, or 2:1) is 70$\pm$15\% for young pairs, 24$\pm$8\% for adolescent pairs, and 15$\pm$2\% for mature pairs. The fraction of systems with at least one nearly commensurable pair (either first or second-order) is 86$\pm13$\% among young systems, 38$\pm12$\% for adolescent systems, and 23$\pm3$\% for mature systems. First-order commensurabilities prevail across all age groups, with an admixture of second-order commensurabilities. Commensurabilities are more common in systems with high planet multiplicity and low mutual inclinations. Observed period ratios often deviate from perfect commensurability by $\sim$1\% even among young planets, too large to be explained by resonant repulsion with equilibrium eccentricity tides. We also find that super-Earths in the radius gap ($1.5-1.9R_\oplus$) are less likely to be near-resonant (11.9$\pm2.0\%$) compared to Earth-sized planets ($R_p<1R_\oplus$; 25.3$\pm4.4\%$) or mini-Neptunes ($1.9R_\oplus \leq R_p<2.5R_\oplus$; 14.4$\pm1.8\%$).
△ Less
Submitted 21 August, 2024; v1 submitted 10 June, 2024;
originally announced June 2024.
-
Spectral Energy Distributions of Disc-Embedded Accreting Protoplanets
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
Many dozens of protoplanetary discs show signatures of sculpting by planets. To help find these protoplanets by direct imaging, we compute their broadband spectral energy distributions, which overlap with the JWST (James Webb Space Telescope) and ALMA (Atacama Large Millimeter Array) passbands. We consider circumplanetary spherical envelopes, and combinations of circumplanetary discs and envelopes…
▽ More
Many dozens of protoplanetary discs show signatures of sculpting by planets. To help find these protoplanets by direct imaging, we compute their broadband spectral energy distributions, which overlap with the JWST (James Webb Space Telescope) and ALMA (Atacama Large Millimeter Array) passbands. We consider circumplanetary spherical envelopes, and combinations of circumplanetary discs and envelopes, heated by protoplanet accretion and external irradiation. Searches with JWST's NIRCam (Near-Infrared Camera) and the blue portion of MIRI (Mid-Infrared Instrument) are most promising since $\sim$300--1000 K protoplanets outshine their $\sim$20--50 K circumstellar environs at wavelengths of $\sim$2--20 $μ$m. Detection is easier if circumplanetary dust settles into discs (more likely for more massive planets) or is less abundant per unit mass gas (because of grain growth or filtration). At wavelengths longer than 20 $μ$m, circumplanetary material is difficult to see against the circumstellar disc's surface layers that directly absorb starlight and reprocess it to the far-infrared. Such contaminating circumstellar emission can be serious even within the evacuated gaps observed by ALMA. Only in strongly depleted regions, like the cavity of the transitional disc PDS 70 where two protoplanets have been confirmed, may long-wavelength windows open for protoplanet study. We compile a list of candidate protoplanets and identify those with potentially the highest accretion luminosities, all peaking in the near-infrared.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Damping Obliquities of Hot Jupiter Hosts by Resonance Locking
Authors:
J. J. Zanazzi,
Janosz Dewberry,
Eugene Chiang
Abstract:
When orbiting hotter stars, hot Jupiters are often highly inclined relative to their host star equator planes. By contrast, hot Jupiters orbiting cooler stars are more aligned. Prior attempts to explain this correlation between stellar obliquity and effective temperature have proven problematic. We show how resonance locking -- the coupling of the planet's orbit to a stellar gravity mode (g mode)…
▽ More
When orbiting hotter stars, hot Jupiters are often highly inclined relative to their host star equator planes. By contrast, hot Jupiters orbiting cooler stars are more aligned. Prior attempts to explain this correlation between stellar obliquity and effective temperature have proven problematic. We show how resonance locking -- the coupling of the planet's orbit to a stellar gravity mode (g mode) -- can solve this mystery. Cooler stars with their radiative cores are more likely to be found with g-mode frequencies increased substantially by core hydrogen burning. Strong frequency evolution in resonance lock drives strong tidal evolution; locking to an axisymmetric g mode damps semi-major axes, eccentricities, and as we show for the first time, obliquities. Around cooler stars, hot Jupiters evolve into spin-orbit alignment and may avoid engulfment. Hotter stars lack radiative cores, and therefore preserve congenital spin-orbit misalignments. We focus on resonance locks with axisymmetric modes, supplementing our technical results with simple physical interpretations, and show that non-axisymmetric modes also damp obliquity. Outstanding issues regarding the dissipation of tidally-excited modes and the disabling of resonance locks are discussed quantitatively.
△ Less
Submitted 2 May, 2024; v1 submitted 8 March, 2024;
originally announced March 2024.
-
Sweeping Secular Resonances and Giant Planet Inclinations in Transition Discs
Authors:
J. J. Zanazzi,
Eugene Chiang
Abstract:
The orbits of some warm Jupiters are highly inclined (20$^\circ$-50$^\circ$) to those of their exterior companions. Comparable misalignments are inferred between the outer and inner portions of some transition discs. These large inclinations may originate from planet-planet and planet-disc secular resonances that sweep across interplanetary space as parent discs disperse. The maximum factor by whi…
▽ More
The orbits of some warm Jupiters are highly inclined (20$^\circ$-50$^\circ$) to those of their exterior companions. Comparable misalignments are inferred between the outer and inner portions of some transition discs. These large inclinations may originate from planet-planet and planet-disc secular resonances that sweep across interplanetary space as parent discs disperse. The maximum factor by which a seed mutual inclination can be amplified is of order the square root of the angular momentum ratio of the resonant pair. We identify those giant planet systems (e.g. Kepler-448 and Kepler-693) which may have crossed a secular resonance, and estimate the required planet masses and semimajor axes in transition discs needed to warp their innermost portions (e.g. in CQ Tau). Passage through an inclination secular resonance could also explain the hypothesized large mutual inclinations in apsidally-orthogonal warm Jupiter systems (e.g. HD 147018).
△ Less
Submitted 3 October, 2023; v1 submitted 20 July, 2023;
originally announced July 2023.
-
The maximum accretion rate of a protoplanet: how fast can runaway be?
Authors:
Nick Choksi,
Eugene Chiang,
Jeffrey Fung,
Zhaohuan Zhu
Abstract:
The hunt is on for dozens of protoplanets hypothesised to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations we calculate maximum gas accretion rates for planet masses $M_{\rm p}$ from 1$\,M_{\oplus}$ to $10\,M_{\rm J}$. When the planet is…
▽ More
The hunt is on for dozens of protoplanets hypothesised to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations we calculate maximum gas accretion rates for planet masses $M_{\rm p}$ from 1$\,M_{\oplus}$ to $10\,M_{\rm J}$. When the planet is small enough that its sphere of influence is fully embedded in the disc, with a Bondi radius $r_{\rm Bondi}$ smaller than the disc's scale height $H_{\rm p}$ -- such planets have thermal mass parameters $q_{\rm th} \equiv (M_{\rm p}/M_{\star}) / (H_{\rm p}/R_{\rm p})^3 \lesssim 0.3$, for host stellar mass $M_{\star}$ and orbital radius $R_{\rm p}$ -- the maximum accretion rate follows a Bondi scaling, with $\max \dot{M}_{\rm p} \propto ρ_{\rm g} M_{\rm p}^2 / (H_{\rm p}/R_{\rm p})^3$ for ambient disc density $ρ_{\rm g}$. For more massive planets with $0.3 \lesssim q_{\rm th} \lesssim 10$, the Hill sphere replaces the Bondi sphere as the gravitational sphere of influence, and $\max \dot{M}_{\rm p} \propto ρ_{\rm g} M_{\rm p}^1$, with no dependence on $H_{\rm p}/R_{\rm p}$. In the strongly superthermal limit when $q_{\rm th} \gtrsim 10$, the Hill sphere pops well out of the disc, and $\max \dot{M}_{\rm p} \propto ρ_{\rm g} M_{\rm p}^{2/3} (H_{\rm p}/R_{\rm p})^1$. Applied to the two confirmed protoplanets PDS 70b and c, our numerically calibrated maximum accretion rates imply their Jupiter-like masses may increase by up to a factor of $\sim$2 before their parent disc dissipates.
△ Less
Submitted 25 July, 2023; v1 submitted 2 May, 2023;
originally announced May 2023.
-
Giant Impacts and Debris Disk Morphology
Authors:
Joshua W. Jones,
Eugene Chiang,
Gaspard Duchene,
Paul Kalas,
Thomas M. Esposito
Abstract:
Certain debris disks have non-axisymmetric shapes in scattered light which are unexplained. The appearance of a disk depends on how its constituent Keplerian ellipses are arranged. The more the ellipses align apsidally, the more non-axisymmetric the disk. Apsidal alignment is automatic for fragments released from a catastrophic collision between solid bodies. We synthesize scattered light images,…
▽ More
Certain debris disks have non-axisymmetric shapes in scattered light which are unexplained. The appearance of a disk depends on how its constituent Keplerian ellipses are arranged. The more the ellipses align apsidally, the more non-axisymmetric the disk. Apsidal alignment is automatic for fragments released from a catastrophic collision between solid bodies. We synthesize scattered light images, and thermal emission images, of such giant impact debris. Depending on the viewing geometry, and if and how the initial apsidal alignment is perturbed, the remains of a giant impact can appear in scattered light as a one-sided or two-sided "fork", a lopsided "needle", or a set of "double wings". Double wings are difficult to reproduce in other scenarios involving gravitational forcing or gas drag, which do not align orbits as well. We compare our images with observations and offer a scorecard assessing whether the scattered light asymmetries in HD 15115, HD 32297, HD 61005, HD 111520, HD 106906, beta Pic, and AU Mic are best explained by giant impacts, gravitational perturbations, or sculpting by the interstellar medium.
△ Less
Submitted 1 May, 2023; v1 submitted 17 March, 2023;
originally announced March 2023.
-
Chaotic winds from a dying world: a one-dimensional map for evolving atmospheres
Authors:
Joshua Bromley,
Eugene Chiang
Abstract:
Planets which are smaller than Mercury and heated to sublimation temperatures of $\sim$2000 K lose mass catastrophically in dusty evaporative winds. The winds are observed to gust and recede largely without pattern; transit depths from the Kepler mission vary randomly from orbit to orbit by up to a factor of 10 or more. We explain how chaotic outflows may arise by constructing a map for the wind m…
▽ More
Planets which are smaller than Mercury and heated to sublimation temperatures of $\sim$2000 K lose mass catastrophically in dusty evaporative winds. The winds are observed to gust and recede largely without pattern; transit depths from the Kepler mission vary randomly from orbit to orbit by up to a factor of 10 or more. We explain how chaotic outflows may arise by constructing a map for the wind mass-loss rate as a function of time. The map is built on three statements: (1) The wind mass-loss rate scales in proportion to the surface equilibrium vapor pressure, rising exponentially with ground temperature. (2) Because the wind takes a finite time to escape the planet's gravity well, the surface mass-loss rate at any time determines the wind optical depth at a later time -- the atmosphere has hysteresis. (3) The ground temperature increases with optical depth (greenhouse effect) when the atmosphere is optically thin, and decreases with optical depth when the atmosphere is optically thick (nuclear winter). Statement (3) follows from how dust condenses in the face of intense stellar irradiation. As discussed recently, condensates initially naked before the star must be silicate-rich and iron-poor, staying cool enough for condensation by absorbing weakly in the visible and emitting strongly in the infrared. Later, when grains are numerous enough to self-shield from starlight, they may accrete more iron and reverse their visible-to-infrared opacity ratio. Depending on parameters, the map for the wind can regularly boom and bust between a greenhouse and a nuclear winter, or erupt into chaos. Lyapunov times are measured in orbital periods, the time for the wind to turn by Coriolis forces away from the planet's dayside, out of the Hill sphere.
△ Less
Submitted 27 March, 2023; v1 submitted 9 February, 2023;
originally announced February 2023.
-
Testing the Interaction Between a Substellar Companion and a Debris Disk in the HR 2562 System
Authors:
Stella Yimiao Zhang,
Gaspard Duchêne,
Robert J. De Rosa,
Megan Ansdell,
Quinn Konopacky,
Thomas Esposito,
Eugene Chiang,
Malena Rice,
Brenda Matthews,
Paul Kalas,
Bruce Macintosh,
Franck Marchis,
Stan Metchev,
Jenny Patience,
Julien Rameau,
Kimberly Ward-Duong,
Schuyler Wolff,
Michael P. Fitzgerald,
Vanessa P. Bailey,
Travis S. Barman,
Joanna Bulger,
Christine H. Chen,
Jeffrey K. Chilcotte,
Tara Cotten,
René Doyon
, et al. (29 additional authors not shown)
Abstract:
The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the…
▽ More
The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the orbital fit and, in combination with absolute astrometry from GAIA, place a 3$σ$ upper limit of 18.5 $M_J$ on the companion's mass. To interpret the ALMA observations, we used radiative transfer modeling to determine the disk properties. We find that the disk is well resolved and nearly edge on. While the misalignment angle between the disk and the orbit is weakly constrained due to the short orbital arc available, the data strongly support a (near) coplanar geometry for the system. Furthermore, we find that the models that describe the ALMA data best have an inner radius that is close to the companion's semi-major axis. Including a posteriori knowledge of the system's SED further narrows the constraints on the disk's inner radius and place it at a location that is in reasonable agreement with, possibly interior to, predictions from existing dynamical models of disk truncation by an interior substellar companion. HR\,2562 has the potential over the next few years to become a new testbed for dynamical interaction between a debris disk and a substellar companion.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
Exciting the TTV Phases of Resonant Sub-Neptunes
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
There are excesses of sub-Neptunes just wide of period commensurabilities like the 3:2 and 2:1, and corresponding deficits narrow of them. Any theory that explains this period ratio structure must also explain the strong transit timing variations (TTVs) observed near resonance. Besides an amplitude and a period, a sinusoidal TTV has a phase. Often overlooked, TTV phases are effectively integration…
▽ More
There are excesses of sub-Neptunes just wide of period commensurabilities like the 3:2 and 2:1, and corresponding deficits narrow of them. Any theory that explains this period ratio structure must also explain the strong transit timing variations (TTVs) observed near resonance. Besides an amplitude and a period, a sinusoidal TTV has a phase. Often overlooked, TTV phases are effectively integration constants, encoding information about initial conditions or the environment. Many TTVs near resonance exhibit non-zero phases. This observation is surprising because dissipative processes that capture planets into resonance also damp TTV phases to zero. We show how both the period ratio structure and the non-zero TTV phases can be reproduced if pairs of sub-Neptunes capture into resonance in a gas disc while accompanied by a third eccentric non-resonant body. Convergent migration and eccentricity damping by the disc drives pairs to orbital period ratios wide of commensurability; then, after the disc clears, secular forcing by the third body phase-shifts the TTVs. The scenario predicts that resonant planets are apsidally aligned and possess eccentricities up to an order of magnitude larger than previously thought.
△ Less
Submitted 19 March, 2023; v1 submitted 28 November, 2022;
originally announced November 2022.
-
Multiwavelength Vertical Structure in the AU Mic Debris Disk: Characterizing the Collisional Cascade
Authors:
David Vizgan,
A. Meredith Hughes,
Evan S. Carter,
Kevin M. Flaherty,
Margaret Pan,
Eugene Chiang,
Hilke Schlichting,
David J. Wilner,
Sean M. Andrews,
John M. Carpenter,
Atilla Moór,
Meredith A. MacGregor
Abstract:
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…
▽ More
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.
△ Less
Submitted 11 July, 2022;
originally announced July 2022.
-
A likely flyby of binary protostar Z CMa caught in action
Authors:
Ruobing Dong,
Hauyu Baobab Liu,
Nicolas Cuello,
Christophe Pinte,
Peter Abraham,
Eduard Vorobyov,
Jun Hashimoto,
Agnes Kospal,
Eugene Chiang,
Michihiro Takami,
Lei Chen,
Michael Dunham,
Misato Fukagawa,
Joel Green,
Yasuhiro Hasegawa,
Thomas Henning,
Yaroslav Pavlyuchenkov,
Tae-Soo Pyo,
Motohide Tamura
Abstract:
Close encounters between young stellar objects in star forming clusters are expected to dramatically perturb circumstellar disks. Such events are witnessed in numerical simulations of star formation, but few direct observations of ongoing encounters have been made. Here we report sub-0".1 resolution Atacama Large Millimeter Array (ALMA) and Jansky Very Large Array (JVLA) observations towards the m…
▽ More
Close encounters between young stellar objects in star forming clusters are expected to dramatically perturb circumstellar disks. Such events are witnessed in numerical simulations of star formation, but few direct observations of ongoing encounters have been made. Here we report sub-0".1 resolution Atacama Large Millimeter Array (ALMA) and Jansky Very Large Array (JVLA) observations towards the million year old binary protostar Z CMa in dust continuum and molecular line emission. A point source ~4700 au from the binary has been discovered at both millimeter and centimeter wavelengths. It is located along the extension of a ~2000 au streamer structure previously found in scattered light imaging, whose counterpart in dust and gas emission is also newly identified. Comparison with simulations shows signposts of a rare flyby event in action. Z CMa is a "double burster", as both binary components undergo accretion outbursts, which may be facilitated by perturbations to the host disk by flybys.
△ Less
Submitted 20 February, 2022; v1 submitted 14 January, 2022;
originally announced January 2022.
-
Mysterious Dust-emitting Object Orbiting TIC 400799224
Authors:
Brian P. Powell,
Veselin Kostov,
Saul Rappaport,
Andrei Tokovinin,
Avi Shporer,
Karen Collins,
Hank Corbett,
Tamas Borkovits,
Bruce Gary,
Eugene Chiang,
Joseph Rodriguez,
Nicholas Law,
Thomas Barclay,
Robert Gagliano,
Andrew Vanderburg,
Greg Olmschenk,
Ethan Kruse,
Joshua Schlieder,
Alan Soto,
Erin Goeke,
Thomas Jacobs,
Martti Kristiansen,
Daryll LaCourse,
Mark Omohundro,
Hans Schwengeler
, et al. (2 additional authors not shown)
Abstract:
We report the discovery of a unique object of uncertain nature -- but quite possibly a disintegrating asteroid or minor planet -- orbiting one star of the widely separated binary TIC 400799224. We initially identified the system in data from TESS Sector 10 via an abnormally-shaped fading event in the light curve (hereafter 'dips'). Follow-up speckle imaging determined that TIC 400799224 is actuall…
▽ More
We report the discovery of a unique object of uncertain nature -- but quite possibly a disintegrating asteroid or minor planet -- orbiting one star of the widely separated binary TIC 400799224. We initially identified the system in data from TESS Sector 10 via an abnormally-shaped fading event in the light curve (hereafter 'dips'). Follow-up speckle imaging determined that TIC 400799224 is actually two stars of similar brightness at 0.62" separation, forming a likely bound binary with projected separation of ~300 au. We cannot yet determine which star in the binary is host to the dips in flux. ASAS-SN and Evryscope archival data show that there is a strong periodicity of the dips at ~19.77 days, leading us to believe that an occulting object is orbiting the host star, though the duration, depth, and shape of the dips vary substantially. Statistical analysis of the ASAS-SN data shows that the dips only occur sporadically at a detectable threshold in approximately one out of every three to five transits, lending credence to the possibility that the occulter is a sporadically-emitted dust cloud. The cloud is also fairly optically thick, blocking up to 37% or 75% of the light from the host star, depending on the true host. Further observations may allow for greater detail to be gleaned as to the origin and composition of the occulter, as well as to a determination of which of the two stars comprising TIC 400799224 is the true host star of the dips.
△ Less
Submitted 3 October, 2021;
originally announced October 2021.
-
Testing planet formation from the ultraviolet to the millimeter
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
Gaps imaged in protoplanetary discs are suspected to be opened by planets. We compute the present-day mass accretion rates $\dot{M}_{\rm p}$ of seven hypothesized gap-embedded planets, plus the two confirmed planets in the PDS 70 disc. The accretion rates are based on disc gas surface densities $Σ_{\rm gas}$ from C$^{18}$O observations, and planet masses $M_{\rm p}$ from simulations fitted to obse…
▽ More
Gaps imaged in protoplanetary discs are suspected to be opened by planets. We compute the present-day mass accretion rates $\dot{M}_{\rm p}$ of seven hypothesized gap-embedded planets, plus the two confirmed planets in the PDS 70 disc. The accretion rates are based on disc gas surface densities $Σ_{\rm gas}$ from C$^{18}$O observations, and planet masses $M_{\rm p}$ from simulations fitted to observed gaps. Assuming accretion is Bondi-like, we find in eight out of nine cases that $\dot{M}_{\rm p}$ is consistent with the time-averaged value given by the current planet mass and system age, $M_{\rm p}/t_{\rm age}$. As system ages are comparable to circumstellar disc lifetimes, these gap-opening planets may be undergoing their last mass doublings, reaching final masses of $M_{\rm p} \sim 10-10^2 \, M_\oplus$ for the non-PDS 70 planets, and $M_{\rm p} \sim 1-10 \, M_{\rm J}$ for the PDS 70 planets. For another fifteen gaps without C$^{18}$O data, we predict $Σ_{\rm gas}$ by assuming their planets are accreting at their time-averaged $\dot{M}_{\rm p}$. Bondi accretion rates for PDS 70b and c are orders of magnitude higher than accretion rates implied by measured U-band and H$α$ fluxes, suggesting most of the accretion shock luminosity emerges in as yet unobserved wavebands, or that the planets are surrounded by dusty, highly extincting, quasi-spherical circumplanetary envelopes. Thermal emission from such envelopes or from circumplanetary discs, on Hill sphere scales, peaks at wavelengths in the mid-to-far-infrared and can reproduce observed mm-wave excesses.
△ Less
Submitted 29 November, 2021; v1 submitted 30 September, 2021;
originally announced October 2021.
-
Eccentric millisecond pulsars by resonant convection
Authors:
Sivan Ginzburg,
Eugene Chiang
Abstract:
Eccentric millisecond pulsars (eMSPs) with white dwarf companions exhibit orbital eccentricities orders of magnitude larger than predicted by turbulent convection in the white dwarfs' red giant progenitors. The orbital periods of eMSPs cluster around $P=20-30$ d, remarkably close to the red giant convective eddy turnover time $t_{\rm eddy}$. We propose that the anomalously large eccentricities are…
▽ More
Eccentric millisecond pulsars (eMSPs) with white dwarf companions exhibit orbital eccentricities orders of magnitude larger than predicted by turbulent convection in the white dwarfs' red giant progenitors. The orbital periods of eMSPs cluster around $P=20-30$ d, remarkably close to the red giant convective eddy turnover time $t_{\rm eddy}$. We propose that the anomalously large eccentricities are resonantly driven by convective flows somehow made coherent when the turnover time matches the tidally locked red giant's spin period, which is also the tidal forcing period. Numerical simulations of rotating red giants and magnetic field studies of stars show some evidence for especially ordered flow patterns when the convective Rossby number $P/t_{\rm eddy}$ is of order unity. We show that resonant convection boosts eccentricities by a factor of $(t_{\rm nuc}/P)^{1/2}\approx 3\times 10^3$ over the random-walk values that characterize conventional MSPs, in good agreement with observations ($t_{\rm nuc}$ is the giant's nuclear burning time-scale). We also show how variations in the eddy turnover time arising from red giant metallicity variations can reproduce the observed effective width of the resonance, $ΔP/P\approx 0.4$.
△ Less
Submitted 21 September, 2021;
originally announced September 2021.
-
Obliquity Constraints on the Planetary-mass Companion HD 106906 b
Authors:
Marta L. Bryan,
Eugene Chiang,
Caroline V. Morley,
Gregory N. Mace,
Brendan P. Bowler
Abstract:
We constrain the angular momentum architecture of HD 106906, a 13 $\pm$ 2 Myr old system in the ScoCen complex composed of a compact central binary, a widely separated planetary-mass tertiary HD 106906 b, and a debris disk nested between the binary and tertiary orbital planes. We measure the orientations of three vectors: the companion spin axis, companion orbit normal, and disk normal. Using near…
▽ More
We constrain the angular momentum architecture of HD 106906, a 13 $\pm$ 2 Myr old system in the ScoCen complex composed of a compact central binary, a widely separated planetary-mass tertiary HD 106906 b, and a debris disk nested between the binary and tertiary orbital planes. We measure the orientations of three vectors: the companion spin axis, companion orbit normal, and disk normal. Using near-IR high-resolution spectra from Gemini/IGRINS, we obtain a projected rotational velocity of $v\sin{i_p}$ = 9.5 $\pm$ 0.2 km/s for HD 106906 b. This measurement together with a published photometric rotation period implies the companion is viewed nearly pole-on, with a line-of-sight spin axis inclination of $i_p$ = 14 $\pm$ 4 degrees or 166 $\pm$ 4 degrees. By contrast, the debris disk is known to be viewed nearly edge-on. The likely misalignment of all three vectors suggests HD 106906 b formed by gravitational instability in a turbulent environment, either in a disk or cloud setting.
△ Less
Submitted 30 August, 2021;
originally announced August 2021.
-
Primordial Obliquities of Brown Dwarfs and Super-Jupiters from Fragmenting Gravito-Turbulent Discs
Authors:
R. Michael Jennings,
Eugene Chiang
Abstract:
Super-Jupiters, brown dwarfs, and stars can form from the collapse of self-gravitating discs. Such discs are turbulent, with flocculent spiral arms accelerating gas to transonic speeds horizontally and vertically. Objects that fragment from gravito-turbulent discs should spin with a wide range of directions, reflecting the random orientations of their parent eddies. We show by direct numerical sim…
▽ More
Super-Jupiters, brown dwarfs, and stars can form from the collapse of self-gravitating discs. Such discs are turbulent, with flocculent spiral arms accelerating gas to transonic speeds horizontally and vertically. Objects that fragment from gravito-turbulent discs should spin with a wide range of directions, reflecting the random orientations of their parent eddies. We show by direct numerical simulation that obliquities of newly collapsed fragments can range up to 45$^\circ$. Subsequent collisions between fragments can further alter the obliquity distribution, up to 90$^\circ$ or down to near-zero. The large obliquities of newly discovered super-Jupiters on wide orbits around young stars may be gravito-turbulent in origin. Obliquely spinning fragments are born on orbits that may be inclined relative to their parent discs by up to 20$^\circ$, and gravitationally stir leftover material to many times the pre-fragmentation disc thickness.
△ Less
Submitted 21 August, 2021; v1 submitted 25 May, 2021;
originally announced May 2021.
-
A coplanar circumbinary protoplanetary disk in the TWA 3 triple M dwarf system
Authors:
Ian Czekala,
Álvaro Ribas,
Nicolás Cuello,
Eugene Chiang,
Enrique Macías,
Gaspard Duchêne,
Sean M. Andrews,
Catherine C. Espaillat
Abstract:
We present sensitive ALMA observations of TWA 3, a nearby, young ($\sim$10 Myr) hierarchical system composed of three pre-main sequence M3--M4.5 stars. For the first time, we detected ${}^{12}$CO and ${}^{13}$CO $J$=2-1 emission from the circumbinary protoplanetary disk around TWA 3A. We jointly fit the protoplanetary disk velocity field, stellar astrometric positions, and stellar radial velocitie…
▽ More
We present sensitive ALMA observations of TWA 3, a nearby, young ($\sim$10 Myr) hierarchical system composed of three pre-main sequence M3--M4.5 stars. For the first time, we detected ${}^{12}$CO and ${}^{13}$CO $J$=2-1 emission from the circumbinary protoplanetary disk around TWA 3A. We jointly fit the protoplanetary disk velocity field, stellar astrometric positions, and stellar radial velocities to infer the architecture of the system. The Aa and Ab stars ($0.29\pm0.01\,M_\odot$ and $0.24\pm0.01\,M_\odot$, respectively) comprising the tight ($P=35$ days) eccentric ($e=0.63\pm0.01$) spectroscopic binary are coplanar with their circumbinary disk (misalignment $< 6^{\circ}$ with 68% confidence), similar to other short-period binary systems. From models of the spectral energy distribution, we found the inner radius of the circumbinary disk ($r_\mathrm{inner} = 0.50 - 0.75$ au) to be consistent with theoretical predictions of dynamical truncation $r_\mathrm{cav}/a_\mathrm{inner} \approx 3$. The outer orbit of the tertiary star B ($0.40\pm0.28\,M_\odot$, $a\sim65 \pm 18$ au, $e=0.3\pm0.2$) is not as well constrained as the inner orbit, however, orbits coplanar with the A system are still preferred (misalignment $ < 20^{\circ}$). To better understand the influence of the B orbit on the TWA 3A circumbinary disk, we performed SPH simulations of the system and found that the outer edge of the gas disk ($r_\mathrm{outer}=8.5\pm0.2$ au) is most consistent with truncation from a coplanar, circular or moderately eccentric orbit, supporting the preference from the joint orbital fit.
△ Less
Submitted 2 March, 2021; v1 submitted 23 February, 2021;
originally announced February 2021.
-
Resolving Structure in the Debris Disk around HD 206893 with ALMA
Authors:
Ava Nederlander,
A. Meredith Hughes,
Anna J. Fehr,
Kevin M. Flaherty,
Kate Y. L. Su,
Attila Moor,
Eugene Chiang,
Sean M. Andrews,
David J. Wilner,
Sebastian Marino
Abstract:
Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and waveleng…
▽ More
Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and wavelength of 1.3 mm with the Atacama Large Millimeter/submillimeter Array (ALMA). We observe a broad disk extending from a radius of <51 au to 194^{+13}_{-2} au. We model the disk with a continuous, gapped, and double power-law model of the surface density profile, and find strong evidence for a local minimum in the surface density distribution near a radius of 70 au, consistent with a gap in the disk with an inner radius of 63^{+8}_{-16} au and width 31^{+11}_{-7} au. Gapped structure has been observed in four other debris disks -- essentially every other radially resolved debris disk observed with sufficient angular resolution and sensitivity with ALMA -- and could be suggestive of the presence of an additional planetary-mass companion.
△ Less
Submitted 21 January, 2021;
originally announced January 2021.
-
As the Worlds Turn: Constraining Spin Evolution in the Planetary-Mass Regime
Authors:
Marta L. Bryan,
Sivan Ginzburg,
Eugene Chiang,
Caroline Morley,
Brendan P. Bowler,
Jerry W. Xuan,
Heather A. Knutson
Abstract:
To understand how planetary spin evolves and traces planet formation processes, we measure rotational line broadening in eight planetary-mass objects (PMOs) of various ages (1--800 Myr) using near-infrared high-resolution spectra from NIRSPEC/Keck. Combining these with published rotation rates, we compile 27 PMO spin velocities, 16 of which derive from our NIRSPEC/Keck program. Our data are consis…
▽ More
To understand how planetary spin evolves and traces planet formation processes, we measure rotational line broadening in eight planetary-mass objects (PMOs) of various ages (1--800 Myr) using near-infrared high-resolution spectra from NIRSPEC/Keck. Combining these with published rotation rates, we compile 27 PMO spin velocities, 16 of which derive from our NIRSPEC/Keck program. Our data are consistent with spin velocities $v$ scaling with planetary radius $R$ as $v \propto 1/R$. We conclude that spin angular momentum is conserved as objects cool and contract over the sampled age range. The PMOs in our sample spin at rates that are approximately an order of magnitude below their break-up values, consistent with the hypothesis that they were spun down by magnetized circum-PMO disks (CPDs) during the formation era at ages $\lesssim$ a few Myr. There is a factor of 4--5 variation in spin velocity that has yet to be understood theoretically. It also remains to be seen whether spin evolves on timescales $\gtrsim$ 1 Gyr for PMOs, as it does for stars and high-mass brown dwarfs emitting magnetized winds.
△ Less
Submitted 14 October, 2020;
originally announced October 2020.
-
An ALMA survey of $λ$ Orionis disks: from supernovae to planet formation
Authors:
Megan Ansdell,
Thomas J. Haworth,
Jonathan P. Williams,
Stefano Facchini,
Andrew Winter,
Carlo F. Manara,
Alvaro Hacar,
Eugene Chiang,
Sierk van Terwisga,
Nienke van der Marel,
Ewine F. van Dishoeck
Abstract:
Protoplanetary disk surveys by the Atacama Large Millimeter/sub-millimeter Array (ALMA) are now probing a range of environmental conditions, from low-mass star-forming regions like Lupus to massive OB clusters like $σ$ Orionis. Here we conduct an ALMA survey of protoplanetary disks in $λ$ Orionis, a ~5 Myr old OB cluster in Orion, with dust mass sensitivities comparable to the surveys of nearby re…
▽ More
Protoplanetary disk surveys by the Atacama Large Millimeter/sub-millimeter Array (ALMA) are now probing a range of environmental conditions, from low-mass star-forming regions like Lupus to massive OB clusters like $σ$ Orionis. Here we conduct an ALMA survey of protoplanetary disks in $λ$ Orionis, a ~5 Myr old OB cluster in Orion, with dust mass sensitivities comparable to the surveys of nearby regions (~0.4 $M_\oplus$). We assess how massive OB stars impact planet formation, in particular from the supernova that may have occurred ~1 Myr ago in the core of $λ$ Orionis; studying these effects is important as most planetary systems, including our Solar System, are likely born in cluster environments. We find that the effects of massive stars, in the form of pre-supernova feedback and/or a supernova itself, do not appear to significantly reduce the available planet-forming material otherwise expected at the evolved age of $λ$ Orionis. We also compare a lingering massive "outlier" disk in $λ$ Orionis to similar systems in other evolved regions, hypothesizing that these outliers host companions in their inner disks that suppress disk dispersal to extend the lifetimes of their outer primordial disks. We conclude with numerous avenues for future work, highlighting how $λ$ Orionis still has much to teach us about perhaps one of the most common types of planet-forming environments in the Galaxy.
△ Less
Submitted 30 September, 2020;
originally announced October 2020.
-
Chondrules from high-velocity collisions: thermal histories and the agglomeration problem
Authors:
Nick Choksi,
Eugene Chiang,
Harold C. Connolly Jr.,
Zack Gainsforth,
Andrew J. Westphal
Abstract:
We assess whether chondrules, once-molten mm-sized spheres filling the oldest meteorites, could have formed from super-km/s collisions between planetesimals in the solar nebula. High-velocity collisions release hot and dense clouds of silicate vapor which entrain and heat chondrule precursors. Thermal histories of CB chondrules are reproduced for colliding bodies $\sim$10--100 km in radius. The sl…
▽ More
We assess whether chondrules, once-molten mm-sized spheres filling the oldest meteorites, could have formed from super-km/s collisions between planetesimals in the solar nebula. High-velocity collisions release hot and dense clouds of silicate vapor which entrain and heat chondrule precursors. Thermal histories of CB chondrules are reproduced for colliding bodies $\sim$10--100 km in radius. The slower cooling rates of non-CB, porphyritic chondrules point to colliders with radii $\gtrsim$ 500 km. How chondrules, collisionally dispersed into the nebula, agglomerated into meteorite parent bodies remains a mystery. The same orbital eccentricities and inclinations that enable energetic collisions prevent planetesimals from re-accreting chondrules efficiently and without damage; thus the sedimentary laminations of the CB/CH chondrite Isheyevo are hard to explain by direct fallback of collisional ejecta. At the same time, planetesimal surfaces may be littered with the shattered remains of chondrules. The micron-sized igneous particles recovered from comet 81P/Wild-2 may have originated from in-situ collisions and subsequent accretion in the proto-Kuiper belt, obviating the need to transport igneous solids across the nebula. Asteroid sample returns from Hayabusa2 and OSIRIS-REx may similarly contain chondrule fragments.
△ Less
Submitted 28 March, 2021; v1 submitted 21 September, 2020;
originally announced September 2020.
-
Dynamical Evidence of a Spiral Arm--Driving Planet in the MWC 758 Protoplanetary Disk
Authors:
Bin Ren,
Ruobing Dong,
Rob G. van Holstein,
Jean-Baptiste Ruffio,
Benjamin A. Calvin,
Julien H. Girard,
Myriam Benisty,
Anthony Boccaletti,
Thomas M. Esposito,
Élodie Choquet,
Dimitri Mawet,
Laurent Pueyo,
Tomas Stolker,
Eugene Chiang,
Jozua de Boer,
John H. Debes,
Antonio Garufi,
Carol A. Grady,
Dean C. Hines,
Anne-Lise Maire,
François Ménard,
Maxwell Millar-Blanchaer,
Marshall D. Perrin,
Charles A. Poteet,
Glenn Schneider
Abstract:
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 Larg…
▽ More
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.
△ Less
Submitted 29 July, 2020; v1 submitted 9 July, 2020;
originally announced July 2020.
-
Heavy-metal Jupiters by major mergers: metallicity vs. mass for giant planets
Authors:
Sivan Ginzburg,
Eugene Chiang
Abstract:
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 orb…
▽ More
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.
△ Less
Submitted 22 June, 2020;
originally announced June 2020.
-
Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign
Authors:
Thomas M. Esposito,
Paul Kalas,
Michael P. Fitzgerald,
Maxwell A. Millar-Blanchaer,
Gaspard Duchene,
Jennifer Patience,
Justin Hom,
Marshall D. Perrin,
Robert J. De Rosa,
Eugene Chiang,
Ian Czekala,
Bruce Macintosh,
James R. Graham,
Megan Ansdell,
Pauline Arriaga,
Sebastian Bruzzone,
Joanna Bulger,
Christine H. Chen,
Tara Cotten,
Ruobing Dong,
Zachary H. Draper,
Katherine B. Follette,
Li-Wei Hung,
Ronald Lopez,
Brenda C. Matthews
, et al. (40 additional authors not shown)
Abstract:
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. Observ…
▽ More
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.
△ Less
Submitted 23 June, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
How Consumption and Repulsion Set Planetary Gap Depths and the Final Masses of Gas Giants
Authors:
Mickey Rosenthal,
Eugene Chiang,
Sivan Ginzburg,
Ruth Murray-Clay
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 June, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
The Gemini Planet Imager view of the HD 32297 debris disk
Authors:
Gaspard Duchene,
Malena Rice,
Justin Hom,
Joseph Zalesky,
Thomas M. Esposito,
Maxwell A. Millar-Blanchaer,
Bin Ren,
Paul Kalas,
Michael Fitzgerald,
Pauline Arriaga,
Sebastian Bruzzone,
Joanna Bulger,
Christine H. Chen,
Eugene Chiang,
Tara Cotten,
Ian Czekala,
Robert J. De Rosa,
Ruobing Dong,
Zachary H. Draper,
Katherine B. Follette,
James R. Graham,
Li-Wei Hung,
Ronald Lopez,
Bruce Macintosh,
Brenda C. Matthews
, et al. (38 additional authors not shown)
Abstract:
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 ana…
▽ More
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.
△ Less
Submitted 13 April, 2020;
originally announced April 2020.
-
Sub-Neptune Formation: The View from Resonant Planets
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
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 betwee…
▽ More
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.
△ Less
Submitted 20 May, 2020; v1 submitted 6 March, 2020;
originally announced March 2020.
-
Obliquity Constraints on an Extrasolar Planetary-Mass Companion
Authors:
Marta L. Bryan,
Eugene Chiang,
Brendan P. Bowler,
Caroline V. Morley,
Sarah Millholland,
Sarah Blunt,
Katelyn B. Ashok,
Eric Nielsen,
Henry Ngo,
Dimitri Mawet,
Heather A. Knutson
Abstract:
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 orbi…
▽ More
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.
△ Less
Submitted 27 May, 2020; v1 submitted 25 February, 2020;
originally announced February 2020.
-
The First Habitable Zone Earth-Sized Planet From TESS II: $Spitzer$ Confirms TOI-700 d
Authors:
Joseph E. Rodriguez,
Andrew Vanderburg,
Sebastian Zieba,
Laura Kreidberg,
Caroline V. Morley,
Jason D. Eastman,
Stephen R. Kane,
Alton Spencer,
Samuel N. Quinn,
Ryan Cloutier,
Chelsea X. Huang,
Karen A. Collins,
Andrew W. Mann,
Emily Gilbert,
Joshua E. Schlieder,
Elisa V. Quintana,
Thomas Barclay,
Gabrielle Suissa,
Ravi kumar Kopparapu,
Courtney D. Dressing,
George R. Ricker,
Roland K. Vanderspek,
David W. Latham,
Sara Seager,
Joshua N. Winn
, et al. (21 additional authors not shown)
Abstract:
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 ho…
▽ More
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).
△ Less
Submitted 10 June, 2020; v1 submitted 3 January, 2020;
originally announced January 2020.
-
Breaking the centrifugal barrier to giant planet contraction by magnetic disc braking
Authors:
Sivan Ginzburg,
Eugene Chiang
Abstract:
During the runaway phase of their formation, gas giants fill their gravitational spheres of influence out to Bondi or Hill radii. When runaway ends, planets shrink several orders of magnitude in radius until they are comparable in size to present-day Jupiter; in 1D models, the contraction occurs on the Kelvin-Helmholtz time-scale $t_{\rm KH}$, which is initially a few thousand years. However, if a…
▽ More
During the runaway phase of their formation, gas giants fill their gravitational spheres of influence out to Bondi or Hill radii. When runaway ends, planets shrink several orders of magnitude in radius until they are comparable in size to present-day Jupiter; in 1D models, the contraction occurs on the Kelvin-Helmholtz time-scale $t_{\rm KH}$, which is initially a few thousand years. However, if angular momentum is conserved, contraction cannot complete, as planets are inevitably spun up to their breakup periods $P_{\rm break}$. We consider how a circumplanetary disc (CPD) can de-spin a primordially magnetized gas giant and remove the centrifugal barrier, provided the disc is hot enough to couple to the magnetic field, a condition that is easier to satisfy at later times. By inferring the planet's magnetic field from its convective cooling luminosity, we show that magnetic spin-down times are shorter than contraction times throughout post-runaway contraction: $t_{\rm mag}/t_{\rm KH}\sim(P_{\rm break}/t_{\rm KH})^{1/21}\lesssim 1$. Planets can spin down until they corotate with the CPD's magnetospheric truncation radius, at a period $P_{\rm max}/P_{\rm break} \sim (t_{\rm KH}/P_{\rm break})^{1/7}$. By the time the disc disperses, $P_{\rm max}/P_{\rm break}\sim 20-30$; further contraction at fixed angular momentum can spin planets back up to $\sim$$10 P_{\rm break}$, potentially explaining observed rotation periods of giant planets and brown dwarfs.
△ Less
Submitted 25 October, 2019;
originally announced October 2019.
-
The Endgame of Gas Giant Formation: Accretion Luminosity and Contraction Post-Runaway
Authors:
Sivan Ginzburg,
Eugene Chiang
Abstract:
Giant planets are thought to form by runaway gas accretion onto solid cores. Growth must eventually stop running away, ostensibly because planets open gaps (annular cavities) in their surrounding discs. Typical models stop runaway by artificially capping the accretion rate and lowering it to zero over an arbitrarily short time-scale. In reality, post-runaway accretion persists as long as the disc…
▽ More
Giant planets are thought to form by runaway gas accretion onto solid cores. Growth must eventually stop running away, ostensibly because planets open gaps (annular cavities) in their surrounding discs. Typical models stop runaway by artificially capping the accretion rate and lowering it to zero over an arbitrarily short time-scale. In reality, post-runaway accretion persists as long as the disc remains. During this final and possibly longest phase of formation, when the planet is still emerging from the disc, its mass can more than double, and its radius contracts by orders of magnitude. By drawing from the theory of how gaps clear, we find that post-runaway accretion luminosities diverge depending on disc viscosity: luminosities fall in low-viscosity discs but continue to rise past runaway in high-viscosity discs. This divergence amounts to a factor of $10^2$ by the time the disc disperses. Irrespective of the specifics of how planets interact with discs, the observed luminosity and age of an accreting planet can be used to calculate its instantaneous mass, radius, and accretion rate. We perform this exercise for the planet candidates embedded within the discs orbiting PDS 70, HD 163296, and MWC 758, inferring masses of $1-10\, M_{\rm J}$, accretion rates of $0.1-10\, M_{\rm J}$/Myr, and radii of $1-10\, R_{\rm J}$. Our radii are computed self-consistently from the planet's concurrent contraction and accretion and do not necessarily equal the value of $2R_{\rm J}$ commonly assumed; in particular, the radius depends on the envelope opacity as $R \propto κ^{0.5}$.
△ Less
Submitted 11 October, 2019;
originally announced October 2019.
-
Circumplanetary Disk Dynamics in the Isothermal and Adiabatic Limits
Authors:
Jeffrey Fung,
Zhaohuan Zhu,
Eugene Chiang
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 November, 2019; v1 submitted 20 September, 2019;
originally announced September 2019.
-
Sculpting eccentric debris disks with eccentric gas rings
Authors:
Jonathan W. Lin,
Eugene Chiang
Abstract:
Many debris disks seen in scattered light have shapes that imply their dust grains trace highly eccentric, apsidally aligned orbits. Apsidal alignment is surprising, especially for dust. Even when born from an apse-aligned ring of parent bodies, dust grains have their periastra dispersed in all directions by stellar radiation pressure. The periastra cannot be re-oriented by planets within the shor…
▽ More
Many debris disks seen in scattered light have shapes that imply their dust grains trace highly eccentric, apsidally aligned orbits. Apsidal alignment is surprising, especially for dust. Even when born from an apse-aligned ring of parent bodies, dust grains have their periastra dispersed in all directions by stellar radiation pressure. The periastra cannot be re-oriented by planets within the short dust lifetimes at the bottom of the collisional cascade. We propose that what re-aligns dust orbits is drag exerted by second-generation gas. Gas is largely immune to radiation pressure, and when released by photodesorption or collisions within an eccentric ring of parent bodies should occupy a similarly eccentric, apse-aligned ring. Dust grains launched onto misaligned orbits cross the eccentric gas ring supersonically and can become dragged into alignment within collisional lifetimes. The resultant dust configurations, viewed nearly but not exactly edge-on, with periastra pointing away from the observer, appear moth-like, with kinked wings and even doubled pairs of wings, explaining otherwise mysterious features in HD 61005 ("The Moth") and HD 32297, including their central bulbs when we account for strong forward scattering from irregularly shaped particles. Around these systems we predict gas at Kuiper-belt-like distances to move on highly elliptical streamlines that owe their elongation, ultimately, to highly eccentric planets. Unresolved issues and an alternative explanation for apsidal alignment are outlined.
△ Less
Submitted 24 July, 2019;
originally announced July 2019.
-
The Degree of Alignment Between Circumbinary Disks and Their Binary Hosts
Authors:
Ian Czekala,
Eugene Chiang,
Sean M. Andrews,
Eric L. N. Jensen,
Guillermo Torres,
David J. Wilner,
Keivan G. Stassun,
Bruce Macintosh
Abstract:
All four circumbinary (CB) protoplanetary disks orbiting short-period ($P < 20$ day) double-lined spectroscopic binaries (SB2s)---a group that includes UZ Tau E, for which we present new ALMA data---exhibit sky-plane inclinations $i_{\rm disk}$ which match, to within a few degrees, the sky-plane inclinations $i_\star$ of their stellar hosts. Although for these systems the true mutual inclinations…
▽ More
All four circumbinary (CB) protoplanetary disks orbiting short-period ($P < 20$ day) double-lined spectroscopic binaries (SB2s)---a group that includes UZ Tau E, for which we present new ALMA data---exhibit sky-plane inclinations $i_{\rm disk}$ which match, to within a few degrees, the sky-plane inclinations $i_\star$ of their stellar hosts. Although for these systems the true mutual inclinations $θ$ between disk and binary cannot be directly measured because relative nodal angles are unknown, the near-coincidence of $i_{\rm disk}$ and $i_\star$ suggests that $θ$ is small for these most compact of systems. We confirm this hypothesis using a hierarchical Bayesian analysis, showing that 68% of CB disks around short-period SB2s have $θ< 3.0^\circ$. Near co-planarity of CB disks implies near co-planarity of CB planets discovered by Kepler, which in turn implies that the occurrence rate of close-in CB planets is similar to that around single stars. By contrast, at longer periods ranging from $30-10^5$ days (where the nodal degeneracy can be broken via, e.g., binary astrometry), CB disks exhibit a wide range of mutual inclinations, from co-planar to polar. Many of these long-period binaries are eccentric, as their component stars are too far separated to be tidally circularized. We discuss how theories of binary formation and disk-binary gravitational interactions can accommodate all these observations.
△ Less
Submitted 17 June, 2019; v1 submitted 7 June, 2019;
originally announced June 2019.
-
The End of Runaway: How Gap Opening Limits the Final Masses of Gas Giants
Authors:
Sivan Ginzburg,
Eugene Chiang
Abstract:
Gas giants are thought to form by runaway accretion: an instability driven by the self-gravity of growing atmospheres that causes accretion rates to rise super-linearly with planet mass. Why runaway should stop at a Jupiter or any other mass is unknown. We consider the proposal that final masses are controlled by circumstellar disc gaps (cavities) opened by planetary gravitational torques. We deve…
▽ More
Gas giants are thought to form by runaway accretion: an instability driven by the self-gravity of growing atmospheres that causes accretion rates to rise super-linearly with planet mass. Why runaway should stop at a Jupiter or any other mass is unknown. We consider the proposal that final masses are controlled by circumstellar disc gaps (cavities) opened by planetary gravitational torques. We develop a fully time-dependent theory of gap formation and couple it self-consistently to planetary growth rates. When gaps first open, planetary torques overwhelm viscous torques, and gas depletes as if it were inviscid. In low-viscosity discs, of the kind motivated by recent observations and theory, gaps stay predominantly in this inviscid phase and planet masses finalize at $M_{\rm final}/M_\star\sim(Ωt_{\rm disc})^{0.07}(H/a)^{2.73}(Gρ_0/Ω^2)^{1/3}$, with $M_\star$ the host stellar mass, $Ω$ the planet's orbital angular velocity, $t_{\rm disc}$ the gas disc's lifetime, $H/a$ its aspect ratio, and $ρ_0$ its unperturbed density. This final mass is independent of the dimensionless viscosity $α$ and applies to large orbital distances, typically beyond $\sim$10 AU, where disc scale heights exceed planet radii. It evaluates to a few Jupiter masses at 10-100 AU, increasing gradually with distance as gaps become harder to open.
△ Less
Submitted 9 May, 2019;
originally announced May 2019.
-
The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics From 10-100 AU
Authors:
Eric L. Nielsen,
Robert J. De Rosa,
Bruce Macintosh,
Jason J. Wang,
Jean-Baptiste Ruffio,
Eugene Chiang,
Mark S. Marley,
Didier Saumon,
Dmitry Savransky,
S. Mark Ammons,
Vanessa P. Bailey,
Travis Barman,
Celia Blain,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
Ian Czekala,
Rene Doyon,
Gaspard Duchene,
Thomas M. Esposito,
Daniel Fabrycky,
Michael P. Fitzgerald,
Katherine B. Follette,
Jonathan J. Fortney,
Benjamin L. Gerard
, et al. (40 additional authors not shown)
Abstract:
We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation be…
▽ More
We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^αa^β$, finding $α$ = -2.4 $\pm$ 0.8 and $β$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.
△ Less
Submitted 10 April, 2019;
originally announced April 2019.
-
The Mass of Stirring Bodies in the AU Mic Debris Disk Inferred from Resolved Vertical Structure
Authors:
Cail Daley,
A. Meredith Hughes,
Evan S. Carter,
Kevin Flaherty,
Zachary Lambros,
Margaret Pan,
Hilke Schlichting,
Eugene Chiang,
Mark Wyatt,
David Wilner,
Sean Andrews,
John Carpenter
Abstract:
The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but…
▽ More
The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution.
△ Less
Submitted 8 April, 2019; v1 submitted 29 March, 2019;
originally announced April 2019.
-
Dynamical Constraints on the HR 8799 Planets with GPI
Authors:
Jason J. Wang,
James R. Graham,
Rebekah Dawson,
Daniel Fabrycky,
Robert J. De Rosa,
Laurent Pueyo,
Quinn Konopacky,
Bruce Macintosh,
Christian Marois,
Eugene Chiang,
S. Mark Ammons,
Pauline Arriaga,
Vanessa P. Bailey,
Travis Barman,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
René Doyon,
Gaspard Duchêne,
Thomas M. Esposito,
Michael P. Fitzgerald,
Katherine B. Follette,
Benjamin L. Gerard,
Stephen J. Goodsell,
Alexandra Z. Greenbaum
, et al. (30 additional authors not shown)
Abstract:
The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assu…
▽ More
The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.
△ Less
Submitted 11 September, 2018;
originally announced September 2018.
-
Multiple Disk Gaps and Rings Generated by a Single Super-Earth: II. Spacings, Depths, and Number of Gaps, with Application to Real Systems
Authors:
Ruobing Dong,
Shengtai Li,
Eugene Chiang,
Hui Li
Abstract:
ALMA has found multiple dust gaps and rings in a number of protoplanetary disks in continuum emission at millimeter wavelengths. The origin of such structures is in debate. Recently, we documented how one super-Earth planet can open multiple (up to five) dust gaps in a disk with low viscosity ($α\lesssim10^{-4}$). In this paper, we examine how the positions, depths, and total number of gaps opened…
▽ More
ALMA has found multiple dust gaps and rings in a number of protoplanetary disks in continuum emission at millimeter wavelengths. The origin of such structures is in debate. Recently, we documented how one super-Earth planet can open multiple (up to five) dust gaps in a disk with low viscosity ($α\lesssim10^{-4}$). In this paper, we examine how the positions, depths, and total number of gaps opened by one planet depend on input parameters, and apply our results to real systems. Gap locations (equivalently, spacings) are the easiest metric to use when making comparisons between theory and observations, as positions can be robustly measured. We fit the locations of gaps empirically as functions of planet mass and disk aspect ratio. We find that the locations of the double gaps in HL Tau and TW Hya, and of all three gaps in HD 163296, are consistent with being opened by a sub-Saturn mass planet. This scenario predicts the locations of other gaps in HL Tau and TW Hya, some of which appear consistent with current observations. We also show how the Rossby wave instability may develop at the edges of several gaps and result in multiple dusty vortices, all caused by one planet. A planet as low in mass as Mars may produce multiple dust gaps in the terrestrial planet forming region.
△ Less
Submitted 14 September, 2018; v1 submitted 20 August, 2018;
originally announced August 2018.
-
Direct Imaging of the HD 35841 Debris Disk: A Polarized Dust Ring from Gemini Planet Imager and an Outer Halo from HST/STIS
Authors:
Thomas M. Esposito,
Gaspard Duchêne,
Paul Kalas,
Malena Rice,
Élodie Choquet,
Bin Ren,
Marshall D. Perrin,
Christine H. Chen,
Pauline Arriaga,
Eugene Chiang,
Eric L. Nielsen,
James R. Graham,
Jason J. Wang,
Robert J. De Rosa,
Katherine B. Follette,
S. Mark Ammons,
Megan Ansdell,
Vanessa P. Bailey,
Travis Barman,
Juan Sebastián Bruzzone,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
Rene Doyon,
Michael P. Fitzgerald
, et al. (33 additional authors not shown)
Abstract:
We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending…
▽ More
We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.
△ Less
Submitted 7 June, 2018;
originally announced June 2018.
-
A Balanced Budget View on Forming Giant Planets by Pebble Accretion
Authors:
Jonathan W. Lin,
Eve J. Lee,
Eugene Chiang
Abstract:
Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with th…
▽ More
Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with the largest accretion cross-sections also have the fastest radial drift speeds and are the most easily drained out of discs. We present a global model of planet formation by pebble accretion that keeps track of the disc's mass budget. Cores, each initialized with a lunar mass, grow from discs whose finite stores of mm-cm sized pebbles drift inward across all radii in viscously accreting gas. For every 1 $M_\oplus$ netted by a core, at least 10 $M_\oplus$ and possibly much more are lost to radial drift. Core growth rates are typically exponentially sensitive to particle Stokes number, turbulent Mach number, and solid surface density. This exponential sensitivity, when combined with disc migration, tends to generate binary outcomes from 0.1-30 AU: either sub-Earth cores remain sub-Earth, or explode into Jupiters, with the latter migrating inward to varying degrees. When Jupiter-breeding cores assemble from mm-cm sized pebbles, they do so in discs where such particles drain out in $\sim$10$^5$ yr or less; such fast-draining discs do not fit mm-wave observations.
△ Less
Submitted 6 August, 2018; v1 submitted 1 June, 2018;
originally announced June 2018.
-
The Eccentric Cavity, Triple Rings, Two-Armed Spirals, and Double Clumps of the MWC 758 Disk
Authors:
Ruobing Dong,
Sheng-yuan Liu,
Josh Eisner,
Sean Andrews,
Jeffrey Fung,
Zhaohuan Zhu,
Eugene Chiang,
Jun Hashimoto,
Hauyu Baobab Liu,
Simon Casassus,
Thomas Esposito,
Yasuhiro Hasegawa,
Takayuki Muto,
Yaroslav Pavlyuchenkov,
David Wilner,
Eiji Akiyama,
Motohide Tamura,
John Wisniewski
Abstract:
Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution…
▽ More
Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $μ$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.
△ Less
Submitted 30 May, 2018;
originally announced May 2018.
-
A Decade of MWC 758 Disk Images: Where Are the Spiral-Arm-Driving Planets?
Authors:
Bin Ren,
Ruobing Dong,
Thomas M. Esposito,
Laurent Pueyo,
John H. Debes,
Charles A. Poteet,
Élodie Choquet,
Myriam Benisty,
Eugene Chiang,
Carol A. Grady,
Dean C. Hines,
Glenn Schneider,
Rémi Soummer
Abstract:
Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and co-rotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10-yr timescale, we are able to provide dynamical constraint…
▽ More
Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and co-rotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10-yr timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 $L'$-band observations in 2017. MWC~758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of ${0.6^\circ}^{+3.3^\circ}_{-0.6^\circ}\, \mathrm{yr}^{-1}$ at $3σ$ for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best fit orbit is at $89$ au ($0.59"$), with a $3σ$ lower limit of 30 au ($0.20"$). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.
△ Less
Submitted 19 June, 2018; v1 submitted 18 March, 2018;
originally announced March 2018.
-
Secular Dynamics of an Exterior Test Particle: The Inverse Kozai and Other Eccentricity-Inclination Resonances
Authors:
Benjamin R. Vinson,
Eugene Chiang
Abstract:
The behavior of an interior test particle in the secular 3-body problem has been studied extensively. A well-known feature is the Lidov-Kozai resonance in which the test particle's argument of periapse librates about $\pm 90^\circ$ and large oscillations in eccentricity and inclination are possible. Less explored is the inverse problem: the dynamics of an exterior test particle and an interior per…
▽ More
The behavior of an interior test particle in the secular 3-body problem has been studied extensively. A well-known feature is the Lidov-Kozai resonance in which the test particle's argument of periapse librates about $\pm 90^\circ$ and large oscillations in eccentricity and inclination are possible. Less explored is the inverse problem: the dynamics of an exterior test particle and an interior perturber. We survey numerically the inverse secular problem, expanding the potential to hexadecapolar order and correcting an error in the published expansion. Four secular resonances are uncovered that persist in full $N$-body treatments (in what follows, $\varpi$ and $Ω$ are the longitudes of periapse and of ascending node, $ω$ is the argument of periapse, and subscripts 1 and 2 refer to the inner perturber and outer test particle): (i) an orbit-flipping quadrupole resonance requiring a non-zero perturber eccentricity $e_1$, in which $Ω_2-\varpi_1$ librates about $\pm 90^\circ$; (ii) a hexadecapolar resonance (the "inverse Kozai" resonance) for perturbers that are circular or nearly so and inclined by $I \simeq 63^\circ/117^\circ$, in which $ω_2$ librates about $\pm 90^\circ$ and which can vary the particle eccentricity by $Δe_2 \simeq 0.2$ and lead to orbit crossing; (iii) an octopole "apse-aligned" resonance at $I \simeq 46^\circ/107^\circ$ wherein $\varpi_2 - \varpi_1$ librates about $0^\circ$ and $Δe_2$ grows with $e_1$; and (iv) an octopole resonance at $I \simeq 73^\circ/134^\circ$ wherein $\varpi_2 + \varpi_1 - 2 Ω_2$ librates about $0^\circ$ and $Δe_2$ can be as large as 0.3 for small $e_1 \neq 0$. The more eccentric the perturber, the more the particle's eccentricity and inclination vary; also, more polar orbits are more chaotic. Our inverse solutions may be applied to the Kuiper belt and debris disks, circumbinary planets, and stellar systems.
△ Less
Submitted 20 December, 2017; v1 submitted 28 November, 2017;
originally announced November 2017.
-
Optically Thin Core Accretion: How Planets Get Their Gas in Nearly Gas-Free Disks
Authors:
Eve J. Lee,
Eugene Chiang,
Jason W. Ferguson
Abstract:
Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such en…
▽ More
Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such environments can be optically thin to both incident starlight and internally generated thermal radiation. We construct time-dependent models of such atmospheres, showing that super-Earths/sub-Neptunes can accrete their $\sim$1%-by-mass gas envelopes, and super-puffs/sub-Saturns their $\sim$20%-by-mass envelopes, over a wide range of nebular depletion histories requiring no fine tuning. Although nascent atmospheres can exhibit stratospheric temperature inversions effected by atomic Fe and various oxides that absorb strongly at visible wavelengths, the rate of gas accretion remains controlled by the radiative-convective boundary (rcb) at much greater pressures. For dusty envelopes, the temperature at the rcb $T_{\rm rcb} \simeq 2500$ K is still set by ${\rm H}_2$ dissociation; for dust-depleted envelopes, $T_{\rm rcb}$ tracks the temperature of the visible or thermal photosphere, whichever is deeper, out to at least $\sim$5 AU. The rate of envelope growth remains largely unchanged between the old radiative diffusion models and the new optically thin models, reinforcing how robustly super-Earths form as part of the endgame chapter in disk evolution.
△ Less
Submitted 6 February, 2018; v1 submitted 6 October, 2017;
originally announced October 2017.
-
Stellar Winds and Dust Avalanches in the AU Mic Debris Disk
Authors:
Eugene Chiang,
Jeffrey Fung
Abstract:
We explain the fast-moving, ripple-like features in the edge-on debris disk orbiting the young M dwarf AU Mic. The bright features are clouds of sub-micron dust repelled by the host star's wind. The clouds are produced by avalanches: radial outflows of dust that gain exponentially more mass as they shatter background disk particles in collisional chain reactions. The avalanches are triggered from…
▽ More
We explain the fast-moving, ripple-like features in the edge-on debris disk orbiting the young M dwarf AU Mic. The bright features are clouds of sub-micron dust repelled by the host star's wind. The clouds are produced by avalanches: radial outflows of dust that gain exponentially more mass as they shatter background disk particles in collisional chain reactions. The avalanches are triggered from a region a few AU across -- the "avalanche zone" -- located on AU Mic's primary "birth" ring, at a true distance of $\sim$35 AU from the star but at a projected distance more than a factor of 10 smaller: the avalanche zone sits directly along the line of sight to the star, on the side of the ring nearest Earth, launching clouds that disk rotation sends wholly to the southeast, as observed. The avalanche zone marks where the primary ring intersects a secondary ring of debris left by the catastrophic disruption of a progenitor up to Varuna in size, less than tens of thousands of years ago. Only where the rings intersect are particle collisions sufficiently violent to spawn the sub-micron dust needed to seed the avalanches. We show that this picture works quantitatively, reproducing the masses, sizes, and velocities of the observed escaping clouds. The Lorentz force exerted by the wind's magnetic field, whose polarity reverses periodically according to the stellar magnetic cycle, promises to explain the observed vertical undulations. The timescale between avalanches, about 10 yr, might be set by time variability of the wind mass-loss rate or, more speculatively, by some self-regulating limit cycle.
△ Less
Submitted 30 August, 2017; v1 submitted 27 July, 2017;
originally announced July 2017.
-
A Three-Dimensional View of Turbulence: Constraints on Turbulent Motions in the HD 163296 Protoplanetary Disk using DCO$^+$
Authors:
Kevin M. Flaherty,
A. Meredith Hughes,
Sanaea C. Rose,
Jacob B. Simon,
Chunhua Qi,
Sean M. Andrews,
Agnes Kospal,
David J. Wilner,
Eugene Chiang,
Philip J. Armitage,
Xue-ning Bai
Abstract:
Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combine…
▽ More
Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $α\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment.
△ Less
Submitted 14 June, 2017;
originally announced June 2017.
-
Multiple Disk Gaps and Rings Generated by a Single Super-Earth
Authors:
Ruobing Dong,
Shengtai Li,
Eugene Chiang,
Hui Li
Abstract:
We investigate the observational signatures of super-Earths (i.e., Earth-to-Neptune mass planets) in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of sub-mm-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the mm-wave thermal c…
▽ More
We investigate the observational signatures of super-Earths (i.e., Earth-to-Neptune mass planets) in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of sub-mm-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the mm-wave thermal continuum for direct comparison with observations. In low viscosity gas ($α\lesssim10^{-4}$), a super-Earth opens two annular gaps to either side of its orbit by the action of Lindblad torques. This double gap and its associated gas pressure gradients cause dust particles to be dragged by gas into three rings: one ring sandwiched between the two gaps, and two rings located at the gap edges farthest from the planet. Depending on system parameters, additional rings may manifest for a single planet. A double gap located at tens of AUs from a host star in Taurus can be detected in the dust continuum by the Atacama Large Millimeter Array (ALMA) at an angular resolution of ~0".03 after two hours of integration. Ring and gap features persist in a variety of background disk profiles, last for thousands of orbits, and change their relative positions and dimensions depending on the speed and direction of planet migration. Candidate double gaps have been observed by ALMA in systems like HL Tau (D5 and D6) and TW Hya (at 37 and 43 AU); we submit that each double gap is carved by one super-Earth in nearly inviscid gas.
△ Less
Submitted 15 May, 2017; v1 submitted 12 May, 2017;
originally announced May 2017.