# A Hydrodynamical Mechanism for Generating Astrophysical Jets

## Abstract

Whenever in a classical accretion disk the thin disk approximation fails interior to a certain radius, a transition from Keplerian to radial infalling trajectories should occur. We show that this transition is actually expected to occur interior to a certain critical radius, provided surface density profiles are steeper than , and further, that it probably corresponds to the observationally inferred phenomena of thick hot walls internally limiting the extent of many stellar accretion disks. Infalling trajectories will lead to the convergent focusing and concentration of matter towards the very central regions, most of which will simply be swallowed by the central object. We show through a perturbative hydrodynamical analysis, that this will naturally develop a well collimated pair of polar jets. A first analytic treatment of the problem described is given, proving the feasibility of purely hydrodynamical mechanisms of astrophysical jet generation. The purely hydrodynamic jet generation mechanism explored here complements existing ideas focusing on the role of magnetic fields, which probably account for the large-scale collimation and stability of jets.

X. Hernandez: Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 70264, CD Universitaria, CP 04510 , México DF, México () Pablo. L. Rendón: Centro de Ciencias Aplicadas y Desarrollo Tecnológico, Universidad Nacional Autónoma de México, AP 70–186, México D.F. 04510, México Rosa G. Rodríguez-Mota: Centro de Ciencias Aplicadas y Desarrollo Tecnológico, Universidad Nacional Autónoma de México, AP 70–186, México D.F. 04510, México A. Capella:Instituto de Matematicas, Universidad Nacional Autónoma de México, México D.F. 04510, México. \listofauthorsHernandez, Rendón, Rodríguez-Mota & Capella \indexauthorHernandez, X. \indexauthorRendón, P.L. \indexauthorRodríguez-Mota, R.G. \indexauthorCapella, A. \resumenWhenever in a classical accretion disk the thin disk approximation fails interior to a certain radius, a transition from Keplerian to radial infalling trajectories should occur. We show that this transition is actually expected to occur interior to a certain critical radius, provided surface density profiles are steeper than , and further, that it probably corresponds to the observationally inferred phenomena of thick hot walls internally limiting the extent of many stellar accretion disks. Infalling trajectories will lead to the convergent focusing and concentration of matter towards the very central regions, most of which will simply be swallowed by the central object. We show through a perturbative hydrodynamical analysis, that this will naturally develop a well collimated pair of polar jets. A first analytic treatment of the problem described is given, proving the feasibility of purely hydrodynamical mechanisms of astrophysical jet generation. The purely hydrodynamic jet generation mechanism explored here complements existing ideas focusing on the role of magnetic fields, which probably account for the large-scale collimation and stability of jets. \addkeywordaccretion, accretion discs \addkeywordhydrodynamics \addkeywordgalaxies:jets \addkeywordISM:jets and outflows

## 1 Introduction

Astrophysical jets occur over a large range of astronomical scales, from the stellar Newtonian scales of HH objects (e.g. Reipurth & Bally 2001), to the relativistic cases of microquasars and gamma ray bursts (e.g. Mirabel & Rodriguez 1999, Mészáros 2002), to the Mega-parsec extragalactic scales of AGN jets (e.g. Marscher et al. 2002). Although the processes of propagation and collimation appear to be relatively well understood in terms of the interplay between the hydrodynamics of the problem (Scheuer 1997) and magnetic fields (e.g. Blandford 1990), the precise mechanism of jet generation remains as the most uncertain part of the system.

Existing jet generation mechanisms focus on the interaction between the rotating material in the inner regions of the accretion disk and the magnetic field of the central (rotating or otherwise) star, or perhaps an external magnetic field threading the disk. See for example Blandford & Payne (1982), Henriksen & Valls-Gabaud (1994), Meier et al. (2001), Price et al. (2003). Although generically successful, detailed comparisons with observations do dot always yield consistent results, e.g. Ferreira et al. (2006). Further, the details of relative orientation of stellar spin, accretion disk and magnetic field configurations, have to be supplied in somewhat highly specific manners.

As an alternative, we explore the possibility of hydrodynamical jet generation mechanism where magnetic fields play no rôle, but only the intrinsic hydrodynamical physics of the interplay between the accretion and the central star wholly determine the characteristics of the jet. The main obstacle to such a scheme is the presence of a centrifugal barrier associated with the angular momentum content of the material in the accretion disk. We show however, that the thin disk approximation for accretion disks, which is equivalent to the assumption of quasi-circular orbits for the disk (e.g. Pringle 1981), is expected to break down internal to a certain critical radius, resulting in a transition to quasi-radial flow for the disk. Then, an analytic first order perturbation treatment serves to demonstrate the feasibility of purely hydrodynamical jets. We show also that many generic features of observed astrophysical jets across a wide range of scales, can be naturally accounted for in the general model presented. The presence of magnetic fields in jet phenomena is evident empirically, however, it is not impossible that their main contribution to the problem could be restricted to the radiation of the jet material and to its long range collimation, with purely hydrodynamical physics playing a part in the actual jet generation mechanisms.

The problem of jet formation has been studied extensively in the context of classical hydrodynamics, most often regarding fluid-body interactions. The appearance of stable coaxial jets resulting from radially-symmetric velocity fields over thin fluid sheets has been established by, among others, Taylor (1960). The rôle played by both the magnitude and the direction of velocity in the formation of this type of jet is the subject of a theoretical study by Glauert (1956), where it is shown that at the point at which the jet forms, a large velocity gradient is observed and momentum flux is constant, with horizontal momentum being transformed into vertical momentum. Similar results were obtained by King & Needham (1994), who provide an asymptotic study of a jet formed by a vertical plate accelerating into a semi-infinite expanse of stationary fluid of finite depth with a free surface and a gravitational restoring force. It is found that as the fluid approaches the plate in a horizontal direction, a gradual rise in free-surface elevation occurs. Eventually, a thin region is reached where vertical velocity dominates horizontal velocity as a consequence of the fluid finding it more difficult to overcome the inertia it would meet by continuing its horizontal motion than escaping vertically towards the low-pressure free surface. In essence, this same mechanism can allow for jet formation in the context of the problem we study here.

In section 2 we develop the criterion for transition to radial flow in a standard accretion disk. The perturbative solution to the resulting problem of a radially infalling disk is then developed in section 3. Section 4 presents trajectories for particular cases of the solution obtained in the previous section, in dimensionless units. Finally, we give our conclusions in section 5.

## 2 The transition to Radial Flow

We start from a standard thin accretion disk where material orbits on quasi-circular Keplerian orbits around a central star of mass M. Assuming axial symmetry and cylindrical coordinates, the total angular momentum of a ring at radius of width will be:

(1) |

where and are the surface density and orbital frequency profiles of the disk, respectively. If the breaking torque on a given ring of the disk, due to its exterior ring is denoted by , the total breaking torque on the ring at radius will be:

(2) |

The rate of change of the angular momentum of the ring in question can now be calculated as:

(3) |

If at any radius a substantial fraction of the angular momentum of the ring is lost through viscous torques over an orbital period, the assumption of quasi-circular orbits will break down, and the disk will make a transition to mostly radial orbits. This condition can be stated as:

(4) |

Substitution of equations (1) and (2) into the above yields:

(5) |

as the condition for the onset of radial flow in the disk.

We can now introduce a model for the torques in terms of the rate of shear and the effective viscosity (e.g. Von Wiezsäcker 1943, Pringle 1981) as:

(6) |

Taking to substitute the above equation into condition (5) gives

(7) |

To proceed further we can take for example, constant and a model for the disk surface density profile of the form:

(8) |

of the type often used in models of accretion disks when fitted to observations (e.g. Hughes et al. 2009). Use of the two above forms for and reduces condition(7) to:

(9) |

The above condition will always be met in any accretion disk with , interior to a transition radius given by:

(10) |

It is interesting that directly observed accretion disks have spectra which when modelled typically yield , in general , e.g. Hartmann et al. (1998), Lada (2006), Hughes et al. (2009). This leads one to expect the transition to radial flow to occur in many of the stellar accretion disks, interior to radii given by equation (10).

Finally, we can explore the consequences of introducing an prescription (Shakura & Sunyaev 1973) into condition (9), , where and are respectively the sound speed and height of the disk, with a dimensionless number, yielding the dimensionless condition:

(11) |

where is the ratio of the Keplerian orbital velocity in the disk to the sound speed in the disk. Although the above condition is only valid for constant , it illustrates the equivalence between the assumption of quasi-circular orbits for the material in the disk, and the assumption that the disk is thin. We see that the breakdown of the assumption of quasi-circular orbits (condition 9), corresponds to the point where disk becomes fat.

Writing in astrophysical units as:

(12) |

we see that for values typical of what is observationally inferred for the stellar mass and position and temperature of accretion disk walls in T Tauri protoplanetary disks, those indicated in the above equation (e.g. D’Alessio et al. 2005, Espaillat et al. 2007, Hughes et al 2009), one approaches the breakdown of the thin disk approximation, , and consequently the transition to radial flow.

Taking typical inferred values for at the wall of (e.g. D’Alessio et al. 2005, Espaillat et al. 2007, Hughes et al 2009), we can now write the condition for the transition to radial flow, locally at the wall, as:

(13) |

We see that for a standard value of the above equation implies a reasonable value of at the thick wall, significantly higher than the values of and lower, which apply for the body of the thin accretion disk beyond this radius. A substantial increase of as decreases is expected in any turbulence driven viscosity model for accretion disks, e.g. Firmani, Hernandez & Gallagher (1996), in the context of galactic disks.

In terms of the debate surrounding the inference of inner holes in observed accretion disks, many solutions have been proposed in terms of disk clearing mechanisms; grain growth (e.g. Strom et al. 1989, Dullemond & Dominik 2005), photoevaporation (e.g. Clarke et al. 2001), magnetorotational instability inside-out clearing (e.g. Chiang & Murray-Clay 2007, Dutrey et al. 2008), binarity (e.g. Ireland & Kraus 2008) and planet-disk interactions (e.g. Rice et al. 2003). None of the above is entirely satisfactory, as noted by Hughes et al. (2009), mostly due to their incompatibility with a steady state solution. An alternative solution under the proposed scenario, is that there is no actual disk material clearing, only a transition to predominantly radial flow at the thick wall, and consequently a shear-free flow interior to this point. Once the disk heating mechanism is removed, as one should expect from the analysis presented in this section, the inner disk disappears from sight.

This transition at a thickened inner boundary is qualitatively what ADAF models propose. In those models, substantial turbulent and thermal pressure remain within the inner flow. Still, given the present lack of a definitive jet generation mechanism, we consider it interesting to explore the consequences of a model where cooling drives the flow to a condition. As we show in the following sections, we prove that such a situation will quite naturally yield purely hydrodynamical jet solutions, which we find encouraging.

## 3 Hydrodynamical Jet Solutions

We shall now model the physical situation resulting from the scenario described above as an axially symmetrical distribution of gas in free fall towards a central star of mass . Taking a spherical coordinate system with the angle between the positive vertical direction and the position vector we have:

(14) |

(15) |

(16) |

for the continuity equation, and the radial and angular components of Euler’s equation. In the above, , , and give the radial velocity, angular velocity, matter density and pressure, respectively. We have neglected temporal derivatives, as we are interested at this point, in the characteristics of steady state solutions, although the temporal derivatives are reinstated in the Appendix for the purpose of stability analysis. We take as a background state a free-falling axially symmetrical distribution of gas described by , and all components of negligible, a consistent solution to eqs.(15) and (16), used only in the description of this background. Having ignored the inclusion of the radius at which the transition to radial flow takes place in the choice of limits the validity of the analysis to radial scales along the plane of the disk much smaller than . This is justified by the fact that jets appear as phenomena extremely localised towards . We now take a density profile given by:

(17) |

were is a dimensionless function of and describes the polar angle dependence of the infalling material, for example, one can ask for , diminishing symmetrically towards the poles. The choice of this last function will determine the details of the problem, and can be thought of as something of the type

(18) |

with a normalisation constant and a form constant describing the flattening of the disk of infalling material. However, we shall mostly leave results indicated in terms of . The continuity equation (14) now fixes through:

(19) |

and hence , which completes the description of the background state through:

(20) |

with a constant which determines the point at which .

We now study the first order departures from the assumed background state through a perturbative approach, to explore the possibility that such departures might naturally give rise to a hydrodynamical jet. Regarding these first order perturbations, we shall also be interested in steady state solutions given by:

(21) |

(22) |

(23) |

where quantities with subscript (1) denote the perturbation on the background solution. Notice that is, as required by the assumed background state, a function only of the radial coordinate . For the background state assumed, the perturbations on it will be solved for in a fully self consistent manner. Writing eqs.(14), (15) and (16) to first order in the perturbation one obtains after rearranging terms:

(24) |

(25) |

(26) |

In the above we have assumed an isothermal equation of state for the perturbation , an assumption often used in the modeling of astrophysical jets, e.g. the T Tauri jets observed and modelled by Hartigan et al. (2004). This idealised case serves to illustrate clearly the consequences of the physical setup being considered, as it allows for an analytic solution. The generalisation to more realistic adiabatic, polytropic or otherwise equations of state can be performed numerically, and can be expected to yield qualitatively similar results, although interesting differences in the details can be expected to emerge, which will be considered latter. In the above three equations we have introduced the constants and .

To make further progress we can attempt a solution through the method of separation of variables, proposing a solution of the form: , , . This ansatz yields two independent systems of three equations each, one for the radial, and one for the angular dependences of the perturbations. The radial equations become:

(27) |

(28) |

(29) |

Where we have used the results of the angular ones:

(30) |

(31) |

(32) |

In splitting the radial and angular dependences of the perturbed continuity equation, eq.(24), we have used the result of the angular equation of the perturbed radial Euler equation, eq.(30). The constants , and are the separation constants of the problem. Firstly we turn to the angular system, which allows an exact solution.

We can take eq.(32) and substitute into it the product from eq.(31), and the product from eq.(30) to obtain an equation involving only :

(33) |

having solution:

(34) |

In the above equation and are normalisation constants, and and are the Legendre polynomials of the first and second kinds, respectively. The index of these functions is given by the relation , where . As typical of separation of variables problems, we see the solution as an eigenvalue problem, with the separation constants determining the order of the solution function. The requirement of axial symmetry forces to be even. Any such desired angular distribution for can now be constructed as an infinite series of the above functions. For simplicity we analyse the case of , , :

(35) |

This will result in slightly less material along the plane of the disk, and slightly more along the poles, for the total density , with respect to the background state . At this point eqs.(30) and (31) can be used to obtain:

(36) |

(37) |

We see the potential for jet formation already in the two above equations; if the background state is relatively devoid of material towards the poles, could be small for , and eq.(36) will lead to large radial velocities along the poles. Further, eq.(37) shows that the angular velocity will always tend to zero both on the plane of the disk, and along the poles, where movement will necessarily be radial, although not necessarily positive.

The case of the radial system is more cumbersome, as the linear operator which appears is of third order. However, in the Appendix we show through a first order stability analysis which terms can be safely ignored in the high-frequency approximation, and for small , the region where jets develop. Equation (58) is used to show that the right-hand side of eq.(29) can be ignored to leading order, so that the resulting equation can be readily integrated to give:

(38) |

In the above we have taken the integration constant which appears as being equal to zero, from requiring for . Substituting the above relation for into eq.(27) leads to:

(39) |

We see from equation (59) in the Appendix that the second term on the left-hand side of the above equation can be dismissed, to leading order, in the regime of validity there established. We then obtain:

(40) |

(41) |

where is a characteristic radius at which . Now from equation (28), where all the terms are of the same order, as is shown in equation (62), we have

(42) |

In the above equation we have also taken the integration constant as zero, from requiring for . Choosing without loss of generality the two characteristic radii and both equal to , we can now write the full solution to the perturbation as:

(43) |

(44) |

(45) |

where we have introduced and as the angular part of , . The full solution to the perturbation can be seen to depend only on the two parameters and , normalisation constants for the densities of the background state and the perturbation, with the velocities of the perturbation solution depending only, and linearly, on the ratio , which determines the validity regime of the perturbative approach through the condition .

As was already evident from eq.(37), we see form eq.(45) that the angular velocity will be zero only for , and . Thus, movement along the plane of the disk will remain along the plane, but also, along the poles movement will be exclusively radial. This last point, together with the positive sign of the radial velocity along the poles (c.f. eq.(44) for ), provides for a well collimated jet along the poles. Hence, the background state proposed is seen to be unstable towards jet formation.

From eq.(44) we see that one has only to ask for a background state where matter density does not grow towards the poles, in order to obtain ejection velocities along them, which could become very large for flattened disks with relatively empty poles. The axial symmetry condition imposed on eq.(34) guarantees both axial symmetry and symmetry above and below the plane of the disk for the full solution. We see also that if one takes higher orders for , the index of the Legendre polynomial solution to eq.(33), one obtains increasingly more critical angles at positions intermediary between and at which the angular velocity goes to zero. In fact, more complex geometries and asymmetric jets appear, as inferred observationally by e.g. Ferreira et al. (2006), if is taken as an arbitrary real number. However, modelling a situation where a polar jet dominates the ejection identifies as the leading order.

Notice that the qualitative behaviour of the solution is guaranteed by the exactness of the angular solution, the approximation used for solving the radial problem will only introduce an error in the magnitudes of the velocities outside of , but will not change the fact that velocities will be of radial infall along the plane of the disk, , and of radial outflow along the poles, the jet solution for . This is reinforced by the stability analysis shown in the appendix, where we show that the angular system of equations does not respond to a temporal perturbation to the first order solution treated here, while the response of the radial system is merely the introduction of periodic and bound factors multiplying the radial solutions given in this section.

The two constants of the problem, and , can now be calculated once a choice of is specified, from the two conditions:

(46) |

(47) |

where is a suitable angle defining the opening of the jet, in all likelihood very small, as will be see in the following section. In the above equations gives the matter accretion rate onto the central star, and the matter ejection rate due to the jet. Dimensionally, the two quantities above will scale as:

(48) |

(49) |

where and are two dimensionless constants which will depend on the choice of , and which would be expected to be of order 1.

Qualitatively, this type of model naturally furnishes a tight disk-jet connection (c.f. eq. 44) e.g., as now firmly established in microquasars and AGN jets (see e.g. Marscher et al. 2002, Chatterjee et al. 2009). In the above systems bursts of enhanced jet activity are seen to follow temporal dips in disk luminosity output after small characteristic delay times. In the present model, such a situation would be expected if the critical radius for transition to radial flow in the disk made a sudden transition to higher values. Again, the drop in disk output might not reflect the disk material disappearing (in this case being swallowed by the central black hole, as sometimes proposed), but simply fading from view as heating mechanisms shut down, then naturally enhancing jet activity as the effective increases.

## 4 Particular Solutions

In order to present a sample of the trajectories expected in the model, we turn to the full solution to the problem for , eqs.(21) and (22), but written in dimensionless form:

(50) |

(51) |

The above remain in spherical coordinates, where and . Notice that the threshold for jet development can now be stated as the condition , . A choice of and now allows to numerically integrate trajectories. We take:

(52) |

as a model for the background state, where small values of will result in thin disks, while more spherical mass distributions correspond to large values of . As a first example we take () and , a fairly thick disk, with only very small changes in density from the plane to the poles of a factor of 1.6, an almost spherically symmetric density configuration compatible with what appears in some ADAF models. This value of is far from representing a flattened disk, and hence one which does not force the jet solution, see eq.(44). With these parameters, and initial conditions specified in dimensionless cylindrical coordinates as and ranging from 1.0 to 1.4, we solve eqs.(50) and (51) through a finite differences scheme to plot figure 1.

In spite of having taken , the validity of the perturbative approach could be broken by the appearance of large velocities in the perturbation, breaking the validity condition . We explicitly keep track of this condition, and show the appearance of perturbation radial velocities violating it with a change to thin lines in the figure, angular velocities for the perturbation will always remain of order , provided .

For this case, the two lowermost curves present trajectories which all turn downwards to converge onto the central star. These are solutions which essentially follow the background state, infalling onto the bottom of the potential well. As one raises the initial value of however, the following 4 curves begin to increasingly deviate from the free fall trajectories and, although still ending up at , clearly show a change in behaviour. A threshold is eventually crossed and curves of a very different type ensue, the six upper jet trajectories shown in figure 1. We see the pressure gradients associated with the distribution of matter in the background solution acting to break the fall of the incoming material, turn it back, and then accelerate it vertically through the vertical density gradients. These jet trajectories rapidly converge with height to eventually yield a well collimated structure.

Note that the region where the perturbation condition ceases to be valid appears only on the ”jet” trajectories, and only after these have turned upwards, at points where in fact, velocities already exceed the local escape velocities. This point guarantees the appearance of a jet, even though the details of the flow in this jet region could differ somewhat from what is shown by the thin lines in the figure. The same applies to the validity region of the radial equation: the jet solution appears within the validity regime of all the approximations taken. the qualitative form of the full solution will not deviate much from what is shown in figure 1, due to the exact character of the angular solution.

Notice that the generic validity of the solution proposed does not require that all of the disk material should loose all of its angular momentum at the critical radius , only that some of the disk material should loose most of its angular momentum at that point. Any minor, residual angular momentum remaining on the disk fraction which forms the jet will only establish a finite jet cross-sectional area through the appearance of a centrifugal barrier, which will present only a small correction on the scenario given here. Still, the empirical presence of spectroscopically studied accretion disks truncated interior to certain critical radii suggests the very substantial reduction of the shears in that region, as could happen if a transition to mostly radial flow takes place.

In going back to eqs.(50) and (51) we see that one expects

(53) |

It is reassuring that for the jets associated with T Tauri stars for example, values of of between 0.1 and 0.01 on average are observationally inferred and therefore of order (e.g Hartigan et al. 1996, Gullbring et al. 1998, Hartigan et al. 2004, Ferreira et al. 2006). These are compatible with the value used to plot fig. (1), and hence ones which will readily yield jet solutions.

Figure 2 is analogous to fig. 1, but presents an example for a much more flattened disk having , a significantly flattened disk where the density contrast between plane and pole at fixed radius is now of a factor of 2.5, orders of magnitude larger than in the last example. We take this time . Again, we see the main features of the solution being well established in the region interior to , with a jet solution again established within the validity regime of all approximations used.

The long range stability and coherence of these structures lies outside the scope of this work, and is in all probability furnished by a series of mechanisms extensively explored in the literature including angular momentum, magnetic fields and pressure containment of the surrounding medium, e.g. Begelman et al. (1984), Blandford (1990), Falle (1991), Kaiser & Alexander (1997).

In going to the more extreme jet phenomena associated with stellar black holes (e.g Mirabel & Rodriguez 1999), quasars (e.g. Marscher 2002) and gamma ray bursts (e.g. Mészáros 2002), it is natural to expect the ideas presented here to apply, but amplified to a much more extreme regime by the appearance of corresponding relativistic and general relativistic effects, to first order, the shift in the divergence in the potential from in the Newtonian case to in the general relativistic one. It is therefore natural to expect purely hydrodynamic jet generation mechanism to apply across all classes of astrophysical objects, specially given the qualitative scalings and similarities which appear over all astrophysical jet classes (e.g. Miarabel & Rodriguez 1999, Mendoza, Hernandez & Lee 2005), in addition to the magnetically driven processes traditionally found in the literature.

Notice from eq.(44) the intimate link between the jet velocity and the physical state of the infalling material, and through . This implies that temporal variations in the parameters of the infalling material will result in temporal variations in the density and velocity of the jet material, in a way described by eqs.(43) and (44). The above can serve as a physical description of the key processes relevant to the formation of internal shocks in astrophysical jets, the main ingredient behind phenomena such as HH objects and gamma ray bursts.

Although the solution to the problem presented is rigorous, and the appearance of the resulting hydrodynamical jets solid, the relevance of our model is clearly constrained by the validity, or otherwise, of the starting hypothesis; most critically, the appearance of an inner region in free fall. Although plausibility arguments for the appearance of this regime have been given, the reality of it is certainly not assured. Still, given the absence of a fully satisfactory and general solution to the problem of jet generation, and given the clear and transparent physics of the one shown here, we believe it is interesting to present such a simple hydrodynamical jet solution.

## 5 Conclusions

We show that given a radially infalling accretion, a purely hydrodynamical jet ensues.

We calculate the condition for the transition from quasi-circular to quasi-radial flow in a standard accretion disk, and show it will always occur for power law surface density profiles of the form , interior to a critical radius, provided .

Comparison with inferred inner holes in observed accretion disks yields results consistent with our estimates for the above transition radius, the point where shears in the flow would be substantially reduced.

Well collimated jets readily appear, proving the existence of purely hydrodynamical mechanisms for the generation of astrophysical jets.

## 6 Acknowledgments

Xavier Hernández acknowledges the hospitality of the Observatoire de Paris for the duration of a sabbatical stay during which many of the ideas presented here were first developed, and financial support from project PAPIIT IN103011 DGAPA UNAM. Pablo L. Rendón and Rosa Rodríguez acknowledge financial support from project PAPIIT IN110411 DGAPA UNAM. Antonio Capella acknowledges financial support from project PAPIIT IN101410 DGAPA UNAM.

## Appendix A Appendix: Stability Analysis

We study the stability of the system by means of a standard linear perturbation approach. Furthermore, we develop here a first order analysis that also serves to justify the approximations made in section 3 of the main body of the paper. We follow the approach used by Papaloizou & Pringle (1984) to study the stability of this system by means of linear perturbation of the time-dependent equations of the fluid, as opposed to the steady-state versions stated previously as equations (24)-(26). The time-dependence of the perturbed quantities is written as , where is a temporal frequency. We then obtain

(54) | |||||

for the continuity equation, and

(55) |

(56) |

for the radial and angular components of the Euler equation, respectively. We now proceed as before, separating variables in such a way that we require the angular components , and to be dimensionless, while the radial components , and have the same dimensions as the original functions. Encouragingly, the equations for the angular components are exactly the same as those obtained for the steady-state analysis, equations (30)-(32). We conclude that the angular components are not time-dependent and are thus unconditionally stable, reinforcing the overall character of the solution presented in section 3. The resulting equations for the radial components do include an extra term due to the temporal dependence. In order to establish the relative magnitude of the different terms in these equations, we shall rewrite them in dimensionless form, using the following scalings for that effect:

(57) |

where is a characteristic density, , and . We now consider the high-frequency limit, as Papaloizou & Pringle (1984) have also done. Notice then that . Equation (54) may now be rewritten as

(58) |

where . If we now require and , we observe that all the terms on the left-hand side of the equation are , whereas the term on the right-hand side is . Thus, at leading order the term on the right-hand side vanishes.

Using this same scaling, equation (55) is now

(59) |

where we observe that the term on the right-hand side is . Again, at , the right-hand side of the equation vanishes, and we may solve the equation to obtain

(60) |

(61) | |||||

where is a constant, and is the incomplete gamma function. The real parts of these solutions are the dimensionless forms of the solutions we obtained for and in section 3 of this paper, modulated by an oscillating and bounded function.

Finally, in dimensionless form, equation (56) is

(62) |

where . If we now require , all three terms are . Notice that the orders of magnitude chosen here for the separation constants give as a result , so that the choice of value given to in section 3 is consistent with our analysis. Taking the real part of the solution given in (61) and substituting in the above equation, we may solve to obtain

(63) | |||||

where is a constant. Considering the real part of the above solution, we again find that we have the dimensionless form of the solution obtained for earlier in this paper, modulated by an oscillating, bounded function.

Thus, the effect of linear perturbations on the time-dependent system is to provoke bounded oscillations around the steady-state solutions for density and both angular and radial velocities. We conclude that the system is neutrally stable.

Our choice of scalings has permitted us to write the equations of the fluid in dimensionless form in terms of only three dimensionless parameters, , and . In this form, as is expected, it becomes rather more straightforward to establish the leading-order terms in the equations. The hypothesis of small radial distance translates into the smallness of , and and , formally very similar, are related to the separation constants. It should be pointed out that the dimensional constants and are independently required to be of the same order, and by simply letting the dimensionless constant be , the sole restriction placed upon these three constants in section 3, that also be , is fulfilled. We can then consistently fix the order of the Legendre polynomial solutions presented in section 3 as a function of the particular values given to the separation constants, a behaviour typically associated with the method of separation of variables.

### Footnotes

- affiliation: Instituto de Astronomía, UNAM, México.
- affiliationmark:
- affiliationmark:
- affiliationmark:

### References

- Begelman, M. Blandford, R. D., Rees, M., 1984, Rev. Mod. Phys., 56, 256
- Blandford, R. D., Payne, D. G., 1982, MNRAS, 199, 883
- Blandford, R. D., 1990, in “Active Galactic Nuclei, eds. Courvoisier, T. L., Mayor, M., Saas-Fee Advanced Course 20 (Les Diablerets: Springer-Verlag), 161-275
- Chatterjee, R., et al., 2009, ApJ, 704, 1689
- Chiang, E., Murray-Clay, R., 2007, Nat. Phys., 3, 604
- Clarke, C. J., Gendrin, A., Sotomayor, M., 2001, MNRAS, 328, 485
- D’Alessio, P., et al., 2005, ApJ, 621, 461
- Dullemond, C. P., Dominik, C., 2005, A&A, 434, 971
- Dutrey, A., et al., 2008, A&A, 490, L15
- Falle, S. A. E. G., 1991, MNRAS, 250, 581
- Ferreira, J., Dougados, C., Cabrit, S., 2006, A&A, 453, 785
- Firmani, C., Hernandez, X., Gallagher, J., 1996, A&A, 308 403
- Glauert, M. B., 1956, J. Fluid Mech., 1, 625
- Gullbring, E., Hartmann, L., Briceño, C., Calvet, N., 1998, ApJ, 492, 323
- Hartmann, L., Calvet, N., Gullbring, E., D’Alessio, P., 1998, ApJ, 495, 385
- Hartigan, P., Edwards, S., Gandhour, L., 1995, ApJ, 452, 736
- Hartigan, P., Edwards, S., Pierson, R., 2004, ApJ, 609, 261
- Henriksen, R. N., Valls-Gabaud, D., 1994, MNRAS, 266, 681
- Hughes, A. M., Wilner, D. J., Calvet, N., D’Alessio, P., Claussen, M. J., Hogerheijde, M. R., ApJ, 2007, 664, 536
- Hughes, A. M., et al., 2009, ApJ, 698, 131
- Ireland, M. J., Kraus, A. L., 2008, ApJ, 678, L59
- Keiser, C. R., Alexander, P., 1997, MNRAS, 286, 215
- King, A. C., Needham, D. J., 1994, J. Fluid Mech., 268, 89
- Lada, C., et al., 2006, AJ, 131, 1574
- Marscher, A. P., Jorstad, S. G., Gomez, J., Aller, M. F., Teräsranta, H., Lister, M. L., Stirling, A. M., 2002, Nature, 417, 625
- Meier, D. L., Koide, S., Uchida, Y., 2001, Science, 291, 84
- Mendoza, S., Hernandez, X., Lee, W. H., 2005, Rev. Mex. Astron. Astrofis., 41, 61
- Mészáros, P., 2002, ARA&A, 40, 137
- Mirabel, I. F., Rodriguez, L. F., 1999, ARA&A, 37, 409
- Papaloizou, J. C. B., Pringle, J. E., 1984, MNRAS, 208, 721
- Price, D. J., Pringle, J. E., King, A. R., 2003, MNRAS, 339, 1223
- Pringle, J. E., 1981, ARA&A, 19, 137
- Reipurth, B., Bally, J., 2001, ARA&A, 39, 403
- Rice, W. K. M., Wood, K., Armitage, P. J., Whitney, B. A., Bjorkman, J. E., 2003, MNRAS, 324, 79
- Shakura, N. I., Sunyaev, R. A., 1973, A&A, 24, 337
- Scheuer, P. A. G., 1974, MNRAS, 166, 513
- Strom, K. M.,Strom, S. E., Edwards, S., Cabrit, S., Skrutskie, M. F., 1989, AJ, 97, 1451
- Taylor, G., 1960, Proc. R. Soc. Lond. A, 259, 1
- von Weizsäcker, C. F., 1943, Z. Astrophys., 22, 319