Spiral arms in protoplanetary discs are thought to be linked to the presence of companions. We test the hypothesis that the double spiral arm morphology observed in the transition disc MWC 758 can be generated by an $\approx 10$ M$_{\rm Jup}$ companion on an eccentric orbit internal to the spiral arms. Previous studies on MWC 758 have assumed an external companion. We compare simulated observations from three dimensional hydrodynamics simulations of disc-companion interaction to scattered light, infrared and CO molecular line observations, taking into account observational biases. The inner companion hypothesis is found to explain the double spiral arms, as well as several additional features seen in MWC 758 -- the arc in the northwest, substructures inside the spiral arms, the cavity in CO isotopologues, and the twist in the kinematics. Testable predictions include detection of fainter spiral structure, detection of a point source south-southeast of the primary, and proper motion of the spiral arms.

In this work we demonstrate that the inner spiral structure observed in AB Aurigae can be created by a binary star orbiting inside the dust cavity. We find that a companion with a mass-ratio of 0.25, semi-major axis of 40 au, eccentricity of 0.5, and inclination of 90{\deg} produces gaseous spirals closely matching the ones observed in $^{12}$CO (2-1) line emission. Based on dust dynamics in circumbinary discs (Poblete, Cuello, and Cuadra 2019), we constrain the inclination of the binary with respect to the circumbinary disc to range between 60{\deg} and 90{\deg}. We predict that the stellar companion is located roughly 0.18 arcsec from the central star towards the east-southeast, above the plane of the disc. Should this companion be detected in the near future, our model indicates that it should be moving away from the primary star at a rate of 6 mas/yr on the plane of the sky. Since our companion is inclined, we also predict that the spiral structure will appear to change with time, and not simply co-rotate with the companion.

We present observations of the young multiple system UX Tauri to look for circumstellar disks and for signs of dynamical interactions. We obtained SPHERE/IRDIS deep differential polarization images in the J and H bands. We also used ALMA archival CO data. Large extended spirals are well detected in scattered light coming out of the disk of UX Tau A. The southern spiral forms a bridge between UX Tau A and C. These spirals, including the bridge connecting the two stars, all have a CO (3-2) counterpart seen by ALMA. The disk of UX Tau C is detected in scattered light. It is much smaller than the disk of UX Tau A and has a major axis along a different position angle, suggesting a misalignment. We performed PHANTOM SPH hydrodynamical models to interpret the data. The scattered light spirals, CO emission spirals and velocity patterns of the rotating disks, and the compactness of the disk of UX Tau C all point to a scenario in which UX Tau A has been perturbed very recently (about 1000 years) by the close passage of UX Tau C.

Close-in giant planets represent the most significant evidence of planetary migration. If large exomoons form around migrating giant planets which are more stable (e.g. those in the Solar System), what happens to these moons after migration is still under intense research. This paper explores the scenario where large regular exomoons escape after tidal-interchange of angular momentum with its parent planet, becoming small planets by themselves. We name this hypothetical type of object a \textit{ploonet}. By performing semi-analytical simulations of tidal interactions between a large moon with a close-in giant, and integrating numerically their orbits for several Myr, we found that in $\sim$50 per cent of the cases a young ploonet may survive ejection from the planetary system, or collision with its parent planet and host star, being in principle detectable. Volatile-rich ploonets are dramatically affected by stellar radiation during both planetocentric and siderocentric orbital evolution, and their radius and mass change significantly due to the sublimation of most of their material during time-scales of hundred of Myr. We estimate the photometric signatures that ploonets may produce if they transit the star during the phase of evaporation, and compare them with noisy lightcurves of known objects ("Kronian" stars and non-periodical dips in dusty lightcurves). Additionally, the typical transit timing variations (TTV) induced by the interaction of a ploonet with its planet are computed. We find that present and future photometric surveys' capabilities can detect these effects and distinguish them from those produced by other nearby planetary encounters.

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'', 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.

We now have several observational examples of misaligned broken protoplanetary discs, where the disc inner regions are strongly misaligned with respect to the outer disc. Current models suggest that this disc structure can be generated with an internal misaligned companion (stellar or planetary), but the occurrence rate of these currently unobserved companions remains unknown. Here we explore whether a strong misalignment between the inner and outer disc can be formed without such a companion. We consider a disc that has an existing gap --- essentially separating the disc into two regions --- and use a flyby to disturb the discs, leading to a misalignment. Despite considering the most optimistic parameters for this scenario, we find maximum misalignments between the inner and outer disc of $\sim$45$^{\circ}$ and that these misalignments are short-lived. We thus conclude that the currently observed misaligned discs must harbour internal, misaligned companions.

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 test the hypothesis that the disc cavity in the `transition disc' Oph IRS 48 is carved by an unseen binary companion. We use 3D dust-gas smoothed-particle hydrodynamics simulations to demonstrate that marginally coupled dust grains concentrate in the gas over-density that forms in in the cavity around a low binary mass ratio binary. This produces high contrast ratio dust asymmetries at the cavity edge similar to those observed in the disc around IRS 48 and other transition discs. This structure was previously assumed to be a vortex. However, we show that the observed velocity map of IRS 48 displays a peculiar asymmetry that is not predicted by the vortex hypothesis. We show the unusual kinematics are naturally explained by the non-Keplerian flow of gas in an eccentric circumbinary cavity. We further show that perturbations observed in the isovelocity curves of IRS 48 may be explained as the product of the dynamical interaction between the companion and the disc. The presence of a $\sim$0.4 M$_{\odot}$ companion at a $\sim$10 au separation can qualitatively explain these observations. High spatial resolution line and continuum imaging should be able to confirm this hypothesis.

Recent observations have revealed that protoplanetary discs often exhibit cavities and azimuthal asymmetries such as dust traps and clumps. The presence of a stellar binary system in the inner disc regions has been proposed to explain the formation of these structures. Here, we study the dust and gas dynamics in circumbinary discs around eccentric and inclined binaries. This is done through two-fluid simulations of circumbinary discs, considering different values of the binary eccentricity and inclination. We find that two kinds of dust structures can form in the disc: a single horseshoe-shaped clump, on top of a similar gaseous over-density; or numerous clumps, distributed along the inner disc rim. The latter features form through the complex interplay between the dust particles and the gaseous spirals caused by the binary. All these clumps survive between one and several tens of orbital periods at the feature location. We show that their evolution strongly depends on the gas-dust coupling and the binary parameters. Interestingly, these asymmetric features could in principle be used to infer or constrain the orbital parameters of a stellar companion - potentially unseen - inside the inner disc cavity. Finally, we apply our findings to the disc around AB Aurigae.

Context. Dynamical studies suggest that most circumbinary discs (CBDs) should be coplanar. However, some theoretical works show that under certain conditions a CBD can become polar, which means that its rotation vector is orthogonal with respect to the binary orbital plane. Interestingly, very recent observations show that polar CBDs exist in nature (e.g. HD 98800). Aims. We test the predictions of CBD alignment around eccentric binaries based on linear theory. In particular, we compare prograde and retrograde CBD configurations. Then, assuming planets form in these systems, we thoroughly characterise the orbital behaviour and stability of misaligned (P-type) particles. This is done for massless and massive particles. Methods. The evolution of the CBD alignment for various configurations was modelled through three-dimensional hydrodynamical simulations. For the orbital characterisation and the analysis stability, we relied on long-term N-body integrations and structure and chaos indicators, such as $\Delta e$ and MEGNO. Results. We confirm previous analytical predictions on CBD alignment, but find an unexpected symmetry breaking between prograde and retrograde configurations. Therefore, the likelihood of becoming polar for a highly misaligned CBD is higher than previously thought. Regarding the stability of circumbinary P-type planets (also know as Tatooines), polar orbits are stable over a wide range of binary parameters. In particular, for binary eccentricities below 0.4 the orbits are stable for any value of the binary mass ratio. In the absence of gas, planets with masses below $10^{-5}\,M_{\odot}$ have negligible effects on the binary orbit. Finally, we suggest that mildly eccentric equal-mass binaries should be searched for polar Tatooines.

Dynamical studies suggest that most of the circumbinary discs (CBDs) should be coplanar. However, under certain initial conditions, the CBD can evolve toward polar orientation. Here we extend the parametric study of polar configurations around detached close-in binaries through $N$-body simulations. For polar configurations around binaries with mass ratios $q$ below $0.7$, the nominal location of the mean motion resonance (MMR) $1~:~4$ predicts the limit of stability for $e_{B} > 0.1$. Alternatively, for $e_{B} < 0.1$ or $q \sim 1$, the nominal location of the MMR $1~:~3$ is the closest stable region. The presence of a} giant planet increases the region of forbidden polar configurations around low mass ratio binaries with eccentricities $e_B\sim0.4$ with respect to rocky earth-like planets. For equal mass stars, the eccentricity excitation $\Delta e$ of polar orbits smoothly increases with decreasing distance to the binary. For $q<1$, $\Delta e$ can reach values as high as $0.4$. Finally, we studied polar configurations around $HD~98800BaBb$ and show that the region of stability is strongly affected by the relative positions of the nodes. The most stable configurations in the system correspond to polar particles, which are not expected to survive on longer time-scales due to the presence of the external perturber HD~$98800AaAb$.

As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

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.

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.

We present a methodology to interpret observations of protoplanetary discs where a flyby, also called a tidal encounter, is suspected. In case of a flyby, protoplanetary discs can be significantly disturbed. The resulting dynamical and kinematical signatures can last for several thousands of years after the flyby and hence deeply affect the evolution of the disc. These effects are stronger for closer encounters and more massive perturbers. For the very same flyby parameters, varying the inclination of the perturber's orbit produces a broad range of disc structures: spirals, bridges, warps and cavities. We study this kind of features both in the gas and in the dust for grains ranging from 1 {\mu}m to 10 cm in size. Interestingly, the dust exhibits a different dynamical behaviour compared to the gas because of gas-drag effects. Finally, flybys can also trigger high accretion events in the disc-hosting star, readily similar to FU Orionis-type outbursts. All this information can be used to infer the flyby parameters from an incomplete set of observations at different wavelengths. Therefore, the main scope of our flyby scene investigation (FSI) methodology is to help to interpret recent "puzzling" disc observations.

We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions.

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.

The size and density of dust grains determine their response to gas drag in protoplanetary discs. Aerodynamical (size x density) sorting is one of the proposed mechanisms to explain the grain properties and chemical fractionation of chondrites. However, the efficiency of aerodynamical sorting and the location in the disc in which it could occur are still unknown. Although the effects of grain sizes and growth in discs have been widely studied, a simultaneous analysis including dust composition is missing. In this work we present the dynamical evolution and growth of multicomponent dust in a protoplanetary disc using a 3D, two-fluid (gas+dust) Smoothed Particle Hydrodynamics (SPH) code. We find that the dust vertical settling is characterised by two phases: a density-driven phase which leads to a vertical chemical sorting of dust and a size-driven phase which enhances the amount of lighter material in the midplane. We also see an efficient radial chemical sorting of the dust at large scales. We find that dust particles are aerodynamically sorted in the inner disc. The disc becomes sub-solar in its Fe/Si ratio on the surface since the early stage of evolution but sub-solar Fe/Si can be also found in the outer disc-midplane at late stages. Aggregates in the disc mimic the physical and chemical properties of chondrites, suggesting that aerodynamical sorting played an important role in determining their final structure.

Calcium Aluminium Inclusions (CAIs) formed in the Solar System, some 4,567 million years ago. CAIs are almost always surrounded by Wark-Lovering Rims (WLRs), which are a sequence of thin, mono/bi-mineralic layers of refractory minerals, with a total thickness in the range of 1 to 100 microns. Recently, some CAIs have been found that have tektite-like bowl-shapes. To form such shapes, the CAI must have travelled through a rarefied gas at hypersonic speeds. We show how CAIs may have been ejected from the inner solar accretion disc via the centrifugal interaction between the solar magnetosphere and the inner disc rim. They subsequently punched through the hot, inner disc rim wall at hypersonic speeds. This re-entry heating partially or completely evaporated the CAIs. Such evaporation could have significantly increased the metal abundances of the inner disc rim. High speed movement through the inner disc produced WLRs. To match the observed thickness of WLRs required metal abundances at the inner disc wall that are of order ten times that of standard solar abundances. The CAIs cooled as they moved away from the protosun, the deduced CAI cooling rates are consistent with the CAI cooling rates obtained from experiment and observation. The speeds and gas densities required to form bowl-shaped CAIs are also consistent with the expected speeds and gas densities for larger, ~ 1 cm, CAIs punching through an inner accretion disc wall.

Photophoresis is a physical process based on momentum exchange between an illuminated dust particle and its gaseous environment. Its net effect in protoplanetary discs (PPD) is the outward transport of solid bodies from hot to cold regions. This process naturally leads to the formation of ring-shaped features where dust piles up. In this work, we study the dynamical effects of photophoresis in PPD by including the photophoretic force in the two-fluid (gas+dust) smoothed particle hydrodynamics (SPH) code developed by Barri\`ere-Fouchet et al. (2005). We find that the conditions of pressure and temperature encountered in the inner regions of PPD result in important photophoretic forces, which dramatically affect the radial motion of solid bodies. Moreover, dust particles have different equilibrium locations in the disc depending on their size and their intrinsic density. The radial transport towards the outer parts of the disc is more efficient for silicates than for iron particles, which has important implications for meteoritic composition. Our results indicate that photophoresis must be taken into account in the inner regions of PPD to fully understand the dynamics and the evolution of the dust composition.

It is widely accepted that rocky planets form in the inner regions of protoplanetary disks (PPD) about 1 - 10 AU from the star. However, theoretical calculations show that when particles reach the size for which the radial migration is the fastest they tend to be accreted very efficiently by the star. This is known as the radial-drift barrier. We explore the photophoresis in the inner regions of PPD as a possible mechanism for preventing the accretion of solid bodies onto the star. Photophoresis is the thermal creep induced by the momentum exchange of an illuminated solid particle with the surrounding gas. Recent laboratory experiments predict that photophoresis would be able to stop the inward drift of macroscopic bodies (from 1 mm to 1 m in size). This extra force has been included in our two-fluid (gas+dust) SPH code in order to study its efficiency. We show that the conditions of pressure and temperature encountered in the inner regions of PPD result in strong dynamical effects on the dust particles due to photophoresis. Our simulations show that there is a radial and a vertical sorting of the dust grains according to their sizes and their intrinsic densities. Thus, our calculations support the fact that photophoresis is a mechanism which can have a strong effect on the morphology of the inner regions of PPD, ultimately affecting the fate of planetesimals.