Search

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 (, 17 pages, 14 figures

Many stars are in binaries or higher-order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the dynamics of gas and dust in discs around each component of a binary system under the hypothesis that the evolution of the circumprimary and the circumsecondary discs is independent. It is known from previous studies that the secular evolution of the gas in binary discs is hastened due to the tidal interactions with their hosting stars. Here we prove that binarity affects dust dynamics too, possibly in a more dramatic way than the gas. In particular, the presence of a stellar companion significantly reduces the amount of solids retained in binary discs because of a faster, more efficient radial drift, ultimately shortening their lifetime. We prove that how rapidly discs disperse depends both on the binary separation, with discs in wider binaries living longer, and on the disc viscosity. Although the less-viscous discs lose high amounts of solids in the earliest stages of their evolution, they are dissipated slowly, while those with higher viscosities show an opposite behaviour. The faster radial migration of dust in binary discs has a striking impact on planet formation, which seems to be inhibited in this hostile environment, unless other disc substructures halt radial drift further in. We conclude that if planetesimal formation were viable in binary discs, this process would take place on very short time scales., 18 pages, 11 figures; accepted for publication in MNRAS

The recently discovered exoplanets in binary or higher-order multiple stellar systems sparked a new interest in the study of proto-planetary discs in stellar aggregations. Here we focus on disc solids, as they make up the reservoir out of which exoplanets are assembled and dominate (sub-)millimetre disc observations. These observations suggest that discs in binary systems are fainter and smaller than in isolated systems. In addition, disc dust sizes are consistent with tidal truncation only if they orbit very eccentric binaries. In a previous study we showed that the presence of a stellar companion hastens the radial migration of solids, shortening disc lifetime and challenging planet formation. In this paper we confront our theoretical and numerical results with observations: disc dust fluxes and sizes from our models are computed at ALMA wavelengths and compared with Taurus and $\rho$ Ophiuchus data. A general agreement between theory and observations is found. In particular, we show that the dust disc sizes are generally smaller than the binary truncation radius due to the combined effect of grain growth and radial drift: therefore, small disc sizes do not require implausibly high eccentricities to be explained. Furthermore, the observed binary discs are compatible within $1\sigma$ with a quadratic flux-radius correlation similar to that found for single-star discs and show a close match with the models. However, the observational sample of resolved binary discs is still small and additional data are required to draw more robust conclusions on the flux-radius correlation and how it depends on the binary properties., Comment: 20 pages, 10 figures; accepted for publication in MNRAS

The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M, Comment: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

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)., 11 pages, 5 figures, accepted for publication in MNRAS

Tidal disruption events (TDEs) are among the brightest transients in the optical, ultraviolet, and X-ray sky. These flares are set into motion when a star is torn apart by the tidal field of a massive black hole, triggering a chain of events which is -- so far -- incompletely understood. However, the disruption process has been studied extensively for almost half a century, and unlike the later stages of a TDE, our understanding of the disruption itself is reasonably well converged. In this Chapter, we review both analytical and numerical models for stellar tidal disruption. Starting with relatively simple, order-of-magnitude physics, we review models of increasing sophistication, the semi-analytic ``affine formalism,'' hydrodynamic simulations of the disruption of polytropic stars, and the most recent hydrodynamic results concerning the disruption of realistic stellar models. Our review surveys the immediate aftermath of disruption in both typical and more unusual TDEs, exploring how the fate of the tidal debris changes if one considers non-main sequence stars, deeply penetrating tidal encounters, binary star systems, and sub-parabolic orbits. The stellar tidal disruption process provides the initial conditions needed to model the formation of accretion flows around quiescent massive black holes, and in some cases may also lead to directly observable emission, for example via shock breakout, gravitational waves or runaway nuclear fusion in deeply plunging TDEs., Comment: Review chapter in book: 'The Tidal Disruption of Stars by Massive Black Holes', Space Science Reviews, Springer. Comments welcome

Recent observations have shown that circumbinary discs can be misaligned with respect to the binary orbital plane.The lack of spherical symmetry, together with the non-planar geometry of these systems, causes differential precession which might induce the propagation of warps. While gas dynamics in such environments is well understood, little is known about dusty discs. In this work, we analytically study the problem of dust traps formation in misaligned circumbinary discs. We find that pile-ups may be induced not by pressure maxima, as the usual dust traps, but by a difference in precession rates between the gas and dust. Indeed, this difference makes the radial drift inefficient in two locations, leading to the formation of two dust rings whose position depends on the system parameters. This phenomenon is likely to occur to marginally coupled dust particles $(\text{St}\gtrsim1)$ as both the effect of gravitational and drag force are considerable. We then perform a suite of three-dimensional SPH numerical simulations to compare the results with our theoretical predictions. We explore the parameter space, varying stellar mass ratio, disc thickness, radial extension, and we find a general agreement with the analytical expectations. Such dust pile-up prevents radial drift, fosters dust growth and may thus promote the planet formation in circumbinary discs., Accepted for publication in MNRAS, 12 pages, 12 figures, 2 tables, 2 appendices

Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars., 57 pages, 17 figures. Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Corrected several typos from previous submission

Understanding how giant planets form requires observational input from directly imaged protoplanets. We used VLT/NACO and VLT/SPHERE to search for companions in the transition disc of 2MASS J19005804-3645048 (hereafter CrA-9), an accreting M0.75 dwarf with an estimated age of 1-2 Myr. We found a faint point source at $\sim$0.7'' separation from CrA-9 ($\sim$108 au projected separation). Our 3-epoch astrometry rejects a fixed background star with a $5\sigma$ significance. The near-IR absolute magnitudes of the object point towards a planetary-mass companion. However, our analysis of the 1.0-3.8$\mu$m spectrum extracted for the companion suggests it is a young M5.5 dwarf, based on both the 1.13-$\mu$m Na index and comparison with templates of the Montreal Spectral Library. The observed spectrum is best reproduced with high effective temperature ($3057^{+119}_{-36}$K) BT-DUSTY and BT-SETTL models, but the corresponding photometric radius required to match the measured flux is only $0.60^{+0.01}_{-0.04}$ Jovian radius. We discuss possible explanations to reconcile our measurements, including an M-dwarf companion obscured by an edge-on circum-secondary disc or the shock-heated part of the photosphere of an accreting protoplanet. Follow-up observations covering a larger wavelength range and/or at finer spectral resolution are required to discriminate these two scenarios., Comment: 24 pages, 14 figures, 4 tables, to be published in MNRAS

IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$\alpha$ imaging and Keck astrometric measures of the binary to develop a coherent dynamical model of this system. The system features an azimuthal asymmetry detected at the western edge of the cavity in Band~7 observations and a wiggling outflow. Dust emission in ALMA Band 4 from the proximity of the individual stars suggests the presence of marginally resolved circumstellar discs. We estimate the binary orbital parameters from the measured arc of the orbit from Keck and ALMA astrometry. We further constrain these estimates using considerations from binary-disc interaction theory. We finally perform three SPH gas + dust simulations based on the theoretical constraints; we post-process the hydrodynamic output using radiative transfer Monte Carlo methods and directly compare the models with observations. Our results suggest that a highly eccentric $e\sim 0.5\textrm{--}0.7$ equal mass binary, with a semi-major axis of $\sim 55$ au, and small/moderate orbital plane vs. circumbinary disc inclination $\theta\lesssim 30^\circ$ provides a good match with observations. A dust mass of $\sim 1.5\times 10^{-4} {\rm M_\odot}$ best reproduces the flux in Band 7 continuum observations. Synthetic CO line emission maps qualitatively capture both the emission from the central region and the non-Keplerian nature of the gas motion in the binary proximity., Comment: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is particularly interesting if we consider that gravitational instabilities could be a fundamental first step for the formation of planetesimals and planets. In this Letter, we model the rotation curve obtained by CO data of Elias 2-27 with a theoretical rotation curve including both the disk self-gravity and the star contribution to the gravitational potential. We compare this model with a purely Keplerian one and with a simple power-law function. We find that (especially for the $^{13}$CO isotopologue) the rotation curve is better described by considering not only the star, but also the disk self-gravity. We are thus able to obtain for the first time a dynamical estimate of the disk mass of $0.08\pm0.04\,M_{\odot}$ and the star mass of $0.46\pm0.03\,M_{\odot}$ (in the more general case), the latter being comparable with previous estimates. From these values, we derive that the disk is 17$\%$ of the star mass, meaning that it could be prone to gravitational instabilities. This result would strongly support the hypothesis that the two spiral arms are generated by gravitational instabilities., Comment: 14 pages, 9 figures, 2 tables, accepted for publication on ApJL

It is likely that young protostellar disks undergo a self-gravitating phase. Such systems are characterized by the presence of a spiral pattern that can be either in a quasi-steady state or in a nonlinear unstable condition. This spiral wave affects both the gas dynamics and kinematics, resulting in deviations from the Keplerian rotation. Recently, a lot of attention has been devoted to kinematic studies of planet-forming environments, and we are now able to measure even small perturbations of velocity field ($\lesssim 1 \%$ of the Keplerian speed) thanks to high spatial and spectral resolution observations of protostellar disks. In this work, we investigate the kinematic signatures of gravitational instability: we perform an analytical study of the linear response of a self-gravitating disk to a spiral-like perturbation, focusing our attention on the velocity field perturbations. We show that unstable disks have clear kinematic imprints into the gas component across the entire disk extent, due to the GI spiral wave perturbation, resulting in deviations from Keplerian rotation. The shape of these signatures depends on several parameters, but they are significantly affected by the cooling factor: by detecting these features, we can put constraints on protoplanetary disk cooling., Comment: Accepted for publication in APJL, 13 pages, 4 figures

Accretion onto black holes is an efficient mechanism in converting the gas mass-energy into energetic outputs as radiation, wind and jet. Tidal disruption events, in which stars are tidally torn apart and then accreted onto supermassive black holes, offer unique opportunities of studying the accretion physics as well as the wind and jet launching physics across different accretion regimes. In this review, we systematically describe and discuss the models that have been developed to study the accretion flows and jets in tidal disruption events. A good knowledge of these physics is not only needed for understanding the emissions of the observed events, but also crucial for probing the general relativistic space-time around black holes and the demographics of supermassive black holes via tidal disruption events., Comment: Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review, "The Tidal Disruption of Stars by Massive Black Holes," vol. 79

Observations with the Atacama Large Millimeter/Submillimeter array (ALMA) have dramatically improved our understanding of the site of exoplanet formation: protoplanetary discs. However, many basic properties of these discs are not well-understood. The most fundamental of these is the total disc mass, which sets the mass budget for planet formation. Discs with sufficiently high masses can excite gravitational instability and drive spiral arms that are detectable with ALMA . Although spirals have been detected in ALMA observations of the dust , their association with gravitational instability, and high disc masses, is far from clear. Here we report a prediction for kinematic evidence of gravitational instability. Using hydrodynamics simulations coupled with radiative transfer calculations, we show that a disc undergoing such instability has clear kinematic signatures in molecular line observations across the entire disc azimuth and radius which are independent of viewing angle. If these signatures are detected, it will provide the clearest evidence for the occurrence of gravitational instability in planet-forming discs, and provide a crucial way to measure disc masses., 10 pages, 7 figures, accepted to The Astrophysical Journal

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($\alpha_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($\tau\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation., Comment: accepted for publication in ApJ

Determining the mechanisms that drive the evolution of protoplanetary disks is a necessary step to understand how planets form. Here we measured the mass accretion rate for young stellar objects at age >5 Myr, a critical test for the current models of disk evolution. We present the analysis of the spectra of 36 targets in the ~5-10 Myr old Upper Scorpius region for which disk masses were measured with ALMA. We find that the mass accretion rates in this sample of old but still survived disks are similarly high as those of the younger (10^-9 Msun/yr while having low disk masses. Furthermore, the median values of the measured mass accretion rates in the disk mass ranges where our sample is complete at a level ~60-80% are compatible in these three regions. At the same time, the spread of mass accretion rates at any given disk mass is still >0.9 dex even at age>5 Myr. These results are in contrast with simple models of viscous evolution, which would predict that the values of the mass accretion rate diminish with time, and a tighter correlation with disk mass at age>5 Myr. Similarly, simple models of internal photoevaporation cannot reproduce the observed mass accretion rates, while external photoevaporation might explain the low disk masses and high accretion rates. A partial possible solution to the discrepancy with the viscous models is that the gas-to-dust ratio of the disks at >5 Myr is significantly different and higher than the canonical 100, as suggested by some dust and gas disk evolution models. The results shown here require the presence of several inter-playing processes, such as detailed dust evolution, external photoevaporation and possibly MHD winds, to explain the secular evolution of protoplanetary disks., Comment: Accepted on Astronomy & Astrophysics. Language edited version

Context. A large portion of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question of whether and how planets in binary systems form. Protoplanetary disks are the birthplaces of planets, and characterizing them is crucial in order to understand the planet formation process. Aims. Our goal is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks. We also aim to trace evidence for binary-disk interactions. Methods. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyzed the observed disk morphology and substructures. We ran 2D hydrodynamical models to simulate the evolution of the circumbinary ring over the lifetime of the disk. Results. The disk and also the cavity and the inner region are highly structured, with several shadowed regions, spiral structures, and streamer-like filaments. Some of these are detected here for the first time. The streamer-like filaments appear to connect the outer ring with the northern arc. Their azimuthal spacing suggests that they may be generated through periodic perturbations by the binary, which tear off material from the inner edge of the outer disk once during each orbit. By comparing observations to hydrodynamical simulations, we find that the main features, in particular, the gap size, but also the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with a semimajor axis of ~35 au on an orbit coplanar with the circumbinary ring., Astronomy & Astrophysics, 639, ISSN:0004-6361, ISSN:1432-0746

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability., Comment: Accepted for publication in the Astrophysical Journal. The resolution of several figures has been downgraded to comply with the size limit of arXiv submissions

Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events., A chapter in the ISSI review book "The Tidal Disruption of Stars by Massive Black Holes", to be published in Space Science Reviews

Numerical simulations have historically played a major role in understanding the hydrodynamics of the tidal disruption process. Given the complexity of the geometry of the system, the challenges posed by the problem have indeed stimulated much work on the numerical side. Smoothed Particles Hydrodynamics methods, for example, have seen their very first applications in the context of tidal disruption and still play a major role to this day. Likewise, initial attempts at simulating the evolution of the disrupted star with the so-called affine method have been historically very useful. In this Chapter, we provide an overview of the numerical techniques used in the field and of their limitations, and summarize the work that has been done to simulate numerically the tidal disruption process., 34 pages, Accepted chapter for Springer Space Science Reviews book on Tidal Disruption Events

Recent mm-wavelength surveys performed with the Atacama Large Millimeter Array (ALMA) have revealed protoplanetary discs characterized by rings and gaps. A possible explanation for the origin of such rings is the tidal interaction with an unseen planetary companion. The protoplanetary disc around DS Tau shows a wide gap in the ALMA observation at 1.3 mm. We construct a hydrodynamical model for the dust continuum observed by ALMA assuming the observed gap is carved by a planet between one and five Jupiter masses. We fit the shape of the radial intensity profile along the disc major axis varying the planet mass, the dust disc mass, and the evolution time of the system. The best fitting model is obtained for a planet with $M_{\rm p}=3.5\,M_{\rm Jup}$ and a disc with $M_{\rm dust}= 9.6\cdot10^{-5}\,M_{\odot}$. Starting from this result, we also compute the expected signature of the planet in the gas kinematics, as traced by CO emission. We find that such a signature (in the form of a "kink" in the channel maps) could be observed by ALMA with a velocity resolution between $0.2-0.5\,\rm{kms}^{-1}$ and a beam size between 30 and 50 mas., 16 pages, 15 figures, accepted for publication on MNRAS

Binary systems exert a gravitational torque on misaligned discs orbiting them, causing differential precession which may produce disc warping and tearing. While this is well understood for gas-only discs, misaligned cirumbinary discs of gas and dust have not been thoroughly investigated. We perform SPH simulations of misaligned gas and dust discs around binaries to investigate the different evolution of these two components. We choose two different disc aspect ratios: A thin case for which the gas disc always breaks, and a thick one where a smooth warp develops throughout the disc. For each case, we run simulations of five different dust species with different degrees of coupling with the gas component, varying in Stokes number from 0.002 (strongly coupled dust) to 1000 (effectively decoupled dust). We report two new phenomena: First, large dust grains in thick discs pile up at the warp location, forming narrow dust rings, due to a difference in precession between the gas and dust components. These pile ups do not form at gas pressure maxima, and hence are different from conventional dust traps. This effect is most evident for St ~ 10 - 100. Second, thin discs tear and break only in the gas, while dust particles with St > 10 form a dense dust trap due to the steep pressure gradient caused by the break in the gas. We find that dust with St < 0.02 closely follow the gas particles, for both thin and thick discs, with radial drift becoming noticeable only for the largest grains in this range., Accepted for publication in MNRAS, 11 pages, 9 figures

A tidal disruption event (TDE) occurs when a star passes too close to a supermassive black hole and gets torn apart by its gravitational tidal field. After the disruption, the stellar debris form an expanding gaseous stream. The morphology and evolution of this stream is particularly interesting as it ultimately determines the observational properties of the event itself. In this work we perform 3D hydrodynamical simulations of the TDE of a star modelled as a polytropic sphere of index {\gamma} = 5/3, and study the gravitational stability of the resulting gas stream. We provide an analytical solution for the evolution of the stream in the bound, unbound and marginally bound case, that allows us to describe the stream properties and analyse the time-scales of the physical processes involved, applying a formalism developed in star formation context. Our results are that, when fragmentation occurs, it is fueled by the failure of pressure in supporting the gas against its self-gravity. We also show that a stability criterion that includes also the stream gas pressure proves to be far more accurate than one that only considers the black hole tidal forces, giving analytical predictions of the time evolution of the various forces associated to the stream. Our results point out that fragmentation occurs on timescales longer compared with the observational windows of these events and is thus not expected to give rise to significant observational features., Comment: 12 pages,15 figures, accepted for publication in MNRAS

At the final stages of a supermassive black hole coalescence, the emission of gravitational waves will efficiently remove energy and angular momentum from the binary orbit, allowing the separation between the compact objects to shrink. In the scenario where a circumprimary disc is present, a squeezing phase will develop, in which the tidal interaction between the disc and the secondary black hole could push the gas inwards, enhancing the accretion rate on to the primary and producing what is known as an electromagnetic precursor. In this context, using 3D hydrodynamic simulations, we study how an adiabatic circumprimary accretion disc responds to the varying gravitational potential as the secondary falls onto the more massive object. We included a cooling prescription controlled by the parameter \beta = \Omega t_{cool}, which will determine how strong the final accretion rate is: a hotter disc is thicker, and the tidal interaction is suppressed for the gas outside the binary plane. Our main results are that for scenarios where the gas cannot cool fast enough (\beta>30) the disc becomes thick and renders the system invisible, while for \beta, Comment: to be published in MNRAS

Recent ALMA observations of the protoplanetary disc around HD~169142 reveal a peculiar structure made of concentric dusty rings: a main ring at $\sim$20 au, a triple system of rings at $\sim 55-75$ au in millimetric continuum emission and a perturbed gas surface density from the $^{12}$CO,$^{13}$CO and C$^{18}$O$(J=2-1)$ surface brightness profile. In this Letter, we perform three-dimensional numerical simulations and radiative transfer modeling exploring the possibility that two giant planets interacting with the disc and orbiting in resonant locking can be responsible for the origin of the observed dust inner rings structure. We find that in this configuration the dust structure is actually long lived while the gas mass of the disc is accreted onto the star and the giant planets, emptying the inner region. In addition, we also find that the innermost planet is located at the inner edge of the dust ring, and can accrete mass from the disc, generating a signature in the dust ring shape that can be observed in mm ALMA observations., 8 pages, 5 figures, accepted Apj Letters

We present 3 mm ALMA continuum and line observations at resolutions of 6.5 au and 13 au respectively, toward the Class 0 system IRAS 16293-2422 A. The continuum observations reveal two compact sources towards IRAS 16293-2422 A, coinciding with compact ionized gas emission previously observed at radio wavelengths (A1 and A2), confirming the long-known radio sources as protostellar. The emission towards A2 is resolved and traces a dust disk with a FWHM size of ~12 au, while the emission towards A1 sets a limit to the FWHM size of the dust disk of ~4 au. We also detect spatially resolved molecular kinematic tracers near the protostellar disks. Several lines of the J=5-4 rotational transition of HNCO, NH2CHO and t-HCOOH are detected, with which we derived individual line-of-sight velocities. Using these together with the CS (J=2-1), we fit Keplerian profiles towards the individual compact sources and derive masses of the central protostars. The kinematic analysis indicates that A1 and A2 are a bound binary system. Using this new context for the previous 30 years of VLA observations, we fit orbital parameters to the relative motion between A1 and A2 and find the combined protostellar mass derived from the orbit is consistent with the masses derived from the gas kinematics. Both estimations indicate masses consistently higher (0.5< M1, Comment: Accepted for publication in ApJ. 23 pages, 14 figures, 4 tables

We report the discovery of a low-mass stellar companion around the young Herbig Be star MWC 297. We performed multi-epoch high-contrast imaging in the near infrared (NIR) with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. The companion is found at projected separation of 244.7$\pm$13.2 au and a position angle of 176.4$\pm$0.1 deg. The large separation supports formation via gravitational instability. From the spectrum, we estimate a mass of 0.1-0.5 M$_{\odot}$, the range conveying uncertainties in the extinction of the companion and in evolutionary models at young ages. The orbit coincides with a gap in the dust disk inferred from the Spectral Energy Distribution (SED). The young age ($\lesssim$ 1 Myr) and mass ratio with the central star ($\sim 0.01$) makes the companion comparable to PDS 70~b, suggesting a relation between formation scenarios and disk dynamics., Comment: 4 figures

The disc surrounding PDS 70, with two directly imaged embedded giant planets, is an ideal laboratory to study planet-disc interaction. We present three-dimensional smoothed particle hydrodynamics simulations of the system. In our simulations, planets, which are free to migrate and accrete mass, end up in a locked resonant configuration that is dynamically stable. We show that features observed at infrared (scattered light) and millimetre (thermal continuum) wavelengths are naturally explained by the accretion stream onto the outer planet, without requiring a circumplanetary disc around planet c. We post-processed our near-infrared synthetic images in order to account for observational biases known to affect high-contrast images. Our successful reproduction of the observations indicates that planet-disc dynamical interactions alone are sufficient to explain the observations of PDS 70., 15 pages, 11 figures

In this paper we derive the gravitational wave stochastic background from tidal disruption events (TDEs). We focus on both the signal emitted by main sequence stars disrupted by super-massive black holes (SMBHs) in galaxy nuclei, and on that from disruptions of white dwarfs by intermediate mass black holes (IMBHs) located in globular clusters. We show that the characteristic strain $h_{\rm c}$'s dependence on frequency is shaped by the pericenter distribution of events within the tidal radius, and under standard assumptions $h_{\rm c} \propto f^{-1/2}$. This is because the TDE signal is a burst of gravitational waves at the orbital frequency of the closest approach. In addition, we compare the background characteristic strains with the sensitivity curves of the upcoming generation of space-based gravitational wave interferometers: the Laser Interferometer Space Antenna (LISA), TianQin, ALIA, the DECI-hertz inteferometer Gravitational wave Observatory (DECIGO) and the Big Bang Observer (BBO). We find that the background produced by main sequence stars might be just detected by BBO in its lowest frequency coverage, but it is too weak for all the other instruments. On the other hand, the background signal from TDEs with white dwarfs will be within reach of ALIA, and especially of DECIGO and BBO, while it is below the LISA and TianQin sensitive curves. This background signal detection will not only provide evidence for the existence of IMBHs up to redshift $z\sim 3$, but it will also inform us on the number of globular clusters per galaxy and on the occupation fraction of IMBHs in these environments., Comment: Accepted for Publications in MNRAS. 11 pages, 5 figures

Protoplanetary discs are the site of star and planet formation, and their evolution and consequent dispersal deeply affect the formation of planetary systems. In the standard scenario they evolve on timescales ~Myr due to the viscous transport of angular momentum. The analytical self-similar solution for their evolution predicts also specific disc isochrones in the accretion rate - disc mass plane. However, photoevaporation by radiation emitted by the central star is likely to dominate the gas disc dispersal of the innermost region, introducing another (shorter) timescale for this process. In this paper, we include the effect of internal (X and EUV) photoevaporation on the disc evolution, finding numerical solutions for a population of protoplanetary discs. Our models naturally reproduce the expected quick dispersal of the inner region of discs when their accretion rates match the rate of photoevaporative mass loss, in line with previous studies. We find that photoevaporation preferentially removes the lightest discs in the sample. The net result is that, counter-intuitively, photoevaporation increases the average disc mass in the sample, by dispersing the lightest discs. At the same time, photoevaporation also reduces the mass accretion rate by cutting the supply of material from the outer to the inner disc. In a purely viscous framework, this would be interpreted as the result of a longer viscous evolution, leading to an overestimate of the disc age. Our results thus show that photoevaporation is a necessary ingredient to include when interpreting observations of large disc samples with measured mass accretion rates and disc masses. Photoevaporation leaves a characteristic imprint on the shape of the isochrone. Accurate data in the accretion rate - disc mass plane in the low disc mass region therefore give clues on the typical photoevaporation rate., 8 pages, 7 figures. Accepted for publication on MNRAS

This paper reports on a new analysis of archival ALMA $870\,\mu$m dust continuum observations. Along with the previously observed bright inner ring ($r \sim 20-40\,$au), two addition substructures are evident in the new continuum image: a wide dust gap, $r \sim 40-150\,$au, and a faint outer ring ranging from $r \sim 150\,$au to $r \sim 250\,$au and whose presence was formerly postulated in low-angular-resolution ALMA cycle 0 observations but never before observed. Notably, the dust emission of the outer ring is not homogeneous, and it shows two prominent azimuthal asymmetries that resemble an eccentric ring with eccentricity $e = 0.07 $. The characteristic double-ring dust structure of HD 100546 is likely produced by the interaction of the disk with multiple giant protoplanets. This paper includes new smoothed-particle-hydrodynamic simulations with two giant protoplanets, one inside of the inner dust cavity and one in the dust gap. The simulations qualitatively reproduce the observations, and the final masses and orbital distances of the two planets in the simulations are 3.1 $M_{J}$ at 15 au and 8.5 $M_{J}$ at 110 au, respectively. The massive outer protoplanet substantially perturbs the disk surface density distribution and gas dynamics, producing multiple spiral arms both inward and outward of its orbit. This can explain the observed perturbed gas dynamics inward of 100 au as revealed by ALMA observations of CO. Finally, the reduced dust surface density in the $\sim 40-150\,$au dust gap can nicely clarify the origin of the previously detected H$_2$O gas and ice emission., Comment: Accepted for publication

Spatially resolved observations of protoplanetary discs are revealing that their inner regions can be warped or broken from the outer disc. A few mechanisms are known to lead to such 3D structures; among them, the interaction with a stellar companion. We perform a 3D SPH simulation of a circumbinary disc misaligned by $60^\circ$ with respect to the binary orbital plane. The inner disc breaks from the outer regions, precessing as a rigid body, and leading to a complex evolution. As the inner disc precesses, the misalignment angle between the inner and outer discs varies by more than $100^\circ$. Different snapshots of the evolution are post-processed with a radiative transfer code, in order to produce observational diagnostics of the process. Even though the simulation was produced for the specific case of a circumbinary disc, most of the observational predictions hold for any disc hosting a precessing inner rim. Synthetic scattered light observations show strong azimuthal asymmetries, where the pattern depends strongly on the misalignment angle between inner and outer disc. The asymmetric illumination of the outer disc leads to azimuthal variations of the temperature structure, in particular in the upper layers, where the cooling time is short. These variations are reflected in asymmetric surface brightness maps of optically thick lines, as CO $J$=3-2. The kinematical information obtained from the gas lines is unique in determining the disc structure. The combination of scattered light images and (sub-)mm lines can distinguish between radial inflow and misaligned inner disc scenarios., 17 pages, 17 figures. Accepted for publication in MNRAS

Mergers are fundamental to the standard paradigm of galaxy evolution, and provide a natural formation mechanism for supermassive black hole binaries. The formation process of such a binary can have a direct impact on the rate at which stars are tidally disrupted by one or the other black hole, and the luminous signature of the tidal disruption itself can have distinct imprints of a binary companion. In this chapter we review our current understanding of the influence of black hole binarity on the properties of tidal disruption events. We discuss the rates of tidal disruption by supermassive black hole binaries, the impact of a second black hole on the fallback of debris and the formation of an accretion flow, and the prospects for detection of tidal disruption events by supermassive black hole binaries., 22 pages, 5 figures. Accepted chapter for Springer Space Science Reviews book on tidal disruption events

While planets are commonly discovered around main-sequence stars, the processes leading to their formation are still far from being understood. Current planet population synthesis models, which aim to describe the planet formation process from the protoplanetary disk phase to the time exoplanets are observed, rely on prescriptions for the underlying properties of protoplanetary disks where planets form and evolve. The recent development in measuring disk masses and disk-star interaction properties, i.e., mass accretion rates, in large samples of young stellar objects demand a more careful comparison between the models and the data. We performed an initial critical assessment of the assumptions made by planet synthesis population models by looking at the relation between mass accretion rates and disk masses in the models and in the currently available data. We find that the currently used disk models predict mass accretion rate in line with what is measured, but with a much lower spread of values than observed. This difference is mainly because the models have a smaller spread of viscous timescales than what is needed to reproduce the observations. We also find an overabundance of weakly accreting disks in the models where giant planets have formed with respect to observations of typical disks. We suggest that either fewer giant planets have formed in reality or that the prescription for planet accretion predicts accretion on the planets that is too high. Finally, the comparison of the properties of transition disks with large cavities confirms that in many of these objects the observed accretion rates are higher than those predicted by the models. On the other hand, PDS70, a transition disk with two detected giant planets in the cavity, shows mass accretion rates well in line with model predictions., Accepted for publication on A&A Letter, 5 pages + 4 of appendix, 3 figures - revised version after language editing

We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~\mu$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths., Comment: ApJ accepted, 25 pages, 11 figures

Recent observations of protoplanetary discs reveal disc substructures potentially caused by embedded planets. We investigate how the gas surface density in discs changes the observed morphology in scattered light and dust continuum emission. Assuming that disc substructures are due to embedded protoplanets, we combine hydrodynamical modelling with radiative transfer simulations of dusty protoplanetary discs hosting planets. The response of different dust species to the gravitational perturbation induced by a planet depends on the drag stopping time - a function of the generally unknown local gas density. Small dust grains, being stuck to the gas, show spirals. Larger grains decouple, showing progressively more axisymmetric (ring-like) substructure as decoupling increases with grain size or with the inverse of the gas disc mass. We show that simultaneous modelling of scattered light and dust continuum emission is able to constrain the Stokes number, ${\rm St}$. Hence, if the dust properties are known, this constrains the local gas surface density, $\Sigma_{\rm gas}$, at the location of the structure, and hence the total gas mass. In particular, we found that observing ring-like structures in mm-emitting grains requires ${\rm St} \gtrsim 0.4$ and therefore $\Sigma_{\rm gas} \lesssim 0.4\,\textrm{g/cm}^{2}$. We apply this idea to observed protoplanetary discs showing substructures both in scattered light and in the dust continuum., Comment: Accepted for publication on MNRAS, 12 pages, 6 figures

Gravitational waves can be emitted by accretion discs if they undergo instabilities that generate a time varying mass quadrupole. In this work we investigate the gravitational signal generated by a thick accretion disc of $1 M_{\odot}$ around a static super-massive black hole of $10^{6}M_{\odot}$, assumed to be formed after the tidal disruption of a solar type star. This torus has been shown to be unstable to a global non-axisymmetric hydrodynamic instability, the Papaloizou-Pringle instability, in the case where it is not already accreting and has a weak magnetic field. We start by deriving analytical estimates of the maximum amplitude of the gravitational wave signal, with the aim to establish its detectability by the Laser Interferometer Space Antenna (LISA). Then, we compare these estimates with those obtained through a numerical simulation of the torus, made with a 3D smoothed particle hydrodynamics code. Our numerical analysis shows that the measured strain is two orders of magnitude lower than the maximum value obtained analytically. However, accretion discs affected by the Papaloizou-Pringle instability may still be interesting sources for LISA, if we consider discs generated after deeply penetrating tidal disruptions of main sequence stars of higher mass., Accepted for Publications in MNRAS. 9 pages, 4 figures

The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and spatially resolve the disks around all primary stars, and those around eight secondary and one tertiary star. The dust radii of disks around multiple stellar systems are smaller than those around single stars in the same stellar mass range and in the same region. The disks in multiple stellar systems also show a steeper decay of the millimeter continuum emission at the outer radius than disks around single stars, suggestive of the impact of tidal truncation on the shape of the disks in multiple systems. However, the observed ratio between the dust disk radii and the observed separation of the stars in the multiple systems is consistent with analytic predictions of the effect of tidal truncation only if the eccentricities of the binaries are rather high (typically >0.5), or if the observed dust radii are a factor of two smaller than the gas radii, as is typical for isolated systems. Similar high-resolution studies targeting the gaseous emission from disks in multiple stellar systems are required to resolve this question., Accepted on Astronomy & Astrophysics; 15 pages plus appendix, revised version after language editing

A correlation between proto-planetary disc radii and sub-mm fluxes has been recently reported. In this Letter we show that the correlation is a sensitive probe of grain growth processes. Using models of grain growth and drift, we have shown in a companion paper that the observed disc radii trace where the dust grains are large enough to have a significant sub-mm opacity. We show that the observed correlation emerges naturally if the maximum grain size is set by radial drift, implying relatively low values of the viscous $\alpha$ parameter $ \lesssim 0.001$. In this case the relation has an almost universal normalisation, while if the grain size is set by fragmentation the flux at a given radius depends on the dust-to-gas ratio. We highlight two observational consequences of the fact that radial drift limits the grain size. The first is that the dust masses measured from the sub-mm could be overestimated by a factor of a few. The second is that the correlation should be present also at longer wavelengths (e.g. 3mm), with a normalisation factor that scales as the square of the observing frequency as in the optically thick case., Comment: 6 pages, 3 figures; accepted on MNRAS letters

Aims. R Coronae Australis (R CrA) is the brightest star of the Coronet nebula of the Corona Australis (CrA) star forming region. It has very red colors, probably due to dust absorption and it is strongly variable. High contrast instruments allow for an unprecedented direct exploration of the immediate circumstellar environment of this star. Methods. We observed R CrA with the near-IR channels (IFS and IRDIS) of SPHERE at VLT. In this paper, we used four different epochs, three of them from open time observations while one is from the SPHERE guaranteed time. The data were reduced using the DRH pipeline and the SPHERE Data Center. On the reduced data we implemented custom IDL routines with the aim to subtract the speckle halo.We have also obtained pupil-tracking H-band (1.45-1.85 micron) observations with the VLT/SINFONI near-infrared medium-resolution (R~3000) spectrograph. Results. A companion was found at a separation of 0.156" from the star in the first epoch and increasing to 0.18400 in the final one. Furthermore, several extended structures were found around the star, the most noteworthy of which is a very bright jet-like structure North-East from the star. The astrometric measurements of the companion in the four epochs confirm that it is gravitationally bound to the star. The SPHERE photometry and the SINFONI spectrum, once corrected for extinction, point toward an early M spectral type object with a mass between 0.3 and 0.55 M?. The astrometric analyis provides constraints on the orbit paramenters: e~0.4, semi-major axis at 27-28 au, inclination of ~ 70{\deg} and a period larger than 30 years. We were also able to put constraints of few MJup on the mass of possible other companions down to separations of few tens of au., Comment: 11 pages, 13 figures, accepted by A&A

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet., Accepted for publication in A&A, 11 pages, 8 figures

ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters., Comment: 9 pages, 5 figures, re-submitted to MNRAS after referee's comments

The combination of high-resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission (λ = 1.3 mm) and CO isotopologues (12CO, 13CO, C18O, J = 2 - 1) integrated intensity from the disc around the nearby (d = 162 pc), intermediate-mass (M_{\star }=1.67 M_{\odot }) pre-main-sequence star CQ Tau. The data show an inner depression in continuum and in both 13CO and C18O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc spectral energy distribution. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between ~10-1 and ~10-3. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centred at R_dust=53 au and with a width of σ =13 au. Three-dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of {~ }20 au and with a mass of {~ } 6 - 9 M_Jup reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one-planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.

Proto-planetary disc surveys conducted with ALMA are measuring disc radii in multiple star forming regions. The disc radius is a fundamental quantity to diagnose whether discs undergo viscous spreading, discriminating between viscosity or angular momentum removal by winds as drivers of disc evolution. Observationally, however, the sub-mm continuum emission is dominated by the dust, which also drifts inwards, complicating the picture. In this paper we investigate, using theoretical models of dust grain growth and radial drift, how the radii of dusty viscous proto-planetary discs evolve with time. Despite the existence of a sharp outer edge in the dust distribution, we find that the radius enclosing most of the dust $\textit{mass}$ increases with time, closely following the evolution of the gas radius. This behaviour arises because, although dust initially grows and drifts rapidly onto the star, the residual dust retained on Myr timescales is relatively well coupled to the gas. Observing the expansion of the dust disc requires using definitions based on high fractions of the disc $\textit{flux}$ (e.g. 95 per cent) and very long integrations with ALMA, because the dust grains in the outer part of the disc are small and have a low sub-mm opacity. We show that existing surveys lack the sensitivity to detect viscous spreading. The disc radii they measure do not trace the mass radius or the sharp outer edge in the dust distribution, but the outer limit of where the grains have significant sub-mm opacity. We predict that these observed radii should shrink with time., 18 pages, 14 figures; accepted on MNRAS

The process of tidal disruption of stars by a supermassive black hole (SMBH) provides luminous UV and soft X-ray flares with peak luminosities of $\approx 10^{46}$ ergs/sec and duration of a few months. As part of a wider exploration of the effects of stellar rotation on the outcome of a TDE, we have performed hydrodynamical simulations of the disruption of a rotating star whose spin axis is opposite to the orbital axis. Such a retrograde rotation makes the star more resilient to tidal disruption, so that, even if its orbit reaches the formal tidal radius, it actually stays intact after the tidal encounter. However, the outer layers of the star are initially stripped away from the core, but then fall back onto the star itself, producing a newly formed accretion disc around the star. We estimate that the accretion rate onto the star would be strongly super-Eddington (for the star) and would result in an X-ray flare with luminosity of the order of $\approx 10^{40}$ ergs/sec and duration of a few months. We speculate that such events might be responsible for the known X-ray flares from Sgr A* in the recent past., Comment: 8 pages, 10 figures, accepted for publication in MNRAS

Stars are spheres of gas held together by self-gravity. When flying by a black hole, however, the star self-binding force can be overwhelmed by the black hole tides and the star can be torn apart. This is a physically rich and fascinating event which will be described by first introducing the concept of black hole from a mathematical point of view. We will then dive into the physics of the tidal disruption and proceed describing the accompanying electromagnetic flare and gravitational wave burst in the frequency range of the Laser Interferometer Space Antenna. This empowers such events to discover the elusive black holes with mass intermediate between the solar and the million/billion solar masses.

We estimate the accretion rates produced when a circumprimary gas disc is pushed into the primary supermassive black hole (SMBH) by the tidal force of the decaying secondary during a SMBH merger. Using the 3D Smoothed Particle Hydrodynamics (SPH) code {\sc phantom}, we extend previous investigations of co-planar discs to the case where the disc and binary orbital planes are misaligned. We consider a geometrically thin disc with inclination angles varying from 1 to 180 degrees and a binary with mass ratio $q$=$10^{-3}$. We find that discs with small inclination angles (< 10 degrees) produce an increase in luminosity exceeding the Eddington rate. By contrast, discs with inclinations between 20 and 30 degrees show a less pronounced rise in the accretion rate, whilst discs inclined by 180 degrees show no peak in the mass accretion rate. While previous analytic work predicted that the effective tidal torque drops with increasing inclination angle, we show that the misaligned snowplough effect remains important even for angles larger than the disc aspect ratio. The rise in the accretion rate produced by discs inclined at small angles to the binary orbit can produce an electromagnetic counterpart to the gravitational wave signal emitted from final stages of the binary orbital decay., Accepted 2018 December 11. Received 2018 December 4; in original form 2018 September 18

In this paper we analyse giant gap-opening planet migration in protoplanetary discs, focusing on the type II migration regime. According to standard type II theory, planets migrate at the same rate as the gas in the disc, as they are coupled to the disc viscous evolution; however, recent studies questioned this paradigm, suggesting that planets migrate faster than the disc material. We study the problem through 2D long-time simulations of systems consistent with type II regime, using the hydrodynamical grid code FARGO3D. Even though our simulations confirm the presence of an initial phase characterised by fast migration, they also reveal that the migration velocity slows down and eventually reaches the theoretical prediction if we allow the system to evolve for enough time. We find the same tendency to evolve towards the theoretical predictions at later times when we analyse the mass flow through the gap and the torques acting on the planet. This transient is related to the initial conditions of our (and previous) simulations, and is due to the fact that the shape of the gap has to adjust to a new profile, once the planet is set into motion. Secondly, we test whether the type II theory expectation that giant planet migration is driven by viscosity is consistent with our simulation by comparing simulations with the same viscosity and different disc mass (or viceversa). We find a good agreement with the theory, since when the discs are characterised by the same viscosity, the migration properties are the same., Comment: 12 pages, 8 figures, Accepted for publication in MNRAS

In the past years, high angular resolution observations have revealed that circumstellar discs appear in a variety of shapes with diverse substructures being ubiquitous. This has given rise to the question of whether these substructures are triggered by planet-disc interactions. Besides direct imaging, one of the most promising methods to distinguish between different disc shaping mechanisms is to study the kinematics of the gas disc. In particular, the deviations of the rotation profile from Keplerian velocity can be used to probe perturbations in the gas pressure profile that may be caused by embedded planets. In this paper we aim to analyze the gas brightness temperature and kinematics of the transitional disc around the CQ Tau star in order to resolve and characterize substructure in the gas, caused by possible perturbers. For our analysis we use spatially resolved ALMA observations of 12CO, 13CO and C18O (J=2-1). We further extract robust line centroids for each channel map and fit a number of Keplerian disc models to the velocity field. The gas kinematics of the CQ Tau disc present non-Keplerian features, showing bent and twisted iso-velocity curves in 12CO and 13CO. Significant spiral structures are detected between 10-180 au in both the brightness temperature and the rotation velocity of 12CO after subtraction of an azimuthally symmetric model, which may be tracing planet-disc interactions with an embedded planet or low-mass companion. We identify three spirals, two in the brightness temperature and one in the velocity residuals, spanning a large azimuth and radial extent. The brightness temperature spirals are morphologically connected to spirals observed in NIR scattered light in the same disc, indicating a common origin. Together with the observed large dust and gas cavity, the spirals support the hypothesis of a massive embedded companion in the CQ Tau disc., Comment: 16 pages, 15 figures; accepted for publication in A&A

A large number of tidal disruption event (TDE) candidates have been observed recently, often differing in their observational features. Two classes appear to stand out: X-ray and optical TDEs, the latter featuring lower effective temperatures and luminosities. These differences can be explained if the radiation detected from the two categories of events originates from different locations. In practice, this location is set by the evolution of the debris stream around the black hole and by the energy dissipation associated with it. In this paper, we build an analytical model for the stream evolution, whose dynamics is determined by both magnetic stresses and shocks. Without magnetic stresses, the stream always circularizes. The ratio of the circularization timescale to the initial stream period is $t_{\rm ev}/t_{\rm min} = 8.3 (M_{\rm h}/10^6 M_{\odot})^{-5/3} \beta^{-3}$, where $M_{\rm h}$ is the black hole mass and $\beta$ is the penetration factor. If magnetic stresses are strong, they can lead to the stream ballistic accretion. The boundary between circularization and ballistic accretion corresponds to a critical magnetic stresses efficiency $v_{\rm A}/v_{\rm c} \approx 10^{-1}$, largely independent of $M_{\rm h}$ and $\beta$. However, the main effect of magnetic stresses is to accelerate the stream evolution by strengthening self-crossing shocks. Ballistic accretion therefore necessarily occurs on the stream dynamical timescale. The shock luminosity associated to energy dissipation is sub-Eddington but decays as $t^{-5/3}$ only for a slow stream evolution. Finally, we find that the stream thickness rapidly increases if the stream is unable to cool completely efficiently. A likely outcome is its fast evolution into a thick torus, or even an envelope completely surrounding the black hole., Comment: 17 pages, 11 figures, accepted for publication in MNRAS. 3D visualizations of the model can be found at http://home.strw.leidenuniv.nl/~bonnerot/research.html