Engineering topics

Beloborodov, Andrei M.

Publications and source records attributed to Beloborodov, Andrei M..

16 records

DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Comptonization by reconnection plasmoids in black hole coronae II: Electron-ion plasma

ABSTRACT We perform 2D particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization 1 ≤ σ ≤ 3 (defined as the ratio of magnetic tension to plasma rest-mass energy density), which is expected in the coronae of accretion discs around black holes. We find that magnetic dissipation proceeds with inefficient energy exchange between the heated ions and the Compton-cooled electrons. As a result, most electrons are kept at a low temperature in Compton equilibrium with radiation, and so thermal Comptonization cannot reach photon energies $\sim 100\,$ keV observed from accreting black holes. Nevertheless, magnetic reconnection efficiently generates $\sim 100\,$ keV photons because of mildly relativistic bulk motions of the plasmoid chain formed in the reconnection layer. Comptonization by the plasmoid motions dominates the radiative output and controls the peak of the radiation spectrum Epk. We find Epk ∼ 40 keV for σ = 1 and Epk ∼ 100 keV for σ = 3. In addition to the X-ray peak around 100 keV, the simulations show a non-thermal MeV tail emitted by a non-thermal electron population generated near X-points of the reconnection layer. The results are consistent with the typical hard state of accreting black holes. In particular, we find that the spectrum of Cygnus X-1 is well explained by electron-ion reconnection with σ ∼ 3.

Astronomy & Astrophysics↗

DOE OSTI2022

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-transport Two-temperature GRMHD Simulations

Multiwavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nucleus transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole (BH). However, the formation, structure, and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semiopaque corona reach, respectively, $T_i$ ≳ 10 10 K and $T_e$ ≳ 5 × 10 8 K. The system produces radiation, powerful collimated jets, and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning BH by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Magnetar Bursts Due to Alfvén Wave Nonlinear Breakout

The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from 10 36 –10 43 erg s -1 . Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to magnetar bursts. Using fully 3D force-free electrodynamics simulations, we show that such magnetar bursts may be produced by Alfvén waves launched from localized magnetar quakes: a wave packet propagates to the outer magnetosphere, becomes nonlinear, and escapes the magnetosphere, forming an ultra-relativistic ejecta. The ejecta pushes open the magnetospheric field lines, creating current sheets behind it. Magnetic reconnection can happen at these current sheets, leading to plasma energization and X-ray emission. The angular size of the ejecta can be compact, ≲1 sr if the quake launching region is small, ≲0.01 sr at the stellar surface. We discuss implications for the FRBs and the coincident X-ray burst from SGR 1935+2154.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Relativistic Alfvén Waves Entering Charge-starvation in the Magnetospheres of Neutron Stars

Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current j A . We derive an analytic expression for the evolution of the wavevector k and the growth of j A . In the strongly sheared regime, j A may exceed the maximum current j0 that can be supported by the background e ± plasma. We investigate these charge-starved waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave is able to propagate successfully even when κ ≡ j A /j 0 $\ll$ 1. It sustains j A by compressing and advecting the plasma along the magnetic field lines with an increasing Lorentz factor, γ ≳ κ 1/2 . The simulations show how plasma instabilities lead to gradual dissipation of the wave energy. Our results suggest that an extremely high charge-starvation parameter κ ≳ 10 4 may be required in order for this mechanism to power the observed fast radio bursts (FRBs) from SGR 1935+2154. However, cosmological FRBs with much higher luminosities are unlikely to be a result of charge-starvation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Comptonization by reconnection plasmoids in black hole coronae I: Magnetically dominated pair plasma

We perform 2D particle-in-cell simulations of reconnection in magnetically dominated electron–positron plasmas subject to strong Compton cooling. We vary the magnetization σ $\gg$ 1, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of σ, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives ~25 percent of the dissipated reconnection power for σ = 10 and ~40 percent for σ = 40. We complement our particle-in-cell studies with a Monte Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for σ ~ 10, yields a spectrum consistent with the typical hard state of accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Coherent Electromagnetic Emission from Relativistic Magnetized Shocks

Relativistic magnetized shocks are a natural source of coherent emission, offering a plausible radiative mechanism for fast radio bursts (FRBs). We present first-principles 3D simulations that provide essential information for the FRB models based on shocks: the emission efficiency, spectrum, and polarization. The simulated shock propagates in an e ± plasma with magnetization σ > 1 . The measured fraction of shock energy converted to coherent radiation is ≃ 10 -3 σ -1 , and the energy-carrying wave number of the wave spectrum is ≃ 4ω c /c, where ω c is the upstream gyrofrequency. The ratio of the O-mode and X-mode energy fluxes emitted by the shock is ≃ 0.4σ -1 . The dominance of the X mode at σ >> 1 is particularly strong, approaching 100% in the spectral band around 2ω c . We also provide a detailed description of the emission mechanism for both X and O modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Fast Dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object's luminosity. We describe a dissipation mechanism that works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields B 1 and B 2 that propagate in opposite directions along the background magnetic field B 0 and collide. The collision forms a thin current sheet perpendicular to B 0 , which absorbs the incoming waves. The current sheet is sustained by an electric field E breaking the magnetohydrodynamic condition E < B and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude A = B 1 / B 0 = B 2 /B 0 propagating in a uniform B 0 . The mechanism is activated when A > 1/2. It dissipates a large fraction of the wave energy, † = (2A - 1)/A 2 , reaching 100% when A = 1. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate A of main interest, the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at A $\gg$ 1. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Synchrotron self-Compton radiation from magnetically dominated turbulent plasmas in relativistic jets

ABSTRACT Relativistic jets launched by rotating black holes are powerful emitters of non-thermal radiation. Extraction of the rotational energy via electromagnetic stresses produces magnetically dominated jets, which may become turbulent. Studies of magnetically dominated plasma turbulence from first principles show that most of the accelerated particles have small pitch angles, i.e. the particle velocity is nearly aligned with the local magnetic field. We examine synchrotron self-Compton radiation from anisotropic particles in the fast cooling regime. The small pitch angles reduce the synchrotron cooling rate and promote the role of inverse Compton (IC) cooling, which can occur in two different regimes. In the Thomson regime, both synchrotron and IC components have soft spectra, νFν ∝ ν1/2. In the Klein–Nishina regime, synchrotron radiation has a hard spectrum, typically νFν ∝ ν, over a broad range of frequencies. Our results have implications for the modelling of BL Lacertae objects (BL Lacs) and gamma-ray bursts (GRBs). BL Lacs produce soft synchrotron and IC spectra, as expected when Klein–Nishina effects are minor. The observed synchrotron and IC luminosities are typically comparable, which indicates a moderate anisotropy with pitch angles θ ≳ 0.1. Rare orphan gamma-ray flares may be produced when θ ≪ 0.1. The hard spectra of GRBs may be consistent with synchrotron radiation when the emitting particles are IC cooling in the Klein–Nishina regime, as expected for pitch angles θ ∼ 0.1. Blazar and GRB spectra can be explained by turbulent jets with a similar electron plasma magnetization parameter, σe ∼ 104, which for electron–proton plasmas corresponds to an overall magnetization σ = (me/mp)σe ∼ 10.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Kinetic Simulations of Radiative Magnetic Reconnection in the Coronae of Accreting Black Holes

We perform 2D and 3D particle-in-cell simulations of reconnection in magnetically dominated e ± plasmas subject to strong Compton cooling. Magnetic reconnection under such conditions can operate in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically dominated plasmoids with a small thermal pressure. Compton drag clears cavities inside plasmoids and also affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in nonradiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two distinct components: (1) A mildly relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the dissipated reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives ~20% of the dissipated reconnection power. It is populated by particles accelerated impulsively at X-points or "picked up" by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support magnetic reconnection as a mechanism powering hard X-ray emission from magnetically dominated regions of accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

Magnetars as Astrophysical Laboratories of Extreme Quantum Electrodynamics: The Case for a Compton Telescope

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.

Wadiasingh, Zorawar↗

NASA NTRS2014

NuSTAR Observations of X-Ray Bursts from the Magnetar 1E 1048.1-5937

We report the detection of eight bright X-ray bursts from the 6.5 s magnetar 1E 1048.1−5937, during a 2013 July observation campaign with the Nuclear Spectroscopic Telescope Array. We study the morphological and spectral properties of these bursts and their evolution with time. The bursts resulted in count rate increases by orders of magnitude, sometimes limited by the detector dead time, and showed blackbody spectra with kT is approx. 6-8 keV in the T90 duration of 1-4 s, similar to earlier bursts detected from the source. We find that the spectra during the tail of the bursts can be modeled with an absorbed blackbody with temperature decreasing with flux. The burst flux decays followed a power law of index 0.8-0.9. In the burst tail spectra, we detect a is approx. 13 keV emission feature, similar to those reported in previous bursts from this source as well as from other magnetars observed with the Rossi X-ray Timing Explorer.We explore possible origins of the spectral feature such as proton cyclotron emission, which implies a magnetic field strength of B is approx. 2×10(exp15) G in the emission region. However, the consistency of the energy of the feature in different objects requires further explanation.

MAGNETAR 1E 1048.1-5937↗

NASA NTRS2014

NuSTAR Observations of the Magnetar 1E 2259+586

We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data show pulsations above 20 keV for the first time and we report evidence that one of the pulses in the double-peaked pulse profile shifts position with energy. The pulsed fraction of the magnetar is shown to increase strongly with energy. Our spectral analysis reveals that the soft X-ray spectrum is well characterized by an absorbed double blackbody or blackbody plus power-law model in agreement with previous reports. Our new hard X-ray data, however, suggest that an additional component, such as a power law, is needed to describe the NuSTAR and Swift spectrum. We also fit the data with the recently developed coronal outflow model by Beloborodov for hard X-ray emission from magnetars. The outflow from a ring on the magnetar surface is statistically preferred over outflow from a polar cap.

X-rays: individual↗

NASA NTRS2014

Timing and Flux Evolution of the Galactic Center Magnetar SGR J1745-2900

We present the X-ray timing and spectral evolution of the Galactic Center magnetar SGR J1745-2900 (SGR*) for the first 4 months post-discovery using data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift observatories. Our timing analysis reveals a large increase in the magnetar spin-down rate by a factor of 2.6 plus or minus 0.07 over our data span. We further show that the change in spin evolution was likely coincident with a bright X-ray burst observed in 2013 June by Swift, and if so, there was no accompanying discontinuity in the frequency. We find that the source 3 to 10 kiloelectronvolt flux has declined monotonically by a factor of approximately 2 over an 80-day period post-outburst accompanied by an approximately 20 percent decrease in the source's blackbody temperature, although there is evidence for both flux and kiloteslas having leveled off. We argue that the torque variations are likely to be magnetospheric in nature and will dominate over any dynamical signatures of orbital motion around Sgr A*.

Stars: magnetic field↗

NASA NTRS2014

NuStar and Swift Observations of the Fast Rotating Magnetized White Dwarf AE Aquarii

AE Aquarii is a cataclysmic variable with the fastest known rotating magnetized white dwarf (P(sub spin) = 33.08 s). Compared to many intermediate polars, AE Aquarii shows a soft X-ray spectrum with a very low luminosity (LX (is) approximately 10(exp 31) erg per second). We have analyzed overlapping observations of this system with the NuSTAR and the Swift X-ray observatories in 2012 September. We find the 0.5-30 keV spectra to be well fitted by either an optically thin thermal plasma model with three temperatures of 0.75(+0.18 / −0.45), 2.29(+0.96 / −0.82), and 9.33 (+6.07 / −2.18) keV, or an optically thin thermal plasma model with two temperatures of 1.00 (+0.34 / −0.23) and 4.64 (+1.58 / −0.84) keV plus a power-law component with photon index of 2.50 (+0.17 / −0.23). The pulse profile in the 3-20 keV band is broad and approximately sinusoidal, with a pulsed fraction of 16.6% +/- 2.3%. We do not find any evidence for a previously reported sharp feature in the pulse profile.

cataclysmic variables↗

NASA NTRS2014

NuSTAR AND SWIFT Observations of the Fast Rotating Magnetized White Dwarf AE Aquarii

AE Aquarii is a cataclysmic variable with the fastest known rotating magnetized white dwarf (P(sub spin) = 33.08 s). Compared to many intermediate polars, AE Aquarii shows a soft X-ray spectrum with a very low luminosity (LX (is) approximately 10(exp 31) erg per second). We have analyzed overlapping observations of this system with the NuSTAR and the Swift X-ray observatories in 2012 September. We find the 0.5-30 keV spectra to be well fitted by either an optically thin thermal plasma model with three temperatures of 0.75(+0.18 / -0.45), 2.29(+0.96 / -0.82), and 9.33 (+6.07 / -2.18) keV, or an optically thin thermal plasma model with two temperatures of 1.00 (+0.34 / -0.23) and 4.64 (+1.58 / -0.84) keV plus a power-law component with photon index of 2.50 (+0.17 / -0.23). The pulse profile in the 3-20 keV band is broad and approximately sinusoidal, with a pulsed fraction of 16.6% +/- 2.3%. We do not find any evidence for a previously reported sharp feature in the pulse profile.

power-law component↗