Black holes are known to launch powerful relativistic jets and emit highly variable gamma radiation. How these jets are loaded with plasma remains poorly understood. Spark gaps are thought to drive particle acceleration and pair creation in the black-hole magnetosphere. In this paper, we perform 2D axisymmetric general-relativistic particle-in-cell simulations of a monopole black-hole magnetosphere with a realistic treatment of inverse Compton scattering and pair production. We find that the magnetosphere can self-consistently fill itself with plasma and activate the Blandford-Znajek mechanism. A highly time-dependent spark gap opens near the inner light surface which injects pair plasma into the magnetosphere. These results may account for the high-energy activity observed in active galactic nuclei and explain the origin of plasma at the base of the jet.

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

Black holes drive powerful plasma jets to relativistic velocities. This plasma should be collisionless, and self-consistently supplied by pair creation near the horizon. We present general-relativistic collisionless plasma simulations of Kerr-black-hole magnetospheres which begin from vacuum, inject electron-positron pairs based on local unscreened electric fields, and reach steady states with electromagnetically powered Blandford-Znajek jets and persistent current sheets. Particles with negative energy-at-infinity are a general feature, and can contribute significantly to black-hole rotational-energy extraction in a variant of the Penrose process. The generated plasma distribution depends on the pair-creation environment, and we describe two distinct realizations of the force-free electrodynamic solution. This sensitivity suggests that plasma kinetics will be useful in interpreting future horizon-resolving submillimeter and infrared observations.

Accreting neutron stars can power a wide range of astrophysical phenomena including short- and long-duration gamma-ray bursts, ultra-luminous X-ray sources, and X-ray binaries. Numerical simulations are a valuable tool for studying the accretion-disk--magnetosphere interaction that is central to these problems, most clearly for the recently discovered transitional millisecond pulsars. However, magnetohydrodynamic (MHD) methods, widely used for simulating accretion, have difficulty in highly magnetized stellar magnetospheres, while force-free methods, suitable for such regions, cannot include the accreting gas. We present an MHD method that can stably evolve essentially force-free, highly magnetized regions, and describe the first time-dependent relativistic simulations of magnetized accretion onto millisecond pulsars. Our axisymmetric general-relativistic MHD simulations for the first time demonstrate how the interaction of a turbulent accretion flow with a pulsar's electromagnetic wind can lead to the transition of an isolated pulsar to the accreting state. This transition naturally leads to the formation of relativistic jets, whose power can greatly exceed the power of the isolated pulsar's wind. If the accretion rate is below a critical value, the pulsar instead expels the accretion stream. More generally, our simulations produce for the first time the four possible accretion regimes, in order of decreasing mass accretion rate: (a) crushed magnetosphere and direct accretion; (b) magnetically channeled accretion onto the stellar poles; (c) the propeller state, where material enters through the light cylinder but is prevented from accreting by the centrifugal barrier; (d) almost perfect exclusion of the accretion flow from the light cylinder by the pulsar wind.

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 alpha-disc model. We study a range of disc alpha-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.

The interaction of a rotating star's magnetic field with a surrounding plasma disk lies at the heart of many questions posed by neutron stars in X-ray binaries. We consider the opening of stellar magnetic flux due to differential rotation along field lines coupling the star and disk, using a simple model for the disk-opened flux, the torques exerted on the star by the magnetosphere, and the power extracted by the electromagnetic wind. We examine the conditions under which the system enters an equilibrium spin state, in which the accretion torque is instantaneously balanced by the pulsar wind torque alone. For magnetic moments, spin frequencies, and accretion rates relevant to accreting millisecond pulsars, the spin-down torque from this enhanced pulsar wind can be substantially larger than that predicted by existing models of the disk-magnetosphere interaction, and is in principle capable of maintaining spin equilibrium at frequencies less than 1 kHz. We speculate that this mechanism may account for the non-detection of frequency increases during outbursts of SAX J1808.4-3658 and XTE J1814-338, and may be generally responsible for preventing spin-up to sub-millisecond periods. If the pulsar wind is collimated by the surrounding environment, the resulting jet can satisfy the power requirements of the highly relativistic outflows from Cir X-1 and Sco X-1. In this framework, the jet power scales relatively weakly with accretion rate, L_j ~ Mdot^{4/7}, and would be suppressed at high accretion rates only if the stellar magnetic moment is sufficiently low.

Using an analytic, axisymmetric approach that includes general relativity, coupled to a condition for pair production deduced from simulations, we derive general results about the spatial distribution of pair-producing field lines over the pulsar polar cap. In particular, we show that pair production on magnetic field lines operates over only a fraction of the polar cap for an aligned rotator for general magnetic field configurations, assuming the magnetic field varies spatially on a scale that is larger than the size of the polar cap. We compare our result to force-free simulations of a pulsar with a dipole surface field and find excellent agreement. Our work has implications for first-principles simulations of pulsar magnetospheres, and for explaining observations of pulsed radio and high-energy emission.

The equatorial current sheet in pulsar magnetospheres is often regarded as an ideal site for particle acceleration via relativistic reconnection. Using 2D spherical particle-in-cell simulations, we investigate particle acceleration in the axisymmetric pulsar magnetosphere as a function of the injected plasma multiplicity and magnetization. We observe a clear transition from a highly charge-separated magnetosphere for low plasma injection with little current and spin-down power, to a nearly force-free solution for high plasma multiplicity characterized by a prominent equatorial current sheet and high spin-down power. We find significant magnetic dissipation in the current sheet, up to 30% within 5 light-cylinder radii in the high-multiplicity regime. The simulations unambiguously demonstrate that the dissipated Poynting flux is efficiently channeled to the particles in the sheet, close to the Y-point within about 1-2 light cylinder radii from the star. The mean particle energy in the sheet is given by the upstream plasma magnetization at the light cylinder. The study of particle orbits shows that all energetic particles originate from the boundary layer between the open and the closed field lines. Energetic positrons always stream outward, while high-energy electrons precipitate back towards the star through the sheet and along the separatrices, which may result in auroral-like emission. Our results suggest that the current sheet and the separatrices may be the main source of high-energy radiation in young pulsars.

We propose a scenario for launching relativistic jets from rotating black holes, in which small-scale magnetic flux loops, sustained by disc turbulence, are forced to inflate and open by differential rotation between the black hole and the accretion flow. This mechanism does not require a large-scale net magnetic flux in the accreting plasma. Estimates suggest that the process could operate effectively in many systems, and particularly naturally and efficiently when the accretion flow is retrograde. We present the results of general-relativistic force-free electrodynamic simulations demonstrating the time evolution of the black hole's magnetosphere, the cyclic formation of jets, and the effect of magnetic reconnection. The jets are highly variable on timescales ~ 10-10^3 r_ g/c, where r_g is the black hole's gravitational radius. The reconnecting current sheets observed in the simulations may be responsible for the hard X-ray emission from accreting black holes.

Magnetar magnetospheres are believed to be strongly twisted, due to shearing of the stellar crust by internal magnetic stresses. We present time-dependent axisymmetric simulations showing in detail the evolution of relativistic force-free magnetospheres subjected to slow twisting through large angles. When the twist amplitude is small, the magnetosphere moves quasi-statically through a sequence of equilibria of increasing free energy. At some twist amplitude the magnetosphere becomes tearing-mode unstable to forming a resistive current sheet, initiating large-scale magnetic reconnection in which a significant fraction of the magnetic free energy can be dissipated. This "critical" twist angle is insensitive to the resistive length scale. Rapid shearing temporarily stabilizes the magnetosphere beyond the critical angle, allowing the magnetosphere of a rapidly differentially rotating star to store, and dissipate, more free energy. In addition to these effects, shearing the surface of a rotating star increases the spindown torque applied to the star. If shearing is much slower than rotation, the resulting spikes in spindown rate can occur on timescales anywhere from the long twisting timescale to the stellar spin period or shorter, depending both on the stellar shear distribution and the existing distribution of magnetospheric twists. A model in which energy is stored in the magnetosphere, and released by a magnetospheric instability, therefore predicts large changes in the measured spindown rate before SGR giant flares.

We present the first simulations of evolving, strongly twisted magnetar magnetospheres. Slow shearing of the magnetar crust is seen to lead to a series of magnetospheric expansion and reconnection events, corresponding to X-ray flares and bursts. The axisymmetric simulations include rotation of the neutron star and the magnetic wind through the light cylinder. We study how the increasing twist affects the spindown rate of the star, finding that a dramatic increase in spindown occurs. Particularly spectacular are explosive events caused by the sudden opening of large amounts of overtwisted magnetic flux, which may be associated with the observed giant flares. These events are accompanied by a short period of ultra-strong spindown, resulting in an abrupt increase in spin period, such as was observed in the giant flare of SGR 1900+14.

We describe a new scheme for evolving the equations of force-free electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This pseudospectral code uses global orthogonal basis function expansions to take accurate spatial derivatives, allowing the use of an unstaggered mesh and the complete force-free current density. The method has low numerical dissipation and diffusion outside of singular current sheets. We present a range of one- and two-dimensional tests, and demonstrate convergence to both smooth and discontinuous analytic solutions. As a first application, we revisit the aligned rotator problem, obtaining a steady solution with resistivity localised in the equatorial current sheet outside the light cylinder.

We derive a general expression for the large-scale halo bias, in theories with a scale-dependent linear growth, using the excursion set formalism. Such theories include modified gravity models, and models in which the dark energy clustering is non-negligible. A scale dependence is imprinted in both the formation and evolved biases by the scale-dependent growth. Mergers are accounted for in our derivation, which thus extends earlier work which focused on passive evolution. There is a simple analytic form for the bias for those theories in which the nonlinear collapse of perturbations is approximately the same as in general relativity. As an illustration, we apply our results to a simple Yukawa modification of gravity, and use SDSS measurements of the clustering of luminous red galaxies to constrain the theory's parameters.

Fry (1996) showed that galaxy bias has the tendency to evolve towards unity, i.e. in the long run, the galaxy distribution tends to trace that of matter. Generalizing slightly Fry's reasoning, we show that his conclusion remains valid in theories of modified gravity (or equivalently, complex clustered dark energy). This is not surprising: as long as both galaxies and matter are subject to the same force, dynamics would drive them towards tracing each other. This holds, for instance, in theories where both galaxies and matter move on geodesics. This relaxation of bias towards unity is tempered by cosmic acceleration, however: the bias tends towards unity but does not quite make it, unless the formation bias were close to unity. Our argument is extended in a straightforward manner to the case of a stochastic or nonlinear bias. An important corollary is that dynamical evolution could imprint a scale dependence on the large scale galaxy bias. This is especially pronounced if non-standard gravity introduces new scales to the problem: the bias at different scales relaxes at different rates, the larger scales generally more slowly and retaining a longer memory of the initial bias. A consistency test of the current (general relativity + uniform dark energy) paradigm is therefore to look for departure from a scale independent bias on large scales. A simple way is to measure the relative bias of different populations of galaxies which are at different stages of bias relaxation. Lastly, we comment on the possibility of directly testing the Poisson equation on cosmological scales, as opposed to indirectly through the growth factor.

Solar active regions, produced by the emergence of tubes of strong magnetic field in the photosphere, are restricted to within 35 degrees of the solar equator. The nature of the dynamo processes that create and renew these fields, and are therefore responsible for solar magnetic phenomena, are not well understood. We analyze the magneto-rotational stability of the solar tachocline for general field geometry. This thin region of strong radial and latitudinal differential rotation, between the radiative and convective zones, is unstable at latitudes above 37 degrees, yet is stable closer to the equator. We propose that small-scale magneto-rotational turbulence prevents coherent magnetic dynamo action in the tachocline except in the vicinity of the equator, thus explaining the latitudinal restriction of active regions. Tying the magnetic dynamo to the tachocline elucidates the physical conditions and processes relevant to solar magnetism.