In 2007 September, the Crab Nebula exhibited a bright gamma-ray flare in the GeV energy range that was detected by AGILE. The observed emission at >160 MeV indicates that the radiating electrons had energies above the classical synchrotron radiation-reaction limit, thus presenting a serious challenge to classical models for electron acceleration in astrophysical environments. In this paper, we apply our recently developed time-dependent self-similar analytical model describing electrostatic acceleration in the explosive reconnection region around the pulsar wind termination shock to the 2007 September flare. This event was unique in that it displayed both long-duration "wave" and short-duration "sub-flare" features. The unusual temporal variation makes this flare an especially interesting test for our model. We demonstrate that our model can reproduce the time-dependent gamma-ray spectrum for this event, as well as the associated gamma-ray light curve, obtained by integrating the spectrum for photon energies >100 MeV. This establishes that our time-dependent electrostatic acceleration model can explain both wave and sub-flare transients, which lends further support to the theoretical framework we have developed. We also further examine the validity of the self-similar electric and magnetic field evolution implied by our model. We conclude that strong electrostatic acceleration driven by shock-induced magnetic reconnection is able to power the Crab Nebula gamma-ray flares by energizing the electrons on sub-Larmor timescales.

The $\gamma$-ray flares from the Crab nebula observed by {\it AGILE} and {\it Fermi}-LAT between 2007-2013 reached GeV photon energies and lasted several days. The strongest emission, observed during the 2011 April "super-flare," exceeded the quiescent level by more than an order of magnitude. These observations challenge the standard models for particle acceleration in pulsar wind nebulae, because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Particle-in-cell simulations have suggested that the classical synchrotron limit can be exceeded if the electrons also experience electrostatic acceleration due to shock-driven magnetic reconnection. In this paper, we revisit the problem using an analytic approach based on solving a fully time-dependent electron transport equation describing the electrostatic acceleration, synchrotron losses, and escape experienced by electrons in a magnetically confined plasma "blob" as it encounters and passes through the pulsar-wind termination shock. We show that our model can reproduce the $\gamma$-ray spectra observed during the rising and decaying phases of each of the two sub-flare components of the 2011 April super-flare. We integrate the spectrum for photon energies $\ge 100\,$MeV to obtain the light curve for the event, which agrees with the observations. We find that strong electrostatic acceleration occurs on both sides of the termination shock, driven by magnetic reconnection. We also find that the dominant mode of particle escape changes from diffusive escape to advective escape as the blob passes through the shock.

The {\gamma}-ray flares from the Crab nebula observed by AGILE and Fermi-LAT reaching GeV energies and lasting several days challenge the standard models for particle acceleration in pulsar wind nebulae, because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Previous modeling has suggested that the synchrotron limit can be exceeded if the electrons experience electrostatic acceleration, but the resulting spectra do not agree very well with the data. As a result, there are still some important unanswered questions about the detailed particle acceleration and emission processes occurring during the flares. We revisit the problem using a new analytical approach based on an electron transport equation that includes terms describing electrostatic acceleration, stochastic wave-particle acceleration, shock acceleration, synchrotron losses, and particle escape. An exact solution is obtained for the electron distribution, which is used to compute the associated {\gamma}-ray synchrotron spectrum. We find that in our model the {\gamma}-ray flares are mainly powered by electrostatic acceleration, but the contributions from stochastic and shock acceleration play an important role in producing the observed spectral shapes. Our model can reproduce the spectra of all the Fermi-LAT and AGILE flares from the Crab nebula, using magnetic field strengths in agreement with the multi-wavelength observational constraints. We also compute the spectrum and duration of the synchrotron afterglow created by the accelerated electrons, after they escape into the region on the downstream side of the pulsar wind termination shock. The afterglow is expected to fade over a maximum period of about three weeks after the {\gamma}-ray flare.

Many accreting black holes manifest time lags during outbursts, in which the hard Fourier component typically lags behind the soft component. Despite decades of observations of this phenomenon, the underlying physical explanation for the time lags has remained elusive, although there are suggestions that Compton reverberation plays an important role. However, the lack of analytical solutions has hindered the interpretation of the available data. In this paper, we investigate the generation of X-ray time lags in Compton scattering coronae using a new mathematical approach based on analysis of the Fourier-transformed transport equation. By solving this equation, we obtain the Fourier transform of the radiation Green's function, which allows us to calculate the exact dependence of the time lags on the Fourier frequency, for both homogeneous and inhomogeneous coronal clouds. We use the new formalism to explore a variety of injection scenarios, including both monochromatic and broadband (bremsstrahlung) seed photon injection. We show that our model can successfully reproduce both the observed time lags and the time-averaged (quiescent) X-ray spectra for Cyg~X-1 and GX~339-04, using a single set of coronal parameters for each source. The time lags are the result of impulsive bremsstrahlung injection occurring near the outer edge of the corona, while the time-averaged spectra are the result of continual distributed injection of soft photons throughout the cloud.

We examine the stability of an Einstein-Maxwell perfect fluid configuration with a privileged direction of symmetry by means of a $1+1+2$-tetrad formalism. We use this formalism to cast, in a quasi linear symmetric hyperbolic form the equations describing the evolution of the system. This hyperbolic reduction is used to discuss the stability of solutions of the linear perturbation. By restricting the analysis to isotropic fluid configurations, we made use of a constant electrical conductivity coefficient for the fluid (plasma), and the nonlinear stability for the case of an infinitely conducting plasma is also considered. As a result of this analysis we provide a complete classification and characterization of various stable and unstable configurations. We found in particular that in many cases the stability conditions is strongly determined by the constitutive equations by means of the square of the velocity of sound and the electric conductivity, and a threshold for the emergence of the instability appears in both contracting and expanding systems.

Accreting black hole sources show a wide variety of rapid time variability, including the manifestation of time lags during X-ray transients, in which a delay (phase shift) is observed between the Fourier components of the hard and soft spectra. Despite a large body of observational evidence for time lags, no fundamental physical explanation for the origin of this phenomenon has been presented. We develop a new theoretical model for the production of X-ray time lags based on an exact analytical solution for the Fourier transform describing the diffusion and Comptonization of seed photons propagating through a spherical corona. The resulting Green's function can be convolved with any source distribution to compute the associated Fourier transform and time lags, hence allowing us to explore a wide variety of injection scenarios. We show that thermal Comptonization is able to self-consistently explain both the X-ray time lags and the steady-state (quiescent) X-ray spectrum observed in the low-hard state of Cyg X-1. The reprocessing of bremsstrahlung seed photons produces X-ray time lags that diminish with increasing Fourier frequency, in agreement with the observations for a wide range of sources.

We consider a complex scalar field minimally coupled to gravity and to a U(1) gauge symmetry and we construct of a first order symmetric hyperbolic evolution system for the Einstein-Maxwell-Klein-Gordon system. Our analysis is based on a 1+3 tetrad formalism which makes use of the components of the Weyl tensor as one of the unknowns. In order to ensure the symmetric hyperbolicity of the evolution equations, implied by the Bianchi identity, we introduce a tensor of rank 3 corresponding to the covariant derivative of the Faraday tensor, and two tensors of rank 2 for the covariant derivative of the vector potential and the scalar field.

We examine the problem of the construction of a first order symmetric hyperbolic evolution system for the Einstein-Maxwell-Euler system. Our analysis is based on a 1+3 tetrad formalism which makes use of the components of the Weyl tensor as one of the unknowns. In order to ensure the symmetric hyperbolicity of the evolution equations implied by the Bianchi identity, we introduce a tensor of rank 3 corresponding to the covariant derivative of the Faraday tensor. Our analysis includes the case of a perfect fluid with infinite conductivity (ideal magnetohydrodynamics) as a particular subcase.