Source author record

Indrek Vurm

Indrek Vurm 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

21works
5topics
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

21 published item(s)

preprint2022arXiv

The first nova eruption in a novalike variable: YZ Ret as seen in X-rays and gamma-rays

Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT and NuSTAR observations complemented by XMM-Newton X-ray grating spectroscopy to probe the physical parameters of the shocked ejecta and the nova-hosting white dwarf. The XMM-Newton observations revealed a super-soft X-ray emission which is dominated by emission lines of CV, CVI, NVI, NVII, and OVIII rather than a blackbody-like continuum, suggesting CO-composition of the white dwarf in a high-inclination binary system. Fermi/LAT detected YZ Ret for 15 days with the gamma-ray spectrum best described by a power law with an exponential cut-off at 1.9 +/-0.6 GeV. In stark contrast with theoretical predictions and in keeping with previous NuSTAR observations of Fermi-detected classical novae (V5855 Sgr and V906 Car), the 3.5-78 keV X-ray emission is found to be two orders of magnitude fainter than the GeV emission. The X-ray emission observed by NuSTAR is consistent with a single-temperature thermal plasma. We detect no non-thermal tail of the GeV emission expected to extend down to the NuSTAR band. NuSTAR observations continue to challenge theories of high-energy emission from shocks in novae.

preprint2021arXiv

High-Energy Neutrinos and Gamma-Rays from Non-Relativistic Shock-Powered Transients

Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discovery of time-correlated optical and gamma-ray emission in classical novae has revealed the important role of radiative shocks in powering these events, enabling an unprecedented view of the properties of ion acceleration, including its efficiency and energy spectrum, under similar physical conditions to shocks in extragalactic transients. Here we introduce a model for connecting the radiated optical fluence of non-relativistic transients to their maximal neutrino and gamma-ray fluence. We apply this technique to a wide range of extragalactic transient classes in order to place limits on their contributions to the cosmological high-energy gamma-ray and neutrino backgrounds. Based on a simple model for diffusive shock acceleration at radiative shocks, calibrated to novae, we demonstrate that several of the most luminous transients can accelerate protons up to energies $E_{\rm max} \gtrsim 10^{16}$ eV, sufficient to contribute to the IceCube astrophysical background. Furthermore, several of the considered sources$-$particularly hydrogen-poor supernovae$-$may serve as "hidden" gamma-ray sources due to the high gamma-ray opacity of their ejecta, evading constraints imposed by the non-blazar Fermi-LAT background. However, adopting an ion acceleration efficiency $\sim$ 0.3-1$\%$ motivated by nova observations, we find that currently known classes of non-relativistic, potentially shock-powered transients contribute at most a few percent of the total IceCube background.

preprint2020arXiv

Direct evidence for shock-powered optical emission in a nova

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

preprint2020arXiv

X-ray Spectroscopy of the gamma-ray Brightest Nova V906 Car (ASASSN-18fv)

Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.5-78 keV emission of V906 Car is much weaker than the simultaneously observed >100 MeV emission. No non-thermal X-ray emission is detected, and our deep limits imply that the gamma-rays are likely hadronic. After correcting for substantial absorption (N_H ~ 2 x 10^23 cm^-2), the thermal X-ray luminosity (from a 9 keV optically-thin plasma) is just ~2% of the gamma-ray luminosity. We consider possible explanations for the low thermal X-ray luminosity, including the X-rays being suppressed by corrugated, radiative shock fronts or the X-rays from the gamma-ray producing shock are hidden behind an even larger absorbing column (N_H >10^25 cm^-2). Adding XMM-Newton and Swift/XRT observations to our analysis, we find that the evolution of the intrinsic X-ray absorption requires the nova shell to be expelled 24 days after the outburst onset. The X-ray spectra show that the ejecta are enhanced in nitrogen and oxygen, and the nova occurred on the surface of a CO-type white dwarf. We see no indication of a distinct super-soft phase in the X-ray lightcurve, which, after considering the absorption effects, may point to a low mass of the white dwarf hosting the nova.

preprint2016arXiv

Polarization of gamma-ray bursts in the dissipative photosphere model

The MeV spectral peak of gamma-ray bursts (GRBs) is best explained as photospheric emission from a dissipative relativistic jet. The observed non-blackbody spectrum shows that sub-photospheric dissipation involves both thermal plasma heating and injection of nonthermal particles, which quickly cool through inverse Compton scattering and emission of synchrotron radiation. Synchrotron photons emitted around and above the photosphere are predicted to dominate the low-energy part of the GRB spectrum, starting from roughly a decade in energy below the MeV peak. We show that this leads to a unique polarization signature: a rise in GRB polarization toward lower energies. We compute the polarization degree of GRB radiation as a function of photon energy for a generic jet model, and show the predictions for GRBs 990123, 090902B and 110721A. The expected polarization is significant in the X-ray band, in particular for bursts similar to GRB 090902B. Radiation in the MeV peak (and at higher energies) is unpolarized as long as the jet is approximately uniform on angular scales $δθ\gtrsim Γ^{-1}$ where $Γ$ is the bulk Lorentz factor of the jet.

preprint2016arXiv

Radiative Transfer Models for Gamma-Ray Bursts

We present global radiative transfer models for heated relativistic jets. The simulations include all relevant radiative processes, starting deep in the opaque zone and following the evolution of radiation to and beyond the photosphere of the jet. The transfer models are compared with three gamma-ray bursts GRB 990123, GRB 090902B, and GRB 130427A, which have well-measured and different spectra. The models provide good fits to the observed spectra in all three cases, and we obtain estimates for the jet magnetization parameter $\varepsilon_{\rm B}$ and the Lorentz factor $Γ$. In the small sample of three bursts, $\varepsilon_{\rm B}$ varies between 0.01 and 0.05, and $Γ$ varies between 400 and 1200.

preprint2016arXiv

Shocks in nova outflows. II. Synchrotron radio emission

The discovery of GeV gamma-rays from classical novae indicates that shocks and relativistic particle acceleration are energetically key in these events. Further evidence for shocks comes from thermal keV X-ray emission and an early peak in the radio light curve on a timescale of months with a brightness temperature which is too high to result from freely expanding photo-ionized gas. Paper I developed a one dimensional model for the thermal emission from nova shocks. This work concluded that the shock-powered radio peak cannot be thermal if line cooling operates in the post-shock gas at the rate determined by collisional ionization equilibrium. Here we extend this calculation to include non-thermal synchrotron emission. Applying our model to three classical novae, we constrain the amplification of the magnetic field $ε_B$ and the efficiency $ε_e$ of accelerating relativistic electrons of characteristic Lorentz factor $γ\sim 100$. If the shocks are radiative (low velocity $v_{\rm sh} \lesssim 1000$ km s$^{-1}$) and cover a large solid angle of the nova outflow, as likely characterize those producing gamma-rays, then values of $ε_e \sim 0.01-0.1$ are required to achieve the peak radio brightness for $ε_B = 10^{-2}$. Such high efficiencies exclude secondary pairs from pion decay as the source of the radio-emitting particles, instead favoring the direct acceleration of electrons at the shock. If the radio-emitting shocks are instead adiabatic (high velocity), as likely characterize those responsible for the thermal X-rays, then much higher brightness temperatures are possible, allowing the radio-emitting shocks to cover a smaller outflow solid angle.

preprint2015arXiv

Gamma-ray novae as probes of relativistic particle acceleration at non-relativistic shocks

The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface.

preprint2015arXiv

Measuring Ambient Densities and Lorentz Factors of Gamma-Ray Bursts from GeV and Optical Observations

Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accompanied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash has been detected. Recent work suggests that the GeV+optical flash is emitted by the external blast wave from the explosion in a medium loaded with copious $e^\pm$ pairs. The full light curve of the flash is predicted by a first-principle radiative transfer simulation and can be tested against observations. Here we examine a sample of 7 bursts with best GeV+optical data and test the model. We find that the observed light curves are in agreement with the theoretical predictions and allow us to measure three parameters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, and the external density. With one possible exception of GRB 090510 (which is the only short burst in the sample) the ambient medium is consistent with a wind from a Wolf-Rayet progenitor. The wind density parameter $A=ρr^2$ varies in the sample around $10^{11}$g/cm. The initial Lorentz factor of the blast wave varies from 200 to 540 and correlates with the burst luminosity. Radiative efficiency of the prompt emission in the sample is between 0.1 and 0.8. For the two bursts with detected optical flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the external blast wave. Remarkably, the model reproduces the entire optical light curve of GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV data. The spectrum of GeV flashes is predicted to extend above 0.1 TeV, where they can be detected by ground-based Cherenkov telescopes.

preprint2015arXiv

Novae as Tevatrons: Prospects for CTA and IceCube

The discovery of novae as sources of ~GeV gamma-rays highlights the key role of shocks and relativistic particle acceleration in these transient systems. Although there is evidence for a spectral cut-off above energies ~1-100 GeV at particular epochs in some novae, the maximum particle energy achieved in these accelerators has remained an open question. The high densities of the nova ejecta (~10 orders of magnitude larger than in supernova remnants) render the gas far upstream of the shock neutral and shielded from ionizing radiation. The amplification of the magnetic field needed for diffusive shock acceleration requires ionized gas, thus confining the acceleration process to a narrow photo-ionized layer immediately ahead of the shock. Based on the growth rate in this layer of the hybrid non-resonant cosmic ray current-driven instability (considering also ion-neutral damping), we quantify the maximum particle energy, Emax, across the range of shock velocities and upstream densities of interest. We find values of Emax ~ 10 GeV - 10 TeV, which are broadly consistent with the inferred spectral cut-offs, but which could also in principle lead to emission extending to higher energies >100 GeV accessible to atmosphere Cherenkov telescopes, such as the planned Cherenkov Telescope Array (CTA). Detecting TeV neutrinos with IceCube in hadronic scenarios appears to be more challenging, although the prospects are improved if the shock power during the earliest, densest phases of the nova outburst is higher than is implied by the observed GeV light curves, due to downscattering of the gamma-rays by electrons within the ejecta. Novae provide ideal nearby laboratories to study magnetic field amplification and the onset of cosmic ray acceleration, because other time-dependent sources (e.g. radio supernovae) typically occur too distant to detect as gamma-ray sources.

preprint2014arXiv

On the origin of GeV emission in gamma-ray bursts

The most common progenitors of gamma-ray bursts (GRBs) are massive stars with strong stellar winds. We show that the GRB blast wave in the wind should emit a bright GeV flash. It is produced by inverse Compton cooling of the thermal plasma behind the forward shock wave. The main part of the flash is shaped by scattering of the prompt MeV radiation (emitted at smaller radii) which streams through the external blast wave. The inverse-Compton flash is bright due to the huge e+- enrichment of the external medium. At late times, the blast wave switches to normal synchrotron-self-Compton cooling. The mechanism is demonstrated by a detailed transfer simulation. The observed prompt MeV radiation is taken as an input of the simulation; we use GRB 080916C as an example. The result reproduces the GeV flash observed by the Fermi telescope. It explains the delayed onset, the steep rise, the peak flux, the time of the peak, the long smooth decline, and the spectral slope of GeV emission. The wind density required to reproduce all these features is typical of Wolf-Rayet stars. Our simulation predicts strong TeV emission 1 min after the burst trigger; then a cutoff in the observed high-energy spectrum is expected from absorption by extragalactic background light. In addition, a bright optical counterpart of the GeV flash is predicted for plausible values of the magnetic field; such a double (optical+GeV) flash has been observed in GRB 130427A.

preprint2014arXiv

Pair-dominated GeV-optical flash in GRB 130427A

We show that the light curve of the double GeV+optical flash in GRB 130427A is consistent with radiation from the blast wave in a wind-type medium with density parameter $A=ρr^2\sim 5\times 10^{10}$ g cm$^{-1}$. The peak of the flash is emitted by copious $e^\pm$ pairs created and heated in the blast wave; our first-principle calculation determines the pair-loading factor and temperature of the shocked plasma. Using detailed radiative transfer simulations we reconstruct the observed double flash. The optical flash is dominated by synchrotron emission from the thermal plasma behind the forward shock, and the GeV flash is produced via inverse Compton (IC) scattering by the same plasma. The seed photons for IC scattering are dominated by the prompt MeV radiation during the first tens of seconds, and by the optical to X-ray afterglow thereafter. IC cooling of the thermal plasma behind the forward shock reproduces all GeV data from a few seconds to $\sim 1$ day. We find that the blast wave Lorentz factor at the peak of the flash is $Γ\approx 200$, and the forward shock magnetization is $ε_B\sim 2\times 10^{-4}$. An additional source is required by the data in the optical and X-ray bands at times $>10^2$ s; we speculate that this additional source may be a long-lived reverse shock in the explosion ejecta.

preprint2014arXiv

Shocks in nova outflows. I. Thermal emission

Evidence for shocks in nova outflows include (1) multiple velocity components in the optical spectra; (2) keV X-ray emission weeks to months after the outburst; (3) early radio flare on timescales of months, in excess of that predicted from the freely expanding photo-ionized gas; and (4) ~ GeV gamma-rays. We present a 1D model for the shock interaction between the fast nova outflow and a dense external shell (DES) and its associated thermal X-ray, optical, and radio emission. The forward shock is radiative initially when the density of shocked gas is highest, at which times radio emission originates from the dense cooling layer immediately downstream of the shock. The radio light curve is characterized by sharper rises to maximum and later peak times at progressively lower frequencies, with a peak brightness temperature that is approximately independent of frequency. We apply our model to the recent gamma-ray classical nova V1324 Sco, obtaining an adequate fit to the early radio maximum for reasonable assumptions about the fast nova outflow and assuming the DES possesses a velocity ~1e3 km/s and mass ~ 2e-4 M_sun; the former is consistent with the velocities of narrow line absorption systems observed previously in nova spectra, while the total ejecta mass of the DES and fast outflow is consistent with that inferred independently by modeling the late radio peak. Rapid evolution of the early radio light curves require the DES possess a steep outer density profile, which may indicate that the onset of mass loss from the white dwarf was rapid, providing indirect evidence that the DES was expelled by the thermonuclear runaway event. Reprocessed X-rays from the shock absorbed by the DES at early times may contribute significantly to the optical/UV emission, which we speculate is responsible for the previously unexplained `plateaus' and secondary maxima in nova optical light curves.

preprint2013arXiv

Hot accretion flow in black hole binaries: a link connecting X-rays to the infrared

Multiwavelength observations of Galactic black hole transients have opened a new path to understanding the physics of the innermost parts of the accretion flows. While the processes giving rise to their X-ray continuum have been studied extensively, the emission in the optical and infrared (OIR) energy bands was less investigated and remains poorly understood. The standard accretion disc, which may contribute to the flux at these wavelengths, is not capable of explaining a number of observables: the infrared excesses, fast OIR variability and a complicated correlation with the X-rays. It was suggested that these energy bands are dominated by the jet emission, however, this scenario does not work in a number of cases. We propose here an alternative, namely that most of the OIR emission is produced by the extended hot accretion flow. In this scenario, the OIR bands are dominated by the synchrotron radiation from the non-thermal electrons. An additional contribution is expected from the outer irradiated part of the accretion disc heated by the X-rays. We discuss properties of the model and compare them to the data. We show that the hot-flow scenario is consistent with many of the observed spectral data, at the same time naturally explaining X-ray timing properties, fast OIR variability and its correlation with the X-rays.

preprint2013arXiv

Ionization Break-Out from Millisecond Pulsar Wind Nebulae: an X-ray Probe of the Origin of Superluminous Supernovae

Magnetic spin-down of a millisecond neutron star has been proposed as the power source of hydrogen-poor "superluminous" supernovae (SLSNe-I). However, producing an unambiguous test that can distinguish this model from alternatives, such as circumstellar interaction, has proven challenging. After the supernova explosion, the pulsar wind inflates a hot cavity behind the expanding stellar ejecta: the nascent millisecond pulsar wind nebula. Electron/positron pairs injected by the wind cool through inverse Compton scattering and synchrotron emission, producing a pair cascade and hard X-ray spectrum inside the nebula. These X-rays ionize the inner exposed side of the ejecta, driving an ionization front that propagates outwards with time. Under some conditions this front can breach the ejecta surface within months after the optical supernova peak, allowing ~0.1-1 keV photons to escape the nebula unattenuated with a characteristic luminosity L_X ~ 1e43-1e45 erg/s. This "ionization break-out" may explain the luminous X-ray emission observed from the transient SCP 06F, providing direct evidence that this SLSN was indeed engine-powered. Luminous break-out requires a low ejecta mass and that the spin-down time of the pulsar be comparable to the photon diffusion timescale at optical maximum, the latter condition similar to that required for a supernova with a high optical fluence. These relatively special requirements may explain why most SLSNe-I are not accompanied by detectable X-ray emission. Global asymmetry of the supernova ejecta increases the likelihood of an early break-out along the direction of lowest density. Atomic states with lower threshold energies are more readily ionized at earlier times near optical maximum, allowing UV break-out across a wider range of pulsar and ejecta properties than X-ray break-out, possibly contributing to the blue/UV colors of SLSNe-I.

preprint2012arXiv

On thermalization in gamma-ray burst jets and the peak energies of photospheric spectra

The low energy spectral slopes of the prompt emission of most gamma-ray bursts (GRBs) are difficult to reconcile with radiatively efficient optically thin emission models irrespective of the radiation mechanism. An alternative is to ascribe the radiation around the spectral peak to a thermalization process occurring well inside the Thomson photosphere. This quasi-thermal spectrum can evolve into the observed non-thermal shape by additional energy release at moderate to small Thomson optical depths, which can readily give rise to the hard spectral tail. The position of the spectral peak is determined by the temperature and Lorentz factor of the flow in the termalization zone, where the total number of photons carried by the jet is established. To reach thermalization, dissipation alone is not sufficient and photon generation requires an efficient emission/absorption process in addition to scattering. We perform a systematic study of all relevant photon production mechanisms searching for possible conditions in which thermalization can take place. We find that a significant fraction of the available energy should be dissipated at intermediate radii, $\sim 10^{10}$-- a few$\times 10^{11}$ cm and the flow there should be relatively slow: the bulk Lorentz factor could not exceed a few tens for all but the most luminous bursts with the highest $\Epk$-s. The least restrictive constraint for successful thermalization, $Γ\lesssim 20$, is obtained if synchrotron emission acts as the photon source. This requires, however, a non-thermal acceleration deep below the Thomson photosphere transferring a significant fraction of the flow energy to relativistic electrons with Lorentz factors between 10 and 100. Other processes require bulk flow Lorentz factors of order of a few for typical bursts. We examine the implications of these results to different GRB photospheric emission models.

preprint2011arXiv

A self-consistent hybrid Comptonization model for broad-band spectra of accreting supermassive black holes

The nature of the broad-band spectra of supermassive accreting black holes in active galactic nuclei (AGNs) is still unknown. The hard X-ray spectra of Seyferts as well as of Galactic stellar-mass black holes (GBHs) are well represented by thermal Comptonization, but the origin of the seed photons is less certain. The MeV tails observed in GBHs provide evidence in favour of non-thermal electron tails and it is possible that such electrons are also present in the X-ray emitting regions of AGNs. Using simulations with the kinetic code that self-consistently models electron and photon distributions, we find that the power-law-like X-ray spectra in AGNs can be explained in terms of the synchrotron self-Compton radiation of hybrid thermal/non-thermal electrons, similarly to the hard/low state of GBHs. Under a very broad range of parameters the model predicts a rather narrow distribution of photon spectral slopes consistent with that observed from LINERs and Seyferts at luminosities less than 3 per cent of the Eddington luminosity. The entire infrared to X-ray spectrum of these objects can be described in terms of our model, suggesting a tight correlation between the two energy bands. We show that the recently found correlation between slope and the Eddington ratio at higher luminosities can be described by the increasing fraction of disc photons in the emitting region, which may be associated with the decreasing inner radius of the optically thick accretion disc. The increasing flux of soft photons is also responsible for the transformation of the electron distribution from nearly thermal to almost completely non-thermal. The softer X-ray spectra observed in narrow-line Seyfert galaxies may correspond to non-thermal Comptonization of the disc photons, predicting that no cutoff should be observed up to MeV energies in these sources, similarly to the soft-state GBHs.

preprint2011arXiv

A synchrotron self-Compton -- disk reprocessing model for optical/X-ray correlation in black hole X-ray binaries

Physical picture of the emission mechanisms operating in the X-ray binaries was put under question by the simultaneous optical/X-ray observations with high time resolution. The light curves of the two energy bands appeared to be connected and the cross-correlation functions observed in three black hole binaries exhibited a complicated shape. They show a dip of the optical emission a few seconds before the X-ray peak and the optical flare just after the X-ray peak. This behavior could not be explained in terms of standard optical emission candidates (e.g., emission from the cold accretion disk or a jet). We propose a novel model, which explains the broadband optical to the X-ray spectra and the variability properties. We suggest that the optical emission consists of two components: synchrotron radiation from the non-thermal electrons in the hot accretion flow and the emission produced by reprocessing of the X-rays in the outer part of the accretion disk. The first component is anti-correlated with the X-rays, while the second one is correlated, but delayed and smeared relative to the X-rays. The interplay of the components explains the complex shape of the cross-correlation function, the features in the optical power spectral density as well as the time lags.

preprint2011arXiv

Gamma-ray bursts from magnetized collisionally-heated jets

Jets producing gamma-ray bursts (GRBs) are likely to carry a neutron component that drifts with respect to the proton component. The neutron-proton collisions strongly heat the jet and generate electron-positron pairs. We investigate radiation produced by this heating using a new numerical code. Our results confirm the recent claim that collisional heating generates the observed Band-type spectrum of GRBs. We extend the model to study the effects of magnetic fields on the emitted spectrum. We find that the spectrum peak remains near 1 MeV for the entire range of the magnetization parameter $0<ε_B<2$ that is explored in our simulations. The low-energy part of the spectrum softens with increasing $ε_B$, and a visible soft excess appears in the keV band. The high-energy part of the spectrum extends well above the GeV range and can contribute to the prompt emission observed by Fermi/LAT. Overall, the radiation spectrum created by the collisional mechanism appears to agree with observations, with no fine-tuning of parameters.

preprint2010arXiv

Theory of Compton scattering by anisotropic electrons

Compton scattering plays an important role in various astrophysical objects such as accreting black holes and neutron stars, pulsars, and relativistic jets, clusters of galaxies as well as the early Universe. In most of the calculations it is assumed that the electrons have isotropic angular distribution in some frame. However, there are situations where the anisotropy may be significant due to the bulk motions, or anisotropic cooling by synchrotron radiation, or anisotropic source of seed soft photons. We develop here an analytical theory of Compton scattering by anisotropic distribution of electrons that can simplify significantly the calculations. Assuming that the electron angular distribution can be represented by a second order polynomial over cosine of some angle (dipole and quadrupole anisotropy), we integrate the exact Klein-Nishina cross-section over the angles. Exact analytical and approximate formulae valid for any photon and electron energies are derived for the redistribution functions describing Compton scattering of photons with arbitrary angular distribution by anisotropic electrons. The analytical expressions for the corresponding photon scattering cross-section on such electrons as well as the mean energy of scattered photons, its dispersion and radiation pressure force are also derived. We applied the developed formalism to the accurate calculations of the thermal and kinematic Sunyaev-Zeldovich effects for arbitrary electron distributions.

preprint2009arXiv

Time-dependent modeling of radiative processes in hot magnetized plasmas

Numerical simulations of radiative processes in magnetized compact sources such as hot accretion disks around black holes, relativistic jets in active galaxies and gamma-ray bursts are complicated because the particle and photon distributions span many orders of magnitude in energy, they also strongly depend on each other, the radiative processes behave significantly differently depending on the energy regime, and finally due to the enormous difference in the time-scales of the processes. We have developed a novel computer code for the time-dependent simulations that overcomes these problems. The processes taken into account are Compton scattering, electron-positron pair production and annihilation, Coulomb scattering as well as synchrotron emission and absorption. No approximation has been made on the corresponding rates. For the first time, we solve coupled integro-differential kinetic equations for photons and electrons/positrons without any limitations on the photon and lepton energies. A numerical scheme is proposed to guarantee energy conservation when dealing with synchrotron processes in electron and photon equations. We apply the code to model non-thermal pair cascades in the blackbody radiation field, to study the synchrotron self-absorption as particle thermalization mechanism, and to simulate time evolution of stochastically heated pairs and corresponding synchrotron self-Compton photon spectra which might be responsible for the prompt emission of gamma-ray bursts. Good agreement with previous works is found in the parameter regimes where comparison is feasible, with the differences attributable to our improved treatment of the microphysics.