We present, for the first time, the structure of the axisymmetric force-free magnetosphere of an aligned rotating magnetic dipole, in the case in which there exists a sufficiently large charge density (whose origin we do not question) to satisfy the ideal MHD condition, ${\bf E\cdot B}=0$, everywhere. The unique distribution of electric current along the open magnetic field lines which is required for the solution to be continuous and smooth is obtained numerically. With the geometry of the field lines thus determined we compute the dynamics of the associated MHD wind. The main result is that the relativistic outflow contained in the magnetosphere is not accelerated to the extremely relativistic energies required for the flow to generate gamma rays. We expect that our solution will be useful as the starting point for detailed studies of pulsar magnetospheres under more general conditions, namely when either the force-free and/or the ideal MHD condition ${\bf E\cdot B}=0$ are not valid in the entire magnetosphere. Based on our solution, we consider that the most likely positions of such an occurrence are the polar cap, the crossings of the zero space charge surface by open field lines, and the return current boundary, but not the light cylinder.

We have applied axisymmetric MHD simulations to investigate the impact of the accretion disk magnetic flux profile on the jet collimation. Using the ZEUS-3D code modified for magnetic diffusivity, our simulations evolve from an initial hydrostatic equilibrium state in a force-free magnetic field. Considering a power law for the disk poloidal magnetic field profile Bp ~ r^{-mu} and for the disk wind density profile rho ~ r^{-mu_rho} we performed a systematic study over a wide parameter range mu and mu_rho. We find a degree of collimation (ratio of mass flow rates in axial and lateral direction) decreasing for steeper disk magnetic field profiles (increasing mu). Varying the total magnetic flux doesn't change the degree of jet collimation substantially, it only affects the time scale of outflow evolution and the terminal jet speed. As our major result we find a general relation between the collimation degree with the disk wind magnetization power law exponent. Outflows with high collimation degree resulting from a flat disk magnetic field profile tend to be unsteady, producing axially propagating knots as discussed earlier. Depending slightly on the inflow density profile this unsteady behavior sets in for mu < 0.4. We also performed simulations of jet formation with artificially enhanced decay of the toroidal magnetic field in order to investigate the idea of a purely "poloidal collimation" discussed in the literature. These outflows remain weakly collimated and propagate with lower velocity. Thanks to our large numerical grid size (7x14 AU for protostars), we may apply our results to recently observed hints of jet rotation (DG Tau) indicating a relatively flat disk magnetic field profile, mu ~ 0.5. In general, our results are applicable to both stellar and extragalactic sources of MHD jets.Comment: accepted by ApJ, high resolution version under www.mpia-hd.mpg.de/homes/fendt

We investigate the launching mechanism of relativistic jets from black hole sources, in particular the strong winds from the surrounding accretion disk. Numerical investigations of the disk wind launching—the simulation of the accretion–ejection transition—have so far almost only been done for nonrelativistic systems. From these simulations we know that resistivity, or magnetic diffusivity, plays an important role for the launching process. Here we extend this treatment to general relativistic magnetohydrodynamics (GR-MHD), applying the resistive GR-MHD code rHARM. Our model setup considers a thin accretion disk threaded by a large-scale open magnetic field. We run a series of simulations with different Kerr parameter, field strength, and diffusivity level. Indeed, we find strong disk winds with, however, mildly relativistic speed, the latter most probably due to our limited computational domain. Further, we find that magnetic diffusivity lowers the efficiency of accretion and ejection, as it weakens the efficiency of the magnetic lever arm of the disk wind. As a major driving force of the disk wind we disentangle the toroidal magnetic field pressure gradient; however, magnetocentrifugal driving may also contribute. Black hole rotation in our simulations suppresses the accretion rate owing to an enhanced toroidal magnetic field pressure that seems to be induced by frame dragging. Comparing the energy fluxes from the Blandford–Znajek-driven central spine and the surrounding disk wind, we find that the total electromagnetic energy flux is dominated by the total matter energy flux of the disk wind (by a factor of 20). The kinetic energy flux of the matter outflow is comparatively small and comparable to the Blandford–Znajek electromagnetic energy flux.