Gravitational coupling between young planets and their parent disks is often explored using numerical simulations, which typically treat the disk thermodynamics in a highly simplified manner. In particular, many studies adopt the locally isothermal approximation, in which the disk temperature is a fixed function of the stellocentric distance. We explore the dynamics of planet-driven density waves in disks with more general thermodynamics, in which the temperature is relaxed towards an equilibrium profile on a finite cooling timescale $t_{\rm c}$. We use both linear perturbation theory and direct numerical simulations to examine the global structure of density waves launched by planets in such disks. A key diagnostic used in this study is the behavior of the wave angular momentum flux (AMF), which directly determines the evolution of the underlying disk. The AMF of free waves is constant for slowly cooling (adiabatic) disks, but scales with the disk temperature for rapidly cooling (and locally isothermal) disks. However, cooling must be extremely fast, with $\beta = \Omega t_{\rm c} \lesssim 10^{-3}$ for the locally isothermal approximation to provide a good description of density wave dynamics in the linear regime (relaxing to $\beta \lesssim 10^{-2}$ when nonlinear effects are important). For intermediate cooling timescales, density waves are subject to a strong linear damping. This modifies the appearance of planet-driven spiral arms and the characteristics of axisymmetric structures produced by massive planets: in disks with $\beta \approx 0.1$ -- $1$, a near-thermal mass planet opens only a single wide gap around its orbit, in contrast to the several narrow gaps produced when cooling is either faster or slower.

High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings. When the ring edges are relatively sharp and the dust surface density becomes comparable to the gas surface density, we find that dust feedback enhances the radial gradients of both the azimuthal velocity profile and the potential vorticity profile at the ring edges. This eventually leads to instabilities on meso-scales (spatial scales of several disk scale heights), causing dusty rings to be populated with many compact regions with highly concentrated dust densities on meso-scales. We also produce synthetic dust emission images using our simulation results and discuss the comparison between simulations and observations.

We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite'', steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite'' disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2\Omega_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $\Omega_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($\Omega_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue'' $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2\Omega_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite'' disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.

It has been recently suggested that the multiple concentric rings and gaps discovered by ALMA in many protoplanetary disks may be produced by a single planet, as a result of the complex propagation and dissipation of the multiple spiral density waves it excites in the disk. Numerical efforts to verify this idea have largely utilized the so-called locally isothermal approximation with a prescribed disk temperature profile. However, in protoplanetary disks this approximation does not provide an accurate description of the density wave dynamics on scales of tens of au. Moreover, we show that locally isothermal simulations tend to overestimate the contrast of ring and gap features, as well as misrepresent their positions, when compared to simulations in which the energy equation is evolved explicitly. This outcome is caused by the non-conservation of the angular momentum flux of linear perturbations in locally isothermal disks. We demonstrate this effect using simulations of locally isothermal and adiabatic disks (with essentially identical temperature profiles) and show how the dust distributions, probed by mm wavelength observations, differ between the two cases. Locally isothermal simulations may thus underestimate the masses of planets responsible for the formation of multiple gaps and rings on scales of tens of au observed by ALMA. We suggest that caution should be exercised in using the locally isothermal simulations to explore planet-disk interaction, as well as in other studies of wave-like phenomena in astrophysical disks.

The quiet Sun exhibits a wealth of magnetic activities that are fundamental for our understanding of solar and astrophysical magnetism. The magnetic fields in the quiet Sun are observed to evolve coherently, interacting with each other to form distinguished structures as they are advected by the horizontal photospheric flows. We study coherent structures in photospheric flows in a region of quiet Sun consisted of supergranules. Supergranular turbulence is investigated by detecting hyperbolic and elliptic Lagrangian coherent structures (LCS) using the horizontal velocity fields derived from Hinode intensity maps. Repelling/attracting LCS are found by computing the forward/backward finite-time Lyapunov exponent (FTLE). The Lagrangian centre of a supergranular cell is given by the local maximum of the forward FTLE; the Lagrangian boundaries of supergranular cells are given by the ridges of the backward FTLE. Objective velocity vortices are found by calculating the Lagrangian-averaged vorticity deviation, and false vortices are filtered by applying a criterion given by the displacement vector. The Lagrangian centres of neighboring supergranular cells are interconnected by ridges of the repelling LCS, which provide the transport barriers that allow the formation of vortices and the concentration of strong magnetic fields in the valleys of the repelling LCS. The repelling LCS also reveal the most likely sites for magnetic reconnection. We show that the ridges of the attracting LCS expose the locations of the sinks of photospheric flows at supergranular junctions, which are the preferential sites for the formation of kG magnetic flux tubes and persistent vortices.

Recent observations of protoplanetary disks, as well as simulations of planet-disk interaction, have suggested that a single planet may excite multiple spiral arms in the disk, in contrast to the previous expectations based on linear theory (predicting a one-armed density wave). We re-assess the origin of multiple arms in the framework of linear theory, by solving for the global two-dimensional response of a non-barotropic disk to an orbiting planet. We show that the formation of a secondary arm in the inner disk, at about half of the orbital radius of the planet, is a robust prediction of linear theory. This arm becomes stronger than the primary spiral at several tenths of the orbital radius of the planet. Several additional, weaker spiral arms may also form in the inner disk. On the contrary, a secondary spiral arm is unlikely to form in the outer disk. Our linear calculations, fully accounting for the global behavior of both the phases and amplitudes of perturbations, generally support the recently proposed WKB phase argument for the secondary arm origin (as caused by the intricacy of constructive interference of azimuthal harmonics of the perturbation at different radii). We provide analytical arguments showing that the process of a single spiral wake splitting up into multiple arms is a generic linear outcome of wave propagation in differentially rotating disks. It is not unique to planet-driven waves and occurs also in linear calculations of spiral wakes freely propagating with no external torques. These results are relevant for understanding formation of multiple rings and gaps in protoplanetary disks.

We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the AREPO code. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-mass binaries with arbitrary eccentricities. We evolve the flow over long (viscous) timescales until a quasi-steady is reached, in which the mass supply rate at large distances $\dot{M}_0$ (assumed constant) equals the time-averaged mass transfer rate across the disk and the total mass accretion rate onto the binary components. This quasi-steady state allows us to compute the secular angular momentum transfer rate onto the binary, $\langle\dot{J}_b\rangle$, and the resulting orbital evolution. Through direct computation of the gravitational and accretion torques on the binary, we find that $\langle\dot{J}_b\rangle$ is consistently positive (i.e., the binary gains angular momentum), with $l_0\equiv\langle\dot{J}_b\rangle/\dot M_0$ in the range of $(0.4-0.8)a_b^2\Omega_b$, depending on the binary eccentricity (where $a_b,~\Omega_b$ are the binary semi-major axis and angular frequency); we also find that this $\langle\dot{J}_b\rangle$ is equal to the net angular momentum current across the CBD, indicating that global angular momentum balance is achieved in our simulations. We compute the time-averaged rate of change of the binary orbital energy for eccentric binaries, and thus obtain the secular rates $\langle\dot a_b\rangle$ and $\langle \dot{e}_b\rangle$. In all cases, $\langle\dot{a}_b\rangle$ is positive, i.e., the binary expands while accreting. We discuss the implications of our results for the merger of supermassive binary black holes and for the formation of close stellar binaries.

Spectroscopic observations of some metal-rich white dwarfs (WDs), believed to be polluted by planetary material, reveal the presence of compact gaseous metallic disks orbiting them. The observed variability of asymmetric, double-peaked emission line profiles in about half of such systems could be interpreted as the signature of precession of an eccentric gaseous debris disk. The variability timescales --- from decades down to $1.4$ yr (recently inferred for the debris disk around HE 1349--2305) --- are in rough agreement with the rate of general relativistic (GR) precession in the test particle limit. However, it has not been demonstrated that this mechanism can drive such a fast, coherent precession of a radially extended (out to $1 R_\odot$) gaseous disk mediated by internal stresses (pressure). Here we use the linear theory of eccentricity evolution in hydrodynamic disks to determine several key properties of eccentric modes in gaseous debris disks around WDs. We find a critical dependence of both the precession period and radial eccentricity distribution of the modes on the inner disk radius, $r_\mathrm{in}$. For small inner radii, $r_\mathrm{in} \lesssim (0.2 - 0.4) R_\odot$, the modes are GR-driven, with periods of $\approx 1 - 10$ yr. For $r_\mathrm{in} \gtrsim (0.2 - 0.4) R_\odot$, the modes are pressure-dominated, with periods of $\approx 3 - 20$ yr. Correspondence between the variability periods and inferred inner radii of the observed disks is in general agreement with this trend. In particular, the short period of HE 1349--2305 is consistent with its small $r_\mathrm{in}$. Circum-WD debris disks may thus serve as natural laboratories for studying the evolution of eccentric gaseous disks.

We investigate the migration of a low-mass ($\lesssim 10 M_\oplus$) planet near the inner edge of a protoplanetary disc using two-dimensional viscous hydrodynamics simulations. We employ an inner boundary condition representing the truncation of the disc at the stellar corotation radius. As described by Tsang (2011), wave reflection at the inner disc boundary modifies the Type I migration torque on the planet, allowing migration to be halted before the planet reaches the inner edge of the disc. For low-viscosity discs ($\alpha \lesssim 10^{-3}$), planets may be trapped with semi-major axes as large as $3-5$ times the inner disc radius. In general, planets are trapped closer to the inner edge as either the planet mass or the disc viscosity parameter $\alpha$ increases, and farther from the inner edge as the disc thickness is increased. This planet trapping mechanism may impact the formation and migration history of close-in compact multiplanet systems.

A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after $10^4$ orbits, or $2.5$ Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

We carry out numerical simulations of circumbinary discs, solving the viscous hydrodynamics equations on a polar grid covering an extended disc outside the binary co-orbital region. We use carefully controlled outer boundary conditions and long-term integrations to ensure that the disc reaches a quasi-steady state, in which the time-averaged mass accretion rate onto the binary, $\langle\dot{M}\rangle$, matches the mass supply rate at the outer disc. We focus on binaries with comparable masses and a wide range of eccentricities ($e_\mathrm{B}$). For $e_\mathrm{B} \lesssim 0.05$, the mass accretion rate of the binary is modulated at about $5$ times the binary period; otherwise it is modulated at the binary period. The inner part of the circumbinary disc ($r \lesssim 6 a_\mathrm{B}$) generally becomes coherently eccentric. For low and high $e_\mathrm{B}$, the disc line of apsides precesses around the binary, but for intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), it instead becomes locked with that of the binary. By considering the balance of angular momentum transport through the disc by advection, viscous stress, and gravitational torque, we determine the time-averaged net angular momentum transfer rate to the binary, $\langle\dot{J}\rangle$. The specific angular momentum, $l_0 = \langle\dot{J}\rangle/\langle\dot{M}\rangle$, depends non-monotonically on $e_\mathrm{B}$. Contrary to previous claims, we find that $l_0$ is positive for most $e_\mathrm{B}$, implying that the binary receives net angular momentum, which may cause its separation to grow with time. The minimum $l_0$ occurs at intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), corresponding to the regime where the inner eccentric disc is apsidally aligned with the binary.

The physical mechanism of angular momentum transport in poorly ionized regions of protoplanetary discs, the dead zones (DZs), is not understood. The presence of a DZ naturally leads to conditions susceptible to the Rossby wave instability (RWI), which produces vortices and spiral density waves that may revive the DZ and be responsible for observed large-scale disc structures. We present a series of two-dimensional hydrodynamic simulations to investigate the role of the RWI in DZs, including its impact on the long-term evolution of the disc and its morphology. The nonlinear RWI can generate Reynolds stresses (effective $\alpha$ parameter) as large as $0.01 - 0.05$ in the DZ, helping to sustain quasi-steady accretion throughout the disc. It also produces novel disc morphologies, including azimuthal asymmetries with $m = 1, 2$, and atypical vortex shapes. The angular momentum transport strength and morphology are most sensitive to two parameters: the radial extent of the DZ and the disc viscosity. The largest Reynolds stresses are produced when the radial extent of the DZ is less than its distance to the central star. Such narrow DZs lead to a single vortex or two coherent antipodal vortices in the quasi-steady state. The edges of wider DZs evolve separately, resulting in two independent vortices and reduced angular momentum transport efficiency. In either case, we find that, because of the Reynolds stresses generated by the nonlinear RWI, gravitational instability is unlikely to play a role in angular momentum transport across the DZ, unless the accretion rate is sufficiently high.

Recent observations have shown that circumstellar and circumbinary discs in young stellar binaries are often misaligned with respect to the binary orbital plane. We analyze the tidal truncation of such misaligned discs due to torques applied to the disc at the Lindblad resonances from the tidal forcings of the binary. We consider eccentric binaries with arbitrary binary-disc inclination angles. We determine the dependence of the tidal forcing strengths on the binary parameters and show that they are complicated non-monotonic functions of eccentricity and inclination. We adopt a truncation criterion determined by the balance between resonant torque and viscous torque, and use it to calculate the outer radii of circumstellar discs and the inner radii of circumbinary discs. Misaligned circumstellar discs have systematically larger outer radii than aligned discs, and are likely to fill their Roche lobes if inclined by more than $45^\circ - 90^\circ$, depending on the binary mass ratio and disc viscosity parameter. Misaligned circumbinary discs generally have smaller inner radii than aligned discs, but the details depend sensitively on the binary and disc parameters.

We examine the role played by viscosity in the excitation of global oscillation modes (both axisymmetric and non-axisymmetric) in accretion discs around black holes using two-dimensional hydrodynamic simulations. The turbulent viscosity is modeled by the $\alpha$-ansatz, with different equations of state. We consider both discs with transonic radial inflows across the innermost stable circular orbit, and stationary discs truncated by a reflecting wall at their inner edge, representing a magnetosphere. In transonic discs, viscosity can excite several types of global oscillation modes. These modes are either axisymmetric with frequencies close to multiples of the maximum radial epicyclic frequency $\kappa_\mathrm{max}$, non-axisymmetric with frequencies close to multiples of of the innermost stable orbit frequency $\Omega_\mathrm{ISCO}$, or hybrid modes whose frequencies are linear combinations of these two frequencies. Small values of the viscosity parameter $\alpha$ primarily produce non-axisymmetric modes, while axisymmetric modes become dominant for large $\alpha$. The excitation of these modes may be related to an instability of the sonic point, at which the radial infall speed is equal to the sound speed of the gas. In discs with a reflective inner boundary, we explore the effect of viscosity on trapped $p$-modes which are intrinsically overstable due to the corotation resonance effect. The effect of viscosity is either to reduce the growth rates of these modes, or to completely suppress them and excite a new class of higher frequency modes. The latter requires that the dynamic viscosity scales positively with the disc surface density, indicating that it is a result of the classic viscous overstability effect.

We study the development of coherent structures in local simulations of the magnetorotational instability in accretion discs in regimes of on-off intermittency. In a previous paper [Chian et al., Phys. Rev. Lett. 104, 254102 (2010)], we have shown that the laminar and bursty states due to the on-off spatiotemporal intermittency in a one-dimensional model of nonlinear waves correspond, respectively, to nonattracting coherent structures with higher and lower degrees of amplitude-phase synchronization. In this paper we extend these results to a three-dimensional model of magnetized Keplerian shear flows. Keeping the kinetic Reynolds number and the magnetic Prandtl number fixed, we investigate two different intermittent regimes by varying the plasma beta parameter. The first regime is characterized by turbulent patterns interrupted by the recurrent emergence of a large-scale coherent structure known as two-channel flow, where the state of the system can be described by a single Fourier mode. The second regime is dominated by the turbulence with sporadic emergence of coherent structures with shapes that are reminiscent of a perturbed channel flow. By computing the Fourier power and phase spectral entropies in three-dimensions, we show that the large-scale coherent structures are characterized by a high degree of amplitude-phase synchronization.