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 $\gamma$-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 $\gamma$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $\gamma$-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 $\gamma$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques.

We report the first detection of hard (>10 keV) X-ray emission simultaneous with gamma rays in a nova eruption. Observations of the nova V5855 Sgr carried out with the NuSTAR satellite on Day 12 of the eruption revealed faint, highly absorbed thermal X-rays. The extreme equivalent hydrogen column density towards the X-ray emitting region (~3 x 10$^{24}$ cm$^{-2}$) indicates that the shock producing the X-rays was deeply embedded within the nova ejecta. The slope of the X-ray spectrum favors a thermal origin for the bulk of the emission, and the constraints of the temperature in the shocked region suggest a shock velocity compatible with the ejecta velocities inferred from optical spectroscopy. While we do not claim the detection of non-thermal X-rays, the data do not allow us to rule out an additional, fainter component dominating at energy above 20 keV, for which we obtained upper limits. The inferred luminosity of the thermal X-rays is too low to be consistent with the gamma-ray luminosities if both are powered by the same shock under standard assumptions regarding the efficiency of non-thermal particle acceleration and the temperature distribution of the shocked gas.

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.

On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A.

Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large Area Telescope have shown that many classical novae are accompanied by gigaelectronvolt gamma-ray emission. This emission likely originates from strong shocks, providing new insights into the properties of nova outflows and allowing them to be used as laboratories to study the unknown efficiency of particle acceleration in shocks. Here we report gamma-ray and optical observations of the Milky Way nova ASASSN-16ma, which is among the brightest novae ever detected in gamma-rays. The gamma-ray and optical light curves show a remarkable correlation, implying that the majority of the optical light comes from reprocessed emission from shocks rather than the white dwarf. The ratio of gamma-ray to optical flux in ASASSN-16ma directly constrains the acceleration efficiency of non-thermal particles to be ~0.005, favouring hadronic models for the gamma-ray emission. The need to accelerate particles up to energies exceeding 100 gigaelectronvolts provides compelling evidence for magnetic field amplification in the shocks.

Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair production and annihilation. We simulate shocks with properties relevant for GRBs and study the steady-state solutions, which are accurate deep below the jet photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upwards from the typical upstream photon energy. Strong shocks (for which the downstream pressure is much larger than the upstream pressure) have rising $\nu F_\nu$ shock spectra. The spectrum extends up to $\epsilon_{max} \equiv E_{max}/m_e c^2 \sim v^2$ for non-relativistic shocks, where $m_e$ is the electron rest mass and $v$ is the relative speed between the upstream and downstream in units of the speed of light $c$. For mildly relativistic shocks the power law softens at $\epsilon \gtrsim 10^{-1}$ due to Klein-Nishina effects, and shocks with $v\gamma \gtrsim 1$, where $\gamma \equiv (1-v^2)^{-1/2}$, produce electron-positron pairs. As an example, a strong shock with $v\gamma = 3$ and a photon-to-proton ratio of $n_\gamma/n_p = 2 \times 10^5$ has a peak pair-to-proton ratio of $Z_\pm \approx 225$. The main effect of pairs in a steady-state shock is to decrease its spatial width by a factor of $\sim Z_\pm$. The post-shock spectrum thermalizes in the downstream. In absence of emission and absorption processes, kinetic equilibrium at temperature $\theta_d \equiv kT_d/m_e c^2 \approx \epsilon_d/3$ is reached at an optical depth of $\tau \gg \theta_d^{-1}$ behind the shock, where $\epsilon_d$ is the average downstream photon energy.

A gamma-ray burst (GRB) jet running into an external medium is expected to generate luminous GeV-TeV emission lasting from minutes to several hours. The high-energy emission results from inverse Compton upscattering of prompt and afterglow photons by shock-heated {\it thermal} plasma. At its peak the high-energy radiation carries a significant fraction of the power dissipated at the forward shock. We discuss in detail the expected TeV luminosity, using a robust "minimal" emission model. Then, using the statistical properties of the GRB population (luminosity function, redshift distribution, afterglow energy) we simulate the expected detection rates of GRBs by current and upcoming atmospheric Cherenkov instruments. We find that GRBs exploding into a low-density interstellar medium must produce TeV emission that would have already been detected by the currently operating Cherenkov telescopes. The absence of detections is consistent with explosions into a dense wind of the GRB progenitor. If the typical environment of long GRBs is a Wolf-Rayet progenitor wind, as suggested by the recent analysis of Fermi LAT data, then the upcoming Cherenkov Telescope Array (CTA) should be able to detect a fraction of GRBs that trigger the space-borne detectors. Since absorption by the extragalactic background light limits the detectability above 0.1 TeV for all but the most nearby bursts ($z\lesssim 1$), the reduced energy threshold of CTA is the key improvement over current instruments, which should boost the detection rate by at least a factor of 3.

Multiwavelength radiation from relativistic particles accelerated at shocks in novae and other astrophysical sources carries a wealth of information about the outflow properties and the microphysical processes at work near the shocks. The observation of GeV gamma-rays from novae by Fermi/LAT demonstrates that the shocks in these systems can accelerate particles to energies of at least $\sim 10$ GeV. The low-energy extension of the same non-thermal particle distribution inevitably gives rise to emission extending into the X-ray band. Above $\gtrsim 10$ keV this radiation can escape the system without significant absorption/attenuation, and can potentially be detected by NuSTAR. We present theoretical models for hard X-ray and gamma-ray emission from radiative shocks in both leptonic and hadronic scenarios, accounting for the rapid evolution of the downstream properties due to the fast cooling of thermal plasma. Due to strong Coulomb cooling of the mildly relativistic electrons nominally responsible for producing hard X-ray emission, only a fraction of $10^{-4} - 10^{-3}$ of the gamma-ray luminosity is radiated in the NuSTAR band; nevertheless, this emission could be detectable simultaneous with the LAT emission in bright gamma-ray novae with a $\sim 50$ ks exposure. Our work demonstrates how combined hard X-ray and gamma-ray observations can be used to constrain properties of the nova outflow (velocity, density and mass outflow rate) and particle acceleration at the shock. A very low X-ray to gamma-ray luminosity ratio ($L_{\rm X}/L_{\gamma} \lesssim 5 \times 10^{-4}$) would disfavor leptonic models for the gamma-ray emission. Our model can also be applied to other astrophysical environments with radiative shocks, including Type IIn supernovae and colliding winds in massive star binaries.

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 $\delta\theta \gtrsim \Gamma^{-1}$ where $\Gamma$ is the bulk Lorentz factor of the jet.

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 $\epsilon_B$ and the efficiency $\epsilon_e$ of accelerating relativistic electrons of characteristic Lorentz factor $\gamma \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 $\epsilon_e \sim 0.01-0.1$ are required to achieve the peak radio brightness for $\epsilon_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.

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 $\Gamma$. In the small sample of three bursts, $\varepsilon_{\rm B}$ varies between 0.01 and 0.05, and $\Gamma$ varies between 400 and 1200.

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.

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=\rho 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.

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.

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=\rho 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 $\Gamma\approx 200$, and the forward shock magnetization is $\epsilon_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.

This paper introduces a kinetic code that simulates gamma-ray burst (GRB) afterglow emission from the external forward shock and presents examples of some of its applications. One interesting research topic discussed in the paper is the high-energy radiation produced by Compton scattering of the prompt GRB photons against the shock-accelerated electrons. The difference between the forward shock emission in a wind-type and a constant-density medium is also studied, and the emission due to Maxwellian electron injection is compared to the case with pure power-law electrons. The code calculates the time-evolving photon and electron distributions in the emission region by solving the relativistic kinetic equations for each particle species. For the first time, the full relativistic equations for synchrotron emission/absorption, Compton scattering, and pair production/annihilation were applied to model the forward shock emission. The synchrotron self-absorption thermalization mechanism, which shapes the low-energy end of the electron distribution, was also included in the electron equation. The simulation results indicate that inverse Compton scattering of the prompt GRB photons can produce a luminous TeV emission component, even when pair production in the emission region is taken into account. This very high-energy radiation may be observable in low-redshift GRBs. The test simulations also show that the low-energy end of a pure power-law distribution of electrons can thermalize owing to synchrotron self-absorption in a wind-type environment, but without an observable impact on the radiation spectrum. Moreover, a flattening in the forward shock X-ray light curve may be expected when the electron injection function is assumed to be purely Maxwellian instead of a power law.

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.

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.

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.

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.

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, $\Gamma\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.

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.

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.

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<\epsilon_B<2$ that is explored in our simulations. The low-energy part of the spectrum softens with increasing $\epsilon_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.

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.

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.

The origin of dramatically different electron distributions responsible for Comptonization in black hole X-ray binaries (BHBs) in their various states is discussed. We solve the coupled kinetic equations for photons and electrons without approximations on the relevant cross-sections accounting for Compton scattering, synchrotron radiation, and Coulomb collisions. In the absence of external soft photons, the electrons are efficiently thermalized by synchrotron self-absorption and Coulomb scattering even for pure nonthermal electron injection. The resulting quasi-thermal synchrotron self-Compton spectra have very stable slopes and electron temperatures similar to the hard states of BHBs. The hard spectral slopes observed in the X-rays, the cutoff at 100 keV, and the MeV tail together require low magnetic fields ruling out the magnetic dissipation mechanism. The motion of the accretion disk toward the black hole results in larger Compton cooling and lower equilibrium electron temperature. Our self-consistent simulations show that in this case both electron and photon distributions attain a power-law-dominated shape similar to what is observed in the soft state. The electron distribution in the Cyg X-1 soft state might require a strong magnetic field, being consistent with the magnetically dominated corona.

We present detailed calculations of the electron thermalization by synchrotron self-absorption accounting for cooling by Compton scattering. For the first time, we solve coupled kinetic equations for electrons and photons without any approximations on the relevant cross-sections and compute self-consistently the resulting electron and photon distributions. The results presented in this contribution may be applied to the magnetized coronae around Galactic black hole accretion disks.