No Arabic abstract
Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa).
We present 3D smoothed particle hydrodynamics simulations of protoplanetary discs undergoing a flyby by a stellar perturber on a parabolic orbit lying in a plane inclined relative to the disc mid-plane. We model the disc as a mixture of gas and dust, with grains ranging from 1 {mu}m to 10 cm in size. Exploring different orbital inclinations, periastron distances and mass ratios, we investigate the disc dynamical response during and after the flyby. We find that flybys induce evolving spiral structure in both gas and dust which can persist for thousands of years after periastron. Gas and dust structures induced by the flyby differ because of drag-induced effects on the dust grains. Variations in the accretion rate by up to an order of magnitude occur over a time-scale of order 10 years or less, inducing FU Orionis-like outbursts. The remnant discs are truncated and warped. The dust disc is left more compact than the gas disc, both because of disc truncation and accelerated radial drift of grains induced by the flyby.
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.
Misalignments between the orbital planes of planets and the equatorial planes of their host stars have been observed in our solar system, in transiting exoplanets, and in the orbital planes of debris discs. We present a mechanism that causes such a spin-orbit misalignment for a protoplanetary disc due to its movement through an ambient medium. Our physical explanation of the mechanism is based on the theoretical solutions to the Stark problem. We test this idea by performing self-consistent hydrodynamical simulations and simplified gravitational $N$-body simulations. The $N$-body model reduces the mechanism to the relevant physical processes. The hydrodynamical simulations show the mechanism in its full extent, including gas-dynamical and viscous processes in the disc which are not included in the theoretical framework. We find that a protoplanetary disc embedded in a flow changes its orientation as its angular momentum vector tends to align parallel to the relative velocity vector. Due to the force exerted by the flow, orbits in the disc become eccentric, which produces a net torque and consequentially changes the orbital inclination. The tilting of the disc causes it to contract. Apart from becoming lopsided, the gaseous disc also forms a spiral arm even if the inclination does not change substantially. The process is most effective at high velocities and observational signatures are therefore mostly expected in massive star-forming regions and around winds or supernova ejecta. Our $N$-body model indicates that the interaction with supernova ejecta is a viable explanation for the observed spin-orbit misalignment in our solar system.
We study the three-dimensional evolution of a viscous protoplanetary disc which is perturbed by a passing star on a parabolic orbit. The aim is to test whether a single stellar flyby is capable to excite significant disc inclinations which would favour the formation of so-called misaligned planets. We use smoothed particle hydrodynamics to study inclination, disc mass and angular momentum changes of the disc for passing stars with different masses. We explore different orbital configurations for the perturbers orbit to find the parameter spaces which allow significant disc inclination generation. Prograde inclined parabolic orbits are most destructive leading to significant disc mass and angular momentum loss. In the remaining disc, the final disc inclination is only below $20^circ$. This is due to the removal of disc particles which have experienced the strongest perturbing effects. Retrograde inclined parabolic orbits are less destructive and can generate disc inclinations up to $60^circ$. The final disc orientation is determined by the precession of the disc angular momentum vector about the perturbers orbital angular momentum vector and by disc orbital inclination changes. We propose a sequence of stellar flybys for the generation of misalignment angles above $60^circ$. The results taken together show that stellar flybys are promising and realistic for the explanation of misaligned Hot Jupiters with misalignment angles up to 60degr.
A key problem in protoplanetary disc evolution is understanding the efficiency of dust radial drift. This process makes the observed dust disc sizes shrink on relatively short timescales, implying that discs started much larger than what we see now. In this paper we use an independent constraint, the gas radius (as probed by CO rotational emission), to test disc evolution models. In particular, we consider the ratio between the dust and gas radius, $R_{rm CO}/R_{rm dust}$. We model the time evolution of protoplanetary discs under the influence of viscous evolution, grain growth, and radial drift. Then, using the radiative transfer code RADMC with approximate chemistry, we compute the dust and gas radii of the models and investigate how $R_{rm CO}/R_{rm dust}$ evolves. Our main finding is that, for a broad range of values of disc mass, initial radius, and viscosity, $R_{rm CO}/R_{rm dust}$ becomes large (>5) after only a short time (<1 Myr) due to radial drift. This is at odds with measurements in young star forming regions such as Lupus, which find much smaller values, implying that dust radial drift is too efficient in these models. Substructures, commonly invoked to stop radial drift in large, bright discs, must then be present, although currently unresolved, in most discs.