Researcher profile

Mitchell C. Begelman

Mitchell C. Begelman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2022arXiv

A radio, optical, UV and X-ray view of the enigmatic changing look Active Galactic Nucleus 1ES~1927+654 from its pre- to post-flare states

The nearby type-II AGN 1ES1927+654 went through a violent changing-look (CL) event beginning December 2017 during which the optical and UV fluxes increased by four magnitudes over a few months, and broad emission lines newly appeared in the optical/UV. By July 2018 the X-ray coronal emission had completely vanished, only to reappear a few months later. In this work we report the evolution of the radio, optical, UV and X-rays from the pre-flare state through mid-2021 with new and archival data from the Very Long Baseline Array (VLBA), the European VLBI Network, the Very Large Array (VLA), the Telescopio Nazionale Galileo (TNG), Gran Telescopio Canarias (GTC), The Neil Gehrels Swift observatory and XMM-Newton. The main results from our work are: (1) The source has returned to its pre-CL state in optical, UV, and X-ray; the disk-corona relation has been re-established as has been in the pre-CL state, with an $α_{\rm OX}\sim 1.02$. The optical spectra are dominated by narrow emission lines. (2) The UV light curve follows a shallower slope of $\propto t^{-0.91\pm 0.04}$ compared to that predicted by a tidal disruption event. We conjecture that a magnetic flux inversion event is the possible cause for this enigmatic event. (3) The compact radio emission which we tracked in the pre-CL (2014), during CL (2018) and post-CL(2021) at spatial scales $<1$ pc was at its lowest level during the changing look event in 2018, nearly contemporaneous with a low $2-10$ keV emission. The radio to X-ray ratio of the compact source $L_{\rm Radio}/L_{\rm X-ray}\sim 10^{-5.5}$, follows the Gudel-Benz relation, typically found in coronally active stars, and several AGN. (4) We do not detect any presence of nascent jets at the spatial scales of $\sim 5-10$ pc.

preprint2022arXiv

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma. Generalizing recent treatments by Alves et al. (2018) and Davelaar et al. (2020), we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. (2018) is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

preprint2022arXiv

Magnetization of Relativistic Current-Carrying Jets with Radial Velocity Shear

Astrophysical jets, launched from the immediate vicinity of accreting black holes, carry away large amounts of power in a form of bulk kinetic energy of jet particles and electromagnetic flux. Here we consider a simple analytical model for relativistic jets at larger distances from their launching sites, assuming a cylindrical axisymmetric geometry with a radial velocity shear, and purely toroidal magnetic field. We argue that, as long as the jet plasma is in magnetohydrostatic equilibrium, such outflows tend to be particle dominated, i.e. the ratio of the electromagnetic to particle energy flux, integrated over the jet cross-sectional area, is typically below unity, $σ< 1$. At the same time, for particular magnetic and radial velocity profiles, magnetic pressure may still dominate over particle pressure for certain ranges of the jet radius, i.e. the local jet plasma parameter $β_{pl} < 1$, and this may be relevant in the context of particle acceleration and production of high-energy emission in such systems. The jet magnetization parameter can be elevated up to the modest values $σ\lesssim \mathcal{O}(10)$ only in the case of extreme gradients or discontinuities in the gaseous pressure, and a significantly suppressed velocity shear. Such configurations, which consist of a narrow, unmagnetized jet spine surrounded by an extended, force-free layer, may require an additional poloidal field component to stabilize them against current-driven oscillations, but even this will not elevate substantially their $σ$ parameter.

preprint2022arXiv

Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc

Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of $10-100$ for $\approx 60$ minutes. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from PIC simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies $γ_e \lesssim 10^6$ and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simpified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard ($\lesssim 100$ keV) X-ray energies.

preprint2022arXiv

What really makes an accretion disc MAD

Magnetically arrested accretion discs (MADs) around black holes (BH) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here we propose instead that it is mainly a dynamically important {\it toroidal} field -- the result of dynamo action triggered by the significant but still relatively weak vertical field -- that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotatonal instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (strongly magnetized but not &#34;arrested&#34;) discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.

preprint2020arXiv

First-principles demonstration of diffusive-advective particle acceleration in kinetic simulations of relativistic plasma turbulence

Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio kinetic simulations of driven magnetized turbulence in relativistic pair plasma. By statistically analyzing the motion of tracked particles, we demonstrate the diffusive nature of NTPA and measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $γm_e c^2$. We find that $D(γ)$ scales as $γ^2$ in the high-energy nonthermal tail, in line with 2nd-order Fermi acceleration theory, but has a much weaker scaling at lower energies. We also find that $A$ is not negligible and reduces NTPA by tending to pull particles towards the peak of the particle energy distribution. This study provides strong support for the FP picture of turbulent NTPA, thereby enhancing our understanding of space and astrophysical plasmas.

preprint2020arXiv

Kinetic turbulence in shining pair plasma: intermittent beaming and thermalization by radiative cooling

High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization ($0.04 \lesssim σ\lesssim 11$). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of nonthermal energetic particles (extending the tail by a factor of $\sim 2$). This indicates that nonthermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g., second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, nonthermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*).

preprint2020arXiv

Magnetic flux inversion in a peculiar changing look AGN

We argue that the changing-look event in the active galactic nucleus 1ES 1927+654, followed by a dip of 3 orders of magnitude in the X-ray luminosity, is controlled by a change in the accretion rate and an inversion of magnetic flux in a magnetically arrested disk (MAD). Before the changing-look event, strong magnetic flux on the black hole powers X-ray emission via the Blandford-Znajek process while the UV emission is produced by a radiatively inefficient magnetized disk. An advection event, bringing flux of the opposite polarity, propagates inward leading, first, to a rise in the UV/optical luminosity and, then, to a dip in the X-ray luminosity when it reaches the black hole. We estimate the timescale for magnetic flux advection and find that the observed timescale between the beginning of the changing-look event and the minimum in the X-ray luminosity, $\approx200$ days, is in agreement with the time needed to cancel the magnetic flux in a MAD extending to $\approx180\:r_g$. Although flux inversion events might be rare due to the large ratio of flux-to-mass that is needed, we argue that AGN showing an unusually high ratio of X-ray to UV luminosity are prime candidates for such events. We also suggest that similar events may lead to jet interruptions in radio-loud objects.

preprint2020arXiv

Strongly magnetized accretion disks: structure and accretion from global magnetohydrodynamic simulations

We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r \approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $β_0 = 1000,\, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks&#39; vertical structure in all three cases, with accretion occurring at $z/R\approx 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $β_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically &#34;elevated&#34; state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).

preprint2008arXiv

Hot self-similar relativistic MHD flows

We consider axisymmetric relativistic jets with a toroidal magnetic field and an ultrarelativistic equation of state, with the goal of studying the lateral structure of jets whose pressure is matched to the pressure of the medium through which they propagate. We find all self-similar steady-state solutions of the relativistic MHD equations for this setup. One of the solutions is the case of a parabolic jet being accelerated by the pressure gradient as it propagates through a medium with pressure declining as p(z)\propto z^{-2}. As the jet material expands due to internal pressure gradients, it runs into the ambient medium resulting in a pile-up of material along the jet boundary, while the magnetic field acts to produce a magnetic pinch along the axis of the jet. Such jets can be in a lateral pressure equilibrium only if their opening angle θ_j at distance z is smaller than about 1/γ, where γis the characteristic bulk Lorentz-factor at this distance; otherwise, different parts of the jet cannot maintain causal contact. We construct maps of optically thin synchrotron emission from our self-similar models. We suggest that the boundary pile-up may be the reason for the limb-brightening of the sub-parsec jet of M87. We find that if the synchrotron emissivity falls with the distance from the jet axis, the polarization fraction rises toward the edge, as seen in 3C273 and Mkn501. Projection effects and the emissivity pattern of the jet have a strong effect on the observed polarization signal, so the interpretation of the polarization data in terms of the geometry of magnetic fields is rather uncertain. For example, jets with toroidal magnetic fields display the `spine-sheath&#39; polarization angle pattern seen in some BL Lac objects.

preprint2007arXiv

Temporal and Angular Properties of GRB Jets Emerging from Massive Stars

We study the long-term evolution of relativistic jets in collapsars and examine the effects of viewing angle on the subsequent gamma ray bursts. We carry out a series of high-resolution simulations of a jet propagating through a stellar envelope in 2D cylindrical coordinates using the FLASH relativistic hydrodynamics module. For the first time, simulations are carried out using an adaptive mesh that allows for a large dynamic range inside the star while still being efficient enough to follow the evolution of the jet long after it breaks out from the star. Our simulations allow us to single out three phases in the jet evolution: a precursor phase in which relativistic material turbulently shed from the head of the jet first emerges from the star, a shocked jet phase where a fully shocked jet of material is emerging, and an unshocked jet phase where the jet consists of a free-streaming, unshocked core surrounded by a thin boundary layer of shocked jet material. The appearance of these phases will be different to observers at different angles. The precursor has a wide opening angle and would be visible far off axis. The shocked phase has a relatively narrow opening angle that is constant in time. During the unshocked jet phase the opening angle increases logarithmically with time. As a consequence, some observers see prolonged dead times of emission even for constant properties of the jet injected in the stellar core. We also present an analytic model that is able to reproduce the overall properties of the jet and its evolution. We finally discuss the observational implications of our results, emphasizing the possible ways to test progenitor models through the effects of jet propagation in the star. In an appendix, we present 1D and 2D tests of the FLASH relativistic hydrodynamics module.

preprint2006arXiv

The physics of GRB jets and their interaction with the progenitor star

It is now generally accepted that long gamma-ray bursts are associated with the final evolutionary stages of massive stars. As a consequence, their jets must propagate through the stellar progenitor and break out on their surface, before they can reach the photospheric radius and produce the gamma-ray photons. We investigate the role of the progenitor star in shaping the jet properties. We show that even a jet powered by a steady engine can develop a rich phenomenology at the stellar surface. We present special-relativistic simulations and compare the results to analytic considerations. We show that the jet is complex in the time as well as in the angular domain, so that observers located along different lines of sight detect significantly different bursts.

preprint2002arXiv

Heating, conduction and minimum temperatures in cooling flows

There is mounting observational evidence from Chandra for strong interaction between keV gas and AGN in cooling flows. It is now widely accepted that the temperatures of cluster cores are maintained at a level of 1 keV and that the mass deposition rates are lower than earlier ROSAT/Einstein values. Recent theoretical results suggest that thermal conduction can be very efficient even in magnetized plasmas. Motivated by these discoveries, we consider a ``double heating model&#39;&#39; which incorporates the effects of simultaneous heating by both the central AGN and thermal conduction from the hot outer layers of clusters. Using hydrodynamical simulations, we demonstrate that there exists a family of solutions that does not suffer from the cooling catastrophe. In these cases, clusters relax to a stable final state, which is characterized by minimum temperatures of order 1 keV and density and temperature profiles consistent with observations. Moreover, the accretion rates are much reduced, thereby reducing the need for excessive mass deposition rates required by the standard cooling flow models.

preprint1998arXiv

RXTE Observation of Cygnus X--1: III. Implications for Compton Corona and ADAF Models

We have recently shown that a `sphere+disk&#39; geometry Compton corona model provides a good description of Rossi X-ray Timing Explorer (RXTE) observations of the hard/low state of Cygnus X--1. Separately, we have analyzed the temporal data provided by RXTE. In this paper we consider the implications of this timing analysis for our best-fit `sphere+disk&#39; Comptonization models. We focus our attention on the observed Fourier frequency-dependent time delays between hard and soft photons. We consider whether the observed time delays are: created in the disk but are merely reprocessed by the corona; created by differences between the hard and soft photon diffusion times in coronae with extremely large radii; or are due to `propagation&#39; of disturbances through the corona. We find that the time delays are most likely created directly within the corona; however, it is currently uncertain which specific model is the most likely explanation. Models that posit a large coronal radius [or equivalently, a large Advection Dominated Accretion Flow (ADAF) region] do not fully address all the details of the observed spectrum. The Compton corona models that do address the full spectrum do not contain dynamical information. We show, however, that simple phenomenological propagation models for the observed time delays for these latter models imply extremely slow characteristic propagation speeds within the coronal region.

preprint1998arXiv

RXTE Observation of Cygnus X-1: II. Timing Analysis

We present timing analysis for a Rossi X-ray Timing Explorer observation of Cygnus X-1 in its hard/low state. This was the first RXTE observation of Cyg X-1 taken after it transited back to this state from its soft/high state. RXTE&#39;s large effective area, superior timing capabilities, and ability to obtain long, uninterrupted observations have allowed us to obtain measurements of the power spectral density (PSD), coherence function, and Fourier time lags to a decade lower in frequency and half a decade higher in frequency than typically was achieved with previous instruments. Notable aspects of our observations include a weak 0.005 Hz feature in the PSD coincident with a coherence recovery; a `hardening&#39; of the high-frequency PSD with increasing energy; a broad frequency range measurement of the coherence function, revealing rollovers from unity coherence at both low and high frequency; and an accurate determination of the Fourier time lags over two and a half decades in frequency. As has been noted in previous similar observations, the time delay is approximately proportional to f^(-0.7), and at a fixed Fourier frequency the time delay of the hard X-rays compared to the softest energy channel tends to increase logarithmically with energy. Curiously, the 0.01-0.2 Hz coherence between the highest and lowest energy bands is actually slightly greater than the coherence between the second highest and lowest energy bands. We carefully describe all of the analysis techniques used in this paper, and we make comparisons of the data to general theoretical expectations. In a companion paper, we make specific comparisons to a Compton corona model that we have successfully used to describe the energy spectral data from this observation.