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<M2<2 Msun) than previous estimations using lower resolution observations of the gas kinematics. The ALMA high-resolution data provides a unique insight into the gas kinematics and masses of a young deeply embedded bound binary system.

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.

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.

A large fraction 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 if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks, and trace evidence for binary-disk interactions. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyze the observed disk morphology and substructures. We run 2D hydrodynamical models simulating the evolution of the circumbinary ring over the lifetime of the disk. The disk, as well as the cavity and the inner region are highly structured with several shadowed regions, spiral structures, and streamer-like filaments, some of them detected 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 by periodic perturbations by the binary, tearing 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, as well as the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with semi-major axis of $\sim$35 au on an orbit coplanar with the circumbinary ring.

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.

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 (<3 Myr old) star-forming regions of Lupus and Cha I, when considering the dependence on stellar and disk mass. In particular, several disks show high mass accretion rates >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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

GrailQuest (Gamma-ray Astronomy International Laboratory for Quantum Exploration of Space-Time) is an ambitious astrophysical mission concept that uses a fleet of small satellites, whose scientific objectives are discussed below. Within Quantum Gravity theories, different models for space-time quantisation predict an energy dependent speed for photons. Although the predicted discrepancies are minuscule, Gamma-Ray Bursts, occurring at cosmological distances, could be used to detect this signature of space-time granularity with a new concept of modular observatory of huge overall collecting area consisting in a fleet of small satellites in low orbits, with sub-microsecond time resolution and wide energy band (keV-MeV). The enormous number of collected photons will allow to effectively search these energy dependent delays. Moreover, GrailQuest will allow to perform temporal triangulation of high signal-to-noise impulsive events with arc-second positional accuracies: an extraordinary sensitive X-ray/Gamma all-sky monitor crucial for hunting the elusive electromagnetic counterparts of Gravitational Waves. A pathfinder of GrailQuest is already under development through the HERMES (High Energy Rapid Modular Ensemble of Satellites) project: a fleet of six 3U cube-sats to be launched by the end of 2021.

We perform 3D general relativistic smoothed particle hydrodynamics (GRSPH) simulations of tidal disruption events involving 1 $M_\odot$ stars and $10^6 M_\odot$ rotating supermassive black holes. We consider stars on initially elliptical orbits both in, and inclined to, the black hole equatorial plane. We confirm that stream-stream collisions caused by relativistic apsidal precession rapidly circularise the disrupted material into a disc. For inclined trajectories we find that nodal precession induced by the black hole spin (i.e. Lense-Thirring precession) inhibits stream-stream collisions only in the first orbit, merely causing a short delay in forming a disc, which is inclined to the black hole equatorial plane. We also investigate the effect of radiative cooling on the remnant disc structure. We find that with no cooling a thick, extended, slowly precessing torus is formed, with a radial extent of 5 au (for orbits with a high penetration factor). Radiatively efficient cooling produces a narrow, rapidly precessing ring close to pericentre. We plot the energy dissipation rate, which tracks the pancake shock, stream-stream collisions and viscosity. We compare this to the effective luminosity due to accretion onto the black hole. We find energy dissipation rates of $\sim10^{45}$ erg s$^{-1}$ for stars disrupted at the tidal radius, and up to $\sim10^{47}$ erg s$^{-1}$ for deep encounters.

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).

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.

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.

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.

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.

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.

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.

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 ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $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 $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. 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 $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm 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.

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.

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.

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.

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.

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.

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.

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.

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.

The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes.

Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($\alpha=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases.

The binary system GG Tau A is observed to have a circumbinary disc with a dust ring located further out than expected, assuming a co-planar disc and a corresponding semi-major axis of 34 AU. Given the binary separation, this large cavity can be explained by relaxing the assumption of a co-planar disc and instead fit the observations with a mis-aligned circumbinary disc around an eccentric binary with a wider semi-major axis of 60 AU, consistent with fitting the proper motion data for the system. We run SPH simulations to check this possibility and indeed we find that a misalignment angle of 30 degrees and a binary eccentricity of 0.45 fit both the astrometric data and the disc cavity. However, such configuration could in principle be unstable to polar, rather that planar alignment. We investigate the secular evolution of this configuration and show that it is indeed stable throughout the disc lifetime.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

The disruption of a star by a supermassive black hole generates a sudden bright flare. Previous studies have focused on the disruption by single black holes, for which the fallback rate decays as~$\propto t^{-5/3}$. In this paper, we generalise the study to the case of a supermassive black hole binary (SMBHB), using both analytical estimates and hydrodynamical simulations, looking for specific observable signatures. The range of binary separation for which it is possible to distinguish between the disruption created by a single or a binary black hole concerns typically separations of order a few milliparsecs for a primary of mass $\sim 10^6M_{\odot}$. When the fallback rate is affected by the secondary, it undergoes two types interruptions, depending on the initial inclination $\theta$ of the orbit of the star relative to the plane of the SMBHB. For $\theta \lesssim 70^\circ$, periodic sharp interruptions occur and the time of first interruption depends on the distance of the secondary black hole with the debris. If $\theta \gtrsim 70^\circ$, a first smooth interruption occurs, but not always followed by a further recovery of the fallback rate. This implies that most of the TDEs around a SMBHB will undergo periodic sharp interruptions of their lightcurve.

We describe a new implementation of the one-fluid method in the SPH code Phantom to simulate the dynamics of dust grains in gas protoplanetary discs. We revise and extend previously developed algorithms by computing the evolution of a new fluid quantity that produces a more accurate and numerically controlled evolution of the dust dynamics. Moreover, by limiting the stopping time of uncoupled grains that violate the assumptions of the terminal velocity approximation, we avoid fatal numerical errors in mass conservation. We test and validate our new algorithm by running 3D SPH simulations of a large range of disc models with tightly- and marginally-coupled grains.

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.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.

During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.

We demonstrate that the compact, thick disc formed in a tidal disruption event may be unstable to non-axisymmetric perturbations in the form of the Papaloizou-Pringle instability. We show this can lead to rapid redistribution of angular momentum that can be parameterised in terms of an effective Shakura-Sunyaev $\alpha$ parameter. For remnants that have initially weak magnetic fields, this may be responsible for driving mass accretion prior to the onset of the magneto-rotational instability. For tidal disruptions around a $10^6$ M$_{\odot}$ black hole, the measured accretion rate is super-Eddington but is not sustainable over many orbits. We thus identify a method by which the torus formed in tidal disruption event may be significantly accreted before the magneto-rotational instability is established.

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.

Type-C QPOs in X-ray binaries have been often interpreted as a consequence of relativistic Lense-Thirring precession around a spinning black hole and they potentially offer a way to measure black hole spins and masses. The connection between relativistic precession and the resulting QPO has been made either in terms of a simplified model involving a single test particle producing the QPO, or in terms of a global model where a geometrically thick accretion flow precesses coherently as a rigid body. In this paper, we analyse similarities and differences between these two models, sometimes considered as in opposition to each other. We demonstrate that the former is the limiting case of the latter when the radial extent of the precessing flow is very small, and that solid lower limits to the black hole spin can be obtained by considering the test particle model alone. We also show that the global precession model naturally accounts for the range of frequencies observed for Type-C QPOs without the need to invoke a truncation of the inner accretion flow before it reaches the innermost stable circular orbit. Finally, we show that, in order to maintain rigid precession, the thick accretion flow should be radially narrow, and that if it extends beyond $10-10^2$ gravitational radii it aligns with the black hole spin too fast to produce a coherent QPO.

In this paper, we compare simple viscous diffusion models for the disc evolution with the results of recent surveys of the properties of young protoplanetary discs. We introduce the useful concept of `disc isochrones' in the accretion rate - disc mass plane and explore a set of Montecarlo realization of disc initial conditions. We find that such simple viscous models can provide a remarkable agreement with the available data in the Lupus star forming region, with the key requirement that the average viscous evolutionary timescale of the discs is comparable to the cluster age. Our models produce naturally a correlation between mass accretion rate and disc mass that is shallower than linear, contrary to previous results and in agreement with observations. We also predict that a linear correlation, with a tighter scatter, should be found for more evolved disc populations. Finally, we find that such viscous models can reproduce the observations in the Lupus region only in the assumption that the efficiency of angular momentum transport is a growing function of radius, thus putting interesting constraints on the nature of the microscopic processes that lead to disc accretion.

We investigate the relation between disk mass and mass accretion rate to constrain the mechanism of angular momentum transport in protoplanetary disks. Dust mass and mass accretion rate in Chamaeleon I are correlated with a slope close to linear, similar to the one recently identified in Lupus. We investigate the effect of stellar mass and find that the intrinsic scatter around the best-fit Mdust-Mstar and Macc-Mstar relations is uncorrelated. Disks with a constant alpha viscosity can fit the observed relations between dust mass, mass accretion rate, and stellar mass, but over-predict the strength of the correlation between disk mass and mass accretion rate when using standard initial conditions. We find two possible solutions. 1) The observed scatter in Mdust and Macc is not primoridal, but arises from additional physical processes or uncertainties in estimating the disk gas mass. Most likely grain growth and radial drift affect the observable dust mass, while variability on large time scales affects the mass accretion rates. 2) The observed scatter is primordial, but disks have not evolved substantially at the age of Lupus and Chamaeleon I due to a low viscosity or a large initial disk radius. More accurate estimates of the disk mass and gas disk sizes in a large sample of protoplanetary disks, either through direct observations of the gas or spatially resolved multi-wavelength observations of the dust with ALMA, are needed to discriminate between both scenarios or to constrain alternative angular momentum transport mechanisms such as MHD disk winds.

We carried out a Bayesian homogeneous determination of the orbital parameters of 231 transiting giant planets (TGPs) that are alone or have distant companions; we employed DE-MCMC methods to analyse radial-velocity (RV) data from the literature and 782 new high-accuracy RVs obtained with the HARPS-N spectrograph for 45 systems over 3 years. Our work yields the largest sample of systems with a transiting giant exoplanet and coherently determined orbital, planetary, and stellar parameters. We found that the orbital parameters of TGPs in non-compact planetary systems are clearly shaped by tides raised by their host stars. Indeed, the most eccentric planets have relatively large orbital separations and/or high mass ratios, as expected from the equilibrium tide theory. This feature would be the outcome of high-eccentricity migration (HEM). The distribution of $\alpha=a/a_R$, where $a$ and $a_R$ are the semi-major axis and the Roche limit, for well-determined circular orbits peaks at 2.5; this also agrees with expectations from the HEM. The few planets of our sample with circular orbits and $\alpha >5$ values may have migrated through disc-planet interactions instead of HEM. By comparing circularisation times with stellar ages, we found that hot Jupiters with $a < 0.05$ au have modified tidal quality factors $10^{5} < Q'_p < 10^{9}$, and that stellar $Q'_s > 10^{6}-10^{7}$ are required to explain the presence of eccentric planets at the same orbital distance. As a by-product of our analysis, we detected a non-zero eccentricity for HAT-P-29; we determined that five planets that were previously regarded to have hints of non-zero eccentricity have circular orbits or undetermined eccentricities; we unveiled curvatures caused by distant companions in the RV time series of HAT-P-2, HAT-P-22, and HAT-P-29; and we revised the planetary parameters of CoRoT-1b.

We present ALMA observations of the 2M1207 system, a young binary made of a brown dwarf with a planetary-mass companion at a projected separation of about 40 au. We detect emission from dust continuum at 0.89 mm and from the $J = 3 - 2$ rotational transition of CO from a very compact disk around the young brown dwarf. The small radius found for this brown dwarf disk may be due to truncation from the tidal interaction with the planetary-mass companion. Under the assumption of optically thin dust emission, we estimated a dust mass of 0.1 $M_{\oplus}$ for the 2M1207A disk, and a 3$\sigma$ upper limit of $\sim 1~M_{\rm{Moon}}$ for dust surrounding 2M1207b, which is the tightest upper limit obtained so far for the mass of dust particles surrounding a young planetary-mass companion. We discuss the impact of this and other non-detections of young planetary-mass companions for models of planet formation, which predict the presence of circum-planetary material surrounding these objects.

When a star gets tidally disrupted by a supermassive black hole, its magnetic field is expected to pervade its debris. In this paper, we study this process via smoothed particle magnetohydrodynamical simulations of the disruption and early debris evolution including the stellar magnetic field. As the gas stretches into a stream, we show that the magnetic field evolution is strongly dependent on its orientation with respect to the stretching direction. In particular, an alignment of the field lines with the direction of stretching induces an increase of the magnetic energy. For disruptions happening well within the tidal radius, the star compression causes the magnetic field strength to sharply increase by an order of magnitude at the time of pericentre passage. If the disruption is partial, we find evidence for a dynamo process occurring inside the surviving core due to the formation of vortices. This causes an amplification of the magnetic field strength by a factor of $\sim 10$. However, this value represents a lower limit since it increases with numerical resolution. For an initial field strength of 1 G, the magnetic field never becomes dynamically important. Instead, the disruption of a star with a strong 1 MG magnetic field produces a debris stream within which magnetic pressure becomes similar to gas pressure a few tens of hours after disruption. If the remnant of one or multiple partial disruptions is eventually fully disrupted, its magnetic field could be large enough to magnetically power the relativistic jet detected from Swift J1644+57. Magnetized streams could also be significantly thickened by magnetic pressure when it overcomes the confining effect of self-gravity.

Black hole X-ray transients show a variety of state transitions during their outburst phases, characterized by changes in their spectral and timing properties. In particular, power density spectra (PDS) show quasi periodic oscillations (QPOs) that can be related to the accretion regime of the source. We looked for type-C QPOs in the disc-dominated state (i.e. the high soft state) and in the ultra-luminous state in the RXTE archival data of 12 transient black hole X-ray binaries known to show QPOs during their outbursts. We detected 6 significant QPOs in the soft state that can be classified as type-C QPOs. Under the assumption that the accretion disc in disc-dominated states extends down or close to the innermost stable circular orbit (ISCO) and that type-C QPOs would arise at the inner edge of the accretion flow, we use the relativistic precession model (RPM) to place constraints on the black hole spin. We were able to place lower limits on the spin value for all the 12 sources of our sample while we could place also an upper limit on the spin for 5 sources.

The protoplanetary disk around the GGTau A binary system is so far one of the most studied young circumbinary disk. Observations of the dust continuum emission at sub-mm/mm wavelengths detected a dust ring located between 200AU and 300AU from the center of mass of the system. If the disk and the binary orbit are coplanar, given the observed projected separation of the stars in the binary system, the classical theory of tidal interaction between a binary system and its circumbinary disk predicts a truncation radius which is significantly smaller than the measured inner radius of the mm-sized dust ring. We investigate the origin of this dust ring structure in the GG Tau A disk, test whether the interaction between the binary and the disk can produce a dust trap at the location of the observed ring, and discuss an alternative scenario which invokes a misalignment between the disk and the stellar orbital planes. We run a set of hydrodynamical simulations for an orbit consistent with the astrometric solutions for the GG Tau A stellar proper motions, and for different disk temperature and viscosities. We then apply a dust evolution model in post-processing in order to to retrieve the distribution of mm-sized grains. Comparing models and observations, we show that, if the binary orbit and the disk were coplanar,the tidal truncation of the circumbinary disk would occur at a radius that is too small compared to the inner edge inferred by the dust observations, and that the pressure bump and the dust ring in the models would be located at <150AU from the center of mass of the stellar system. This shows that the GG Tau A disk cannot be coplanar with the binary orbital plane. We also discuss the viability of the misaligned disk scenario, suggesting that in order for dust trapping to occur at the observed radius, the disk and orbital plane must be misaligned by an angle of about 25 degrees.

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.

Most massive galaxies are thought to contain a supermassive black hole in their centre surrounded by a tenuous gas environment, leading to no significant emission. In these quiescent galaxies, tidal disruption events represent a powerful detection method for the central black hole. Following the disruption, the stellar debris evolves into an elongated gas stream, which partly falls back towards the disruption site and accretes onto the black hole producing a luminous flare. Using an analytical treatment, we investigate the interaction between the debris stream and the gas environment of quiescent galaxies. Although we find dynamical effects to be negligible, we demonstrate that Kelvin-Helmholtz instability can lead to the dissolution of the stream into the ambient medium before it reaches the black hole, likely dimming the associated flare. This result is robust against the presence of a typical stellar magnetic field and fast cooling within the stream. Furthermore, we find this effect to be enhanced for disruptions involving more massive black holes and/or giant stars. Consequently, although disruptions of evolved stars have been proposed as a useful probe of black holes with masses $\gtrsim 10^8 \, M_{\odot}$, we argue that the associated flares are likely less luminous than expected.

We investigate whether the rings, lopsided features and horseshoes observed at millimetre wavelengths in transitional discs can be explained by the dynamics of gas and dust at the edge of the cavity in circumbinary discs. We use 3D dusty smoothed particle hydrodynamics calculations to show that binaries with mass ratio $q \gtrsim 0.04$ drive eccentricity in the central cavity, naturally leading to a crescent-like feature in the gas density, which is accentuated in the mm dust grain population with intensity contrasts in mm-continuum emission of 10 or higher. We perform mock observations to demonstrate that these features closely match those observed by ALMA, suggesting that the origin of rings, dust horseshoes and other non-axisymmetric structures in transition discs can be explained by the presence of massive companions.

We present three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations investigating the dependence of the accretion rate on the disc thickness around an equal-mass, circular black hole binary system. We find that for thick/hot discs, with $H/R\gtrsim 0.1$, the binary torque does not prevent the gas from penetrating the cavity formed in the disc by the binary (in line with previous investigations). The situation drastically changes for thinner discs, in this case the mass accretion rate is suppressed, such that only a fraction (linearly dependent on $H/R$) of the available gas is able to flow within the cavity and accrete on to the binary. Extrapolating this result to the cold and thin accretion discs expected around supermassive black hole binary systems implies that this kind of systems accretes less material than predicted so far, with consequences not only for the electromagnetic and gravitational waves emissions during the late inspiral phase but also for the recoil speed of the black hole formed after binary coalescence, thus influencing also the evolutionary path both of the binary and of the host galaxy. Our results, being scale-free, are also applicable to equal mass, circular binaries of stellar mass black holes, such as the progenitor of the recently discovered gravitational wave source GW150914.

[Abridged] Star and planet formation are the complex outcomes of gravitational collapse and angular momentum transport mediated by protostellar and protoplanetary disks. In this review we focus on the role of gravitational instability in this process. We begin with a brief overview of the observational evidence for massive disks that might be subject to gravitational instability, and then highlight the diverse ways in which the instability manifests itself in protostellar and protoplanetary disks: the generation of spiral arms, small scale turbulence-like density fluctuations, and fragmentation of the disk itself. We present the analytic theory that describes the linear growth phase of the instability, supplemented with a survey of numerical simulations that aim to capture the non-linear evolution. We emphasize the role of thermodynamics and large scale infall in controlling the outcome of the instability. Despite apparent controversies in the literature, we show a remarkable level of agreement between analytic predictions and numerical results. We highlight open questions related to (1) the development of a turbulent cascade in thin disks, and (2) the role of mode-mode coupling in setting the maximum angular momentum transport rate in thick disks.

We identify two distinct physical mechanisms for dust gap opening by embedded planets in protoplanetary discs based on the symmetry of the drag-induced motion around the planet: I) A mechanism where low mass planets, that do not disturb the gas, open gaps in dust by tidal torques assisted by drag in the inner disc, but resisted by drag in the outer disc; and II) The usual, drag assisted, mechanism where higher mass planets create pressure maxima in the gas disc which the drag torque then acts to evacuate further in the dust. The first mechanism produces gaps in dust but not gas, while the second produces partial or total gas gaps which are deeper in the dust phase. Dust gaps do not necessarily indicate gas gaps.

We study accretion rates during the gravitational wave-driven merger of a binary supermassive black hole embedded in an accretion disc, formed by gas driven to the centre of the galaxy. We use 3D simulations performed with PHANTOM, a Smoothed Particle Hydrodynamics code. Contrary to previous investigations, we show that there is evidence of a "squeezing phenomenon", caused by the compression of the inner disc gas when the secondary black hole spirals towards the primary. This causes an increase in the accretion rates that always exceed the Eddington rate. We have studied the main features of the phenomenon for a mass ratio $q = 10^{-3}$ between the black holes, including the effects of numerical resolution, the secondary accretion radius and the disc thickness. With our disc model with a low aspect ratio, we show that the mass expelled from the orbit of the secondary is negligible ($< 5\%$ of the initial disc mass), different to the findings of previous 2D simulations with thicker discs. The increase in the accretion rates in the last stages of the merger leads to an increase in luminosity, making it possible to detect an electromagnetic precursor of the gravitational wave signal emitted by the coalescence.

A tidal disruption event (TDE) occurs when a star wanders close enough to a black hole to be disrupted by its tidal force. The debris of a tidally disrupted star are expected to form an accretion disc around the supermassive black hole. The light curves of these events sometimes show a quasi-periodic modulation of the flux that can be associated with the precession of the accretion disc due to the Lense-Thirring ("frame-dragging") effect. Since the initial star orbit is in general inclined with respect to the black hole spin, this misalignment combined with the Lense-Thirring effect leads to a warp in the disc. In this paper we provide a simple model of the system composed by a thick and narrow accretion disc surrounding a spinning supermassive black hole, with the aim to: (a) compute the expected precession period as a function of the system parameters, (b) discuss the conditions that have to be satisfied in order to have rigid precession, (c) investigate the alignment process, highlighting how different mechanisms play a role leading the disc and the black hole angular momenta into alignment.

The potential of tidal disruption of stars to probe otherwise quiescent supermassive black holes cannot be exploited, if their dynamics is not fully understood. So far, the observational appearance of these events has been derived from analytical extrapolations of the debris dynamical properties just after disruption. By means of hydrodynamical simulations, we investigate the subsequent fallback of the stream of debris towards the black hole for stars already bound to the black hole on eccentric orbits. We demonstrate that the debris circularize due to relativistic apsidal precession which causes the stream to self-cross. The circularization timescale varies between 1 and 10 times the period of the star, being shorter for more eccentric and/or deeper encounters. This self-crossing leads to the formation of shocks that increase the thermal energy of the debris. If this thermal energy is efficiently radiated away, the debris settle in a narrow ring at the circularization radius with shock-induced luminosities of $\sim 10-10^3 \, L_{\rm Edd}$. If instead cooling is impeded, the debris form an extended torus located between the circularization radius and the semi-major axis of the star with heating rates $\sim 1-10^2 \, L_{\rm Edd}$. Extrapolating our results to parabolic orbits, we infer that circularization would occur via the same mechanism in $\sim 1$ period of the most bound debris for deeply penetrating encounters to $\sim 10$ for grazing ones. We also anticipate the same effect of the cooling efficiency on the structure of the disc with associated luminosities of $\sim 1-10 \, L_{\rm Edd}$ and heating rates of $\sim 0.1-1 \, L_{\rm Edd}$. In the latter case of inefficient cooling, we deduce a viscous timescale generally shorter than the circularization timescale. This suggests an accretion rate through the disc tracing the fallback rate, if viscosity starts acting promptly.

We explain the axisymmetric gaps seen in recent long-baseline observations of the HL Tau protoplanetary disc with the Atacama Large Millimetre/Submillimetre Array (ALMA) as being due to the different response of gas and dust to embedded planets in protoplanetary discs. We perform global, three dimensional dusty smoothed particle hydrodynamics calculations of multiple planets embedded in dust/gas discs which successfully reproduce most of the structures seen in the ALMA image. We find a best match to the observations using three embedded planets with masses of 0.2, 0.27 and 0.55 $M_{\rm J}$ in the three main gaps observed by ALMA, though there remain uncertainties in the exact planet masses from the disc model.

We determine the true and the projected obliquity of HAT-P-36 and WASP-11/HAT-P-10 systems, respectively, which are both composed of a relatively cool star and a hot-Jupiter planet. Thanks to the high-resolution spectrograph HARPS-N, we observed the Rossiter-McLaughlin effect for both the systems by acquiring precise radial-velocity measurements during planetary transit events. We also present photometric observations comprising six light curves covering five transit events, obtained using three medium-class telescopes and the telescope-defocussing technique. One transit of WASP-11/HAT-P-10 was followed simultaneously from two observatories. The three transit light curves of HAT-P-36b show anomalies that are attributable to starspot complexes on the surface of the parent star, in agreement with the analysis of its spectra that indicate a moderate activity. By analysing the complete HATNet data set of HAT-P-36, we estimated the stellar rotation period by detecting a periodic photometric modulation in the light curve caused by star spots, obtaining Prot=15.3 days, which implies that the inclination of the stellar rotational axis with respect to the line of sight is 65 degree. We used the new spectroscopic and photometric data to revise the main physical parameters and measure the sky-projected misalignment angle of the two systems. We found \lambda=-14 degree for HAT-P-36 and \lambda=7 degree for WASP-11/HAT-P-10, indicating in both cases a good spin-orbit alignment. In the case of HAT-P-36, we also measured its real obliquity, which turned out to be 25 degrees.

Interactions between a supermassive black hole binary and the surrounding accretion disc can both assist the binary inspiral and align the black hole spins to the disc angular momentum. While binary migration is due to angular-momentum transfer within the circumbinary disc, the spin-alignment process is driven by the mass accreting on to each black hole. Mass transfer between different disc components thus couples the inspiral and the alignment process together. Mass is expected to leak through the cavity cleared by the binary, and preferentially accretes on to the lighter (secondary) black hole which orbits closer to the disc edge. Low accretion rate on to the heavier (primary) black hole slows the alignment process down. We revisit the problem and develop a semi-analytical model to describe the coupling between gas-driven inspiral and spin alignment, finding that binaries with mass ratio q<~0.2 approach the gravitational-wave driven inspiral in differential misalignment: light secondaries prevent primaries from aligning. Binary black holes with misaligned primaries are ideal candidates for precession effects in the strong-gravity regime and may suffer from moderately large (~1500 km/s) recoil velocities.

Tidal disruption events occur when a star passes too close to a massive black hole and it is totally ripped apart by tidal forces. It may also happen that the star is not close enough to the black hole to be totally disrupted and a less dramatic event might happen. If the stellar orbit is bound and highly eccentric, just like some stars in the centre of our own Galaxy, repeated flares should occur. When the star approaches the black hole tidal radius at periastron, matter might be stripped resulting in lower intensity outbursts recurring once every orbital period. We report on Swift observations of a recent bright flare from the galaxy IC 3599 hosting a middle-weight black hole, where a possible tidal disruption event was observed in the early 1990s. By light curve modelling and spectral fitting we can consistently account for the events as the non-disruptive tidal stripping of a star into a highly eccentric orbit. The recurrence time is 9.5 yr. IC 3599 is also known to host a low-luminosity active galactic nucleus. Tidal stripping from this star over several orbital passages might be able to spoon-feed also this activity.

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.

Tidal disruption events (TDEs) occur when a star, passing too close to a massive black hole, is ripped apart by tidal forces. A less dramatic event occurs if the star orbits just outside the tidal radius, resulting in a mild stripping of mass. Thus, if a star orbits a central black hole on one of these bound eccentric orbits, weaker outbursts will occur recurring every orbital period. Thanks to five Swift observations, we observed a recent flare from the close by (92 Mpc) galaxy IC 3599, where a possible TDE was already observed in December 1990 during the Rosat All-Sky Survey. By light curve modeling and spectral fitting, we account for all these events as the non-disruptive tidal stripping of a single star into a 9.5 yr highly eccentric bound orbit. This is the first example of periodic partial tidal disruptions, possibly spoon-feeding the central black hole.

In this paper, we revisit the issue of estimating the "fossil" disc mass in the circumprimary disc, during the merger of a supermassive black hole binary. As the binary orbital decay speeds up due to the emission of gravitational waves, the gas in the circumprimary disc might be forced to accrete rapidly and could in principle provide a significant electromagnetic counterpart to the gravitational wave emission. Since the luminosity of such flare is proportional to the gaseous mass in the circumprimary disc, estimating such mass accurately is important. Previous investigations of this issue have produced contradictory results, with some authors estimating super-Eddington flares and large disc mass, while others suggesting that the "fossil" disc mass is very low, even less than a Jupiter mass. Here, we perform simple 1D calculations to show that such very low estimates of the disc mass are an artifact of the specific implementation of the tidal torque in 1D models. In particular, for moderate mass ratios of the binary, the usual formula for the torque used in 1D models significantly overestimates the width of the gap induced by the secondary and this artificially leads to a very small leftover circumprimary disc. Using a modified torque, calibrated to reproduce the correct gap width as estimated by 3D models, leads to fossil disc masses of the order of one solar mass. The rapid accretion of the whole circumprimary disc would produce peak luminosities of the order of 1-20 times the Eddington luminosity. Even if a significant fraction of the gas escapes accretion by flowing out the secondary orbit during the merger (an effect not included in our calculations), we would still predict close to Eddington luminosities that might be easily detected.

We revisit the TrES-4 system parameters based on high-precision HARPS-N radial-velocity measurements and new photometric light curves. A combined spectroscopic and photometric analysis allows us to determine a spectroscopic orbit with an amplitude $K=51\pm3$ m s$^{-1}$. The derived mass of TrES-4b is found to be $M_{\rm p} = 0.49\pm0.04 \rm M_{Jup}$, significantly lower than previously reported. Combined with the large radius ($R_{\rm p} = 1.84_{-0.09}^{+0.08} \rm R_{Jup}$) inferred from our analysis, TrES-4b becomes the second-lowest density transiting hot Jupiter known. We discuss several scenarios to explain the puzzling discrepancy in the mass of TrES-4b in the context of the exotic class of highly inflated transiting giant planets.

XO-2 is the first confirmed wide stellar binary system where the almost twin components XO-2N and XO-2S have planets. This stimulated a detailed characterization study of the stellar and planetary components based on new observations. We collected high-resolution spectra with the HARPS-N spectrograph and multi-band light curves. Spectral analysis led to an accurate determination of the stellar atmospheric parameters and characterization of the stellar activity. We collected 14 transit light curves of XO-2Nb used to improve the transit parameters. Photometry provided accurate magnitude differences between the stars and a measure of their rotation periods. The iron abundance of XO-2N was found to be +0.054 dex greater, within more than 3-sigma, than that of XO-2S. We confirm a long-term variation in the radial velocities of XO-2N, and we detected a turn-over with respect to previous measurements. We suggest the presence of a second massive companion in an outer orbit or the stellar activity cycle as possible causes of the observed acceleration. The latter explanation seems more plausible with the present dataset. We obtained an accurate value of the projected spin-orbit angle for the XO-2N system (lambda=7+/-11 degrees), and estimated the real 3-D spin-orbit angle (psi=27 +12/-27 degrees). We measured the XO-2 rotation periods, and found a value of P=41.6 days in the case of XO-2N, in excellent agreement with the predictions. The period of XO-2S appears shorter, with an ambiguity between 26 and 34.5 days that we cannot solve with the present dataset alone. XO-2N appears to be more active than the companion, and this could be due to the fact that we sampled different phases of their activity cycle, or to an interaction between XO-2N and its hot Jupiter that we could not confirm.

In this paper we present simulated Atacama Large Millimetre/sub-millimetre Array (ALMA) observations of self-gravitating circumstellar discs with different properties in size, mass and inclination, located in four of the most extensively studied and surveyed star-forming regions. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation and representative dust opacities, we have initially constructed maps of the expected emission at sub-mm wavelengths of a large sample of discs with different properties. We have then simulated realistic observations of discs as they may appear with ALMA using the Common Astronomy Software Application ALMA simulator. We find that, with a proper combination of antenna configuration and integration time, the spiral structure characteristic of self-gravitating discs is readily detectable by ALMA over a wide range of wavelengths at distances comparable to TW Hydrae ($\sim 50 \,$pc), Taurus - Auriga and Ophiucus ($\sim 140 \,$pc) star-forming regions. However, for discs located in Orion complex ($\sim 400 \,$pc) only the largest discs in our sample (outer radius of 100 au) show a spatially resolved structure while the smaller ones (outer radius of 25 au) are characterized by a spiral structure that is not conclusively detectable with ALMA.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

On 2011 August 11, INTEGRAL discovered the hard X-ray source IGR J17361-4441 near the centre of the globular cluster NGC 6388. Follow up observations with Chandra showed the position of the transient was inconsistent with the cluster dynamical centre, and thus not related to its possible intermediate mass black hole. The source showed a peculiar hard spectrum (Gamma \approx 0.8) and no evidence of QPOs, pulsations, type-I bursts, or radio emission. Based on its peak luminosity, IGR J17361-4441 was classified as a very faint X-ray transient, and most likely a low-mass X-ray binary. We re-analysed 200 days of Swift/XRT observations, covering the whole outburst of IGR J17361-4441 and find a t^{-5/3} trend evident in the light curve, and a thermal emission component that does not evolve significantly with time. We investigate whether this source could be a tidal disruption event, and for certain assumptions find an accretion efficiency epsilon \approx 3.5E-04 (M_{Ch}/M) consistent with a massive white dwarf, and a disrupted minor body mass M_{mb}=1.9E+27(M/M_{Ch}) g in the terrestrial-icy planet regime. These numbers yield an inner disc temperature of the order kT_{in} \approx 0.04 keV, consistent with the blackbody temperature of kT_{in} \approx 0.08 keV estimated by spectral fitting. Although the density of white dwarfs and the number of free-floating planets are uncertain, we estimate the rate of planetary tidal disruptions in NGC 6388 to be in the range 3E-06 to 3E-04 yr^{-1}. Averaged over the Milky Way globular clusters, the upper limit value corresponds to 0.05 yr^{-1}, consistent with the observation of a single event by INTEGRAL and Swift.

We performed an intensive radial velocity monitoring of XO-2S, the wide companion of the transiting planet-host XO-2N, using HARPS-N at TNG in the framework of the GAPS programme. The radial velocity measurements indicate the presence of a new planetary system formed by a planet that is slightly more massive than Jupiter at 0.48 au and a Saturn-mass planet at 0.13 au. Both planetary orbits are moderately eccentric and were found to be dynamically stable. There are also indications of a long-term trend in the radial velocities. This is the first confirmed case of a wide binary whose components both host planets, one of which is transiting, which makes the XO-2 system a unique laboratory for understanding the diversity of planetary systems.

We are entering the era in which observations of protoplanetary discs properties can indirectly probe the presence of massive planets or low mass stellar companions interacting with the disc. In particular, the detection of warped discs can provide important clues to the properties of the star-disc system. In this paper we show how observations of warped discs can be used to infer the dynamical properties of the systems. We concentrate on circumbinary discs, where the mass of the secondary can be planetary. First, we provide some simple relations that link the amplitude of the warp in the linear regime to the parameters of the system. Secondly, we apply our method to the case of TW Hya, a transition disc for which a warp has been proposed based on spectroscopic observations. Assuming values for the disc and stellar parameters from observations, we conclude that, in order for a warp induced by a planetary companion to be detectable, the planet mass should be large ($M_{\rm p} \approx 10 - 14M_{\rm J}$) and the disc should be viscous ($\alpha \approx 0.15 - 0.25$). We also apply our model to LkCa 15 and T Cha, where a substellar companion has been detected within the central cavity of the transition discs.

We present new 3 mm ATCA data of two Class I Young Stellar Objects in the Ophiucus star forming region: Elias29 and WL12. For our analysis we compare them with archival 1.1 mm SMA data. In the (u,v) plane the two sources present a similar behavior: a nearly constant non-zero emission at long baselines, which suggests the presence of an unresolved component and an increase of the fluxes at short baselines, related to the presence of an extended envelope. Our data analysis leads to unusually low values of the spectral index $\alpha_{\rm 1.1-3mm}$, which may indicate that mm-sized dust grains have already formed both in the envelopes and in the disk-like structures at such early stages. To explore the possible scenarios for the interpretation of the sources we perform a radiative transfer modeling using a Monte Carlo code, in order to take into account possible deviations from the Rayleigh-Jeans and optically thin regimes. Comparison between the model outputs and the observations indicates that dust grains may form aggregates up to millimeter size already in the inner regions of the envelopes of Class I YSOs. Moreover, we conclude that the embedded disk-like structures in our two Class Is are probably very compact, in particular in the case of WL12, with outer radii down to tens of AU.

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

KH 15D is a protostellar binary system that shows a peculiar light curve. In order to model it, a narrow circumbinary precessing disc has been invoked, but a proper dynamical model has never been developed. In this paper, we analytically address the issue of whether such a disc can rigidly precess around KH 15D, and we relate the precessional period to the main parameters of the system. Then, we simulate the disc's dynamics by using a 1D model developed in a companion paper, such that the warp propagates into the disc as a bending wave, which is expected to be the case for protostellar discs. The validity of such an approach has been confirmed by comparing its results with full 3D SPH simulations on extended discs. In the present case, we use this 1D code to model the propagation of the warp in a narrow disc. If the inner truncation radius of the disc is set by the binary tidal torques at {\sim} 1 AU, we find that the disc should extend out to 6-10 AU (depending on the models), and is therefore wider than previously suggested. Our simulations show that such a disc does reach an almost steady state, and then precesses as a rigid body. The disc displays a very small warp, with a tilt inclination that increases with radius in order to keep the disc in equilibrium against the binary torque. However, for such wider discs, the presence of viscosity leads to a secular decay of the tilt on a timescale of {\approx} 3000 ({\alpha}/0.05)^(-1) years, where {\alpha} is the disc viscosity parameter. The presence of a third body (such as a planet), orbiting at roughly 10 AU might simultaneously explain the outer truncation of the disc and the maintenance of the tilt for a prolonged time.

In this paper we analyse the propagation of warps in protostellar circumbinary discs. We use these systems as a test environment in which to study warp propagation in the bending-wave regime, with the addition of an external torque due to the binary gravitational potential. In particular, we want to test the linear regime, for which an analytic theory has been developed. In order to do so, we first compute analytically the steady state shape of an inviscid disc subject to the binary torques. The steady state tilt is a monotonically increasing function of radius. In the absence of viscosity, the disc does not present any twist. Then, we compare the time-dependent evolution of the warped disc calculated via the known linearised equations both with the analytic solutions and with full 3D numerical simulations, which have been performed with the PHANTOM SPH code using 2 million particles. We find a good agreement both in the tilt and in the phase evolution for small inclinations, even at very low viscosities. Moreover, we have verified that the linearised equations are able to reproduce the diffusive behaviour when {\alpha} > H/R, where {\alpha} is the disc viscosity parameter. Finally, we have used the 3D simulations to explore the non-linear regime. We observe a strongly non-linear behaviour, which leads to the breaking of the disc. Then, the inner disc starts precessing with its own precessional frequency. This behaviour has already been observed with numerical simulations in accretion discs around spinning black holes. The evolution of circumstellar accretion discs strongly depends on the warp evolution. Therefore the issue explored in this paper could be of fundamental importance in order to understand the evolution of accretion discs in crowded environments, when the gravitational interaction with other stars is highly likely, and in multiple systems.

We discuss the importance and potential contribution of Athena+ to the science questions related to gamma-ray bursts, tidal disruption events and supernova shock break-out. Athena+ will allow breakthrough observations involving high-resolution X-ray spectroscopic observations of high-z gamma-ray bursts, observations of tidal disruption events tailored to determine the mass and potentially the spin of the black hole responsible for the tidal disruption and observations of X-rays from the supernova shock breakout providing a measure of the radius of the exploding star or of the companion in the case of type Ia supernovae. We briefly discuss survey facilities that will provide triggers to these events envisaged to be operational around 2028.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We present SCUBA-2 850um observations for 7 very low mass stars (VLMS) and brown dwarfs (BDs): 3 in Taurus, 4 in the TWA, and all classical T Tauri (cTT) analogs. We detect 2 of the 3 Taurus disks, but none of the TWA ones. Our 3sigma limits correspond to a dust mass of 1.2 MEarth in Taurus and a mere 0.2 MEarth in the TWA (3--10x deeper than previous work). We combine our data with other sub-mm/mm surveys of Taurus, rho Oph and the TWA to investigate trends in disk mass and grain growth during the cTT phase. We find : (1) The minimum disk outer radius required to explain the upper envelope of sub-mm/mm fluxes is 100 AU for intermediate-mass stars, solar-types and VLMS, and 20 AU for BDs. (2) While the upper envelope of disk masses increases with Mstar from BDs to VLMS to solar-types, no increase is seen from solar-type to intermediate-mass stars. We propose this is due to enhanced photoevaporation around intermediate masses. (3) Many disks around Taurus and rho Oph intermediate-mass and solar-type stars evince an opacity index beta of 0--1, indicating large grains. Of the only four VLMS/BDs in these regions with multi-wavelength data, three are consistent with large grains, though optically thick disks are not ruled out. (4) For the TWA VLMS (TWA 30A, B), combining our fluxes with Mdot and ages suggests substantial grain growth by 10 Myr. The degree of grain growth in the TWA BDs (2M1207A, SSPM1102) remains largely unknown. (5) A Bayesian analysis shows that mean(log[Mdisk/Mstar]) = -2.4, roughly constant all the way from intermediate-mass stars to VLMS/BDs, and (6) the disk mass in close solar-type Taurus binaries is significantly lower than in singles (by a factor of 10), while that in wide solar-type Taurus binaries is closer to that in singles (lower by a factor of 3). (7) We discuss the implications for planet formation, and for the dependence of Mdot on Mstar.

A long-term multi-purpose observational programme has started with HARPS-N@TNG aimed to characterise the global architectural properties of exoplanetary systems. In this first paper we fully characterise the transiting system Qatar-1. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter-McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric transit light-curves are analysed and a spectroscopic characterisation of the host star atmospheric parameters is performed based on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment of the system; this information, combined with the spectroscopic determination of the host star properties, allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of lambda = -8.4+-7.1 deg. The planet, with a mass of 1.33+-0.05 M_J, is found to be significantly more massive than previously reported. The host star is confirmed to be metal-rich ([Fe/H]= 0.20+-0.10) and slowly rotating (vsinI = 1.7+-0.3 km/s), though moderately active, as indicated by strong chromospheric emission in the Ca II H&K line cores (logR'_HK about -4.60). The system is well aligned and fits well within the general lambda vs Teff trend. We definitely rule out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions.

In this contribution, I review the recent developments on the modelling of the lightcurve of tidal disruption events. Our understanding has evolved significantly from the earlier seminal results that imply a simple power-law decay of the bolometric light curve as $t^{-5/3}$. We now know that the details of the rise to the peak of the lightcurve is determined mainly by the internal structure of the disrupted star. We also have improved models for the disc thermal emission, showing that in this case the decline of the luminosity with time should be much flatter than the standard $t^{-5/3}$ law, especially in optical and UV wavelengths, while the X-ray lightcurve is generally best suited to track the bolometric one. Finally, we are just starting to explore the interesting general relativistic effects that might arise for such events, for which the tidal radius lies very close to the black hole event horizon. In particular, I describe here some possible evidences for relativistic Lense-Thirring precession from the light curve of the event Swift J1644.

In this Letter we revisit arguments suggesting that the Bardeen-Petterson effect can coalign the spins of a central supermassive black hole binary accreting from a circumbinary (or circumnuclear) gas disc. We improve on previous estimates by adding the dependence on system parameters, and noting that the nonlinear nature of warp propagation in a thin viscous disc affects alignment. This reduces the disc's ability to communicate the warp, and can severely reduce the effectiveness of disc-assisted spin alignment. We test our predictions with a Monte Carlo realization of random misalignments and accretion rates and we find that the outcome depends strongly on the spin magnitude. We estimate a generous upper limit to the probability of alignment by making assumptions which favour it throughout. Even with these assumptions, about 40% of black holes with $a \gtrsim 0.5$ do not have time to align with the disc. If the residual misalignment is not small and it is maintained down to the final coalescence phase this can give a powerful recoil velocity to the merged hole. Highly spinning black holes are thus more likely of being subject to strong recoils, the occurrence of which is currently debated.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

The present observational understanding of the evolution of the mass accretion rates (Macc) in pre-main sequence stars is limited by the lack of accurate measurements of Macc over homogeneous and large statistical samples of young stars. Such observational effort is needed to properly constrain the theory of star formation and disk evolution. Based on HST/WFPC2 observations, we present a study of Macc for a sample of \sim 700 sources in the Orion Nebula Cluster, ranging from the Hydrogen-burning limit to M\ast \sim 2M\odot. We derive Macc from both the U-band excess and the H{\alpha} luminosity (LH{\alpha}), after determining empirically both the shape of the typical accretion spectrum across the Balmer jump and the relation between the accretion luminosity (Lacc) and LH{\alpha}, that is Lacc/L\odot = (1.31\pm0.03)\cdotLH{\alpha}/L\odot + (2.63\pm 0.13). Given our large statistical sample, we are able to accurately investigate relations between Macc and the parameters of the central star such as mass and age. We clearly find Macc to increase with stellar mass, and decrease over evolutionary time, but we also find strong evidence that the decay of Macc with stellar age occurs over longer timescales for more massive PMS stars. Our best fit relation between these parameters is given by: log(Macc/M\odot\cdotyr)=(-5.12 \pm 0.86) -(0.46 \pm 0.13) \cdot log(t/yr) -(5.75 \pm 1.47)\cdot log(M\ast/M\odot) + (1.17 \pm 0.23)\cdot log(t/yr) \cdot log(M\ast/M\odot). These results also suggest that the similarity solution model could be revised for sources with M\ast > 0.5M\odot. Finally, we do not find a clear trend indicating environmental effects on the accretion properties of the sources.

We investigate the evolution of the surface density of a circumbinary accretion disc after the mass loss induced by the merger of two supermassive black holes. We first introduce an analytical model, under the assumption of a disc composed of test particles, to derive the surface density evolution of the disc following the mass loss. The model predicts the formation of sharp density peaks in the disc; the model also allows us to compute the typical timescale for the formation of these peaks. To test and validate the model, we run numerical simulations of the process using the Smoothed Particle Hydrodynamics (SPH) code PHANTOM, taking fluid effects into account. We find good agreement in the shape and position of the peaks between the model and the simulations. In a fluid disc, however, the epicyclic oscillations induced by the mass loss can dissipate, and only some of the predicted peaks form in the simulation. To quantify how fast this dissipation proceeds, we introduce an appropriate parameter, and we show that it is effective in explaining the differences between the analytical, collisionless model and a real fluid disc.

We have revised our analysis of the SPIRE observations of 2MASSW J1207334-393254 (2M1207). Recent PACS observations show a bright source located ~25" east of 2M1207. There are issues in terms of the detection/non-detection of the bright source when comparing the Spitzer, WISE, and PACS observations. It is apparently inconsistent, perhaps due to variability or low signal-to-noise of the data. We have looked into the possible misidentification of the target, and have revised the measured SPIRE fluxes and the disc parameters for 2M1207. We have also reviewed which among the various formation mechanisms of this system would still be valid.

The tidal disruption of a solar mass star around a supermassive black hole has been extensively studied analytically and numerically. In these events the star develops into an elongated banana-shaped structure. After completing an eccentric orbit, the bound debris fall onto the black hole, forming an accretion disk and emitting radiation. The same process may occur on planetary scales, if a minor body orbits too close to its star. In the Solar System, comets fall directly onto our Sun or onto planets. If the star is a compact object, the minor body can become tidally disrupted. Indeed, one of the first mechanisms invoked to produce strong gamma-ray emission involved accretion of comets onto neutron stars in our Galaxy. Here we report that the peculiarities of the `Christmas' gamma-ray burst (GRB 101225A) can be explained by a tidal disruption event of a minor body around a Galactic isolated neutron star. This result would indicate either that minor bodies can be captured by compact stellar remnants more frequently than it occurs in our Solar System or that minor body formation is relatively easy around millisecond radio pulsars. A peculiar Supernova associated to a GRB may provide an alternative explanation.

It has been recently proposed that giant protoplanets migrating inward through the disc more rapidly than they contract could be tidally disrupted when they fill their Roche lobes ~ 0.1 AU away from their parent protostars. Here we consider the process of mass and angular momentum exchange between the tidally disrupted planet and the surrounding disc in detail. We find that the planet's adiabatic mass-radius relation and its ability to open a deep gap in the disc determine whether the disruption proceeds as a sudden runaway or a balanced quasi-static process. In the latter case the planet feeds the inner disc through its Lagrangian L1 point like a secondary star in a stellar binary system. As the planet loses mass it gains specific angular momentum and normally migrates in the outward direction until the gap closes. Numerical experiments show that planet disruption outbursts are preceded by long "quiescent" periods during which the disc inward of the planet is empty. The hole in the disc is created when the planet opens a deep gap, letting the inner disc to drain onto the star while keeping the outer one stalled behind the planet. We find that the mass-losing planet embedded in a realistic protoplanetary disc spawns an extremely rich set of variability patterns. In a subset of parameter space, there is a limit cycle behaviour caused by non-linear interaction between the planet mass loss and the disc hydrogen ionisation instability. We suggest that tidal disruptions of young massive planets near their stars may be responsible for the observed variability of young accreting protostars such as FU Ori, EXor and T Tauri stars in general.

Unambiguous detection of the tidal disruption of a star would allow an assessment of the presence and masses of supermassive black holes in quiescent galaxies. It would also provide invaluable information on bulge scale stellar processes (such as two-body relaxation) via the rate at which stars are injected into the tidal sphere of influence of the black holes. This rate, in turn, is essential to predict gravitational radiation emission by compact object inspirals. The signature of a tidal disruption event is thought to be a fallback rate for the stellar debris onto the black hole that decreases as $t^{-5/3}$. This mass flux is often assumed to yield a luminous signal that decreases in time at the same rate. In this paper, we calculate the monochromatic lightcurves arising from such an accretion event. Differently from previous studies, we adopt a more realistic description of the fallback rate and of the super-Eddigton accretion physics. We also provide simultaneous lightcurves in optical, UV and X-rays. We show that, after a few months, optical and UV lightcurves scale as $t^{-5/12}$, and are thus substantially flatter than the $t^{-5/3}$ behaviour, which is a prerogative of the bolometric lightcurve, only. At earlier times and for black hole masses $< 10^7 M_{\sun}$, the wind emission dominates: after reaching a peak of $10^{41}-10^{43}$ erg/s at roughly a month, the lightcurve decreases steeply as $\sim t^{-2.6}$, until the disc contribution takes over. The X-ray band, instead, is the best place to detect the $t^{-5/3}$ "smoking gun" behaviour, although it is displayed only for roughly a year, before the emission steepens exponentially.

We examine the tidal disruption event scenario to explain Sw 1644+57, a powerful and persistent X-ray source which suddenly became active as GRB 110328A. The precise localization at the center of a z=0.35 galaxy argues for activity of the central engine as the underlying cause. We look at the suggestion by Bloom et al of the possibility of a tidal disruption event (TDE). We argue that Sw 1644+57 cannot be explained by the traditional TDE model in which the periastron distance is close to the tidal disruption radius - three independent lines of argument indicate the orbit must be deeply plunging or else the powerful jet we are observing could not be produced. These arguments stem from (i) comparing the early X-ray light curve to the expected theoretical fallback rate, (ii) looking at the time of transition to disk-dominated decay, and (iii) considering the TDE rate. Due to the extreme excess in the tidal force above that which would be required minimally to disrupt the star in a deeply plunging orbit at periastron, we suggest this scenario might be referred to more descriptively as a TOE (tidal obliteration event) rather than a TDE.

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

Stimulated by recent results by Meru and Bate (2010a,b), we revisit the issue of resolution requirements for simulating self-gravitating accretion discs with Smoothed Particle Hydrodynamics (SPH). We show that the results by Meru and Bate (2010a) are consistent with those of Meru and Bate (2010b) if they are both interpreted as driven by resolution effects, therefore implying that the resolution criterion for cooling gaseous discs is a function of the imposed cooling rate. We discuss two possible numerical origins of such dependence, which are both consistent with the limited number of available data. Our results tentatively indicate that convergence for current simulations is being reached for a number of SPH particles approaching 10 millions (for a disc mass of order 10 per cent of the central object mass), which would set the critical cooling time for fragmentation at about $15\Omega^{-1}$, roughly a factor two larger than previously thought. More in general, we discuss the extent to which the large number of recent numerical results are reliable or not. We argue that those results that pertain to the dynamics associated with gravitational instabilities (such as the locality of angular momentum transport, and the relationship between density perturbation and induced stress) are robust, while those pertaining to the thermodynamics of the system (such as the determination of the critical cooling time for fragmentation) can be affected by poor resolution.

Semi-analytic models of self-gravitating discs often approximate the angular momentum transport generated by the gravitational instability using the phenomenology of viscosity. This allows the employment of the standard viscous evolution equations, and gives promising results. It is, however, still not clear when such an approximation is appropriate. This paper tests this approximation using high resolution 3D smoothed particle hydrodynamics (SPH) simulations of self-gravitating protostellar discs with radiative transfer. The nature of angular momentum transport associated with the gravitational instability is characterised as a function of both the stellar mass and the disc-to-star mass ratio. The effective viscosity is calculated from the Reynolds and gravitational stresses in the disc. This is then compared to what would be expected if the effective viscosity were determined by assuming local thermodynamic equilibrium or, equivalently, that the local dissipation rate matches the local cooling rate. In general, all the discs considered here settle into a self-regulated state where the heating generated by the gravitational instability is modulated by the local radiative cooling. It is found that low-mass discs can indeed be represented by a local "alpha-parametrisation", provided that the disc aspect ratio is small (H/R < 0.1) which is generally the case when the disc-to-star mass ratio q <0.5. However, this result does not extend to discs with masses approaching that of the central object. These are subject to transient burst events and global wave transport, and the effective viscosity is not well modelled by assuming local thermodynamic equilibrium. In spite of these effects, it is shown that massive (compact) discs can remain stable and not fragment, evolving rapidly to reduce their disc-to-star mass ratios through stellar accretion and radial spreading.

In this paper we present simulated observations of massive self-gravitating circumstellar discs using the Atacama Large Millimetre/sub-millimetre Array (ALMA). Using a smoothed particle hydrodynamics model of a $0.2M_{\odot}$ disc orbiting a $1M_{\odot}$ protostar, with a cooling model appropriate for discs at temperatures below $\sim 160$K and representative dust opacities, we have constructed maps of the expected emission at sub-mm wavelengths. We have then used the CASA ALMA simulator to generate simulated images and visibilities with various array configurations and observation frequencies, taking into account the expected thermal noise and atmospheric opacities. We find that at 345 GHz (870 $\mu$m) spiral structures at a resolution of a few AU should be readily detectable in approximately face-on discs out to distances of the Taurus-Auriga star-forming complex.

In this paper we revisit the issue of the propagation of warps in thin and viscous accretion discs. In this regime warps are know to propagate diffusively, with a diffusion coefficient approximately inversely proportional to the disc viscosity. Previous numerical investigations of this problem (Lodato & Pringle 2007) did not find a good agreement between the numerical results and the predictions of the analytic theories of warp propagation, both in the linear and in the non-linear case. Here, we take advantage of a new, low-memory and highly efficient SPH code to run a large set of very high resolution simulations (up to 20 million SPH particles) of warp propagation, implementing an isotropic disc viscosity in different ways, to investigate the origin of the discrepancy between the theory and the numerical results. Our new and improved analysis now shows a remarkable agreement with the analytic theory both in the linear and in the non-linear regime, in terms of warp diffusion coefficient and precession rate. It is worth noting that the resulting diffusion coefficient is inversely proportional to the disc viscosity only for small amplitude warps and small values of the disc $\alpha$ coefficient ($\alpha < 0.1$). For non-linear warps, the diffusion coefficient is a function of both radius and time, and is significantly smaller than the standard value. Warped accretion discs are present in many contexts, from protostellar discs to accretion discs around supermassive black holes. In all such cases, the exact value of the warp diffusion coefficient may strongly affect the evolution of the system and therefore its careful evaluation is critical in order to correctly estimate the system dynamics (abridged).

In this paper we consider the effects of opacity regimes on the stability of self-gravitating protoplanetary discs to fragmentation into bound objects. Using a self-consistent 1-D viscous disc model, we show that the ratio of local cooling to dynamical timescales Omega*tcool has a strong dependence on the local temperature. We investigate the effects of temperature-dependent cooling functions on the disc's gravitational stability through controlled numerical experiments using an SPH code. We find that such cooling functions raise the susceptibility of discs to fragmentation through the influence of temperature perturbations - the average value of Omega*tcool has to increase to prevent local variability leading to collapse. We find the effects of temperature dependence to be most significant in the "opacity gap" associated with dust sublimation, where the average value of Omega*tcool at fragmentation is increased by over an order of magnitude. We then use this result to predict where protoplanetary discs will fragment into bound objects, in terms of radius and accretion rate. We find that without temperature dependence, for radii < ~10AU a very large accretion rate ~10^-3 Msun/yr is required for fragmentation, but that this is reduced to 10^-4 Msun/yr with temperature-dependent cooling. We also find that the stability of discs with accretion rates < ~10^-7 Msun/yr at radii > ~50AU is enhanced by a lower background temperature if the disc becomes optically thin.

The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.

We investigate the evolution of the relative angle between the stellar rotation axis and the circumstellar disc axis of a star that forms in a stellar cluster from the collapse of a turbulent molecular cloud. This is an inherently chaotic environment with variable accretion, both in terms of rate and the angular momentum of the material, and dynamical interactions between stars. We find that the final stellar rotation axis and disc spin axis can be strongly misaligned, but this occurs primarily when the disc is truncated by a dynamical encounter so that the final disc rotation axis depends simply on what fell in last. This may lead to planetary systems with orbits that are misaligned with the stellar rotation axis, but only if the final disc contains enough mass to form planets. We also investigate the time variability of the inner disc spin axis, which is likely to determine the direction of a protostellar jet. We find that the jet direction varies more strongly for lighter discs, such as those that have been truncated by dynamical interactions or have suffered a period of rapid accretion. Finally, we note that variability of the angular momentum of the material accreting by a star implies that the internal velocity field of such stars may be more complicated than that of aligned differential rotation.

In this Letter we show that if planetesimals form in spiral features in self-gravitating discs, as previously suggested by the idealised simulations of Rice et al, then in realistic protostellar discs, this process will be restricted to the outer regions of the disc (i.e. at radii in excess of several tens of A.U.). This restriction relates to the requirement that dust has to be concentrated in spiral features on a timescale that is less than the (roughly dynamical) lifetime of such features, and that such rapid accumulation requires spiral features whose fractional amplitude is not much less than unity. This in turn requires that the cooling timescale of the gas is relatively short, which restricts the process to the outer disc. We point out that the efficient conversion of a large fraction of the primordial dust in the disc into planetesimals could rescue this material from the well known problem of rapid inward migration at a $\sim$ metre size scale and that in principle the collisional evolution of these objects could help to re-supply small dust to the protostellar disc. We also point out the possible implications of this scenario for the location of planetesimal belts inferred in debris discs around main sequence stars, but stress that further dynamical studies are required in order to establish whether the disc retains a memory of the initial site of planetesimal creation.

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre--existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework. The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high--photon--energy variability on the relatively long orbital timescale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary. We also find that the disc can drive the binary to merger on a reasonable timescale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e. $a_0 \lesssim 0.05$ pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is $\la 0.001$. Gas discs thus appear not to be an effective solution to the `last parsec' problem for a significant class of mergers.

We perform numerical analyses of the structure induced by gravitational instabilities in cooling gaseous accretion discs. For low enough cooling rates a quasi-steady configuration is reached, with the instability saturating at a finite amplitude in a marginally stable disc. We find that the saturation amplitude scales with the inverse square root of the cooling parameter beta = t_cool / t_dyn, which indicates that the heating rate induced by the instability is proportional to the energy density of the induced density waves. We find that at saturation the energy dissipated per dynamical time by weak shocks due is of the order of 20 per cent of the wave energy. From Fourier analysis of the disc structure we find that while the azimuthal wavenumber is roughly constant with radius, the mean radial wavenumber increases with radius, with the dominant mode corresponding to the locally most unstable wavelength. We demonstrate that the density waves excited in relatively low mass discs are always close to co-rotation, deviating from it by approximately 10 per cent. This can be understood in terms of the flow Doppler-shifted phase Mach number -- the pattern speed self-adjusts so that the flow into spiral arms is always sonic. This has profound effects on the degree to which transport through self-gravity can be modelled as a viscous process. Our results thus provide (a) a detailed description of how the self-regulation mechanism is established for low cooling rates, (b) a clarification of the conditions required for describing the transport induced by self-gravity through an effective viscosity, (c) an estimate of the maximum amplitude of the density perturbation before fragmentation occurs, and (d) a simple recipe to estimate the density perturbation in different thermal regimes.

In this paper we revisit the arguments for the basis of the time evolution of the flares expected to arise when a star is disrupted by a supermassive black hole. We present a simple analytic model relating the lightcurve to the internal density structure of the star. We thus show that the standard lightcurve proportional to $t^{-5/3}$ only holds at late times. Close to the peak luminosity the lightcurve is shallower, deviating more strongly from $t^{-5/3}$ for more centrally concentrated (e.g. solar--type) stars. We test our model numerically by simulating the tidal disruption of several stellar models, described by simple polytropic spheres with index $\gamma$. The simulations agree with the analytical model given two considerations. First, the stars are somewhat inflated on reaching pericentre because of the effective reduction of gravity in the tidal field of the black hole. This is well described by a homologous expansion by a factor which becomes smaller as the polytropic index becomes larger. Second, for large polytropic indices wings appear in the tails of the energy distribution, indicating that some material is pushed further away from parabolic orbits by shocks in the tidal tails. In all our simulations, the $t^{-5/3}$ lightcurve is achieved only at late stages. In particular we predict that for solar type stars, this happens only after the luminosity has dropped by at least two magnitudes from the peak. We discuss our results in the light of recent observations of flares in otherwise quiescent galaxies and note the dependence of these results on further parameters, such as the star/hole mass ratio and the stellar orbit.

I review recent progresses in the dynamics and the evolution of self-gravitating accretion discs. Accretion discs are a fundamental component of several astrophysical systems on very diverse scales, and can be found around supermassive black holes in Active Galactic Nuclei (AGN), and also in our Galaxy around stellar mass compact objects and around young stars. Notwithstanding the specific differences arising from such diversity in physical extent, all these systems share a common feature where a central object is fed from the accretion disc, due to the effect of turbulence and disc instabilities, which are able to remove the angular momentum from the gas and allow its accretion. In recent years, it has become increasingly apparent that the gravitational field produced by the disc itself (the disc's self-gravity) is an important ingredient in the models, especially in the context of protostellar discs and of AGN discs. Indeed, it appears that in many cases (and especially in the colder outer parts of the disc) the development of gravitational instabilities can be one of the main agents in the redistribution of angular momentum. In some cases, the instability can be strong enough to lead to the formation of gravitationally bound clumps within the disc, and thus to determine the disc fragmentation. As a result, progress in our understanding of the dynamics of self-gravitating discs is essential to understand the processes that lead to the feeding of both young stars and of supermassive black holes in AGN. At the same time, understanding the fragmentation conditions is important to determine under which conditions AGN discs would fragment and form stars and whether protostellar discs might form giant gaseous planets through disc fragmentation.

We investigate the orbital evolution of planetesimals in a self-gravitating circumstellar disc in the size regime ($\sim 1-5000$ km) where the planetesimals behave approximately as test particles in the disc's non-axisymmetric potential. We find that the particles respond to the stochastic, regenerative spiral features in the disc by executing large random excursions (up to a factor of two in radius in $\sim 1000$ years), although typical random orbital velocities are of order one tenth of the Keplerian speed. The limited time frame and small number of planetesimals modeled does not permit us to discern any {\it net} direction of planetesimal migration. Our chief conclusion is that the high eccentricities ($\sim 0.1$) induced by interaction with spiral features in the disc is likely to be highly unfavourable to the collisional growth of planetesimals in this size range while the disc is in the self-gravitating regi me. Thus {\it if}, as recently argued by Rice et al 2004, 2006, the production of planetesimals gets under way when the disc is in the self-gravitating regime (either at smaller planetesimal size scales, where gas drag is important, or via gravitational fragmentation of the solid component), then the planetesimals thus produced would not be able to grow collisionally until the disc ceased to be self-gravitating. It is unclear, however,given the large amplitude excursions undergone by planetesimals in the self-gravitating disc, whether they would be retained in the disc throughout this period, or whether they would instead be lost to the central star.

We investigate the evolution of high redshift seed black hole masses at late times and their observational signatures. The massive black hole seeds studied here form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this way, we follow the mass assembly of these black holes to the present time using a Monte-Carlo merger tree. Using this machinery we predict the black hole mass function at high redshifts and at the present time; the integrated mass density of black holes and the luminosity function of accreting black holes as a function of redshift. These predictions are made for a set of three seed models with varying black hole formation efficiency. Given the accuracy of current observational constraints, all 3 models can be adequately fit. Discrimination between the models appears predominantly at the low mass end of the present day black hole mass function which is not observationally well constrained. However, all our models predict that low surface brightness, bulgeless galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at z=0. This effect is more pronounced for low mass galaxies. This is the key discriminant between the models studied here and the Population III remnant seed model. We find that there exists a population of low mass galaxies that do not host nuclear black holes. Our prediction of the shape of the black hole mass - velocity dispersion relation at the low mass end is in agreement with the recent observational determination from the census of low mass galaxies in the Virgo cluster.

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the sub-stellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarfs of mass 0.05M_{\odot}, our models predict a maximum planetary mass of ~5M_{\oplus}, orbiting with semi-major axis ~1AU. However, we note that the predictions for the mass - semi-major axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the sub-stellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasise the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type-II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.

In this paper we explore numerically the evolution of a warped accretion disc. Here, we focus here on the regime where the warp evolves diffusively. By comparing the numerical results to a simple diffusion model, we are able to determine the diffusion coefficient of the warp, $\alpha_2$, as a function of the relevant disc parameters. We find that while in general the disc behaviour is well reproduced by the diffusion model and for relatively large viscosities the warp diffusion is well described by the linear theory (in particular confirming that the warp diffusion coefficient is inversely proportional to viscosity), significant non-linear effects are present as the viscosity becomes smaller, but still dominates over wave-propagation effects. In particular, we find that the inverse dependence of the diffusion coefficient on viscosity breaks down at low viscosities, so that $\alpha_2$ never becomes larger than a saturation value $\alpha_{\rm max}$ of order unity. This can have major consequences in the evolution of systems where a warped disc is present. In particular, it affects the location of the warp radius in the Bardeen-Petterson effect and therefore the spin up (or spin down) of supermassive black holes in the nuclei of galaxies. Additionally, we also find that while the rate of warp diffusion does not depend significantly on the detailed viscosity formulation, the rate of internal precession generated by the warp is strongly affected by it. Such effects should be considered with care when modeling the evolution of warped discs. This emphasises the need to test the above results using different numerical schemes, and with higher resolution, in order to investigate the degree to which numerical simulations are able to provide accurate modeling of the complex fluid dynamics of warped discs. (Abridged)

A number of previous studies of the fragmentation of self-gravitating protostellar discs have modeled radiative cooling with a cooling timescale (t_{cool}) parameterised as a simple multiple (beta_{cool}) of the local dynamical timescale. Such studies have delineated the `fragmentation boundary' in terms of a critical value of beta_{cool} (beta_{crit}), where the disc fragments if beta_{cool} < beta_{crit}. Such an approach however begs the question of how in reality a disc could ever be assembled with beta_{cool} < beta_{crit}. Here we adopt the more realistic approach of gradually reducing beta_{cool}, as might correspond to changes in thermal regime due to secular changes in the disc density profile. We find that when beta_{cool} is gradually reduced (on a timescale longer than t_{cool}), the disc is stabilised against fragmentation, compared with models in which beta_{cool} is reduced rapidly. We therefore conclude that a disc's ability to remain in a self-regulated, self-gravitating state (without fragmentation) is partly dependent on its thermal history, as well as its current cooling rate. Nevertheless, a slow reduction in t_{cool} appears only to lower the fragmentation boundary by about a factor two in t_{cool} and thus only permits maximum alpha values (parameterising the efficiency of angular momentum transfer in the disc) that are about a factor two higher than determined hitherto. Our results therefore do not undermine the notion of a fundamental upper limit to the heating rate that can be delivered by gravitational instabilities before the disc is subject to fragmentation. An important implication of this work, therefore, is that self-gravitating discs can enter into the regime of fragmentation via secular evolution and it is not necessary to invoke rapid (impulsive) events to trigger fragmentation.

In this paper we derive the mass function of seed black holes that result from the central mass concentrated via disc accretion in collapsed haloes at redshift $z\approx 15$. Using standard arguments including stability, we show that these pre-galactic discs can assemble a significant mass concentration in the inner regions, providing fuel for the formation and initial growth of super-massive black holes. Assuming that these mass concentrations do result in central seed black holes, we determine the mass distribution of these seeds as a function of key halo properties. The seed mass distribution determined here turns out to be asymmetric and skewed to higher masses. Starting with these initial seeds, building up to $10^9$ solar masses by $z = 6$ to power the bright quasars is not a problem in the standard LCDM cosmogony. These seed black holes in gas rich environments are likely to grow into the supermassive black holes at later times via mergers and accretion. Gas accretion onto these seeds at high redshift will produce miniquasars that likely play an important role in the reionization of the Universe. Some of these seed black holes on the other hand could be wandering in galaxy haloes as a consequence of frequent mergers, powering the off-nuclear ultra-luminous X-ray sources detected in nearby galaxies.

In this paper, we use high-resolution smoothed particle hydrodynamics (SPH) simulations to investigate the response of a marginally stable self-gravitating protostellar disc to a close parabolic encounter with a companion discless star. Our main aim is to test whether close brown dwarfs or massive planets can form out of the fragmentation of such discs. We follow the thermal evolution of the disc by including the effects of heating due to compression and shocks and a simple prescription for cooling and find results that contrast with previous isothermal simulations. In the present case we find that fragmentation is inhibited by the interaction, due to the strong effect of tidal heating, which results in a strong stabilization of the disc. A similar behaviour was also previously observed in other simulations involving discs in binary systems. As in the case of isolated discs, it appears that the condition for fragmentation ultimately depends on the cooling rate.

The spectral energy distributions (SEDs) of some T Tauri stars display a deficit of near-IR flux that could be a consequence of an embedded Jupiter-mass planet partially clearing an inner hole in the circumstellar disc. Here, we use two-dimensional numerical simulations of the planet-disc interaction, in concert with simple models for the dust dynamics, to quantify how a planet influences the dust at different radii within the disc. We show that pressure gradients at the outer edge of the gap cleared by the planet act as a filter - letting particles smaller than a critical size through to the inner disc while holding back larger particles in the outer disc. The critical particle size depends upon the disc properties, but is typically of the order of 10 microns. This filtration process will lead to discontinuous grain populations across the planet's orbital radius, with small grains in the inner disc and an outer population of larger grains. We show that this type of dust population is qualitatively consistent with SED modelling of systems that have optically thin inner holes in their circumstellar discs. This process can also produce a very large gas-to-dust ratio in the inner disc, potentially explaining those systems with optically thin inner cavities that still have relatively high accretion rates.

An unsolved issue in the standard core accretion model for gaseous planet formation is how kilometre-sized planetesimals form from, initially, micron-sized dust grains. Solid growth beyond metre sizes can be difficult both because the sticking efficiency becomes very small, and because these particles should rapidly migrate into the central star. We consider here how metre-sized particles evolve in self-gravitating accretion discs using simulations in which the gravitational influence of the solid particles is also included. Metre-sized particles become strongly concentrated in the spiral structures present in the disc and, if the solid to gas density ratio is sufficiently high, can fragment due to their own self-gravity to form planetesimals directly. This result suggests that planetesimal formation may occur very early in the star formation process while discs are still massive enough to be self-gravitating. The dependence of this process on the surface density of the solids is also consistent with the observation that extrasolar planets are preferentially found around high metallicity stars.

In this paper we discuss the evolution of gravitationally unstable pre-galactic discs that result from the collapse of haloes at high redshift $z \approx 10$ or so, which have not yet been enriched by metals. In cases where molecular hydrogen formation is suppressed the discs are maintained at a temperature of a few thousand degrees Kelvin. However, when molecular hydrogen is present cooling can proceed down to a few hundred degrees Kelvin. Analogous to the case of the larger scale proto-galactic discs, we assume that the evolution of these discs is mainly driven by angular momentum redistribution induced by the development of gravitational instabilities in the disc. We also properly take into account the possibility of disc fragmentation. We thus show that this simple model naturally predicts the formation of supermassive black holes in the nuclei of such discs and provides a robust determination of their mass distribution as a function of halo properties. We estimate that roughly 5% of discs resulting from the collapse of haloes with $M\approx 10^7 M_{\odot}$ should host a massive black hole with a mass $M_{\rm BH}\approx 10^5 M_{\odot}$. We confirm our arguments with time-dependent calculations of the evolution of the surface density and of the accretion rate in these primordial discs. This mechanism offers an efficient way to form seed black holes at high redshift. The predicted masses for our black hole seeds enable the comfortable assembly of $10^9 M_{\odot}$ black holes powering the luminous quasars detected by the Sloan Digital Sky Survey at $z = 6$ for a concordance cosmology. (abridged)

Disks are a natural outcome of the star formation process in which they play a crucial role. Luminous, massive stars of spectral type earlier than B4 are likely to be those that benefit most from the existence of accretion disks, which may significantly reduce the effect of radiation pressure on the accreting material. The scope of the present contribution is to review the current knowledge about disks in young high-mass (proto)stars and discuss their implications. The issues of disk stability and lifetime are also discussed. We conclude that for protostars of less than ~20 solar masses, disks with mass comparable to that of the central star are common. Above this limit the situation is unclear and there are no good examples of proto O4-O8 stars surrounded by accretion disks: in these objects only huge, massive, toroidal, non-equilibrium rotating structures are seen. It is clear on the other hand that the observed disks in stars of 10-20 solar masses are likely to be unstable and with short lifetimes.

In this paper we consider the process of alignment of a spinning black hole and a surrounding misaligned accretion disc. We use a simplified set of equations, that describe the evolution of the system in the case where the propagation of warping disturbances in the accretion disc occurs diffusively, a situation likely to be common in the thin discs in Active Galactic Nuclei (AGN). We also allow the direction of the hole spin to move under the action of the disc torques. In such a way, the evolution of the hole-disc system is computed self-consistently. We consider a number of different situations and we explore the relevant parameter range, by varying the location of the warp radius $R_{\rm w}$ and the propagation speed of the warp. We find that the dissipation associated with the twisting of the disc results in a large increase in the accretion rate through the disc, so that AGN accreting from a misaligned disc are likely to be significantly more luminous than those accreting from a flat disc. We compute explicitly the time-scales for the warping of the disc and for the alignment process and compare our results with earlier estimates based on simplified steady-state solutions. We also confirm earlier predictions that, under appropriate circumstances, accretion can proceed in a counter-aligned fashion, so that the accreted material will spin-down the hole, rather than spinning it up. Our results have implication in a number of different observational features of AGN such as the orientation and shape of jets, the shape of X-ray iron lines, and the possibility of obscuration and absorption of X-ray by the outer disc as well as the general issue of the spin history of black holes during their growth.

In this paper we discuss the nature and the possible formation scenarios of the companion of the brown dwarf 2MASS 1207334-393254. We initially discuss the basic physical properties of this object and conclude that, although from its absolute mass ($5M_{\rm Jup}$), it is a planetary object, in terms of its mass ratio $q$ and of its separation $a$ with respect to the primary brown dwarf, it is consistent with the statistical properties of binaries with higher primary mass. We then explore the possible formation mechanism for this object. We show that the standard planet formation mechanism of core accretion is far too slow to form this object within 10 Myr, the observed age of the system. On the other hand, the alternative mechanism of gravitational instability (proposed both in the context of planet and of binary formation) may, in principle, work and form a system with the observed properties.

The issue of fragmentation in self-gravitating gaseous accretion discs has implications both for the formation of stars in discs in the nuclei of active galaxies, and for the formation of gaseous planets or brown dwarfs in circumstellar discs. It is now well established that fragmentation occurs if the disc is cooled on a timescale smaller than the local dynamical timescale, while for longer cooling times the disc reaches a quasi-steady state in thermal equilibrium, with the cooling rate balanced by the heating due to gravitational stresses. We investigate here how the fragmentation boundary depends on the assumed equation of state. We find that the cooling time required for fragmentation increases as the specific heat ratio, gamma, decreases, exceeding the local dynamical timescale for gamma = 7/5. This result can be easily interpreted as a consequence of there being a maximum stress (in units of the local disc pressure) that can be sustained by a self-gravitating disc in quasi-equilibrium. Fragmentation occurs if the cooling time is such that the stress required to reach thermal equilibrium exceeds this value, independent of gamma. This result suggest that a quasi-steady, self-gravitating disc can never produce a stress that results in the viscous alpha parameter exceeding ~ 0.06.

We present the results of a photometric monitoring campaign of three well studied FU Orionis systems (FU Orionis, V1057 Cygni and V1515 Cygni) undertaken at Maidanak Observatory between 1981 and 2003. When combined with photometric data in the literature, this database provides a valuable resource for searching for short timescale variability - both periodic and aperiodic - as well as for studying the secular evolution of these systems. In the case of V1057 Cyg (which is the system exhibiting the largest changes in brightness since it went into outburst) we compare the photometric data with time dependent models. We show that prior to the end of the `plateau' stage in 1996, the evolution of V1057 Cyg in the $V-(B-V)$ colour-magnitude diagram is well represented by disc instability models in which the outburst is triggered by some agent - such as an orbiting planet - in the inner disc. Following the end of the plateau phase in 1996, the dimming and irregular variations are consistent with occultation of the source by a variable dust screen, which has previously been interpreted in terms of dust condensation events in the observed disc wind. Here we instead suggest that this effect results from the interaction between the wind and an infalling dusty envelope, the existence of this envelope having been previously invoked in order to explain the mid infrared emissio of FU Orionis systems. We discuss how this model may explain some of the photometric and spectroscopic characteristics of FU Orionis systems in general.

In this paper, we extend our previous analysis (Lodato & Rice 2004) of the transport properties induced by gravitational instabilities in cooling, gaseous accretion discs to the case where the disc mass is comparable to the central object. In order to do so, we have performed global, three-dimensional smoothed particle hydrodynamics simulations of massive discs. These new simulations show a much more complex temporal evolution with respect to the less massive case. Whereas in the low disc mass case a self-regulated, marginally stable state (characterized by an approximately constant radial profile of the stability parameter $Q$) is easily established, in the high disc mass case we observe the development of an initial transient and subsequent settling down in a self-regulated state in some simulations, or a series or recurrent spiral episodes, with low azimuthal wave number $m$, in others. Accretion in this last case can therefore be a highly variable process. On the other hand, we find that the secular evolution of the disc is relatively slow. In fact, the time-average of the stress induced by self-gravity results in accretion time-scales much longer than the dynamical timescale, in contrast with previous isothermal simulations of massive accretion discs. We have also compared the resulting stress tensor with the expectations based on a local theory of transport, finding no significant evidence for global wave energy transport.

In this paper we consider the evolution of small planetesimals in marginally stable, self-gravitating protoplanetary discs. The drag force between the disc gas and the embedded planetesimals generally causes the planetesimals to drift inwards through the disc at a rate that depends on the particle size. In a marginally stable, self-gravitating disc, however, the planetesimals are significantly influenced by the non-axisymmetric spiral structures resulting from the growth of the gravitational instability. For small particles, that are strongly coupled to the disc gas, and for large particles, that have essentially decoupled from the disc gas, the effect is not particularly significant. Intermediate sized particles, which would generally have the largest radial drift rates, do, however, become significantly concentrated at the peaks of the spiral arms. These high density regions may persist for, of order, an orbital period and may attain densities comparable to that of the disc gas. Although at the end of the simulation only ~ 25 % of the planetesimal particles lie in regions of enhanced density, during the course of the simulation at least 75 % of the planetesimal particles have at some stage been in a such a region. We find that the concentration of particles in the spiral arms results in an increased collision rate, an effect that could significantly accelerate planetesimal growth. The density enhancements may also be sufficient for the growth of planetesimals through direct gravitational collapse. The interaction between small planetesimals and self-gravitating spiral structures may therefore play an important role in the formation of large planetesimals that will ultimately coagulate to form terrestrial planets or the cores of gas/ice giant planets.

FU Orionis are young stellar objects undergoing episodes of enhanced luminosity, which are generally ascribed to a sudden increase of mass accretion rate in the surrounding protostellar disc. Models invoking a thermal instability in the disc are able to reproduce many features of the outburst, but cannot explain the rapid rise time-scale observed in many cases. Here we explore the possibility (originally suggested by Clarke & Syer 1996) that the thermal instability is triggered away from the disc inner edge (at a distance of $\approx 10R_{\odot}$ from the central protostar) due to the presence of a massive planet embedded in the disc. We have constructed simple, one-dimensional, time-dependent models of the disc evolution, taking into account both the interaction between the disc and the planet, and the thermal evolution of the disc. We are indeed able to reproduce rapid rise outbursts (with rise time-scale $\approx 1$ yr), with a planet mass $M_{\mathrm{s}}=10-15M_{\mathrm{Jupiter}}$. We show that the luminosity and the duration of the outbursts are increasing functions of planet mass. We also show that the inward migration of the planet is significantly slowed once it reaches the radius where it is able to trigger FU Orionis outbursts, thus suggesting that a single planet may be involved in triggering several outbursts.

In this paper we examine the issue of characterising the transport associated with gravitational instabilities in relatively cold discs, discussing in particular the conditions under which it can be described within a local, viscous framework. We present the results of global, three-dimensional, SPH simulations of self-gravitating accretion discs, in which the disc is cooled using a simple parametrisation for the cooling function. Our simulations show that the disc settles in a ``self-regulated'' state, where the axisymmetric stability parameter $Q\approx 1$ and where transport and energy dissipation are dominated by self-gravity. We have computed the gravitational stress tensor and compared our results with expectations based on a local theory of transport. We find that, as long as the disc mass is smaller than $0.25M_{\star}$ and the aspect ratio $H/R\lesssim 0.1$, transport is determined locally, thus allowing for a viscous treatment of the disc evolution.

In this paper we examine the issue of characterizing the transport associated with gravitational instabilities in relatively cold discs, discussing in particular under which condition it can be described within a local, viscous framework. We present the results of global, three-dimensional, SPH simulations of self-gravitating accretion discs, in which the disc is cooled using a simple parameterization for the cooling function. Our simulations show that the disc settles in a ``self-regulated'' state, where the axisymmetric stability parameter $Q\approx 1$ and where transport and energy dissipation are dominated by self-gravity. We have computed the gravitational stress tensor and compared our results with expectations based on a local theory of transport. We find that, for disc masses smaller than $0.25M_{\star}$ and aspect ratio $H/r\lesssim 0.1$, transport is determined locally, thus allowing for a viscous treatment of the disc evolution.

Studies of the Spectral Energy Distribution of Young Stellar Objects suggest that the outer disk of FU Orionis objects might be self-gravitating. In this paper we propose a method to test directly whether, in these objects, significant deviations from Keplerian rotation occur. In a first approach, we have used a simplified model of the disk vertical structure that allows us to quickly bring out effects related to the disk self-gravity. We find that the often studied optical and near-infrared line profiles are produced too close to the central object to provide significant evidence for non-Keplerian rotation. Based on parameters relevant for the case of FU Ori, we show that high-resolution long-wavelength spectroscopy, of the far-infrared H$_2$ pure rotational lines (sometimes observed in ``passive'' protostellar disks) and sub-mm CO lines, should be well suited to probe the rotation curve in the outer disk, thus measuring to what extent it is affected by the disk self-gravity. The results of the present exploratory paper should be extended soon to a more realistic treatment of the disk vertical structure.

The nucleus of the Seyfert galaxy NGC 1068 is believed to host a supermassive black hole. Evidence for the presence of a massive central object is provided by water maser emission, which displays a linear pattern in the sky, suggestive of a rotating disk. The rotating disk hypothesis is further strengthened by the declining shape of the derived rotation curve. Similar maser emission from NGC 4258 has led to a reliable estimate of the mass of the central black hole, because in this case the rotation curve is Keplerian. In the case of NGC 1068 the rotation curve traced by the water maser is non-Keplerian. In this paper we provide an interpretation of the non-Keplerian rotation in NGC 1068 by means of a self-gravitating accretion disk model. We obtain a good fit to the available data and derive a black hole mass M_{bh}=(8.0\pm 0.3) 10^6M_{sun}. The resulting disk mass is comparable to the black hole mass. As an interesting by-product of our fitting procedure, we are able to estimate the viscosity parameter, which turns out to be \alpha\approx 10^{-2}, in line with some theoretical expectations.

The long wavelength emission of protostellar objects is commonly attributed to a disk of gas and dust around the central protostar. In the first stages of disk accretion or in the case of high mass protostars, the disk mass is likely to be sufficiently large, so that the disk self-gravity may have an impact on the dynamics and the emission properties of the disk. In this paper we describe the spectral energy distribution (SED) produced by a simple, non-flaring, self-gravitating accretion disk model. Self-gravity is included in the calculation of the rotation curve of the disk and in the energy balance equation, as a term of effective heating related to Jeans instability. In order to quantify in detail the requirements on the mass of the disk and on the accretion rate posed on the models by realistic situations, we compare the SEDs produced by these models with the observed SEDs of a small sample of well-studied protostellar objects. We find that relatively modest disks - even lighter than the central star - can lead to an interesting fit to the infrared SED of the FU Orionis objects considered, while in the case of T Tauri stars the required parameters fall outside the range suggested as acceptable by the general theoretical and observational scenario. On the basis of the present results, we may conclude that the contribution of a self-gravitating disk is plausible in several cases (in particular, for FU Orionis objects) and that, in the standard irradiation dominated disk scenario, it would help softening the requirements encountered by Keplerian accretion models.

In the dynamics of accretion disks, the presence of collective effects associated with the self-gravity of the disk is expected to affect not only the momentum transport, but also the relevant energy balance equations, which could differ substantially from the non-self-gravitating case. Here we follow the picture that, when the disk is sufficiently cold, the stirring due to Jeans-related instabilities acts as a source of effective heating. The corresponding reformulation of the energy equations allows us to: ({\it i}) demonstrate how self-regulation can be established, so that the stability parameter $Q$ is maintained close to a threshold value, with weak dependence on radius; ({\it ii}) rediscuss the opacity properties in the self-gravitating regime. In particular, we show that, if cooling is dominated by {\it bremsstrahlung}, an optically thin stationary accretion solution is thermally stable, even in the non-advective case, provided the disk is self-gravitating. The details of the cooling function have little effect on the structure of such accretion disk, which is in any case brought by self-gravity to self-regulate. This condition of self-gravitating accretion is expected to be appropriate for the outer regions of many disks of astrophysical interest. With the reformulation of the energy equations described in this paper we have also secured: ({\it iii}) a starting point for the study of the emission properties of self-gravitating accretion disks; ({\it iv}) a tool to analyze the structure of the transition region, where the disk becomes self-gravitating.

We consider a class of steady-state self-gravitating accretion disks for which efficient cooling mechanisms are assumed to operate so that the disk is self-regulated at a condition of approximate marginal Jeans stability. In an earlier paper, this scenario had been shown to lead naturally, in the absence of a central point mass, to a self-similar solution characterized by a flat rotation curve. In this article we investigate the entire parameter space available for such self-regulated accretion disks and provide two non-trivial extensions of the model. The first extension is that of a bimodal disk, obtained by partially relaxing the self-regulation constraint, so that full matching with an inner "standard" Keplerian accretion disk takes place. The second extension is the construction of self-regulated accretion disks embedded in a diffuse spherical "halo". The analysis is further strengthened by a careful discussion of the vertical structure of the disk, in such a way that the transition from self-gravity dominated to non-gravitating disks is covered uniformly.