We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale.

We aim to unveil the observational imprint of physical mechanisms that govern planetary formation in the young, multiple system GG Tau A. We present ALMA observations of $^{12}$CO and $^{13}$CO 3-2 and 0.9 mm continuum emission with 0.35" resolution. The $^{12}$CO 3-2 emission, found within the cavity of the circumternary dust ring (at radius $< 180$ AU) where no $^{13}$CO emission is detected, confirms the presence of CO gas near the circumstellar disk of GG Tau Aa. The outer disk and the recently detected hot spot lying at the outer edge of the dust ring are mapped both in $^{12}$CO and $^{13}$CO. The gas emission in the outer disk can be radially decomposed as a series of slightly overlapping Gaussian rings, suggesting the presence of unresolved gaps or dips. The dip closest to the disk center lies at a radius very close to the hot spot location at $\sim250-260$~AU. The CO excitation conditions indicate that the outer disk remains in the shadow of the ring. The hot spot probably results from local heating processes. The two latter points reinforce the hypothesis that the hot spot is created by an embedded proto-planet shepherding the outer disk.

Mass models of galactic disks traditionally rely on axisymmetric density and rotation curves, paradoxically acting as if their most remarkable asymmetric features, such as lopsidedness or spiral arms, were not important. In this article, we relax the axisymmetry approximation and introduce a methodology that derives 3D gravitational potentials of disk-like objects and robustly estimates the impacts of asymmetries on circular velocities in the disk midplane. Mass distribution models can then be directly fitted to asymmetric line-of-sight velocity fields. Applied to the grand-design spiral M99, the new strategy shows that circular velocities are highly nonuniform, particularly in the inner disk of the galaxy, as a natural response to the perturbed gravitational potential of luminous matter. A cuspy inner density profile of dark matter is found in M99, in the usual case where luminous and dark matter share the same center. The impact of the velocity nonuniformity is to make the inner profile less steep, although the density remains cuspy. On another hand, a model where the halo is core dominated and shifted by 2.2-2.5 kpc from the luminous mass center is more appropriate to explain most of the kinematical lopsidedness evidenced in the velocity field of M99. However, the gravitational potential of luminous baryons is not asymmetric enough to explain the kinematical lopsidedness of the innermost regions, irrespective of the density shape of dark matter. This discrepancy points out the necessity of an additional dynamical process in these regions: possibly a lopsided distribution of dark matter.

We report a simple method to generate potential/surface density pairs in flat axially symmetric finite size disks. Potential/surface density pairs consist of a ``homogeneous'' pair (a closed form expression) corresponding to a uniform disk, and a ``residual'' pair. This residual component is converted into an infinite series of integrals over the radial extent of the disk. For a certain class of surface density distributions (like power laws of the radius), this series is fully analytical. The extraction of the homogeneous pair is equivalent to a convergence acceleration technique, in a matematical sense. In the case of power law distributions, the convergence rate of the residual series is shown to be cubic inside the source. As a consequence, very accurate potential values are obtained by low order truncation of the series. At zero order, relative errors on potential values do not exceed a few percent typically, and scale with the order N of truncation as 1/N**3. This method is superior to the classical multipole expansion whose very slow convergence is often critical for most practical applications.

We report an analytical expression for the locations of Lindblad resonances induced by a perturbing protoplanet, including the effect of disk gravity. Inner, outer and differential torques are found to be enhanced compared to situations where a keplerian velocity field for the dynamics of both the disk and the planet is assumed. Inward migration is strongly accelerated when the disk gravity is only accounted for in the planet orbital motion. The addition of disk self-gravity slows down the planet drift but not enough to stop it.

We demonstrate the high accuracy of the density splitting method to compute the gravitational potential and field in the plane of razor-thin, axially symmetric discs, as preliminarily outlined in Pierens & Hure (2004). Because residual kernels in Poisson integrals are not C^infinity-class functions, we use a dynamical space mapping in order to increase the efficiency of advanced quadrature schemes. In terms of accuracy, results are better by orders of magnitude than for the classical FFT-methods.

Structure of a quasi-stationary stellar cluster is modelled assuming that it is embedded in the gravitational field of a super-massive black hole. Gradual orbital decay of stellar trajectories is caused by the dissipative interaction with an accretion disc. Gravitational field of the disc is constructed and its effect on the cluster structure is taken into account as an axially symmetric perturbation. Attention is focused on a circumnuclear region (r<10^4 gravitational radii) where the effects of the central black hole and the disc dominate over the influence of an outer galaxy. It is shown how the stellar system becomes gradually flattened towards the disc plane. For certain combinations of the model parameters, a toroidal structure is formed by a fraction of stars. Growing anisotropy of stellar velocities as well as their segregation occur. The mass function of the inner cluster is modified and it progressively departs from the asymptotic form assumed in the outer cluster. A new stationary distribution can be characterized in terms of velocity dispersion of the stellar sample in the central region of the modified cluster.

We investigate the analogy between circumstellar disks and the Taylor-Couette flow. Using the Reynolds similarity principle, the analogy results in a number of parameter-free predictions about stability of the disks, and their turbulent transport properties, provided the disk structure is available. We discuss how the latter can be deduced from interferometric observations of circumstellar material. We use the resulting disk structure to compute the molecular transport coefficients, including the effect of ionization by the central object. The resulting control parameter indicates that the disk is well into the turbulent regime. The analogy is also used to compute the effective accretion rate, as a function of the disk characteristic parameters (orbiting velocity, temperature and density). These values are in very good agreement with experimental, parameter-free predictions derived from the analogy. The turbulent viscosity is also computed and found to correspond to an $\alpha$-parameter $2\times 10^{-4}<\alpha<2\times 10^{-2}$. Predictions regarding fluctuations are also checked: luminosity fluctuations in disks do obey the same universal distribution as energy fluctuations observed in a laboratory turbulent flow. Radial velocity dispersion in the outer part of the disk is predicted to be of the order of 0.1 km/s, in agreement with available observations. All these issues provide a proof of the turbulent character of the circumstellar disks, as well as a parameter-free theoretical estimate of effective accretion rates.

Fluid discs and tori around black holes are discussed within different approaches and with the emphasis on the role of disc gravity. First reviewed are the prospects of investigating the gravitational field of a black hole--disc system by analytical solutions of stationary, axially symmetric Einstein's equations. Then, more detailed considerations are focused to middle and outer parts of extended disc-like configurations where relativistic effects are small and the Newtonian description is adequate. Within general relativity, only a static case has been analysed in detail. Results are often very inspiring, however, simplifying assumptions must be imposed: ad hoc profiles of the disc density are commonly assumed and the effects of frame-dragging and completely lacking. Astrophysical discs (e.g. accretion discs in active galactic nuclei) typically extend far beyond the relativistic domain and are fairly diluted. However, self-gravity is still essential for their structure and evolution, as well as for their radiation emission and the impact on the environment around. For example, a nuclear star cluster in a galactic centre may bear various imprints of mutual star--disc interactions, which can be recognised in observational properties, such as the relation between the central mass and stellar velocity dispersion.

Self-gravitating accretion discs have only been studied in a few nearby objects using maser spots at the parsec-scale. We find a new spectral window for observing the self-gravitating accretion disc in super-Eddington accreting Active Galactic Nuclei (AGNs). This window is determined by calculating the outermost radius (rsg) of a non self-gravitating disc and the corresponding emission wavelength (lsg) as a function of various disc parameters. We find that lsg reaches ~ 4000 AA for alpha=0.1, when Mdot > 70 (M_BH / 10^7 M_Sun)^-1 L_Edd/c^2 (where alpha, Mdot, M_BH and L_Edd are, respectively, the viscosity parameter, gas accretion rate onto the central black hole (BH), the BH mass and the Eddington luminosity). Moreover, lsg is as small as 1500 AA for alpha = 0.001, which is the smallest alpha case in this study. Therefore, the window for observing the self-gravitating part of an AGN accretion disc is from 2 micron to lsg. Incidentally, rsg can be less than the photon trapping radius for Mdot >~ 10^3.3 L_Edd/c^2. Namely, a self-gravitating, optically-thick, advection-dominated accretion disc is expected to appear in the extremely high accretion rate regime. Next, we demonstrate that the Mid-Infrared to X-ray spectrum of a bright, well-studied Narrow-Line Seyfert 1 galaxy, Ton S 180, is indeed well fitted by the spectrum arising from the following three components: an inner slim disc (with a corona), an outer, self-gravitating non-Keplerian disc and a dusty torus. The total mass, BH mass plus the entire disc mass, is found to be about (1.4 - 8.0) M_BH. If the surface density varies with radius r in proportion to r^-0.6, the total mass is consistent with the central mass estimated by H_beta and [O III] widths.

For a few nearby active nuclei, the disc orbiting the black hole is traced by water maser emission. By combining a simple model together with observed velocity profiles, we show that it is possible to put constraints on the black hole mass and on the distribution of matter (shape, density, size) in the outer disc. We then report possible parameters for the non-keplerian disc and for the black hole in NGC 1068, and mention an uncertainty of at least 25% on the central mass in NGC 4258.

We demonstrate that the ``sub-keplerian'' rotation curve of maser spots in NGC 1068 can be explained by the gravitational attraction of the disc orbiting the central black hole. Possible parameters matching observations are: black hole mass of 12 million solar masses, disc outer edge > 1.3 pc, aspect ratio in the range 0.003 to 0.3, surface density varying approximately as 1/R, and disc mass of about 9.4 million solar masses. The physical conditions required for the excitation of masers are fulfilled, and the outer disc would stand in a gravitationally marginally stable state.

Current layered accretion models neglect the properties of the ``dead zone''. However, as argued here from simple considerations, the thickness of this zone is a critical quantity when the disc is in hydrostatic equilibrium. It controls not only the structure of the superficial, active layers, but also the mid-plane density and the total disc mass, and should therefore be introduced in models of that kind, steady or not. But in the absence of intrinsic heating, the dead zone must have a tiny size which, given the non-stationary and turbulent character of the global flow, makes very likely its mixing together with the two active layers.

We address the question of the relations between the black hole's mass, the accretion rate, the bolometric luminosity, the optical luminosity and the size of the Broad Line Region (BLR) in Active Galactic Nuclei, using recent observational data obtained from monitoring campaigns. We show that a standard accretion disc cannot account for the observed optical luminosity, unless it radiates at super-Eddington rates. This implies the existence of another, dominant emission mechanism in the optical range, or a non standard disc (non stationary, ADAF and/or strong outflows). Narrow Line Seyfert 1 galaxies (NLS1s) are most extreme in this context: they have larger bolometric to Eddington luminosity ratios than Broad Line Seyfert 1 (BLS1s), and most likely a larger "non disc" component in the optical range. From realistic simulations of self-gravitating $\alpha$-discs, we have systematically localized the gravitationally unstable disc and shown that, given uncertainties on both the model and observations, it coincides quite well with the size of the BLR. We therefore suggest that the gravitationally unstable disc is the source which releases BLR clouds in the medium. However the influence of the ionization parameter is also required to explain the correlation found between the size of the BLR and the luminosity. In this picture the size of the BLR in NLS1s (relative to the black hole size) is larger (and the emission line width smaller) than in BLS1s simply because their Eddington ratio is larger.

The equations governing the vertical structure of a stationary keplerian accretion disc are presented. The model is based on the alpha-viscosity, includes self-gravity, convective transport and turbulent pressure. A few properties of the model are discussed for circumstellar and AGN discs. We show the strong sensitivity of the disc structure to the viscous energy deposition towards the vertical axis, specially when entering inside the self-gravitating part of the disc. The local version of the alpha-prescription leads to a "singular" behavior which is also predicted by the vertically averaged model. With respect, a much softer transition is observed with the "alpha-P" formalism. Turbulent pressure is important only for alpha > 0.1. It lowers vertical density gradients, significantly thickens the disc, tends to wash out density inversions and pushes the self-gravitating region to slightly larger radii. Curves localizing the inner edge of the self-gravitating disc as functions of the viscosity parameter and accretion rate are given. The lower alpha, the closer to the center the self-gravitating regime, and the sensitivity to the accretion rate is generally weak, except for alpha < 0.1. This study suggests that models aiming to describe T-Tauri discs beyond about a few to a few tens astronomical units from the central protostar using the alpha-theory should consider vertical self-gravity. The Primitive Solar Nebula was probably a bit (if not strongly) self-gravitating at the actual orbit of giant planets. Alpha-discs hosted by active galaxies are self-gravitating beyond about a thousand Schwarzchild radii. The inferred surface density remains too high to lower the accretion time scale. More efficient mechanisms driving accretion are required.

We have developed 1D time-dependent numerical models of accretion discs, using an adaptive grid technique and an implicit numerical scheme, in which the disc size is allowed to vary with time. The code fully resolves the cooling and heating fronts propagating in the disc. We show that models in which the radius of the outer edge of the disc is fixed produce incorrect results, from which probably incorrect conclusions about the viscosity law have been inferred. In particular we show that outside-in outbursts are possible when a standard bimodal behaviour of the Shakura-Sunyaev viscosity parameter alpha is used. We also discuss to what extent insufficient grid resolutions have limited the predictive power of previous models. We find that the global properties (magnitudes, etc. ...) of transient discs can be addressed by codes using a high, but reasonable, number of fixed grid points. However, the study of the detailed physical properties of the transition fronts generally requires resolutions which are out of reach of fixed grid codes. It appears that most time-dependent models of accretion discs published in the literature have been limited by resolution effects, improper outer boundary conditions, or both.

We present a model for the dwarf nova WZ Sge which does not require assuming unusually and unexplained low values of the viscosity alpha-parameter during exceptionally long quiescent states of this system. We propose that the inner parts of the accretion disc are disrupted by either a magnetic field or evaporation, so that the disc is stable (or very close to being stable) in quiescence, as the mass transfer rate is very low and the disc can sit on the cool, lower branch of the thermal equilibrium curve. Outbursts are triggered by an enhanced mass transfer, which brings the disc into the unstable regime of the standard dwarf nova disc instability model. The resulting eruptions are strongly affected by the irradiation of the secondary star. Our model reproduces very well the recurrence time and the characteristics of the light curve in outburst.

We show that if the inner regions of accretion discs in quiescent dwarf nova systems are removed by magnetic disruption or evaporated by siphon flows, the remaining disc is globally stable for mass transfer rates $\lta 10^{15}$ g s$^{-1}$. This implies that (super)outbursts in these systems have to be triggered by an enhanced mass transfer form the companion. We suggest that the lack of normal outbursts in WZ Sge results only from its low mass transfer rate and that the viscosity in its disc does not have to be few orders of magnitude lower than in other quiescent dwarf novae.