Radio emission from SN 1181 hosting a white dwarf merger product

The remnant of the historical supernova 1181 is claimed to be associated with a white dwarf merger remnant J005311. The supernova remnant (SNR) shock, and a termination shock expected to be formed by the intense wind of J005311, are potential sites for radio emission via synchrotron emission from shock-accelerated electrons. In this paper, we estimate the radio emission from these two shocks, and find the peak radio flux to be 0.1--10 mJy (at 0.01--1 GHz) in the outer SNR shock and 0.01--0.1 mJy (at 1--10 GHz) in the inner termination shock. We also search for radio emission from this source in the archival data of the Karl G. Jansky Very Large Array (VLA) Sky Survey at 3 GHz, NRAO VLA Sky Survey at 1.4 GHz and the Canadian Galactic Plane Survey at 408 MHz, resulting in no significant detection. While targeted observations with higher sensitivity are desired, we particularly encourage those at higher frequency and angular resolution to probe the inner termination shock and its evolution.


Introduction
Mergers of binary white dwarfs ( WDs) evolve in different paths depending on the compositions and masses of the two WDs.Such mergers are thought to result in Type Ia supernovae ( SNe) ( Webbink 1984 ;Iben & Tutukov 1984 ) and/or neutron stars via accretion-induced collapse ( Saio & Nomoto 1985 ;Taam & van den Heuvel 1986 ) , but some can leave behind massive rapidly rotating WDs ( Pshirkov et al. 2020 ;Caiazzo et al. 2021 ) .The evolutionary paths of such WD merger products have started to be explored in detail ( Schwab et al. 2016 ;Zhang et al. 2019 ;Wu et al. 2023 ;Yao et al. 2023 ;Zhong et al. 2024 ) .
In 2019, an infrared nebula IRAS 00500 + 6713, hosting a massive WD J005311 at the center, was discovered ( Gvaramadze et al. 2019 ) .Spectroscopic observations revealed an intense wind blowing from the WD, with velocity v wind ∼ 16000 km s −1 and mass-loss rate Ṁwind ∼ 10 −6 M yr −1 ( Gvaramadze et al. 2019 ;Lykou et al. 2023 ) , supporting that this WD is a merger product ( see also Kashiyama et al. 2019 ) .Interestingly, the expansion rate and extent of the nebula, as well as its position on the celestial sphere and distance of d = 2.3 kpc ( Bailer-Jones et al. 2021 ) , makes J005311 the prime candidate for the remnant of the historical SN 1181 ( Ritter et al. 2021 ) .
Furthermore, X-ray data of this object obtained by XMM-Newton and Chandra and analyzed by Oskinova et al. ( 2020 ) and Ko et al. ( 2023 ) revealed that there is a more intense X-ray source at the center ( of intensity ∼10 −13 erg cm −2 s −1 arcsec −2 ) in addition to the diffuse emission ( of ∼10 −17 erg cm −2 s −1 arcsec −2 ) likely powered by the merger ejecta.While the XMM observations did not resolve the compact inner source, Chandra observations with a better angular resolution found that the inner source has a finite extent ( 0.0087-0.018pc) .Ko et al. ( 2023 ) constructed a theoretical model that comprehensibly explains the multiple layers observed in various wavelengths ( see their schematic figure and figure 1 ) .They have newly suggested that the inner X-ray emission is powered by a termination shock, formed by the interaction between the unshocked SN ejecta and the fast wind from the central WD.They also confirm that the outer diffuse region is powered by the interaction between the SN 1181 ejecta and the surrounding interstellar medium ( ISM) , like SN remnants ( SNRs) .Using their model, they have succeeded in explaining the multi-wavelength observation properties, if the wind started blowing recently within the past 30 years.
SNRs and young SNe within years of explosion emit synchrotron radiation in radio, as the collisionless shocks formed by the interaction can accelerate particles and amplify the ambient magnetic field ( e.g., Weiler et al. 2002 ;Dubner & Giacani 2015 ;Chevalier & Fransson 2017 ) .This is likely the case for SN 1181 as well, but radio emission from this source has yet to be reported.Thus existing and future radio observations may offer a strong test for the model.
In this work, we consider radio emission from both the outer and inner shocked regions, which are schematically  Ko et al. ( 2023 ) .We expect synchrotron radiation from the outer SN shock sweeping the interstellar medium ( as in typical SNRs) , and the WD wind termination shock sweeping the inner SN ejecta.
shown in figure 1 .The radio source of interest is the region swept by the outer forward shock ( as in typical SNRs) at θ out ≈ 2 , and the inner region swept by the wind termination shock at θ in ≈ 1 -2 .We first search for such signals in archival radio data that had covered this object.We next model the expected radio signal, and find that the radio flux in these regions may be bright enough to be detectable with deeper targeted follow-up by current telescopes.
This work is structured as follows.In section 2 we present the archival radio observations and their results.Since no significant radio emission was found, we present upper limits to the radio flux for the inner and outer shocked regions.In section 3 we describe our model of the radio synchrotron emission.In section 4 we summarize our results, and mention future observability for the inner region, based on the time evolution of the emitting region and its radio luminosity.

Constraints from archival data
To constrain the radio emission, we search for radio observational data in the archive of the NRAO VLA Sky Survey ( NVSS; Condon et al. 1998 ) , the Very Large Array Sky Survey ( VLASS; Lacy et al. 2020 ) , and the Canadian Galactic Plane Survey ( CGPS; Taylor et al. 2003 ;Tung et al. 2017 ) .We specifically search for emission from the inner region with VLASS that has a smaller beam size of a few arcsec, while we mainly probe the outer region for the other two surveys with larger beam sizes.NVSS is a 1.4 GHz continuum survey covering the entire sky north of ∼40 • declination.The region containing SN 1181 was observed in 1993 November with a synthesized beam size of 45 ( full width at half maximum; FWHM) .The rms noise level around the region is 0.48 mJy beam −1 .We find an emission with a peak flux density of 1.8 mJy at ∼50 east of the location of SN 1181 ( figure 2 ) .The source has a size of a few arcmin, most consistent with the outer shocked region ( SNR) with a diameter of ∼240 ( Oskinova et al. 2020 ;Ko et al. 2023 ) .Since the emission is marginal ( signal-to-noise ratio of 3.7) , we treat it as a nondetection and discuss it using the 3 σ upper limit.VLASS is the third radio survey project of the VLA designed to cover the sky visible at the VLA with a frequency bandwidth of 2-4 GHz.The region containing SN 1181 was observed three times: 2017 November, 2020 August, and 2023 March, with a synthesized beam size of 3 .1 × 2 .2, 3 .4 × 2 .2, and 3 .1 × 2 .4, respectively.The typical rms noise level of the Quick Look images is 0.14 mJy beam −1 ( Gordon et al. 2021 ) .We do not find any significant emission.Since the largest angular scale structure of the S -band B-array observations is approximately 58 , which is smaller than the spatial extents of the outer shocked region ( 100 -120 ) , it is possible that the emission was resolved out in the observations.CGPS observations at 408 MHz were conducted from 1995 to 2009 with the Synthesis Telescope at the Dominion Radio Astrophysical Observatory ( Tung et al. 2017 ) .The typical angular resolution is 2 .8 × 2 .8 cosec δ, and the typical sensitivity is 3 mJy beam −1 rms.No significant emission is found.

Radio emission modeling
In a collisionless shock, we expect that magnetic fields are amplified and electrons are accelerated through diffusive shock acceleration.In this section we model the expected radio emission from the inner and outer shocked regions.

Synchrotron emission/absorption modeling
We model each shocked region as one-zone and parameterize the characteristic magnetic field strength B in and B out at each region, using an amplification efficiency B as Hereafter i = {in, out} refers to quantities in the inner and outer shocked region, respectively.Here ρ i and v i are the upstream density and upstream velocity ( in the shock rest frame) in region i , respectively.We fix the value of B as B = 0.01, a typical value inferred for both SNe and SNRs ( e.g., Chevalier & Fransson 2006 ;Maeda 2012 ;Murase et al. 2019 ;Reynolds et al. 2021 ) .
The density of the WD wind powering the inner shocked region can be expressed as ρ wind = Ṁwind / (4 π r 2 in v wind ) , where Ṁwind and v wind are, respectively, the mass-loss rate and velocity of the wind, and r in is the radius of the inner shock front.For the outer shocked region, the relevant parameters ρ out and v out correspond to the ISM density and the SN forward shock velocity.From spectral analysis ( Gvaramadze et al. 2019 ;Lykou et al. 2023 ) and dynamical modeling ( Ko et al. 2023 ) of this object, the values of these parameters in 2021 A.D. were estimated to be those shown in table 1 .Here, since the termination shock velocity ( ∼800 km s −1 from our modeling) is much smaller than the wind velocity of v wind ≈ 1.5 × 10 4 km s −1 , we set v in = v wind .Here we adopt the fiducial case of the model; (t age , in , M ej ) = (20 yr , 0 .35 M ) .
Substituting values from table 1 into equation ( 1) , the magnetic field strength at each region can be estimated as ( 2) The maximum energy E max = γ max m e c 2 of the electrons, where γ max is the corresponding Lorentz factor, is set following Tsuna and Kawanaka ( 2019 ) by requiring the acceleration timescale ( assuming the Bohm limit) t acc ,i ≈ (20 / 3) γ i m e c 3 / (eB i v 2 i ) to be smaller than the age t age,i and the synchrotron cooling timescale Here m e , e , c , and σ T are, respectively, the electron mass, the elementary charge, the speed of light, and the Thomson crosssection.A plausible range of t age,in = 12-30 yr is obtained from modeling of the inner shocked region ( Ko et al. 2023 ) , while t age,out ≈ 840 yr is estimated from the SN date of 1181.The corresponding maximum frequencies of synchrotron radiation ν max = 3 γ 2 max eB i / 16 m e c are ν max,in ∼ 2 × 10 19 Hz and ν max,out ∼ 7 × 10 16 Hz, which are much higher than the radio frequency of interest.Hence we expect that synchrotron emission spectra from accelerated electrons can easily be extended to the GHz band.
Downloaded from https://academic.oup.com/pasj/advance-article/doi/10.1093/pasj/psae023/7642730 by University of Tokyo Library user on 16 April 2024 Fig. 3. Radio spectra for selected parameter sets in table 2 .The blue and orange lines refer to the emissions from the inner and outer shocked regions, respectively.The solid, dashed, and dotted lines represent, respectively, the optimistic, fiducial, and pessimistic parameter sets.The triangles show the 3 σ upper limits ( with NVSS being a tentative detection) , color-coded as constraints on the outer/inner region ( for details see section 2 ) .The horizontal gray line shows the 5 σ sensitivity for Jansky VLA in the X-band ( 8-12 GHz) , for an integration time of one hour.Note that the sensitivity applies to the total flux density when not spatially resolved, and a higher sensitivity is required when spatially resolved.
The relativistic electrons are assumed to follow a power-law energy distribution N( E) ∝ E −p ( p > 2) , expected for diffusive shock acceleration.For an isotropic electron velocity distribution, the synchrotron power per frequency for each shocked region i is ( Rybicki & Lightman 1979 ) : where is the gamma function.Here C i is the normalization constant determined by the fraction of total dissipated energy given to relativistic electrons e as where we assume v i is constant for simplicity and adopt the present-day value.We adopt a minimum electron energy E min = 2 m e c 2 .For the power-law index p , while a range of p = 2-3 is adopted for an SNR in general, p ≈ 3 is found from radio modeling for young SNRs ( e.g., Chevalier & Fransson 2006 ;Maeda 2012 ) .Here we consider p = 2.5 and p = 3.It can be shown as follows that synchrotron cooling does not affect the radio emission.Finding the Lorentz factor γ cool,i , where t cool,i is equal to the age, the frequencies ν c,i = 3 γ 2 cool ,i eB i / (16 m e c ) beyond which cooling is important can be estimated for both regions as At lower frequencies, the synchrotron emission is subject to synchrotron self-absorption.The absorption coefficient for a power-law distribution of electron energy [equation ( 5) ] is given as ( Rybicki & Lightman 1979 ) : and the optical depth is obtained as τ sa,i ( ν) = α sa,i ( ν) r i .
Taking into account these absorption processes, the observed specific luminosity can be expressed as The luminosities and fluxes ( assuming the distance of 2.3 kpc) for sets of model parameters ( p , e ) are shown in figure 3 .The solid, dashed and dotted lines represent the optimistic, fiducial and pessimistic cases, respectively, which vary the values of ( p , e ) as in table 2 .The values for the fiducial cases are inspired from previous modeling of SNe and SNRs ( e.g., Maeda 2012 ;Reynolds et al. 2021 ) , and those for the optimistic cases are set to roughly explain the flux of the marginal finding from the NVSS image in section 2 .The upper limits obtained from the previous surveys introduced in section 2 are shown as triangles.We find that for the outer shocked region the peak flux is about 0.1-10 mJy at Downloaded from https://academic.oup.com/pasj/advance-article/doi/10.1093/pasj/psae023/7642730 by University of Tokyo Library user on 16 April 2024  frequencies around 0.1-1 GHz.For the inner shocked region, the peak flux is about 0.01-1 mJy around 1-10 GHz.Synchrotron self-absorption determines the rise at frequencies lower than the peak.The flux becomes maximum when the optical depth is around unity.The frequency at the peak ( ν max,i ) can be described as ( Substituting them into equations ( 9) and ( 10 ) , we obtain the maximum flux: (p = 3) .

Free-free emission/absorption modeling
As the shocked regions contain charged particles due to shock heating/cooling, they may contribute to both free-free emission and absorption at radio frequencies.Here we infer free-free emission/absorption from three regions-the outer shocked region, the inner shocked WD wind and the inner shocked ejecta-and find that they are likely negligible in our situation.
For the shocked wind region, the flux of free-free emission is given ( in cgs units) as ( Rybicki & Lightman 1979 ) : where n e,in,wind is the electron number density and Z ion,in and n ion,in,wind are the charge and number density of each ion, and T e,in,wind is the electron temperature of the inner shocked wind.
As the inner shocked wind is responsible for the inner X-ray emission ( Ko et al. 2023 ) also negligible compared to the synchrotron component.Here, we adopt Z ion,out = 1, since the abundance of the outer shocked region is near the solar abundance and hydrogen and helium are expected to be fully ionized ( Ko et al. 2023 ) .
The free-free absorption coefficient can be expressed as ( Mezger & Henderson 1967 ) where n e is the electron number density and Z ion and n ion are the charge and the number density of each ion.Evaluating the equation as was done for the free-free emission, we find that free-free absorption from these two regions are unimportant, with optical depths ( ∼α ff,i r i ) of both 1 in the radio band.The estimation for the inner shocked ejecta is less straightforward than the other two regions, as the electron temperature is expected to not exceed 1-a few eV, and thus its cooling emission would not contribute to the observed X-rays ( see appendix A of Ko et al. 2023 ) .The cooling timescale in this region is 3-4 orders of magnitude shorter than the age, and the region is expected to further compress into a thin shell by radiative cooling.This runaway cooling should continue until the shocked region recombines, and we may thus expect the free-free emission/absorption in this region to be negligible for an extremely small ionization fraction of ∼10 −4 -10 −3 .
So far we have neglected the possibility that the shocked ejecta may become photo-ionized by the central WD that has a high effective temperature of ∼2 × 10 5 K ( Gvaramadze et al. 2019 ;Lykou et al. 2023 ) .The optically thick wind from the WD may reprocess the ionizing photons from the central star and strongly suppress its flux, as found in models of optically thick winds from Wolf-Rayet stars ( e.g., Smith et al. 2002 ;Götberg et al. 2017 ;Sander & Vink 2020 ) .However, the situation can be complicated by the possible clumpiness of the wind and the ejecta, as well as the uncertainies in their composition and geometry.
Here we adopt a more empirical approach, and use the nondetection from the VLASS survey in section 2 to infer the ionization of the shocked ejecta.The radius, mass and density of the shocked ejecta are obtained by one-zone modeling ( Ko et al. 2023 ) of a WD wind, with mass-loss rate of Ṁw , sweeping the bound part of the SN ejecta with a total mass of M ej and following a density profile ρ ∝ r −1.5 at the innermost part ( Tsuna et al. 2021 ) .We adopt the model parameters in Ko et al. ( 2023 ) of ( Ṁw , t age , in , M ej ) = (8 × 10 −7 M yr −1 , 20 yr , 0 .35 M ) that reproduce the observed parameters in table 1 , although the discussions below are not sensitive to the assumed values.By solving for the pressure in the thin shell at 2021 A.D. following subsection 3.3 of Ko et al. ( 2023 ) , we find the corresponding gas density to be where m u is the atomic mass, T e,in,ej is the electron temperature of the inner shocked ejecta, and μ ion is the mean molecular weight.We assume an oxygen-dominated composition for simplicity as in the Appendix of Ko et al. ( 2023 ) , and that the ionized region is composed of neutral and singly-ionized oxygen with μ ion = 16( 1 − 0.5 x OII ) , where x OII is the ionization fraction.The corresponding number density of electrons and oxygen ions is The flux of free-free emission from the inner shocked ejecta is thus Here, we have adopted the inner shocked ejecta volume of where M in,ej is the mass of the inner shocked ejecta.
The archive VLASS data constrains the emission to be less than 0.35 mJy, which requires the shocked ejecta to be almost neutral with x OII 2 × 10 −2 , consistent with the above expectation. 1 However, the above scaling also shows that the free-free flux is sensitive to the value of x OII .Here we aim to obtain a conservative estimate for the radio detectability, and assume that free-free emission is negligible as expected for a plausible range of x OII = 10 −4 -10 −3 .On the other hand, the free-free absorption optical depth is given as and is thus expected to be negligible regardless of the exact value of x OII .
To summarize, we consider the free-free emission/ absorption to be negligible in our radio modeling for all the shocked regions.While detailed photoionization calculations would be desired for better estimates of the free-free emission from the inner shocked region, future radio observations in multiple frequency bands may also be able to probe the relative contribution, as the spectral shape of free-free emission is different from that of synchrotron emission.

Discussion
We modeled radio signals from the shocked regions of the SN 1181 remnant, assuming that they are powered by synchrotron emission from electrons accelerated by the shocks ( see also figure 1 ) .Our search using archival radio data resulted in no significant signal, although one targeting the outer region may be a marginal detection.
Nevertheless, the radio signals we estimated may be bright enough to be detectable by targeted observations with current and future radio facilities.In particular, the possibility of detecting radio emissions from the inner shocked region is discussed here, since clarifying the radio emissions from the inner termination shock can constrain the nature of the fast-blowing wind.The inner region is generally bright at ≈10 GHz regardless of the electron spectral index p , and for p = 3 the peak lies around 10 GHz.
In this section, we estimate the future time evolution of the radio emission from the termination shock and discuss the observation strategy.

Future observation prospects
We calculate the time evolution of the inner shocked region, and discuss prospects for future observations.We fix Ṁw and adopt several parameter sets of ( t w , M ej ) which satisfy the constraints reported by Ko et al. ( 2023 ) , including that for r in obtained by high-resolution Chandra observations.Then we calculate the time evolution of r in by solving the equations governing the shock dynamics in Ko et al. ( 2023 ) , and also calculate the time evolution of L ν, in at 10 GHz, adopting the fiducial value of e for p = 2.5 as an example.Here we have also neglected the free-free emission/absorption.
These results are shown in figure 4 .The left-hand panel shows the time evolution of the radius for the parameter set (t age , in , M ej ) = (15 yr , 0 .25 M ) , (20 yr , 0 .35 M ) which correspond to the models we adopted for figure 3 , and (25 yr , 0 .55 M ) .This figure suggests that the source can be spatially resolved with a spatial resolution of at least 1. 5 .
For the light curve at 10 GHz in the right-hand panel, we generally find that the radio emission is gradually decaying at present.We also see that the earlier the wind was launched, the longer the rise time.This is because, for a nearly fixed r in , the earlier the wind launches, the slower the expansion rate, and the slower the rate of change of the absorption optical depth that governs the rise.
For a fiducial value of e = 5 × 10 −4 , the 10 GHz flux in the 2020s will typically be ≈30 μJy, and will very weakly decrease over time in the decades that follow.This is the sensitivity required for current and future radio observations to detect the central radio emission.
For example, the Jansky VLA satisfies the required spatial resolution, so targeted observations may spatially resolve the inner shocked region for a sufficiently long exposure.In particular, next-generation telescopes with even higher sensitivity and angular resolution, such as the Square Kilometer Array 2 and the Next Generation VLA, 3 will achieve these requirements.These projects allow detailed studies of spatially resolved structures of the termination shock region, as well as long-term monitoring to probe the activity of the remnant WD over timescales as long as decades.These high-resolution observations offer a unique probe for the formation mechanisms of the merger ejecta and the remnant wind.Throughout this paper, we assumed spherical symmetry, since the multiwavelength images of the inner and outer shocked regions are overall consistent with a spherical one.Nevertheless, numerical simulations of WD mergers aimed to reproduce SN Iax find the ejecta significantly departing from spherical symmetry ( Kashyap et al. 2018 , figure 4 ) .Magnetohydrodynamical simulations modeling the central wind from a WD merger product also find asymmetry in the wind velocity profile ( Zhong et al. 2024 ) , which can lead to significant angular dependence in the radio power.If free-free emission from the shocked ejecta is important, multi-dimensional effects ( e.g., mixing via Rayleigh-Taylor instabilities) can become important for characterizing the radio emission profile ( D. Tsuna & X. Huang in preparation) .With good spatial resolution ( see figure 4 ) , future observations may probe the asphericity of the innermost region and test these models.

Fig. 1 .
Fig. 1.Schematic picture of SNR 1181.The left-hand panel shows the observed features of this object, and the right-hand panel shows the explanations of each region based on the model proposed by Ko et al. (2023) .We expect synchrotron radiation from the outer SN shock sweeping the interstellar medium ( as in typical SNRs) , and the WD wind termination shock sweeping the inner SN ejecta.

Fig. 2 .
Fig. 2. NVSS 1.4 GHz continuum map centered at SN 1181 ( cross mark) .The contours are ±2 σ and ±3 σ .The beam size is shown in the lower lefthand corner of the map.

Fig. 4 .
Fig. 4. Time evolution of the inner shock radius ( left) and luminosity at 10 GHz ( right) .We fix Ṁwind and v in as values in table 1 , and adopt the fiducial value of e for p = 2.5.The red range in the left-hand panel shows the constraint of the shock radius obtained by Chandra in 2021 A.D.

Table 1 .
Parameters relevant for the WD remnant, the X-ray nebula, and the radio emission modeling in this work.

Table 2 .
Values of p , e adopted in this work. mJy B T e,in,wind 1 as valid in the radio band, where h and k B denotes the Planck constant and the Boltzmann constant, respectively, we obtainF ff ν, in , wind ≈ 6 .8 × 10 −38 T e , in , wind Z 2 ( Ko et al. 2023 )1 .2keV,(18)which is negligible compared to the synchrotron component in the frequency of interest.Here, we adopt Z ion,in = 8, since the abundance of the inner shocked wind is oxygen-dominant( Ko et al. 2023 ).A similar calculation for the outer shocked region yields F ff ν, out ≈ 6 .8 × 10 −38 T e , out Z 2