• We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.
• We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.
• ### In-flight Calibration of Hitomi Soft X-ray Spectrometer (3) Effective Area(1801.02104)

Jan. 7, 2018 astro-ph.IM, astro-ph.HE
We present the result of the in-flight calibration of the effective area of the Soft X-ray Spectrometer (SXS) onboard the Hitomi X-ray satellite using an observation of the Crab nebula. We corrected for the artifacts when observing high count rate sources with the X-ray microcalorimeter. We then constructed a spectrum in the 0.5-20 keV band, which we modeled with a single power-law continuum attenuated by an interstellar extinction. We evaluated the systematic uncertainty upon the spectral parameters by various calibration items. In the 2-12 keV band, the SXS result is consistent with the literature values in flux (2.20 $\pm$ 0.08) $\times$10$^{-8}$ erg s$^{-1}$ cm$^{-2}$ with a 1$\sigma$ statistical uncertainty) but is softer in the power-law index (2.19 $\pm$ 0.11). The discrepancy is attributable to the systematic uncertainty of about $+$6/$-$7% and $+$2/$-$5% respectively for the flux and the power-law index. The softer spectrum is affected primarily by the systematic uncertainty of the Dewar gate valve transmission and the event screening.
• ### Inflight Calibration of the Hitomi Soft X-ray Spectrometer (2) Point Spread Function(1712.08321)

Dec. 22, 2017 astro-ph.IM
We present results of inflight calibration of the point spread function (PSF) of the Soft X-ray Telescope (SXT-S) that focuses X-ray onto the pixel array of the Soft X-ray Spectrometer system (SXS). We make a full array image of a point-like source by extracting a pulsed component of the Crab nebula emission. Within the limited statistics afforded by an exposure time of only 6.9~ksec and the limited knowledge of the systematic uncetainties, we find that the raytracing model of 1'.2 half-power-diameter (HPD) is consistent with an image of the observed event distributions across pixels. The ratio between the Crab pulsar image and the raytracing shows scatter from pixel to pixel that is 40% or less in all except one pixel. The pixel-to-pixel ratio has a spread of 20%, on average, for the 15 edge pixels, with an averaged statistical error of 17% (1 sigma). In the central 16 pixels, the corresponding ratio is 15% with an error of 6%.
• We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.
• ### X-ray Measurements of the Particle Acceleration Properties at Inward Shocks in Cassiopeia A(1710.06992)

Dec. 6, 2017 astro-ph.HE
We present new evidence that the bright non-thermal X-ray emission features in the interior of the Cassiopeia A supernova remnant (SNR) are caused by inward moving shocks based on Chandra and NuSTAR observations. Several bright inward-moving filaments were identified using monitoring data taken by Chandra in 2000-2014. These inward-moving shock locations are nearly coincident with hard X-ray (15-40 keV) hot spots seen by NuSTAR. From proper motion measurements, the transverse velocities were estimated to be in the range $\sim$2,100-3,800 km s$^{-1}$ for a distance of 3.4 kpc. The shock velocities in the frame of the expanding ejecta reach values of $\sim$5,100-8,700 km s$^{-1}$, slightly higher than the typical speed of the forward shock. Additionally, we find flux variations (both increasing and decreasing) on timescales of a few years in some of the inward-moving shock filaments. The rapid variability timescales are consistent with an amplified magnetic field of $B \sim$ 0.5-1 mG. The high speed and low photon cut-off energy of the inward-moving shocks are shown to imply a particle diffusion coefficient that departs from the Bohm regime ($k_0 = D_0/D_{\rm 0,Bohm} \sim$ 3-8) for the few simple physical configurations we consider in this study. The maximum electron energy at these shocks is estimated to be $\sim$8-11 TeV, smaller than the values of $\sim$15-34 TeV inferred for the forward shock. Cassiopeia A is dynamically too young for its reverse shock to appear to be moving inward in the observer frame. We propose instead that the inward-moving shocks are a consequence of the forward shock encountering a density jump of $\gtrsim$ 5-8 in the surrounding material.
• Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.
• The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.
• ### Smoothness-constrained model for nonparametric item response theory(1704.07736)

April 25, 2017 stat.AP
This paper is concerned with the nonparametric item response theory (NIRT) for estimating item characteristic curves (ICCs) and latent abilities of examinees on educational and psychological tests. In contrast to parametric models, NIRT models can estimate various forms of ICCs under mild shape restrictions, such as the constraints of monotone homogeneity and double monotonicity. However, NIRT models frequently suffer from estimation instability because of the great flexibility of nonparametric ICCs, especially when there is only a small amount of item-response data. To improve the estimation accuracy, we propose a novel NIRT model constrained by monotone homogeneity and smoothness based on ordered latent classes. Our smoothness constraints avoid overfitting of nonparametric ICCs by keeping them close to logistic curves. We also implement a tailored expectation--maximization algorithm to calibrate our smoothness-constrained NIRT model efficiently. We conducted computational experiments to assess the effectiveness of our smoothness-constrained model in comparison with the common two-parameter logistic model and the monotone-homogeneity model. The computational results demonstrate that our model obtained more accurate estimation results than did the two-parameter logistic model when the latent abilities of examinees for some test items followed bimodal distributions. Moreover, our model outperformed the monotone-homogeneity model because of the effect of the smoothness constraints.
• ### Investigating consumers' store-choice behavior via hierarchical variable selection(1704.00665)

April 3, 2017 stat.AP
This paper is concerned with a store-choice model for investigating consumers' store-choice behavior based on scanner panel data. Our store-choice model enables us to evaluate the effects of the consumer/product attributes not only on the consumer's store choice but also on his/her purchase quantity. Moreover, we adopt a mixed-integer optimization (MIO) approach to selecting the best set of explanatory variables with which to construct a store-choice model. We devise two MIO models for hierarchical variable selection in which the hierarchical structure of product categories is used to enhance the reliability and computational efficiency of the variable selection. We assess the effectiveness of our MIO models through computational experiments on actual scanner panel data. These experiments are focused on the consumer's choice among three types of stores in Japan: convenience stores, drugstores, and grocery supermarkets. The computational results demonstrate that our method has several advantages over the common methods for variable selection, namely, the stepwise method and $L_1$-regularized regression. Furthermore, our analysis reveals that convenience stores tend to be chosen because of accessibility, drugstores are chosen for the purchase of specific products at low prices, and grocery supermarkets are chosen for health food products by women with families.
• ### Multi-Year X-ray Variations of Iron-K and Continuum Emissions in the Young Supernova Remnant Cassiopeia A(1701.07905)

Jan. 26, 2017 astro-ph.HE
We found simultaneous decrease of Fe-K line and 4.2-6 keV continuum of Cassiopeia A with the monitoring data taken by Chandra in 2000-2013. The flux change rates in the whole remnant are $-$0.65$\pm$0.02 \% yr$^{-1}$ in the 4.2--6.0 keV continuum and $-$0.6$\pm$0.1 \% yr$^{-1}$ in the Fe-K. In the eastern region where the thermal emission is considered to dominate, the variations show the largest values: $-$1.03$\pm$0.05 \% yr$^{-1}$ (4.2-6 keV band) and $-$0.6$\pm$0.1 \% yr$^{-1}$ (Fe-K line). In this region, the time evolution of the emission measure and the temperature have a decreasing trend. This could be interpreted as the adiabatic cooling with the expansion of $m = 0.66$. On the other hand, in the non-thermal emission dominated regions, the variations of the 4.2--6 keV continuum show the smaller rates: $-$0.60$\pm$0.04 \% yr$^{-1}$ in the southwestern region, $-$0.46$\pm$0.05 \% yr$^{-1}$ in the inner region and $+$0.00$\pm$0.07 \% yr$^{-1}$ in the forward shock region. In particular, the flux does not show significant change in the forward shock region. These results imply that a strong braking in the shock velocity has not been occurring in Cassiopeia A ($<$ 5 km s$^{-1}$ yr$^{-1}$). All of our results support that the X-ray flux decay in the remnant is mainly caused by the thermal components.