Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Constraining the number of compact remnants near Sgr A*

Due to dynamical friction stellar mass black holes and neutron stars are expected to form high-density cusps in the inner parsec of our Galaxy. These compact remnants, expected to number around 20 000, may be accreting cold dense gas present there, and give rise to potentially observable X-ray emission. Here we build a simple but detailed time-dependent model of such emission. The possibility that these accretion flows are radiatively inefficient is taken into account and brings in some uncertainty in the conclusions. Despite this uncertainty, we find that at least several X-ray sources of this nature should be detectable with Chandra at any one time. Turning this issue around, we also ask a question of what current observational constraints might be telling us about the total number of compact remnants. In our 'best guess' model, a cusp of ∼40 000 remnants overpredicts the number of discrete sources and the total X-ray luminosity of the inner parsec, and is hence ruled out. In the most radiatively inefficient scenario that we consider, the radiative efficiency is set to be as small as  ɛ= 10⁻⁵  . In this rather unlikely scenario, a cusp of ∼40 000 black holes would be allowed by the data, but several individual sources should still be visible. Future observations of the distribution and orbits of the cold ionized gas in the inner parsec of our Galaxy will put tighter constraints on the cusp of compact remnants.     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

Accretion in stellar clusters and the initial mass function

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses M ∗≳0.5 M_⊙ . The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a ρ  ∝  R ⁻² distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of γ =−3/2 (where Salpeter is γ =−2.35) . Once the stars dominate the potential and are virialized, which occurs first in the cluster core, the accretion radius is the Bondi–Hoyle radius. The resultant mass spectrum has an asymptotic limit of γ =−2 with slightly steeper slopes ( γ ≈−2.5) if the stars are already mass-segregated. Simulations of accretion on to clusters containing 1000 stars show that, as expected, the low-mass stars accumulate the majority of their masses during the gas-dominated phase whereas the high-mass stars accumulate the majority of their masses during the stellar-dominated phase. This results in a mass spectrum with a relatively shallow γ ≈3/2 power law for low-mass stars and a steeper power law for high-mass stars −2.5≲ γ ≤−2 . This competitive accretion model also results in a mass-segregated cluster.     

The impact of a supernova explosion in a very massive binary

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than  100 M_⊙  . In a Population III binary system, a more massive star can result in the formation of a black hole (BH) and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early Universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of  10³ M_⊙  that results from a direct collapse of a very massive star and a companion star of  100 M_⊙  that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early Universe.     

Strong helium 10 830-Å absorption in Sakurai's object (V4334 Sgr)

We report the appearance and evolution during 1998 of strong neutral helium ³S–³P^o absorption at ∼10 830 Å in Sakurai's Object (V4334 Sgr), which is believed to be a planetary nebula nucleus (PNN) undergoing a final helium shell-flash. First detected on 1998 March 18, the profile of the He  i feature is P Cygni-like. The absorption depth has increased in three subsequent spectra in 1998. If this is owing to a wind, the profile indicates a wind velocity of ∼670±50 km s⁻¹. The strong C  i 10 690-Å line seen prior to the appearance of the helium feature has disappeared; however Sr  ii and CN absorption features remain present. We tentatively identify several new features as Si  i . Taken together with other observations we suggest that the data are consistent with Sakurai's Object entering a phase in which it seems to have become a member of the R Coronae Borealis-type class of stars.     

The eccentricities of the barium stars

We investigate the eccentricities of barium (Ba  ii ) stars formed via a stellar wind accretion model. We carry out a series of Monte Carlo simulations using a rapid binary evolution algorithm, which incorporates full tidal evolution, mass loss and accretion, and nucleosynthesis and dredge-up on the thermally pulsing asymptotic giant branch. We follow the enhancement of barium in the envelope of the accreting main-sequence companion and dilution into its convective envelope once the star ascends the giant branch.
The observed eccentricities of Ba  ii stars are significantly smaller than those of an equivalent set of normal red giants but are nevertheless non-zero. We show that such a distribution of eccentricities is consistent with a wind accretion model for Ba  ii star production with weak viscous tidal dissipation in the convective envelopes of giant stars. We successfully model the distribution of orbital periods and the number of observed Ba  ii stars. The actual distribution of eccentricities is quite sensitive to the strength of the tides, so that we are able to confirm that this strength is close to, but less than, what is expected theoretically and found with alternative observational tests. Two systems – one very short-period but eccentric, and one long-period and highly eccentric – still lie outside the envelope of our models, and so require a more exotic formation mechanism. All our models, even those which were a good fit to the observed distributions, overproduced the number of high-period barium stars, a problem that could not be solved by some combination of the three parameters: tidal strength, tidal enhancement and wind accretion efficiency.     

Inflow and outflow from the accretion disc of the microquasar SS 433: UKIRT spectroscopy

A succession of near-infrared (near-IR) spectroscopic observations, taken nightly throughout an entire cycle of SS 433's orbit, reveal (i) the persistent signature of SS 433's accretion disc, having a rotation speed of  ∼500 km s⁻¹  , (ii) the presence of circumbinary disc recently discovered at optical wavelengths by Blundell, Bowler & Schmidtobreick (2008) and (iii) a much faster outflow than has previously been measured for the disc wind, with a terminal velocity of  ∼1500 km s⁻¹  . The increased wind terminal velocity results in a mass-loss rate of  ∼10⁻⁴ M_⊙ yr⁻¹  . These, together with the newly (upwardly) determined masses for the components of the SS 433 system, result in an accurate diagnosis of the extent to which SS 433 has super-Eddington flows. Our observations imply that the size of the companion star is comparable with the semiminor axis of the orbit which is given by     , where e is the eccentricity. Our relatively spectral resolution at these near-IR wavelengths has enabled us to deconstruct the different components that comprise the Brackett-γ (Brγ) line in this binary system, and their physical origins. With this line being dominated throughout our series of observations by the disc wind, and the accretion disc itself being only a minority (∼15 per cent) contribution, we caution against use of the unresolved Brγ line intensity as an 'accretion signature' in X-ray binaries or microquasars in any quantitative way.     

The Galactic Center: a laboratory for AGN?

Summary. This paper reviews the physical state of stars and Interstellar Matter in the Galactic Bulge (radius kpc from the dynamical center of the Galaxy), in the Nuclear Bulge (kpc) and in the Sgr A Radio and GMC Complex, i.e. the central \,pc of our Galaxy. The Galactic Bulge is devoid of cold Interstellar Matter and consists mainly of old stars, while the Nuclear Bulge accounts for of the mass of all of the Interstellar Matter in the Galaxy. A similar ratio holds for the formation rate of medium and high mass stars in Bulge and Disk. The metal abundance of the Interstellar Matter in the Galactic Bulge is found to be . The H-to-CO conversion factors to be applied to molecular gas in the Central Region are by factors 3 (Arimoto et al. 1996) to 10 (Sodroski et al. 1995) lower than in the solar vicinity. Hence, most H masses derived for the Central Region appear to be considerably overestimated. The Nuclear Bulge is pervaded by a thermal plasma (K) which is responsible for the diffuse free-free emission. Lyman continuum photon and dust IR luminosity of the Nuclear Bulge again account for of the respective total luminosities of the Galaxy. Magnetic fields in the Nuclear Bulge are strong (up to mG) as compared with the Galactic Disk (a few tens of G). The field lines are oriented parallel to the galactic plane inside giant molecular clouds and perpendicular to the plane in the intercloud medium. The compact source Sgr A* is close to or at the dynamical center of the Galaxy. Its radio spectrum with a high frequency cut-off at GHz, a low frequency turnover at GHz and a flux density dependence in between can be explained by synchrotron emission from quasi-monoenergetic relativistic electrons. Due to an extinction between Sun and Galactic Center corresponding to , an intrinsic weakness of this source in the near infrared, and a strong background emission from warm dust there are only upper limits available for the flux density of Sgr A* in the far, mid and near infrared and X-ray regime. The size of Sgr A* in the radio regime is cm, its dereddened K-band flux density is mJy, its luminosity has upper limits of (if radiation comes from an Accretion Disk) and (if black-body radiation from an object with a single temperature of K is assumed). If anyone of the soft X-ray sources detected by ROSAT actually coincides with Sgr A*, its X-ray luminosity would be less than a few . With a dark mass of Sgr A* is the best candidate for a starving black hole, although there are no observational indications for the presence of a (Standard) Accretion Disk. While the radio/IR spectrum of Sgr A* is purely nonthermal, the spectrum integrated over the central parsec resembles that of a Seyfert galaxy. Sgr A* is embedded in the Hii region Sgr A West with part of the ionized gas forming a minispiral. Sgr A West is surrounded by the Circum Nuclear Disk, an irregular shaped assembly of molecular gas which extends from pc and rotates around the Galactic Center with an estimated dynamical time scale of \,yr. The total luminosity of of the central parsec is due to the radiation of early-type stars of which have now been directly identified as luminous blue supergiants. It is still debated, however, if these stars can also account for all of the ionization of Sgr A West. In addition, the central parsec contains red giants, AGB stars, and a few super giants of which the brightest are now identified by direct imaging. These stars – together with a few million low mass main sequence stars – account for the bulk of the 2.2\,m emission. The spatial distributions of the three stellar populations in the central pc are remarkably different. Sgr A* is – along the line-of-sight – presumably located close to the center of the Hii region Sgr A West, which in turn is located in front of the extended (pc) synchrotron source Sgr A East, which appears to be the remnant of a gigantic explosion (of the order of the energy of a single supernova explosion) which took place yr ago inside the GMC Sgr A East Core. X-ray observations show within pc a pervasive hot (keV) plasma of expansion age of yr. Both phenomena – as well as the formation of the Circum Nuclear Disk – may have the same origin. Influx of material is observed within the Nuclear Bulge on all distance scales. In the Nuclear Bulge (pc) as well as in the Circum Nuclear Disk (pc) inflow towards the Galactic Center occurs primarily in the galactic plane and amounts to a few . The accretion rate into the central Black Hole, deduced from the luminosity of Sgr A*, however, appears to be lower by at least five orders of magnitude (assuming standard disk accretion). But in an equilibrium state only part of the infalling mass which is not accreted by the Black Hole can be consumed by star formation. A mass inflow rate varying with time is a more natural explanation. Comparing the physical state of the Center of our Galaxy with that of Active Galactic Nuclei derived from observations and modelling, we find that most of the basic characteristics of an AGN are also present in the Galactic Center. Lacking are, however, both the evidence for a standard Accretion Disk and a hard UV spectrum with accompanying high excitation emission lines in the Galactic Center which are characteristic for AGN. The luminosity of the central parsec, , amounts to only of the total luminosity of the Galaxy of . Seen from a distance of M31 (kpc) with an angular resolution of (corresponding to a linear size of pc) the Center of our Galaxy would appear as a mildly active nucleus with some starburst activity and would probably be classified as a weak Seyfert galaxy. The synchrotron spectrum of Sgr A*, however, would be completely masked by reprocessed stellar light (i.e. free-free and dust emission). Received: October 21, 1996     

Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.     

X-rays from cusps of compact remnants near galactic centres

Compact remnants – stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low-mass stars, forming a high concentration cusp. Same physical region may also contain very high-density molecular clouds and accretion discs that are needed to fuel supermassive black hole (SMBH) activity. Here we estimate gas capture rates on to the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This 'diffuse' emission may be of importance for local moderately bright active galactic nuclei (AGNs), especially low-luminosity AGNs. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.     

SPH simulations of negative (nodal) superhumps: a parametric study

Negative superhumps in cataclysmic variable systems result when the accretion disc is tilted with respect to the orbital plane. The line of nodes of the tilted disc precesses slowly in the retrograde direction, resulting in a photometric signal with a period slightly less than the orbital period. We use the method of smoothed particle hydrodynamics to simulate a series of models of differing mass ratio and effective viscosity to determine the retrograde precession period and superhump period deficit  ɛ₋  as a function of system mass ratio q . We tabulate our results and present fits to both  ɛ₋  and  ɛ₊  versus q , as well as compare the numerical results with those compiled from the literature of negative superhump observations. One surprising result is that while we find negative superhumps most clearly in simulations with an accretion stream present, we also find evidence for negative superhumps in simulations in which we shut off the mass transfer stream completely, indicating that the origin of the photometric signal is more complicated than previously believed.     

A possible EGRET source with a two-sided radio jet structure on a parsec scale – 3EG J0808+5114

Recent observations have revealed that damped Lyα clouds (DLAs) host star formation activity. In order to examine if such star formation activity can be triggered by ionization fronts, we perform high-resolution hydrodynamics and radiative transfer simulations of the effect of radiative feedback from propagating ionization fronts on high-density clumps. We examine two sources of ultraviolet (UV) radiation field to which high-redshift ( z ∼ 3) galaxies could be exposed: one corresponding to the UV radiation originating from stars within the DLA, itself, and the other corresponding to the UV background radiation. We find that, for larger clouds, the propagating I-fronts created by local stellar sources can trigger cooling instability and collapse of significant part, up to 85 per cent, of the cloud, creating conditions for star formation in a time-scale of a few Myr. The passage of the I-front also triggers collapse of smaller clumps (with radii below ∼4 pc), but in these cases the resulting cold and dense gas does not reach conditions conducive to star formation. Assuming that 85 per cent of the gas initially in the clump is converted into stars, we obtain a star formation rate of  ∼0.25 M_⊙ yr⁻¹ kpc⁻²  . This is somewhat higher than the value derived from recent observations. On the other hand, the background UV radiation which has harder spectrum fails to trigger cooling and collapse. Instead, the hard photons which have long mean free-path heat the dense clumps, which as a result expand and essentially dissolve in the ambient medium. Therefore, the star formation activity in DLAs is strongly regulated by the radiative feedback, both from the external UV background and internal stellar sources and we predict quiescent evolution of DLAs (not starburst-like evolution).     

Constraints on the ionizing flux emitted by T Tauri stars

We present the results of an analysis of ultraviolet observations of T Tauri stars (TTs). By analysing emission measures taken from the literature, we derive rates of ionizing photons from the chromospheres of five classical TTs in the range  ∼10⁴¹–10⁴⁴  photon s⁻¹, although these values are subject to large uncertainties. We propose that the He  ii /C  iv line ratio can be used as a reddening-independent indicator of the hardness of the ultraviolet spectrum emitted by TTs. By studying this line ratio in a much larger sample of objects, we find evidence for an ionizing flux which does not decrease, and may even increase, as TTs evolve. This implies that a significant fraction of the ionizing flux from TTs is not powered by the accretion of disc material on to the central object, and we discuss the significance of this result and its implications for models of disc evolution. The presence of a significant ionizing flux in the later stages of circumstellar disc evolution provides an important new constraint on disc photoevaporation models.     

Evolution of Class 0 protostars: models versus observations

The rates at which mass accumulates into protostellar cores can now be predicted in numerical simulations. Our purpose here is to develop methods to compare the statistical properties of the predicted protostars with the observable parameters. This requires (1) an evolutionary scheme to convert numerically derived mass accretion rates into evolutionary tracks and (2) a technique to compare the tracks to the observed statistics of protostars. Here, we use a 3D Kolmogorov–Smirnov test to quantitatively compare model evolutionary tracks and observations of Class 0 protostars.
We find that the wide range of accretion functions and time-scales associated with gravoturbulent simulations naturally overcome difficulties associated with schemes that use a fixed accretion pattern. This implies that the location of a protostar on an evolutionary track does not precisely determine the present age or final accrued mass. Rather, we find that predictions of the final mass for protostars from observed   T _bol– L _bol  values are uncertain by a factor of 2 and that the bolometric temperature is not always a reliable measure of the evolutionary stage. Furthermore, we constrain several parameters of the evolutionary scheme and estimate a lifetime of Class 0 sources of  2–6 × 10⁴ yr  , which is related to the local free-fall time and thus to the local density at the onset of the collapse. Models with Mach numbers smaller than six are found to best explain the observational data. Generally, only a probability of 70 per cent was found that our models explain the current observations. This is caused by not well-understood selection effects in the observational sample and the simplified assumptions in the models.     

X‐shooting Herbig Ae/Be stars: Accretion probed by near‐infrared He I emission

The Herbig Ae/Be stars are intermediate mass pre‐main sequence stars that bridge the gap between the low mass T Tauri stars and the Massive Young Stellar Objects. In this mass range, the acting star forming mechanism switches from magnetically controlled accretion to an as yet unknown mechanism, but which is likely to be direct disk accretion onto the star. We observed a large sample of Herbig Ae/Be stars with X‐shooter to address this issue from a multi‐wavelength perspective. It is the largest such study to date, not only because of the number of objects involved, but also because of the large wavelength coverage from the blue to the near‐infrared. This allows many accretion diagnostics to be studied simultaneously. By correlating the various properties with mass, temperature and age, we aim to determine where and whether the magnetically controlled mass accretion mechanism halts and the proposed direct disk accretion takes over. Here, we will give an overview of the background, present some observations and discuss our initial results. We will introduce a new accretion diagnostic for the research of Herbig Ae/Be stars, the HeI 1.083 μm line (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Possible evidence for the disc origin for the powering of jets in Sgr A* and nearby elliptical galaxies

Recent VLBA observation indicates the existence of an elongated (jet) structure in the compact radio source Sgr A*. This is hard to explain in the context of advection-dominated accretion flow (ADAF) model for this source. On the other hand, the mass accretion rate favoured by ADAF is 10–20 times smaller than that favoured by the hydrodynamical simulation based on Bondi capture. If the latter were adopted, the predicted radio flux would significantly exceed the observation. A similar situation exists in the case of nearby giant ellipticals, where the canonical ADAF model – the widely assumed standard model for these sources – also significantly overpredicts the radio flux. Based on these facts, in this paper we propose a truncated ADAF model for Sgr A* and three ellipticals M87, NGC 4649 and NGC 4636. We assume that the accretion disc is truncated at a certain radius R _tr within which the jet forms by extracting the energy of the disc. The radio flux is greatly suppressed owing to the radiative truncation of the disc and the fits to the observational data are excellent. For example, for Sgr A*, the model fits the observational spectrum very well from radio including the 'excess' below the break frequency to hard X-ray under a high accretion rate near the simulation value, and the predicted size-frequency relationship is also in excellent agreement with the observation; for M87, the predicted upper limit of the jet location is 24 R _g, in excellent agreement with the observational result that the jet is formed on scales smaller than 30 R _g, and the ≈20 per cent variability at ∼1 keV – which is hard to explain in another model that succeeded in explaining the low radio flux of M87 – is also marginally interpreted. The success of the model supplies possible evidence for the disc rather than the hole origin for the powering of jets.     

Numerical simulations of astrophysical jets from Keplerian discs – III. The effects of mass loading

We present 2.5D time-dependent simulations of the non-linear evolution of non-relativistic outflows from the surface of Keplerian accretion discs. The gas is accelerated from the surface of the disc (which is a fixed platform in these simulations) into a cold corona in stable hydrostatic equilibrium. We explore the dependence of the resulting jet characteristics upon the mass loading of the winds. Two initial configurations of the threading disc magnetic field are studied: a potential field and a uniform vertical field configuration.
We show that the nature of the resulting highly collimated, jet-like outflows (steady or episodic) is determined by the mass load of the disc wind. The mass load controls the interplay between the collimating effects of the toroidal field and the kinetic energy density in the outflow. In this regard, we demonstrate that the onset of episodic behaviour of jets appears to be determined by the quantity     which compares the speed for a toroidal Alfvén wave to cross the diameter of the jet, with the flow speed v _p along the jet. This quantity decreases with increasing load. For sufficiently large N (small mass loads), disturbances appear to grow leading to instabilities and shocks. Knots are then generated and the outflow becomes episodic. These effects are qualitatively independent of the initial magnetic configuration that we employed and are probably generic to a wide variety of magnetized accretion disc models.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.