Relativistic Jets from Accretion Disks

The jets observed to emanate from many compact accreting objects may arise from the twisting of a magnetic field threading a differentially rotating accretion disk which acts to magnetically extract angular momentum and energy from the disk. Two main regimes have been discussed, hydromagnetic jets, which have a significant mass flux and have energy and angular momentum carried by both matter and electromagnetic field and, Poynting jets, where the mass flux is small and energy and angular momentum are carried predominantly by the electromagnetic field. Here, we describe recent theoretical work on the formation of relativistic Poynting jets from magnetized accretion disks. Further, we describe new relativistic, fully electromagnetic, particle-in-cell (PIC) simulations of the formation of jets from accretion disks. Analog Z-pinch experiments may help to understand the origin of astrophysical jets.     

An electric-current description of solar flares

The presently prevailing theories of solar flares rely on the hypothetical presence of magnetic flux tubes beneath the photosphere and the two subsequent hypotheses, their emergence above the photosphere and explosive magnetic reconnection, converting magnetic energy carried by the flux tubes to solar flare energy. In this paper, we discuss solar flares from an entirely different point of view, namely in terms of power supply by a dynamo process in the photosphere. By this process, electric currents flowing along the magnetic field lines are generated and the familiar ‘force-free’ fields or the ‘sheared’ magnetic fields are produced. Upward field-aligned currents thus generated are carried by downward streaming electrons; these electrons can excite hydrogen atoms in the chromosphere, causing the optical Hα flares or ‘low temperature flares’. It is thus argued that as the ‘force-free’ fields are being built up for the magnetic energy storage, a flare must already be in progress.     

Linkage between accretion disks and blazars

The magnetic field in an accretion disk is estimated assuming that all of the angular momentum within prescribed accretion disk radii is removed by a jet. The magnetic field estimated at the base of the jet is extrapolated to the blazar emission region using a model for a relativistic axisymmetric jet combined with some simplifying assumptions based on the relativistic nature of the flow. The extrapolated magnetic field is compared with estimates based upon the synchrotron and inverse Compton emission from three blazars, MKN 501, MKN 421 and PKS 2155-304. The magnetic fields evaluated from pure synchrotron self-Compton models are inconsistent with the magnetic fields extrapolated in this way. However, in two cases inverse Compton models in which a substantial part of the soft photon field is generated locally agree well, mainly because these models imply magnetic field strengths consistent with an important Poynting Flux component. This comparison is based on estimating the mass accretion rate from the jet energy flux. Further comparisons along these lines will be facilitated by independent estimates of the mass accretion rate in blazars and by more detailed models for jet propagation near the black hole.     

The Dynamics of Gas around Active Galactic Nuclei

This paper examines the relationship between the accretion flows into the accretion disk around the black holes of active galactic nuclei (AGN), and the jets which they produce. A large-scale accretion flow around the nucleus is proposed as the thick dusty toroid structure of the unified model. Physically, this is similar to the Kahn ‘cocoon star’ model. High (super-Eddington) accretion rates into the accretion disk are needed. It is assumed that all black holes produce relativistic jets. However, entrainment of thermal (wind) material into the jet determines the AGN class, radio-loud or radio quiet. The jet interacts with the ambient medium as it expands, generating a cocoon of fast radiative or partially-radiative shocks around it. Such a model can explain both the radio properties and the emission line properties over a wide variety of AGN classes. This revised version was published online in July 2006 with corrections to the Cover Date.     

An essay on sunspots and solar flares

The presently prevailing theories of sunspots and solar flares rely on the hypothetical presence of magnetic flux tubes beneath the photosphere and the two subsequent hypotheses, their emergence above the photosphere and explosive magnetic reconnection, converting magnetic energy carried by the flux tubes for solar flare energy.In this paper, we pay attention to the fact that there are large-scale magnetic fields which divide the photosphere into positive and negative (line-of-sight) polarity regions and that they are likely to be more fundamental than sunspot fields, as emphasized most recently by McIntosh (1981). A new phenomenological model of the sunspot pair formation is then constructed by considering an amplification process of these largescale fields near their boundaries by shear flows, including localized vortex motions. The amplification results from a dynamo process associated with such vortex flows and the associated convergence flow in the largescale fields.This dynamo process generates also some of the familiar “force-free” fields or the “sheared” magnetic fields in which the magnetic field-aligned currents are essential. Upward field-aligned currents generated by the dynamo process are carried by downward streaming electrons which are expected to be accelerated by an electric potential structure; a similar structure is responsible for accelerating auroral electrons in the magnetosphere. Depending on the magnetic field configuration and the shear flows, the current-carrying electrons precipitate into different geometrical patterns, causing circular flares, umbral flares, two-ribbon flares, etc. Thus, it is suggested that “low temperature flares” are directly driven by the photospheric dynamo process.     

Isolated versus common envelope dynamos in planetary nebula progenitors

The origin, evolution and role of magnetic fields in the production and shaping of proto-planetary nebulae (PPNe) and planetary nebulae (PNe) are a subject of active research. Most PNe and PPNe are axisymmetric with many exhibiting highly collimated outflows; however, it is important to understand whether such structures can be generated by isolated stars or require the presence of a binary companion. Towards this end, we study a dynamical, large-scale α−Ω interface dynamo operating in a 3.0 M_⊙ Asymptotic Giant Branch (AGB) star in both an isolated setting and a setting in which a low-mass companion is embedded inside the envelope. The back reaction of the fields on the shear is included and differential rotation and rotation deplete via turbulent dissipation and Poynting flux. For the isolated star, the shear must be resupplied in order to sufficiently sustain the dynamo. Furthermore, we investigate the energy requirements that convection must satisfy to accomplish this by analogy to the Sun. For the common envelope case, a robust dynamo results, unbinding the envelope under a range of conditions. Two qualitatively different types of explosion may arise: (i) magnetically induced, possibly resulting in collimated bipolar outflows and (ii) thermally induced from turbulent dissipation, possibly resulting in quasi-spherical outflows. A range of models is presented for a variety of companion masses.     

Magnetic Energy Storage and Current Density Distributions for Different Force-Free Models

In the last decades, force-free-field modelling has been used extensively to describe the coronal magnetic field and to better understand the physics of solar eruptions at different scales. Especially the evolution of active regions has been studied by successive equilibria in which each computed magnetic configuration is subject to an evolving photospheric distribution of magnetic field and/or electric-current density. This technique of successive equilibria has been successful in describing the rate of change of the energetics for observed active regions. Nevertheless the change in magnetic configuration due to the increase/decrease of electric current for different force-free models (potential, linear and nonlinear force-free fields) has never been studied in detail before. Here we focus especially on the evolution of the free magnetic energy, the location of the excess of energy, and the distribution of electric currents in the corona. For this purpose, we use an idealised active region characterised by four main polarities and a satellite polarity, allowing us to specify a complex topology and sheared arcades to the coronal magnetic field but no twisted flux bundles. We investigate the changes in the geometry and connectivity of field lines, the magnetic energy and current-density content as well as the evolution of null points. Increasing the photospheric current density in the magnetic configuration does not dramatically change the energy-storage processes within the active region even if the magnetic topology is slightly modified. We conclude that for reasonable values of the photospheric current density (the force-free parameter α<0.25 Mm⁻¹), the magnetic configurations studied do change but not dramatically: i) the original null point stays nearly at the same location, ii) the field-line geometry and connectivity are slightly modified, iii) even if the free magnetic energy is significantly increased, the energy storage happens at the same location. This extensive study of different force-free models for a simple magnetic configuration shows that some topological elements of an observed active region, such as null points, can be reproduced with confidence only by considering the potential-field approximation. This study is a preliminary work aiming at understanding the effects of electric currents generated by characteristic photospheric motions on the structure and evolution of the coronal magnetic field.     

Magnetohydrodynamic instabilities and turbulence in accretion disks

In this paper we review the possibilities for magnetohydrodynamic processes to handle the angular momentum transport in accretion disks. Traditionally the angular momentum transport has been considered to be the result of turbulent viscosity in the disk, although the Keplerian flow in accretion disks is linearly stable towards hydrodynamic perturbations. It is on the other hand linearly unstable to some magnetohydrodynamic (MHD) instabilities. The most important instabilities are the Parker and Balbus-Hawley instabilities that are related to the magnetic buoyancy and the shear flow, respectively. We discuss these instabilities not only in the traditional MHD framework, but also in the context of slender flux tubes, that reduce the complexity of the problem while keeping most of the stability properties of the complete problem. In the non-linear regime the instabilities produce turbulence. Recent numerical simulations describe the generation of magnetic fields by a dynamo in the resulting turbulent flow. Eventually such a dynamo may generate a global magnetic field in the disk. The relation of the MHD-turbulence to observations of accretion disks is still obscure. It is commonly believed that magnetic fields can be highly efficient in transporting the angular momentum, but emission lines, short-time scale variability and non-thermal radiation, which a stellar astronomer would take as signs of magnetic variability, are more commonly observed during periods of low accretion rates. Received October 12, 1995 / Accepted November 16, 1995     

Viscous-resistive ADAF with a general large-scale magnetic field

We have studied the structure of hot accretion flow bathed in a general large-scale magnetic field. We have considered magnetic parameters , where are the Alfvén sound speeds in three direction of cylindrical coordinate (r,φ,z). The dominant mechanism of energy dissipation is assumed to be the magnetic diffusivity due to turbulence and viscosity in the accretion flow. Also, we adopt a more realistic model for kinematic viscosity (ν=αc _s H), with both c _s and H as a function of magnetic field. As a result in our model, the kinematic viscosity and magnetic diffusivity (η=η ₀ c _s H) are not constant. In order to solve the integrated equations that govern the behavior of the accretion flow, a self-similar method is used. It is found that the existence of magnetic resistivity will increase the radial infall velocity as well as sound speed and vertical thickness of the disk. However the rotational velocity of the disk decreases by the increase of magnetic resistivity. Moreover, we study the effect of three components of global magnetic field on the structure of the disk. We found out that the radial velocity and sound speed are Sub-Keplerian for all values of magnetic field parameters, but the rotational velocity can be Super-Keplerian by the increase of toroidal magnetic field. Also, Our numerical results show that all components of magnetic field can be important and have a considerable effect on velocities and vertical thickness of the disk.     

A new expression for cylindrically-symmetric magnetic fields in coronal magnetic loops

We present a set of cylindrically-symmetric force-free magnetic fields with non-constant scalar function scalar. We found that the kink instability of the fields can be suppressed by reducing the length of the flux tube. By using the pressure profile in coronal magnetic loops obtained on the basis of the observational data, and by neglecting the effect of gravity, these force-free fields ars modified to non-force-free ones. For the plasma of finite conductivity the time and space dependent magnetic fields are obtained, and the ohmic dissipation per unit volume per second is calculated. For the magnetic fields, presented in the investigation, it is also found that, due to the large electrical conductivity of the plasma, the ohmic dissipation is negligable in comparison to the conduction and the radiation loss. Hence, for the energy equilibrium in a coronal loop, the contribution of ohmic dissipation is insignificant.     

Can self-organized critical accretion disks generate a log-normal emission variability in AGN?

Active Galactic Nuclei (AGN), such as Seyfert galaxies, quasars, etc., show light variations in all wavelength bands, with various amplitude and in many time scales. The variations usually look erratic, not periodic nor purely random. Many of these objects also show lognormal flux distribution and RMS–flux relation and power law frequency distribution. So far, the lognormal flux distribution of black hole objects is only observational facts without satisfactory explanation about the physical mechanism producing such distribution in the accretion disk. One of the most promising models based on cellular automaton mechanism has been successful in reproducing PSD (Power Spectral Density) of the observed objects but could not reproduce lognormal flux distribution. Such distribution requires the existence of underlying multiplicative process while the existing SOC models are based on additive processes. A modified SOC model based on cellular automaton mechanism for producing lognormal flux distribution is presented in this paper. The idea is that the energy released in the avalanche and diffusion in the accretion disk is not entirely emitted instantaneously as in the original cellular automaton model. Some part of the energy is kept in the disk and thus increase its energy content so that the next avalanche will be in higher energy condition and will release more energy. The later an avalanche occurs, the more amount of energy is emitted to the observers. This can provide multiplicative effects to the flux and produces lognormal flux distribution.     

Conservation of both current and helicity in a quadrupolar model for solar flares

A model for a solar flare, involving magnetic reconnection transferring flux and current between current-carrying magnetic loops connecting two pairs of footpoints, is generalized to include conservation of magnetic helicity during reconnection, as well as conservation of current at all four footpoints. For a set of force-free loops, with the ith loop having flux F _i and current I _i, the self and mutual helicities are proportional to the self and mutual inductances with the constant of proportionality determined by α_i=F _i/μ₀ I _i. In a constant-α model, the change in magnetic energy is proportional to the change in helicity, and conservation of helicity implies conservation of magnetic energy, so that a flare cannot occur. In a quadrupolar model, with α₁>α₂ initially, α₁ increases and α₂ decreases when flux and current are transferred from loops 1 and 2 to loops 3 and 4. A model that conserves both current and helicity is constructed; it depends on the initial αs, and otherwise is somewhat simpler than when helicity is neglected.     

Coronal activity from dynamos in astrophysical rotators

We show that a steady mean-field dynamo in astrophysical rotators leads to an outflow of relative magnetic helicity and thus magnetic energy available for particle and wind acceleration in a corona. The connection between energy and magnetic helicity arises because mean-field generation is linked to an inverse cascade of magnetic helicity. To maintain a steady state in large magnetic Reynolds number rotators, there must then be an escape of relative magnetic helicity associated with the mean field, accompanied by an equal and opposite contribution from the fluctuating field. From the helicity flow, a lower limit on the magnetic energy deposited in the corona can be estimated. Steady coronal activity including the dissipation of magnetic energy, and formation of multi-scale helical structures therefore necessarily accompanies an internal dynamo. This highlights the importance of boundary conditions which allow this to occur for non-linear astrophysical dynamo simulations. Our theoretical estimate of the power delivered by a mean-field dynamo is consistent with that inferred from observations to be delivered to the solar corona, the Galactic corona, and Seyfert 1 AGN coronae.     

Jet-driven disk accretion in low luminosity AGN?

We explore an accretion model for low luminosity AGN (LLAGN) that attributes the low radiative output to a low mass accretion rate, , rather than a low radiative efficiency. In this model, electrons are assumed to drain energy from the ions as a result of collisionless plasma microinstabilities. Consequently, the accreting gas collapses to form a geometrically thin disk at small radii and is able to cool before reaching the black hole. The accretion disk is not a standard disk, however, because the radial disk structure is modified by a magnetic torque which drives a jet and which is primarily responsible for angular momentum transport. We also include relativistic effects. We apply this model to the well known LLAGN M87 and calculate the combined disk-jet steady-state broadband spectrum. A comparison between predicted and observed spectra indicates that M87 may be a maximally spinning black hole accreting at a rate of ∼10⁻³ M _⊙ yr⁻¹. This is about 6 orders of magnitude below the Eddington rate for the same radiative efficiency. Furthermore, the total jet power inferred by our model is in remarkably good agreement with the value independently deduced from observations of the M87 jet on kiloparsec scales.     

An unconventional accretion disc dynamo

In the light of recent results from numerical simulations of accretion disc MHD turbulence, we revisit the problem of the configuration of large-scale magnetic fields resulting from an α Ω dynamo operating in a thin accretion disc. In particular, we analyse the consequences of the peculiar sign of the α -effect suggested by numerical simulations . We determine the symmetry of the fastest-growing modes in the kinematic dynamo approximation and, in the framework of an ' α -quenched' dynamo model, study the evolution of the magnetic field. We find that the resulting field for this negative polarity of the α -effect generally has dipole symmetry with respect to the disc midplane, although the existence of an equilibrium configuration depends on the properties of the turbulence. The role of magnetic field dragging is discussed and, finally, the presence of an external uniform magnetic field is included to address the issue of magneto centrifugal wind launching from accretion discs.     

Jet-driven disk accretion in low luminosity AGN?*

We explore an accretion model for low luminosity AGN (LLAGN) that attributes the low radiative output to a low mass accretion rate, , rather than a low radiative efficiency. In this model, electrons are assumed to drain energy from the ions as a result of collisionless plasma microinstabilities. Consequently, the accreting gas collapses to form a geometrically thin disk at small radii and is able to cool before reaching the black hole. The accretion disk is not a standard disk, however, because the radial disk structure is modified by a magnetic torque which drives a jet and which is primarily responsible for angular momentum transport. We also include relativistic effects. We apply this model to the well known LLAGN M87 and calculate the combined disk-jet steady-state broadband spectrum. A comparison between predicted and observed spectra indicates that M87 may be a maximally spinning black hole accreting at a rate of ∼10⁻³  M _⊙ yr⁻¹. This is about 6 orders of magnitude below the Eddington rate for the same radiative efficiency. Furthermore, the total jet power inferred by our model is in remarkably good agreement with the value independently deduced from observations of the M87 jet on kiloparsec scales. ^* This paper has previously been published in Astrophysics and Space Science, vol. 310:3–4.     

The Alpha Effect and the Observed Twist and Current Helicity of Solar Magnetic Fields

We present a straightforward comparison of model calculations for the α-effect, helicities, and magnetic field line twist in the solar convection zone with magnetic field observations at atmospheric levels. The model calculations are carried out in a mixing-length approximation for the turbulence with a profile of the solar internal rotation rate obtained from helioseismic inversions. The magnetic field data consist of photospheric vector magnetograms of 422 active regions for which spatially-averaged values of the force-free twist parameter and of the current helicity density are calculated, which are then used to determine latitudinal profiles of these quantities. The comparison of the model calculations with the observations suggests that the observed twist and helicity are generated in the bulk of the convection zone, rather than in a layer close to the bottom. This supports two-layer dynamo models where the large-scale toroidal field is generated by differential rotation in a thin layer at the bottom while the α-effect is operating in the bulk of the convection zone. Our previous observational finding was that the moduli of the twist factor and of the current helicity density increase rather steeply from zero at the equator towards higher latitudes and attain a certain saturation at about 12 – 15^∘. In our dynamo model with algebraic nonlinearity, the increase continues, however, to higher latitudes and is more gradual. This could be due to the neglect of the coupling between small-scale and large-scale current and magnetic helicities and of the latitudinal drift of the activity belts in the model.     

Magnetic flares in Active Galactic Nuclei: modeling the iron Kα line

The X‐ray spectra of Active Galactic Nuclei (AGN) are complex and vary rapidly in time as seen in recent observations. Magnetic flares above the accretion disk can account for the extreme variability of AGN. They also explain the observed iron Kα fluorescence lines. We present radiative transfer modeling of the X‐ray reflection due to emission from magnetic flares close to the marginally stable orbit. The hard X‐ray primary radiation coming from the flare source illuminates the accretion disk. A Compton reflection/reprocessed component coming from the disk surface is computed for different emission directions. We assume that the density structure remains adjusted to the hydrostatic equilibrium without external illumination because the flare duration is only a quarter‐orbit. The model takes into account the variations of the incident radiation across the hot spot underneath the flare source. The integrated spectrum seen by a distant observer is computed for flares at different orbital phases close to the marginally stable orbit of a Schwarzschild black hole and of a maximally rotating Kerr black hole. The calculations include relativistic and Doppler corrections of the spectra using a ray tracing technique. We explore the practical possibilities to map out the azimuthal irradiation pattern of the inner accretion disks and conclude that the next generation of X‐ray satellites should reveal this structure from iron Kα line profiles and X‐ray lightcurves. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)