skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The effects of disc self-gravity and radiative cooling on the formation of gaps and spirals by young planets
ABSTRACT We have carried out 2D hydrodynamical simulations to study the effects of disc self-gravity and radiative cooling on the formation of gaps and spirals. (1) With disc self-gravity included, we find stronger, more tightly wound spirals and deeper gaps in more massive discs. The deeper gaps are due to the larger Angular Momentum Flux (AMF) of the waves excited in more massive discs, as expected from the linear theory. The position of the secondary gap does not change, provided that the disc is not extremely massive (Q ≳ 2). (2) With radiative cooling included, the excited spirals become monotonically more open (less tightly wound) as the disc’s cooling time-scale increases. On the other hand, the amplitude and strength of the spirals decrease when the cooling time increases from a small value to ∼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 Tcool ∼ 1. Consequently, the induced primary gap is narrower and the secondary gap becomes significantly shallower when the cooling time becomes ∼1/Ω. When the secondary gap is present, the position of it moves to the inner disc from the fast cooling cases to the slow cooling cases. The dependence of gap properties on the cooling time-scale (e.g. in AS 209) provides a new way to constrain the disc optical depth and thus disc surface density.  more » « less
Award ID(s):
1753168
PAR ID:
10158146
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
493
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2287 to 2305
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract 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.1Mth ≤ Mp ≤ 1Mth, where Mth is the thermal mass), dimensionless stopping time (10−3 ≤ τs ≤ 1), and solid abundance (0 < Z ≤ 1). We find that when the dust particles are tightly coupled to the gas (τ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 ≲ τs ≲ 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 (Mp < Mth) 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. 
    more » « less
  2. ABSTRACT Spirals in protoplanetary discs have been used to locate the potential planet in discs. Since only the spiral shape from a circularly orbiting perturber is known, most previous works assume that the planet is in a circular orbit. We develop a simple semi-analytical method to calculate the shape of the spirals launched by an eccentric planet. We assume that the planet emits wavelets during its orbit, and the wave fronts of these propagating wavelets form the spirals. The resulting spiral shape from this simple method agrees with numerical simulations exceptionally well. The spirals excited by an eccentric planet can detach from the planet, bifurcate, or even cross each other, which are all reproduced by this simple method. The spiral’s bifurcation point corresponds to the wavelet that is emitted when the planet’s radial speed reaches the disc’s sound speed. Multiple spirals can be excited by an eccentric planet (more than five spirals when e ≳ 0.2). The pitch angle and pattern speed are different between different spirals and can vary significantly across one spiral. The spiral wakes launched by high-mass eccentric planets steepen to spiral shocks and the crossing of spiral shocks leads to distorted or broken spirals. With the same mass, a more eccentric planet launches weaker spirals and induces a shallower gap over a long period of time. The observed unusually large/small pitch angles of some spirals, the irregular multiple spirals, and the different pattern speeds between different spirals may suggest the existence of eccentric perturbers in protoplanetary discs. 
    more » « less
  3. ABSTRACT Stars form from the gravitational collapse of turbulent, magnetized molecular cloud cores. Our non-ideal MHD simulations reveal that the intrinsically anisotropic magnetic resistance to gravity during the core collapse naturally generates dense gravomagneto sheetlets within inner protostellar envelopes – disrupted versions of classical sheet-like pseudo-discs. They are embedded in a magnetically dominant background, where less dense materials flow along the local magnetic field lines and accumulate in the dense sheetlets. The sheetlets, which feed the disc predominantly through its upper and lower surfaces, are the primary channels for mass and angular momentum transfer from the envelope to the disc. The protostellar disc inherits a small fraction (up to 10 per cent) of the magnetic flux from the envelope, resulting in a disc-averaged net vertical field strength of 1–10 mG and a somewhat stronger toroidal field, potentially detectable through ALMA Zeeman observations. The inherited magnetic field from the envelope plays a dominant role in disc angular momentum evolution, enabling the formation of gravitationally stable discs in cases where the disc field is relatively well-coupled to the gas. Its influence remains significant even in marginally gravitationally unstable discs formed in the more magnetically diffusive cases, removing angular momentum at a rate comparable to or greater than that caused by spiral arms. The magnetically driven disc evolution is consistent with the apparent scarcity of prominent spirals capable of driving rapid accretion in deeply embedded protostellar discs. The dense gravomagneto sheetlets observed in our simulations may correspond to the ‘accretion streamers’ increasingly detected around protostars. 
    more » « less
  4. null (Ed.)
    ABSTRACT Recent ALMA molecular line observations have revealed 3D gas velocity structure in protoplanetary discs, shedding light on mechanisms of disc accretion and structure formation. (1) By carrying out viscous simulations, we confirm that the disc’s velocity structure differs dramatically using vertical stress profiles from different accretion mechanisms. Thus, kinematic observations tracing flows at different disc heights can potentially distinguish different accretion mechanisms. On the other hand, the disc surface density evolution is mostly determined by the vertically integrated stress. The sharp disc outer edge constrained by recent kinematic observations can be caused by a radially varying α in the disc. (2) We also study kinematic signatures of a young planet by carrying out 3D planet–disc simulations. The relationship between the planet mass and the ‘kink’ velocity is derived, showing a linear relationship with little dependence on disc viscosity, but some dependence on disc height when the planet is massive (e.g. 10MJ). We predict the ‘kink’ velocities for the potential planets in DSHARP discs. At the gap edge, the azimuthally averaged velocities at different disc heights deviate from the Keplerian velocity at similar amplitudes, and its relationship with the planet mass is consistent with that in 2D simulations. After removing the planet, the azimuthally averaged velocity barely changes within the viscous time-scale, 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 ‘kink’ velocity is needed to probe young planets in protoplanetary discs. 
    more » « less
  5. ABSTRACT Recent high angular resolution ALMA observations have revealed numerous gaps in protoplanetary discs. A popular interpretation has been that planets open them. Most previous investigations of planet gap-opening have concentrated on viscous discs. Here, we carry out 2D (axisymmetric) global simulations of gap opening by a planet in a wind-launching non-ideal MHD disc with consistent thermochemistry. We find a strong concentration of poloidal magnetic flux in the planet-opened gap, where the gas dynamics are magnetically dominated. The magnetic field also drives a fast (nearly sonic) meridional gas circulation in the denser disc regions near the inner and outer edges of the gap, which may be observable through high-resolution molecular line observations. The gap is more ionized than its denser surrounding regions, with a better magnetic field–matter coupling. In particular, it has a much higher abundance of molecular ion HCO+, consistent with ALMA observations of the well-studied AS 209 protoplanetary disc that has prominent gaps and fast meridional motions reaching the local sound speed. Finally, we provide fitting formulae for the ambipolar and Ohmic diffusivities as a function of the disc local density, which can be used for future 3D simulations of planet gap-opening in non-ideal MHD discs where thermochemistry is too computationally expensive to evolve self-consistently with the magneto-hydrodynamics. 
    more » « less