-
Spectroastrometric Survey of Protoplanetary Disks with Inner Dust Cavities
Authors:
Stanley K. Jensen Jr,
Sean D. Brittain,
Andrea Banzatti,
Joan R. Najita,
John S. Carr,
Joshua Kern,
Janus Kozdon,
Jonathan Zrake,
Jeffrey Fung
Abstract:
We present high-resolution spectra and spectroastrometric (SA) measurements of fundamental rovibrational CO emission from nine nearby ($\lesssim$300 pc) protoplanetary disks where large inner dust cavities have been observed. The emission line profiles and SA signals are fit with a slab disk model that allows the eccentricity of the disk and intensity of the emission to vary as power laws. Six of…
▽ More
We present high-resolution spectra and spectroastrometric (SA) measurements of fundamental rovibrational CO emission from nine nearby ($\lesssim$300 pc) protoplanetary disks where large inner dust cavities have been observed. The emission line profiles and SA signals are fit with a slab disk model that allows the eccentricity of the disk and intensity of the emission to vary as power laws. Six of the sources are well fit with our model, and three of these sources show asymmetric line profiles that can be fit by adopting a non-zero eccentricity. The three other sources have components in either their line profile or SA signal that are not captured by our disk model. Two of these sources (V892 Tau, CQ Tau) have multi-epoch observations that reveal significant variability. CQ Tau and AB Aur have CO line profiles with centrally-peaked components that are similar to line profiles that have been interpreted as evidence of molecular gas arising from a wide-angle disk wind. Alternatively, emission from a circumplanetary disk (CPD) could also account for this component. The interpretations of these results can be clarified in the future with additional epochs that will test the variability timescale of these SA signals. We discuss the utility of using high-resolution spectroscopy for probing the dynamics of gas in the disk and the scenarios that can give rise to profiles that are not fit with a simple disk model.
△ Less
Submitted 30 January, 2024; v1 submitted 28 December, 2023;
originally announced December 2023.
-
A Systematic Study of Planetary Envelope Growth with 3D Radiation-Hydrodynamics Simulations
Authors:
Avery Bailey,
Jim Stone,
Jeffrey Fung
Abstract:
In the core accretion model of planet formation, envelope cooling regulates the accretion of material and ultimately sets the timescale to form a giant planet. Given the diversity of planet-forming environments, opacity uncertainties, and the advective transport of energy by 3-dimensional recycling flows, it is unclear whether 1D models can adequately describe envelope structure and accretion in a…
▽ More
In the core accretion model of planet formation, envelope cooling regulates the accretion of material and ultimately sets the timescale to form a giant planet. Given the diversity of planet-forming environments, opacity uncertainties, and the advective transport of energy by 3-dimensional recycling flows, it is unclear whether 1D models can adequately describe envelope structure and accretion in all regimes. Even in 3D models, it is unclear whether approximate radiative transfer methods sufficiently model envelope cooling particularly at the planetary photosphere. To address these uncertainties, we present a suite of 3D radiation hydrodynamics simulations employing methods that directly solve the transfer equation. We perform a parameter space study, formulated in terms of dimensionless parameters, for a variety of envelope optical depths and cooling times. We find that the thermodynamic structure of the envelope ranges from adiabatic to isothermal based on the cooling time and by extension, the background disk temperature and density. Our models show general agreement with 1D static calculations, suggesting a limited role of recycling flows in determining envelope structure. By adopting a dimensionless framework, these models can be applied to a wide range of formation conditions and assumed opacities. In particular, we dimensionalize them to the case of a super-Earth and proto-Jupiter and place upper limits on the 3D mass accretion rates prior to runaway growth. Finally, we evaluate the fidelity of approximate radiative transfer methods and find that even in the most challenging cases, more approximate methods are sufficiently accurate and worth their savings in computational cost.
△ Less
Submitted 4 October, 2023;
originally announced October 2023.
-
Ro-vibrational Spectroscopy of CI Tau -- Evidence of a Multi-Component Eccentric Disk Induced by a Planet
Authors:
Janus Kozdon,
Sean Brittain,
Jeffrey Fung,
Josh Kern,
Stanley Jensen,
John Carr,
Joan Najita,
Andrea Banzatti
Abstract:
CI Tau is currently the only T Tauri star with an inner protoplanetary disk that hosts a planet, CI Tau b, that has been detected by a radial velocity survey. This provides the unique opportunity to study disk features that were imprinted by that planet. We present multi-epoch spectroscopic data, taken with NASA IRTF in 2022, of the ${}^{12}$CO and hydrogen Pf$β$ line emissions spanning 9 consecut…
▽ More
CI Tau is currently the only T Tauri star with an inner protoplanetary disk that hosts a planet, CI Tau b, that has been detected by a radial velocity survey. This provides the unique opportunity to study disk features that were imprinted by that planet. We present multi-epoch spectroscopic data, taken with NASA IRTF in 2022, of the ${}^{12}$CO and hydrogen Pf$β$ line emissions spanning 9 consecutive nights, which is the proposed orbital period of CI Tau b. We find that the star's accretion rate varied according to that 9~d period, indicative of companion driven accretion. Analysis of the ${}^{12}$CO emission lines reveals that the disk can be described with an inner and outer component spanning orbital radii 0.05-0.13~au and 0.15-1.5~au, respectively. Both components have eccentricities of about 0.05 and arguments of periapses that are oppositely aligned. We present a proof-of-concept hydrodynamic simulation that shows a massive companion on a similarly eccentric orbit can recreate a similar disk structure. Our results allude to such a companion being located around an orbital distance of 0.14~au. However, this planet's orbital parameters may be inconsistent with those of CI Tau b whose high eccentricity is likely not compatible with the low disk eccentricities inferred by our model.
△ Less
Submitted 7 September, 2023; v1 submitted 24 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.
-
Planet-Disk Interactions
Authors:
Sijme-Jan Paardekooper,
Ruobing Dong,
Paul Duffell,
Jeffrey Fung,
Frederic S. Masset,
Gordon Ogilvie,
Hidekazu Tanaka
Abstract:
Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments.…
▽ More
Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments. The back-reaction of perturbations induced in the disk, both wave-like and non-wavelike, is a change in orbital elements of the planet. The efficiency of orbital migration is a long-standing problem in planet formation theory. We discuss recent progress in planet-disk interactions for different planet masses and disk parameters, in particular the level of turbulence, and progress in modeling observational signatures of embedded planets.
△ Less
Submitted 8 November, 2022; v1 submitted 17 March, 2022;
originally announced March 2022.
-
The Triangulum Extended (TREX) Survey: The Stellar Disk Dynamics of M33 as a Function of Stellar Age
Authors:
A. C. N. Quirk,
P. Guhathakurta,
K. Gilbert,
L. Chemin,
J. Dalcanton,
B. Williams,
A. Seth,
E. Patel,
J. Fung,
P. Tangirala,
I. Yusufali
Abstract:
Triangulum, M33, is a low mass, relatively undisturbed spiral galaxy that offers a new regime in which to test models of dynamical heating. In spite of its proximity, the dynamical heating history of M33 has not yet been well constrained. In this work, we present the TREX Survey, the largest stellar spectroscopic survey across the disk of M33. We present the stellar disk kinematics as a function o…
▽ More
Triangulum, M33, is a low mass, relatively undisturbed spiral galaxy that offers a new regime in which to test models of dynamical heating. In spite of its proximity, the dynamical heating history of M33 has not yet been well constrained. In this work, we present the TREX Survey, the largest stellar spectroscopic survey across the disk of M33. We present the stellar disk kinematics as a function of age to study the past and ongoing dynamical heating of M33. We measure line of sight velocities for ~4,500 disk stars. Using a subset, we divide the stars into broad age bins using Hubble Space Telescope and Canada-France-Hawaii-Telescope photometric catalogs: massive main sequence stars and helium burning stars (~80 Myr), intermediate mass asymptotic branch stars (~1 Gyr), and low mass red giant branch stars (~4 Gyr). We compare the stellar disk dynamics to that of the gas using existing HI, CO, and Halpha kinematics. We find that the disk of M33 has relatively low velocity dispersion (~16 km/s), and unlike in the Milky Way and Andromeda galaxies, there is no strong trend in velocity dispersion as a function of stellar age. The youngest disk stars are as dynamically hot as the oldest disk stars and are dynamically hotter than predicted by most M33 like low mass simulated analogs in Illustris. The velocity dispersion of the young stars is highly structured, with the large velocity dispersion fairly localized. The cause of this high velocity dispersion is not evident from the observations and simulated analogs presented here.
△ Less
Submitted 9 February, 2022;
originally announced February 2022.
-
Dust Dynamics in Transitional Disks: Clumping and Disk Recession
Authors:
Jiaqing Bi,
Jeffrey Fung
Abstract:
The role of radiation pressure in dust migration and the opening of inner cavities in transitional disks is revisited in this paper. Dust dynamics including radiation pressure is often studied in axisymmetric models, but in this work, we show that highly non-axisymmetric features can arise from an instability at the inner disk edge. Dust grains clump into high density features there, allowing radi…
▽ More
The role of radiation pressure in dust migration and the opening of inner cavities in transitional disks is revisited in this paper. Dust dynamics including radiation pressure is often studied in axisymmetric models, but in this work, we show that highly non-axisymmetric features can arise from an instability at the inner disk edge. Dust grains clump into high density features there, allowing radiation to leak around them and penetrate deeper into the disk, changing the course of dust migration. Our proof-of-concept, two-dimensional, vertically-averaged simulations show that the combination of radiation pressure, shadowing, and gas drag can produce a net outward migration, or recession, of the dust component of the disk. The recession speed of the inner disk edge is on the order of $10^{-5}$ times Keplerian speed in our parameter space, which is faster than the background viscous flow, assuming a Shakura & Sunyaev viscosity alpha $\lesssim 10^{-3}$. This speed, if sustained over the lifetime of the disk, can result in a dust cavity as large as tens of au.
△ Less
Submitted 14 February, 2022; v1 submitted 4 February, 2022;
originally announced February 2022.
-
The TREX Survey: Kinematical Complexity Throughout M33's Stellar Disk and Evidence for a Stellar Halo
Authors:
Karoline M. Gilbert,
Amanda C. N. Quirk,
Puragra Guhathakurta,
Erik Tollerud,
Jennifer Wojno,
Julianne J. Dalcanton,
Meredith J. Durbin,
Anil Seth,
Benjamin F. Williams,
Justin T. Fung,
Pujita Tangirala,
Ibrahim Yusufali
Abstract:
We present initial results from a large spectroscopic survey of stars throughout M33's stellar disk. We analyze a sample of 1667 red giant branch (RGB) stars extending to projected distances of $\sim 11$ kpc from M33's center ($\sim 18$ kpc, or $\sim 10$ scale lengths, in the plane of the disk). The line-of-sight velocities of RGB stars show the presence of two kinematical components. One componen…
▽ More
We present initial results from a large spectroscopic survey of stars throughout M33's stellar disk. We analyze a sample of 1667 red giant branch (RGB) stars extending to projected distances of $\sim 11$ kpc from M33's center ($\sim 18$ kpc, or $\sim 10$ scale lengths, in the plane of the disk). The line-of-sight velocities of RGB stars show the presence of two kinematical components. One component is consistent with rotation in the plane of M33's HI disk and has a velocity dispersion ($\sim 19$ km s$^{-1}$) consistent with that observed in a comparison sample of younger stars, while the second component has a significantly higher velocity dispersion. A two-component fit to the RGB velocity distribution finds that the high dispersion component has a velocity dispersion of $59.3^{+2.6}_{-2.5}$ km s$^{-1}$ and rotates very slowly in the plane of the disk (consistent with no rotation at the $<1.5σ$ level), which favors interpreting it as a stellar halo rather than a thick disk population. A spatial analysis indicates that the fraction of RGB stars in the high-velocity-dispersion component decreases with increasing radius over the range covered by the spectroscopic sample. Our spectroscopic sample establishes that a significant high-velocity-dispersion component is present in M33's RGB population from near M33's center to at least the radius where M33's HI disk begins to warp at 30$'$ ($\sim 7.5$ kpc) in the plane of the disk. This is the first detection and spatial characterization of a kinematically hot stellar component throughout M33's inner regions.
△ Less
Submitted 29 October, 2021;
originally announced October 2021.
-
Cooling-Induced Vortex Decay in Keplerian Disks
Authors:
Jeffrey Fung,
Tomohiro Ono
Abstract:
Vortices are readily produced by hydrodynamical instabilities, such as the Rossby wave instability, in protoplanetary disks. However, large-scale asymmetries indicative of dust-trapping vortices are uncommon in sub-millimeter continuum observations. One possible explanation is that vortices have short lifetimes. In this paper, we explore how radiative cooling can lead to vortex decay. Elliptical v…
▽ More
Vortices are readily produced by hydrodynamical instabilities, such as the Rossby wave instability, in protoplanetary disks. However, large-scale asymmetries indicative of dust-trapping vortices are uncommon in sub-millimeter continuum observations. One possible explanation is that vortices have short lifetimes. In this paper, we explore how radiative cooling can lead to vortex decay. Elliptical vortices in Keplerian disks go through adiabatic heating and cooling cycles. Radiative cooling modifies these cycles and generates baroclinicity that changes the potential vorticity of the vortex. We show that the net effect is typically a spin down, or decay, of the vortex for a sub-adiabatic radial stratification. We perform a series of two-dimensional shearing box simulations, varying the gas cooling (or relaxation) time, $t_{\rm cool}$, and initial vortex strength. We measure the vortex decay half-life, $t_{\rm half}$, and find that it can be roughly predicted by the timescale ratio $t_{\rm cool}/t_{\rm turn}$, where $t_{\rm turn}$ is the vortex turnaround time. Decay is slow in both the isothermal ($t_{\rm cool}\ll t_{\rm turn}$) and adiabatic ($t_{\rm cool}\gg t_{\rm turn}$) limits; it is fastest when $t_{\rm cool}\sim0.1\,t_{\rm turn}$, where $t_{\rm half}$ is as short as $\sim300$ orbits. At tens of au where disk rings are typically found, $t_{\rm turn}$ is likely much longer than $t_{\rm cool}$, potentially placing vortices in the fast decay regime.
△ Less
Submitted 13 August, 2021;
originally announced August 2021.
-
Observational Signatures of Planets in Protoplanetary Disks: Temperature structures in spiral arms
Authors:
Dhruv Muley,
Ruobing Dong,
Jeffrey Fung
Abstract:
High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with $β$-cooling to a vertically-stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a c…
▽ More
High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with $β$-cooling to a vertically-stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a contribution from stellar irradiation, 0.3-3% inside the planet radius but always <0.5% outside. When cooling time equals or exceeds dynamical time, however, this is overwhelmed by hydrodynamic PdV work, which introduces a 10-20% perturbation within a factor of 2 from the planet's orbital radius. We devise an empirical fit of the spiral amplitude $Δ(T)$ to take into account both effects. Where cooling is slow, we find also that temperature perturbations from buoyancy spirals -- a strictly 3D, non-isothermal phenomenon -- become nearly as strong as those from Lindblad spirals, which are amenable to 2D and isothermal studies. Our findings may help explain observed thermal features in disks like TW Hydrae and CQ Tauri, and underscore that 3D effects have a qualitatively important effect on disk structure.
△ Less
Submitted 13 July, 2021;
originally announced July 2021.
-
The Three Dimensional Flow Field Around Planets on Eccentric Orbits
Authors:
Avery Bailey,
Jim Stone,
Jeffrey Fung
Abstract:
We investigate the properties of the hydrodynamic flow around eccentric protoplanets and compare them with the often assumed case of a circular orbit. To this end, we perform a set of 3D hydrodynamic simulations of protoplanets with small eccentricities ($e\leq 0.1$). We adopt an isothermal equation of state and concentrate resolution on the protoplanet to investigate flows down to the scale of th…
▽ More
We investigate the properties of the hydrodynamic flow around eccentric protoplanets and compare them with the often assumed case of a circular orbit. To this end, we perform a set of 3D hydrodynamic simulations of protoplanets with small eccentricities ($e\leq 0.1$). We adopt an isothermal equation of state and concentrate resolution on the protoplanet to investigate flows down to the scale of the protoplanet's circumplanetary disk (CPD). We find enhanced prograde rotation exterior to the CPD for low planet masses undergoing subsonic eccentric motion. If the eccentricity is made large enough to develop a bow shock, this trend reverses and rotation becomes increasingly retrograde. The instantaneous eccentric flow field is dramatically altered compared to circular orbits. Whereas the latter exhibit a generic pattern of polar inflow and midplane outflow, the flow geometry depends on orbital phase in the eccentric case. For even the modest eccentricities tested here, the dominant source of inflow can come from the midplane instead of the poles. We find that the amount of inflow and outflow increases for higher $e$ and lower protoplanet masses, thereby recycling more gas through the planet's Bondi radius. These increased fluxes may increase the pebble accretion rate for eccentric planets up to several times that of the circular orbit rate. In response to eccentric motion, the structure and rotation of the planet's bound CPD remains unchanged. Because the CPD regulates the eventual accretion of gas onto the planet, we predict little change to the gas accretion rates between eccentric and circular planets.
△ Less
Submitted 30 October, 2020;
originally announced November 2020.
-
Visualizing the Kinematics of Planet Formation
Authors:
Disk Dynamics Collaboration,
Philip J. Armitage,
Jaehan Bae,
Myriam Benisty,
Edwin A. Bergin,
Simon Casassus,
Ian Czekala,
Stefano Facchini,
Jeffrey Fung,
Cassandra Hall,
John D. Ilee,
Miriam Keppler,
Aleksandra Kuznetsova,
Romane Le Gal,
Ryan A. Loomis,
Wladimir Lyra,
Natascha Manger,
Sebastian Perez,
Christophe Pinte,
Daniel J. Price,
Giovanni Rosotti,
Judit Szulagyi,
Kamber Schwarz,
Jacob B. Simon,
Richard Teague
, et al. (1 additional authors not shown)
Abstract:
A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on…
▽ More
A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the 'Visualizing the Kinematics of Planet Formation' workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.
△ Less
Submitted 15 September, 2020; v1 submitted 9 September, 2020;
originally announced September 2020.
-
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.
-
A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag
Authors:
Jeffrey Fung,
Dhruv Muley
Abstract:
Numerical simulations of dust-gas dynamics are one of the fundamental tools in astrophysical research, such as the study of star and planet formation. It is common to find tightly coupled dust and gas in astrophysical systems, which demands that any practical integration method be able to take time steps $Δt$ much longer than the stopping time $t_{\rm s}$ due to drag. A number of methods have been…
▽ More
Numerical simulations of dust-gas dynamics are one of the fundamental tools in astrophysical research, such as the study of star and planet formation. It is common to find tightly coupled dust and gas in astrophysical systems, which demands that any practical integration method be able to take time steps $Δt$ much longer than the stopping time $t_{\rm s}$ due to drag. A number of methods have been developed to ensure stability in this stiff ($Δt\gg t_{\rm s}$) regime, but there remains large room for improvement in terms of accuracy. In this paper, we describe an easy-to-implement method, the "staggered semi-analytic method" (SSA), and conduct numerical tests to compare it to other implicit and semi-analytic methods, including the $2^{\rm nd}$ order implicit method and the Verlet method. SSA makes use of a staggered step to better approximate the terminal velocity in the stiff regime. In applications to protoplanetary disks, this not only leads to orders-of-magnitude higher accuracy than the other methods, but also provides greater stability, making it possible to take time steps 100 times larger in some situations. SSA is also $2^{\rm nd}$ order accurate and symplectic when $Δt \ll t_{\rm s}$. More generally, the robustness of SSA makes it applicable to linear dust-gas drag in virtually any context.
△ Less
Submitted 18 September, 2019; v1 submitted 4 September, 2019;
originally announced September 2019.
-
PDS 70: A transition disk sculpted by a single planet
Authors:
Dhruv Muley,
Jeffrey Fung,
Nienke van der Marel
Abstract:
The wide, deep cavities of transition disks are often believed to have been hollowed out by nascent planetary systems. PDS 70, a ${\sim}5$ Myr old transition disk system in which a multi-Jupiter-mass planet candidate at 22 au coexists with a ${\sim}30$ au gas and ${\sim}60$ au dust-continuum gap, provides a valuable case study for this hypothesis. Using the PEnGUIn hydrodynamics code, we simulate…
▽ More
The wide, deep cavities of transition disks are often believed to have been hollowed out by nascent planetary systems. PDS 70, a ${\sim}5$ Myr old transition disk system in which a multi-Jupiter-mass planet candidate at 22 au coexists with a ${\sim}30$ au gas and ${\sim}60$ au dust-continuum gap, provides a valuable case study for this hypothesis. Using the PEnGUIn hydrodynamics code, we simulate the orbital evolution and accretion of PDS 70b in its natal disk over the lifetime of the system. When the accreting planet reaches about 2.5 Jupiter masses, it spontaneously grows in eccentricity and consumes material from a wide swathe of the PDS 70 disk; radiative transfer post-processing with DALI shows that this accurately reproduces the observed gap profile. Our results demonstrate that super-Jupiter planets can single-handedly carve out transition disk cavities, and indicate that the high eccentricities measured for such giants may be a natural consequence of disk-planet interaction.
△ Less
Submitted 29 May, 2019; v1 submitted 19 February, 2019;
originally announced February 2019.
-
Observational Signatures of Planets in Protoplanetary Disks: Planet-Induced Line Broadening in Gaps
Authors:
Ruobing Dong,
Sheng-yuan Liu,
Jeffrey Fung
Abstract:
Protoplanets can produce structures in protoplanetary disks via gravitational disk-planet interactions. Once detected, such structures serve as signposts of planet formation. Here we investigate the kinematic signatures in disks produced by multi-Jupiter mass ($M_{\rm J}$) planets using 3D hydrodynamics and radiative transfer simulations. Such a planet opens a deep gap, and drives transonic vertic…
▽ More
Protoplanets can produce structures in protoplanetary disks via gravitational disk-planet interactions. Once detected, such structures serve as signposts of planet formation. Here we investigate the kinematic signatures in disks produced by multi-Jupiter mass ($M_{\rm J}$) planets using 3D hydrodynamics and radiative transfer simulations. Such a planet opens a deep gap, and drives transonic vertical motions inside. Such motions include both a bulk motion of the entire half-disk column, and turbulence on scales comparable to and smaller than the scale height. They significantly broaden molecular lines from the gap, producing double-peaked line profiles at certain locations, and a kinematic velocity dispersion comparable to thermal after azimuthal averaging. The same planet does not drive fast vertical motions outside the gap, except at the inner spiral arms and the disk surface. Searching for line broadening induced by multi-$M_{\rm J}$ planets inside gaps requires an angular resolution comparable to the gap width, an assessment of the gap gas temperature to within a factor of 2, and a high sensitivity needed to detect line emission from the gap.
△ Less
Submitted 14 May, 2019; v1 submitted 23 November, 2018;
originally announced November 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.
-
Subaru/HiCIAO $HK_{\rm s}$ imaging of LkH$α$ 330 - multi-band detection of the gap and spiral-like structures
Authors:
Taichi Uyama,
Jun Hashimoto,
Takayuki Muto,
Eiji Akiyama,
Ruobing Dong,
Jerome de Leon,
Itsuki Sakon,
Tomoyuki Kudo,
Nobuhiko Kusakabe,
Masayuki Kuzuhara,
Mickael Bonnefoy,
Lyu Abe,
Wolfgang Brandner,
Timothy D. Brandt,
Joseph C. Carson,
Thayne Currie,
Sebastian Egner,
Markus Feldt,
Jeffrey Fung,
Miwa Goto,
Carol A. Grady,
Olivier Guyon,
Yutaka Hayano,
Masahiko Hayashi,
Saeko S. Hayashi
, et al. (34 additional authors not shown)
Abstract:
We present $H$- and $K_{\rm s}$-bands observations of the LkH$α$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investiga…
▽ More
We present $H$- and $K_{\rm s}$-bands observations of the LkH$α$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investigated the azimuthal profiles at the ring and the outer-disk regions as well as radial profiles in the directions of the spiral-like structures and semi-major axis. Azimuthal analysis shows a large variety in wavelength and implies that the disk has non-axisymmetric dust distributions. The radial profiles in the major-axis direction (PA=$271^\circ$) suggest that the outer region ($r\geq0\farcs25$) may be influenced by shadows of the inner region of the disk. The spiral-like directions (PA=10$^\circ$ and 230$^\circ$) show different radial profiles, which suggests that the surfaces of the spiral-like structures are highly flared and/or have different dust properties. Finally, a color-map of the disk shows a lack of an outer eastern region in the $H$-band disk, which may hint the presence of an inner object that casts a directional shadow onto the disk.
△ Less
Submitted 8 June, 2018; v1 submitted 16 April, 2018;
originally announced April 2018.
-
Inner Super-Earths, Outer Gas Giants: How Pebble Isolation and Migration Feedback Keep Jupiters Cold
Authors:
Jeffrey Fung,
Eve Lee
Abstract:
The majority of gas giants (planets of masses $\gtrsim10^2 M_\oplus$) are found to reside at distances beyond $\sim1$ au from their host stars. Within 1 au, the planetary population is dominated by super-Earths of $2-20 M_\oplus$. We show that this dichotomy between inner super-Earths and outer gas giants can be naturally explained should they form in nearly inviscid disks. In laminar disks, a pla…
▽ More
The majority of gas giants (planets of masses $\gtrsim10^2 M_\oplus$) are found to reside at distances beyond $\sim1$ au from their host stars. Within 1 au, the planetary population is dominated by super-Earths of $2-20 M_\oplus$. We show that this dichotomy between inner super-Earths and outer gas giants can be naturally explained should they form in nearly inviscid disks. In laminar disks, a planet can more easily repel disk gas away from its orbit. The feedback torque from the pile-up of gas inside the planet's orbit slows down and eventually halts migration. A pressure bump outside the planet's orbit traps pebbles and solids, starving the core. Gas giants are born cold and stay cold: more massive cores are preferentially formed at larger distances, and they barely migrate under disk feedback. We demonstrate this using 2D hydrodynamical simulations of disk-planet interaction lasting up to $10^5$ years: we track planet migration and pebble accretion until both come to an end by disk feedback. Whether cores undergo runaway gas accretion to become gas giants or not is determined by computing 1D gas accretion models. Our simulations show that in an inviscid minimum mass solar nebula, gas giants do not form inside $\sim$0.5 au, nor can they migrate there while the disk is present. We also explore the dependence on disk mass, and find that gas giants form further out in less massive disks.
△ Less
Submitted 16 May, 2018; v1 submitted 9 March, 2018;
originally announced March 2018.
-
The Orbit of the Companion to HD 100453A: Binary-Driven Spiral Arms in a Protoplanetary Disk
Authors:
Kevin Wagner,
Ruobing Dong,
Patrick Sheehan,
Daniel Apai,
Markus Kasper,
Melissa McClure,
Katie M. Morzinski,
Laird Close,
Jared Males,
Phil Hinz,
Sascha P. Quanz,
Jeffrey Fung
Abstract:
HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, whic…
▽ More
HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD 100453AB between 2015-2017 utilizing extreme adaptive optics systems on the Very Large Telescope and Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a=1.06"+/-0.09", e=0.17+/-0.07, and i=32.5+/- 6.5 degrees. We utilized publicly available ALMA CO data to constrain the inclination of the disk to i~28 degrees, which is relatively co-planar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamical and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-dirven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material.
△ Less
Submitted 11 January, 2018;
originally announced January 2018.
-
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.
-
Complex Spiral Structure in the HD 100546 Transitional Disk as Revealed by GPI and MagAO
Authors:
Katherine B. Follette,
Julien Rameau,
Ruobing Dong,
Laurent Pueyo,
Laird M. Close,
Gaspard Duchene,
Jeffrey Fung,
Clare Leonard,
Bruce Macintosh,
Jared R. Males,
Christian Marois,
Maxwell A. Millar-Blanchaer,
Katie M. Morzinski,
Wyatt Mullen,
Marshall Perrin,
Elijah Spiro,
Jason Wang,
S. Mark Ammons,
Vanessa P. Bailey,
Travis Barman,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
Robert J. De Rosa,
Rene Doyon
, et al. (31 additional authors not shown)
Abstract:
We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at Hα. The new GPI H -band total intensity data represent a significant enhanc…
▽ More
We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at Hα. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.
△ Less
Submitted 19 April, 2017;
originally announced April 2017.
-
Save the Planet, Feed the Star: How Super-Earths Survive Migration and Drive Disk Accretion
Authors:
Jeffrey Fung,
Eugene Chiang
Abstract:
Two longstanding problems in planet formation include (1) understanding how planets survive migration, and (2) articulating the process by which protoplanetary disks disperse---and in particular how they accrete onto their central stars. We can go a long way toward solving both problems if the disk gas surrounding planets has no intrinsic diffusivity ("viscosity"). In inviscid, laminar disks, a pl…
▽ More
Two longstanding problems in planet formation include (1) understanding how planets survive migration, and (2) articulating the process by which protoplanetary disks disperse---and in particular how they accrete onto their central stars. We can go a long way toward solving both problems if the disk gas surrounding planets has no intrinsic diffusivity ("viscosity"). In inviscid, laminar disks, a planet readily repels gas away from its orbit. On short timescales, zero viscosity gas accumulates inside a planet's orbit to slow Type I migration by orders of magnitude. On longer timescales, multiple super-Earths (distributed between, say, $\sim$0.1--10 AU) can torque inviscid gas out of interplanetary space, either inward to feed their stars, or outward to be blown away in a wind. We explore this picture with 2D hydrodynamics simulations of Earths and super-Earths embedded in inviscid disks, confirming their slow/stalled migration even under gas-rich conditions, and showing that disk transport rates range up to $\sim$$10^{-7} M_\odot~{\rm yr^{-1}}$ and scale as $\dot{M} \propto ΣM_{\rm p}^{3/2}$, where $Σ$ is the disk surface density and $M_{\rm p}$ is the planet mass. Gas initially sandwiched between two planets is torqued past both into the inner and outer disks. In sum, sufficiently compact systems of super-Earths can clear their natal disk gas, in a dispersal history that may be complicated and non-steady, but which conceivably leads over Myr timescales to large gas depletions similar to those characterizing transition disks.
△ Less
Submitted 24 March, 2017; v1 submitted 27 January, 2017;
originally announced January 2017.
-
What is the Mass of a Gap-Opening Planet?
Authors:
Ruobing Dong,
Jeffrey Fung
Abstract:
High contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their ma…
▽ More
High contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their masses, from the observed gap properties, such as their depths and widths. Combing 2D and 3D hydrodynamics simulations with 3D radiative transfer simulations, we investigate the morphology of planet-opened gaps in near-infrared scattered light images. Quantitatively, we obtain correlations that directly link intrinsic gap depths and widths in the gas surface density to observed depths and widths in images of disks at modest inclinations under finite angular resolution. Subsequently, the properties of the surface density gaps enable us to derive the disk scale height at the location of the gap $h$, and to constrain the quantity $M_{\rm p}^2/α$, where $M_{\rm p}$ is the mass of the gap-opening planet and $α$ characterizes the viscosity in the gap. As examples, we examine the gaps recently imaged by VLT/SPHERE, Gemini/GPI, and Subaru/HiCIAO in HD 97048, TW Hya, HD 169142, LkCa 15, and RX J1615.3-3255. Scale heights of the disks and possible masses of the gap-opening planets are derived assuming each gap is opened by a single planet. Assuming $α=10^{-3}$, the derived planet mass in all cases are roughly between 0.1-1 $M_{\rm J}$.
△ Less
Submitted 24 July, 2017; v1 submitted 14 December, 2016;
originally announced December 2016.
-
How Bright are Planet-Induced Spiral Arms in Scattered Light?
Authors:
Ruobing Dong,
Jeffrey Fung
Abstract:
Recently, high angular resolution imaging instruments such as SPHERE and GPI have discovered many spiral-arm-like features in near-infrared scattered light images of protoplanetary disks. Theory and simulations have suggested that these arms are most likely excited by planets forming in the disks; however, a quantitative relation between the arm-to-disk brightness contrast and planet mass is still…
▽ More
Recently, high angular resolution imaging instruments such as SPHERE and GPI have discovered many spiral-arm-like features in near-infrared scattered light images of protoplanetary disks. Theory and simulations have suggested that these arms are most likely excited by planets forming in the disks; however, a quantitative relation between the arm-to-disk brightness contrast and planet mass is still missing. Using 3D hydrodynamics and radiative transfer simulations, we examine the morphology and contrast of planet-induced arms in disks. We find a power-law relation for the face-on arm contrast ($δ_{\rm max}$) as a function of planet mass ($M_{\rm p}$) and disk aspect ratio ($h/r$): $δ_{\rm max}\approx((M_{\rm p}/M_{\rm J})/(h/r)^{1.38})^{0.22}$. With current observational capability, at a 30 AU separation, the minimum planet mass for driving detectable arms in a disk around a 1 Myr 1$M_\odot$ star at 140 pc at low inclinations is around Saturn mass. For planets more massive than Neptune masses, they typically drive multiple arms. Therefore in observed disks with spirals, it is unlikely that each spiral arm originates from a different planet. We also find only massive perturbers with at least multi-Jupiter masses are capable of driving bright arms with $δ_{\rm max}\gtrsim2$ as found in SAO 206462, MWC 758, and LkH$α$~330, and these arms do not follow linear wave propagation theory. Additionally, we find the morphology and contrast of the primary and secondary arms are largely unaffected by a modest level of viscosity with $α\lesssim0.01$. Finally, the contrast of the arms in the SAO 206462 disk suggests that the perturber SAO 206462 b at $\sim100$ AU is about $5-10M_{\rm J}$ in mass.
△ Less
Submitted 24 July, 2017; v1 submitted 1 December, 2016;
originally announced December 2016.
-
Planetary Torque in 3D Isentropic Disks
Authors:
Jeffrey Fung,
Frederic Masset,
Elena Lega,
David Velasco
Abstract:
Planet migration is inherently a three-dimensional (3D) problem, because Earth-size planetary cores are deeply embedded in protoplanetary disks. Simulations of these 3D disks remain challenging due to the steep requirement in resolution. Using two different hydrodynamics code, FARGO3D and PEnGUIn, we simulate disk-planet interaction for a 1 to 5 Earth-mass planet embedded in an isentropic disk. We…
▽ More
Planet migration is inherently a three-dimensional (3D) problem, because Earth-size planetary cores are deeply embedded in protoplanetary disks. Simulations of these 3D disks remain challenging due to the steep requirement in resolution. Using two different hydrodynamics code, FARGO3D and PEnGUIn, we simulate disk-planet interaction for a 1 to 5 Earth-mass planet embedded in an isentropic disk. We measure the torque on the planet and ensure that the measurements are converged both in resolution and between the two codes. We find that the torque is independent of the smoothing length of the planet's potential ($r_{\rm s}$), and that it has a weak dependence on the adiabatic index of the gaseous disk ($γ$). The torque values correspond to an inward migration rate qualitatively similar to previous linear calculations. We perform additional simulations with explicit radiative transfer using FARGOCA, and again find agreement between 3D simulations and existing torque formulae. We also present the flow pattern around the planets, and show that active flow is present within the planet's Hill sphere, and meridional vortices are shed downstream. The vertical flow speed near the planet is faster for a smaller $r_{\rm s}$ or $γ$, up to supersonic speeds for the smallest $r_{\rm s}$ and $γ$ in our study.
△ Less
Submitted 11 January, 2017; v1 submitted 28 October, 2016;
originally announced October 2016.
-
Gap Opening in 3D: Single Planet Gaps
Authors:
Jeffrey Fung,
Eugene Chiang
Abstract:
Giant planets can clear deep gaps when embedded in 2D (razor-thin) viscous circumstellar disks. We show by direct simulation that giant planets are just as capable of carving out gaps in 3D. Surface density maps are similar between 2D and 3D, even in detail. In particular, the scaling $Σ_{\rm gap} \propto q^{-2}$ of gap surface density with planet mass, derived from a global "zero-dimensional" bal…
▽ More
Giant planets can clear deep gaps when embedded in 2D (razor-thin) viscous circumstellar disks. We show by direct simulation that giant planets are just as capable of carving out gaps in 3D. Surface density maps are similar between 2D and 3D, even in detail. In particular, the scaling $Σ_{\rm gap} \propto q^{-2}$ of gap surface density with planet mass, derived from a global "zero-dimensional" balance of Lindblad and viscous torques, applies equally well to results obtained at higher dimensions. Our 3D simulations reveal extensive, near-sonic, meridional flows both inside and outside the gaps; these large-scale circulations might bear on disk compositional gradients, in dust or other chemical species. At high planet mass, gap edges are mildly Rayleigh unstable and intermittently shed streams of material into the gap - less so in 3D than in 2D.
△ Less
Submitted 19 September, 2016; v1 submitted 7 June, 2016;
originally announced June 2016.
-
How Spirals and Gaps Driven by Companions in Protoplanetary Disks Appear in Scattered Light at Arbitrary Viewing Angles
Authors:
Ruobing Dong,
Jeffrey Fung,
Eugene Chiang
Abstract:
Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a di…
▽ More
Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a disk by a companion varies with inclination and position angle in NIR scattered light. We compare the cases of a $3M_{\rm J}$ and a $0.1M_{\odot}$ companion, and make predictions suitable for testing with Gemini/GPI, VLT/NACO/SPHERE, and Subaru/HiCIAO/SCExAO. We find that the two trailing arms produced by an external perturber can have a variety of morphologies in inclined systems - they may appear as one trailing arm; two trailing arms on the same side of the disk; or two arms winding in opposite directions. The disk ring outside a planetary gap may also mimic spiral arms when viewed at high inclinations. We suggest potential explanations for the features observed in HH 30, HD 141569 A, AK Sco, HD 100546, and AB Aur. We emphasize that inclined views of companion-induced features cannot be converted into face-on views using simple and commonly practiced image deprojections.
△ Less
Submitted 5 May, 2016; v1 submitted 15 February, 2016;
originally announced February 2016.
-
An M Dwarf Companion and Its Induced Spiral Arms in the HD 100453 Protoplanetary Disk
Authors:
Ruobing Dong,
Zhaohuan Zhu,
Jeffrey Fung,
Roman Rafikov,
Eugene Chiang,
Kevin Wagner
Abstract:
Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to exami…
▽ More
Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.
△ Less
Submitted 15 December, 2015;
originally announced December 2015.
-
Inferring Planet Mass from Spiral Structures in Protoplanetary Disks
Authors:
Jeffrey Fung,
Ruobing Dong
Abstract:
Recent observations of protoplanetary disk have reported spiral structures that are potential signatures of embedded planets, and modeling efforts have shown that a single planet can excite multiple spiral arms, in contrast to conventional disk-planet interaction theory. Using two and three-dimensional hydrodynamics simulations to perform a systematic parameter survey, we confirm the existence of…
▽ More
Recent observations of protoplanetary disk have reported spiral structures that are potential signatures of embedded planets, and modeling efforts have shown that a single planet can excite multiple spiral arms, in contrast to conventional disk-planet interaction theory. Using two and three-dimensional hydrodynamics simulations to perform a systematic parameter survey, we confirm the existence of multiple spiral arms in disks with a single planet, and discover a scaling relation between the azimuthal separation of the primary and secondary arm, $φ_{\rm sep}$, and the planet-to-star mass ratio $q$: $φ_{\rm sep} = 102^{\circ} (q/0.001)^{0.2}$ for companions between Neptune mass and 16 Jupiter masses around a 1 solar mass star, and $φ_{\rm sep} = 180^{\circ}$ for brown dwarf mass companions. This relation is independent of the disk's temperature, and can be used to infer a planet's mass to within an accuracy of about 30% given only the morphology of a face-on disk. Combining hydrodynamics and Monte-Carlo radiative transfer calculations, we verify that our numerical measurements of $φ_{\rm sep}$ are accurate representations of what would be measured in near-infrared scattered light images, such as those expected to be taken by Gemini/GPI, VLT/SPHERE, or Subaru/SCExAO in the future. Finally, we are able to infer, using our scaling relation, that the planet responsible for the spiral structure in SAO 206462 has a mass of about 6 Jupiter masses.
△ Less
Submitted 9 December, 2015; v1 submitted 3 November, 2015;
originally announced November 2015.
-
The 3D Flow Field Around an Embedded Planet
Authors:
Jeffrey Fung,
Pawel Artymowicz,
Yanqin Wu
Abstract:
3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our multi-GPU hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar str…
▽ More
3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our multi-GPU hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet's Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the "transient" horseshoe flow. The flow continues as it rolls up into a pair of up-down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a $Σ\sim r^{-3/2}$ radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of 3 slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.
△ Less
Submitted 18 September, 2015; v1 submitted 12 May, 2015;
originally announced May 2015.
-
Irradiation Instability at the Inner Edges of Accretion Disks
Authors:
Jeffrey Fung,
Pawel Artymowicz
Abstract:
An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around AGNs. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also pr…
▽ More
An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around AGNs. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also present the linear growth rate as a function of beta, the ratio between radiation force and gravity, and c_s, the sound speed of the disk, obtained using two methods: a semi-analytic analysis of the linearized equations and a numerical simulation using the GPU-accelerated hydrodynamical code PEnGUIn. In particular, we find that IRI occurs at beta~0.1 if the transition region extends as wide as ~0.05r, and at higher beta values if it is wider. This threshold value applies to c_s ranging from 3% of the Keplerian orbital speed to 5%, and becomes higher if c_s is lower. Furthermore, in the nonlinear evolution of the instability, disks with a large beta and small c_s exhibit 'clumping', extreme local surface density enhancements that can reach over ten times the initial disk surface density.
△ Less
Submitted 1 July, 2014; v1 submitted 17 March, 2014;
originally announced March 2014.
-
How Empty are Disk Gaps Opened by Giant Planets?
Authors:
Jeffrey Fung,
Ji-Ming Shi,
Eugene Chiang
Abstract:
Gap clearing by giant planets has been proposed to explain the optically thin cavities observed in many protoplanetary disks. How much material remains in the gap determines not only how detectable young planets are in their birth environments, but also how strong corotation torques are, which impacts how planets can survive fast orbital migration. We determine numerically how the average surface…
▽ More
Gap clearing by giant planets has been proposed to explain the optically thin cavities observed in many protoplanetary disks. How much material remains in the gap determines not only how detectable young planets are in their birth environments, but also how strong corotation torques are, which impacts how planets can survive fast orbital migration. We determine numerically how the average surface density inside the gap, sigma_gap, depends on planet-to-star mass ratio q, Shakura-Sunyaev viscosity parameter alpha, and disk height-to-radius aspect ratio h/r. Our results are derived from our new GPU-accelerated Lagrangian hydrodynamical code PEnGUIn, and are verified by independent simulations with ZEUS90. For Jupiter-like planets, we find sigma_gap \propto q^-2.2 alpha^1.4 (h/r)^6.6, and for near brown dwarf masses, sigma_gap \propto q^-1 alpha^1.3 (h/r)^6.1. Surface density contrasts inside and outside gaps can be as large as 10^4, even when the planet does not accrete. We derive a simple analytic scaling, sigma_gap \propto q^-2 alpha^1 (h/r)^5, that compares reasonably well to empirical results, especially at low Neptune-like masses, and use discrepancies to highlight areas for progress.
△ Less
Submitted 29 December, 2013; v1 submitted 1 October, 2013;
originally announced October 2013.
-
Radiative transfer effects during primordial helium recombination
Authors:
Jens Chluba,
Jeffrey Fung,
Eric R. Switzer
Abstract:
In this paper we refine the calculation of primordial helium recombination, accounting for several additional effects that were neglected or treated more approximately in previous studies. These include consideration of (i) time-dependent radiative transfer interaction between the 2^1 P_1 - 1^1 S_0 and 2^3 P_1 - 1^1 S_0 resonances; (ii) time-dependent radiative transfer for the partially overlappi…
▽ More
In this paper we refine the calculation of primordial helium recombination, accounting for several additional effects that were neglected or treated more approximately in previous studies. These include consideration of (i) time-dependent radiative transfer interaction between the 2^1 P_1 - 1^1 S_0 and 2^3 P_1 - 1^1 S_0 resonances; (ii) time-dependent radiative transfer for the partially overlapping n^1 P_1 - 1^1 S_0, n^1 D_2 - 1^1 S_0 and n^3 P_1 - 1^1 S_0 series with 3 <= n <= 10; (iii) electron scattering within a kernel approach. We also briefly discuss the effect of electron scattering and HI quadrupole lines on the recombination of hydrogen. Although the physics of all the considered processes is interesting and subtle, for the standard cosmology, with Y_p \sim 0.24, the overall correction to the ionization history during helium recombination with respect to the previous implementation of CosmoRec remains smaller than |DeltaNe/Ne| \sim 0.05%. The dominant improvement is caused by consistent inclusion of resonance scattering for the 2^1 P_1 - 1^1 S_0 resonance. For cosmologies with a large helium fraction, Y_p \sim 0.4, the difference reaches |D N_e/N_e|\sim 0.22% at z\sim1800, however, the overall correction to the CMB power spectra is small, exceeding |D C_l/C_l| \sim 0.05 % only at l>\sim 3000. In comparison to the current version of Recfast the difference reaches |D C_l/C_l|\sim0.4% at l\sim 3000 for Y_p \sim 0.4, and also for the standard value Y_p \sim 0.24 we find differences |DC_l/C_l|>\sim 0.1% at l\sim 2500. The new processes are now included by the cosmological recombination code CosmoRec and can be activated as needed for most settings without affecting its runtime significantly.
△ Less
Submitted 10 May, 2012; v1 submitted 2 October, 2011;
originally announced October 2011.