Source author record

Peter A. Becker

Peter A. Becker appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2020arXiv

A two-fluid model for black-hole accretion flows: particle acceleration and disc structure

Hot, tenuous advection-dominated accretion flows around black holes are ideal sites for the Fermi acceleration of relativistic particles at standing shock waves in the accretion disc. Previous work has demonstrated that the shock-acceleration process can be efficient enough to power the observed, strong outflows in radio-loud active galaxies such as M87. However, the dynamical effect (back-reaction) on the flow, exerted by the pressure of the relativistic particles, has not been previously considered, and this effect can have a significant influence on the disc structure. We reexamine the problem by developing a new, two-fluid model for the structure of the accretion disc that includes the dynamical effect of the relativistic particle pressure, combined with the pressure of the background (thermal) gas. The new model is analogous to the two-fluid model of cosmic ray acceleration in supernova-driven shock waves. As part of the model, we also develop a new set of shock jump conditions, which are solved along with the hydrodynamic conservation equations to determine the structure of the accretion disc. The solutions include the formation of a mildly relativistic outflow (jet) at the shock radius, driven by the relativistic particles accelerated in the disc. One of our main conclusions is that in the context of the new two-fluid accretion model, global smooth (shock-free) solutions do not exist, and the disc must always contain a standing shock wave, at least in the inviscid case considered here.

preprint2020arXiv

A two-fluid model for black-hole accretion flows: Particle acceleration, outflows, and TeV emission

The multi-wavelength spectrum observed from M87 extends from radio wavelengths up to TeV gamma-ray energies. The radio through GeV components have been interpreted successfully using SSC models based on misaligned blazar jets, but the origin of the intense TeV emission detected during flares in 2004, 2005, and 2010 remains puzzling. It has been previously suggested that the TeV flares are produced when a relativistic proton jet originating in the core of M87 collides with a molecular cloud (or stellar atmosphere) located less than one parsec from the central black hole. We explore this scenario in detail here using a self-consistent model for the acceleration of relativistic protons in a shocked, two-fluid ADAF accretion disc. The relativistic protons accelerated in the disc escape to power the observed jet outflows. The distribution function for the jet protons is used to compute the TeV emission produced when the jet collides with a cloud or stellar atmosphere. The simulated broadband radiation spectrum includes radio, X-ray, and GeV components generated via synchrotron, as well as TeV emission generated via the production and decay of muons, positrons, and electrons. The self-consistency of the model is verified by computing the relativistic particle pressure using the distribution function, and comparing it with the relativistic particle pressure obtained from the hydrodynamical model. We demonstrate that the model is able to reproduce the multi-wavelength spectrum from M87 observed by VERITAS and HESS during the high-energy flares in 2004, 2005, and 2010.

preprint2019arXiv

The giant outburst of 4U 0115+634 in 2011 with Suzaku and RXTE

We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

preprint2016arXiv

An Integrated Model for the Production of X-Ray Time Lags and Quiescent Spectra from Homogeneous and Inhomogeneous Black Hole Accretion Coronae

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.

preprint2016arXiv

Standing Shock Instability in Advection-Dominated Accretion Flows

Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura linearization method and the Nakayama instability boundary conditions, we discover that there are regions of parameter space where disks/shocks with outflows can be stable or unstable. In regions of instability, we find that preshock deceleration is always unstable to the zeroth mode with zero frequency of oscillation, but is always stable to the fundamental mode and overtones. Furthermore, we also find that preshock acceleration is always unstable to the zeroth mode and that the fundamental mode and overtones become increasingly less stable as the shock location moves away from the horizon when the disk half-height expands above $\sim 12$ gravitational radii at the shock radius. In regions of stability, we demonstrate the zeroth mode to be stable for the velocity profiles that exhibit preshock acceleration and deceleration. Moreover, for models that are linearly unstable, our model suggests the possible existence of quasi-periodic oscillations (QPOs) with ratios 2:3 and 3:5. These ratios are believed to occur in stellar and supermassive black hole candidates, for example, in GRS 1915+105 and Sgr A*, respectively. We expect that similar QPO ratios also exist in regions of stable shocks.

preprint2016arXiv

The NuSTAR X-ray Spectrum of Hercules X-1: A Radiation-Dominated Radiative Shock

We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of the Her X-1 35-day accretion disk precession period. This model allows us to estimate the accretion rate, the Comptonizing temperature of the radiating plasma, the radius of the magnetic polar cap, and the average scattering opacity parameters in the accretion column. This is in contrast to previous phenomenological models that characterized the shape of the X-ray spectrum but could not determine the physical parameters of the accretion flow. We describe the spectral fitting details and discuss the interpretation of the accretion flow physical parameters.

preprint2016arXiv

Time-Dependent Electron Acceleration in Blazar Transients: X-Ray Time Lags and Spectral Formation

Electromagnetic radiation from blazar jets often displays strong variability, extending from radio to $γ$-ray frequencies. In a few cases, this variability has been characterized using Fourier time lags, such as those detected in the X-rays from Mrk~421 using BeppoSAX. The lack of a theoretical framework to interpret the data has motivated us to develop a new model for the formation of the X-ray spectrum and the time lags in blazar jets based on a transport equation including terms describing stochastic Fermi acceleration, synchrotron losses, shock acceleration, adiabatic expansion, and spatial diffusion. We derive the exact solution for the Fourier transform of the electron distribution, and use it to compute the Fourier transform of the synchrotron radiation spectrum and the associated X-ray time lags. The same theoretical framework is also used to compute the peak flare X-ray spectrum, assuming that a steady-state electron distribution is achieved during the peak of the flare. The model parameters are constrained by comparing the theoretical predictions with the observational data for Mrk~421. The resulting integrated model yields, for the first time, a complete first-principles physical explanation for both the formation of the observed time lags and the shape of the peak flare X-ray spectrum. It also yields direct estimates of the strength of the shock and the stochastic MHD wave acceleration components in the Mrk~421 jet.

preprint2015arXiv

Fourier Analysis of Blazar Variability: Klein-Nishina Effects and the Jet Scattering Environment

The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing the PSDs and time lags produced via a combination of stochastic particle injection and emission via the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism used the Thomson cross section and simple $δ$-function approximations to model the synchrotron and Compton emissivities. Here we expand upon this work, using the full Compton cross section and detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed using the $δ$-function approximations and those computed using the accurate expressions, provided the observed photons are produced primarily by electrons with energies exceeding the lower limit of the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the cooling timescales of the electrons primarily responsible for the observed emission, and the associated time lags are related to the difference in electron cooling timescales between the two energy channels, as expected. If the electron cooling timescales can be determined from the observed time lags and/or the observed EC PSDs, then one could in principle use the method developed here to determine the energy of the external seed photon source for EC, which is an important unsolved problem in blazar physics.

preprint2015arXiv

Looking into the Theory of Pulsar Accretion: Cen X-3 and XTE J1946+274

This is an overview of pulsar accretion modeling. The physics of pulsar accretion, i.e., the process of plasma flow onto the neutron star surface, can be constrained from the spectral properties of the X-ray source. We discuss a new implementation of the physical continuum model developed by Becker and Wolff (2007, ApJ 654, 435). The model incorporates Comptonized blackbody, bremsstrahlung, and cyclotron emission. We discuss preliminary results of applying the new tool to the test cases of Suzaku data of Cen X-3 and XTE J1946+274. Cen X-3 is a persistent accreting pulsar with an O-star companion observed during a bright period. XTE J1946+274 is a transient accreting pulsar with a Be companion observed during a dim period. Both sources show spectra that are well described with an empirical Fermi Dirac cutoff power law model. We extend the spectral analysis by making the first steps towards a physical description of Cen X-3 and XTE J1946+274.

preprint2015arXiv

The Transient Accereting X-Ray Pulsar XTE J1946+274: Stability of the X-Ray Properties at Low Flux and Updated Orbital Solution

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

preprint2014arXiv

A Fourier Transformed Bremsstrahlung Flash Model for the Production of X-Ray Time Lags in Accreting Black Hole Sources

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.

preprint2014arXiv

Fourier Analysis of Blazar Variability

Blazars display strong variability on multiple timescales and in multiple radiation bands. Their variability is often characterized by power spectral densities (PSDs) and time lags plotted as functions of the Fourier frequency. We develop a new theoretical model based on the analysis of the electron transport (continuity) equation, carried out in the Fourier domain. The continuity equation includes electron cooling and escape, and a derivation of the emission properties includes light travel time effects associated with a radiating blob in a relativistic jet. The model successfully reproduces the general shapes of the observed PSDs and predicts specific PSD and time lag behaviors associated with variability in the synchrotron, synchrotron self-Compton (SSC), and external Compton (EC) emission components, from sub-mm to gamma-rays. We discuss applications to BL Lacertae objects and to flat-spectrum radio quasars (FSRQs), where there are hints that some of the predicted features have already been observed. We also find that FSRQs should have steeper PSD power-law indices than BL Lac objects at Fourier frequencies < 10^{-4} Hz, in qualitative agreement with previously reported observations by the Fermi Large Area Telescope.

preprint2012arXiv

No anti-correlation between cyclotron line energy and X-ray flux in 4U 0115+634

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3--50\,keV band. In agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7keV. Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10keV, the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the continuum.

preprint2012arXiv

TeV blazar variability: the firehose instability?

Recently observed minute timescale variability of blazar emission at TeV energies has imposed severe constraints on jet models and TeV emission mechanisms. We focus on a robust jet instability to explain this variability. As a consequence of the bulk outflow of the jet plasma, the pressure is likely to be anisotropic, with the parallel pressure $P_{||}$ in the forward jet direction exceeding the perpendicular pressure $P_{\perp}$. Under these circumstances, the jet is susceptible to the firehose instability, which can cause disruptions in the large scale jet structure and result in variability of the observed radiation. For a realistic range of parameters, we find that the growth timescale of the firehose instability is $\approx$ a few minutes, in good agreement with the observed TeV variability timescales for Mrk 501 (Albert et al. 2007) and PKS 2155-304 (Aharonian et al. 2007).

preprint2012arXiv

Understanding the Last Mile - Physics of the Accretion Column

Accreting X-ray pulsars are among the best observed objects of X-ray astronomy with a rich data set of observational phenomena in the spectral and timing domain. While the general picture for these sources is well established, the detailed physics behind the observed phenomena are often subject of debate. We present recent observational, theoretical and modeling results for the structure and dynamics of the accretion column in these systems. Our results indicate the presence of different accretion regimes and possible explanations for observed variations of spectral features with luminosity.

preprint2009arXiv

Particle Acceleration and the Formation of Relativistic Outflows in Viscous Accretion Disks with Shocks

In this Letter, we present a new self-consistent theory for the production of the relativistic outflows observed from radio-loud black hole candidates and active galaxies as a result of particle acceleration in hot, viscous accretion disks containing standing, centrifugally-supported isothermal shocks. This is the first work to obtain the structure of such disks for a relatively large value of the Shakura-Sunyaev viscosity parameter ($α=0.1$), and to consider the implications of the shock for the acceleration of relativistic particles in viscous disks. In our approach, the hydrodynamics and the particle acceleration are coupled and the solutions are obtained self-consistently based on a rigorous mathematical method. We find that particle acceleration in the vicinity of the shock can provide enough energy to power the observed relativistic jet in M87.

preprint2009arXiv

Particle Acceleration in Advection-Dominated Accretion Disks with Shocks: Green's Function Energy Distribution

The distribution function describing the acceleration of relativistic particles in an advection-dominated accretion disk is analyzed using a transport formalism that includes first-order Fermi acceleration, advection, spatial diffusion, and the escape of particles through the upper and lower surfaces of the disk. When a centrifugally-supported shock is present in the disk, the concentrated particle acceleration occurring in the vicinity of the shock channels a significant fraction of the binding energy of the accreting gas into a population of relativistic particles. These high-energy particles diffuse vertically through the disk and escape, carrying away both energy and entropy and allowing the remaining gas to accrete. The dynamical structure of the disk/shock system is computed self-consistently using a model previously developed by the authors that successfully accounts for the production of the observed relativistic outflows (jets) in M87 and \SgrA. This ensures that the rate at which energy is carried away from the disk by the escaping relativistic particles is equal to the drop in the radial energy flux at the shock location, as required for energy conservation. We investigate the influence of advection, diffusion, and acceleration on the particle distribution by computing the nonthermal Green's function, which displays a relatively flat power-law tail at high energies. We also obtain the energy distribution for the particles escaping from the disk, and we conclude by discussing the spectrum of the observable secondary radiation produced by the escaping particles.

preprint2009arXiv

Study of the accreting pulsar 4U 0115+634 with a bulk and thermal Comptonization model

Highly magnetized pulsars accreting matter in a binary system are bright sources in the X-ray band (0.1-100 keV). Despite the early comprehension of the basic emission mechanism, their spectral energy distribution is generally described by phenomenological or simplified models. We propose a study of the spectral emission from the high mass X-ray binary pulsar 4U 0115+634 by means of thermal and bulk Comptonization models based on the physical properties of such objects. For this purpose, we analyze the BeppoSAX data in the energy range 0.7-100 keV of the 1999 giant outburst, 12 days after the maximum. We model the spectral energy distribution of the system using a two-component continuum. At higher energy, above ~7 keV, the emission is due to thermal and bulk Comptonization of the seed photons produced by cyclotron cooling of the accretion column, and at lower energy, the emission is due to thermal Comptonization of a blackbody source in a diffuse halo close to the stellar surface. From the best fit parameters, we argue that the cyclotron emission is produced ~1.7 km above the stellar surface, and escapes from the column near its base, where the absorption features are generated by the interaction with the magnetic field in a surrounding halo. We find that in 4U 0115+634, the observed spectrum is dominated by reprocessed cyclotron radiation, whereas in other bright sources with stronger magnetic fields such as Her X-1, the spectrum is dominated by reprocessed bremsstrahlung.

preprint2003arXiv

On the integration of products of Whittaker functions with respect to the second index

Several new formulas are developed that enable the evaluation of a family of definite integrals containing the product of two Whittaker W-functions. The integration is performed with respect to the second index, and the first index is permitted to have any complex value, within certain restrictions required for convergence. The method utilizes complex contour integration along with various symmetry relations satisfied by the Whittaker functions. The new results derived in this paper are complementary to the previously known integrals of products of Whittaker functions, which generally treat integration with respect to either the first index or the primary argument. A physical application involving radiative transport is discussed.