Force-free Electrodynamics Research Articles

Novel scheme for simulating the force-free equations: Boundary conditions and the evolution of solutions towards stationarity

Force-Free Electrodynamics (FFE) describes a particular regime of magnetically dominated relativistic plasmas, which arises on several astrophysical scenarios of interest such as pulsars or active galactic nuclei. In this article, we present a full 3D numerical implementation of the FFE evolution around a Kerr black hole. The novelty of our approach is three-folded: i) We use the "multi-block" technique to represent a domain with $S^2 \times \mathbb{R}^{+}$ topology within a stable finite-differences scheme. ii) We employ as evolution equations those arising from a covariant hyperbolization of the FFE system. iii) We implement stable and constraint-preserving boundary conditions to represent an outer region given by a uniform magnetic field aligned or misaligned respect to the symmetry axis. We find stationary jet solutions which reach equilibrium --through boundary conditions-- with the outer numerical surface. This is so, even when the outer boundary is located very close to the central region (i.e. $r_{out}\sim 10M $). These numerical solutions reproduce most of the known results for analogue astrophysical settings.

Expanded solutions of force-free electrodynamics on general Kerr black holes

In this work, expanded solutions of force-free magnetospheres on general Kerr black holes are derived through a radial distance expansion method. From the regular conditions both at the horizon and at spatial infinity, two previously known asymptotical solutions (one of them is actually an exact solution) are identified as the only solutions that satisfy the same conditions at the two boundaries. Taking them as initial conditions at the boundaries, expanded solutions up to the first few orders are derived by solving the stream equation order by order. It is shown that our extension of the exact solution can (partially) cure the problems of the solution: it leads to magnetic domination and a mostly timelike current for restricted parameters.

Simulations of the magnetospheres of accreting millisecond pulsars

Accreting pulsars power relativistic jets and display a complex spin phenomenology. These behaviours may be closely related to the large-scale configuration of the star's magnetic field, shaped by its interaction with the surrounding accretion disc. Here, we present the first relativistic simulations of the interaction of a pulsar magnetosphere with an accretion flow. Our axisymmetric simulations treat the magnetospheric, or coronal, regions using a resistive extension of force-free electrodynamics. The magnetic field is also evolved inside the disc, which is a defined volume with a specified velocity field and conductivity profile, found using an α-disc model. We study a range of disc α-parameters, thicknesses, magnetic Prandtl numbers and inner truncation radii. We find that a large fraction of the magnetic flux in the pulsar's closed zone is opened by the intrusion of the disc, leading to an enhancement of the power extracted by the pulsar wind and the spin-down torque applied to the pulsar. In our simulations, most of the spin-down contribution to the stellar torque acts on open field lines. The efficiency of field-line opening is high in the simulations’ long-term quasi-steady states, which implies that a millisecond pulsar's electromagnetic wind could be strong enough to power the observed neutron-star radio jets, and may significantly affect the pulsar's spin evolution.

Metric Independence of Vacuum and Force-Free Electromagnetic Fields.

Electromagnetic fields which solve the vacuum Maxwell equations in one spacetime are well known to also be solutions in all spacetimes with conformally related metrics. This provides a sense in which electromagnetism alone cannot be used to measure certain aspects of geometry. We show that there is actually much more which cannot be so measured; relatively little of a spacetime's geometry is in fact imprinted in any particular electromagnetic field. This is demonstrated by finding a much larger class of metric transformations-involving five free functions-which preserve Maxwell solutions both in vacuum, without local currents, and also for the force-free electrodynamics associated with a tenuous plasma. One consequence of this is that many of the exact force-free fields which have previously been found around Schwarzschild and Kerr black holes are also solutions in appropriately identified flat backgrounds. As a more direct application, we use our metric transformations to write down a large class of electromagnetic waves which remain unchanged by a large class of gravitational waves propagating "in the same direction."

Ill posedness of force-free electrodynamics in Euler potentials

We prove that the initial value problem for Force-free Electrodynamics in Euler variables is not well posed. We establish this result by showing that a well-posedness criterion provided by Kreiss fails to hold for this theory, and using a theorem provided by Strang. To show the nature of the problem we display a particular bounded (in Sobolev norms) sequence of initial data for the Force-free equations such that at any given time as close to zero as one wishes, the corresponding evolution sequence is not bounded. Thus, the Force-free evolution is non continuous in that norm with respect to the initial data. We furthermore prove that this problem is also ill-posed in the Leray-Ohya sense.

Accumulative coupling between magnetized tenuous plasma and gravitational waves

We explicitly compute the plasma wave (PW) induced by a plane gravitational wave (GW) travelling through a region of strongly magnetized plasma, governed by force-free electrodynamics. The PW co-moves with the GW and absorbs its energy to grow over time, creating an essentially force-free counterpart to the inverse-Gertsenshtein effect. The time-averaged Poynting flux of the induced PW is comparable to the vacuum case, but the associated current may offer a more sensitive alternative to photodetection when designing experiments for detecting/constraining high frequency gravitational waves. Aside from the exact solutions, we also offer an analysis of the general properties of the GW to PW conversion process, which should find use when evaluating electromagnetic counterparts to astrophysical gravitational waves, that are generated directly by the latter as a second order phenomenon.

Force-Free Electromagnetic Fields within Spinor Framework

The article deals with spinor representation of the force-free electrodynamics. The equations of the force-free electromagnetic field describe the physics of pulsars and black holes whose magnetospheres are filled with magnetically dominated relativistic plasma. The paper is a brief pedagogical introduction to the mathematics of the subject, based on 2-spinor calculi. The objective is to present the nonlinear theory of force-free fields in a compact and elegant form that the spinor framework provides. First, the algebraic classification of the Maxwell tensor is presented. Then, the reduced system of differential equations is obtained for two types of electromagnetic field and the basic features of the solutions are described. The null force-free field is connected with the shear-free geodesic null congruence in a space-time and is derived from a linear equation for a complex function. The magnetic force-free field is associated with the time-like 2-surface that represents the world-sheet of magnetic field line. The simplified system includes 4 linear differential equations for a real function. The article is educational in nature and there are no new solutions of force-free equations obtained.

CRAB FLARES DUE TO TURBULENT DISSIPATION OF THE PULSAR STRIPED WIND

ABSTRACT We interpret γ-ray flares from the Crab Nebula as the signature of turbulence in the pulsar’s electromagnetic outflow. Turbulence is triggered upstream by dynamical instability of the wind’s oscillating magnetic field and accelerates non-thermal particles. On impacting the wind-termination shock, these particles emit a distinct synchrotron component , which is constantly modulated by intermittency of the upstream plasma flow. Flares are observed when the high-energy cutoff of emerges above the fast-declining nebular emission around 0.1–1 GeV. Simulations carried out in the force-free electrodynamics approximation predict the striped wind to become fully turbulent well ahead of the wind-termination shock, provided its terminal Lorentz factor is .

QED plasma and magnetars

Magnetars are surrounded by diffuse plasma in magnetic field strengths well above the quantum electrodynamic critical value. We derive equations of ``quantum force-free electrodynamics'' for this plasma using effective field theory arguments. We argue that quantum effects do not modify the large scale structure of the magnetosphere, and in particular that the spin-down rate does not deviate significantly from the classical result. We provide definite evolution equations that can be used to explore potentially important small-scale corrections, such as shock formation, which has been proposed as a mechanism for both burst and quiescent emission from magnetars.

Covariant hyperbolization of force-free electrodynamics

Force-Free Flectrodynamics (FFE) is a non-linear system of equations modeling the evolution of the electromagnetic field, in the presence of a magnetically dominated relativistic plasma. This configuration arises on several astrophysical scenarios, which represent exciting laboratories to understand physics in extreme regimes. We show that this system, when restricted to the correct constraint submanifold, is symmetric hyperbolic. In numerical applications is not feasible to keep the system in that submanifold, and so, it is necessary to analyze its structure first in the tangent space of that submanifold and then in a whole neighborhood of it. As already shown by Pfeiffer, a direct (or naive) formulation of this system (in the whole tangent space) results in a weakly hyperbolic system of evolution equations for which well-possednes for the initial value formulation does not follows. Using the generalized symmetric hyperbolic formalism due to Geroch, we introduce here a covariant hyperbolization for the FFE system. In fact, in analogy to the usual Maxwell case, a complete family of hyperbolizers is found, both for the restricted system on the constraint submanifold as well as for a suitably extended system defined in a whole neighborhood of it. A particular symmetrizer among the family is then used to write down the pertaining evolution equations, in a generic ($3+1$)-decomposition on a background spacetime. Interestingly, it turns out that for a particular choice of the lapse and shift functions of the foliation, our symmetrized system reduces to the one found by Pfeiffer. Finally we analyze the characteristic structure of the resulting evolution system.

Hertz potential formalism for force-free electrodynamics and its application to Brennan–Gralla–Jacobson solutions

The Hertz potential is a powerful tool for the source-free electrodynamics. Especially for the algebraically special spacetime background, the Hertz potential formalism simplifies the Maxwell equations quite much. In astrophysics, strong electric–magnetic field is very common. Force-free electrodynamics is a good approximation for strong enough electric–magnetic field compared to the inertial energy of the involved plasma. For example, the force-free model has been extensively used to describe the magnetosphere of stars in the universe. In this paper, we extend the Hertz potential formalism to the force-free electrodynamics. The Hertz potential formalism simplifies the force-free dynamical equations as much as that in the source-free case. As an application, we use the Hertz potential formalism to the Schwarzschild background. And the Brennan–Gralla–Jacobson solutions are recovered straightforwardly.

Near-horizon extreme Kerr magnetospheres

Analytical solutions to force-free electrodynamics around black holes are fundamental for building simple models of accretion disk and jet dynamics. We present a (nonexhaustive) classification of complex highest-weight solutions to the force-free equations in the near-horizon region of the extremal Kerr black hole. Several classes of real magnetically dominated or null solutions, either axisymmetric or nonaxisymmetric, are described which admit finite energy and angular momentum with respect to the asymptotically flat observer. Subtleties related to the velocity of light surface in the near-horizon region are discussed.

Structure of Aristotelian electrodynamics

Aristotelian electrodynamics (AE) describes the regime of a plasma with a very strong electric field that is not shorted out, with the charge current determined completely by pair production and the balance of the Lorentz 4-force against the curvature radiation reaction. Here it is shown how the principal null directions and associated eigenvalues of the field tensor govern AE, and how force-free electrodynamics arises smoothly from AE when the eigenvalues (and therefore the electric field in some frame) vanish. A criterion for validity of AE and force-free electrodynamics is proposed in terms of a pair of ``field curvature scalars'' formed from the first derivative of the principal null directions.

Exact solutions for extreme black hole magnetospheres

We present new exact solutions of Force-Free Electrodynamics (FFE) in the Near-Horizon region of an Extremal Kerr black hole (NHEK) and offer a complete classification of the subset that form highest-weight representations of the spacetime's SL(2,R) x U(1) isometry group. For a natural choice of spacetime embedding of this isometry group, the SL(2,R) highest-weight conditions lead to stationary solutions with non-trivial angular dependence, as well as axisymmetry when the U(1)-charge vanishes. In addition, we unveil a hidden SL(2,C) symmetry of the equations of FFE that stems from the action of a complex automorphism group, and enables us to generate an SL(2,C) family of (generically time-dependent) solutions. We then obtain still more general solutions with less symmetry by appealing to a principle of linear superposition that holds for solutions with collinear currents and allows us to resum the highest-weight primaries and their SL(2,R)-descendants.

Force-free electrodynamics around extreme Kerr black holes

Plasma-filled magnetospheres can extract energy from a spinning black hole and provide the power source for a variety of observed astrophysical phenomena. These magnetospheres are described by the highly nonlinear equations of force-free electrodynamics, or FFE. Typically these equations can only be solved numerically. In this paper we consider the FFE equations very near the horizon of a maximally spinning black hole, where the energy extraction takes place. Thanks to an enhanced conformal symmetry which appears in this near-horizon region, we are able to analytically obtain several infinite families of exact solutions of the full nonlinear equations.

Spacetime approach to force-free magnetospheres

Force-Free Electrodynamics (FFE) describes magnetically dominated relativistic plasma via non-linear equations for the electromagnetic field alone. Such plasma is thought to play a key role in the physics of pulsars and active black holes. Despite its simple covariant formulation, FFE has primarily been studied in 3+1 frameworks, where spacetime is split into space and time. In this article we systematically develop the theory of force-free magnetospheres taking a spacetime perspective. Using a suite of spacetime tools and techniques (notably exterior calculus) we cover 1) the basics of the theory, 2) exact solutions that demonstrate the extraction and transport of the rotational energy of a compact object (in the case of a black hole, the Blandford-Znajek mechanism), 3) the behavior of current sheets, 4) the general theory of stationary, axisymmetric magnetospheres and 5) general properties of pulsar and black hole magnetospheres. We thereby synthesize, clarify and generalize known aspects of the physics of force-free magnetospheres, while also introducing several new results.

On the magnetosphere of an accelerated pulsar

We report on a remarkable class of exact solutions to force-free electrodynamics that has four-current along the light cones of an arbitrary timelike worldline in flat spacetime. No symmetry is assumed, and the solutions are given in terms of a free function of three variables. The field configuration should describe the outer magnetosphere of a pulsar moving on the worldline. The power radiated is the sum of an acceleration (Larmor-type) term and a pulsar-type term.

Pulsar spin-down luminosity: Simulations in general relativity

Adopting our new method for matching general relativistic, ideal magnetohydrodynamics to its force-free limit, we perform the first systematic simulations of force-free pulsar magnetospheres in general relativity. We endow the neutron star with a general relativistic dipole magnetic field, model the interior with ideal magnetohydrodynamics, and adopt force-free electrodynamics in the exterior. Comparing the spin-down luminosity to its corresponding Minkowski value, we find that general relativistic effects give rise to a modest enhancement: the maximum enhancement for $n=1$ polytropes is $\ensuremath{\sim}23%$. Evolving a rapidly rotating $n=0.5$ polytrope we find an even greater enhancement of $\ensuremath{\sim}35%$. Using our simulation data, we derive fitting formulas for the pulsar spin-down luminosity as a function of the neutron star compaction, angular speed, and dipole magnetic moment. We expect stiffer equations of state and more rapidly spinning neutron stars to lead to even larger enhancements in the spin-down luminosity.

Exact solutions to force-free electrodynamics in black hole backgrounds

A shared property of several of the known exact solutions to the equations of force-free electrodynamics is that their charge-current four-vector is null. We examine the general properties of null-current solutions and then focus on the principal congruences of the Kerr black hole spacetime. We obtain a large class of exact solutions, which are in general time-dependent and non-axisymmetric. These solutions include waves that, surprisingly, propagate without scattering on the curvature of the black hole’s background. They may be understood as generalizations to Robinson’s solutions to vacuum electrodynamics associated with a shear-free congruence of null geodesics. When stationary and axisymmetric, our solutions reduce to those of Menon and Dermer, the only previously known solutions in Kerr. In Kerr, all of our solutions have null electromagnetic fields ( and E2 = B2). However, in Schwarzschild or flat spacetime there is freedom to add a magnetic monopole field, making the solutions magnetically dominated (B2 > E2). This freedom may be used to reproduce the various flat-spacetime and Schwarzschild-spacetime (split) monopole solutions available in the literature (due to Michel and later authors), and to obtain a large class of time-dependent, non-axisymmetric generalizations. These generalizations may be used to model the magnetosphere of a conducting star that rotates with arbitrary prescribed time-dependent rotation axis and speed. We thus significantly enlarge the class of known exact solutions, while organizing and unifying previously discovered solutions in terms of their null structure.

The Force-Free Electrodynamics Method for the Extrapolation of Coronal Magnetic Fields from Vector Magnetograms

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.