BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Jeffrey Fung

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.COastro-ph.GAastro-ph.SRphysics.comp-ph

Catalog footprint

What is connected

12works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

12 published item(s)

preprint2022arXiv

Dust Dynamics in Transitional Disks: Clumping and Disk Recession

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Visualizing the Kinematics of Planet Formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Circumplanetary Disk Dynamics in the Isothermal and Adiabatic Limits

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gap Opening in 3D: Single Planet Gaps

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inferring Planet Mass from Spiral Structures in Protoplanetary Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The 3D Flow Field Around an Embedded Planet

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Irradiation Instability at the Inner Edges of Accretion Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How Empty are Disk Gaps Opened by Giant Planets?

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative transfer effects during primordial helium recombination

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.

0 citations0 reviewsNo code linkedOpen access