Non-thermal Electron Population Research Articles

Sunyaev-Zeldovich signatures from non-thermal, relativistic electrons using CMB maps

Relativistically energetic, non-thermal population of electrons can give rise to unique CMB spectral distortion signatures which can be significantly different from thermal Sunyaev-Zeldovich signal or y-distortion. These signatures depend upon the spectrum of non-thermal electrons, therefore, a detection can inform us about the existence and abundance of non-thermal electrons in our universe. Using public CMB maps and data, we derive upper limits on non-thermal y-parameter for a relativistic, power-law electron distribution. With future CMB experiments, we may be in a position to detect or put significantly tighter constraints on these signals which can affect our understanding of non-thermal electron distributions in our universe.

Impact of Electron Precipitation on Brown Dwarf Atmospheres and the Missing Auroral H3+ Emission

Recent observations have demonstrated that very low-mass stars and brown dwarfs are capable of sustaining strong magnetic fields despite their cool and neutral atmospheres. These kilogauss field strengths are inferred based on strong, highly circularly polarized gigahertz radio emission, a consequence of the electron cyclotron maser instability. Crucially, these observations imply the existence of energetic nonthermal electron populations, associated with strong current systems, as are found in the auroral regions of the magnetized planets of the solar system. Intense auroral electron precipitation will lead to electron collisions with the H2 gas that should generate the ion H3+ . With this motivation, we targeted a sample of ultracool dwarfs, known to exhibit signatures associated with aurorae, in search of the K-band emission features of H3+ using the Keck telescopes on Maunakea. From our sample of nine objects, we found no clear indication of H3+ emission features in our low-to-medium-resolution spectra (R ∼ 3600). We also modeled the impact of an auroral electron beam on a brown dwarf atmosphere, determining the depth at which energetic beams deposit their energy and drive particle impact ionization. We find that the H3+ nondetections can be explained by electron beams of typical energies ≳2–10 keV, which penetrate deeply enough that any H3+ produced is chemically destroyed before radiating energy through its infrared transitions. Strong electron beams could further explain the lack of UV auroral detections and suggest that most or nearly all of the precipitating auroral energy must ultimately emerge as thermal emissions deep in brown dwarf atmospheres.

Synchrotron Polarization Signatures of Surface Waves in Supermassive Black Hole Jets

Supermassive black holes in active galactic nuclei are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of nonthermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The telltale polarization signatures of the waves could be observable by future very long baseline interferometric observations, e.g., the next-generation Event Horizon Telescope.

Electrostatic solitary waves in a bi-ion plasma with two suprathermal electron populations – application to Saturn’s magnetosphere

ABSTRACT Non-thermal particle distributions characterized by a high-energy tail are ubiquitous in space plasmas. They are usually described by a kappa distribution function, that has been shown to be an excellent fit in most real circumstances. Among other space missions, Cassini and Voyager have both recorded evidence of a coexistence of non-thermal electron populations (with different characteristics) in Saturn’s magnetosphere, and subsequent studies showed that these are well-described by using different tailor-fit realizations of the (parametrized) kappa distribution. Motivated by these observations we have formulated a multifluid plasma model incorporating two types of (positive) ions and two distinct kappa-distributed electron populations, in order to study electrostatic solitary waves (ESWs) in Saturn’s magnetosphere from first principles. Our analysis reveals that the spectral index (in fact, the κ parameter value related to the cold electron population mainly) is vital in explaining the difference among different types of non-linear structures. A comparison with spacecraft observations suggests that our theoretical model provides an efficient framework for the interpretation of ESW observations in Saturn’s magnetosphere. Our qualitative predictions may also apply to other planetary magnetospheres, where a similar multispecies plasma composition may be present.

Understanding the nature of the optical emission in gamma-ray bursts: analysis from TAROT, COATLI, and RATIR observations

ABSTRACT We collected the optical light-curve data of 227 gamma-ray bursts (GRBs) observed with the TAROT, COATLI, and RATIR telescopes. These consist of 133 detections and 94 upper limits. We constructed average light curves in the observer and rest frames in both X-rays (from Swift/X-Ray Telescope) and the optical. Our analysis focused on investigating the observational and intrinsic properties of GRBs. Specifically, we examined observational properties, such as the optical brightness function of the GRBs at T = 1000 s after the trigger, as well as the temporal slope of the afterglow. We also estimated the redshift distribution for the GRBs within our sample. Of the 227 GRBs analysed, we found that 116 had a measured redshift. Based on these data, we calculated a local rate of ρ0 = 0.2 Gpc−3 yr−1 for these events with z < 1. To explore the intrinsic properties of GRBs, we examined the average X-ray and optical light curves in the rest frame. We use the afterglowpy library to generate synthetic curves to constrain the parameters typical of the bright GRB jet, such as energy (〈E0〉 ∼ 1053.6 erg), opening angle (〈θcore〉 ∼ 0.2 rad), and density (〈n0〉 ∼ 10−2.1 cm−3). Furthermore, we analyse microphysical parameters, including the fraction of thermal energy in accelerated electrons (〈ϵe〉 ∼ 10−1.37) and in the magnetic field (〈ϵB〉 ∼ 10−2.26), and the power-law index of the population of non-thermal electrons (〈p〉 ∼ 2.2).

Deciphering Faint Gyrosynchrotron Emission from a Coronal Mass Ejection Using Spectropolarimetric Radio Imaging

Measurements of the plasma parameters of coronal mass ejections (CMEs), particularly the magnetic field and nonthermal electron population entrained in the CME plasma, are crucial to understand their propagation, evolution, and geo-effectiveness. Spectral modeling of gyrosynchrotron (GS) emission from CME plasma has been regarded as one of the most promising remote-sensing techniques for estimating spatially resolved CME plasma parameters. Imaging the very low flux density CME GS emission in close proximity to the Sun with orders of magnitude higher flux density has, however, proven to be rather challenging. This challenge has only recently been met using the high dynamic range imaging capability of the Murchison Widefield Array (MWA). Although routine detection of GS is now within reach, the challenge has shifted to constraining the large number of free parameters in GS models, a few of which are degenerate, using the limited number of spectral points at which the observations are typically available. These degeneracies can be broken using polarimetric imaging. For the first time, we demonstrate this using our recently developed capability of high-fidelity polarimetric imaging on the data from the MWA. We show that spectropolarimetric imaging, even when only sensitive upper limits on circularly polarization flux density are available, is not only able to break the degeneracies but also yields tighter constraints on the plasma parameters of key interest than possible with total intensity spectroscopic imaging alone.

Modelling kilonova afterglows: Effects of the thermal electron population and interaction with GRB outflows

AbstractGiven an increasing number of gamma-ray bursts accompanied by potential kilonovae, there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distribution affect kilonova afterglow light curves. We employ semi-analytic afterglow model, PyBlastAfterglow. We consider kilonova ejecta profiles from ab-initio numerical relativity binary neutron star merger simulations, targeted to GW170817. We do not perform model selection. We find that the emission from thermal electrons dominates at early times. If the interstellar medium density is high (${\simeq }0.1\, \, \text{cm}^{-3}$), it adds an early time peak to the light curve. As ejecta decelerates, the spectral and temporal indexes change in a characteristic way that, if observed, can be used to reconstruct the ejecta velocity distribution. For the low interstellar medium density, inferred for GRB 170817A, the emission from the non-thermal electron population generally dominates. We also assess how kilonova afterglow light curves change if the interstellar medium has been partially removed and pre-accelerated by laterally expanding gamma-ray burst ejecta. For the latter, we consider properties informed by observations of GRB170817A. We find that the main effect is the emission suppression at early time ${\lesssim }10^{3}\,$ days, and at its maximum it reaches ${\sim }40{{\ \rm per\ cent}}$ when the fast tail of the kilonova ejecta moves subsonically through the wake of laterally spreading gamma-ray burst ejecta. The subsequent rebrightening, when these ejecta break through and shocks form, is very mild (${\lesssim }10{{\ \rm per\ cent}}$) and may not be observable.

Probing Plasma Composition with the Next Generation Event Horizon Telescope (ngEHT)

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models, including a constant electron to magnetic pressure (constant βe model) with a population of non-thermal electrons as well as an R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations are enhanced with every pair added. On the contrary, we show that, at higher frequencies, a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of a large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function.

Highly Energetic Electrons Accelerated in Strong Solar Flares as a Preferred Driver of Sunquakes

Sunquakes are enhanced seismic waves excited in some energetic solar flares. Up to now, their origin has still been controversial. In this Letter, we select and study 20 strong flares in Solar Cycle 24, whose impulse phase is fully captured by the Reuven Ramaty High Energy Solar Spectroscopic Imager. For 11 out of 12 sunquake-active flares in our sample, the hard X-ray emission shows a good temporal and spatial correlation with the white-light enhancement and the sunquake. Spectral analysis also reveals a harder photon spectrum that extends to several hundred keV, implying a considerable population of flare-accelerated nonthermal electrons at high energies. Quantitatively, the total energy of electrons above 300 keV in sunquake-active flares is systematically different from that in sunquake-quiet flares, while the difference is marginal for electrons above 50 keV. All these facts support highly energetic electrons as a preferred driver of the sunquakes. Such an electron-driven scenario can be reasonably accommodated in the framework of a recently proposed selection rule for sunquake generation. For the remaining one event, the sunquake epicenter is cospatial with a magnetic imprint, i.e., a permanent change of magnetic field on the photosphere. Quantitative calculation shows that the flare-induced downward Lorentz force can do enough work to power the sunquake, acting as a viable sunquake driver for this specific event.

Tailoring surface topographies on solids with Mid-IR femtosecond laser pulses

Irradiation of solids with ultrashort pulses using laser sources in the mid-infrared (mid-IR) spectral region is a yet predominantly unexplored field that opens broad possibilities for efficient and precise surface texturing for a wide range of applications. In the present work, we investigate both experimentally and theoretically the impact of laser sources on the generation of surface modification related effects and on the subsequent surface patterning of metallic and semiconducting materials. Through a parametric study we correlate the mid-IR pulsed laser parameters with the onset of material damage and the formation of a variety of periodic surface structures at a laser wavelength of λL = 3200 nm and a pulse duration of τp = 45 fs. Results for nickel and silicon indicate that the produced topographies comprise both high and low spatial frequency induced periodic structures, similar to those observed at lower wavelengths, while groove formation is absent. The investigation of the damage thresholds suggests the incorporation of nonlinear effects generated from three-photon-assisted excitation (for silicon) and the consideration of the role of the non-thermal excited electron population (for nickel) at very short pulse durations. The results demonstrate the potential of surface structuring with mid-IR pulses, which can constitute a systematic novel engineering approach with strong fields at long-wavelength spectral regions that can be used for advanced industrial laser applications.

Ultrasensitive probing of plasmonic hot electron occupancies

Non-thermal and thermal carrier populations in plasmonic systems raised significant interest in contemporary fundamental and applied physics. Although the theoretical description predicts not only the energies but also the location of the generated carriers, the experimental justification of these theories is still lacking. Here, we demonstrate experimentally that upon the optical excitation of surface plasmon polaritons, a non-thermal electron population appears in the topmost domain of the plasmonic film directly coupled to the local fields. The applied all-optical method is based on spectroscopic ellipsometric determination of the dielectric function, allowing us to obtain in-depth information on surface plasmon induced changes of the directly related electron occupancies. The ultrahigh sensitivity of our method allows us to capture the signatures of changes induced by electron-electron scattering processes with ultrafast decay times. These experiments shed light on the build-up of plasmonic hot electron population in nanoscale media.

Polarization-sensitive Compton Scattering by Accelerated Electrons

We describe upgrades to a numerical code that computes synchrotron and inverse-Compton emission from relativistic plasma including full polarization. The introduced upgrades concern a scattering kernel that is now capable of scattering the polarized and unpolarized photons on a nonthermal population of electrons. We describe the scheme to approach this problem and we test the numerical code against a known analytic solution. Finally, using the upgraded code, we predict the polarization of light that is scattered off subrelativistic thermal or relativistic thermal and nonthermal free electrons. The upgraded code enables more realistic simulations of emissions from plasma jets associated with accreting compact objects.

Observation of non-thermal electrons outside the SOL in the Wendelstein 7-X stellarator

• Non-thermal electrons can be observed far outside the main plasma in W7-X. • Occurrence is limited to electron-cyclotron resonance heated plasmas at low density. • Electrons cause hot spots on plasma-facing components with considerable heat fluxes. Hot spots and bright patterns on plasma-facing components (PFC) in visible light are observed in electron cyclotron resonance heated (ECRH) plasmas at low densities in the stellarator Wendelstein 7-X. The events are often located far outside of any convective plasma loads and led in some cases to damages of diagnostic components, where the interaction zones qualitatively agree with fast particle loss simulations. Reciprocating electric probes indicate that this phenomenon is related to a non-thermal electron population that can have a beam like character, being directed in one parallel direction. Further, the electrons can be trapped on rational flux surfaces and be used to map magnetic islands.

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.

Call and Response: A Time-resolved Study of Chromospheric Evaporation in a Large Solar Flare

We studied an X1.6 solar flare produced by NOAA Active Region 12602 on 2014 October 22. The entirety of this event was covered by RHESSI, IRIS, and Hinode/EIS, allowing analysis of the chromospheric response to a nonthermal electron driver. We derived the energy contained in nonthermal electrons via RHESSI spectral fitting and linked the time-dependent parameters of this call to the response in Doppler velocity, density, and nonthermal width across a broad temperature range. The total energy injected was 4.8 × 1030 erg and lasted 352 s. This energy drove explosive chromospheric evaporation, with a delineation in both Doppler and nonthermal velocities at the flow reversal temperature, between 1.35 and 1.82 MK. The time of peak electron injection (14:06 UT) corresponded to the time of highest velocities. At this time, we found 200 km s−1 blueshifts in the core of Fe xxiv, which is typically assumed to be at rest. Shortly before this time, the nonthermal electron population had the shallowest spectral index (≈6), corresponding to the peak nonthermal velocity in Si iv and Fe xxi. Nonthermal velocities in Fe xiv, formed near the flow reversal temperature, were low and not correlated with density or Doppler velocity. Nonthermal velocities in ions with similar temperatures were observed to increase and correlate with Doppler velocities, implying unresolved flows surrounding the flow reversal point. This study provides a comprehensive, time-resolved set of chromospheric diagnostics for a large X-class flare, along with a time-resolved energy injection profile, ideal for further modeling studies.

Plasma dynamics in the flaring loop observed by RHESSI

Context. Hard X-rays (HXRs) contain the most direct information about the non-thermal electron population in solar flares. The approximation of the HXR emission mechanism (bremsstrahlung), known as the thick-target model, is well developed. It allows one to diagnose the physical conditions within a flaring structure. The thick-target model predicts that in flare foot points, we should observe lowering of HXR sources’ altitude with increasing energy. Aims. The foot point of HXR sources result from the direct interaction of non-thermal electron beams with plasma in the lower part of the solar atmosphere, where the density increases rapidly. Therefore, we can estimate the plasma density distribution along the non-thermal electron beam directly from the observations of the altitude-energy relation obtained for the HXR foot point sources. However, the relation is not only density-dependent. Its shape is also determined by the power-law distribution of non-thermal electrons. Additionally, during the impulsive phase, the plasma density and a degree of ionisation within foot points may change dramatically due to heating and chromospheric evaporation. For this reason, the interpretation of observed HXR foot point sources’ altitudes is not straightforward and needs detailed numerical modelling of the electron precipitation process. Methods. We present the results of numerical modelling of one well-observed solar flare. We used HXR observations obtained by RHESSI. The numerical model was calculated using the hydrodynamic 1D model with an application of the Fokker-Planck formalism for non-thermal beam precipitation. Results. HXR data were used to trace chromospheric density changes during a non-thermal emission burst, in detail. We have found that the amount of mass that evaporated from the chromosphere is in the range of 2.7 × 1013 − 4.0 × 1014 g. This is in good agreement with the ranges obtained from hydrodynamical modelling of a flaring loop (2.3 × 1013 − 3.3 × 1013 g), and from an analysis of observed emission measure in the loop top (3.9 × 1013 − 5.3 × 1013 g). Additionally, we used specific scaling laws which gave another estimation of the evaporated mass around 2 × 1014 g. Conclusions. Consistency between the obtained values shows that HXR images may provide an important constraint for models – a mass of plasma that evaporated due to a non-thermal electron beam depositing energy in the chromosphere. High-energy, non-thermal sources’ (above 20 keV in this case) positions fit the column density changes obtained from the hydrodynamical model perfectly. Density changes seem to be less affected by the electrons’ spectral index. The obtained results significantly improve our understanding of non-thermal electron beam precipitation and allow us to refine the energy balance in solar flare foot points during the impulsive phase.

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

Electrons accelerated on Earth by a rich variety of wave scattering or stochastic processes generate hard non-thermal X-ray bremsstrahlung up to >~ 1 MeV and power Earth's various types of aurorae. Although Jupiter's magnetic field is an order of magnitude larger than Earth's, space-based telescopes have previously detected X-rays only up to ~7 keV. On the basis of theoretical models of the Jovian auroral X-ray production, X-ray emission in the ~2-7 keV band has been interpreted as thermal (arising from electrons characterized by a Maxwell-Boltzmann distribution) bremsstrahlung. Here we report the observation of hard X-rays in the 8-20 keV band from the Jovian aurorae, obtained with the NuSTAR X-ray observatory. The X-rays fit to a flat power-law model with slope 0.60+/-0.22 - a spectral signature of non-thermal, hard X-ray bremsstrahlung. We determine the electron flux and spectral shape in the keV to MeV energy range using coeval in situ measurements by the Juno spacecraft's JADE and JEDI instruments. Jovian electron spectra of the form we observe have previously been interpreted to arise in stochastic acceleration, rather than coherent acceleration by electric fields. We reproduce the X-ray spectral shape and approximate flux observed by NuSTAR, and explain the non-detection of hard X-rays by Ulysses, by simulating the non-thermal population of electrons undergoing precipitating electron energy loss, secondary electron generation and bremsstrahlung emission in a model Jovian atmosphere. The results highlight the similarities between the processes generating hard X-ray auroras on Earth and Jupiter, which may be occurring on Saturn, too.

Observation and origin of non-thermal hard X-rays from Jupiter

Electrons accelerated on Earth by a rich variety of wave-scattering or stochastic processes1,2 generate hard, non-thermal X-ray bremsstrahlung up to ~1 MeV (refs. 3,4) and power Earth’s various types of aurorae. Although Jupiter’s magnetic field is an order of magnitude larger than Earth’s, space-based telescopes have previously detected X-rays only up to ~7 keV (ref. 5). On the basis of theoretical models of the Jovian auroral X-ray production6,7,8, X-ray emission in the ~2–7 keV band has been interpreted as thermal (arising from electrons characterized by a Maxwell–Boltzmann distribution) bremsstrahlung5,9. Here we report the observation of hard X-rays in the 8–20 keV band from the Jovian aurorae, obtained with the NuSTAR X-ray observatory. The X-rays fit to a flat power-law model with slope of 0.60 ± 0.22—a spectral signature of non-thermal, hard X-ray bremsstrahlung. We determine the electron flux and spectral shape in the kiloelectronvolt to megaelectronvolt energy range using coeval in situ measurements taken by the Juno spacecraft’s JADE and JEDI instruments. Jovian electron spectra of the form we observe have previously been interpreted as arising in stochastic acceleration, rather than coherent acceleration by electric fields10. We reproduce the X-ray spectral shape and approximate flux observed by NuSTAR, and explain the non-detection of hard X-rays by Ulysses11, by simulating the non-thermal population of electrons undergoing precipitating electron energy loss, secondary electron generation and bremsstrahlung emission in a model Jovian atmosphere. The results highlight the similarities between the processes generating hard X-ray aurorae on Earth and Jupiter, which may be occurring on Saturn, too.

The Intrinsic Structure of Sagittarius A* at 1.3 cm and 7 mm

Abstract Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets in which to resolve the innermost region of an SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multiwavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ∼704 ± 102 μas (axial ratio ∼ 1.19 − 0.19 + 0.24 ) and ∼300 ± 25 μas (axial ratio ∼1.28 ± 0.2) at 1.349 cm and 6.95 mm, respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power law, with an index ∼1.2 ± 0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes.

Thermal Electrons in Mildly Relativistic Synchrotron Blast Waves

Abstract Numerical models of collisionless shocks robustly predict an electron distribution composed of both thermal and nonthermal electrons. Here, we explore in detail the effect of thermal electrons on the emergent synchrotron emission from subrelativistic shocks. We present a complete “thermal + nonthermal” synchrotron model and derive properties of the resulting spectrum and light curves. Using these results, we delineate the relative importance of thermal and nonthermal electrons for subrelativistic shock-powered synchrotron transients. We find that thermal electrons are naturally expected to contribute significantly to the peak emission if the shock velocity is ≳0.2c, but would be mostly undetectable in nonrelativistic shocks. This helps explain the dichotomy between typical radio supernovae and the emerging class of “AT2018cow-like” events. The signpost of thermal electron synchrotron emission is a steep optically-thin spectral index and a ν 2 optically-thick spectrum. These spectral features are also predicted to correlate with a steep postpeak light-curve decline rate, broadly consistent with observed AT2018cow-like events. We expect that thermal electrons may be observable in other contexts where mildly relativistic shocks are present and briefly estimate this effect for gamma-ray burst afterglows and binary–neutron-star mergers. Our model can be used to fit spectra and light curves of events and accounts for both thermal and nonthermal electron populations with no additional physical degrees of freedom.