Appendix

For the fits discussed within this paper, we have used a combination of binning by signal-to-noise and uniform channel binning. Binning was employed to allow the use of $ \chi ^2$ statistics, which are faster to minimize (this is a non-trivial concern when using the warmabs model, which is an extremely slow to evaluate code). Furthermore, adding of data and binning of channels can serve to average over systematic calibration uncertainties. However, any non-uniform binning, e.g., a signal-to-noise criterion, can introduce biases in line fits.

For our data, the signal-to-noise per channel rapidly varies at energies $ < 1$keV, and there is no (small) set of uniform channel binnings that achieves adequate signal-to-noise in these channels. Thus, we choose mixed criteria (using the ISIS group  function) that ensures both a minimum signal-to-noise and a minimum number of channels in the binned data. In practice for these particular data, only the minimum channel criterion applies above 1keV, as this is sufficient to ensure signal-to-noise $ \ge 4$ in all of these channels.

Figure: Left: Results of the blind line search, utilizing 51 components. The first 24 are presented in Table 7, along with possible line IDs, presuming an outflowing wind (in the cosmological rest frame of PG1211+143 ) with outflow velocity $ v_{\rm out} = -17\,300~\rm km\,s^{-1}$ ( $ z_{\rm out}\sim -0.0561$). Lines labeled with question marks are those found within the first 24 components for which we have no suggested line ID. Right: Chandra 1st-order detector counts stacked in velocity bins centered on the nine identified lines in Table 7. The dashed lines are the velocities of absorber components suggested by Pounds et al. (2016a), as well as the $ v_{\rm out} = -17\,300~\rm km\,s^{-1}$ ( $ z_{\rm out}=-0.0561$) absorber component discussed in this work.
\includegraphics[width=1.0\textwidth, viewport=35 155 575 750, clip, angle=0]{figures/fig10_spectral_ids.ps}
\includegraphics[width=1.0\textwidth, angle=0]{figures/fig10_vshift.ps}


Table: Results of Blind Line Search

$ \lambda_{\rm obs}$
$ \lambda_{\rm rest}$ $ \sigma$ $ \Delta C$ $ -$EW Component ID $ \lambda_{\rm lab}$ $ z_{\rm out}$
(Å) (Å) (Å)   (mÅ)     (Å) (c)

4.142
3.832 0.013 $ -$11.5 $ -$12.6 11 - - -
6.173 5.830 0.002 $ -$13.1 $ -$6.9 6 SiXIV (Ly$ \alpha $) 6.182 $ -0.0586$
6.793 6.285 0.014 $ -$12.4 $ -$10.9 4 SiXIII (He$ \alpha $) 6.648 $ -0.0565$
8.219 7.603 0.009 $ -$7.0 $ -$7.6 21 - - -
8.613 7.968 0.026 $ -$10.9 $ -$16.8 10 MgXII (Ly$ \alpha $) 8.421 $ -0.0555$
8.913 8.418 0.008 $ -$13.1 $ -$12.8 3 - - -
9.160 8.651 0.002 $ -$9.8 $ -$9.0 13 MgXI (He$ \alpha $) 9.169 $ -0.0576$
10.504 9.920 0.003 $ -$7.2 $ -$9.8 20 - - -
12.089 11.417 0.003 $ -$9.9 22.7 19 - - -
12.133 11.459 0.022 $ -$30.0 $ -$55.8 0 NeX (Ly$ \alpha $) 12.133 $ -0.0576$
13.453 12.705 0.013 $ -$8.3 $ -$31.7 16 NeIX (He$ \alpha $) 13.447 $ -0.0565$
13.909 13.136 0.002 $ -$10.1 35.7 15 - - -
15.140 14.007 0.017 $ -$13.6 $ -$50.4 5 - - -
16.338 15.430 0.009 $ -$12.5 $ -$52.3 9 - - -
16.518 15.600 0.009 $ -$12.5 $ -$50.0 8 - - -
16.691 15.763 0.002 $ -$8.7 115.5 17 - - -
17.739 16.753 1.132 $ -$15.6 141.7 2 - - -
18.963 17.544 0.010 $ -$10.1 $ -$61.2 12 OVII (He$ \beta$) 18.627 $ -0.0597$
22.080 20.428 0.025 $ -$7.5 $ -$131.6 23 OVII (He$ \alpha $) 21.602 $ -0.0555$
23.933 22.603 0.004 $ -$8.2 117.5 22 - - -
24.431 23.073 0.600 $ -$22.4 $ -$50.5 1 - - -
24.855 23.473 0.129 $ -$10.8 $ -$354.3 14 NVII (Ly$ \alpha $) 24.781 $ -0.0544$
25.584 24.161 0.130 $ -$7.7 $ -$437.8 18 - - -
29.610 27.964 0.099 $ -$13.2 2631.8 7 - - -

\begin{tablenotes}
\item[1]\textbf{Notes.} Results from a blind search to the un...
...hift in the
rest frame of PG\,1211+143~\ if this ID is correct.
\end{tablenotes}

Our fits indicating an outflowing wind with a velocity in the PG1211+143 frame of $ -$17300 kms$ ^{-1}$ ( $ z_{\rm out}\sim -0.0561$) are primarily being driven by the $ \alpha $ lines of Helium-like and Hydrogen-like Ne, Mg, and Si. To ensure that these lines are not being biased by our binning (although five of the six primary lines in our fits are within the uniformly binned portion of the spectrum at energies $ >1$keV), we have also performed a “blind line search” of the data. We grid the HEG data to the MEG bins, and combine all spectra, but otherwise do not perform any further channel binning. We use the same continuum model as discussed above, namely an absorbed disk plus powerlaw spectrum with exponential cutoff, and include line emission from the Fe region. We then loop through the spectra, adding one line at a time which is allowed to freely range between emission and absorption. The initial line fit is constrained to a narrow range of wavelengths (16 MEG channels, i.e., $ \approx 0.18$Å), but all possible wavelength bins are searched. The line with the greatest change in fit statistic is retained. After each step, all continuum and line parameters are refit (within the constraints of the existing wavelength region of the added line). This process is repeated (51 times in Figure 10, with the first 24 found lines listed in Table 7). As the goal is to identify candidate lines, we do not calculate confidence intervals for the final lines.

The putative NeX line is our single most significant residual, with the remaining 5 lines from H-like and He-like Ne, Mg, and Si all falling within the 17 most significant residuals. Additionally, there are three other residuals that might be associated with a $ z_{\rm out} \approx -0.0561$ outflow. These, however, fall within a very low signal-to-noise region of the spectrum. Several residuals are broad, and are undoubtedly modifying the continuum fit. Several could be spurious noise features. The rest remain unidentified. However, this blind search highlights those features that are driving the XSTAR and warmabs models to identify an outflowing component in the rest frame of PG1211+143 .

As a further diagnostic of possible absorber components, we take the 1st-order Chandra -HETGS counts, and using the 9 potential H- and He-like lines identified in Table 7, stack the data into cosmological rest frame velocity bins. The results of this stacking are also shown in Figure 10 (note that in this procedure, not every bin is statistically independent from one another, as counts can be reused for different ions). The strong feature at a velocity in the rest frame of PG1211+143  of $ -$17300 kms$ ^{-1}$ ( $ z_{\rm out}\sim -0.0561$) is apparent. We also indicate the velocities of the absorber components suggested by the analysis of Pounds et al. (2016a). There is a feature near the $ \sim-0.19c$ ( $ -57\,000~\rm km\,s^{-1}$) velocity found by Pounds et al. (2016a), but we have not found any single line-like residual at such a velocity. It is possible, however, that such an absorption component, if real, only manifests significantly in a stacked analysis.

Ashkbiz Danehkar
2018-03-28