skip to main content

Title: Flux variability from ejecta in structured relativistic jets with large-scale magnetic fields
Context. Standing and moving shocks in relativistic astrophysical jets are very promising sites for particle acceleration to large Lorentz factors and for the emission from the radio up to the γ -ray band. They are thought to be responsible for at least part of the observed variability in radio-loud active galactic nuclei. Aims. We aim to simulate the interactions of moving shock waves with standing recollimation shocks in structured and magnetized relativistic jets and to characterize the profiles of connected flares in the radio light curve. Methods. Using the relativistic magneto-hydrodynamic code MPI-AMRVAC and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dynamics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks and on the stratified jet morphology. Secondly, we study the interaction of a moving shock wave with the standing shocks. We calculated the synthetic synchrotron maps and radio light curves and analyze the variability at two frequencies 1 and 15.3 GHz and for several observation angles. Finally, we compare the characteristics of our simulated light curves more » with radio flares observed from the blazar 3C 273 with the Owens Valley Radio Observatory and Very Long Baseline Array in the MOJAVE survey between 2008 and 2019. Results. We find that in a structured over-pressured relativistic jet, the presence of the large-scale magnetic field structure changes the properties of the standing shock waves and leads to an opening in the jet. The interaction between waves from inner and outer jet components can produce strong standing shocks. When crossing such standing shocks, moving shock waves accompanying overdensities injected in the base of the jet cause very luminous radio flares. The observation of the temporal structure of these flares under different viewing angles probes the jet at different optical depths. At 1 GHz and for small angles, the self-absorption caused by the moving shock wave becomes more important and leads to a drop in the observed flux after it interacts with the brightest standing knot. A weak asymmetry is seen in the shape of the simulated flares, resulting from the remnant emission of the shocked standing shocks. The characteristics of the simulated flares and the correlation of peaks in the light curve with the crossing of moving and standing shocks favor this scenario as an explanation of the observed radio flares of 3C 273. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Astronomy & Astrophysics
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate non-thermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anticorrelated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the periodmore »of the QPOs is associated with the kink growth time-scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities.« less
  2. ABSTRACT V445 Puppis is the only helium nova observed to date; its eruption in late 2000 showed high velocities up to 8500 km s−1, and a remarkable bipolar morphology cinched by an equatorial dust disc. Here we present multifrequency radio observations of V445 Pup obtained with the Very Large Array (VLA) spanning 1.5–43.3 GHz, and between 2001 January and 2008 March (days ∼89–2700 after eruption). The radio light curve is dominated by synchrotron emission over these 7 yr, and shows four distinct radio flares. Resolved radio images obtained in the VLA’s A configuration show that the synchrotron emission hugs the equatorial disc, and comparisons to near-IR images of the nova clearly demonstrate that it is the densest ejecta – not the fastest ejecta – that are the sites of the synchrotron emission in V445 Pup. The data are consistent with a model where the synchrotron emission is produced by a wind from the white dwarf impacting the dense equatorial disc, resulting in shocks and particle acceleration. The individual synchrotron flares may be associated with density enhancements in the equatorial disc and/or velocity variations in the wind from the white dwarf. This overall scenario is similar to a common picture of shock production in hydrogen-rich classical novae,more »but V445 Pup is remarkable in that these shocks persist for almost a decade, much longer than the weeks or months for which shocks are typically observed in classical novae.« less

    We present one-dimensional hydrodynamical simulations including radiative losses, of internal shocks in the outflows from classical novae, to explore the role of shocks in powering multiwavelength emission from radio to gamma-ray wavelengths. Observations support a picture in which the initial phases of some novae generate a slow, equatorially focused outflow (directly from the outer Lagrange point, or from a circumbinary disc), which then transitions to, or is overtaken by, a faster more isotropic outflow from the white dwarf which collides and shocks the slower flow, powering gamma-ray and optical emission through reprocessing by the ejecta. However, the common occurrence of multiple peaks in nova light curves suggests that the outflow’s acceleration need not be monotonic, but instead can involve successive transitions between ‘fast’ and ‘slow’ modes. Such a time-fluctuating outflow velocity naturally can reproduce several observed properties of nova, such as correlated gamma-ray and optical flares, expansion of the photosphere coincident with (though lagging slightly) the peak flare luminosity, and complex time evolution of spectral lines (including accelerating, decelerating, and merging velocity components). While the shocks are still deeply embedded during the gamma-ray emission, the onset of ∼keV X-ray and ∼10 GHz radio synchrotron emission is typically delayed until themore »forward shock of the outermost monolithic shell (created by merger of multiple internal shock-generated shells) reaches a sufficiently low column through the dense external medium generated by the earliest phase of the outburst.

    « less
  4. Abstract We present a population of 19 radio-luminous supernovae (SNe) with emission reaching L ν ∼ 10 26 –10 29 erg s −1 Hz −1 in the first epoch of the Very Large Array Sky Survey (VLASS) at 2–4 GHz. Our sample includes one long gamma-ray burst, SN 2017iuk/GRB 171205A, and 18 core-collapse SNe detected at ≈1–60 yr after explosion. No thermonuclear explosion shows evidence for bright radio emission, and hydrogen-poor progenitors dominate the subsample of core-collapse events with spectroscopic classification at the time of explosion (79%). We interpret these findings in the context of the expected radio emission from the forward shock interaction with the circumstellar medium (CSM). We conclude that these observations require a departure from the single wind–like density profile (i.e., ρ CSM ∝ r −2 ) that is expected around massive stars and/or from a spherical Newtonian shock. Viable alternatives include the shock interaction with a detached, dense shell of CSM formed by a large effective progenitor mass-loss rate, M ̇ ∼ 10 − 4 – 10 − 1 M ⊙ yr −1 (for an assumed wind velocity of 1000 km s −1 ); emission from an off-axis relativistic jet entering our line of sight; ormore »the emergence of emission from a newly born pulsar-wind nebula. The relativistic SN 2012ap that is detected 5.7 and 8.5 yr after explosion with L ν ∼ 10 28 erg s −1 Hz −1 might constitute the first detections of an off-axis jet+cocoon system in a massive star. However, none of the VLASS SNe with archival data points are consistent with our model off-axis jet light curves. Future multiwavelength observations will distinguish among these scenarios. Our VLASS source catalogs, which were used to perform the VLASS cross-matching, are publicly available at .« less
  5. ABSTRACT During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS’s open magnetic flux – and hence increases its wind luminosity – relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of ∼1 Jy·ms at ∼GHz frequencies lasts ∼1–500 ms following the merger for a source at ∼3 Gpc (Bd/1012 G)8/9, where Bd is the dipole field strength of the more strongly magnetized star.more »Given an outflow geometry concentrated along the binary equatorial plane, the signal may be preferentially observable for high-inclination systems, that is, those least likely to produce a detectable gamma-ray burst.« less