Search

All numerical solutions of the pulsar magnetosphere over the past 25 years show closed-line regions that end a significant distance inside the light cylinder, and manifest thick strongly dissipative separatrix surfaces instead of thin current sheets, with a tip that has a distinct pointed Y shape instead of a T shape. We need to understand the origin of these results which were not predicted by our early theories of the pulsar magnetosphere. In order to gain new intuition on this problem, we set out to obtain the theoretical steady-state solution of the 3D ideal force-free magnetosphere with zero dissipation along the separatrix and equatorial current sheets. In order to achieve our goal, we needed to develop a novel numerical method. We solve two independent non-singular magnetospheric problems in the domains of open and closed field lines, and adjust the shape of their interface (the separatrix) to satisfy pressure balance between the two regions. The solution is obtained with meshless Physics Inspired Neural Networks (PINNs) In this Letter we present preliminary results for an inclined dipole rotator using the new methodology. This is the first time the steady-state 3D problem is addressed directly, and not through a time-dependent simulation that eventually relaxes to a steady-state. Our results demonstrate the potential of the new method to generate the reference solution of the pulsar magnetosphere., Comment: 6 pages, 3 figures, submitted for publication in Astronomy & Astrophysics

The particle-in-cell approach has proven effective at modeling neutron star and black hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where microphysics operates and the system-size scales due to limited numerical resources. A legitimate concern is whether the scale separation currently achieved is large enough, such that results can be safely extrapolated to realistic scales. In this work, our aim is to explore the effect of scaling physical parameters up, and to check whether salient features uncovered by pure kinetic models at smaller scales are still valid, with a special emphasis on particle acceleration and high-energy radiation emitted beyond the light cylinder. To reach this objective, we develop a new hybrid numerical scheme coupling the ideal force-free and the particle-in-cell methods, to optimize the numerical cost of global models. We propose a domain decomposition of the magnetosphere based on the magnetic field topology using the flux function. The force-free model is enforced along open field lines while the particle-in-cell model is restricted to the reconnecting field line region. As a proof of concept, this new hybrid model is applied to simulate a weak millisecond pulsar magnetosphere with realistic scales using high-resolution axisymmetric simulations. Magnetospheric features reported by previous kinetic models are recovered, and strong synchrotron radiation above 100MeV consistent with the Fermi-LAT gamma-ray pulsar population is successfully reproduced. This work further consolidates the shining reconnecting current sheet scenario as the origin of the gamma-ray emission in pulsars, as well as firmly establishes pulsar magnetospheres as at least TeV particle accelerators., Comment: 12 pages, 7 figures, accepted for publication in A&A

We present the first new type of solution of the pulsar equation since 1999. In it, the whole magnetosphere is confined inside the light cylinder and an electrically charged layer wraps around it and holds it together. The reason this new solution has never been obtained before is that all current time-dependent simulations are initialized with a vacuum dipole configuration that extends to infinity, thus their final steady-state solution also extends to infinity. Under special conditions, such a confined configuration may be attained when the neutron star first forms in the interior of a collapsing star during a supernova explosion, or when it accretes from an external wind or disk from a donor star. It is shown that this new maximally closed non-decelerating solution is the limit of a continuous sequence of standard magnetospheres with open and closed field lines when the amount of open field lines gradually drops to zero. The minimum energy solution in this sequence is a standard magnetosphere in which the closed field line region extends up to about 80% of the light cylinder. We estimate that the released energy when the new solution transitions to the minimum energy one is enough to power a fast radio burst., Comment: 15 pages, 6 figures, in press in Journal Universe Special Volume (MDPI) "A New Horizon of Pulsar and Neutron Star: The 55-Year Anniversary" (Eds. C. Zhang, D. Blaschke, J. Petri)

We investigate the generation of gravitational waves from the rotation of an orthogonal pulsar magnetosphere in flat space time. We calculate the first order metric perturbation due to the rotation of the non-axisymmetric distribution of electromagnetic energy density around the central star. We show that gravitational waves from a strong magnetic field pulsar right after its formation within a distance of 1 kpc may be detectable with the new generation of gravitational wave detectors., Comment: 9 pages, 1 figure, accepted for publication in Monthly Notices

In this study, we introduce a novel approach for deriving the solution of the ideal force-free steady-state pulsar magnetosphere in three dimensions. Our method involves partitioning the magnetosphere into the regions of closed and open field lines, and subsequently training two custom Physics Informed Neural Networks (PINNs) to generate the solution within each region. We periodically modify the shape of the boundary separating the two regions (the separatrix) to ensure pressure balance throughout. Our approach provides an effective way to handle mathematical contact discontinuities in Force-Free Electrodynamics (FFE). We present preliminary results in axisymmetry, which underscore the significant potential of our method. Finally, we discuss the challenges and limitations encountered while working with Neural Networks, thus providing valuable insights from our experience., Comment: 12 pages, 7 figures, accepted for publication in Monthly Notices

We develop a classical picture of interference for non-interacting individual classical massive free particles. As long as they remain undetected, particles carry the information of a phase equal to an action integral along their trajectory. At the point of their detection, a special type of detector collects the phases from all individual particles reaching it, adds them up over time as complex numbers, and divides them by the square root of their number. The detector announces a number of detections equal to the square of the amplitude of the resulting complex number. An interference pattern is gradually built from the collection of particle phases in the detection bins of the detector after several repetitions of the experiment. We obtain perfect agreement with three solutions of the Schr\"odinger equation for free particles: a Gaussian wavepacket, two Gaussian wavepackets approaching each other, and a Gaussian wavepacket reflecting off a wall., Comment: 6 figures. arXiv admin note: text overlap with arXiv:2203.10942 Accepted for publication in Particles (MDPI)

We present a statistical model of non-interacting individual classical particles that may lead to a microscopic implementation of quantum mechanics. The model requires the action of a special type of detector that detects and records individual particles. According to our model, microscopic particles are classical entities that follow classical trajectories without nonlocal interactions between them. As long as they remain undetected, particles carry the information of a phase equal to an action integral along their trajectory. At the point of their detection, our special type of detector collects the phases from all particles reaching it, adds them up over time as complex numbers, and divides them by the square root of their number. The detector records a detection each time the square of the magnitude of the ensuing vector crosses an integer value. We test our model with numerical simulations of a double-slit experiment, and obtain an interference pattern analogous to the one observed in an actual physical experiment., Comment: 8 pages, 3 figures, under consideration in Foundations of Physics

We investigate long Gamma-Ray Bursts (GRB) which manifest a sharp linear rise followed by an exponential decay in their gamma-ray prompt emission observed with the BAT instrument on board the Swift satellite. We offer a simple electrodynamic model that may account for these particular characteristics. We associate the sharp rise with the winding of the magnetic field by the fast rotating core that formed in the interior of the stellar precursor. We also associate the subsequent exponential decay with the electromagnetic spin-down of the core following the release of the electromagnetic jet from the stellar interior. Any non-axisymmetric distortion in the rotating core will generate gravitational waves with exponentially decreasing frequency, a so-called "down-chirp". We obtain a detailed estimate of the gravitational wave profile if the distortion of spacetime is due to the winding of a non-axisymmetric component of the magnetic field during that particular phase of the burst. We offer 7 particular time intervals during which one may look into LIGO archival data for the presence of our particular predicted waveforms in order to test our interpretation., Comment: 7 pages, 6 figures, accepted for publication in Monthly Notices

Magnetic fields in galaxy halos are in general very difficult to observe. Most recently, the CHANG-ES collaboration (Continuum HAlos in Nearby Galaxies - an EVLA Survey) investigated in detail the radio halos of 35 nearby edge-on spiral galaxies and detected large scale magnetic fields in 16 of them. We used the CHANG-ES radio polarization data to create Rotation Measure (RM) maps for all galaxies in the sample and stack them with the aim to amplify any underlying universal toroidal magnetic field pattern in the halo above and below the disk of the galaxy. We discovered a large-scale magnetic field in the central region of the stacked galaxy profile, attributable to an axial electric current that universally outflows from the center both above and below the plane of the disk. A similar symmetry-breaking has also been observed in astrophysical jets but never before in galaxy halos. This is an indication that galaxy halo magnetic fields are probably not generated by pure ideal magnetohydrodynamic (MHD) processes in the central regions of galaxies. One such promising physical mechanism is the Cosmic Battery operating in the innermost accretion disk around the central supermassive black hole. We anticipate that our discovery will stimulate a more general discussion on the origin of astrophysical magnetic fields., Comment: 10 pages, 6 figures, accepted for publication in the 4. Extragalactic astronomy section of Astronomy & Astrophysics

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., Comment: 7 pages, 5 figures, Monthly Notices in press

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., Comment: 9 pages, 6 figures, Monthly Notices

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., Comment: 9 pages, 2 figures, International Conference on Black Holes as Cosmic Batteries: UHECRs and Multimessenger Astronomy - BHCB2018, 12-15 September 2018, Foz du Iguazu, Brazil. arXiv admin note: text overlap with arXiv:1901.01404

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., Comment: 18 pages, 4 figures, accepted for publication in the Special Issue "Cosmic Plasmas and Electromagnetic Phenomena" of Journal Galaxies

We investigate the development of the magnetic Rayleigh-Taylor instability at the inner edge of an astrophysical disk around a spinning central black hole. We solve the equations of general relativity that govern small amplitude oscillations of a discontinuous interface in a Keplerian disk threaded by an ordered magnetic field, and we derive a stability criterion that depends on the central black hole spin and the accumulated magnetic field. We also compare our results with the results of GR MHD simulations of black hole accretion flows that reach a magnetically arrested state (MAD). We found that the instability growth timescales that correspond to the simulation parameters are comparable to the corresponding timescales for free-fall accretion from the ISCO onto the black hole. We thus propose that the Rayleigh-Taylor instability disrupts the accumulation of magnetic flux onto the black hole horizon as the disk reaches a MAD state., Comment: 11 pages, 1 figure, accepted for publication in Monthly Notices

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., Comment: 5 pages, 2 figures, accepted for publication in Monthly Notices Letters

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., Comment: accepted to ApJ Letter; 16 pages, 3 figures

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., Comment: 37 pages, 7 figs, 3 tables; accepted to ApJ

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., Comment: 4 pages, 1 figure, Monthly Notices Letters in press

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., Comment: 7 pages, 4 figures, submitted to Monthly Notices

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., Comment: 5 pages, 1 figures, accepted for publication in Monthly Notices Letters

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., Comment: published in 2017 March 6 Nature Astronomy, 23 pages, 4 figures, 4 supplementary figures

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., Comment: 10 pages, accepted for publication in Monthly Notices

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., Comment: 5 pages, 1 figure, Monthly Notices Letters

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., Comment: 12 pages, 4 figures, 1 table, accepted for publication in Astronomy & Astrophysics

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., Comment: 13 pages, 2 figures, accepted for publication in the Journal of Plasma Physics

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., Comment: 8 pages, 6 figures

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., Comment: 5 pages, 4 figures, MNRAS accepted

We obtain an approximate solution for the motion of a charged particle around a Schwarzschild black hole immersed in a weak dipolar magnetic field. We focus on eccentric bound orbits in the equatorial plane of the Schwarzschild black hole and derive an analytic expression for the spectral distribution of the electromagnetic emission from a charged particle on such an orbit. Two sets of harmonic contributions appear, with specific frequency spacing. The expression can be written in compact form, if it is truncated up to the lowest order harmonic contributions., Comment: 13 pages, Journal of Mathematical Physics

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., Comment: v.3 as of 5/6/15 with eliminating extra figs: accepted to ApJ, 28 pages, figs.1-6 (color), 3 tables

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., Comment: 11 pages, 5 figures, acceped for publication in the Astrophysical Journal

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., Comment: 6 pages, 2 figures, presented in "Swift: 10 years of discovery" (Rome 2014)

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., Comment: 14 pages, 5 figures, accepted for publication in the Astrophysical Journal

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., Comment: 10 pages, 4 figures, recommended for publication in the Astrophysical Journal

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., Comment: version 2 (modified), 27 pages, 5 figures, accepted to ApJ

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

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., Comment: 5 pages, 3 figures, accepted for publication in ApJ

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'., Comment: 7 pages, 2 figures, accepted for publication in the Astrophysical Journal

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., Comment: 9 pages, 1 figure, accepted for publication in Solar Physics

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., Comment: submitted to the Astronomical Review, 32pg, 8 figs

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., Comment: 7 pages, 3 figures, accepted for publication in Monthly Notices of the RAS

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., Comment: Submitted to the Astrophysical Journal Letters, 15 pages, 4 figures

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., Comment: 6 pages, 4 figures, 2011 Fermi Symposium proceedings - eConf C110509

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., Comment: 2011 Fermi Symposium proceedings - eConf C110509, 4 pages, 3 figures

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., Comment: 6 pages, 3 figures, accepted in MNRAS