We continue our investigation of particle acceleration in the pulsar equatorial current sheet (ECS) that began with Contopoulos (2019) and Contopoulos & Stefanou (2019). Our basic premise has been that the charge carriers in the current sheet originate in the polar caps as electron-positron pairs, and are carried along field lines that enter the equatorial current sheet beyond the magnetospheric Y-point. In this work we investigate further the charge replenishment of the ECS. We discovered that the flow of pairs from the rims of the polar caps cannot supply both the electric charge and the electric current of the ECS. The ECS must contain an extra amount of positronic (or electronic depending on orientation) electric current that originates in the stellar surface and flows outwards along the separatrices. We develop an iterative hybrid approach that self-consistently combines ideal force-free electrodynamics in the bulk of the magnetosphere with particle acceleration along the ECS. We derive analytic approximations for the orbits of the particles, and obtain the structure of the pulsar magnetosphere for various values of the pair-formation multiplicity parameter kappa. For realistic values kappa >> 1, the magnetosphere is practically indistinguishable from the ideal force-free one, and therefore, the calculation of the spectrum of high-energy radiation must be based on analytic approximations for the distribution of the accelerating electric field in the ECS.

Force-free pulsar magnetospheres develop a large scale poloidal electric current circuit that flows along open magnetic field lines from the neutron star to the termination shock. The electric current closes through the interior of the neutron star where it provides the torque that spins-down the star. In the present work, we study the internal electric current in an axisymmetric rotator. We evaluate the path of the electric current by requiring the minimization of internal Ohmic losses. We find that, in millisecond pulsars, the current reaches the base of the crust, while in pulsars with periods of a few seconds, the bulk of the electric current does not penetrate deeper than about $100$ m. The region of maximum spin-down torque in millisecond pulsars is the base of the crust, while in slowly spinning ones it is the outer crust. We evaluate the corresponding Maxwell stresses and find that, in typical rotation-powered radio pulsars, they are well below the critical stress that can be sustained by the crust. For magnetar-level fields, the Maxwell stresses near the surface are comparable to the critical stress and may lead to the decoupling of the crust from the rest of the stellar rotation.

Contopoulos 2019 proposed that a dissipation zone develops in the magnetosphere of young pulsars at the edge of the closed-line region beyond the light cylinder. This is necessary in order to supply the charge carriers that will establish current closure through the equatorial and separatrix current-sheets. In the present work, we propose to investigate in greater detail this region with a simplified model that we would like to call the `ring-of-fire'. According to this simple model, the dissipation zone is a narrow reconnection layer where electrons and positrons are accelerated inwards and outwards respectively along Speiser orbits that are deflected in the azimuthal direction by the pulsar rotation. After they exit the reconnection layer, the accelerated positrons form the positively charged equatorial current-sheet, and the accelerated electrons form the negatively charged separatrix current-sheet along the boundary of the closed-line region. During their acceleration, particles lose only a small part of their energy to radiation. Most of their energy is lost outside the dissipation region, in the equatorial and separatrix current sheets. Our simple model allows us to obtain high-energy spectra and efficiencies. The radiation emitted by the positrons in the equatorial current-sheet forms a very-high energy tail that extends up to the TeV range.

The detection of the pulsed $\sim 1 $~TeV gamma-ray emission from the Crab pulsar reported by MAGIC and VERITAS collaborations demands a substantial revision of existing models of particle acceleration in the pulsar magnetosphere. In this regard model independent restrictions on the possible production site of the VHE photons become an important issue. In this paper, we consider limitations imposed by the process of conversion of VHE gamma rays into $e^{\pm}$ pairs in the magnetic field of the pulsar magnetosphere. Photons with energies exceeding 1~TeV are effectively absorbed even at large distances from the surface of the neutron star. Our calculations of magnetic absorption in the force-free magnetosphere show that the twisting of the magnetic field due to the pulsar rotation makes the magnetosphere more transparent compared to the dipole magnetosphere. The gamma-ray absorption appears stronger for photons emitted in the direction of rotation than in the opposite direction. There is a small angular cone inside which the magnetosphere is relatively transparent and photons with energy $1.5$~TeV can escape from distances beyond $0.1$~light cylinder radius ($R_{\rm{lc}}$). The emission surface from where photons can be emitted in the observer's direction further restricts the sites of VHE gamma-ray production. For the observation angle $57^{\circ}$ relative to the Crab pulsar axis of rotation and the orthogonal rotation, the emission surface in the open field line region is located as close as $0.4\,R_{\rm{lc}}$ from the stellar surface for a dipole magnetic field, and $0.1\,R_{\rm{lc}}$ for a force-free magnetic field.

Astrophysical accretion flows are associated with energetic emission of radiation and outflows (winds and jets). Extensive observations of these two processes in X-ray binary outbursts are available. A convincing understanding of their dynamics remains, however, elusive. The main agent that controls the dynamics is believed to be a large scale magnetic field that threads the system. We propose that during the quiescent state, the field is held in place by a delicate balance between inward advection and outward diffusion through the accreting matter. We also propose that the source of the field is a growing toroidal electric current generated by the aberrated radiation pressure on the innermost plasma electrons in orbit around the central black hole. This is the astrophysical mechanism of the Cosmic Battery. When the return magnetic field outside the toroidal electric current diffuses through the surrounding disk, the disk magnetic field and its associated accretion rate gradually increase, thus leading the system to an outburst. After the central accretion flow approaches equipartition with radiation, it is disrupted, and the Cosmic Battery ceases to operate. The outward field diffusion is then reversed, magnetic flux reconnects with the flux accumulated around the central black hole and disappears. The magnetic field and the associated accretion rate slowly decrease, and the system is gradualy driven back to quiescence. We conclude that the action (or inaction) of the Cosmic Battery may be the missing key that will allow us to understand the long-term evolution of astrophysical accretion-ejection flows.

We present our view of the immediate environment of an astrophysical black hole. We discuss the origin of the electron-positron and electron-proton jet components, and the origin of the large-scale magnetic field. We show how the aberrated radiation field viewed by the rotating plasma generates an azimuthal induction electric field, and how the curl of this electric field steadily grows the magnetic field to equipartition over a few hundred million dynamical times. That timescale corresponds to weeks for stellar mass black holes, and million years for super-massive black holes. We conclude with numerical and observational confirmation of the action of a Cosmic Battery in accretion flows around astrophysical black holes.

Current closure in the pulsar magnetosphere holds the key to its structure. We must determine not only the global electric circuit, but also the source of its electric charge carriers. We address this issue with the minimum number of assumptions: a) The magnetosphere is everywhere ideal and force-free, except above the polar cap and in some finite part of the current sheet; and b) pairs are produced above the polar cap with multiplicity kappa. We show that a thin region of width delta ~ r_pc/2 kappa << r_pc along the rim of the polar cap provides all the charges that are needed in the equatorial and separatrix electric current sheet. These charges are transferred to the current sheet in a narrow dissipation zone just outside the magnetospheric Y-point. The maximum accelerating potential in this region is equal to the potential drop in the thin polar cap rim, which is approximately equal to 1/kappa times the potential drop from the center to the edge of the polar cap. The dissipated electromagnetic energy is approximately equal to 0.5/kappa times the total pulsar spindown energy loss. Our framework allows to calculate the high energy emission in terms of the pair multiplicity.

Among a number of active galactic nuclei (AGNs) that drive ionized outflows in X-rays, a low-redshift (z = 0.184) quasar, PDS 456, is long known to exhibit one of the exemplary ultra-fast outflows (UFOs). However, the physical process of acceleration mechanism is yet to be definitively constrained. In this work, we model the variations of the Fe K UFO properties in PDS 456 over many epochs in X-ray observations in the context of magnetohydrodynamic (MHD) accretion-disk winds employed in our earlier studies of similar X-ray absorbers. We applied the model to the 2013/2014 XMM-Newton/NuSTAR spectra to determine the UFO's condition; namely, velocity, ionization parameter, column density and equivalent width (EW). Under some provisions on the dependence of X-ray luminosity on the accretion rate applicable to near-Eddington state, our photoionization calculations, coupled to a 2.5-dimensional MHD-driven wind model, can further reproduce the observed correlations of the UFO velocity and the anti-correlation of its EW with X-ray strength of PDS 456. This work demonstrates that UFOs, even without radiative pressure, can be driven as an extreme case purely by magnetic interaction while also producing the observed spectrum and correlations.

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

We analyze a 900-ks stacked Chandra/HETG spectrum of NGC 3783 in the context of magnetically-driven accretion-disk wind models in an effort to provide tight constraints on the global conditions of the underlying absorbers. Motivated by the earlier measurements of its absorption measure distribution (AMD) indicating X-ray-absorbing ionic columns that decrease slowly with decreasing ionization parameter, we employ 2D magnetohydrodynamic (MHD) disk-wind models to describe the global outflow. We compute its photoionization structure along with the wind kinematic properties allowing us to further calculate in a self-consistent fashion the shapes of the major X-ray absorption lines. With the wind radial density profile determined by the AMD, the profiles of the ensemble of the observed absorption features are determined by the two global parameters of the MHD wind; i.e. disk inclination \theta_obs and wind density normalization n_o. Considering the most significant absorption features in the (~1.8A-20A) range, we show that the MHD-wind is best described by n(r)~6.9e11(r/ro)^-1.15 [cm^-3] and \theta_obs=44deg. We argue that winds launched by X-ray heating, radiation pressure or even MHD winds but with steeper radial density profiles are strongly disfavored by data. Considering the properties of Fe K band absorption features (i.e. Fe xxv and Fe xxvi), while typically prominent in the AGN X-ray spectra, they appear to be weak in NGC 3783. For the specific parameters of our model obtained by fitting the AMD and the rest of absorption features, these features are found to be weak in agreement with observation.

In view of the upcoming observations with the Event Horizon Telescope (EHT), we present our thoughts on the immediate environment of an astrophysical black hole. We are concerned that two approximations used in general relativistic magnetohydrodynamic numerical simulations, namely numerical density floors implemented near the base of the black hole jet, and a magnetic field that comes from large distances, may mislead our interpretation of the observations. We predict that three physical processes will manifest themselves in EHT observations, namely dynamic pair formation just above the horizon, electromagnetic energy dissipation along the boundary of the black hole jet, and a region of weak magnetic field separating the black hole jet from the disk wind.

Observations of X-ray absorption lines in magnetically driven disk winds around black hole binaries and active galactic nuclei yield a universal radial density profile rho proportional to r^{-1.2} in the wind. This is in disagreement with the standard Blandford & Payne profile rho_BP proportional to r^{-1.5} expected when the magnetic field is neither advected nor diffusing through the accretion disk. In order to account for this discrepancy, we establish a new paradigm for magnetically driven astrophysical winds according to which the large scale ordered magnetic field that threads the disk is continuously generated by the Cosmic Battery around the inner edge of the disk and continuously diffuses outward. We obtain self-similar solutions of such magnetically advected winds (MAW) and discuss their observational ramifications.

We implement the KORAL code to perform two sets of very long general relativistic radiation magnetohydrodynamic simulations of an axisymmetric optically thin magnetized flow around a non-rotating black hole: one with a new term in the electromagnetic field tensor due to the radiation pressure felt by the plasma electrons on the comoving frame of the electron-proton plasma, and one without. The source of the radiation is the accretion flow itself. Without the new term, the system evolves to a standard accretion flow due to the development of the magneto-rotational instability (MRI). With the new term, however, the system eventually evolves to a magnetically arrested state (MAD) in which a large scale jet-like magnetic field threads the black hole horizon. Our results confirm the secular action of the Cosmic Battery in accretion flows around astrophysical black holes.

The recent discoveries of pulsed X-ray emission from three ultraluminous X-ray (ULX) sources have finally enabled us to recognize a subclass within the ULX class: the great pretenders, neutron stars (NSs) that appear to emit X-ray radiation at isotropic luminosities $L_X = 7\times 10^{39}$~erg~s$^{-1}-1\times 10^{41}$~erg~s$^{-1}$ only because their emissions are strongly beamed toward our direction and our sight lines are offset by only a few degrees from their magnetic-dipole axes. The three known pretenders appear to be stronger emitters than the presumed black holes of the ULX class, such as Holmberg II \& IX X-1, IC10 X-1, and NGC300 X-1. For these three NSs, we have adopted a single reasonable assumption, that their brightest observed outbursts unfold at the Eddington rate, and we have calculated both their propeller states and their surface magnetic-field magnitudes. We find that the results are not at all different from those recently obtained for the Magellanic Be/X-ray pulsars: the three NSs reveal modest magnetic fields of about 0.3-0.4~TG and beamed propeller-line X-ray luminosities of $\sim 10^{36-37}$~erg~s$^{-1}$, substantially below the Eddington limit.

Black hole accretion disks appear to produce invariably plasma outflows that result in blue-shifted absorption features in their spectra. The X-ray absorption-line properties of these outflows are quite diverse, ranging in velocity from non-relativistic ($\sim 300$ km/sec) to sub-relativistic ($\sim 0.1c$ where $c$ is the speed of light) and a similarly broad range in the ionization states of the wind plasma. We report here that semi-analytic, self-similar magnetohydrodynamic (MHD) wind models that have successfully accounted for the X-ray absorber properties of supermassive black holes, also fit well the high-resolution X-ray spectrum of the accreting stellar-mass black hole, GRO J1655-40. This provides an explicit theoretical argument of their MHD origin (aligned with earlier observational claims) and supports the notion of a universal magnetic structure of the observed winds across all known black hole sizes.

In order to investigate the importance of dissipation in the pulsar magnetosphere we decided to combine Force-Free with Aristotelian Electrodynamics. We obtain solutions that are ideal (non-dissipative) everywhere except in an equatorial current sheet where Poynting flux from both hemispheres converges and is dissipated into particle acceleration and radiation. We find significant dissipative losses (up to about 50% of the pulsar spindown luminosity), similar to what is found in global PIC simulations in which particles are provided only on the stellar surface. We conclude that there might indeed exist two types of pulsars, strongly dissipative ones with particle injection only from the stellar surface, and ideal (weakly dissipative) ones with particle injection in the outer magnetosphere and in particular at the Y-point.

This is our first study of the magnetic Rayleigh-Taylor instability at the inner edge of an astrophysical disk around a central back hole. We derive the equations governing small-amplitude oscillations in general relativistic ideal magnetodydrodynamics and obtain a criterion for the onset of the instability. We suggest that static disk configurations where magnetic field is held by the disk material are unstable around a Schwarzschild black hole. On the other hand, we find that such configurations are stabilized by the spacetime rotation around a Kerr black hole. We obtain a crude estimate of the maximum amount of poloidal magnetic flux that can be accumulated around the center, and suggest that it is proportional to the black hole spin. Finally, we discuss the astrophysical implications of our result for the theoretical and observational estimations of the black hole jet power.

The spins of stellar-mass black holes (BHs) and the power outputs of their jets are measurable quantities. Unfortunately, the currently employed methods do not agree and the results are controversial. Two major issues concern the measurements of BH spin and beam (jet) power. The former issue can be resolved by future observations. But the latter issue can be resolved now, if we pay attention to what is expected from theoretical considerations. The question of whether a correlation has been found between the power outputs of few objects and the spins of their BHs is moot because BH beam power does not scale with the square of the spin of the BH. We show that the theoretical BH beam power is a strongly nonlinear function of spin that cannot be approximated by a quadratic relation, as is generally stated when the influence of the magnetic field is not accounted for in the \cite{bla77} model. The BH beam power of ballistic jets should scale a lot more steeply with BH spin irrespective of the magnetic field assumed to thread the horizon and the spin range considered. This behavior may already be visible in the analyses of radio observations by \cite{nar12} and \cite{rus13}. In agreement with previous studies, we also find that the power output that originates in the inner regions of the surrounding accretion disks is higher than that from the BHs and it cannot be ignored in investigations of continuous compact jets from these systems.

Helical magnetic fields embedded in the jets of active galactic nuclei (AGNs) are required by the broad range of theoretical models that advocate for electromagnetic launching of the jets. In most models, the direction of the magnetic field is random, but if the axial field is generated by a Cosmic Battery generated by current in the direction of rotation in the accretion disk, there is a correlation between the directions of the spin of the AGN accretion disk and of the axial field, which leads to a specific direction for the axial electric current, azimuthal magnetic field, and the resulting observed transverse Faraday-rotation (FR) gradient across the jet, due to the systematic change in the line-of-sight magnetic field. We consider new observational evidence for the presence of a nested helical magnetic-field structure such as would be brought about by the operation of the Cosmic Battery, and make predictions about the expected behavior of transverse FR gradients observed on decaparsec and kiloparsec scales. The collected results can be understood if the dominant azimuthal field on parsec scales corresponds to an axial electric current flowing inward along the jet, whereas the (weaker) dominant azimuthal field on kiloparsec scales corresponds to a outward-flowing current in the outer sheath of the jet and/or an extended disk wind.

We want to understand what drives magnetospheric dissipation in the equatorial current sheet. Numerical simulations have limitations and, unless we have a clear a priori understanding of the physical processes involved, their results can be misleading. We argue that the canonical pulsar magnetosphere is strongly dissipative and that a large fraction (up to 30-40% in an aligned rotator) of the spindown luminosity is redirected towards the equator where it is dissipated into particle acceleration and emission of radiation. We show that this is due to the failure of the equatorial electric current to cross the Y-point at the tip of the corotating zone.

Gamma-ray bursts (GRBs) are believed to be powered by the electromagnetic extraction of spin energy from a black hole endowed with a magnetic field supported by electric currents in a surrounding disk (Blandford & Znajek 1977). A generic feature of this mechanism is that, under certain fairly general assumptions, the energy loss rate decays exponentially. In this work, we are looking precisely for such exponential decay in the lightcurves of long duration GRBs observed with the XRT instrument on the Swift satellite. We found out that almost 30 % of XRT lightcurves show such behavior before they reach the afterglow plateau. According to Blandford & Znajek, the duration of the burst depends on the magnetic flux accumulated on the event horizon. This allows us to estimate the surface magnetic field of a possible progenitor. Our estimations are consistent with magnetic fields observed in Wolf-Rayet stars.

Gamma-ray bursts (GRBs) are violent explosions, coming from cosmological distances. They are detected in gamma-rays (also X-rays, UV, optical, radio) almost every day, and have typical durations of a few seconds to a few minutes. Some GRBs have been reported with extraordinary duration of 10^4 sec, the so-called Ultra Long GRBs. It has been debated whether these form a new distinct class of events or whether they are similar to long GRBs. According to Blandford & Znajek (1977), the spin energy of a rotating black hole can be extracted electromagnetically, should the hole be endowed with a magnetic field supported by electric currents in a surrounding disk. We argue that this can be the case for the central engines of GRBs and we show that the duration of the burst depends on the magnetic flux accumulated on the event horizon of the black hole. We thus estimate the surface magnetic field of a possible progenitor star, and we conclude that an Ultra Long GRB may originate from a progenitor star with a relatively low magnetic field.

We present a study of X-ray ionization of magnetohydrodynamic (MHD) accretion-disk winds in an effort to constrain the physics underlying the highly-ionized ultra-fast outflows (UFOs) inferred by X-ray absorbers often detected in various sub-classes of Seyfert active galactic nuclei (AGNs). Our primary focus is to show that magnetically-driven outflows are indeed physically plausible candidates for the observed outflows accounting for the AGN absorption properties of the present X-ray spectroscopic observations. Employing a stratified MHD wind launched across the entire AGN accretion disk, we calculate its X-ray ionization and the ensuing X-ray absorption line spectra. Assuming an appropriate ionizing AGN spectrum, we apply our MHD winds to model the absorption features in an {\it XMM-Newton}/EPIC spectrum of the narrow-line Seyfert, \pg. We find, through identifying the detected features with Fe K$\alpha$ transitions, that the absorber has a characteristic ionization parameter of $\log (\xi_c [erg~cm~s$^{-1}$]) \simeq 5-6$ and a column density on the order of $N_H \simeq 10^{23}$ cm$^{-2}$, outflowing at a characteristic velocity of $v_c/c \simeq 0.1-0.2$ (where $c$ is the speed of light). The best-fit model favors its radial location at $r_c \simeq 200 R_o$ ($R_o$ is the black hole innermost stable circular orbit), with an inner wind truncation radius at $R_{\rm t} \simeq 30 R_o$. The overall K-shell feature in the data is suggested to be dominated by \fexxv\ with very little contribution from \fexxvi\ and weakly-ionized iron, which is in a good agreement with a series of earlier analysis of the UFOs in various AGNs including \pg.

The aberrated radiation pressure at the inner edge of the accretion disk around an astrophysical black hole imparts a relative azimuthal velocity on the electrons with respect to the ions which gives rise to a ring electric current that generates large scale poloidal magnetic field loops. This is the Cosmic Battery established by Contopoulos and Kazanas in 1998. In the present work we perform realistic numerical simulations of this important astrophysical mechanism in advection-dominated accretion flows-ADAF. We confirm the original prediction that the inner parts of the loops are continuously advected toward the central black hole and contribute to the growth of the large scale magnetic field, whereas the outer parts of the loops are continuously diffusing outward through the turbulent accretion flow. This process of inward advection of the axial field and outward diffusion of the return field proceeds all the way to equipartition, thus generating astrophysically significant magnetic fields on astrophysically relevant timescales. We confirm that there exists a critical value of the magnetic Prandtl number between unity and 10 in the outer disk above which the Cosmic Battery mechanism is suppressed.

According to Blandford & Znajek (1977), energy can be extracted electromagnetically from a rotating black hole, should the latter be endowed with a magnetic field supported by electric currents in a surrounding disk. We show that exact models of black hole magnetospheres produce Poynting flux that decreases almost exponentially with time. We went through the Swift BAT-XRT lightcurves and identified a subclass of GRBs that exhibits a clear exponential decay over more than three orders of magnitude in flux (EDOHS GRBs). We estimate the energy given-off in the X-rays and discuss a possible correlation between the peak brightness of the X-ray prompt emission and its decay time. We investigate a possible application of this result in high redshift Cosmology.

We investigate the dynamics of radiation in an astrophysical accretion disk around a Kerr black hole. The source of the radiation is the accretion disk itself, and not a central object as in previous studies of the Poynting-Robertson effect. We generate numerical sky maps from photon trajectories that originate on the surface of the disk as seen from the inner edge of the disk at the position of the innermost stable circular orbit (ISCO). We investigate several accretion disk morphologies with a Shakura-Sunyaev surface temperature distribution. Finally, we calculate the electromotive source of the Cosmic Battery mechanism around the inner edge of the accretion disk and obtain characteristic timescales for the generation of astrophysical magnetic fields.

We investigate the structure of the steady-state force-free magnetosphere around a Kerr black hole in various astrophysical settings. The solution Psi(r,theta) depends on the distributions of the magnetic field line angular velocity omega(Psi) and the poloidal electric current I(Psi). These are obtained self-consistently as eigenfunctions that allow the solution to smoothly cross the two singular surfaces of the problem, the Inner Light Surface (ILS) inside the ergosphere, and the Outer Light Surface (OLS), which is the generalization of the pulsar light cylinder. Magnetic field configurations that cross both singular surfaces (e.g. monopole, paraboloidal) are uniquely determined. Configurations that cross only one light surface e.g. the artificial case of a rotating black hole embedded in a vertical magnetic field) are degenerate. We show that, similarly to pulsars, black hole magnetospheres naturally develop an electric current sheet that potentially plays a very important role in the dissipation of black hole rotational energy and in the emission of high-energy radiation.

In view of recent efforts to probe the physical conditions in the pulsar current sheet, we revisit the standard solution that describes the main elements of the ideal force-free pulsar magnetosphere. The simple physical requirement that the electric current contained in the current layer consists of the local electric charge moving outward at close to the speed of light, yields a new standard solution for the pulsar magnetosphere everywhere ideal force-free except in the current layer. The main elements of the new solution are a) the pulsar spindown rate of the aligned rotator is 23% times larger than that of the orthogonal vacuum rotator, b) only 60% of the magnetic flux that crosses the light cylinder opens up to infinity, c) the electric current closes along the other 40% which gradually converges to the equator, d) this transfers 40% of the total pulsar spindown energy flux in the equatorial current sheet which is then dissipated in the acceleration of particles and in high-energy electromagnetic radiation, e) there is no separatrix current layer. Our solution is a minimum free-parameter solution in that the equatorial current layer is electrostatically supported against collapse and thus does not require a thermal particle population. In that respect, it is one more step toward the development of a new standard solution. We discuss the implications for intermittent pulsars and long duration gamma-ray bursts. We conclude that the physical conditions in the equatorial current layer determine the global structure of the pulsar magnetosphere.

We explore the analogy between a rotating magnetized black hole and an axisymmetric pulsar and derive its electromagnetic spindown after its formation in the core collapse of a supermassive star. The spindown shows two characteristic phases, an early Blandford-Znajek phase that lasts a few hundred seconds, and a late pulsar-like afterglow phase that lasts much longer. During the first phase, the spindown luminosity decreases almost exponentially, whereas during the afterglow phase it decreases as t^{-a} with 1<a<1,5. We associate our findings with long duration gamma-ray bursts (GRB) and compare with observations.

We explore the poloidal structure of two-dimensional (2D) MHD winds in relation to their potential association with the X-ray warm absorbers (WAs) and the highly-ionized ultra-fast outflows (UFOs) in AGN, in a single unifying approach. We present the density $n(r,\theta)$, ionization parameter $\xi(r,\theta)$, and velocity structure $v(r,\theta)$ of such ionized winds for typical values of their fluid-to-magnetic flux ratio, $F$, and specific angular momentum, $H$, for which wind solutions become super-\Alfvenic. We explore the geometrical shape of winds for different values of these parameters and delineate the values that produce the widest and narrowest opening angles of these winds, quantities necessary in the determination of the statistics of AGN obscuration. We find that winds with smaller $H$ show a poloidal geometry of narrower opening angles with their \Alfven\ surface at lower inclination angles and therefore they produce the highest line of sight (LoS) velocities for observers at higher latitudes with the respect to the disk plane. We further note a physical and spatial correlation between the X-ray WAs and UFOs that form along the same LoS to the observer but at different radii, $r$, and distinct values of $n$, $\xi$ and $v$ consistent with the latest spectroscopic data of radio-quiet Seyfert galaxies. We also show that, at least in the case of 3C 111, the winds' pressure is sufficient to contain the relativistic plasma responsible for its radio emission. Stratified MHD disk-winds could therefore serve as a unique means to understand and unify the diverse AGN outflows.

We revisit the Blandford & Znajek (1977) process and solve the fundamental equation that governs the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the magnetic field angular velocity omega and the poloidal electric current I. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem, the inner `light surface' located inside the ergosphere, and the outer `light surface' which is the generalization of the pulsar light cylinder. We find the solution for the simplest possible magnetic field configuration, the split monopole, through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state axisymmetric force-free pulsar magnetosphere (Contopoulos, Kazanas & Fendt 1999). We obtain the rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of previous time dependent simulations. Furthermore, we discuss the physical applicability of magnetic field configurations that do not cross both `light surfaces'.

We present a new improved version of our force-free electrodynamics (FFE) numerical code in spherical coordinates that extrapolates the magnetic field in the inner solar corona from a photospheric vector magnetogram. The code satisfies the photospheric boundary condition and the condition divB=0 to machine accuracy. The performance of our method is evaluated with standard convergence parameters, and is found to be comparable to that of other nonlinear force-free extrapolations.

We present a unified model for the structure and appearance of accretion powered sources across their entire luminosity range from galactic X-ray binaries to luminous quasars, with emphasis on AGN and their phenomenology. Central to this model is the notion of MHD winds launched from the accretion disks that power these objects. These winds provide the matter that manifests as blueshifted absorption features in the UV and X-ray spectra of a large fraction of these sources; furthermore, their density distribution in the poloidal plane determines the "appearance" (i.e. the column and velocity structure of these absorption features) as a function of the observer inclination angle. This work focuses on just the broadest characteristics of these objects; nonetheless, it provides scaling laws that allow one to reproduce within this model the properties of objects spanning a very wide luminosity range and viewed at different inclination angles, and trace them to a common underlying dynamical structure. Its general conclusion is that the AGN phenomenology can be accounted for in terms of three parameters: The wind mass flux in units of the Eddington value, $\dot m$, the observer's inclination angle $\theta$ and the logarithmic slope between the O/UV and X-ray fluxes $\alpha_{OX}$. However, because of a significant correlation between $\alpha_{OX}$ and UV luminosity, we conclude that the AGN structure depends on only two parameters. Interestingly, the correlations implied by this model appear to extend to and consistent with the characteristics of galactic X-ray sources, suggesting the presence of a truly unified underlying structure for accretion powered sources.

The Poynting-Robertson Cosmic Battery proposes that the innermost part of the accretion disk around a black hole is threaded by a large scale dipolar magnetic field generated in situ, and that the return part of the field diffuses outward through the accretion disk. This is different from the scenario that the field originates at large distances and is carried inward by the accretion flow. In view of the importance of large scale magnetic fields in regulating the processes of accretion and outflows, we study the stability of the inner edge of a magnetized disk in general relativity when the distribution of the magnetic field is the one predicted by the Poynting-Robertson Cosmic Battery. We found that as the field grows, the inner edge of the disk gradually moves outward. In a fast spinning black hole with a>0.8M the inner edge moves back in towards the black hole horizon as the field grows beyond some threshold value. In all cases, the inner part of the disk undergoes a dramatic structural change as the field approaches equipartition.

We investigate the shapes of \gamma-ray pulsar light curves using 3D pulsar magnetosphere models of finite conductivity. These models, covering the entire spectrum of solutions between vacuum and force-free magnetospheres, for the first time afford mapping the GeV emission of more realistic, dissipative pulsar magnetospheres. To this end we generate model light curves following two different approaches: (a) We employ the emission patterns of the slot and outer gap models in the field geometries of magnetospheres with different conductivity \sigma. (b) We define realistic trajectories of radiating particles in magnetospheres of different \sigma and compute their Lorentz factor under the influence of magnetospheric electric fields and curvature radiation-reaction; with these at hand we then calculate the emitted radiation intensity. The light curves resulting from these prescriptions are quite sensitive to the value of \sigma, especially in the second approach. While still not self-consistent, these results are a step forward in understanding the physics of pulsar \gamma-radiation.

We present the magnetic and electric field structures as well as the currents and charge densities of pulsar magnetospheres which do not obey the ideal condition, ${\bf E \cdot B =0}$. Since the acceleration of particles and the production of radiation requires the presence of an electric field component parallel to the magnetic field, ${\bf E}_\parallel$, the structure of non-Ideal pulsar magnetospheres is intimately related to the production of pulsar radiation. Therefore, knowledge of the structure of non-Ideal pulsar magnetospheres is important because their comparison (including models for the production of radiation) with observations will delineate the physics and the parameters underlying the pulsar radiation problem. We implement a variety of prescriptions that support nonzero values for ${\bf E}_\parallel$ and explore their effects on the structure of the resulting magnetospheres. We produce families of solutions that span the entire range between the vacuum and the (ideal) Force-Free Electrodynamic solutions. We also compute the amount of dissipation as a fraction of the Poynting flux for pulsars of different angles between the rotation and magnetic axes and conclude that this is at most 20-40% (depending on the non-ideal prescription) in the aligned rotator and 10% in the perpendicular one. We present also the limiting solutions with the property $J=\rho c$ and discuss their possible implication on the determination of the "on/off" states of the intermittent pulsars. Finally, we find that solutions with values of $J$ greater than those needed to null ${\bf E}_\parallel$ locally produce oscillations, potentially observable in the data.

Neutron-star and black-hole X-ray binaries (XRBs) exhibit radio jets, whose properties depend on the X-ray spectral state and history of the source. In particular, black-hole XRBs emit compact, steady radio jets when they are in the so-called hard state, the jets become eruptive as the sources move toward the soft state, disappear in the soft state, and re-appear when the sources return to the hard state. On the other hand, jets from neutron-star X-ray binaries are typically weaker radio emitters than the black-hole ones at the same X-ray luminosity and in some cases radio emission is detected in the soft state. Significant phenomenology has been accumulated so far regarding the spectral states of neutron-star and black-hole XRBs, and there is general agreement about the type of the accretion disk around the compact object in the various spectral states. Our aim is to investigate whether the phenomenology regarding the X-ray emission on one hand and the jet appearance and disappearance on the other can be put together in a consistent physical picture. It has been shown that the so-called Poynting-Robertson Cosmic Battery (PRCB) explains in a natural way the formation of magnetic fields in the disks of AGN and the ejection of jets. We investigate whether the PRCB can also explain the formation, destruction, and variability of jets in XRBs. We find excellent agreement between the conditions under which the PRCB is efficient (i.e., the type of the accretion disk) and the emission or destruction of the radio jet. The disk-jet connection in XRBs is explained in a natural way using the PRCB.

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light curves like Vela, with vacuum dipole fits being more favorable.

The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems. The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA 1 pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models, within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle, maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius, under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry.

We present the structure of the 3D ideal MHD pulsar magnetosphere to a radius ten times that of the light cylinder, a distance about an order of magnitude larger than any previous such numerical treatment. Its overall structure exhibits a stable, smooth, well-defined undulating current sheet which approaches the kinematic split monopole solution of Bogovalov 1999 only after a careful introduction of diffusivity even in the highest resolution simulations. It also exhibits an intriguing spiral region at the crossing of two zero charge surfaces on the current sheet, which shows a destabilizing behavior more prominent in higher resolution simulations. We discuss the possibility that this region is physically (and not numerically) unstable. Finally, we present the spiral pulsar antenna radiation pattern.

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.

We extend our modeling of the ionization structure of magnetohydrodynamic (MHD) accretion-disk winds, previously applied to Seyfert galaxies, to a population of quasi-stellar-objects (QSOs) of much lower X-ray-to-UV flux ratios, i.e. smaller $\alpha_{\rm ox}$ index, motivated by UV/X-ray ionized absorbers with extremely high outflow velocities in UV-luminous QSOs. We demonstrate that magnetically-driven winds ionized by a spectrum with $\alpha_{\rm ox} \simeq -2$ can produce the charge states responsible for \civ ~and \fexxv/\fexxvi ~absorption in wind regions with corresponding maximum velocities of $v$(\civ) $\lsim 0.1c$ and $v({\rm \fexxv}) \lsim 0.6 c$ (where $c$ is the speed of light) and column densities $N_H \sim 10^{23}-10^{24}$ cm$^{-2}$, in general agreement with observations. In contrast to the conventional radiation-driven wind models, {\it high-velocity flows are always present in our MHD-driven winds} but manifest in the absorption spectra only for $\alpha_{\rm ox} \lsim -2$, as larger $\alpha_{\rm ox}$ values ionize the wind completely out to radii too large to demonstrate the presence of these high velocities. We thus predict increasing velocities of these ionized absorbers with decreasing (steeper) $\alpha_{\rm ox}$, a quantity that emerges as the defining parameter in the kinematics of the AGN UV/X-ray absorbers.

We present the 2D ionization structure of self-similar MHD winds off accretion disks around irradiated by a central X-ray source. Based on earlier observational clues and theoretical arguments, we focus our attention on a subset of these winds, namely those with radial density dependence n(r)~1/r. We employ the photoionization code XSTAR to compute the ionic abundances of a large number of ions of different elements and then compile their line-of-sight (LOS) absorption columns. Particular attention is paid to the absorption measure distribution (AMD), namely their Hydrogen-equivalent column per logarithmic ionization parameter \xi interval, d N_H/(d \log \xi), which provides a measure of the winds' radial density profiles. For n(r)~1/r the AMD is found to be independent of \xi, in good agreement with its behavior inferred from the X-ray spectra of several AGNs. For the specific wind structure and X-ray spectrum we also compute detailed absorption line profiles for a number of ions to obtain their LOS velocities, v~100-300 km/sec (at \log \xi~2-3) for Fe XVII and v~1,000-4,000 km/sec (at \log \xi~4-5) for Fe XXV, in good agreement with the observation. Our models describe the X-ray absorption properties of these winds with only two parameters, namely the mass-accretion rate \dot{m} and LOS angle \theta. The probability of obscuration of the X-ray ionizing source in these winds decreases with increasing \dot{m} and increases steeply with \theta. As such, we concur with previous authors that these wind configurations, viewed globally, incorporate all the requisite properties of the parsec scale "torii" invoked in AGN unification schemes. We indicate that a combination of the AMD and absorption line profile observations can uniquely determine these model parameters and their bearing on AGN population demographics.

We investigate the strong electric current sheet that develops at the tip of the pulsar closed line region through time dependent three-dimensional numerical simulations of a rotating magnetic dipole. We show that curvature radiation from relativistic electrons and positrons in the current sheet may naturally account for several features of the high-energy pulsar emission. We obtain light curves and polarization profiles for the complete range of magnetic field inclination angles and observer orientations, and compare our results to recent observations from the Fermi gamma-ray telescope.

The origin of cosmic magnetic (B) fields remains an open question. It is generally believed that very weak primordial B fields are amplified by dynamo processes, but it appears unlikely that the amplification proceeds fast enough to account for the fields presently observed in galaxies and galaxy clusters. In an alternative scenario, cosmic B fields are generated near the inner edges of accretion disks in Active Galactic Nuclei (AGNs) by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with photons. While dynamo processes show no preference for the polarity of the (presumably random) seed field that they amplify, this alternative mechanism uniquely relates the polarity of the poloidal B field to the angular velocity of the accretion disk, resulting in a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry came about by chance being less than 1%. This lends support to the hypothesis that the Universe is seeded by B fields that are generated in AGN via this mechanism and subsequently injected into intergalactic space by the jet outflows.

We propose a simple physical picture for the generation of coherent radio emission in the axisymmetric pulsar magnetosphere that is quite different from the canonical paradigm of radio emission coming from the magnetic polar caps. In this first paper we consider only the axisymmetric case of an aligned rotator. Our picture capitalizes on an important element of the MHD representation of the magnetosphere, namely the separatrix between the corotating closed-line region (the `dead zone') and the open field lines that originate in the polar caps. Along the separatrix flows the return current that corresponds to the main magnetospheric electric current emanating from the polar caps. Across the separatrix, both the toroidal and poloidal components of the magnetic field change discontinuously. The poloidal component discontinuity requires the presence of a significant annular electric current which has up to now been unaccounted for. We estimate the position and thickness of this annular current at the tip of the closed line region, and show that it consists of electrons (positrons) corotating with Lorentz factors on the order of 10^5, emitting incoherent synchrotron radiation that peaks in the hard X-rays. These particles stay in the region of highest annular current close to the equator for a path-length of the order of one meter. We propose that, at wavelengths comparable to that path-length, the particles emit coherent radiation, with radiated power proportional to N^2, where N is the population of particles in the above path-length. We calculate the total radio power in this wavelength regime and its scaling with pulsar period and stellar magnetic field and show that it is consistent with estimates of radio luminosity based on observations.

We investigate the three-dimensional structure of the pulsar magnetosphere through time-dependent numerical simulations of a magnetic dipole that is set in rotation. We developed our own Eulerian finite difference time domain numerical solver of force-free electrodynamics and implemented the technique of non-reflecting and absorbing outer boundaries. This allows us to run our simulations for many stellar rotations, and thus claim with confidence that we have reached a steady state. A quasi-stationary corotating pattern is established, in agreement with previous numerical solutions. We discuss the prospects of our code for future high-resolution investigations of dissipation, particle acceleration, and temporal variability.

Positive and negative pulsar breaking indices suggest that some fraction of the pulsar spindown torque undergoes a cyclic evolution. The observed strong correlation of `anomalous' breaking indices with pulsar age implies that the characteristic periodicity timescale is in the range 100 to 10,000 years depending on the fraction of the spindown torque that undergoes cyclic evolution, 1 to 100% respectively. We argue that the longest variability timescale is consistent with a neutron star magnetic cycle similar to the solar cycle.

We show that toroidal magnetic field annihilation in the equatorial current sheet of the pulsar magnetosphere is related to how fast poloidal magnetic field lines close as we move away from the light cylinder. This determines the reconnection radial electric field which directly accelerates particles in the neutral layer inside the equatorial current sheet. The efficiency of poloidal magnetic field closure near the light cylinder may be measurable through the pulsar braking index. We argue that, the lower the efficiency of pair formation, the higher the braking index. We also argue that synchrotron radiation reaction in the neutral layer does not inhibit the accelerated particles from reaching the maximum energy of about 10^16 eV available in the open pulsar magnetosphere.

We describe the results of numerical "2.5--dimensional" MHD simulations of an initially unmagnetized disk model orbiting a central point--mass and responding to the continual generation of poloidal magnetic field due to a secular source that emulates the Poynting--Robertson (PR) drag on electrons in the vicinity of a luminous stellar or compact accreting object. The fluid in the disk and in the surrounding hotter atmosphere has finite electrical conductivity and allows for the magnetic field to diffuse freely out of the areas where it is generated, while at the same time, the differential rotation of the disk twists the poloidal field and quickly induces a substantial toroidal--field component. The secular PR term has dual purpose in these simulations as the source of the magnetic field and the trigger of a magnetorotational instability (MRI) in the disk. The MRI is especially mild and does not destroy the disk because a small amount of resistivity dampens the instability efficiently. In simulations with moderate resistivities (diffusion timescales up to $\sim$16 local dynamical times) and after $\sim$100 orbits, the MRI has managed to transfer outward substantial amounts of angular momentum and the inner edge of the disk, along with azimuthal magnetic flux, has flowed toward the central point--mass where a new, magnetized, nuclear disk has formed. The toroidal field in this nuclear disk is amplified by differential rotation and it cannot be contained; when it approaches equipartition, it unwinds vertically and produces episodic jet--like outflows. The poloidal field in the inner region cannot diffuse back out if the characteristic diffusion time is of the order of or larger than the dynamical time; it continues to grow linearly in time undisturbed and without saturation, as the outer sections of many poloidal loops are being drawn radially outward.

We propose that cosmological magnetic fields generated in regions of finite spatial dimensions may manifest themselves in the global dynamics of the Universe as `dark energy'. We test our model in the context of spatially flat cosmological models by assuming that the Universe contains non-relativistic matter $\rho_m\propto \alpha^{-3}$, dark energy $\rho_{Q}\propto \alpha^{-3(1+w)}$, and an extra fluid with $\rho_{B} \propto \alpha^{n-3}$ that corresponds to the magnetic field. We place constraints on the main cosmological parameters of our model by combining the recent supernovae type Ia data and the differential ages of passively evolving galaxies. In particular, we find that the model which best reproduces the observational data when $\Omega_m=0.26$ is one with $\Omega_{B}\simeq 0.03$, $n\simeq 7.68$, $\Omega_{Q}\simeq 0.71$ and $w\simeq -0.8$.

We evaluate the result of the recent pioneering numerical simulations in Spitkovsky~2006 on the spindown of an oblique relativistic magnetic dipole rotator. Our discussion is based on our experience from two idealized cases, that of an aligned dipole rotator, and that of an oblique split-monopole rotator. We conclude that the issue of electromagnetic pulsar spindown may not have been resolved yet.

The present work is our first attempt to understand the role of reconnection in the pulsar magnetosphere. Our discussion is based on the observationaly infered fact that, as the pulsar spins down, the region of closed corotating dipolar field lines grows with time. This implies that reconnection must take place in the magnetosphere. We argue that non-dissipative reconnection along the equatorial current sheet allows for the continuous channeling of pulsar spindown energy into particle energy, all the way from the light cylinder to the pulsar wind termination shock, and we propose that this effect may account for the low \sigma values inferred by observations. We present a simple model that allows us to relate the magnetic diffusivity in the equatorial current sheet to an observable pulsar parameter, the braking index n. When n~1, the global structure of the magnetosphere approaches that of a relativistic split monopole where the pulsar spindown energy is carried by the electromagnetic field. However, for values of n>1.5, almost all field lines close inside the pulsar wind termination shock, and thus most of the electromagnetic pulsar spindown energy flux is effectively transformed into particle energy in the equatorial current sheet.

This paper has been withdrawn by the author. It will be posted on astro-ph once the proceedings of the 363rd Heraeus Seminar on Neutron Stars and Pulsars are ready.

We address the issue of electromagnetic pulsar spindown by combining our experience from the two limiting idealized cases which have been studied in great extent in the past: that of an aligned rotator where ideal MHD conditions apply, and that of a misaligned rotator in vacuum. We construct a spindown formula that takes into account the misalignment of the magnetic and rotation axes, and the magnetospheric particle acceleration gaps. We show that near the death line aligned rotators spin down much slower than orthogonal ones. In order to test this approach, we use a simple Monte Carlo method to simulate the evolution of pulsars and find a good fit to the observed pulsar distribution in the P-Pdot diagram without invoking magnetic field decay. Our model may also account for individual pulsars spinning down with braking index n < 3, by allowing the corotating part of the magnetosphere to end inside the light cylinder. We discuss the role of magnetic reconnection in determining the pulsar braking index. We show, however, that n ~ 3 remains a good approximation for the pulsar population as a whole. Moreover, we predict that pulsars near the death line have braking index values n > 3, and that the older pulsar population has preferentially smaller magnetic inclination angles. We discuss possible signatures of such alignment in the existing pulsar data.

Polar magnetospheric gaps consume a fraction of the electric potential that develops accross open field lines. This effect modifies significantly the structure of the axisymmetric pulsar magnetosphere. We present numerical stead-state solutions for various values of the gap potential. We show that a charge starved magnetosphere contains significantly less electric current than one with freely available electric charges. As a result, electromagnetic neutron star braking becomes inefficient. We argue that the magnetosphere may spontaneously rearrange itself to a lower energy configuration through a dramatic release of electromagnetic field energy and magentic flux. Our results might be relevant in understanding the recent December 27, 2004 burst observed in SGR 1806-20.

Using the exact solution of the axisymmetric pulsar magnetosphere derived in a previous publication and the conservation laws of the associated MHD flow, we show that the Lorentz factor of the outflowing plasma increases linearly with distance from the light cylinder. Therefore, the ratio of the Poynting to particle energy flux, generically referred to as $\sigma$, decreases inversely proportional to distance, from a large value (typically $\gsim 10^4$)near the light cylinder to $\s \simeq 1$ at a transistion distance $R_{\rm trans}$. Beyond this distance the inertial effects of the outflowing plasma become important and the magnetic field geometry must deviate from the almost monopolar form it attains between $R_{lc}$ and $R_{\rm trans}$. We anticipate that this is achieved by collimation of the poloidal field lines toward the rotation axis, ensuring that the magnetic field pressure in the equatorial region will fall-off faster than $1/R^2$ ($R$ being the cylindrical radius). This leads both to a value $\s=\ss \ll 1$ at the nebular reverse shock at distance $R_s$ ($R_s \gg R_{\rm trans}$) and to a component of the flow perpendicular to the equatorial component, as required by observation. The presence of the strong shock at $R = R_s$ allows for the efficient conversion of kinetic energy into radiation. We speculate that the Crab pulsar is unique in requiring $\ss \simeq 3 \times 10^{-3}$ because of its small translational velocity, which allowed for the shock distance $R_s$ to grow to values $\gg R_{\rm trans}$.

We examine the effects of a relativistic wind on the spin down of a neutron star and apply our results to the study of Soft Gamma Repeaters (SGRs), thought to be neutron stars with magnetic fields > 10^{14} G. We derive a spin-down formula that includes torques from both dipole radiation and episodic or continuous particle winds. We find that if SGR1806-20 puts out a continuous particle wind of 10^{37} {\rm erg} {\rm s}^{-1}, then the pulsar age is consistent with that of the surrounding supernova remnant, but the derived surface dipole magnetic field is only 3 \times 10^{13} G, in the range of normal radio pulsars. If instead, the particle wind flows are episodic with small duty cycle, then the observed period derivatives imply magnetar-strength fields, while still allowing characteristic ages within a factor of two of the estimated supernova remnant age. Close monitoring of the periods of SGRs will be able to establish or place limits on the wind duty cycle and thus the magnetic field and age of the neutron star.

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 show that the Poynting-Robertson drag effect in an optically thin advection-dominated accretion flow around active gravitating objects generates strong azimuthal electric currents which give rise to astrophysically significant magnetic fields. Although the mechanism is most effective in accreting compact objects, it seems very promising to also account for the generation of stellar dipolar fields during the late protostellar collapse phase, when the star approaches the main sequence.

We obtain self-similar solutions that describe the gravitational collapse of nonrotating, isothermal, magnetic molecular cloud cores. We use simplifying assumptions but explicitly include the induction equation, and the semianalytic solutions we derive are the first to account for the effects of ambipolar diffusion following the formation of a central point mass. Our results demonstrate that, after the protostar first forms, ambipolar diffusion causes the magnetic flux to decouple in a growing region around the center. The decoupled field lines remain approximately stationary and drive a hydromagnetic C-shock that moves outward at a fraction of the speed of sound (typically a few tenths of a kilometer per second), reaching a distance of a few thousand AU at the end of the main accretion phase for a solar-mass star. We also show that, in the absence of field diffusivity, a contracting core will not give rise to a shock if, as is likely to be the case, the inflow speed near the origin is nonzero at the time of point-mass formation. Although the evolution of realistic molecular cloud cores will not be exactly self-similar, our results reproduce the main qualitative features found in detailed core-collapse simulations (Ciolek & Konigl 1998)