We use three dimensional simulations with coupled hydrodynamics and Monte Carlo radiative transfer to show that shadows cast by the inner disc in broken circumbinary discs move within a confined range of position angles on the outer disc. Over time, shadows appear to rock back and forth in azimuth as the inner disc precesses. The effect occurs because the inner disc precesses around a vector that is not the angular momentum vector of the outer disc. We relate our findings to recent observations of shadows in discs.
Context. Despite the recent discovery of spiral-shaped features in protoplanetary discs in the near-infrared and millimetric wavelengths, there is still an active discussion to understand how they formed. In fact, the spiral waves observed in discs around young stars can be due to different physical mechanisms: planet/companion torques, gravitational perturbations or illumination effects. Aims. We study the spirals formed in the gaseous phase due to two diametrically opposed shadows cast at fixed disc locations. The shadows are created by an inclined non-precessing disc inside the cavity, which is assumed to be optically thick. In particular, we analyse the effect of these spirals on the dynamics of the dust particles and discuss their detectability in transition discs. Methods. We perform gaseous hydrodynamical simulations with shadows, then we compute the dust evolution on top of the gaseous distribution, and finally we produce synthetic ALMA observations of the dust emission based on radiative transfer calculations. Results. Our main finding is that mm- to cm-sized dust particles are efficiently trapped inside the shadow-triggered spirals. We also observe that particles of various sizes starting at different stellocentric distances are well mixed inside these pressure maxima. This dynamical effect would favour grain growth and affect the resulting composition of planetesimals in the disc. In addition, our radiative transfer calculations show spiral patterns in the disc at 1.6 {mu}m and 1.3 mm. Due to their faint thermal emission (compared to the bright inner regions of the disc) the spirals cannot be detected with ALMA. Our synthetic observations prove however that shadows are observable as dips in the thermal emission.
In recent years an increasing number of observational studies have hinted at the presence of warps in protoplanetary discs, however a general comprehensive description of observational diagnostics of warped discs was missing. We performed a series of 3D SPH hydrodynamic simulations and combined them with 3D radiative transfer calculations to study the observability of warps in circumbinary discs, whose plane is misaligned with respect to the orbital plane of the central binary. Our numerical hydrodynamic simulations confirm previous analytical results on the dependence of the warp structure on the viscosity and the initial misalignment between the binary and the disc. To study the observational signatures of warps we calculate images in the continuum at near-infrared and sub-millimetre wavelengths and in the pure rotational transition of CO in the sub-millimetre. Warped circumbinary discs show surface brightness asymmetry in near-infrared scattered light images as well as in optically thick gas lines at sub-millimetre wavelengths. The asymmetry is caused by self-shadowing of the disc by the inner warped regions, thus the strength of the asymmetry depends on the strength of the warp. The projected velocity field, derived from line observations, shows characteristic deviations, twists and a change in the slope of the rotation curve, from that of an unperturbed disc. In extreme cases even the direction of rotation appears to change in the disc inwards of a characteristic radius. The strength of the kinematical signatures of warps decreases with increasing inclination. The strength of all warp signatures decreases with decreasing viscosity.
Three-dimensional hydrodynamic numerical simulations have demonstrated that the structure of a protoplanetary disc may be strongly affected by a planet orbiting in a plane that is misaligned to the disc. When the planet is able to open a gap, the disc is separated into an inner, precessing disc and an outer disc with a warp. In this work, we compute infrared scattered light images to investigate the observational consequences of such an arrangement. We find that an inner disc misaligned by a less than a degree to the outer disc is indeed able to cast a shadow at larger radii. In our simulations a planet of around 6 Jupiter masses inclined by around 2 degrees is enough to warp the disc and cast a shadow with a depth of more than 10% of the average flux at that radius. We also demonstrate that warp in the outer disc can cause a variation in the azimuthal brightness profile at large radii. Importantly, this latter effect is a function of the distance from the star and is most prominent in the outer disc. We apply our model to the TW Hya system, where a misaligned, precessing inner disc has been invoked to explain an recently observed shadow in the outer disc. Consideration of the observational constraints suggest that an inner disc precessing due to a misaligned planet is an unlikely explanation for the features found in TW Hya.
Shadows and spirals seem to be common features of transition discs. Among the spiral-triggering mechanisms proposed, only one establishes a causal link between shadows and spirals so far. In fact, provided the presence of shadows in the disc, the combined effect of temperature gradient and differential disc rotation, creates strong azimuthal pressure gradients. After several thousand years, grand-design spirals develop in the gas phase. Previous works have only considered static shadows caused by an inclined inner disc. However, in some cases, the inner regions of circumbinary discs can break and precess. Thus, it is more realistic to consider moving shadow patterns in the disc. In this configuration, the intersection between the inner and the outer discs defines the line of nodes at which the shadows are cast. Here, we consider moving shadows and study the resulting circumbinary disc structure. We find that only static and prograde shadows trigger spirals, in contrast to retrograde ones. Interestingly, if a region of the disc corotates with the shadow, a planet-like signature develops at the co-rotation position. The resulting spirals resemble those caused by a planet embedded in the disc, with similar pitch angles.
As of today ten circumbinary planets orbiting solar type main sequence stars have been discovered. Nearly all of them orbit around the central binary very closely to the region of instability where it is difficult to form them in situ. It is assumed that they formed further out and migrated to their observed position. We extend previous studies to a more realistic thermal disc structure and determine what parameter influence the final parking location of a planet around a binary star. We perform two-dimensional numerical simulations of viscous accretion discs around a central binary that include viscous heating and radiative cooling from the disc surfaces. We vary the binary eccentricity as well as disc viscosity and mass. Concerning the disc evolution we find that it can take well over 100000 binary orbits until an equilibrium state is reached. As seen previously, we find that the central cavity opened by the binary becomes eccentric and precesses slowly in a prograde sense. Embedded planets migrate to the inner edge of the disc. In cases of lower disc viscosity they migrate further in maintaining a circular orbit, while for high viscosity they are parked further out on an eccentric orbit. The final location of an embedded planet is linked to its ability to open a gap in the disc. Gap opening planets separate inner from outer disc, preventing eccentricity excitation in the latter and making it more circular. This allows embedded planets to migrate closer to the binary, in agreement with the observations. The necessary condition for gap opening and the final planet position depend on the planet mass and disc viscosity.