author: Roland Diehl ; Thomas Siegert ; Wolfgang Hillebrandt ; Martin Krause ; Jochen Greiner ; Keiichi Maeda ; Friedrich K. Röpke ; Stuart A. Sim ; Wei Wang ; Xiaoling Zhang
date: Received 15 September 2014 / Accepted 16 December 2014
title: SN2014J gamma-rays from the $^{56}$Ni decay chain

The detection and measurement of gamma-ray lines from the decay chain of $^{56}$Ni provides unique information about the explosion in supernovae. SN2014J at 3.3 Mpc is a sufficiently nearby supernova of type Ia so that such measurements have been feasible with the gamma-ray spectrometer SPI on ESA's INTEGRAL gamma-ray observatory. The $^{56}$Ni freshly-produced in the supernova is understood to power the optical light curve, as it emits gamma-rays upon its radioactive decay first to $^{56}$Co and then to $^{56}$Fe. Gamma-ray lines from $^{56}$Co decay are expected to become directly visible through the overlying white dwarf material several weeks after the explosion, as they progressively penetrate the overlying material of the supernova envelope, diluted as it expands. The lines are expected to be Doppler-shifted or broadened from the kinematics of the $^{56}$Ni ejecta. We aim to exploit high-resolution gamma-ray spectroscopy with the SPI spectrometer on INTEGRAL towards constraining the $^{56}$Ni distribution and kinematics in this supernova. We use the observations with the SPI spectrometer on INTEGRAL together with an improved instrumental background method. We detect the two main lines from $^{56}$Co decay at 847 and 1238 keV, significantly Doppler-broadened, and at intensities (3.65$\pm$`<!-- -->`{=html}1.21) 10$^{-4}$ and (2.27$\pm$`<!-- -->`{=html}0.69) 10$^{-4}$ ph cm$^{-2}$s$^{-1}$, respectively, at brightness maximum. We measure their rise towards a maximum after about 60–100 days and decline thereafter. The intensity ratio of the two lines is found consistent with expectations from $^{56}$Co decay (0.62$\pm$`<!-- -->`{=html}0.28 at brightness maximum, expected is 0.68). We find that the broad lines seen in the late, gamma-ray transparent phase are not representative for the early gamma-ray emission, and rather notice the emission spectrum to be complex and irregular until the supernova is fully transparent to gamma-rays, with progressive uncovering of the bulk of $^{56}$Ni. We infer that the explosion morphology is not spherically symmetric, both in the distribution of $^{56}$Ni and of the unburnt material which occults the $^{56}$Co emission. Comparing light curves from different plausible models, the resulting $^{56}$Ni mass is determined as 0.49$\pm$`<!-- -->`{=html}0.09 M$_\odot$.

# Introduction

Supernovae of type Ia are generally understood to arise from thermonuclear disruption of a CO white dwarf in a binary system, caused by the rapid nuclear energy release from carbon fusion ignited in the central region of the white dwarf . Different scenarios are discussed on how the explosion may be initiated: Either the accretion of material from the companion star causes the white dwarf to reach the Chandrasekhar-mass stability limit, or an external event on the white dwarf such as a major accretion event or a nuclear explosion on the surface makes the white dwarf interior unstable towards runaway nuclear carbon fusion. Once ignited, nuclear fusion at high densities processes the white dwarf material to iron group nuclei, which are the most stable configuration of nuclear matter, with radioactive $^{56}$Ni being a major product of such explosive supernova nucleosynthesis . As the nuclear flame rushes through the star, an explosion is initiated, and nuclear burning then competes with expansion of the material, resulting in some outer parts of the white dwarf not being burnt towards iron group nuclei, but only to intermediate-mass elements, or even left unburnt as carbon and oxygen mainly . Typically, it is expected that about 0.5 M$_\odot$ of $^{56}$Ni are thus produced and embedded in about 0.5–0.9 M$_\odot$ of other material . As the supernova expands, more and more of the $^{56}$Ni gamma-rays from radioactive decay thus should be able to leave the source region where the decay occurs, and be observable with gamma-ray telescopes .

$^{56}$Ni radioactive decay occurs with a first decay to $^{56}$Co after $\tau\sim$ 8 days, i.e. when the supernova is expected to be still opaque to even gamma-rays at MeV energies, converting this radioactivity energy into emission at lower energy photons . When the second decay stage from $^{56}$Co to $^{56}$Fe at $\tau\sim$ 111 days occurs, producing gamma-rays at 846.77 and 1238.29 keV, the supernova envelope begins to be transparent. After a few months, full transparency to gamma-rays will be reached, and the radioactive decay causes the intensity of characteristic gamma-rays to fade away. Between initial gamma-ray leakage and the stage of full transparency, however, the brightness of $^{56}$Co decay gamma-rays is determined by the amount of available $^{56}$Ni and how it is distributed within the expanding supernova. This is the phase, where the rise and fall of the gamma-ray line intensity provides unique information on the type of explosion and the structure of the supernova.

SN2014J was discovered on January 22, 2014 , and recognised as a type Ia explosion from early spectra . It occurred in the nearby starburst galaxy M82 at $\simeq$ 3.3 Mpc distance , and is seen only through major interstellar material in the host galaxy along the line-of-sight, causing rather large reddening with corresponding difficulties to uncover the proprietary supernova brightness and low-energy spectra at the desired precision. The explosion date, therefore, has been inferred as 14 January, UT 14.75, with 0.21 d uncertainty . The light curve reached its blue-band maximum about 20 days after the explosion . Our detection of $^{56}$Ni gamma-rays early on also show that the initiation of the supernova explosion may have been unusual and from a primary nuclear explosion near the surface of the white dwarf, causing additional uncertainty on the structure of the outer supernova envelope. It is therefore of great interest to study the evolution of the gamma-rays from $^{56}$Co decay, as they trace the release of $^{56}$Ni radioactive energy and its occultation. This should complement data obtained from the UV-optical-IR photosphere further outside, which reveals photospheric material composition through a variety of atomic absorption lines, and can be interpreted using radiation transport calculations .

INTEGRAL started observing SN2014J on 31 January, and kept monitoring the supernova over about 7 Msec until 26 June, when visibility constraints terminated this special opportunity. This covers SN2014J emission during days 16.3 – 164.0 after the explosion. Here we describe our analysis of data from the SPI instrument on INTEGRAL, which is unique with keV energy resolution to reveal details of the shape of the characteristic gamma-ray lines from the $^{56}$Ni decay chain, thus providing kinematic constraints, in addition to tracing the evolution of brightness.

# Data and their Analysis

The INTEGRAL space gamma-ray observatory of ESA carries the gamma-ray spectrometer instrument 'SPI' as one of its two main instruments . Both INTEGRAL main telescopes utilize the *coded mask* technique for imaging gamma-ray sources. In this technique, a mask with occulting tungsten blocks and holes in the field of view of the gamma-ray camera imprints a shadowgram of a celestial source signal in the 19-element detector plane. The SPI spectrometer features a camera consisting of 19 high-resolution Ge detectors, which provides a spectrum of celestial gamma-rays attributed to their source through coded-mask shadowgrams, above a large instrumental background. SPI data consist of energy-binned spectra for each of the 15 Ge detectors of the SPI telescope camera which were operational during our observations of SN2014J.

For our analysis, we used exposures of SN2014J accumulated over a special campaign of the INTEGRAL mission, during orbit numbers 1380 to 1428, with one major gap between 23 April and 27 May. The main exposure window of 2.8 Msec was placed in the rising part of the expected gamma-ray line emission, with an additional 1.4 Msec exposure at late times when occultation of gamma-rays should be a minor effect only, and the total $^{56}$Co decay would be observed. After data selections to suppress contaminations, e.g. from solar flare events, our dataset consists of 4.2 Msec of exposure, spread over 1816 telescope pointings. Exposures are taken in typically 3000 s pointings of the telescope towards this sky region, then moving the telescope pointing by 2.1 degrees to shift the shadowgram of the source in the detector plane. A regular pattern of telescope pointings of a 5 by 5 rectangle around the direction of the supernova itself comprises one such cycle of the dithering exposure pattern. The supernova always was in the telescope field of view, with its $\simeq$`<!-- -->`{=html}30 degree opening angle, and variations of sensitivity due to different aspect angles were $<$`<!-- -->`{=html}10%.

Our analysis method is based on a comparison of measured data to models, performed in a data space consisting of the counts per energy bin measured in each of SPI's detectors for each single exposure of the complete observation. We describe data $d_k$ per energy bin $k$ in general terms as a linear combination of the sky model components $M_{ij}$, to which the instrument response matrix $R_{jk}$ is applied, and the background components $B_{jk}$, with parameters $\theta_i$ for $N_I$ sky and $N_B$ background components: $$d_k = \sum_j R_{jk} \sum_{i = 1}^{N_\mathrm{I}} \theta_i M_{ij} +
\sum_{i = N_\mathrm{I} + 1}^{N_\mathrm{I} + N_\mathrm{B}} \theta_i
B_{jk}\label{eq:model-fit}$$ Generally, in our spectroscopy analysis we fit the intensity scaling factor of a model of the sky intensity distribution plus a set of scaling factors for a model of the instrumental background to data in energy bins covering the spectral range of interest. The response matrix encodes the shadowgram effects, i.e., how the occultation by the coded mask affects visibility of the source direction from each of the Ge detectors of the camera. In our case, we use a single sky component for the SN2014J point source ($N_I=1$), and a single background model ($N_B=1$). The latter is derived by a detailed spectroscopic assessment of longterm background and detector behaviour.

Our treatment of instrumental background follows a new approach, which accounts for the physical nature of instrumental background lines and of detector-specific spectral responses, combining data across a broader range of energy and time periods suitably to build a self-consistent description of spectral detector response and background and their variations. Continuum and line backgrounds are treated separately, individual detector responses and their degradations are determined from a combination of spectral lines and their longterm behaviour. Recognising physical processes which cause characteristic instrumental lines has been part of the validation of the background determination.This method is described in detail in , and was applied for the $^{56}$Ni decay lines from SN2014J successfully; we refer to the Supplementary Information therein for details of the data analysis method and its validation.

Our entire dataset for SN2014J consists of 27240 spectra accumulated from single-detector events in 0.5 keV energy bins. We chose to fit five parameters in this analysis: one intensity scaling for our complete background model, and an intensity amplitude for the SN2014J signal itself for four different epochs. This compromise attempts to account for SN201J gamma-ray variation on the scale of $\simeq$ a month, as expected from models , avoiding to prescribe a model/assumption about the rise and fall of the gamma-ray line intensity, beyond a four-element gamma-ray light curve. We also employ an analysis of eleven different time epochs, when we want to investigate evolutions of spectral features on shorter times of $\simeq$ two weeks. Alternatively, we also fit brightness evolutions from a set of candidate models to our measurements, with normalisations of the respective SN2014J light curve model and of instrumental background.

# Results and Discussion

## Detection of $^{56}$Co lines

We determine the spectrum of the gamma-ray signal from SN2014J in two energy bands around the expected $^{56}$Co decay lines, which have energies at rest of 846.77 and 1238.29 keV, the higher-energy line having 68% of the 847 keV line intensity due to the branching ratio of the nuclear de-excitation. The energy bands chosen are 780 to 920 keV (around the 847 keV line) and 1190 to 1290 keV (around the 1238 keV line). We expect Doppler shift and broadening effects, which would be on the order of 15 keV (21 keV) for 5000 km s$^{-1}$ velocity along the line of sight.

Fig.  shows the spectrum for SN2014J, which was derived from the entire observations set covering days 17 – 164 after the supernova explosion, for the two strongest lines emitted in radioactive decay of $^{56}$Co. These integrated time-averaged spectra were derived from fitting a source at the position of SN2014J in four independent epochs (see below for details), thus allowing for time variability of the flux, as expected.

These spectra show significant emission from SN2014J, overall dominated by broadened lines centered near 847 and 1238 keV, as expected. Gaussian profiles as shown were fitted together with an offset accounting for possible continuum. The flux error bars shown per data point were determined from propagating Poissonian uncertainties through our maximum-likelihood fitting of instrument and background model to the measured dataset; horizontal bars indicate the 10 keV wide energy bins. Overall, the significance of line emission detected from SN2014J in these two energy bands is 9.5 and 3.1 $\sigma$, for the 780 to 920 keV and 1190 to 1290 keV bands, respectively. Generally, we detect characteristic $^{56}$Co gamma-ray line emission in agreement with first-order models of SNIa explosions, and also consistent with results reported by . The lines are modestly broadened and somewhat offset as the $^{56}$Ni produced initially in the explosion is partly occulted behind envelope material, uncovering $^{56}$Ni more and more with time and ejecta dilution. The detection has an overall significance of $\simeq$ 10 $\sigma$, limited from the instrumental background at $\simeq$hundred times higher count rate, and Poissonian statistics.

But one expects that the gamma-ray lines from $^{56}$Co decay are gradually emerging from the supernova, as the overlying material becomes transparent to the gamma-rays with successive supernova expansion . At earlier times, spectrum and intensity of gamma-ray emission from the primary $^{56}$Ni and $^{56}$Co energy source may appear different from the late, gamma-ray transparent phase, where $^{56}$Co decay and ejecta kinematics determines the gamma-ray signal. When we analyze the same dataset of our observations in separate epochs and in finer energy bins for the 780–920 keV band, we obtain a cumulative spectrum shown in Fig. . Here we separated our observations in eleven time epochs (details see below), and 2 keV wide energy bins of analysis, close to the instrumental line width of 2.3 keV. This spectrum reveals that the fitted Gaussian profiles in Fig.  may not capture the actual line emission as it evolves.

## Temporal and spectral variations

### Towards a light curve – four epochs

In a first approach towards time-resolved gamma-ray emission, we subdivide our observations into four epochs of post-explosion days 16.3–41.3 (epoch 1), 41.3–66.3 (2), 66.3-99.1 (3), and 134.8-164.0 (epoch 4). Here the first three epochs fall into the phase of gradually rising emission, while after the observation gap between days 100 and 134, the fourth epoch should capture the $^{56}$Co emission in a rather transparent supernova. We again use 10 keV wide energy bins for the analysis, for a high signal-to-noise ratio per bin, which should be adequate as lines are expected to be kinematically broadened (see Fig. ). Fig.  shows the individual spectra measured at the four different epochs, together with Gaussians fitted to the major features.

At the brightest epoch (3), $^{56}$Co decay gamma-ray line intensities are found as (3.65$\pm$`<!-- -->`{=html}1.21) 10$^{-4}$ and (2.27$\pm$`<!-- -->`{=html}0.69) 10$^{-4}$ ph cm$^{-2}$s$^{-1}$ in the 847 and 1238 keV lines, respectively. From the flux measured at this time, we estimate a $^{56}$Ni mass of 0.50$\pm$`<!-- -->`{=html}0.12 M$_\odot$, assuming supernova transparency (see however below, for an alternative $^{56}$Ni mass estimate).

Now we investigate the time evolution of the detected spectral features in more detail. We aim to trace the intensity variations in time steps that are able to discriminate among different models for the explosion, although we may expect that fainter emission before and after maximum flux will result in larger uncertainty of the flux determined at this time.

```latex
\begin{tabular}{@{} ccccccc @{}} % Column formatting, @{} suppresses leading/trailing space
      \toprule
     \multicolumn{7}  {c}  {847 keV line}   {1238 keV line
     } \\
      \cmidrule(r){2-4} % Partial rule. (r) trims the line a little bit on the right; (l) & (lr) also possible
      \cmidrule(l){5-7} % Partial rule. (r) trims the line a little bit on the right; (l) & (lr) also possible
      time  & flux & center    & width & flux & center    & width  \\
      \midrule
        16.3--41.3 & 2.34 & 827.2 & 14.3       & 0.91 & 1259.1$^1$ & 17.7$^1$ \\
        41.3--66.3 & 2.74 & 851.3 & 11.1       & 1.11 & 1259.1$^1$ & 17.7$^1$ \\
        66.3-99.1 & 3.65 & 851.3 & 20.4       & 2.27 & 1259.1 & 17.7 \\
        134.8-164.0 & 1.90 & 846.6 & 12.9       & 0.38 & 1244.9  & 18.8$^2$ \\
      \bottomrule
   \end{tabular}
```

At late epoch (4), we find the expected pair of clear (at least at 847 keV) and broadened lines near the rest energies at 847 and 1238 keV, with a broadening of $\simeq$ 30 and 44 keV (FWHM), respectively. This broadening is equivalent to a velocity spread of $\simeq$ (4570$\pm$`<!-- -->`{=html}1840) km$^{-1}$, determined from the 847 keV line. From Fig. , it is apparent that the 847 keV line varies in position and width among the four different epochs. Also apparent is that there may be multiple emission features in the 780–920 keV energy band (left column in Fig. ). We characterise the spectra through Gaussians on top of a continuum offset, fitting the main, broadened, line features in the spectra plus a second, narrow, Gaussian in the 780–920 keV range.

Tracing the signature for the 847 keV line from the late epoch (4) towards earlier epochs, we can identify a consistent, broadened line, slightly blue-shifted and broader in epochs (3) and (2), while in epoch (1) the broad feature appears red-shifted towards 827 keV centroid energy, i.e. by $\simeq$ 6,920 km$^{-1}$. We chose to identify this broad feature in the spectra representing the emission of the bulk of $^{56}$Co from the 847 keV line, but we note some arbitrariness here in particular for the first epoch. If we chose to identify the narrow line feature near 860 keV with 847 keV $^{56}$Co, we would obtain the expected blue shift at early time, but the line width, and the origin of the emission centered at 827 keV, would be puzzling. The fitted parameters for the broad lines are shown in Table .

The 1238 keV line appears with a centroid at (1245$\pm$`<!-- -->`{=html}5) keV in the last epoch (4). The generally-weaker intensity levels for this line do not allow an independent determination of line shape parameters. We therefore fix the location and width for the early epochs (1) and (2) with their weak signal on the values fitted in epoch (3), and only fit line intensities here. For the weaker signal in epoch (4), we chose to fix the 1238 keV line width to the Doppler broadening determined for the 847 keV line, and only fit intensity and centroid (we assume that the supernova is transparent to $^{56}$Co gamma-rays here and thus both lines should reflect all $^{56}$Co from the supernova in the same way). We thus allow to some extent for different bulk velocities of the observed $^{56}$Co decay at the four epochs, trying to derive as much as we can from the data themselves. We find the 1238 keV line somewhat blue-shifted earlier in (brighter) epochs (2) to (3) at 1259 keV; statistical precision is inadequate to determine its variation at different epochs, in particular the red-shift as indicated in epoch (1) for the 847 keV line cannot be constrained nor confirmed independently.

|    time     |              bulk               |             spread             |
|:-----------:|:-------------------------------:|:------------------------------:|
|  16.3–41.3  | -6920$\pm$`<!-- -->`{=html}1480 | 5060$\pm$`<!-- -->`{=html}1330 |
|  41.3–66.3  | 1600$\pm$`<!-- -->`{=html}1720  | 3940$\pm$`<!-- -->`{=html}1260 |
|  66.3-99.1  | 1600$\pm$`<!-- -->`{=html}1600  | 7250$\pm$`<!-- -->`{=html}1560 |
| 134.8-164.0 |  -80$\pm$`<!-- -->`{=html}1870  | 4570$\pm$`<!-- -->`{=html}1840 |

Velocity values in km s$^{-1}$ as derived from the 847 keV $^{56}$Co line fits in four epochs shown in Fig. , and listed in Table . The 1238 keV line derived centroid and spread are only determined in brightest epoch 3, and are 5050$\pm$`<!-- -->`{=html}1240 and 4290$\pm$`<!-- -->`{=html}1370 km s$^{-1}$. Fixing the 1238 keV line's width to the value determined from the 847 keV line, its centroid also for the weak epoch (4) can be found, at 1610$\pm$`<!-- -->`{=html}960 km s$^{-1}$.

Table  summarises the velocity constraints from this four-epoch analysis. Overall, both lines are consistent in their centroids and broadenings, within uncertainties, although differences are remarkable. The bulk Doppler shifts differ somewhat between the two lines, even at late times (near transparency) we find -80$\pm$`<!-- -->`{=html}1870 km s$^{-1}$ and 1610$\pm$`<!-- -->`{=html}960 km s$^{-1}$ for the 847 and 1238 keV lines, respectively. This may tell us that the $^{56}$Co visible in each of the lines reflects a different subset of the total, with different spatial sampling and thus kinematics. Note that at early times, it is not clear *a priori* which fraction of the $^{56}$Co of the 3-dimensional exploding supernova is visible, and what the magnitude of occultation is, as both $^{56}$Ni and overlying ejecta morphologies are unknown, while towards late times, the true kinematic signature of all of the $^{56}$Ni should be reflected in both $^{56}$Co lines in a consistent way. The transparency of the supernova envelope is expected to vary with energy, transparency being higher at higher gamma-ray energy.

The intensity variation throughout these four epochs of our observations for both lines produce gamma-ray light curves shown in Fig. . Evidently, both lines consistently rise towards a maximum near 60–100 days after explosion, falling off later. The intensity ratios between both lines (see Table ) generally agree (within uncertainties) with the nominal branching ratio of $^{56}$Co decay, in particular for the brightest epoch (3) we find the 1238 keV line being at 62% ($\pm$`<!-- -->`{=html}28%) of the 847 keV line intensity, which compares to a laboratory value of 68%.

For assessment the detection of supernova emission from $^{56}$Co decay, we also check upon our statistical uncertainties. For an independent estimate of uncertainties, we histogram the resulting SN2014J flux values, and compare their distribution near zero flux with expectations from Poissonian statistical uncertainty (these are the error bars shown in Fig. ). Fig.  shows the flux value histogram for all observations in the 780–920 keV band, analysed in the above eleven time bins, and 2 keV energy bins. The distribution of values towards negative fluxes from zero is due to statistical fluctuations only, while a celestial source contributes to the distribution at positive flux values. We thus obtain flux uncertainties from the negative part of the distribution, and total source significance from its positive part. This yields a KS-test p-value of 1.4 10$^{-29}$ for the positive flux values following the same statistics-only distribution, which is equivalent to a probability of 11.3 $\sigma$ that the spectrum contains nonzero signal from SN2014J. Measuring the width of the distribution of fluxes below zero, we find 1.0 10$^{-5}$ ph cm$^{-2}$s$^{-1}$, which compares to a Poissonian error of $\pm$`<!-- -->`{=html}0.54 10$^{-5}$, i.e. 1.08 10$^{-5}$  ph cm$^{-2}$s$^{-1}$ (the error in data points shown in the spectra of Fig. ). We conclude that our Poissonian statistical error estimates are approximately correct. Using this assessed statistical uncertainty, the narrow-line signals from SN2014J shown in Fig.  in the 780–920 keV energy range correspond to a statistical significance of $\simeq$ 2.8–4 $\sigma$.

### Exploring shorter time variations and high spectral resolution

The spectrum in the 847 keV line band of SN2014J (Fig. ) shows, beyond the broadened line discussed above, that apparently additional but narrow lines seem to be present at significance levels exceeding 3 $\sigma$, near 865 keV in epoch (1), 805 keV in epoch(2), and 815 keV in epoch (3). No such features are apparent in the high-energy band 1190–1290 keV.

The large instrumental background may be a concern, creating artificial transient features. We verified adequacy and consistency of our background modeling in these epochs. The orbit-by-orbit spectral fits do not show indications for deviations of fit quality in those energy regions and times. None of the background components show correlations with the sky signal.

```latex
\begin{table*}[htbp]\centering \small
   \begin{tabular}{@{} cccr @{}} % Column formatting, @{} suppresses leading/trailing space
      \toprule
Epoch  &  $\chi^{2}_{best}$ (dof)  &  $\chi^{2}_{alltimes}$ ($\sigma$)  &  $\chi^{2}_{from4thepoch}$ ($\sigma$) \\
    \midrule
1      &  12.62 (6)                &  42.05 (4.28)                      &  42.93 (4.36)                         \\
2      &  17.77 (6)                &  27.48 (1.73)                      &  29.35 (2.04)                         \\
3      &   5.77 (6)                &  15.37 (1.71)                      &  15.79 (1.78)                         \\
4      &  24.78 (9)                &  25.54 (0.41)                      &  - (-)                                \\
     \bottomrule
   \end{tabular}
   \quad\quad
   \begin{tabular}{@{} cccr @{}} % Column formatting, @{} suppresses leading/trailing space
      \toprule
Epoch  &  $\chi^{2}_{best}$ (dof)  &  $\chi^{2}_{alltimes}$ ($\sigma$)  &  $\chi^{2}_{from4thepoch}$ ($\sigma$) \\
    \midrule
1      &  48.36 (62)               &  66.85 (3.04)                      &  66.82 (3.03)                         \\
2      &  66.90 (62)               &  77.73 (1.92)                      &  76.58 (1.72)                         \\
3      &  43.22 (62)               &  60.95 (2.94)                      &  58.81 (2.65)                         \\
4      &  50.66 (65)               &  51.18 (0.29)                      &  - (-)                                \\
     \bottomrule
   \end{tabular}
  \caption{Line shape variations of the 847 keV line. For each of the four observation epochs, we list the $\chi^2$-values of the best fit (column 2; Gaussian with adapted centroid and width), then the intensity-fitted shape of the line as determined in the time-integrated results (Fig.~\ref{fig:SPI-spectrum-847-1}) (column 3), and finally (column 4) the intensity-fitted line shape of the last epoch. The degrees of freedom of the fits are given in column 2 in brackets. The discrepancy of the thus-adopted line shape in each of the epochs is given in $sigma$ units in brackets in columns 3,4. The {\it lefthand} table gives the results for the 10 keV bin data (center column in Fig.~\ref{fig:SPI-lightcurve-847-1238}), the {\it righthand} table gives the results for the 2 keV bin data (left column in Fig.~\ref{fig:SPI-lightcurve-847-1238})}
   \label{table:lineshapefits}
\end{table*}
```

Having confirmed the statistical uncertainties (see above), we now investigate how different resolution in time and in spectral bin width affect our analysis results. Narrow-line features with fluxes at the level of 2 10$^{-4}$ ph cm$^{-2}$s$^{-1}$ are reminiscent of $^{56}$Ni decay, because $^{56}$Ni volume emissivity is much higher due to the shorter decay time than for $^{56}$Co lines, and thus smaller active volume elements with specific kinematic and Doppler shift might contribute transient lines with a typical 1–2 week time scale. Also, if the $^{56}$Ni synthesised in the explosion emerges in clumps, each clump may be expected to have different kinematics and somewhat different turbulence, cooling, and expansion. We analyse the 780–920 keV band therefore in eleven time epochs independently, and in finer energy bins of 2 keV resolving lines as narrow as the instrumental line width. We find that in finer time bins the broad $^{56}$Co line features are difficult to recognise, except towards the late epoch in days 134–164. We also note that those three late-epoch spectra are otherwise featureless and rather flat, which is reassuring. However, all earlier spectra show, more or less, spectral features which are less broad, and which altogether combine to form the broad features seen at coarser energy binning. It appears that between 830 and 860 keV, $\simeq$ 10 keV wide features emerge at times before day 99, but vary in intensity and peak location detail, with a consistent trend towards lower energies. Moreover, around 810 keV, another line feature emerges before days 88, which started to emerge only after days $\simeq$ 60. The feature near 865 keV apparent in the first sample is gone after days 25, while near 820 keV, emission is indicated near days 30 and 55, not seen in adjacent epochs. We conclude that, while being consistent with appearance of the 847 keV line emission, this variability is surprising. When we assume an intrinsic line shape of $^{56}$Co as seen at late times, or integrated and as seen in broad energy bins (Fig. ), we find that the line shapes in the different epochs shown in Fig.  are significantly different in each case (see Table ). This may be understood as different sight lines to embedded $^{56}$Co determining the observed lines at each of our epochs, which in turn implies that the envelope of the supernova includes 3D structure which evolved during the time of our gamma-ray observations.

$^{56}$Ni decay has two strong lines, one at 812 keV, and another one at 158 keV line. We find no emission near 158 keV coincident with the narrow line features in the 800-900 keV band, and therefore we can exclude an origin in late $^{56}$Ni emission (e.g. from fully-ionised $^{56}$Ni regions).

If we take our results as showing three distinct emission components, we may employ a simple model to treat them individually in terms of $^{56}$Ni content, kinematics, and occultation, embedded in an expanding supernova. We roughly identify components at 810$\pm$`<!-- -->`{=html}5 keV, 845$\pm$`<!-- -->`{=html}10 keV, and 865$\pm$`<!-- -->`{=html}7.5 keV. We find that flux time histories of each component are consistent with early occulted and later revealed decaying $^{56}$Co line emission, and we determine optical depths at day 1 of 2000, 600, and 100 for the 810, 845, and 865 keV centered components, respectively. (High opacity for the volume emitting 810 keV $^{56}$Co arises due to early contributions from $^{56}$Ni in the 812 keV line.) The kinematics suggested for such three identified major $^{56}$Co clumps are roughly blue-shifted 6400 km s$^{-1}$, ($\simeq$ 15% of the $^{56}$Ni), red-shifted 13,000 km s$^{-1}$, ($\simeq$ 30% of $^{56}$Ni), and a main, spherically symmetric contribution ($\simeq$ 55% of $^{56}$Ni). While we do not propose this to be the real morphology of SN2014J, this analysis supports our interpretation of aspherical $^{56}$Ni distribution and differences in their respective occulting outer supernova material.

In conclusion, we speculate that occultation and its evolution as the supernova expands may be responsible for the apparent spectral signatures. Clumps or co-moving volume elements carrying $^{56}$Co may lie along a less-occulted line of sight at specific times, thus contributing emission in a particular bulk Doppler shifted energy regime. As the supernova expands, different volume elements may thus contribute at different times, as long as occultation is significant. This may reflect asymmetry in SN2014J, which had been discussed earlier to characterise a subset of SNIa where multiple or staged explosions occur .

# Model comparison and $^{56}$Ni mass

```latex
\begin{table*}[htbp]\centering \small
   \begin{tabular}{@{} lccccccccr @{}} % Column formatting, @{} suppresses leading/trailing space
      \toprule
  \multicolumn{3}  {c} \quad line fit    &  & & 100 keV band  & & & 20 keV band \\
    \midrule
 rank     &   model         & $\chi^2$ & M($^{56}$Ni) &   model & $\chi^2$ & M($^{56}$Ni) &   model & $\chi^2$ & M($^{56}$Ni) \\
      \midrule
1.          &    hed8      &  25.07   & 0.46$\pm$0.06   & hed8  &   7.81    & 0.52$\pm$0.12 & W7 & 16.98 &  0.42$\pm$0.07 \\
2.          &   w7dt               &  25.37  & 0.50$\pm$0.06    & w7dt    & 7.89        & 0.57$\pm$0.13    & det2    & 17.07    & 0.38$\pm$0.06 \\
3.          &   W7A           &   25.57  & 0.51$\pm$0.06    & W7A    & 7.95    & 0.59$\pm$0.13    & m36    & 17.15    & 0.43$\pm$0.07 \\
4.          &   hecd             & 25.61  & 0.52$\pm$0.06    & hecd    & 8.11    & 0.58$\pm$0.13    & dd202c    & 17.26    & 0.42$\pm$0.07 \\
5.          &   hed6             & 26.22  & 0.50$\pm$0.06    & hed6    & 8.54    & 0.56$\pm$0.13    & W7E    & 17.29    & 0.43$\pm$0.07 \\
6.          &   det2             & 27.84  & 0.50$\pm$0.06    & det2    & 9.35    & 0.55$\pm$0.13    & dd4    & 17.33    & 0.45$\pm$0.07 \\
7.          &   w7dn    & 28.27  & 0.54$\pm$0.07        & w7dn    & 9.57    & 0.59$\pm$0.14    & w7dn    & 17.33    & 0.40$\pm$0.06 \\
8.          &   W7E         & 29.48  & 0.56$\pm$0.07        & W7E    & 10.11    & 0.61$\pm$0.15    & det2e2    & 17.67    & 0.44$\pm$0.07 \\
9.          &   m36         & 30.30  & 0.56$\pm$0.07        & m36    & 10.44    & 0.61$\pm$0.15    & hed6    & 17.73    & 0.37$\pm$0.06 \\
10.        &    W7          & 30.52  & 0.54$\pm$0.07        & W7    & 10.50 & 0.59$\pm$0.14    & hecd    & 18.13    & 0.38$\pm$0.06 \\
11.        &    dd202c    & 30.70  & 0.55$\pm$0.07    & dd202c    & 10.51    & 0.61$\pm$0.15    & hed8c    & 18.17    & 0.33$\pm$0.05 \\
12.        &    dd4         & 32.72  & 0.57$\pm$0.08    & dd4    & 11.35    & 0.62$\pm$0.15    & w7dt    & 18.33    & 0.36$\pm$0.06 \\
13.        &    det2e2     & 34.18  & 0.56$\pm$0.07    & det2e2    & 11.80       & 0.61$\pm$0.15    & W7A    & 18.40    & 0.38$\pm$0.06 \\
14.        &    pdd54      & 38.84  & 0.56$\pm$0.08    & pdd54    & 13.52    & 0.61$\pm$0.16    & pdd14    & 18.74    & 0.46$\pm$0.07 \\
15.        &    det2e6     & 50.75  & 0.72$\pm$0.11    & det2e6    & 17.80      & 0.78$\pm$0.25    & det2e6    & 22.94    & 0.66$\pm$0.11 \\
             &  constant      & 35.20  & 0.30$\pm$0.06    & constant    & 8.48    & 0.41$\pm$0.12    & constant    & 26.86    & 0.24$\pm$0.06 \\
             &  bgd only     & 90.92  & -    & bgd only &   27.87    & -    & bgd only    & 57.33    & - \\
     \bottomrule
   \end{tabular}
   \caption{$^{56}$Ni mass fit results for different explosion models, as described in detail by \citet{The:2014}. { We fit the light curves resulting from those models to our gamma-ray data. The fit} quality can be seen from the $\chi^2$ values { (for 9 degrees of freedom, except for the bgd-only case which has 10 d.o.f.)}. Fitted $^{56}$Ni masses are given in M\ensuremath{_\odot} with uncertainties, for the intensity variation time profile of the respective model. We compare results for three analysis approaches, fitting a single broad line \emph{(left)}, for taking the non-zero flux in a broad band \emph{(middle)} (800--900~keV), and for taking nonzero flux in a band where the $^{56}$Co emission could be expected to appear \emph{(right)} (840--860~keV).}
   \label{table:model-parameters}
\end{table*}
```

We may compare the time evolution of our measured fluxes of $^{56}$Co line emission to expectations from different models. For this, we use standard models derived for generic plausible assumptions over the past decades, as presented and discussed by, e.g., . Recenlty, have recompiled a representative set of models specifically for applications to SN2014J, and we use the set discussed there in our aim to capture key properties that may describe SN2014J best.

We return to the issue of how to properly assign $^{56}$Co line emission to earlier epochs, given that the observed emission is not represented by a broad (Gaussian) shape, but rather by evolving appearance of emission across the energy range finally described through such a broadened line (see Fig. ); we also remind that Comptonisation will incur some apparent redshift of $^{56}$Co line emission, while the kinematics of $^{56}$Ni ejected in all directions will imprint its own Doppler shifts, and will be revealed as the overlying supernova material becomes increasingly transparent. This implies that the spectral shape of $^{56}$Co emission is uncertain and depends on the 3-dimensional distributions of $^{56}$Ni and of occulting, overlying material, which is not accounted for in any of these (1-dimensional) models. Therefore, we concentrate on the brighter, more-clearly measured $^{56}$Co line at 847 keV, and then we compare three different approaches to determine $^{56}$Co decay emission originating from the 847 keV line in our different epochs:

Approach (a) uses the line parameters as fitted in late epochs, and fits earlier epochs constraining line width and centroid to within the 2$\sigma$ uncertainty limits of that late-epoch fit. Thus we stabilise the way we infer line emission from this energy region, using the late-epoch emission as the constraint, within earlier revealed emission spots must fall. Model-fitted $^{56}$Ni masses are in the range 0.46–0.52 ($\pm$`<!-- -->`{=html}0.06) M$_\odot$ for most-plausible models with some $^{56}$Ni near the supernova surface (see ). Another approach (b) assumes, as another extreme, that *all* emission above the zero level in the 800–900 keV energy band may be attributed to $^{56}$Co emission originating from the 847 keV line emission. Clearly, this over-estimates early epoch intensties, which, additionally, may carry contributions from the more-intense $^{56}$Ni decay. Correspondingly, model-fitted $^{56}$Ni masses are higher and in the range 0.52–0.59 ($\pm$`<!-- -->`{=html}0.13) M$_\odot$. Approach (c) is an intermediate case, evaluating the energy range of the broad 847 keV line, yet not prescribing a spectral shape constraint. Averaging over the three best-fitting models in our three approaches, we derive a $^{56}$Ni mass of 0.49$\pm$`<!-- -->`{=html}0.09 M$_\odot$.

We adopt the model-predicted $^{56}$Ni masses and $^{56}$Co line fluxes versus time for each of the models, and fit those *light curves* to our data as defined in approaches (a) to (c), scaling their intensity. This scaling factor then determines the $^{56}$Ni mass for SN2014J from each model, as fitted to our observations. The results shown in Table  provide SN2014J $^{56}$Ni masses in the range $\simeq$ 0.4 to 0.8 M$_\odot$. The best-fitting models with their He on the outside also probably give the best $^{56}$Ni mass estimates around (0.5$\pm$`<!-- -->`{=html}0.1) M$_\odot$. The time profiles of even the best-fitting models do not match in detail the fluxes determined for a single broad line, but improve significantly when the spectral constraints are relaxed, thus indirectly confirming our above claim for irregular appearance of $^{56}$Co lines. In Table  we also see that our observations, while providing a clear detection of the $^{56}$Co decay gamma-ray lines, are still inadequate to clearly discriminate among models; most models are in agreement with our observed intensity evolution of the gamma-ray lines, if $^{56}$Ni is a free parameter. Better constraints could be expected, if models with different $^{56}$Ni and ejecta morphologies and corresponding radiation transport would be available, and tested on the broader energy range including all $^{56}$Co lines and continuum emission as it may evolve.

# Conclusions

We analyse the full set of INTEGRAL/SPI observations available for SN2014J, and concentrate on the energy ranges around the brightest gamma-ray line emission expected from $^{56}$Co decay, around the two bright lines at 847 and 1238 keV energy. We employ an analysis that uses the coded-mask shadowgrams as the source aspect is varied across those observations, and searches for emission consistent with these variations, above a large instrumental background. The instrumental background is derived by a detailed spectral analysis and modeling, which accounts for continuum and many instrumental lines arising from cosmic-ray interactions with instrument and spacecraft, in a broader energy range around the lines of interest. In this, we first identify long-time spectral features, then separately fit their appearance in each individual Ge detector with its response characteristics, and finally determine the amplitude variations of the background spectral model on short time scales of our individual pointings. The model background is then extracted for the energy bins where we aim to determine a flux determination from SN2014J. This method overcomes limitations of statistical precision of background measurements in the fine energy bins where we wish to analyse the source, and uses data as measured by the SPI detectors themselves, and the physics extracted from longterm monitoring of background and detector responses.

Our measurements of $^{56}$Co decaying in SN2014J find the two strongest gamma-ray lines at 847 and 1238 keV, clearly distinguished from instrumental background and attributed to the supernova. We analyse separately different epochs of our observations, and derive an intensity time profile in the energy regime around these two lines. The brightness evolution is consistent between both lines, and is overall consistent with expectations from different plausible explosion models. We find that models preferred by the observed gamma-ray light curve also include those with $^{56}$Ni on the outside. We obtain $^{56}$Ni mass values of 0.5–0.6 M$_\odot$ for best-fitting models, here determined from gamma-ray light curve data of the $^{56}$Co decay for the first time. We note, however, that the appearance of the line emission is irregular and not recognised as a smooth emergence of the Doppler broadened emission lines as they are found in the late, gamma-ray transparent observations of SN2014J. We find transient emission features of different widths, but cannot assign them clearly to origins of $^{56}$Ni and $^{56}$Co at specific kinematic properties. We infer that the 3-dimensional structure of the inner supernova, and possibly also of the material overlying the $^{56}$Ni, are not as regular as 1-dimensional models of gamma-ray emergence currently implement.

Our data cannot discriminate among models, although Chandrasekhar-mass models are slightly less favoured. The observed line shapes clearly suggest that the explosion was not spherically symmetric. These results demonstrate the power of gamma-ray line observations from type Ia supernovae, and seed high expectations should a closer event occur during the INTEGRAL lifetime.


abstract: SuperCDMS is currently operating a 10-kg array of cryogenic germanium detectors in the Soudan underground laboratory to search for weakly interacting massive particles, a leading dark matter candidate. These detectors, known as iZIPs, measure ionization and athermal phonons from particle interactions with sensors on both sides of a Ge crystal. The ionization signal can be used to efficiently tag events at high radius and near the top and bottoms surfaces, where diminished charge collection can cause events to mimic WIMP-induced nuclear recoils. Using calibration data taken with a $^{206}$Pb source underground at Soudan, we demonstrate rejection of surface events of $(4.5 \pm 0.9) \times 10^{-4}$ with 46% acceptance of nuclear recoils using the phonon signal only. We also show with $^{133}$Ba calibration data underground that the phonon channels can efficiently identify events near the sidewall. This phonon-based approach can also be extended to lower energies than the ionization-based position reconstruction.
author: Adam Anderson, for the SuperCDMS Collaboration
date: 07.15.2013
title: Phonon-Based Position Determination in SuperCDMS iZIP Detectors

# Introduction

The SuperCDMS experiment is a dark matter search consisting of 15 ultra-pure cylindrical germanium crystals known as interleaved z-sensitive ionization and phonon detectors (iZIPs) each with a mass 0.62 kg and operated at 50 mK at the Soudan Underground Laboratory in Minnesota. Transition-edge sensors (TESs) and aluminum phonon absorbers called QETs are patterned on the top and bottom surfaces of the crystals and grouped in parallel arrays to form 8 phonon channels. These phonon channels measure non-equilibrium phonons produced by particle interactions in the detector. Electron-hole pairs are also produced by particle interactions, and this ionization energy is measured by electrodes that are interleaved with the phonon sensors. When WIMP scatters off a nucleus, it produces a nuclear recoil as opposed to the electronic recoils produced by the majority of background radiation. At a fixed recoil energy, the fraction of energy deposited in ionization, called the ionization yield, is approximately 0.3 for nuclear recoils (NR), while it is 1 for electronic recoils (ER). The simultaneous measurement of ionization and phonons therefore provides powerful discrimination between well-reconstructed electronic and nuclear recoils.

Background events near the top and bottom surfaces and around the outer sidewall of an iZIP can exhibit ionization yield below 1, increasing the chance of leaking into the NR signal band. Betas at 100 keV and gammas at 10 keV stop at about 50 $\mu$m from the surface, where ionization collection can be lower due to back-diffusion of energetic charges into the wrong polarity electrode . Radiogenic nuclei such as $^{206}$Pb from decay of $^{210}$Pb implanted on nearby passive material produce nuclear recoils with ionization yield about 0.2 and a penetration length of microns. Finally, the electric field is weak and fringing near the outer sidewall of the detector, which may cause charge trapping that reduces the yield of all events, including high-energy gammas.

The ionization electrode structure of the iZIP efficiently tags potential background events at both high radius and near the surface. An outer guard electrode forms a ring around the outer edge of each detector face, which measures a signal when the event is at high radius. The phonon and charge sensors are biased at 0V and +2V (-2V), respectively, on side 1 (side 2). This arrangement produces an overall drift field that causes electrons and holes to drift to opposite sides of the detector when an event occurs in the bulk. Electrons and holes from events near the surface, on the other hand, both drift to the same side of the crystal. The difference in charge collection between the two sides is therefore a powerful discriminator of surface events .

Since the phonon channels are sensitive to non-equilibrium phonons, some of the position information of the event is contained in the pulse shapes measured by the phonon sensors and the partitioning of energy between the different channels. Phonon-based identification of surface events and sidewall events is advantageous because it provides some redundancy for the charge-based methods. It also has potential to function at lower energies than the charge, since the phonon resolution is better than the charge resolution in iZIP detectors, and charge-based position estimation is worse, by definition, in regions where charge collection is poor. Extending discrimination to lower energies is especially important in searches for light WIMPs, a topic of recent experimental interest .

# Phonon Physics in the iZIP

When an event occurs in an iZIP, an initial population of high-frequency ($\nu > 1$ THz) prompt phonons are produced. These phonons interact primarily by elastic isotope scattering and inelastic anharmonic decay. The rate for isotope scattering in Ge is given by $(3.67 \times 10^{-41} \textrm{s}^3)\nu^4$, while the rate for anharmonic decay is given by $(6.43 \times 10^{-55} \textrm{s}^4)\nu^5$. At THz frequencies, anharmonic decay causes rapid downconversion of phonons to lower energies. As they down convert, the isotope scattering and anharmonic decay rates decrease quickly and the mean free path of the phonons increases, resulting in a quasidiffusive propagation of phonons away from the initial interaction location. To see how the position information can be probed, consider that a 0.5 THz phonon has a lifetime of about 40 $\mu$s–much longer than the $1.6\mu$s sampling time of the readout electronics–and the radius of the quasidiffusively expanding phonon cloud is expected to be about 0.2 mm These prompt phonons are preferentially absorbed by sensors near the interactions site since the high rate of scattering can produce many interactions with nearby instrumented surfaces. This results in significant position information in the early time component of the pulse, with nearby sensors measuring higher amplitudes and faster risetimes. As time continues, the phonons continue to downconvert until their mean free path exceeds the size of the crystal in the ballistic regime, at which point all position information is lost.

Position information is also carried by Luke-Neganov phonons produced when charges drift across the crystal. When charge carriers in Ge reach their constant drift velocity, additional potential energy released as they drift through the electric field is converted to Luke-Neganov phonons. Since the electric fields are strongest near the top and bottom surfaces of the detector, the Luke phonon power emission is greatest there. Electrons and holes emit Luke phonons at the frequency scales of about 0.3 THz and 0.7 THz, which results in quasidiffusive propagation in a localized cloud near the surface, and strong absorption in a single phonon channel.

# Testing Discrimination

To test the ability of the phonon signal to identify different types of events, data was taken using two iZIP detectors running in the Soudan Underground Laboratory with $^{210}$Pb sources and a $^{133}$Ba source. The former produce low-energy gammas, betas, and $^{206}$Pb nuclei, as shown in Figure , which stop near the surface and are generically classified as "surface events" that can exhibit reduced ionization yield and contaminate the nuclear recoil signal region. They consist of silicon wafers that were implanted with $^{210}$Pb by exposure to $^{222}$Rn . The wafers were positioned above one side of two different detectors in Soudan, and 37.6 live days of unblinded low-background data were analyzed as shown in Figure . The $^{133}$Ba source has strong gamma lines at 302.8 keV, 356 keV, and 383.9 keV, which can produce Compton scatters at high-radius where charge trapping frequently occurs.

For any individual phonon channel or combination of channels, a pulse height estimator is constructed using an optimal filter with a single pulse template. One simple method for quantifying the position-dependence of each event is to compare the relative partitioning of energy between different phonon channels. The resulting phonon energy partitioning will be correlated with the spatial position of the events, so that events near the surface and outer sidewall of the detector can be identified in phonons even if their charge propagation and reconstruction is poor.

To tag events near the surface from the $^{210}$Pb source, we can compare the phonon energy measured on side 1 against side 2 of the detector. Traces from all channels on each side are summed and an optimal filter is used to compute the pulse amplitude. The right panel of Figure shows low-yield data from a detector with a source above side 1. These events are overwhelmingly surface events, and they primarily lie in a band below the line of equal energy partitioning. In contrast, nuclear recoils are centered in a broad band with equal energy collected on each side. A cut is constructed in bins of side 1 phonon energy to optimize the signal significance of nuclear recoils at each energy. The resulting cut values are then fit to the functional form shown by the purple lines in Figure ; events inside the purple bands are accepted, while events outside are rejected. A small spur of events is also visible above the upper cut boundary at low energy. These events are due to a period of abnormally high electronics noise on one channel, rather than true surface events, and they are ignored in the cut optimization.

To tag events at high radius in the $^{133}$Ba calibration data, we sum the pulses from the two outer channels and then compare to the total energy of the event to obtain the fraction of energy in the outer channel. Since the calibration gammas consist of a mixture of events near the sidewall and events in the interior of the detector, one must use the charge information to tag high radius events. Figure shows the fraction of energy measured in the outer phonon channel for non-surface nuclear recoils in $^{252}$Cf calibration data, compared with low-yield ($<0.7$) electronic recoils in $^{133}$Ba that are tagged as being at high radius by the outer charge channel.

# Discrimination Performance

The optimized cut shown in Figure is able to achieve a surface event leakage of $(4.5 \pm 0.9) \times 10^{-4}$ over an energy range from 8.0 to 115 keV recoil energy with a nuclear recoil passage fraction of 46% when measured on $^{252}$Cf nuclear recoils. While this is somewhat worse than the surface event rejection achieved using the charge information, the phonon surface event discrimination is an independent tag which can complement charge-based methods. It is also extensible to lower energies than the charge-based discrimination, which is a topic of active study within SuperCDMS.

The radial cut is additionally able to tag high-radius charge-suppressed gammas similarly to the charge signal. Cutting events with greater than 24% of their energy in the outer channels preserves 85% of the inner nuclear recoils, while rejecting 99.9% of the low-yield outer events that could potentially leak into the signal region. This shows that even for low-yield events with reduced charge collection, the outer phonon channels can still efficiently distinguish gammas near the sidewall from inner nuclear recoils.

# Conclusion and Perspectives

Phonon propagation in SuperCDMS iZIP detectors operating at Soudan provides strong position sensitivity in regions of the detector which are prone to reduced charge collection that produces leakage into the WIMP-induced NR band. Although the iZIP has excellent charge-based defenses against surface and sidewall events, the simple partition estimators of the phonon energy provides even more discrimination against these events. Phonon-only methods are also extensible to lower energies than charge-based methods and will be used in future low-energy threshold analyses of SuperCDMS data. Time-domain pulse-fitting algorithms for complex pulse shape discrimination are also under active development, which may further improve performance over partition-based position estimators.


abstract: In recent times, machine learning (ML) and artificial intelligence (AI) based systems have evolved and scaled across different industries such as finance, retail, insurance, energy utilities, etc. Among other things, they have been used to predict patterns of customer behavior, to generate pricing models, and to predict the return on investments. But the successes in deploying machine learning models at scale in those industries have not translated into the healthcare setting. There are multiple reasons why integrating ML models into healthcare has not been widely successful, but from a technical perspective, general-purpose commercial machine learning platforms are not a good fit for healthcare due to complexities in handling data quality issues, mandates to demonstrate clinical relevance, and a lack of ability to monitor performance in a highly regulated environment with stringent security and privacy needs. In this paper, we describe *Isthmus*, a turnkey, cloud-based platform which addresses the challenges above and reduces time to market for operationalizing ML/AI in healthcare. Towards the end, we describe three case studies which shed light on *Isthmus*' capabilities. These include (1) supporting an end-to-end lifecycle of a model which predicts trauma survivability at hospital trauma centers, (2) bringing in and harmonizing data from disparate sources to create a community data platform for inferring population as well as patient level insights for Social Determinants of Health (SDoH), and (3) ingesting live-streaming data from various IoT sensors to build models, which can leverage real-time and longitudinal information to make advanced time-sensitive predictions.
author: Akshay Arora; Arun Nethi; Priyanka Kharat; Vency Verghese; Grant Jenkins; Steve Miff; Vikas Chowdhry; Xiao Wang
title: Isthmus: Secure, Scalable, Real-time and Robust Machine Learning Platform for Healthcare

# Introduction

Scalable ML/AI driven solutions that are tightly integrated in the clinical decision support system have the potential to improve outcomes, reduce cost, and increase efficiency in healthcare. While healthcare has made steady progress through the years in using data to make evidence-based clinical decisions through randomized clinical trials , \], it has lagged to adopt ML/AI driven solutions at scale. Developments over the last decade have made use of Electronic Health Records (EHRs) widespread in healthcare organizations, which has resulted in an abundance of structured and unstructured clinical and non-clinical data ,but the industry has struggled to generate value from this data. The technical reasons for this slow uptake include uneven data quality, privacy, and confidentiality issues in building workflows, and a lack of interoperability between different systems . To counter the concerns above and deploy ML predictive models in the clinical workflows, we developed a turnkey ML platform called *Isthmus* to seamlessly develop, test, deploy, evaluate, and retrain predictive models.

Developed under clinicians' guidance and oversight from Parkland Health & Hospital System's stakeholders, *Isthmus* balances the domain-rich customization needed for building clinical, predictive models in healthcare with the need of having a turnkey, secure, HIPAA compliant solution. Consequently, the platform reduces time to market to integrate and update AI/ML models to improve the quality and safety of a patient's hospital experience and drive costs down by impacting both clinical and operational outcomes.

# Methodologies

## *Isthmus* Framework

Our vision was to build an end-to-end machine learning framework which would make predictive model development, deployment, evaluation, and retraining seamless by reducing the time to market for integrating clinical predictive insights in clinical workflows to make them actionable. To achieve the goals above we built a cloud hosted, machine learning platform called "*Isthmus*."

*Isthmus*' end-to-end workflow will reduce implementation time when there are multiple models developed and deployed on it,by reusing the environment toolsets, technologies, and data/feature engineering pipelines to either train or score using an already deployed pre-trained model.

Part of the motivation to build such a flexible but scalable and configurable framework was due to the curated set of data transformation techniques that data scientists must perform in terms of imputation, categorical encoding of continuous variables or aggregation of healthcare datasets before building them into features to train a predictive model . Our framework will reuse these techniques as packages, which can be invoked in the deployment flow so that there is consistency in the way features are created for model training and model scoring. Thus, there is standardization of the training and deployment/scoring workflow, which helps in quickly learning through prospective testing the key components which can trigger data or feature drifts as the model runs in a real environment.

Doing this in the same controlled environment, which can ingest either historical or real-time data through the same APIs or secure connection, is key to productionize this within healthcare providers' and payers' systems to drive value. To achieve this, we hosted the entire framework in a secure HIPAA compliant cloud infrastructure to deploy as a turn-key solution as outlined below.

## Cloud Hosted Secure Healthcare Specific Turnkey Platform

*Isthmus* is hosted on a cloud-based infrastructure (Microsoft Azure Cloud Platform ^TM^). This infrastructure is set up with all the state-of-the-art functionalities, including network security, data replication, disaster recovery, and fault tolerance necessary for a robust and enterprise-grade software-as-a-service (SaaS). Cloud resources – compute and storage, both leverage economies of scale to keep the technology costs at a realistic level without having a need to maintain a large number of on-prem resources. Thus, being cost-effective as well as scalable and configurable, *Isthmus* can be adopted by healthcare organizations and systems of all sizes.

### Data Engineering Pipeline

Healthcare data by its very nature is highly complex, high dimensional. and of inconsistent quality. For this data to be useful, it needs a systematic data ingestion approach to collect, store and integrate data-driven insights into clinical and operational processes. To quickly ingest this multi-dimensional data and scale, we developed a configurable and flexible data ingestion pipeline solution, which enables the ingestion of all the relevant health data, including clinical data, claims data, Social Determinants of Health, and streaming IoT data. The data ingestion pipeline is also ready to allow the ingestion of genomics data and high-quality diagnostic imaging data as it becomes available.

The data ingestion pipeline is based on a simplified architecture (Figure 2), which enables user defined transformations for real-time scoring, cleaning, and de-duplication without requiring additional middleware. The raw data is pulled via RESTful API calls to the EHR's API servers or through regular intervals of data fetch using secure file transfer process. Generally, these API servers are the hub for all the API requests, which facilitates the connection between the EHR organizational users and the operational database management system to stream near real-time data seamlessly as a JSON response through the web service APIs upon service request.

The clinical data is mainly comprised of patient demographics, encounter related information, laboratory test result components, medication orders and laboratory vitals. These are pulled through APIs and are parsed through a configurable data munging script where the JSON objects are transformed using user defined templates that are specifically designed for each API response and are aggregated into patient level JSON responses. The pipeline is fully automated, and it ingests the patient level data in batch mode, where the batch size is based on SLA requirements. Thus, the pipeline is scheduled to trigger based on SLA requirements and it continuously pulls the data from the APIs and performs the desired transformation and filtering operations.

This patient level raw JSON data is then preprocessed using an imputation and filter logic, which transforms this data into clinically relevant features that are fed to the machine learning models using a scoring logic script to predict the risk of the acute care condition based on the pre-trained model.The scoring script generates the score response which encompasses the transformed features and the identified risk levels associated with the patient. These responses are aggregated in batch mode, and after cleaning, they are converted into SQL tables using a database operation script and ingressed into the PostgreSQL database. The data is stored in a very secure and reliable manner within the PostgreSQL database. The raw JSON responses are pushed to the Azure Data Lake^TM^ to preserve the raw patient level information for audit purposes.

###  Configurability, Extensibility and Experimentation

Predictive models in the healthcare industry are discrete in nature, especially for acute care systems , where we observe a high complexity of imputation logic and disparate features. The rest of the workflow has unified commonalities for all the models with identical evaluation patterns. If the configuration of these models is distributed and model-specific, it will add an overhead of maintaining and versioning multiple configurations for a given model as well as different environments. This will also increase the cost of additional infrastructure to manage different models.

The *Isthmus* platform provides a unique way of deploying and executing a model workflow for scoring using a single code base, which can support multiple models and versions using a single configuration file as shown in figure 3. It is also designed to use a single infrastructure cluster to execute any number of scoring workflow pipelines in parallel and automate the scoring process using Continuous Integration and Continuous Delivery processes (CI/CD). The use of the configuration methodology facilitates easy upgrading to an existing model or serving a new model in the pipeline workflow as it has a very short delivery cycle.

### Security and governance

All data collected in the context of healthcare is subject to specific privacy and regulatory requirements. Depending upon the type of data, and the context of the data collection and use, it can be Personally Identifiable Information (PII) or Protected Health Information (PHI). This information is often the target of security breaches. This is another reason why every access, analysis, and appending (if any) of this data needs to be audited, monitored, and documented. In the US, compliance with the HIPAA privacy law was mandated among all entities which either create or access PHI. This, in addition to the mandatory business continuity safeguards to avert a possibility of the health system facing outages of care services, puts additional constraint on the traditional feature and data engineering pipelines in a machine learning system. Challenges with data security through traditional on-prem systems and support of life-critical applications are solved by the *Isthmus* platform by adopting a cloud-first strategy. Additionally, the following Azure^TM^ enabled features utilized on *Isthmus* address the security concerns above.

**Access Control** The *Isthmus* platform leverages unique cloud-based security policies, such as the Azure active directory-based service principal for access control as an identity to manage applications and hosted services on the cloud and handle sensitive information (PHI). This eliminates the need for user level login to the cloud applications. Role Based Access Control (RBAC) uses Active Directory policies for managing the authentication. *Isthmus* provides a single role-based access to multi-institutional EHR data.

Additionally, the platform also provides a comprehensive, immutable log management service with easy access across deployed applications using elastic search and Kibana^TM^ dashboard. This ensures a single point of reference to test for any application or system level logs in a responsible manner. With the help of app-insight notifications, the *Isthmus* platform has provisioned real-time alerts for any configured event, including an exception in application or missing data from the source API to provide real-time alerts.

**Server maintenance** In a traditional healthcare setting, on-premise infrastructure solutions are constricted by cost and complexity of maintaining the hardware, firewall, software licensing, and the additional overhead of patch management software and its upgrades. This integration capability is limited to automate any patch management process and to roll out the access policy in a distributed project environment.

The *Isthmus* platform is engineered to easily automate with a cloud-based solution in a limited build and setup time. Moreover, there are prebuilt images developed on this platform which can be customized based on the application requirements. Thus, it makes our infrastructure scalable to meet the minimum performance requirements as well as enables the cloud applications to be easily migrated/replicated to another environment.

### Data replication and fault tolerance

Data replication and fault tolerance with an on-premise infrastructure needs enormous upfront investment. With an increase in processing and data scaling based on utilization, it becomes challenging to manage the business continuity needs. A distributed architecture application like a prediction model has serious challenges related to a single point of failure, which makes it substantially difficult to manage the continuous workflow without interruption when one or more of its components fail. This entails the need of the fault tolerant system design to prevent disruptions due to any single point of failure.

The *Isthmus* application is engineered to overcome these shortcomings and has the capabilities to scale up and accelerate the prediction model workloads to meet the needs of high-performance computing, low-latency, high-bandwidth network communication, and memory-intensive requirements. We adapted this cloud-based solution to resolve the lift and shift the problems like infrastructure upgrade, scalability, transfer and deployment at multiple locations using automated processes and containerization. This has considerably reduced the cost of infrastructure and engendered flexibility for migration/deployment on the cloud environments with minimal application level changes for the code, database, and data model architecture.

*Isthmus* applications have been developed with well-defined replication graphs and disaster recovery strategies for their database and support systems by imposing identical servers running in parallel replication with a mirrored backup of database and system level logs to ensure high levels of data availability. These applications are designed using the micro-services based architecture to reduce the redundancies from all the key components by performing similar activities in each workflow.

### Log management and alerting system

Logging is an essential component for any automated system as it serves as a flashlight on the black box module which is being deployed on a disturbed cluster. This logging information is generated in real-time on the *Isthmus* platform and helps in validating the stability of the system through warning and debug logs. This log data is fed to a high scale analytical engine (elastic search), which is a full-text search engine. This data is then neatly integrated with a visualization dashboard like Kibana^TM^ and provides feeds to self-hosted web front applications using restful APIs. This visualization provides monitors and performance metrics based on application level logs of the automated pipeline for Predictive and analytical applications. This also ensures quality delivery of the model serving on this platform and a quick debugging capability for any production outage.

For any production environment which is automated, having a notification system is critical given the fact that no workflow/infrastructure is perfect. In addition to the log management system, a slack based notification service is also integrated with the *Isthmus* platform to get real-time alerts about the production pipeline. This service has been deployed in real-time in production and is helping the engineering and the data science team to be fully aware of the live status of the pipeline and the patient risk scores. It is found to be useful to behave as the first line of action in the case of anomaly detection. The notification system captures both infrastructure and application failures/exceptions. Thus, this alerting system ensures immediate action and remediation in case of any failed events.

## Model Evaluation

The platform is designed to be a generic multipurpose data science engine. The flexible architecture of this platform has allowed the plug in and use of functional decision-making modules that can run asynchronously without disrupting the integrity of the system. The prediction service on the platform can be leveraged by the model evaluation service where real-time predictions can be interpreted by the models on the fly thereby making it extremely useful for the data scientists and clinicians (or stakeholders) to get actionable insights.

### Model Interpretability

The storage data lakes provisions access of historical data set while the model is running live in production. This data set is pushed to a model explainer script. This script helps in interpreting the model predictions by extracting the top contributing features that helped in making the predictions. This has been proved to be extremely useful by clinicians for making real-time decisions. It has also allowed them to provide necessary feedback to improve the existing models or generate newer models.

###  Silent Mode Testing

Before the model goes live in production, the model is validated in a silent-mode testing environment where the model performance is evaluated by the team of data scientists and clinicians to ensure that the model gives meaningful results and is ready for integration into actual patient care workflows. The *Isthmus* platform provisions this silent mode testing framework to follow all the best practices for the model deployment.

## Model Retraining and Development

In a real-time setting, it is always important that the model is highly stable and easily configurable. Thus, to address this problem, the platform allows retraining capabilities using the same data that was ingested into the model through APIs. The platform leverages this data and helps to generate multiple versions of the model which can be deployed by simply editing the model signature. The platform facilitates the data scientists to perform statistical tests to keep these models updated with the new upcoming data streams.

These features make this platform unique as the same infrastructure with few modifications can be used for model deployment and model development. In addition to the functionalities mentioned above, there is a re-usability of the code and interoperability of the system which makes this platform highly scalable. Thus, this *Isthmus* platform is an end-to-end system which encompasses all the software development processes into a robust function system.

## Multi-disciplinary team

Healthcare is not the only highly regulated industry, but it is an industry with a strong ethical, legal and moral underpinning to its regulations. Every clinical decision made either by a physician using either her judgement or the insights given to her by an ML system can have significant legal, financial, and/or life impacting consequences. Against this backdrop, there has been seminal work \[8\] conducted in demonstrating how individual-level data can be affected by historical prejudices and hence engender systematic biases in the decision of the ML system trained on it. Although there are many statistical methods to adjust for this unfairness, we believe that engaging clinical, operations and financial experts, and other stakeholders from the beginning is critical to ward off such biases and tune the algorithmic features to thresholds acceptable for the environment for which the model is being trained.

Involving clinicians helps in not only understanding the intimate relationship between different features, but it also gives a realistic understanding of issues related to physician cognitive overload \[9\] to improve user experience and help navigate actionable insights. Unlike being a fire-and-forget scenario, healthcare model development processes for clinical interventions are iterative in nature and need both retrospective as well as prospective testing \[10,11\]. To that effect, we adapted traditional statistical model development to a multi-faceted process, which combines pure statistical validation with clinical validation and brings in elements of clinical decision support systems.

# Discussions and Results

## Trauma case study (end-to-end framework)

Trauma surgeons work in a high-stress world of rapid decision-making. One of the most critical decisions they must make is whether a patient is stable enough to operate on, or if the Trauma team needs to focus on providing life-saving care. Trauma surgeons must make this decision quickly and accurately, as operating could save or take the life of their patient. Traditionally, prediction of survivability in trauma patients has always been based on risk assessment algorithms like TRISS, RTS, ISS \[16\]. However, all these algorithms are dependent on the Abbreviated Injury Scale (AIS), an anatomical scoring system that only considers a patient's injuries on arrival . There is not a standard system in place to track the trauma patient's status continuously over time, which is necessary for clinical decision-making as the patient's condition changes while in hospital.

In collaboration with Parkland Trauma Center Team, PCCI created a Trauma predictive model with the aim to monitor and guide clinical decision making during the first 12-72 hours of inpatient stay. This model accurately predicts mortality within the next 48 hours with 93 % accuracy in test data. However, to impact clinical decision-making

PCCI's Trauma predictive model needed a platform to bring the model to fruition with predictions made periodically every hour. Without integration into clinical workflows, PCCI's Trauma model predictions would just add noise to provider decision-making. It was critical that the Trauma Score would be easily accessible for each patient that the trauma surgeons needed to make a timely decision on how to intervene appropriately.

### *Isthmus* as the Solution

Every hour, up-to-the-second data is pulled directly from the EMR via an API and ingested securely into the *Isthmus* platform. In less than a minute, *Isthmus* cleans and transforms the API data, generates the Trauma Score, pushes it back to the EMR, and securely persists the data and results for later analysis. Providers can then view the score directly in the patient chart. *Isthmus* enables continuous, hourly predictions in the most critical window, 12-72 hours post-admission. Providers know they are always seeing a score based on up-to-date labs and vitals, so they can make decisions confidently. Using APIs and custom patient chart builds, the Trauma Score is delivered directly to the front page of the patient chart (an item called Storyboard that provides the most critical information for patient care in a column that persists as providers move through the patient chart). At a glance, providers can view and monitor the score, without a single click in the patient chart.

*Isthmus* has several abilities that have made it a perfect launchpad for the go-live of the Trauma model. These include its (1) ability to run near real-time predictive models end-to-end via data retrieval from APIs, (2)services to perform data cleaning, data transformations, feature engineering, scoring as per the need of the model, (3) functionality to integrate with EHRs to add the scores to patient chart, (4) mechanism to securely persist the data of the scored patients for future analysis, and (4) capability to monitor/ evaluate model's performance.

### Trauma Score Deployed

The Trauma Score has been running live at Parkland Hospital for 1 month, as of the time of publication. It is currently in a clinical validation stage, with providers using the score to validate their own assumptions. So far, the Trauma Score has been generated for 134 patients, correctly predicting 4 out of 4 mortalities. Anecdotally, providers say the score matches their own predictions, formed after a thorough read of the patient chart. Providers have indicated that promising results so far mean that soon they will begin to trust the score and use it for decision-making. This will save them critical minutes, allowing them to focus on operating or life-saving care over reviewing patient charts. These critical minutes could be the difference between life and death for high-risk patients. [^1]

## Community Data Initiative (CDI) on *Isthmus*

The healthcare industry is converging around a common realization that understanding the full context of an individual and the community where s/he lives is critical to improving outcomes and lowering costs. Factors such as transportation access, safe neighborhoods, food availability, school performance, and economic opportunity have always been part of the conversation around healthcare but now, there is an urgency to quantify their impact on the health of individuals and communities. However, data for these social determinants of health is not easily accessible. Documenting social determinants of health is typically not part of the standard documentation workflow during an individual's clinical encounter. And zip code level indicators are not specific enough to make meaningful interventions.

With help from partners across the Dallas community, we utilized the *Isthmus* platform to develop a block group level, longitudinal collection of about 60 key indicators that can be used to measure health, resiliency, and economic vibrancy of neighborhoods. We made these indicators available in an interactive format through a dashboard hosted on the same platform.

By hosting everything on a single platform, we have been able to create a streamlined process for ingesting, cleaning, analyzing, and visualizing data. Through an automated pipeline using Apache NiFi, data is ingested into the blob storage, where data is cleaned according to predefined scripts. This cleaned data is eventually ingested to the PostgreSQL database management system, which is a free and open-source database. The dashboard tool pulls data from the PostgreSQL database management system.

To quantify the impact of CDI data on specific risk conditions, we evaluated its ability to predict the risk of pre-term birth when combined with clinical and claims data. This predictive model was developed at PCCI and uses social determinants of health (SDOH) features extracted from the CDI dataset. These SDOH features proved to be significant for making predictions with this model. Thus, the data from the CDI platform is helping the data scientists and machine learning engineers demonstrate a holistic view of factors that impact preterm birth. By including the data from external sources into a model that uses clinical data, the overall model performance has improved.

## Internet-of-Things on *Isthmus*

The use of internet of things (IoT) applications in the healthcare industry is rapidly growing. The *Isthmus* platform has piloted the IoT data ingestion pipeline with indoor air quality sensors, where real-time data from these sensors' APIs were ingested into the platform. This data was massaged and munged in batch mode using the data cleaning modules developed on the platform.

There was high reusability of the configuration which was developed for the machine learning system as the codebase was made generic to reduce redundancies across multiple projects on the platform. The IoT data is stored and maintained in a PostgreSQL database on the platform with fault tolerance and disaster recovery functionalities. The futuristic goal is to integrate this rich IoT data with the existing machine learning models to have additional features which are believed to improve the model predictions.

# Conclusions

This paper introduced the *Isthmus* Platform, an end-to-end system for developing and deploying ML models. Using *Isthmus*, data scientists can use their familiar ML toolkits and libraries to create models, perform statistical tests, and deploy them.

We addressed the challenges of integrating ML models into application development and model sharing. The proposed architecture supports the sharing of pretrained models across different ML modules run-time environments. As illustrated by the case studies, *Isthmus* provides project level isolation and focuses on code reusability, rather than reinventing the development pieces. The platform demonstrated versatility in terms of serving a prediction service, ingesting IoT data, and integrating an SDOH feature engine.

In the future, we hope to (1) integrate real-time model training with streaming data by tuning the model- coefficients on the fly and (2) perform A/B testing and multi-armed bandit testing to ensure the model stability and reliability over time.

# Acknowledgments

We are extremely grateful to Lyda Hill and the Parkland Foundation for providing financial support to design and build this platform to serve the needs of the New Parkland Hospital and its patients. We are also grateful to the leadership team at the Parkland Health and Hospital System for providing strategic guidance regarding how to best maximize the platform's usefulness to Parkland and to the Parkland IT team for helping with various aspects of the platform build and integration with the EHR.

We would also like to express our deep gratitude to Dr. Manjula Julka and Dr. George Oliver, our clinician supervisors, for their patient guidance, enthusiastic encouragement, and useful critiques of this work. We would also like to thank Dr. Shelley Chang, for her providing her clinical expertise, advice, and assistance. Our grateful thanks are also extended to the PCCI leadership team for their support and encouragement throughout this journey. Finally, we would like to thank everyone on the *Isthmus* team for making *Isthmus* a success.

# Disclaimer

Trademark products mentioned in this paper are properties of respective trademark owners.

[^1]: CDI was supported in part through a grant from Community Council of Greater Dallas (CCGD) and was developed collaboratively with CCGD, DFW Hospital Council Foundation and University of Texas at Dallas