Search

We present JWST observations of the Crab Nebula, the iconic remnant of the historical SN 1054. The observations include NIRCam and MIRI imaging mosaics, plus MIRI/MRS IFU spectra that probe two select locations within the ejecta filaments. We derive a high-resolution map of dust emission and show that the grains are concentrated in the innermost, high-density filaments. These dense filaments coincide with multiple synchrotron bays around the periphery of the Crab's pulsar wind nebula (PWN). We measure synchrotron spectral index changes in small-scale features within the PWN's torus region, including the well-known knot and wisp structures. The index variations are consistent with Doppler boosting of emission from particles with a broken power-law distribution, providing the first direct evidence that the curvature in the particle injection spectrum is tied to the acceleration mechanism at the termination shock. We detect multiple nickel and iron lines in the ejecta filaments and use photoionization models to derive nickel-to-iron abundance ratios that are a factor of 3-8 higher than the solar ratio. We also find that the previously reported order-of-magnitude higher Ni/Fe values from optical data are consistent with the lower values from JWST when we reanalyze the optical emission using updated atomic data and account for local extinction from dust. We discuss the implications of our results for understanding the nature of the explosion that produced the Crab Nebula and conclude that the observational properties are most consistent with a low-mass iron-core-collapse supernova, even though an electron-capture explosion cannot be ruled out., Comment: 32 pages, 3 tables, 20 figures, accepted for publication in ApJL

We provide an overview of various mechanisms, and corresponding powers, of precursor emission to compacts' mergers to be detected by LIGO-Virgo-KAGRA (LVK) collaboration. Expected peak powers, $\leq 10^{43}$ erg s$^{-1}$, are not sufficiently high to be detected by all-sky high-energy satellites (unless beamed). The best chance is the detection of possible coherent radio emission, producing observable signals up to $\sim$ Jansky of flux density. Low-frequency phased array telescopes like LOFAR, the MWA and DSA-2000 are best suited due to their large instantaneous sky coverage. Time-wise, in addition to LIGO early warning alerts up to a minute before the merger, the dispersive delay at lower frequencies of $\sim$ 300 MHz can be of the order of minutes. Optical detections are the most challenging.

J191213.72-441045.1 is a binary system composed of a white dwarf and an M-dwarf in a 4.03-hour orbit. It shows emission in radio, optical, and X-ray, all modulated at the white dwarf spin period of 5.3 min, as well as various orbital sideband frequencies. Like in the prototype of the class of radio-pulsing white dwarfs, AR Scorpii, the observed pulsed emission seems to be driven by the binary interaction. In this work, we present an analysis of far-ultraviolet spectra obtained with the Cosmic Origins Spectrograph at the Hubble Space Telescope, in which we directly detect the white dwarf in J191213.72-441045.1. We find that the white dwarf has an effective temperature of 11485+/-90 K and mass of 0.59+/-0.05 solar masses. We place a tentative upper limit on the magnetic field of ~50 MG. If the white dwarf is in thermal equilibrium, its physical parameters would imply that crystallisation has not started in the core of the white dwarf. Alternatively, the effective temperature could have been affected by compressional heating, indicating a past phase of accretion. The relatively low upper limit to the magnetic field and potential lack of crystallisation that could generate a strong field pose challenges to pulsar-like models for the system and give preference to propeller models with a low magnetic field. We also develop a geometric model of the binary interaction which explains many salient features of the system., Comment: 12 pages, 17 figures. Accepted for publication in MNRAS

Radio emission from magnetars provides a unique probe of the relativistic, magnetized plasma within the near-field environment of these ultra-magnetic neutron stars. The transmitted waves can undergo birefringent and dispersive propagation effects that result in frequency-dependent conversions of linear to circularly polarized radiation and vice-versa, thus necessitating classification when relating the measured polarization to the intrinsic properties of neutron star and fast radio burst (FRB) emission sites. We report the detection of such behavior in 0.7-4 GHz observations of the P = 5.54 s radio magnetar XTE J1810$-$197 following its 2018 outburst. The phenomenon is restricted to a narrow range of pulse phase centered around the magnetic meridian. Its temporal evolution is closely coupled to large-scale variations in magnetic topology that originate from either plastic motion of an active region on the magnetar surface or free precession of the neutron star crust. Our model of the effect deviates from simple theoretical expectations for radio waves propagating through a magnetized plasma. Birefringent self-coupling between the transmitted wave modes, line-of-sight variations in the magnetic field direction and differences in particle charge or energy distributions above the magnetic pole are explored as possible explanations. We discuss potential links between the immediate magneto-ionic environments of magnetars and those of FRB progenitors., Comment: 22 pages (inc. extended data and supplementary materials), 12 figures, 3 tables. Published in Nature Astronomy

Using the ONEDFEL code we perform Free Electron Laser simulations in the astrophysically important guide-field dominated regime. For wigglers' (Alfven waves) wavelengths of tens of meters and beam Lorentz factor $\sim 10^3$, the resulting coherently emitted waves are in the centimeter range. Our simulations show a growth of the wave intensity over fourteen orders of magnitude, over the astrophysically relevant scale of $\sim$ few kilometers. The signal grows from noise (unseeded). The resulting spectrum shows fine spectral sub-structures, reminiscent of the ones observed in Fast Radio Bursts (FRBs)., Comment: accepted by JPP

We reconsider the escape of high brightness coherent emission of Fast Radio Bursts (FRBs) from magnetars' magnetospheres, and conclude that there are numerous ways for the powerful FRB pulse to avoid nonlinear absorption. Sufficiently strong surface fields, $\geq 10\%$ of the quantum field, limit the waves' non-linearity to moderate values. For weaker fields, the electric field experienced by a particle is limited by a combined ponderomotive and parallel-adiabatic forward acceleration of charges by the incoming FRB pulse along the magnetic field lines newly opened during FRB/Coronal Mass Ejection (CME). As a result, particles surf the weaker front part of the pulse, experiencing low radiative losses, and are cleared from the magnetosphere for the bulk of the pulse to propagate. We also find: (i) for propagation across magnetic field, the O-mode suffers much smaller dissipation than the X-mode; (ii) quasi-parallel propagation suffers minimal dissipation; (iii) initial mildly relativistic radial plasma flow further reduces losses; (iv) for oblique propagation of a pulse with limited transverse size, the leading part of the pulse would ponderomotively sweep the plasma aside.

We point out the dominant importance of plasma injection effects for relativistic winds from pulsars and black holes. We demonstrate that outside the light cylinder the magnetically dominated outflows while sliding along the helical magnetic field move in fact nearly radially with very large Lorentz factors $\gamma_0 \gg 1 $, imprinted into the flow during pair production within the gaps. Only at larger distances, $r \geq \gamma_0 (c/\Omega)$, the MHD acceleration $\Gamma \propto r$ takes over. As a result, Blandford-Znajek (BZ) driven outflows would produce spine-brightened images. The best-resolved case of the jet in M87 shows both bright edge-brightened features, as well as weaker spine-brightened feature. Only the spine-brightened component can be BZ-driven/originate from the BH's magnetosphere.

We study dynamics of relativistic Coronal Mass Ejections (CMEs), from launching by shearing of foot-points (either slowly - the ``Solar flare'' paradigm, or suddenly - the ``star quake" paradigm), to propagation in the preceding magnetar wind. For slow shear, most of the energy injected into the CME is first spent on the work done on breaking through the over-laying magnetic field. At later stages, sufficiently powerful CMEs may experience ``detonation" and lead to opening of the magnetosphere beyond some equipartition radius $r_{eq}$, where the energy of the CME becomes larger than the decreasing external magnetospheric energy. Post-CME magnetosphere relaxes via formation of a plasmoid-mediated current sheet, initially at $\sim r_{eq}$ and slowly reaching the light cylinder (this transient stage has much higher spindown rate and may produce an ``anti-glitch''). Both the location of the foot-point shear and the global magnetospheric configuration affect the frequent-and-weak versus rare-and-powerful CME dichotomy - to produce powerful flares the slow shear should be limited to field lines that close near the star. After the creation of a topologically disconnected flux tube, the tube quickly (at $\sim$ the light cylinder) comes into force-balance with the preceding wind, and is passively advected/frozen in the wind afterward. For fast shear (a local rotational glitch), the resulting large amplitude Alfven waves lead to opening of the magnetosphere (which later recovers similarly to the slow shear case). At distances much larger than the light cylinder, the resulting shear Alfven waves propagate through the wind non-dissipatively. Implications to Fast Radio Bursts are discussed.

We first derive a set of equations describing general stationary configurations of relativistic force-free plasma, without assuming any geometric symmetries. We then demonstrate that electromagnetic interaction of merging neutron stars is necessarily dissipative due to the effect of electromagnetic draping - creation of dissipative regions near the star (in the single-magnetized case) or at the magnetospheric boundary (in the double-magnetized case). Our results indicate that even in the single magnetized case we expect that relativistic jets (or ``tongues'') are produced, with correspondingly beamed emission pattern.

In gamma-ray binaries neutron star is orbiting a companion that produces a strong stellar wind. We demonstrate that observed properties of "stellar wind"-"pulsar wind" interaction depend both on the overall wind thrust ratio, as well as more subtle geometrical factors: the relative direction of the pulsar's spin, the plane of the orbit, the direction of motion, and the instantaneous line of sight. Using fully 3D relativistic magnetohydrodynamical simulations we find that the resulting intrinsic morphologies can be significantly orbital phase-dependent: a given system may change from tailward-open to tailward-closed shapes. As a result, the region of unshocked pulsar wind can change by an order of magnitude over a quarter of the orbit. We calculate radiation maps and synthetic light curves for synchrotron (X-ray) and Inverse-Compton emission (GeV-TeV), taking into account $\gamma-\gamma$ absorption. Our modeled light curves are in agreement with the phase-dependent observed light curves of LS5039., Comment: 22 pages, 20 figures

We reconsider the dynamics of accretion flows onto magnetized central star. For dipolar magnetically aligned case, the centrifugal barrier is at $R_{cb} = (2/3)^{1/3} R_c = 0.87 R_c$, where $R_c= ( G M/\Omega^2)^{1/3}$ is the corotation radius. For oblique dipole direct accretion from the corotation radius $R_c$ is possible only for magnetic obliquity satisfying $\tan \theta_\mu \geq 1/( 2 \sqrt{3}) $ ($\theta_\mu \geq 16.1^\circ $). The accretion proceeds in a form of funnel flows - along two streams centered on the $\mu-\Omega$ plane, with azimuthal opening angle $ \cos (\Delta \phi) = { \cot^ 2 {\theta_\mu} }/{12} $. For the magnetosphere distorted by the diamagnetic disk, the centrifugal barrier can be at as small radius as $R_{cb}= 0.719 R_c$ for the fully confined dipole, extending out to $R_{cb} \sim R_c$ for the magnetically balanced case. Type-II X-ray bursts in accreting neutron stars may be mediated by the centrifugal barrier; this requires nearly aligned configuration. Centrifugally-barriered material trapped in the magnetosphere may lead to periodic obscuration ("dips") in the light curve of the host star, e.g., as observed in accreting young stellar objects and X-ray binaries.

Rich observational phenomenology associated with Pulsar B in PSR J0737$-$3039A/B system resembles in many respects phenomena observed in the Earth and Jupiter magnetospheres, originating due to the wind-magnetosphere interaction. We consider particle dynamics in the fast corotating magnetosphere of Pulsar B, when the spin period is shorter than the third adiabatic period. We demonstrate that trapped particles occasionally experience large radial variations of the L-parameter (effective radial distance) due to the parametric interaction of the gyration motion with the large scale electric fields induced by the deformations of the magnetosphere, in what could be called a betatron-induced diffusion. The dynamics of particles from the wind of Pulsar A trapped inside Pulsar B magnetosphere is governed by Mathieu's equation, so that the parametrically unstable orbits are occasionally activated; particle dynamics is not diffusive per se. The model explains the high plasma density on the closed field lines of Pulsar B, and the fact that the observed eclipsing region is several times smaller than predicted by the hydrodynamic models.

We consider propagation of polarization in the inner parts of pair-symmetric magnetar winds, close to the light cylinder. Pair plasmas in magnetic field is birefringent, a $\propto B^2$ effect. As a result, such plasmas work as phase retarders: Stokes parameters follow a circular trajectory on the Poincare sphere. In the highly magnetized regime, $\omega, \, \omega_p \ll \omega_B$, the corresponding rotation rates are independent of the magnetic field. A plasma screen with dispersion measure DM $\sim 10^{-6}$ pc cm$^{-3}$ can induce large polarization changes, including large effective Rotation Measure (RM). The frequency scaling of the (generalized) RM, $ \propto \lambda ^\alpha $, mimics the conventional RM with $\alpha =2$ for small phase shifts, but can be as small as $\alpha =1$. In interpreting observations the frequency scaling of polarization parameters should be fitted independently. The model offers explanations for (i) large circular polarization component observed in FRBs, with right-left switching; (ii) large RM, with possible sign changes; (iii) time-depend variable polarization. Relatively dense and slow wind is needed - the corresponding effect in regular pulsars is small., Comment: arXiv admin note: text overlap with arXiv:2111.00281

Many explosive astrophysical events, like magnetars' bursts and flares, are magnetically driven. We consider dynamics of such magnetic explosions - relativistic expansion of highly magnetized and highly magnetically over-pressurized clouds. The corresponding dynamics is qualitatively different from fluid explosions due to the topological constraint of the conservation of the magnetic flux. Using analytical, relativistic MHD as well as force-free calculations, we find that the creation of a relativistically expanding, causally disconnected flow obeys a threshold condition: it requires sufficiently high initial over-pressure and sufficiently quick decrease of the pressure in the external medium (the pre-explosion wind). In the subcritical case the magnetic cloud just "puffs-up" and quietly expands with the pre-flare wind. We also find a compact analytical solution to the Prendergast's problem - expansion of force-free plasma into vacuum., Comment: 31 pages, 23 figures, submitted to ApJ

Short rise times of Fast Blue Optical Transients (FBOTs) require very light ejected envelopes, $M_{ej} \leq 10^{-1} M_\odot$, much smaller than of a typical supernova. Short peak times also mean that FBOTs should be hydrodynamically, not radioactively powered. The detection by Chandra of X-ray emission in AT2020mrf of $L_X \sim 10^{42} $ erg s$^{-1}$ after 328 days implies total, overall dominant, X-ray energetics at the Gamma Ray Bursts (GRBs) level of $\sim 6 \times 10^{49}$ erg. We further develop a model of Lyutikov & Toonen (2019), whereby FBOTs are the results of a late accretion induced collapse (AIC) of the product of super-Chandrasekhar double white dwarf (WD) merger between ONeMg WD and another WD. Small ejecta mass, and the rarity of FBOTs, result from the competition between mass loss from the merger product to the wind, and ashes added to the core, on time scale of $\sim 10^3-10^4$ years. FBOTs occur only when the envelope mass before AIC is $\leq 10^{-1} M_\odot$. FBOTs proper come from central engine-powered radiation-dominated forward shock as it propagates through ejecta. FBOTs' duration is determined by the diffusion time of photons produced by the NS-driven forward shock within the expanding ejecta. All the photons produced by the central source deep inside the ejecta escape almost simultaneously, producing a short bright event. The high energy emission is generated at the highly relativistic and highly magnetized termination shock, qualitatively similar to Pulsar Wind Nebulae. The X-ray bump observed in AT2020mrf by SRG/eROSITA, predicted by Lyutikov & Toonen (2019), is coming from the break-out of the engine-powered shock from the ejecta into the preceding wind. The model requires total energetics of just few $\times 10^{50}$ ergs, slightly above the observed X-rays. We predict that the system is hydrogen poor.

We find a class of twisted and differentially rotating neutron star magnetospheres that do not have a light cylinder, generate no wind and thus do not spin-down. The magnetosphere is composed of embedded differentially rotating flux surfaces, with the angular velocity decreasing as $\Omega \propto 1/r$ (equivalently, becoming smaller at the foot-points closer to the axis of rotation). For each given North-South self-similar twist profile there is a set of self-similar angular velocity profiles (limited from above) with a "smooth", dipolar-like magnetic field structure extending to infinity. For spin parameters larger than some critical value, the light cylinder appears, magnetosphere opens up, and the wind is generated.

Fast Radio Bursts (FBRs) show highly different polarization properties: high/small RMs, high/small circular/linear fractions. We outline a complicated picture of polarization propagation in the inner parts of the magnetars' winds, at scales $\sim$ few to hundreds of light cylinder radii. The key point is the Faraday rotation of linear polarization in highly magnetized symmetric pair plasma, a $\propto B^2$ effect. Position angle (PA) rotation rate is maximal for propagation across the magnetic field and disappears only for parallel propagation. In the highly magnetized regime, $\omega \ll \omega_B$, it becomes independent of the magnetic field. Very specific properties of PA($\lambda$) (scaling of the rotation angle with the observed wavelength $\lambda$) can help identify/sort out the propagation effects. Two basic regimes in pair plasma predict PA $\propto \lambda$ and $\propto \lambda^3$ (depending on the magnetic dominance); both are different from the conventional plasma's PA = RM $ \lambda^2$. This is the main prediction of the model. A number of effects, all sensitive to the underlying parameters, contribute to the observed complicated polarization patterns: streaming of plasma along magnetic field lines near the light cylinder, Faraday depolarization, effects of limiting polarization, the associated effect of linear-circular conversion, and synchrotron absorption., Comment: This submission is superseded by arXiv:2205.13435

The radiation formation length for relativistic particles, $l_c \sim \gamma^2 \lambda$ ($\gamma$ is the Lorentz factor, $\lambda$ is the emitted wavelength), is much lager than the inter-particle distances in many astrophysical applications. This leads to the importance of plasma effects even for the high energy emission. The consequences are nontrivial: (i) averaging of the phases of the emitting particles reduces the power (a.k.a., a circle current does not emit); (ii) density fluctuations may lead to the sporadic production of coherent emission; (iii) plasma effects during assembly of a photon may lead to the suppression of the emission (Razin-Tsytovich effect for the superluminal modes), or, in the opposite limit of subluminous normal modes, to the newly discussed synchrotron super-radiance. For synchrotron emission the radiation formation length is the same for all emitted waves, $\sim c/\omega_B$ (non-relativistic Larmor length); for curvature emission it is $R/\gamma$ - macroscopically long in pulsar magnetospheres (e.g., kilometers for radio). The popular model of coherent curvature emission by bunches", with kilometers-long radiation formation length, particles swinging-out in a rotating magnetospheres before they finish emitting a wave, extreme requirements on the momentum spreads, and demands on the electric energy needed to keep the electrostatically repulsing charges together, all make that model internally inconsistent. Long radiation formation lengths affect how emission from PIC simulations should be interpreted: phases of the emitted wave should be added over the radiation formation length, not just the powers from the instantaneous acceleration of each particle.

We discuss dissipative processes occurring during production and escape of Fast Radio Bursts (FRBs) from magnetars' magnetospheres, the presumed loci of FRBs. High magnetic fields are required in the emission region, both to account for the overall energetics of FRBs, and in order to suppress ``normal'' (non-coherent) radiative losses of radio emitting particles; this limits the emission radii to $\leq {\rm few} \times 10 R_{NS}$. Radiative losses by particles in the strong FRB pulse may occur in the outer regions of the magnetosphere for longer rotation periods, $P\geq 1$ second. These losses are suppressed by several effects: (i) the ponderomotive pre-acceleration of background plasma along the direction of wave propagation (losses reduced approximately as $\gamma_\parallel^{3}$: smaller frequency, $ \propto \gamma_\parallel^2$ in power, and times scales stretched, $ \propto \gamma_\parallel$); this acceleration is non-dissipative and is reversed on the declining part of the pulse; (ii) Landau-Pomeranchuk-Migdal effects (long radiation formation length and ensuing destructive interference of scattered waves). In some cases an FRB pulse may be dissipated on external perturbations (e.g., an incoming pulse of Alfven waves): this may produce a pulse of UV/soft X-rays, a swan song of an FRB, possibly detectable by Chandra.

Magnetars, the likely sources of Fast Radio Bursts (FRBs), produce both steady highly relativistic magnetized winds, and occasional ejection events. We demonstrate that the requirement of conservation of the magnetic flux dominates the overall dynamics of magnetic explosions. This is missed in conventional hydrodynamic models of the ejections as expanding shell with parametrically added magnetic field, as well as one-dimensional models of magnetic disturbances. Most of the initial free energy of an explosion is actually spent on stretching its own internal magnetic field, while doing minimal $pdV$ work against the surrounding. Magnetic explosions from magnetars come into force balance with the pre-flares wind close to the light cylinder. They are then advected quietly with the wind, or propagate as electromagnetic disturbances. No powerful shock waves are generated in the wind.

We discuss a novel mechanism of axion production during scattering of Alfven waves by a fast moving Schwarzschild black holes. The process couples classical macroscopic objects, and effectively large amplitude electromagnetic (EM) waves, to microscopic axions. The key ingredient is that the motion of a black hole (BH) across magnetic field creates classical non-zero second Poincare invariant, the electromagnetic anomaly (Lyutikov 2011). In the case of magnetized plasma supporting Alfven wave, it is the fluctuating component of the magnetic field that contributes to the anomaly: for sufficiency small BH moving with the super-Alfvenic velocity the plasma does not have enough time to screen the parallel electric field. This creates time-dependent $ {\bf E} \cdot {\bf B} \neq 0$, and production of axions via the axion-EM coupling.

We consider topological configurations of the magnetically coupled spinning stellar binaries (e.g., merging neutron stars or interacting star-planet systems). We discuss conditions when the stellar spins and the orbital motion nearly `compensate' each other, leading to very {\it slow} overall winding of the coupled magnetic fields; slowly winding configurations allow gradual accumulation of magnetic energy, that is eventually released in a flare when the instability threshold is reached. We find that this slow winding can be global and/or local. We describe the topology of the relevant space $\mathbb{F}=T^1S^2$ as the unit tangent bundle of the two-sphere and find conditions for slowly winding configurations in terms of magnetic moments, spins and orbital momentum. These conditions become ambiguous near the topological bifurcation points; in certain cases they also depend on the relative phases of the spin and orbital motions. In the case of merging magnetized neutron stars, if one of the stars is a millisecond pulsar, spinning at $\sim$ 10 msec, the global resonance $\omega_1+\omega_2= 2 \Omega$ (spin-plus beat is two times the orbital period) occurs approximately a second before the merger; the total energy of the flare can be as large as $10\%$ of the total magnetic energy, producing bursts of luminosity $\sim 10^{44}$ erg s$^{-1}$. Higher order local resonances may have similar powers, since the amount of involved magnetic flux tubes may be comparable to the total connected flux.

We discuss constraints that the observed brightness temperatures impose on coherent processes in pulsars and Fast Radio Bursts (FRBs), and in particular on the hypothesis of coherent curvature emission by bunches. We estimate the peak brightness temperature that a bunch of charge $Ze$ can produce via synchrotron and/or curvature emission as $k_B T \sim (Z e)^2/\lambda$, where $\lambda$ is the typical emitted wavelength. We demonstrate that the bunch's electrostatic energy required to produce observed brightness temperature is prohibitively high, of the order of the total {\it bulk } energy. We compare corresponding requirements for the Free Electron Laser mechanism (Lyutikov 2021) and find that in that case the constraints are much easier satisfied.

We develop a model for the radio afterglow of the giant flare of SGR 1806-20 arising due to the interaction of magnetically-dominated cloud, an analogue of Solar Coronal Mass Ejections (CMEs), with the interstellar medium (ISM). The CME is modeled as a spheromak-like configuration. The CME is first advected with the magnetar's wind and later interacts with the ISM, creating a strong forward shock and complicated backwards exhaust flow. Using three-dimensional magnetohydrodynamic simulations, we study various relative configurations of the magnetic field of the CME with respect to the ISM's magnetic field. We show that the dynamics of the forward shock mostly follows the Sedov-Taylor blastwave, while the internal structure of the shocked medium is considerably modified by the back flow, creating a multiple shock configuration. We calculate synthetic synchrotron emissivity maps and light curves using two assumptions: (i) magnetic field compression; (ii) amplification of the magnetic field at the shock.We find that models with magnetic field amplification account better for the observed radio emission., Comment: 20 pages, 20 figures, accepted for publication in MNRAS

We consider magnetospheric structure of rotating neutron stars with internally twisted axisymmetric magnetic fields. The twist-induced and rotation-induced toroidal magnetic fields align/counter-align in different hemispheres. Using analytical and numerical calculations (with PHAEDRA code) we show that as a result the North-South symmetry is broken: the magnetosphere and the wind become "angled", of conical shape. Angling of the magnetosphere affects the spindown (making it smaller for mild twists), makes the return current split unequally at the Y-point, produces anisotropic wind and linear acceleration that may dominate over gravitational acceleration in the Galactic potential and give a total kick up to $\sim 100$ km/s. We also consider analytically the structure of the Y-point in the twisted magnetosphere, and provide estimate of the internal twist beyond which no stable solutions exist: over-twisted magnetospheres must produce plasma ejection events., Comment: Further analysis indicted that angling is a transient effect

We develop a model of the generation of coherent radio emission in the Crab pulsar, magnetars and Fast Radio Bursts (FRBs). Emission is produced by a reconnection-generated beam of particles via a variant of Free Electron Laser (FEL) mechanism, operating in a weakly-turbulent, guide-field dominated plasma. We first consider nonlinear Thomson scattering in a guide-field dominated regime, and apply to model to explain emission bands observed in Crab pulsar and in Fast Radio Bursts. We consider particle motion in a combined fields of the electromagnetic wave and thee lectromagnetic (Alfvenic) wiggler. Charge bunches, created via a ponderomotive force, Compton/Raman scatter the wiggler field coherently. The model is both robust to the underlying plasma parameters and succeeds in reproducing a number of subtle observed features: (i) emission frequencies depend mostly on the length $\lambda_t$ of turbulence and the Lorentz factor of the reconnection generated beam, $\omega \sim \gamma_b^2 ( c/\lambda_t) $ - it is independent of the absolute value of the underlying magnetic field. (ii) The model explains both broadband emission and the presence of emission stripes, including multiple stripes observed in the High Frequency Interpulse of the Crab pulsar. (iii) The model reproduces correlated polarization properties: presence of narrow emission bands in the spectrum favors linear polarization, while broadband emission can have arbitrary polarization. (iv) The mechanism is robust to the momentum spread of the particle in the beam. We also discuss a model of wigglers as non-linear force-free Alfven solitons (light darts).

The compact object in the interacting binary AR Sco has widely been presumed to be a rapidly rotating, magnetized white dwarf (WD), but it has never been detected directly. Isolating its spectrum has proven difficult because the spin-down of the WD generates pulsed synchrotron radiation that far outshines the WD's photosphere. As a result, a previous study of AR Sco was unable to detect the WD in the averaged far-ultraviolet spectrum from a Hubble Space Telescope (HST) observation. In an effort to unveil the WD's spectrum, we reanalyze these HST observations by calculating the average spectrum in the troughs between synchrotron pulses. We identify weak spectral features from the previously unseen WD and estimate its surface temperature to be 11500$\pm$500K. Additionally, during the synchrotron pulses, we detect broad Lyman-$\alpha$ absorption consistent with hot WD spectral models. We infer the presence of a pair of hotspots, with temperatures between 23000K and 28000K, near the magnetic poles of the WD. As the WD is not expected to be accreting from its companion, we describe two possible mechanisms for heating the magnetic poles. The Lyman-$\alpha$ absorption of the hotspots appears relatively undistorted by Zeeman splitting, constraining the WD's field strength to be 100 MG, but the data are insufficient to search for the subtle Zeeman splits expected at lower field strengths., Comment: 10 pages, 7 figures, accepted for publication in the Astrophysical Journal

Afterglows of gamma-ray bursts often show flares, plateaus, and sudden intensity drops: these temporal features are difficult to explain as coming from the forward shock. We calculate radiative properties of early GRB afterglows with the dominant contribution from the reverse shock (RS) propagating in an ultra-relativistic (pulsar-like) wind produced by the long-lasting central engine. RS emission occurs in the fast cooling regime -- this ensures high radiative efficiency and allows fast intensity variations. We demonstrate that: (i) mild wind power, of the order of $\sim 10^{46}$ erg s$^{-1}$, can reproduce the afterglows' plateau phase; (ii) termination of the wind can produce sudden steep decays; (iii) mild variations in the wind luminosity can produce short-duration afterglow flares., Comment: arXiv admin note: substantial text overlap with arXiv:2004.13600

We discuss coherent free electron laser (FEL) operating during explosive reconnection events in magnetized pair plasma of magnetar magnetospheres. The model explains many salient features of Fast Radio Bursts/magnetars' radio emission: temporal coincidence of radio and high energy bursts, high efficiency of conversion of plasma kinetic energy into coherent radiation, presence of variable, narrow-band emission features drifting down in frequency, high degree of linear polarization. The model relies on magnetar-specific drifting $e^\pm$ plasma components (which generate wiggler field due to the development of the firehose instability) and the presence of reconnection-generated particle beam with mild Lorentz factor of $\gamma_b \sim$ few hundred.

We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - trying to lower its magnetic potential energy. As a result a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on few Alfv\'en time scales. In the case of higher order (second order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfv\'en waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetars' giant flares, fast development of tilting instabilities, and no stabilization of the higher order spheromaks, make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds., Comment: 20 pages, 10 figures, accepted for publication in Journal of Plasma Physics

We construct a turbulent model of the Crab Nebula's non-thermal emission. The present model resolves a number of long-standing problems of the Kennel-Coroniti (1984) model: (i) the sigma problem; (ii) the hard spectrum of radio electrons; (iii) the high peak energy of gamma-ray flares; (iv) and the spacial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I accelerated at the wind termination shock via Fermi-I mechanism, and Component-II accelerated in reconnecting turbulence in highly magnetized ($\sigma$ $\gg 1$) plasma in the central part of the Crab Nebula. The reconnecting turbulence Component-II extends from radio to gamma rays: it accelerate radio electrons with a hard spectrum, destroy the large scale magnetic flux (and thus resolves the sigma-problem), and occasionally produces gamma-ray flares (from the largest scale reconnection events). The model reproduces the broad-band spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as spatially resolved evolution of the spectral indices in IR and optical bands.

Lyutikov (2002) predicted "radio emission from soft gamma-ray repeaters (SGRs) during their bursting activity". Detection of a Mega-Jansky radio burst in temporal coincidence with high energy bursts from a Galactic magnetar SGR 1935+2154 confirms that prediction. Similarity of this radio event with Fast Radio Bursts (FRBs) suggests that FRBs are produced within magnetar magnetospheres. We demonstrate that SGR 1935+2154 satisfies the previously derived constraints on the physical parameters at the FRBs' loci. Coherent radio emission is generated in the inner parts of the magnetosphere at $r< 100 R_{\rm NS}$. The radio emission is produced by the yet unidentified plasma emission process, occurring during the initial stages of reconnection events.

We develop a model of early X-ray afterglows of gamma-ray bursts originating from the reverse shock (RS) propagating through ultra-relativistic, highly magnetized pulsar-like winds produced by long-lasting central engines. We first perform fluid and MHD numerical simulations of relativistic double explosions. We demonstrate that even for constant properties of the wind a variety of temporal behaviors can be produced, depending on the energy of the initial explosion and the wind power, the delay time for the switch-on of the wind, and magnetization of the wind. X-ray emission of the highly magnetized RS occurs in the fast cooling regime - this ensures high radiative efficiency and allows fast intensity variations. We demonstrate that: (i) RS emission naturally produces light curves showing power-law temporal evolution with various temporal indices; (ii) mild wind power, of the order of $\sim 10^{46}$ erg s$^{-1}$ (equivalent isotropic), can reproduce the afterglows' plateau phase; (iii) termination of the wind can produce sudden steep decays; (iv) short-duration afterglow flares are due to mild variations in the wind luminosity, with small total injected energy., Comment: 19 pages, 21 figures

We develop a model of the white dwarf (WD) - red dwarf (RD) binaries AR Sco and AE Aqr as systems in a transient propeller stage of highly asynchronous intermediate polars. The WDs are relatively weakly magnetized with magnetic field of $\sim 10^6$ G. We explain the salient observed features of the systems due to the magnetospheric interaction of two stars. Currently, the WD's spin-down is determined by the mass loading of the WD's magnetosphere from the RD's at a mild rate of $\dot{M}_{WD} \sim 10^{-11} M_\odot $/yr. Typical loading distance is determined by the ionization of the RD's wind by the WD's UV flux. The WD was previously spun up by a period of high accretion rate from the RD via Roch lobe overflow with $\dot{M} \sim 10^{-9} M_\odot $/yr, acting for as short a period as tens of thousands of years. The non-thermal X-ray and optical synchrotron emitting particles originate in reconnection events in the magnetosphere of the WD due to the interaction with the flow from the RD. In the case of AR Sco, the reconnection events produce signals at the WD's rotation and beat periods - this modulation is due to the changing relative orientation of the companions' magnetic moments and resulting variable reconnection conditions. Radio emission is produced in the magnetosphere of the RD, we hypothesize, in a way that it is physically similar to the Io-induced Jovian decametric radiation.

Pulsars moving through ISM produce bow shocks detected in hydrogen H$\alpha$ line emission. The morphology of the bow shock nebulae allows one to probe the properties of ISM on scales $\sim 0.01$ pc and smaller. We performed 2D RMHD modeling of the pulsar bow shock and simulated the corresponding H$\alpha$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $\delta\rho/\rho \sim 1$, leads to significant variations of the shape of the shock seen in H$\alpha$ line emission. We successfully reproduce the morphology of the Guitar Nebula. We infer quasi-periodic density variations in the warm component of ISM with a characteristic length of $\sim0.1$~pc. Structures of this scale might be also responsible for the formation of the fine features seen at the forward shock of Tycho SNR in X-rays. Formation of such short periodic density structures in the warm component of ISM is puzzling, and bow-shock nebulae provide unique probes to study this phenomenon., Comment: 12 pages, 7 figures, accepted MNRAS

Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $\eta = L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $\eta \leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars., Comment: 13 pages, 4 figures

Intense radiation field can modify plasma properties, the corresponding refractive index, and lead to such nonlinear propagation effects as self-focusing. We estimate the corresponding effects in pair plasma, both in unmagnetized and strongly magnetically dominated case. First, in the unmagnetized pair plasma the ponderomotive force does not lead to charge separation, but to density depletion. Second, for astrophysically relevant plasmas of pulsar magnetospheres, (and possible loci of Fast Radio Bursts), where cyclotron frequency $\omega_B$ dominates over plasma frequency $\omega_p$ and the frequency of the electromagnetic wave, $\omega_B \gg \omega_p,\, \omega$, we show that (i) there is virtually no nonlinearity due to changing effective mass in the field of the wave; (ii) ponderomotive force is $F_p^{(B)} =- {m_e c^2}/({4 B_0^2}) \nabla E^2$; it is reduced by a factor $(\omega/\omega_B)^2$ if compared to the unmagnetized case ($B_0$ is the external magnetic field and $E$ is the electric field of the wave); (iii) for radiation beam propagating along constant magnetic field in pair plasma with density $n_\pm$, the ponderomotive force leads to appearance of circular currents that lead to the decrease of the field within the beam by a factor $\Delta B/B_0 = 2\pi n_\pm m_e c^2 {E^2}/{B_0^4}$. Applications to the physics of FRBs are discussed; we conclude that for parameters of FRB's the dominant magnetic field completely suppresses nonlinear radiation effects.

We find a new family of solutions for force-free magnetic structures in cylindrical geometry. These solutions have radial power-law dependance and are periodic but non-harmonic in azimuthal direction; they generalize the vacuum $z$-independent potential fields to current-carrying configurations.

Mergers of white dwarfs (WDs) may lead to a variety of transient astrophysical events, SNIa being one possible outcome. Lyutikov & Toonen (2017, 2019) argued that mergers of WDs result, under various parameter regimes, in unusual central engine-powered supernova and a type of short Gamma Ray Bursts that show extended emission tails. Observations by Gvaramadze et al. (2019) of the central star and the nebula J005311 match to the details the model of Lyutikov & Toonen (2017, 2019) for the immediate product of a merger of a heavy ONeMg WD with CO WD (age, luminosity, stellar size, hydrogen deficiency and chemical composition)., Comment: Contribution to "COMPACT WHITE DWARF BINARIES", Yerevan, Armenia

We consider conditions for jet break-out through ejecta following mergers of neutron stars and provide simple relations for the break out conditions. We demonstrate that: (i) break-out requires that the isotropic-equivalent jet energy $E_j$ exceeds the ejecta energy $E_{ej}$ by $E_j \geq E_{ej}/ \beta_0$, where $ \beta_0 = V_{ej}/c$, $V_{ej}$ is the maximum velocity of the ejecta. If the central engine terminates before the break out, the shock approaches the edge of the ejecta slowly $\propto 1/t$; late break out occurs only if at the termination moment the head of the jet was relatively close to the edge. (ii) If there is a substantial delay between the ejecta's and the jet's launching, the requirement on the jet power increases. (iii) The forward shock driven by the jet is mildly strong, with Mach number $M\approx 5/4$ (increasing with time delay $t_d$); (iii) the delay time $t_d$ between the ejecta and the jet's launching is important for $t_d > t_0= ({3 }/{16} ) {c M_{ej} V_{ej}}/{L_j} = 1.01 {\rm sec} M_{ej, -2} L_{j, 51} ^{-1} \left ( {\beta_{ej}} /{0.3} \right)$, where $M_{ej}$ is ejecta mass, $L_j$ is the jet luminosity (isotropic equivalent). For small delays, $t_0 $ is also an estimate of the break-out time.

We build a model of radius-to-frequency mapping in magnetospheres of neutron stars and apply it to frequency drifts observed in Fast Radio Bursts. We assume that an emission patch propagates along the dipolar magnetic field lines producing coherent emission with frequency, direction and polarization defined by the local magnetic field. The observed temporal evolution of the frequency depends on relativistic effects of time contraction and the curvature of the magnetic field lines. The model generically produces linear scaling of the drift rate, $\dot{\omega} \propto - \omega$, matching both numerically and parametrically the rates observed in FBRs; a more complicated behavior of $\dot{\omega} $ is also possible. Fast rotating magnetospheres produce higher drifts rates for similar viewing parameters than the slowly rotating ones. In the case of repeaters same source may show variable drift pattens depending on the observing phase. We expect rotational of polarization position angle through a burst, though by smaller amount than in radio pulsars. All these findings compare favorably with properties of FBRs, strengthening their possible loci in the magnetospheres of neutron stars.

We interpret recent observations of high-to-low frequency drifting features in the spectra of the repeating FRBs as evidence of sharply changing plasma properties in the emission region, presumably the neutron stars magnetospheres. The drifts are then FRBs' analogues of radius-to-frequency mapping in pulsars and Solar type-III radio burst (but not in a sense of a particular emission mechanism). The drifts rates of $\sim 100$ MHz ms$^{-1}$ at frequencies $\sim$ GHz translate to physical size of $\sim {c \omega}/{\dot{\omega}} \sim$ few $\times 10^8$ cm, matching the hypothesis of the FRB origin in the magnetospheres of neutron stars. We suggest that reconnection events result in generation of upward propagating plasma beams that produce radio emission with frequency related to the decreasing local magnetic field and plasma density., Comment: arXiv admin note: substantial text overlap with arXiv:1901.03260

We suggest that narrow, long radio filaments near the Galactic Center arise as kinetic jets - streams of high energy particles escaping from ram-pressure confined pulsar wind nebulae (PWNe). The reconnection between the PWN and interstellar magnetic field allows pulsar wind particles to escape, creating long narrow features. They are the low frequency analogs of kinetic jets seen around some fast-moving pulsars, such as The Guitar and The Lighthouse PWNe. The radio filaments trace a population of pulsars also responsible for the Fermi GeV excess produced by the Inverse Compton scattering by the pulsar wind particles. The magnetic flux tubes are stretched radially by the large scale Galactic winds. In addition to PWNe accelerated particles can be injected at supernovae remnants. The model predicts variations of the structure of the largest filaments on scales of $\sim$ dozens of years - smaller variations can occur on shorter time scales. We also encourage targeted observations of the brightest sections of the filaments and of the related unresolved point sources in search of the powering PWNe and pulsars., Comment: 5 pages, 2 figures

We discuss physical constrains that observations of high brightness temperature coherent radio emission, with brightness temperatures as high as $T_b \sim 10^{35}$ K, impose on the plasma parameters at relativistically moving astrophysical sources. High brightness temperatures imply a minimal plasma energy density at the source. Additional important constraints come from the fact that resonantly emitting particles lose most of their energy to non-resonant inverse Compton and synchrotron processes. We also interpret recent observations of high-to-low frequency drifting features in the spectra of repeating FRBs as analogues of type-III Solar radio bursts produced by reconnection plasma beams within magnetospheres of highly magnetized neutron stars.

We suggest that fast-rising blue optical transients (FBOTs) and the brightest event of the class AT2018cow result from an electron-capture collapse to a \NS\ following a merger of a massive ONeMg white dwarf (WD) with another WD. Two distinct evolutionary channels lead to the disruption of the less massive WD during the merger and the formation of a shell burning non-degenerate star incorporating the ONeMg core. During the shell burning stage a large fraction of the envelope is lost to the wind, while mass and angular momentum are added to the core. As a result, the electron-capture collapse occurs with a small envelope mass, after $\sim 10^2-10^4$ years. During the formation of a neutron star as little as $\sim 10^{-2} M_\odot $ of the material is ejected at the bounce-off with mildly relativistic velocities and total energy $\sim$ few $ 10^{50}$ ergs. This ejecta becomes optically thin on a time scale of days - this is the FBOT. During the collapse, the neutron star is spun up and magnetic field is amplified. The ensuing fast magnetically-dominated relativistic wind from the newly formed neutron star shocks against the ejecta, and later against the wind. The radiation-dominated forward shock produces the long-lasting optical afterglow, while the termination shock of the relativistic wind produces the high energy emission in a manner similar to Pulsar Wind Nebulae. If the secondary WD was of the DA type, the wind will likely have $\sim 10^{-4} M_\odot$ of hydrogen; this explains the appearance of hydrogen late in the afterglow spectrum. The model explains many of the puzzling properties of FBOTs/AT2018cow: host galaxies, a fast and light anisotropic ejecta producing a bright optical peak, afterglow high energy emission of similar luminosity to the optical, and late infra-red features., Comment: arXiv admin note: text overlap with arXiv:1709.02221

We obtained high temporal resolution spectroscopy of the unusual binary system AR Sco covering nearly an orbit. The H$\alpha$ emission shows a complex line structure similar to that seen in some polars during quiescence. Such emission is thought to be due to long-lived prominences originating on the red dwarf. A difference between AR Sco and these other systems is that the white dwarf in AR Sco is rapidly spinning relative to the orbital period. "Slingshot" prominences stable at 3 to 5 stellar radii require surface magnetic fields between 100 and 500 G. This is comparable to the estimated WD magnetic field strength near the surface of the secondary. Our time-resolved spectra also show emission fluxes, line equivalent widths, and continuum color varying over the orbit and the beat/spin periods of the system. During much of the orbit, the optical spectral variations are consistent with synchrotron emission with the highest energy electrons cooling between pulses. On the time-scale of the beat/spin period we detect red and blue-shifted H$\alpha$ emission flashes that reach velocities of 700 km/s. Red-shifted Balmer emission flashes are correlated with the bright phases of the continuum beat pulses while blue-shifted flashes appear to prefer the time of minimum in the beat light curve. We propose that much of the energy generated in AR Sco comes from fast magnetic reconnection events occurring near the inward face of the secondary and we show that the energy generated by magnetic reconnection can account for the observed excess luminosity from the system., Comment: Accepted for publication in ApJ, December 27, 2018, 15 page, 14 figures

We reconsider the structure of thermally driven rotating Parker wind. Rotation, without \Bf, changes qualitatively the structure of the subsonic region: solutions become non-monotonic and do not extend to the origin. For small angular velocities solutions have two critical points - X-point and O-points, which merge at the critical angular velocity of the central star $\Omega_{crit} = G M_\ast/(2 \sqrt{2} c_s R_{b}^2)$ (where $M_\ast$ and $R_{b}$ are mass and radius of the central star, $c_s$ is the sound speed in the wind). For larger spins there is no critical points in the solution. For disk winds (when the base of the wind rotates with Keplerian velocity) launched equatorially the coronal sound speed should be smaller than $\approx 0.22 v_K$ in order to connect to the critical curve ($v_K$ is the Keplerian velocity at a given location on the disk).

We outline a model of the Crab Pulsar Wind Nebula with two different populations of synchrotron emitting particles, arising from two different acceleration mechanisms: (i) Component-I due to Fermi-I acceleration at the equatorial portion of the termination shock, with particle spectral index $p_I \approx 2.2$ above the injection break corresponding to $\gamma_{wind} \sigma_{wind} \sim 10^5$, peaking in the UV ($\gamma_{wind} \sim 10^2$ is the bulk Lorentz factor of the wind, $\sigma _{wind} \sim 10^3$ is wind magnetization); (ii) Component-II due to acceleration at reconnection layers in the bulk of the turbulent Nebula, with particle index $p_{II} \approx 1.6$. The model requires relatively slow but highly magnetized wind. For both components the overall cooling break is in the infra-red at $\sim 0.01$ eV, so that the Component-I is in the fast cooling regime (cooling frequency below the peak frequency). In the optical band Component-I produces emission with the cooling spectral index of $\alpha_o \approx 0.5$, softening towards the edges due to radiative losses. Above the cooling break, in the optical, UV and X-rays, Component-I mostly overwhelms Component-II. We hypothesize that acceleration at large-scale current sheets in the turbulent nebula (Component-II) extends to the synchrotron burn-off limit of $\epsilon_s \approx 100$ MeV. Thus in our model acceleration in turbulent reconnection (Component-II) can produce both hard radio spectra and occasional gamma-ray flares. This model may be applicable to a broader class of high energy astrophysical objects, like AGNe and GRB jets, where often radio electrons form a different population from the high energy electrons., Comment: shortened version of the initial submission, accepted by MNRAS

We consider electromagnetic interaction and precursor emission of merging neutron stars. Orbital motion of the magnetized neutron stars may revive pair production within the common magnetosphere years before the merger, igniting pulsar-like magnetospheric dynamics. We identify two basic scenarios: (i) only one star is magnetized (1M-DNS scenario) and (ii) both stars are magnetized (2M-DNS scenario). Inductively created electric fields can have component along the total magnetic field (gaps) and/or the electric field may exceed the value of the local magnetic field. The key to the detection is orbital modulation of the emission. If only one star is magnetized (1M-DNS scenario) the emission is likely to be produced along the direction of the magnetic field at the location of the secondary; then, if the magnetic axis is misaligned with the orbital spin, this direction is modulated on the orbital period. For the 2M-DNS scenario, the structure of the common magnetosphere of the non-rotating neutron stars is complicated, with gaps, but no $E>B$ regions; there is strong orbital variations for the case of misaligned magnetic moments. For the same parameters of neutron stars the 2M-DNS scenario has intrinsically higher potential than the 1M-DNS one. The overall powers are not very high, $\leq 10^{45} $ erg s$^{-1}$; the best chance to detect electromagnetic precursors to the merging neutron stars is if the interaction of their magnetospheres leads to the production of pulsar-like coherent radio emission modulated at the orbital period, with luminosity of up to $\sim 1$ Jankys at the time the merger., Comment: This ArXive version corrects a sign error in the published paper in Eq. (15). Equations (15)-(18) and figures 2 and 3 are affected. Overall results and conclusions remain the same