Normalized to: Potter, W.
[1]
oai:arXiv.org:1810.04139 [pdf] - 1897854
The Feasibility of Magnetic Reconnection Powered Blazar Flares from
Synchrotron Self-Compton Emission
Submitted: 2018-10-09, last modified: 2019-06-11
Order of magnitude variability has been observed in the blazar sub-class of
Active Galactic Nuclei on minute timescales. These high-energy flares are often
difficult to explain with shock acceleration models due to the small size of
the inferred emitting region, with recent particle-in-cell (PIC) simulations
showing that magnetic reconnection is a promising alternative mechanism. Here,
we present a macroscopic emission model physically motivated by PIC
simulations, where the energy for particle acceleration originates from the
reconnecting magnetic field. We track the radial growth and relative velocity
of a reconnecting plasmoid, modelling particle acceleration and radiative
losses from synchrotron and synchrotron self-Compton (SSC) emission. To test
the viability of magnetic reconnection as the mechanism behind rapid blazar
flares we simultaneously fit our model to the observed light-curve and SED from
the 2016 TeV flare of BL Lacertae. We find generally that, without considering
external photons, reconnecting plasmoids are unable to produce Compton-dominant
TeV flares and so cannot reproduce the observations due to overproduction of
synchrotron emission. Additionally, problematically large plasmoids, comparable
in size to the entire jet radius, are required to emit sufficient SSC
gamma-rays to be observable. However, our plasmoid model can reproduce the
rapid TeV lightcurve of the flare, demonstrating that reconnection is able to
produce rapid, powerful TeV flares on observed timescales. We conclude that
while reconnection can produce SSC flares on the correct timescales, the
primary source of TeV emission cannot be SSC and the size of plasmoids required
may be implausibly large.
[2]
oai:arXiv.org:1709.09188 [pdf] - 1598172
Modelling blazar flaring using a time-dependent fluid jet emission model
- an explanation for orphan flares and radio lags
Submitted: 2017-09-26
Blazar jets are renowned for their rapid violent variability and
multiwavelength flares, however, the physical processes responsible for these
flares are not well understood. In this paper we develop a time-dependent
inhomogeneous fluid jet emission model for blazars. We model optically thick
radio flares for the first time and show that they are delayed with respect to
the prompt optically thin emission by ~ months to decades, with a lag that
increases with the jet power and observed wavelength. This lag is caused by a
combination of the travel time of the flaring plasma to the optically thin
radio emitting sections of the jet and the slow rise time of the radio flare.
We predict two types of flares: symmetric flares - with the same rise and decay
time, which occur for flares whose duration is shorter than both the radiative
lifetime and the geometric path-length delay timescale; extended flares - whose
luminosity tracks the power of particle acceleration in the flare, which occur
for flares with a duration longer than both the radiative lifetime and
geometric delay. Our model naturally produces orphan X-ray and $\gamma$-ray
flares. These are caused by flares which are only observable above the
quiescent jet emission in a narrow band of frequencies. Our model is able to
successfully fit to the observed multiwavelength flaring spectra and
lightcurves of PKS1502+106 across all wavelengths, using a transient flaring
front located within the broad-line region.
[3]
oai:arXiv.org:1404.0013 [pdf] - 1579516
An accretion disc instability induced by a temperature sensitive
$\alpha$ parameter
Submitted: 2014-03-31, last modified: 2017-08-07
In the standard thin disc formalism the dimensionless $\alpha$ parameter is
usually assumed to be constant. However, there are good theoretical reasons for
believing, as well as evidence from simulations, that $\alpha$ is dependent on
intrinsic disc properties. In this paper we analyse the conditions for the
stability of a thin accretion disc in which $\alpha$ is a function of the
magnetic Prandtl number, the ratio of collisional viscosity to resistivity. In
the inner disc, where the free electron opacity and radiation viscosity
dominate, the disc is unstable if $\alpha$ is proportional to the magnetic
Prandtl number with an exponent, $n$, and $6/13<n<10/3$. This is within the
range of values for the power-law index found in MHD simulations with simple
energetics. We calculate the evolution of the unstable disc within the $\alpha$
formalism and show that the physically accessible solutions form a limit cycle,
analogous to the behaviour seen in recurrent dwarf novae. It is noteworthy that
the time-dependent global behaviour of the instability results in cyclic
heating of the inner section of the disc, when parameters appropriate for an
X-ray binary system are used. We calculate a model spectrum of the disc in the
flaring and quiescent states and show that the behaviour is compatible with
X-ray observations of the thermal accretion disc in flaring X-ray binary
systems.
[4]
oai:arXiv.org:1704.02485 [pdf] - 1582435
Demonstration of a magnetic Prandtl number disc instability from first
principles
Submitted: 2017-04-08
Understanding what determines the strength of MHD turbulence in accretion
discs is a question of fundamental theoretical and observational importance. In
this work we investigate whether the dependence of the turbulent accretion disc
stress ($\alpha$) on the magnetic Prandtl number (Pm) is sufficiently sensitive
to induce thermal-viscous instability using 3D MHD simulations. We first
investigate whether the $\alpha$-Pm dependence, found by many previous authors,
has a physical or numerical origin by conducting a suite of local shearing-box
simulations. We find that a definite $\alpha$-Pm dependence persists when
simultaneously increasing numerical resolution and decreasing the absolute
values of both the viscous and resistive dissipation coefficients. This points
to a physical origin of the $\alpha$-Pm dependence. Using a further set of
simulations which include realistic turbulent heating and radiative cooling,
and by giving Pm a realistic physical dependence on the plasma temperature and
density, we demonstrate that the $\alpha$-Pm dependence is sufficiently strong
to lead to a local instability. We confirm that the instability manifests
itself as an unstable limit cycle by mapping the local thermal-equilibrium
curve of the disc. This is the first self-consistent MHD simulation
demonstrating the Pm instability from first principles. This result is
important because a physical Pm instability would lead to the global
propagation of heating and cooling fronts and a transition between disc states
on timescales compatible with the observed hard/soft state transitions in black
hole binaries.
[5]
oai:arXiv.org:1610.08321 [pdf] - 1532255
Using radiative energy losses to constrain the magnetisation and
magnetic reconnection rate at the base of black hole jets
Submitted: 2016-10-26
We calculate the severe radiative energy losses which occur at the base of
black hole jets using a relativistic fluid jet model, including in-situ
acceleration of non-thermal leptons by magnetic reconnection. Our results
demonstrate that including a self-consistent treatment of radiative energy
losses is necessary to perform accurate MHD simulations of powerful jets and
that jet spectra calculated via post-processing are liable to vastly
overestimate the amount of non-thermal emission. If no more than 95% of the
initial total jet power is radiated away by the plasma travels as it travels
along the length of the jet, we can place a lower bound on the magnetisation of
the jet plasma at the base of the jet. For typical powerful jets, we find that
the plasma at the jet base is required to be highly magnetised, with at least
10,000 times more energy contained in magnetic fields than in non-thermal
leptons. Using a simple power-law model of magnetic reconnection, motivated by
simulations of collisionless reconnection, we determine the allowed range of
the large-scale average reconnection rate along the jet, by restricting the
total radiative energy losses incurred and the distance at which the jet first
comes into equipartition. We calculate analytic expressions for the cumulative
radiative energy losses due to synchrotron and inverse-Compton emission along
jets, and derive analytic formulae for the constraint on the initial
magnetisation.
[7]
oai:arXiv.org:1603.06489 [pdf] - 1411424
Surprises in astrophysical gasdynamics
Submitted: 2016-03-21
Much of astrophysics consists of the study of ionised gas under the influence
of gravitational and magnetic fields. Thus, it is not possible to understand
the astrophysical universe without a detailed knowledge of the dynamics of
magnetised fluids. Fluid dynamics is, however, a notoriously tricky subject, in
which it is all too easy for one's a priori intuition to go astray. In this
review, we seek to guide the reader through a series of illuminating yet
deceptive problems, all with an enlightening twist. We cover a broad range of
topics including the instabilities acting in accretion discs, the hydrodynamics
governing the convective zone of the Sun, the magnetic shielding of a cooling
galaxy cluster, and the behaviour of thermal instabilities and evaporating
clouds. The aim of this review is to surprise and intrigue even veteran
astrophysical theorists with an idiosynchratic choice of problems and
counterintuitive results. At the same time, we endeavour to bring forth the
fundamental ideas, to set out important assumptions, and to describe carefully
whatever novel techniques may be appropriate to the problem at hand. By
beginning at the beginning, and analysing a wide variety of astrophysical
settings, we seek not only to make this review suitable for fluid dynamic
veterans, but to engage novice recruits as well with what we hope will be an
unusual and instructive introduction to the subject.
[9]
oai:arXiv.org:1508.00567 [pdf] - 1284908
New constraints on the structure and dynamics of black hole jets
Submitted: 2015-08-03
Accreting black holes produce powerful relativistic plasma jets which emit
radiation across all observable wavelengths but the details of the initial
acceleration and confinement of the jet are uncertain. We apply an innovative
new model that allows us to determine key properties of the acceleration zone
via multi-frequency observations. The central component of the model is a
relativistic steady-state fluid flow, and the emission from physically distinct
regions can be seen to contribute to different energy bands in the overall
spectrum. By fitting with unprecedented accuracy to 42 simultaneous
multiwavelength blazar spectra we are able to constrain the location of the
brightest synchrotron emitting region, and show that there must be a linear
relation between the jet power and the radius of the brightest region of the
jet. We also find a correlation between the length of the accelerating region
and the maximum bulk Lorentz factor of the jet and find evidence for a bimodal
distribution of accretion rates in the blazar population. This allows us to put
constraints on the basic dynamical and structural properties of blazar jets and
to understand the underlying physical differences which result in the blazar
sequence.
[10]
oai:arXiv.org:1502.06479 [pdf] - 1231834
Differential rotation and radiative equilibrium in the Sun: is the
tachocline spreading?
Submitted: 2015-02-23
It is well known that the combination of barotropic rotation and radiative
equilibrium are mutually incompatible in stars. The Sun's internal rotation is
far from barotropic, however, which allows at least the theoretical possibility
that the Sun's thermal balance is one of radiative equilibrium in the region of
the tachocline near the outer boundary of the radiative zone. We show here that
(i) the constraint of radiative equilibrium leads to a straightforward ordinary
differential equation for the Sun's rotation profile, and (ii) solutions of
this equation can be found that, to within current levels of accuracy, closely
resemble the rotation profile deduced from helioseismology. More generally, we
calculate how large a baroclinic deviation from uniform rotation is required to
maintain radiative equilibrium without meridional circulation throughout the
bulk of the radiative zone. Very little deviation is required, well below
detectability. The feasibility of radiative equilibrium for the tachocline
suggests that the issue of a spreading boundary layer may be less severe than
previously thought.
[11]
oai:arXiv.org:1311.2617 [pdf] - 1180645
Uncovering the physics behind the blazar sequence using a realistic
model for jet emission
Submitted: 2013-11-11
Blazar spectra are one of the most important windows into the physical
processes occurring along jets. The spectrum, composed from the different
emitting regions along the jet, allows us to constrain the physical conditions
in the jet. I present my work modelling blazar spectra using an extended
inhomogeneous jet model with an accelerating, magnetically dominated, parabolic
base transitioning to a slowly decelerating, conical section motivated by
observations, simulations and theory. We set the inner geometry of our
multi-zone model using observations of the jet in M87 which transitions from
parabolic to conical at 10^5 Schwarzschild radii. This model is able to
reproduce quiescent blazar spectra very well across all wavelengths (including
radio observations) for a sample of 42 BL Lacs and FSRQs.
Using this inhomogeneous model we are able to constrain the location at which
the synchrotron emission is brightest in these jets by fitting to the optically
thick to thin synchrotron break. We find that the radius of the jet at which
the synchrotron emission is brightest (where the jet first approaches
equipartition) scales approximately linearly with the jet power. We also find a
correlation between the length of the accelerating, parabolic section of the
jet and the maximum bulk Lorentz factor. In agreement with previous work we
find that BL Lacs are low power blazars whereas FSRQs are high power blazars.
Together with our simple jet power-radius relation this leads us to a deeper
understanding of the physics underlying the blazar sequence.
[12]
oai:arXiv.org:1310.0462 [pdf] - 1179641
Synchrotron and inverse-Compton emission from blazar jets - IV. BL Lac
type blazars and the physical basis for the blazar sequence
Submitted: 2013-10-01
In this paper we investigate the properties of a sample of six BL Lacs by
fitting their spectra using our inhomogeneous jet model with an accelerating,
magnetically dominated, parabolic base, which transitions to a slowly
decelerating conical jet with a geometry based on observations of M87. Our
model is able to fit very well to the simultaneous multiwavelength spectra of
all the BL Lacs including radio observations. We find that the BL Lacs have
lower jet powers and bulk Lorentz factors than the sample of Compton-dominant
blazars investigated in Paper III, consistent with the blazar sequence.
Excitingly, we find a correlation between the radius at which the jet first
comes into equipartition and the jet power, in agreement with our prediction
from Paper II. We interpret this result as one of two physical scenarios: a
universal jet geometry which scales linearly with black hole mass or a
dichotomy in Eddington accretion rates between flat spectrum radio quasars
(FSRQs) and BL Lacs. If we assume that the jet geometry of all blazars scales
linearly with black hole mass then we find a plausible range of masses
(10^7-10^10 solar masses). We find that the quiescent gamma-ray spectrum of
Markarian 421 is best fitted by scattering of external CMB photons. We are
unable to fit the spectrum using synchrotron self-Compton emission due to
strong gamma-ray absorption via pair production even using a compact, rapidly
decelerating, jet with a very large bulk Lorentz factor (50), as has been
suggested recently. Finally, we fit to the SEDs of the four high power HSP BL
Lacs recently found by Padovani et al. 2012. We find that their high peak
frequency emission is caused by high maximum electron energies whilst the rest
of their jet properties are typical of relatively high power BL Lacs and
consistent with our predictions.
[13]
oai:arXiv.org:1310.0006 [pdf] - 1179590
The Cosmic Evolution of Fermi BL Lacertae Objects
Ajello, M.;
Romani, R. W.;
Gasparrini, D.;
Shaw, M. S.;
Bolmer, J.;
Cotter, G.;
Finke, J.;
Greiner, J.;
Healey, S. E.;
King, O.;
Max-Moerbeck, W.;
Michelson, P. F.;
Potter, W. J.;
Rau, A.;
Readhead, A. C. S.;
Richards, J. L.;
Schady, P.
Submitted: 2013-09-30
Fermi has provided the largest sample of gamma-ray selected blazars to date.
In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac)
objects detected by it Fermi during its first year of operation. We have
obtained redshift constraints for 206 out of the 211 BL Lacs in our sample
making it the largest and most complete sample of BL Lacs available in the
literature. We use this sample to determine the luminosity function of BL Lacs
and its evolution with cosmic time. We find that for most BL Lac classes, the
evolution is positive with a space density peaking at modest redshift (z~1.2).
The low-luminosity, high-synchrotron peaked (HSP) BL Lacs are an exception,
showing strong negative evolution, with number density increasing for
z$\lesssim$0.5. Since this rise corresponds to a drop-off in the density of
flat-spectrum radio quasars (FSRQs), a possible interpretation is that these
HSPs represent an accretion-starved end-state of an earlier merger-driven
gas-rich phase. We additionally find that the known BL Lac correlation between
luminosity and photon spectral index persists after correction for the
substantial observational selection effects with implications for the so called
`blazar sequence'. Finally, estimating the beaming corrections to the
luminosity function, we find that BL Lacs have an average Lorentz factor of
$\gamma=6.1^{+1.1}_{-0.8}$, and that most are seen within 10$^{\circ}$ of the
jet axis.
[15]
oai:arXiv.org:1303.1182 [pdf] - 1165024
Synchrotron and inverse-Compton emission from blazar jets - III.
Compton-dominant blazars
Submitted: 2013-03-05
In this paper we develop the extended jet model of Potter & Cotter to model
the simultaneous multi-wavelength spectra of six Compton-dominant blazars. We
include an accelerating parabolic base transitioning to a slowly decelerating
conical jet with a geometry set by observations of M87 and consistent with
simulations and theory. We investigate several jet models and find that the
optically thick to thin synchrotron break in the radio spectrum requires the
jet to first come into equipartition at large distances along the jet (10^5
Schwarzschild radii), consistent with the observed transition from parabolic to
conical in the jet of M87. We confirm this result analytically and calculate
the expected frequency core-shift relations for the models under consideration.
We find that a parabolic jet transitioning to a ballistic conical jet, which
starts in equipartition and becomes more particle dominated at larger
distances, fits the multiwavelength data of the six blazars well, whilst an
adiabatic equipartition conical section requires very large bulk Lorentz
factors to reproduce the Compton-dominance of the blazars.
We find that all these blazars require high power, high bulk Lorentz factor
jets observed close to the line of sight as we expect from the blazar sequence
and consistent with the results from Paper II. The inverse-Compton emission in
our fits is due to inverse-Compton scattering of high-redshift CMB photons at
large distances along the jet due to the high bulk Lorentz factors of the jets.
We postulate a new interpretation of the blazar sequence based on the radius of
the transition region of the jet (where the jet is brightest in synchrotron
emission) scaling linearly with black hole mass.
[16]
oai:arXiv.org:1301.0323 [pdf] - 1158783
Spectroscopy of The Largest Ever Gamma-ray Selected BL Lac Sample
Shaw, Michael S.;
Romani, Roger W.;
Cotter, Garret;
Healey, Stephen E.;
Michelson, Peter F.;
Readhead, Anthony C. S.;
Richards, Joseph L.;
Max-Moerbeck, Walter;
King, Oliver G.;
Potter, William J.
Submitted: 2013-01-02
We report on spectroscopic observations covering most of the 475 BL Lacs in
the 2nd Fermi LAT catalog of AGN. Including archival measurements (correcting
several erroneous literature values) we now have spectroscopic redshifts for
44% of the BL Lacs. We establish firm lower redshift limits via intervening
absorption systems and statistical lower limits via searches for host galaxies
for an additional 51% of the sample leaving only 5% of the BL Lacs
unconstrained. The new redshifts raise the median spectroscopic z from 0.23 to
0.33 and include redshifts as large as z=2.471. Spectroscopic redshift minima
from intervening absorbers have ~ z= 0.70, showing a substantial fraction at
large z and arguing against strong negative evolution. We find that detected BL
Lac hosts are bright ellipticals with black hole masses M_\bullet ~ 10^{8.5-9},
substantially larger than the mean of optical AGN and LAT Flat Spectrum Radio
Quasar samples. A slow increase in M_\bullet with z may be due to selection
bias. We find that the power-law dominance of the optical spectrum extends to
extreme values, but this does not strongly correlate with the gamma-ray
properties, suggesting that strong beaming is the primary cause of the range in
continuum dominance.
[17]
oai:arXiv.org:1212.2632 [pdf] - 1158436
Synchrotron and inverse-Compton emission from blazar jets - II. An
accelerating jet model with a geometry set by observations of M87
Submitted: 2012-12-11
In this paper we develop the jet model of Potter & Cotter (2012) to include a
magnetically dominated accelerating parabolic base transitioning to a slowly
decelerating conical jet with a geometry set by recent radio observations of
M87. We conserve relativistic energy-momentum and particle number along the jet
and calculate the observed synchrotron emission from the jet by calculating the
integrated line of sight synchrotron opacity through the jet in the rest frame
of each section of plasma. We calculate the inverse-Compton emission from
synchrotron, CMB, accretion disc, starlight, broad line region, dusty torus and
narrow line region photons by transforming into the rest frame of the plasma
along the jet.
We fit our model to simultaneous multi-wavelength observations of the
Compton-dominant FSRQ type blazar PKS0227-369, with a jet geometry set by M87
and an accelerating bulk Lorentz factor consistent with simulations and theory.
We investigate models in which the jet comes into equipartition at different
distances along the jet and equipartition is maintained via the conversion of
jet bulk kinetic energy into particle acceleration. We find that the jet must
still be magnetically dominated within the BLR and cannot be in equipartition
due to the severe radiative energy losses. The model fits the observations,
including radio data, very well if the jet comes into equipartition outside the
BLR within the dusty torus (1.5pc) or at further distances (34pc). We find that
our fit in which the jet comes into equipartition furthest along the jet, which
has a jet with the geometry of M87 scaled linearly with black hole mass, has an
inferred black hole mass close to previous estimates. This implies that the jet
of PKS0227 might be well described by the same jet geometry as M87.
[18]
oai:arXiv.org:1203.3881 [pdf] - 602216
Synchrotron and inverse-Compton emission from blazar jets I: a uniform
conical jet model
Submitted: 2012-03-17
In the first of a series of papers investigating emission from blazar jets
from radio to high-energy {\gamma}-rays, we revisit the class of models where
the jet has a uniform conical ballistic structure. We argue that by using
simple developments of these models, in the context of new multi-frequency data
extending to gamma-ray energies, valuable insights may be obtained into the
properties that fully realistic models must ultimately have. In this paper we
consider the synchrotron and synchrotron-self-Compton emission from the jet,
modelling the recent simultaneous multi-wavelength observations of BL Lac. This
is the first time these components have been fitted simultaneously for a blazar
using a conical jet model.
In the model we evolve the electron population dynamically along the jet
taking into account the synchrotron and inverse-Compton losses. The
inverse-Compton emission is calculated using the Klein-Nishina cross section
and a relativistic transformation into the jet frame, and we explicitly show
the seed photon population. We integrate synchrotron opacity along the line of
sight through the jet plasma, taking into account the emission and opacity of
each section of the jet. In agreement with previous studies of radio emission,
we find that a conical jet model which conserves magnetic energy produces the
characteristic blazar flat radio spectrum, however, we do not require any
fine-tuning of the model to achieve this. Of particular note, in our model fit
to BL Lac--which at ~10^37W is a relatively low jet-power source--we find no
requirement for significant re-acceleration within the jet to explain the
observed spectrum.
[19]
oai:arXiv.org:1201.0999 [pdf] - 1092729
Spectroscopy of Broad Line Blazars from 1LAC
Shaw, Michael S.;
Romani, Roger W.;
Cotter, Garret;
Healey, Stephen E.;
Michelson, Peter F.;
Readhead, Anthony C. S.;
Richards, Joseph L.;
Max-Moerbeck, Walter;
King, Oliver G.;
Potter, William J.
Submitted: 2012-01-04
We report on optical spectroscopy of 165 Flat Spectrum Radio Quasars (FSRQs)
in the Fermi 1LAC sample, which have helped allow a nearly complete study of
this population. Fermi FSRQ show significant evidence for non-thermal emission
even in the optical; the degree depends on the gamma-ray hardness. They also
have smaller virial estimates of hole mass than the optical quasar sample. This
appears to be largely due to a preferred (axial) view of the gamma-ray FSRQ and
non-isotropic (H/R ~ 0.4) distribution of broad-line velocities. Even after
correction for this bias, the Fermi FSRQ show higher mean Eddington ratios than
the optical population. A comparison of optical spectral properties with Owens
Valley Radio Observatory radio flare activity shows no strong correlation.
[20]
oai:arXiv.org:1108.4414 [pdf] - 406909
A gauge-invariant approach to interactions in the dark sector
Submitted: 2011-08-22
We outline a gauge-invariant framework to calculate cosmological
perturbations in dark energy models consisting of a scalar field interacting
with dark matter via energy and momentum exchanges. Focusing on three
well-known models of quintessence and three common types of dark-sector
interactions, we calculate the matter and dark energy power spectra as well as
the Integrated Sachs-Wolfe (ISW) effect in these models. We show how the
presence of dark-sector interactions can produce a large-scale enhancement in
the matter power spectrum and a boost in the low multipoles of the cosmic
microwave background anisotropies. Nevertheless, we find these enhancements to
be much more subtle than those found by previous authors who model dark energy
using simple ansatz for the equation of state. We also address issues of
instabilities and emphasise the importance of momentum exchanges in the dark
sector.
