BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Zhaohuan Zhu

Zhaohuan Zhu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph.IMastro-ph.HEMachine Learningphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
30works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

30 published item(s)

preprint2025arXiv

Brightness variability in polar circumbinary disks

In binary systems with a strongly misaligned disk, the central binary stars can travel a significant vertical distance above and below the disk's orbital plane. This can cause large changes in illumination of the disk over the course of the binary orbital period. We use both analytic and radiative transfer models to examine the effect of changes in stellar illumination on the appearance of the disk, particularly in the case of the polar disk HD 98800B. We find that the observed flux from the disk can vary significantly over the binary orbital period, producing a periodically varying lightcurve which peaks twice each binary orbit. The amount of flux variation is strongly influenced by the disk geometry. We suggest that these flux variations produce several observable signatures, and that these observables may provide constraints on different properties of the disk such as its vertical structure, geometry, and cooling rate.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning

We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the collisional growth of micron to mm-sized dust particles in early stage planet formation. The simulation data are used to train a U-Net deep learning model to predict gridded three-dimensional representations of the particle density and velocity fields, given as input the corresponding fluid fields. The trained model qualitatively captures the filamentary structure of clustered particles in a highly non-linear regime. We assess model fidelity by calculating metrics of the density field (the radial distribution function) and of the velocity field (the relative velocity and the relative radial velocity between particles). Although trained only on the spatial fields, the model predicts these statistical quantities with errors that are typically <10%. Our results suggest that, given appropriately expanded training data, deep learning could complement direct numerical simulations in predicting particle clustering within turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Structured Distributions of Gas and Solids in Protoplanetary Disks

Recent spatially-resolved observations of protoplanetary disks revealed a plethora of substructures, including concentric rings and gaps, inner cavities, misalignments, spiral arms, and azimuthal asymmetries. This is the major breakthrough in studies of protoplanetary disks since Protostars and Planets VI and is reshaping the field of planet formation. However, while the capability of imaging substructures in protoplanetary disks has been steadily improving, the origin of many substructures are still largely debated. The structured distributions of gas and solids in protoplanetary disks likely reflect the outcome of physical processes at work, including the formation of planets. Yet, the diverse properties among the observed protoplanetary disk population, for example, the number and radial location of rings and gaps in the dust distribution, suggest that the controlling process may differ between disks and/or the outcome may be sensitive to stellar or disk properties. In this review, we (1) summarize the existing observations of protoplanetary disk substructures collected from the literature; (2) provide a comprehensive theoretical review of various processes proposed to explain observed protoplanetary disk substructures; (3) compare current theoretical predictions with existing observations and highlight future research directions to distinguish between different origins; and (4) discuss implications of state-of-the-art protoplanetary disk observations to protoplanetary disk and planet formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

ALMA Discovery of a Disk around the Planetary-mass Companion SR 12 c

We report an Atacama Large Millimeter/submillimeter Array 0.88 mm (Band 7) continuum detection of the accretion disk around SR 12 c, an $\sim$11 $M_{\rm Jup}$ planetary-mass companion (PMC) orbiting its host binary at 980 au. This is the first submillimeter detection of a circumplanetary disk around a wide PMC. The disk has a flux density of $127 \pm14~μ$Jy and is not resolved by the $\sim$0.1&#34; beam, so the dust disk radius is likely less than 5 au and can be much smaller if the dust continuum is optically thick. If, however, the dust emission is optically thin, then the SR 12 c disk has a comparable dust mass to the circumplanetary disk around PDS 70 c but is about five times lower than that of the $\sim$12 $M_{\rm Jup}$ free-floating OTS 44. This suggests that disks around bound and unbound planetary-mass objects can span a wide range of masses. The gas mass estimated with an accretion rate of $10^{-11}~M_\odot$ yr$^{-1}$ implies a gas-to-dust ratio higher than 100. If cloud absorption is not significant, a nondetection of ${}^{12}$CO(3-2) implies a compact gas disk around SR 12 c. Future sensitive observations may detect more PMC disks at 0.88 mm flux densities of $\lesssim$100 $μ$Jy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Direct Formation of Planetary Embryos in Self-Gravitating Disks

Giant planets have been discovered at large separations from the central star. Moreover, a striking number of young circumstellar disks have gas and/or dust gaps at large orbital separations, potentially driven by embedded planetary objects. To form massive planets at large orbital separations through core accretion within disk lifetime, however, an early solid body to seed pebble and gas accretion is desirable. Young protoplanetary disks are likely self-gravitating, and these gravitoturbulent disks may efficiently concentrate solid material at the midplane driven by spiral waves. We run 3D local hydrodynamical simulations of gravitoturbulent disks with Lagrangian dust particles to determine whether particle and gas self-gravity can lead to the formation of dense solid bodies, seeding later planet formation. When self-gravity between dust particles is included, solids of size $\mathrm{St} = 0.1$ to $1$ concentrate within the gravitoturbulent spiral features and collapse under their own self-gravity into dense clumps up to several $M_{\oplus}$ in mass at wide orbits. Simulations with dust that drift most efficiently, $\mathrm{St}=1$, form the most massive clouds of particles, while simulations with smaller dust particles, $\mathrm{St}=0.1$, have clumps with masses an order of magnitude lower. When the effect of dust backreaction onto the gas is included, dust clumps become smaller by a factor of a few but more numerous. The existence of large solid bodies at an early stage of the disk can accelerate the planet formation process, particularly at wide orbital separations, and potentially explain planets distant from the central stars and young protoplanetary disks with substructures.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows

Rings and gaps are commonly observed in the dust continuum emission of young stellar disks. Previous studies have shown that substructures naturally develop in the weakly ionized gas of magnetized, non-ideal MHD disks. The gas rings are expected to trap large mm/cm-sized grains through pressure gradient-induced radial dust-gas drift. Using 2D (axisymmetric) MHD simulations that include ambipolar diffusion and dust grains of three representative sizes (1~mm, 3.3~mm, and 1~cm), we show that the grains indeed tend to drift radially relative to the gas towards the centers of the gas rings, at speeds much higher than in a smooth disk because of steeper pressure gradients. However, their spatial distribution is primarily controlled by meridional gas motions, which are typically much faster than the dust-gas drift. In particular, the grains that have settled near the midplane are carried rapidly inwards by a fast accretion stream to the inner edges of the gas rings, where they are lifted up by the gas flows diverted away from the midplane by a strong poloidal magnetic field. The flow pattern in our simulation provides an attractive explanation for the meridional flows recently inferred in HD 163296 and other disks, including both &#34;collapsing&#34; regions where the gas near the disk surface converges towards the midplane and a disk wind. Our study highlights the prevalence of the potentially observable meridional flows associated with the gas substructure formation in non-ideal MHD disks and their crucial role in generating rings and gaps in dust.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, &#34;famous&#34; disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty &#34;FS CMa&#34; systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

HST/WFC3 H$α$ Direct-Imaging Detection of a Point-like Source in the Disk Cavity of AB Aur

Accreting protoplanets enable the direct characterization of planet formation. As part of a high-contrast imaging search for accreting planets with the Hubble Space Telescope (HST) Wide Field Camera 3, we present H$α$ images of AB Aurigae (AB Aur), a Herbig Ae/Be star harboring a transition disk. The data were collected in two epochs of direct-imaging observations using the F656N narrow-band filter. After subtracting the point spread function of the primary star, we identify a point-like source located at a P.A. of $182.5^{\circ}\pm1.4^{\circ}$ and a separation of $600\pm22$~mas relative to the host star. The position is consistent with the recently identified protoplanet candidate AB Aur b. The source is visible in two individual epochs separated by ${\sim}50$ days and the H$α$ intensities in the two epochs agree. The H$α$ flux density is $F_ν=1.5\pm0.4$~mJy, $3.2\pm0.9$ times of the optical continuum determined by published HST/STIS photometry. In comparison to PDS 70 b and c, the H$α$ excess emission is weak. The central star is accreting and the stellar H$α$ emission has a similar line-to-continuum ratio as seen in AB Aur b. We conclude that both planetary accretion and scattered stellar light are possible sources of the H$α$ emission, and the H$α$ detection alone does not validate AB Aur b as an accreting protoplanet. Disentangling the origin of the emission will be crucial for probing planet formation in the AB Aur disk.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Investigating the Future Potential of an Upgraded ALMA to Image Planet Forming Disks at Sub-au Scales

In recent years, ALMA has been able to observe large-scale substructures within protoplanetary disks. Comparison with the predictions from models of planet-disk interaction has indicated that most of these disk substructures can be explained by the presence of planets with the mass of Neptune or larger at orbital radii of $\approx 5 - 100$ au. Better resolution is needed to observe structures closer to the star, where terrestrial planets are expected to form, as well as structures opened by planets with masses lower than Neptune. We investigate the capabilities of a possible extension to ALMA that would double the longest baseline lengths in the array to detect and resolve disk substructures opened by Earth-mass and Super Earth planets at orbital radii of $1-5$ au. By simulating observations of a family of disk models using this extended configuration in ALMA Bands 6 and 7, we show that an upgraded ALMA would detect gaps in disks formed by super-Earths as close as 1 au, as well as Earth-mass planets down to $2-3$ au from the young host stars in nearby star forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Maximum Temperatures in Evolving Protoplanetary Discs and Composition of Planetary Building Blocks

The maximum temperature and radial temperature profile in a protoplanetary disc are important for the condensation of different elements in the disc. We simulate the evolution of a set of protoplanetary discs from the collapse of their progenitor molecular cloud cores as well as the dust decoupling within the discs as they evolve. We show how the initial properties of the cloud cores affect the thermal history of the protoplanetary discs using a simple viscous disc model. Our results show that the maximum midplane temperature in the disc occurs within 0.5 AU. It increases with the initial cloud temperature and decreases with its angular velocity and the viscosity of the disc. From the observed properties of the molecular cloud cores we find the median value of the maximum temperature is around 1250 K, with roughly 90% of them being less than 1500 K - a value that is lower than the 50% condensation temperatures of most refractory elements. Therefore, only cloud cores with high initial temperatures or low angular velocities and/or low viscosities within the planet-forming discs will result in refractory-rich planetesimals. To reproduce the volatile depletion pattern of CM, CO, and CV chondrites and the terrestrial planets in Solar system, one must either have rare properties of the initial molecular cloud cores like high core temperature, or other sources of energy to heat the disc to sufficiently high temperatures. Alternatively, the volatile depletion observed in these chondrites may be inherited from the progenitor molecular cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves

Local thermal instability can plausibly explain the formation of multiphase gas in many different astrophysical environments, but the theory is only well understood in the optically thin limit of the equations of radiation hydrodynamics (RHD). Here we lay groundwork for transitioning from this limit to a full RHD treatment assuming a gray opacity formalism. We consider a situation where the gas becomes thermally unstable due to the hardening of the radiation field when the main radiative processes are free-free cooling and Compton heating. We identify two ways in which this can happen: (i) when the Compton temperature increases with time, through a rise in either the intensity or energy of a hard X-ray component; and (ii) when attenuation reduces the flux of the thermal component so that the Compton temperature increases with depth through the slab. Both ways likely occur in the broad line region of active galactic nuclei where columns of gas can be ionization bounded. In such instances where attenuation is significant, thermal equilibrium solution curves become position-dependent and it no longer suffices to assess the stability of an irradiated column of gas at all depths using a single equilibrium curve. We demonstrate how to analyze a new equilibrium curve -- the attenuation curve -- for this purpose, and we show that by Field&#39;s instability criterion, a negative slope along this curve indicates that constant density slabs are thermally unstable whenever the gas temperature increases with depth.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Constraining Protoplanetary Disk Accretion and Young Planets Using ALMA Kinematic Observations

Recent ALMA molecular line observations have revealed 3-D gas velocity structure in protoplanetary disks, shedding light on mechanisms of disk accretion and structure formation. 1) By carrying out viscous simulations, we confirm that the disk&#39;s velocity structure differs dramatically using vertical stress profiles from different accretion mechanisms. Thus, kinematic observations tracing flows at different disk heights can potentially distinguish different accretion mechanisms. On the other hand, the disk surface density evolution is mostly determined by the vertically integrated stress. The sharp disk outer edge constrained by recent kinematic observations can be caused by a radially varying $α$ in the disk. 2) We also study kinematic signatures of a young planet by carrying out 3-D planet-disk simulations. The relationship between the planet mass and the &#34;kink&#34; velocity is derived, showing a linear relationship with little dependence on disk viscosity, but some dependence on disk height when the planet is massive, e.g. $10 M_J$. We predict the &#34;kink&#34; velocities for the potential planets in DSHARP disks. At the gap edge, the azimuthally-averaged velocities at different disk heights deviate from the Keplerian velocity at similar amplitudes, and its relationship with the planet mass is consistent with that in 2-D simulations. After removing the planet, the azimuthally-averaged velocity barely changes within the viscous timescale, and thus the azimuthally-averaged velocity structure at the gap edge is due to the gap itself and not directly caused to the planet. Combining both axisymmetric kinematic observations and the residual &#34;kink&#34; velocity is needed to probe young planets in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Kozai-Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star-planet system. We use three-dimensional hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of $H/r\gtrsim 0.05$, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai--Lidov (KL) oscillations. Dissipation in the Kozai--Lidov unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt--KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The ALMA Early Science view of FUor/EXor objects. I. Through the looking-glass of V2775 Ori

As part of an ALMA survey to study the origin of episodic accretion in young eruptive variables, we have observed the circumstellar environment of the star V2775 Ori. This object is a very young, pre-main sequence object which displays a large amplitude outburst characteristic of the FUor class. We present Cycle-2 band 6 observations of V2775 Ori with a continuum and CO (2-1) isotopologue resolution of 0.25\as (103 au). We report the detection of a marginally resolved circumstellar disc in the ALMA continuum with an integrated flux of $106 \pm 2$ mJy, characteristic radius of $\sim$ 30 au, inclination of $14.0^{+7.8}_{-14.5}$ deg, and is oriented nearly face-on with respect to the plane of the sky. The \co~emission is separated into distinct blue and red-shifted regions that appear to be rings or shells of expanding material from quasi-episodic outbursts. The system is oriented in such a way that the disc is seen through the outflow remnant of V2775 Ori, which has an axis along our line-of-sight. The $^{13}$CO emission displays similar structure to that of the \co, while the C$^{18}$O line emission is very weak. We calculated the expansion velocities of the low- and medium-density material with respect to the disc to be of -2.85 km s$^{-1}$ (blue), 4.4 km s$^{-1}$ (red) and -1.35 and 1.15 km s$^{-1}$ (for blue and red) and we derived the mass, momentum and kinetic energy of the expanding gas. The outflow has an hourglass shape where the cavities are not seen. We interpret the shapes that the gas traces as cavities excavated by an ancient outflow. We report a detection of line emission from the circumstellar disc and derive a lower limit of the gas mass of 3 \MJup.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A fast-growing tilt instability of detached circumplanetary disks

Accretion disks in binary systems can exhibit a tilt instability, arising from the interaction between components of the tidal potential and dissipation. Using a linear analysis, we show that the aspect ratios and outer radii of circumplanetary disks provide favorable conditions for tilt growth. We quantify the growth rate of the instability using particle-based ({\sc phantom}) and grid-based ({\sc athena++}) hydrodynamic simulations. For a disk with outer aspect ratio $H/r \simeq 0.1$, initially moderate tilts double on a time scale of about 15-30 binary orbits. Our results imply that detached circumplanetary disks, whose evolution is not entirely controlled by accretion from the circumstellar disk, may commonly be misaligned to the planetary orbital plane. We discuss implications for planetary spin evolution, and possible interactions between the tilt instability and Kozai-Lidov dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced &#34;shoulder&#34; near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $α$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur&#39;s inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as &#34;transitional.&#34;

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Study of Millimeter Variability in FUor Objects

FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object&#39;s jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ALMA 0.88 mm Survey of Disks around Planetary-mass Companions

Characterizing the physical properties and compositions of circumplanetary disks can provide important insights into the formation of giant planets and satellites. We report ALMA 0.88 mm (Band 7) continuum observations of six planetary-mass (10-20 $M_{\rm Jup}$) companions: CT Cha b, 1RXS 1609 b, ROXs 12 b, ROXs 42B b, DH Tau b, and FU Tau b. No continuum sources are detected at the locations of the companions down to 3$σ$ limits of 120-210 $μ$Jy. Given these non-detections, it is not clear whether disks around planetary-mass companions indeed follow the disk flux-host mass trend in the stellar regime. The faint radio brightness of these companion disks may result from a combination of fast radial drift and a lack of dust traps. Alternatively, as disks in binary systems are known to have significantly lower millimeter fluxes due to tidal interactions, these companion disks may instead follow the relationship of moderate-separation binary stars. This scenario can be tested with sensitive continuum imaging at rms levels of $\lesssim$10 $μ$Jy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dust Condensation in Evolving Discs and the Composition of Planetary Building Blocks

Partial condensation of dust from the Solar nebula is likely responsible for the diverse chemical compositions of chondrites and rocky planets/planetesimals in the inner Solar system. We present a forward physical-chemical model of a protoplanetary disc to predict the chemical compositions of planetary building blocks that may form from such a disc. Our model includes the physical evolution of the disc and the condensation, partial advection, and decoupling of the dust within it. The chemical composition of the condensate changes with time and radius. We compare the results of two dust condensation models: one where an element condenses when the midplane temperature in the disc is lower than the 50\% condensation temperature ($\rm T_{50}$) of that element and the other where the condensation of the dust is calculated by a Gibbs free energy minimization technique assuming chemical equilibrium at local disc temperature and pressure. The results of two models are generally consistent with some systematic differences of $\sim 10$\% depending upon the radial distance and an element&#39;s condensation temperature. Both models predict compositions similar to CM, CO, and CV chondrites provided that the decoupling timescale of the dust is on the order of the evolution timescale of the disc or longer. If the decoupling timescale is too short, the composition deviates significantly from the measured values. These models may contribute to our understanding of the chemical compositions of chondrites, and ultimately the terrestrial planets in the solar system, and may constrain the potential chemical compositions of rocky exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori&#39;s Accretion Disk and Observational Signatures of Magnetic Fields

FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori&#39;s accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori&#39;s inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the &#34;surface accretion&#34; mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup

While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolution of $\sim0.3&#39;&#39;$ ($\sim50$ au). The CO isotopologue emission traces four major components of the RU Lup system: a compact Keplerian disk with a radius of $\sim120$ au, a non-Keplerian ``envelope-like&#39;&#39; structure surrounding the disk and extending to $\sim260$ au from the star, at least five blueshifted spiral arms stretching up to 1000 au, and clumps outside the spiral arms located up to 1500 au in projection from RU Lup. We comment on potential explanations for RU Lup&#39;s peculiar gas morphology, including gravitational instability, accretion of material onto the disk, or perturbation by another star. RU Lup&#39;s extended non-Keplerian CO emission, elevated stellar accretion rate, and unusual photometric variability suggest that it could be a scaled-down Class II analog of the outbursting FU Ori systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Near-Infrared Imaging of a Spiral in the CQ Tau Disk

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and &#34;double gap&#34; features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The effects of disk self-gravity and radiative cooling on the formation of gaps and spirals by young planets

We have carried out two-dimensional hydrodynamical simulations to study the effects of disk self-gravity and radiative cooling on the formation of gaps and spirals. (1) With disk self-gravity included, we find stronger, more tightly-wound spirals and deeper gaps in more massive disks. The deeper gaps are due to the larger Angular Momentum Flux (AMF) of the waves excited in more massive disks, as expected from the linear theory. The position of the secondary gap does not change, provided that the disk is not extremely massive ($Q \gtrsim 2$). (2) With radiative cooling included, the excited spirals become monotonically more open (less tightly-wound) as the disk&#39;s cooling timescale increases. On the other hand, the amplitude and strength of the spirals decrease when the cooling time increases from a small value to $\sim 1/Ω$, but then the amplitude starts to increase again when the cooling time continues to increase. This indicates that radiative dissipation becomes important for waves with $T_{cool}\sim$ 1. Consequently, the induced primary gap is narrower and the secondary gap becomes significantly shallower when the cooling time becomes $\sim 1/Ω$. When the secondary gap is present, the position of it moves to the inner disk from the fast cooling cases to the slow cooling cases. The dependence of gap properties on the cooling timescale (e.g. in AS 209) provides a new way to constrain the disk optical depth and thus disk surface density.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Planetary Luminosity Problem: &#34;Missing Planets&#34; and the Observational Consequences of Episodic Accretion

The high occurrence rates of spiral arms and large central clearings in protoplanetary disks, if interpreted as signposts of giant planets, indicate that gas giants form commonly as companions to young stars ($<$ few Myr) at orbital separations of 10--300\,au. However, attempts to directly image this giant planet population as companions to more mature stars ($> 10$\, Myr) have yielded few successes. This discrepancy could be explained if most giant planets form &#34;cold start,&#34; i.e., by radiating away much of their formation energy as they assemble their mass, rendering them faint enough to elude detection at later times. In that case, giant planets should be bright at early times, during their accretion phase, and yet forming planets are detected only rarely through direct imaging techniques. Here we explore the possibility that the low detection rate of accreting planets is the result of episodic accretion through a circumplanetary disk. We also explore the possibility that the companion orbiting the Herbig Ae star HD~142527 may be a giant planet undergoing such an accretion outburst.

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&#39;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 &#34;disks&#34;), 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
preprint2019arXiv

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08&#39;&#39; and 0.11&#39;&#39; respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2&#39;&#39; of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects&#39; circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Morphological signatures induced by dust back reaction in discs with an embedded planet

Recent observations have revealed a gallery of substructures in the dust component of nearby protoplanetary discs, including rings, gaps, spiral arms, and lopsided concentrations. One interpretation of these substructures is the existence of embedded planets. Not until recently, however, most of the modelling effort to interpret these observations ignored the dust back reaction to the gas. In this work, we conduct local-shearing-sheet simulations for an isothermal, inviscid, non-self-gravitating, razor-thin dusty disc with a planet on a fixed circular orbit. We systematically examine the parameter space spanned by planet mass ($0.1M_\text{th} \le M_\text{p} \le 1M_\text{th}$, where $M_\text{th}$ is the thermal mass), dimensionless stopping time ($10^{-3} \le τ_\text{s} \le 1$), and solid abundance ($0 < Z \le 1$). We find that when the dust particles are tightly coupled to the gas ($τ_\text{s} < 0.1$), the spiral arms are less open and the gap driven by the planet becomes deeper with increasing $Z$, consistent with a reduced speed of sound in the approximation of a single dust-gas mixture. By contrast, when the dust particles are marginally coupled ($0.1 \lesssim τ_\text{s} \lesssim 1$), the spiral structure is insensitive to $Z$ and the gap structure in the gas can become significantly skewed and unidentifiable. When the latter occurs, the pressure maximum radially outside of the planet is weakened or even extinguished, and hence dust filtration by a low-mass ($M_\text{p} < M_\text{th}$) planet could be reduced or eliminated. Finally, we find that the gap edges where the dust particles are accumulated as well as the lopsided large-scale vortices driven by a massive planet, if any, are unstable, and they are broken into numerous small-scale dust-gas vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Resolving the FU Ori System with ALMA: Interacting Twin Disks?

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6&#39;&#39;, i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively.

0 citations0 reviewsNo code linkedOpen access