An international x-ray satellite, Hitomi, led by Japan was launched on February 17, 2016, by an H-IIA rocket at the Tanegashima Space Center in Japan and placed in a low-Earth orbit with an altitude of 575 km and an inclination angle of 31 deg.1 Hitomi carries four types of instruments covering a wide energy range from soft x-ray to soft gamma ray. The hard x-ray imaging system composed of two sets of the Hard X-ray Imagers (HXI)2 and two sets of the Hard X-ray Telescopes (HXT)3 is capable of imaging spectroscopy in the hard x-ray band ranging from 5 to 80 keV.2 The two HXI systems are referred to as HXI1 and HXI2, individually paired with HXT1 and HXT2, respectively. Due to the focusing optics, the sensitivity of the hard x-ray imaging system for the point source is 100 times better than those of nonfocusing instruments in the hard x-ray bands, such as hard x-ray detector onboard Suzaku (Suzaku/HXD).4
The HXI is composed of a stacked semiconductor imager5–8 and active shields surrounding the imager. The imager consists of five layers of double-sided strip detectors (DSD) with a strip pitch of and detector area of . The upper four layers are the double-sided Si strips detectors (DSSDs) with a thickness of ,9–11 and the bottom layer is the CdTe double-sided strip detector (CdTe-DSD) with a thickness of .12–16 The DSSDs have p- and n-type strips on the surface of the top and bottom sides of an n-type Si wafer, while the CdTe-DSDs have Pt- and Al-strips on those of a p-type CdTe wafer. By applying positive bias voltages to the n-side of the DSSDs and the Al-side of the CdTe-DSDs, holes and electrons generated by the incident x-ray photon are collected by the p-side/Pt-side strips and the n-side/Al-side strips, respectively. The x-ray-induced charge on the strip electrodes is read-out utilizing dedicated low-noise front-end application-specific integrated circuits (ASICs),17 which are connected to the individual strips. The active shields consist of nine (BGO) scintillators arranged as a well-type structure. Their thicknesses are typically in order to stop protons with energy , trapped at the South Atlantic Anomaly (SAA). The scintillation light of each BGO is read-out by an avalanche photodiode (APD).18 By processing the read-out signals from the APDs in digital filters, veto signals are generated and used for reducing the detector background.19
Although the HXI was lost,2 thorough investigations and evaluations of its in-orbit performance are of great importance for planning and designing future hard x-ray missions. In this paper, we describe the in-orbit performance and calibration results of the HXI. In Sec. 2, in-orbit operations and functionalities are summarized. Standard analysis method of the HXI is described in Sec. 3. Detailed performances on the non-x-ray background (NXB) and energy response are presented in Secs. 4 and 5.
Initial Operations and Observations
After the deployment of the extensible optical bench on February 28, 2016, the temperature of the HXI was gradually cooled to the operation temperature of . On March 8, 2016, a start-up operation of the HXI started. High voltages of the APDs and the DSSD/CdTe-DSD were applied one by one and reached the optimum values on March 12 for HXI1 and March 14 for HXI2.
After being turned on, the HXI performed several observations as listed in Table 1. In spite of the short lifetime of the HXI, x-ray photons from three astronomical objects (IGR J16318-4848, G21.5–0.9, and Crab Nebula) were successfully detected. As well as these data, the HXI observed 164.3 ks (HXI1) and 163.8 ks (HXI2) of blank sky data and 158.7 ks (HXI1) and 160.5 ks (HXI2) of Earth occultation data, including both bright and night Earth. According to the hard x-ray observations by “Swift,”20 count rate on the HXI due to the albedo x-ray/gamma rays is estimated to be less than . Since this count rate is negligible compared with the background rate as shown in Sec. 5, the Earth occultation data are referred to as the NXB in this paper. On the other hand, the summed data of observation sequences named “None2,” “IRU check out,” and “RXJ1856.5–3754,” which includes the cosmic x-ray background (CXB), are referred to as blank sky. These background data provide fruitful information on the in-orbit background properties in the hard x-ray energies as described in Sec. 5.
Observation log of the HXI.
|Start time||Stop time||OBSID||Target name||Notes|
|03-11 21:24||03-13 17:56||10042020–10042030||IGR J16318–4848||HXI1 on|
|03-13 17:56||03-14 16:20||100042040||IGR J16318–4848||Stray light|
|03-14 16:20||03-14 18:00||000007010||None2||—|
|03-14 18:00||03-15 17:56||000007020||None2||HXI2 on|
|03-15 17:56||03-16 19:40||000008010–000008060||IRU check out||—|
|03-16 19:40||03-19 17:00||100043010–100043040||RXJ1856.5–3754||—|
|03-19 17:00||03-23 13:30||100050010–100050050||G21.5–0.9||—|
|03-23 13:30||03-24 11:22||100043050||RXJ1856.5–3754||—|
|03-24 11:22||03-25 11:28||100043060||RXJ1856.5–3754||DTHRa changed|
|03-25 11:28||03-25 18:01||000007010–000007020||Crab||—|
In orbit, all the basic functions of the HXI worked properly. Here, functionalities of the HXI in flight are briefly summarized. For more details, please refer to Nakazawa et al.2 In the imagers, there was no damage or degradation due to the launch. All read-out channels of all the ASICs worked properly, and noise levels were consistent with the ground calibration. The energy resolution was evaluated by fitting the onboard calibration source spectra, which was mounted just above the top layer of the DSSD. From these data, good energy resolutions of 1.0 keV at 13.9 keV and 2.0 keV at 59.5 keV in full-width at half-maximum were obtained. Also, the energy gain was very stable within an uncertainty of bin of the pulse-height invariant (PI) at 59.5 keV, corresponding to 0.1 keV or .
The active shields also showed good performances. The low-energy threshold of each BGO scintillator was the same as the ground calibration results, and anticoincidence rate was consistent with the prelaunch estimation. Light curves of the veto signals from the active shields clearly showed variability corresponding to the geomagnetic cutoff rigidity and decay of the activation component after passages of the SAA. It indicated that the active shields properly monitored the variability of cosmic-ray in the geomagnetic cutoff rigidity and SAA.
In this section, the standard analysis method for the HXI data is summarized. It is composed of four steps: gain correction, event reconstruction, screening, and dead-time correction.
In the gain-correction process, gain-corrected Energy PI (EPI; PI in units of keV, in real number) is calculated. The EPI is generated from the raw analog to digital converter (ADC) value in the gain-correction process and used only in the event reconstruction process. After the event reconstruction, the EPI is converted into PI (in integer number), which is used in further screening and scientific analysis. The relation between PI and EPI is expressed as . For example, ( line) is converted to .
At the first step of the gain correction, all signals from bad channels are excluded from the following processes. In the flight models, only strips located at the edge of the detector, where the leakage current is higher than the other strips, are defined as bad channels. Thus, the detector area within () is available for imaging spectroscopy, which corresponds to . Then, the common mode noise is subtracted from ADC values, before correcting the gain. The common mode noise is a noise where all channels in one read-out ASIC coherently fluctuate. It is estimated in the ASICs by recording the median ADC value, which is 16th smallest in all channels in one ASIC. By the common mode subtraction, the pedestal level of each channel is corrected to zero. Finally, the ADC values of good strips are converted to EPI with third-order polynomial functions.
The gain-correction functions are determined based on the ground calibrations conducted in December and October 2014 for HXI1 and HXI2, respectively. The two HXIs were operated in a low-temperature chamber at Institute of Space and Astronautical Science, where x-ray/gamma-ray photons from radioactive isotopes , , , and irradiated the instrument. From these data, the correspondence between ADC values and photon energies for x-ray and gamma-ray lines is obtained and listed in Table 2. Between these lines, the ADC-energy correspondence is interpolated with third-order spline functions, which are used as the gain-correction function.
X-ray and gamma-ray lines used for gain calibration.
|0||Top||5.9, 13.9, 17.8, 20.8, 26.3, 30.8, 35.0, 59.5, 81.0, 122a|
|Bottom||5.9, 30.8, 59.5, 81.0a, 122a|
|1||Top||13.9, 17.8, 20.8, 26.3, 30.8, 35.0, 59.5, 81.0, 122a|
|Bottom||30.8, 59.5, 81.0a, 122a|
|2||Top||13.9, 17.8, 20.8, 26.3, 30.8, 35.0, 59.5, 81.0, 122a|
|Bottom||30.8, 59.5, 81.0a, 122a|
|3||Top||13.9a, 17.8, 20.8, 30.8, 35.0, 59.5, 81.0, 122a|
|Bottom||30.8, 59.5, 81.0a, 122a|
|4||Top||17.8a, 30.8, 35.0, 59.5, 81.0, 122|
|Bottom||17.8a, 30.8, 35.0, 59.5, 81.0, 122, 136a|
To obtain photon information from the gain-corrected signals in each data acquisition, event reconstruction processes are essential because the HXI imager consists of a stacked double-sided detector. In DSDs, data acquisition for one x-ray photon event usually consists of at least two signals, one each from both sides, and these signals often split into two adjacent strips. Moreover, some signals are detected in multiple layers due to Cd/Te fluorescence lines or Compton scattering. To identify these multisignal events, all the signals exceeding the digital thresholds (DTHR) are read-out simultaneously from all ASICs in all layers even when only one strip in one layer generates the trigger signal. In the standard HXI analysis, the events are reconstructed as shown in Fig. 1: the gain-corrected EPI in each side is obtained at first, then they are combined within one layer, and then combined with information from other layers to finally reconstruct a photon event.
In the first step, all the signals below analysis thresholds are discarded. The analysis thresholds are larger than DTHR for almost all of the strips and are set to individual strips, to be four times (DSSDs) or six times (CdTe-DSDs) the standard deviation of pedestal peaks, which correspond to events with zero energy. Thus, the analysis threshold of the DSSDs is much lower than that of the CdTe-DSDs in order to lower the HXI energy range as far as possible. Under these settings, pedestals from all 126 active strips are below the analysis threshold with probabilities of 99.6% for DSSDs and more than 99.9999% for CdTe-DSDs. It means that the noise contaminates with probabilities of 0.4% for DSSD and for CdTe-DSD. Mean values of the analysis thresholds of all the strips in the bottom sides (-sides) of the top-layer DSSD are 3.56 keV for HXI1 and 3.66 keV for HXI2, and they typically distributes from 3 to 4 keV. These values determine lower limits of the energy range of the HXI because noise levels of the bottom sides are worse than those of the top sides in the DSSDs. Thus, the HXI can observe 5 keV in almost all the strips, while there are two strips with analysis thresholds exceeding 5 keV in HXI2. After applying the analysis thresholds, only signals from the single strip or two adjacent strips are accepted, and signals from more than two strips and those from nonadjacent two strips are discarded. These events are to 3% of all the events in the ground data using radioisotopes.
Since the HXI is composed of double-sided detectors, signals from top and bottom sides of the layer must be combined to obtain two-dimensional (2-D) positional information. EPI from the top side in the DSSDs and that from the bottom side in the CdTe-DSD has better energy resolutions and, thus, is assigned as the EPI value of each layer. Position is simply determined using an intersection point of strips in both sides. When two adjacent strips have signals, the strip with larger pulse height is assigned.
When combining the signals from both sides, consistency between EPI values from both sides is checked using a condition shown in red lines in Fig. 2. Non-x-ray signals by the instrumental noise or a certain cosmic particles can be rejected by this consistency check. The condition for this check is that the pulse heights (EPI) from both sides match within of the energy resolution. Specifically, it is written as21,22 for both Si and CdTe, the second parameters are calculated as (Si) and 0.00044 keV (CdTe), where the pair production energies of (Si) and 4.4 eV (CdTe)12 are used. In addition to these parameters, the low mobility of holes in CdTe is taken into account as a parameter by assuming that every 7.1% of charges are lost during the drift toward the Pt-side strips from the incident position. Thus, for CdTe-DSD and (complete charge collection) for DSSD are assumed. In DSSDs, subpeak events due to the nonuniform electric field,9 which is described in Sec. 4, are also discarded by this process. In the on-ground calibration experiment using radioisotope, 2% to 3% of total events in DSSDs and in CdTe-DSDs are discarded.
After finishing the event reconstruction processes in one layer, hits in five layers are reconstructed as a photon event. In this process, single-hit events detected in a single layer and double-hit events at the combination of one CdTe-DSD and one DSSD with an energy of DSSD consistent with a fluorescence line of Cd or Te are accepted. Otherwise, no values are assigned to PI, hit positions in the final event list and hence discarded in the following processes. In terms of physical processes, this algorithm accepts photoabsorption events and fluorescence escape events, where K-shell fluorescence photons of Cd or Te escaped from CdTe-DSDs are photoabsorbed in DSSDs. Compton scattered events are ignored in current implementation because a fraction of such events composed of Si–Si double hits or Si-CdTe nonfluorescence double hits in total events are less than in the ground data obtained with and radioisotopes.
The event reconstruction algorithm described above must be tested with in-flight data because it was determined based on the ground data analysis. For the purpose of investigating whether this algorithm properly rejects the background data without excluding much of the real x-ray signals, fractions of accepted events (black) and discarded events (red and blue) of the ground calibration data (trigger rate ), Crab data (trigger rate ), and the NXB data (trigger rate ) are shown in Fig. 3. In this figure, the denominator of the fractions is the number of events, in which at least one signal exceeds the analysis thresholds. The nonsignal events account for of ground calibration data, of Crab data, and of the NXB data. These nonsignal events are thought to originate from noise triggers and soft photons below the analysis thresholds because trigger thresholds are set to be as low as possible within a range where the dead-time fraction due to noise triggers does not affect the scientific observations. The reason it is very high in the NXB data is simply because its trigger rate by the external photon and particle background is much lower than those by the instrumental noise (typically a few Hz).
From Fig. 3, the fraction of discarded events over all nonzero signal events is much larger in the NXB data than the ground calibration and Crab data. All inconsistent EPI, nonadjacent 2-strip and -strip are contained in the NXB data, indicating non-x-ray origins of the signal. On the other hand, more than of events are accepted in ground calibration and Crab data, which are presumably dominated by x-ray signals. Thus, in other words, the multilayer nature of the imager and the event reconstruction procedure using their information are effective in reducing the background.
In addition to the event reconstruction process, bad events (e.g., veto events) and bad time intervals (e.g., SAA passages), which are presumably dominated by the instrumental noise or the NXB, are excluded both in the onboard software and the ground analysis software. Basically, the in-flight screening is less stringent than the on-ground screening to flexibly change the screening conditions after observations, on the ground.
The in-flight screening of the HXI is performed in ASICs and HXI Digital Electronics (HXI-DE). In ASICs, only signals exceeding a predefined DTHR are read-out. DTHR is adjustable for each ASIC independently and is also independent from the trigger threshold. On the day before the Crab observation, it was raised up to similar level to the ground-software analysis thresholds for reducing the data size. Read-out data from the ASICs are reduced by further screening in HXI-DE. It assigns “CATEGORY” of high, middle, and low to each event. Assignment of the CATEGORY is performed using the time interval from the previous trigger, number of signals above DTHR, ADC values, flags of active shield coincidence, trigger pattern, and the other flags from the ASIC. This CATEGORY determines priorities to record the event to the data recorder (DR) of the satellite. Since the capacity of the DR is limited, most of the data in middle and low categories are not downloaded to the ground except for those obtained within the interval the satellite is in direct contact with from the operation site at Uchinoura in Japan.
In all observations of the HXI, the screening criteria as listed in Table 3 are applied in HXI-DE. If an event has a flag for the calibration source signal, pseudotrigger, test pulse, or forced trigger, CATEGORY = High is assigned automatically. Then, CATEGORY = Low is assigned to an event if it has a flag of fast BGO, number of signals exceeds 31, which cannot happen with x-ray photons, or any of signals has ADC value of 1023, which is the upper limit of ADC. In the remaining events, if a signal has a flag of HITPAT (hitpattern) BGO, CATEGORY = Middle is assigned. While the fast BGO signal is a veto signal capable to stop ADC conversion if needed, the HITPAT BGO signal is slow but a lower threshold veto signal is to be used for further background rejection in the on-ground data screening (see Ohno et al.19 for more details). Finally, CATEGORY = High is assigned to all the remaining events. In this screening, only events that are almost certainly the background signals are classified to middle or low categories.
Criteria for each CATEGORY.
|High||Calibration source, pseudotrigger, test pulse or forced trigger|
|Not assigned to middle or low|
|Middle||Flagged as HITPAT BGO and not flagged as fast BGO|
|Low||Flagged as fast BGO|
|Number of signals|
|Any of signal has|
The in-flight screening criteria by HXI-DE should also be tested with the Crab and NXB data. In the Crab observation, fractions of high, middle, and low categories are 93.8%, 1.7%, and 4.5%, respectively. Here, these fractions are calculated from the number of events categorized as high/middle/low, which are recorded in the housekeeping data even when the event data are not downloaded. On the other hand, in NXB observations, those are 53.7%, 1.5%, and 44.8%, respectively. The fraction of middle category is very small as expected. Since the low category contains only background-like events, it is expected that more events are classified to low category in the NXB observation than the Crab observation. More detailed screening is applied to the HXI data in the on-ground pipeline analysis. In the standard screening criteria for the scientific observations, time intervals at or around the SAA passages and those in Earth occultations are excluded. After the launch of the satellite, a new screening condition SAA2_HXI==0 is added to the standard screening to reduce the background in top-layer DSSDs, which is described in Sec. 5.
Since the HXI imager is a multilayer detector, an energy-dependent layer selection is applied for maximizing its sensitivity (signal-to-noise ratio). As shown in Fig. 4, the detection efficiency strongly depends on the layer as well as the incident x-ray energies. For example, the top-layer DSSD (layer 0) is capable of detecting only low-energy photons, typically below 30 keV, whereas the CdTe-DSD (layer 4) covers energies above to 30 keV. Although the best detection efficiency can be achieved using all the layers in all energy bands, the background level would inevitably be maximized. Thus, to maximize the sensitivity of the HXI, the energy-dependent layer selection is necessary. As explained in detail in Sec. 5, the NXB of the top-layer DSSD (layer 0) is different from DSSDs in the other layers due to the electron background, while those in the middle layers (layer 1 to 3) are very similar to each other. Thus, the HXI imager is separated into three groups, layer 0, layer 1 to 3, and layer 4, and their sensitivities are estimated using the in-orbit NXB spectra and the effective area. To optimize these sensitivities, only events satisfying conditions of for top-layer DSSD (layer 0), for the DSSDs in lower layers (layer 1 to 3), and for the CdTe-DSD (layer 4) are accepted. These correspond to for the top-layer DSSD, 12 to 80 keV for the DSSDs in lower layers, and for the CdTe-DSD.
The dead-time correction of the spectra and light curves of the HXI is performed by utilizing pseudoevents.23 The pseudoevents are the events triggered by the pseudotrigger, which is generated with a random time interval in the onboard FPGA in HXI Analog Electronics. Frequency of the pseudotrigger is set to be 2 Hz by default. Since the pseudoevents are treated in the same manner as normal events triggered by the ASICs, the number of pseudoevents passing through the data screening divided by the number of input pseudotriggers is a good estimate for the live-time fraction. The process of the dead-time correction is implemented in a dedicated tool hxisgddtime.
The average dead-time per event in the HXI is based on the Crab observation data. It is dominated by reset wait time (), which is a wait time to return to the state for waiting for the next trigger after the previous data acquisition in order to avoid triggering the noise induced by the AD conversion of the previous event. The other components of the dead time are the AD conversion ( to , depending on the pulse height), data transfer from ASICs to the FPGA ( to , depending on the number of signals above the DTHR for one trigger) and so on. In addition to the dead time accompanying each event, events in accidental coincidence with BGO active shields generate the dead time.
In the Crab observation, dead-time fractions estimated by the pseudoevents are 23.4% and 26.2% for HXI1 and HXI2, respectively. Their uncertainties are 1% to 2% point due to the inherent statistical uncertainties of the number of random pseudoevents. Then, we verified this number with independent estimation. Since the trigger rates of HXI1 and HXI2 are and , respectively, fractions of the dead time depending on the event rate are calculated as and 22.9%. Here, is the mean dead-time per event ( for HXI1 and for HXI2) and is the trigger rate. In addition to this, the fraction of the accidental coincidence with the BGO is estimated using the number of events classified to low or middle. By assuming that rates of low or middle events in the Earth occultations do not include the accidental coincidence events, but they are included in the low/middle events in the Crab observation, rates of the accidental coincidence events in the Crab observation are estimated by subtracting the low/middle rates in the Earth occultations (15.2 Hz for HXI1 and 15.5 Hz for HXI2) from those in the Crab observations (38.5 Hz for HXI1 and 38.1 Hz for HXI2). The accidental coincidence rates are calculated to be 23.3 Hz for HXI1 and 22.6 Hz for HXI2, meaning that 4.05% and 3.69% of the trigger rate are discarded, respectively. Thus, by considering the accidental coincidence events, the dead-time fraction is estimated to be 24.2% for HXI1 and 25.7% for HXI2, which match to those estimated by the pseudoevents within the statistical uncertainties.
Energy Response Matrix
Simulations of Detector Devices
The response matrix of the HXI is constructed by running Monte Carlo simulations since Compton scattering and secondary emissions are nonnegligible in hard x-ray bands. The simulation is composed of two steps: (1) Monte Carlo simulations for calculating the interactions of x-ray photons with detectors and passive materials in the HXI system and (2) calculations of charge transportation in the semiconductor detectors. This simulation code is based on an integrated response generator “ComptonSoft,”24 which is available in GitHub repository: https://github.com/odakahirokazu/ComptonSoft.
The Monte Carlo simulation part is based on the Geant4 toolkit library,25,26 which is widely used for the particle tracking in high-energy physics. Figure 5 shows a detailed Geant4 mass model of the HXI implemented in the detector response simulations. Since all the materials in the HXI system affect the detector response, most of the passive materials as well as the imager module and BGO active shields are included in the mass model. The Monte Carlo simulations are performed for each energy bin of the response matrix with monochromatic photons at the central energy of the bin. The incident photons for the simulations are generated in a horizontal plane with the detector size () located above the HXI-S entrance window, with an initial direction to the detector along the optical axis.
After the Monte Carlo simulations, the energy deposits obtained in the simulations are spread by considering the thermal diffusion. In this simulation, it is spread by the 2-D Gaussian with , where is the Boltzmann constant, is the temperature of the detector, and is the drift time from carrier creation to arrival to the electrodes. An additional factor is introduced to reproduce experimental results. Since it only considers the thermal diffusion, this factor represents an effect by the Coulomb repulsion.27,28
Using the spread energy deposits, the induced charge on each read-out strip is calculated with the simulation of the charge transportation. The induced charge is calculated utilizing a weighting potential derived from the Shockley–Ramo theorem29 as
The weighting potential is a solution of Laplace’s equation with the boundary condition of at the read-out electrode and at all the other electrodes. That for the strip detectors is calculated as
The response of the HXI DSSDs is more affected by charge loss due to the complicated electric field structure inside the detector rather than the charge trapping due to finite lifetime of the carriers described in Eq. (3). Since a mobility-lifetime product () of carriers in Si is 2 to 3 orders of magnitude larger than that of CdTe, it is assumed to be infinite in this simulation, which is equivalent to the charge correction efficiency of unity. On the other hand, a significant fraction of the charge is lost by local minimum of the electric potential due to a positive fixed charge on the surface at gaps among strip electrodes.9,30 It makes a subpeak at of the incident x-ray energy in the top side of DSSD, negative peak in the adjacent strip, and no signals in the bottom side of DSSD. Thus, we refer to these events as subpeak events. The subpeak events are unusable for the event reconstruction because they do not have any information of the position in the bottom side. Specifically, they are discarded in the consistency check between EPI values from both sides [Eq. (1) and Fig. 2]. Since the potential local minimum is located on the surface of the top side of the detectors, this effect reduces the detection efficiency at energies below . This effect is simply implemented as rectangular dead regions located at the strip gaps on the surface of the DSSD in our simulation. All induced charges corresponding to energy deposits in these dead regions are set to be zero.
All the parameters of the detector simulations for constructing the response matrix are listed in Table 4. The bias voltage, the mobility-lifetime products of holes and electrons, the diffusion factor, and the noise level of each strip are required for both CdTe-DSD and DSSDs, and one additional parameter for subpeak events is required for DSSDs. In addition to these, trigger efficiency, which reduces the detection efficiency at the lower energy end due to the energy resolution in the shaper for the trigger generation, and absorption by layers on surface of the DSSDs are multiplied to the energy response. The parameters used in the HXI response simulations are determined based on the ground calibration data of flight model and engineering model of the HXI. The mobility-lifetime product of CdTe-DSD is determined by fitting spectra of the single layer experiment of the engineering model detectors. The diffusion factor is set to be 2 as it reproduces the engineering model data. The noise parameter of each channel is estimated from the linewidths of an x-ray line at 59.5 keV from obtained in the ground calibration tests of the flight model detectors. Although only mean values of the noise levels from all the strips are listed in Table 4, the noise level is assigned strip by strip in the simulation.
Simulation parameters for the HXI detector response.
|Layer||0||1 to 3||4||0||1 to 3||4|
|Top-side noise (keV)||0.86||0.90||1.71||0.89||0.90||1.84|
|Bottom-side noise (keV)||2.49||3.04||1.75||2.57||2.92||1.80|
a(μτ)h and (μτ)e are the mobility-lifetime products of holes and electrons, respectively, dsubpeak is a thickness of the dead region where the subpeak events are generated, fdiff is the diffusion factor for spreading carrier clouds, dSiO2 is the thickness of the SiO2 layer on surface of the DSSDs, and Etrig and σtrig are the mean energy and σ of an error function for describing the trigger efficiency, respectively.
Detection efficiency at low energies of DSSDs is affected by three effects: subpeak events, trigger efficiency, and absorption by the inactive layer. Size of the dead region due to the subpeak events is assumed to be a rectangle with a width of , which is the same as the width of the gap among the strip electrodes of the DSSDs. By performing experiments using a single layer of the DSSD engineering model, thickness of the dead region is estimated to be from the energy dependence of the subpeak fraction.30 Trigger efficiency of the top-layer DSSD is assumed to follow the error function, and its mean energy and are estimated by measuring the detected count rate for 5.9-keV line from and its subpeak at 3.2 keV as a function of the trigger threshold in the ASIC. The trigger generation and sample/hold (and then the ADC) are performed in different analog-shaping chains in the ASIC with different shaping times of for the former part and for the latter. The noise level of the trigger is usually worse than that of the EPI values corresponding to the spectral resolution. In the other layers, trigger efficiency is not considered because the trigger threshold is much lower than the analysis threshold applied in the pipeline process. The thickness of the inactive layer, including layers and Al electrodes, on the surface of the DSSDs is estimated to be . Since the difference between photoabsorption cross sections with and Al is negligible, the inactive layer is treated as a layer with a thickness of . Thus, this value is set to the detectors in layers 1 to 3 as listed in the table.
In addition to the relatively rough estimates of the subpeak fraction, trigger efficiency, and the thickness, more detailed tuning of these parameters is needed to reproduce the in-flight data, especially around the lower energy end, in which the photon statistics is the highest in many cases. Among these three parameters, thickness is chosen as a free parameter for adjusting. First, 5 to 12 keV spectra of G21.5–0.9 observed by the SXI and HXI are simultaneously fitted, and the thicknesses of HXI1 and HXI2 are constrained to be 3.7 to and 2.5 to as 90% confidence intervals, respectively. Then, in order to constrain the thickness more tightly, the 5- to 40-keV Crab spectra of the HXI are fitted using response matrices with thickness of 4.0, 4.5, 5.0, and for HXI1 and 2.5, 3.0, 3.5, and for HXI2. As the result, for HXI1 and for HXI2 are found to give the best value. Thus, these values are used for constructing the HXI response matrix. This difference between HXI1 and HXI2 changes the detection efficiency at 5 keV by 7%. Note that this result does not mean that the actual thickness is different between HXI1 and HXI2. The difference would include all of the effect in the lower energies by the subpeak fraction, trigger efficiency, and the thickness.
Validation with Ground Calibration Data
For the purpose of testing the response simulation and its parameters, the simulated spectra and their properties are compared with those of the ground calibration data irradiated by x-rays from and . The simulations are performed in a geometry with a thermostat chamber and passive materials near the radioisotopes. All the x-ray and gamma-ray lines above 10 keV with emission probabilities larger than 0.1% are simulated. The Lund/LBNL Nuclear Data Search31 is referred to for the energies and emission probabilities of all the lines from and 26.3 and 59.5 keV from , while Lépy et al.32 is referred to for the other lines from .
The simulated spectra and the experimental spectra from the ground calibrations of the HXI flight model are shown in Fig. 6. Spectral shapes, including scattered components, tail structures due to the small in CdTe, and energy resolutions are well reproduced. Moreover, detection efficiency is consistent within , and at higher energies above , it matches better than 5%. Here, we should note that the self-absorption effect in the source is considered in this simulation by assuming a -thick Am as an absorber. Figure 7 shows fractions of the 2-strip events of the simulations and the experiments. These fractions have to match with each other for reproducing the spectral shape because energy resolutions of the 2-strip spectra are worse than those of the single-strip spectra. In other words, if the simulation overestimates the 2-strip fraction, the energy resolution of simulated spectra would be worse than those of the experiments. As shown in the upper panels, the 2-strip fractions are reproduced within 2% of the total event number in DSSD p-side and CdTe-DSD Al-side, whose EPI values are used for the spectral analysis. The discrepancy in CdTe-DSD Pt-side has almost no effects on the scientific analysis because the EPI values from CdTe-DSD Pt-side and DSSD n-side are only used for the consistency check between signals in top and bottom sides in the event reconstruction process.
Using the HXI response matrix described above, the Crab spectra are analyzed. All the standard processing and screenings (see Sec. 3) are applied to the data, and the spectra are extracted from circular regions with a radius of 4′. The background spectra are extracted from the blank sky observations (i.e., None2, IRU check out, and RXJ1856.5–3754), which contain the CXB as well as the NXB. The size of the extraction regions for the background spectra is the same as those for the Crab spectra. Net exposures of the Crab spectra after the dead-time correction are 5.92 ks for HXI1 and 6.14 ks for HXI2, respectively, and photons are detected by each of HXI1 and HXI2 in an energy range of 5 to 80 keV.
The Crab spectra observed by the HXI are shown in Fig. 8. They are fitted with an absorbed powerlaw model constant*tbabs*powerlaw convolved with the detector response and the telescope effective area using the spectral analysis software XSPEC.33 As shown in the lower panel in Fig. 8, the deviations between the Crab spectra and the best-fit model are at energies below . The telescope effective area is measured with uncertainties less than on ground34 and confirmed by the Crab observation.3 Although both the telescopes and detectors are well calibrated, residuals of level are seen above 50 keV. This might be due to calibration uncertainties in the telescope effective area or inappropriate modeling of the detector response.
The best-fit parameters and 90% confidence errors for the HXI Crab spectra are listed in Table 5. A difference of normalizations between HXI1 and HXI2, which is expressed by a constant parameter 0.968, is consistent with unity considering the uncertainty in the dead-time correction. For example, we also applied the independent dead-time correction, as discussed in Sec. 3, which are estimated to be 24.2% for HXI1 and 25.7% for HXI2. If we adopt these values, the powerlaw normalization at 1 keV is and a constant parameter is .
Best-fit parameters and 90% confidence errors of the Crab spectra.
|NH (cm−2)||Γ||Normalization at 1 keV (photons cm−2 s−1 keV−1)||Constant||F3−50 keV (10−8 erg cm−2 s−1)|
Our best-fit values of the powerlaw index and a normalization at 1 keV are consistent with the historical values of and proposed by Toor and Seward.35 Also, the spectral slope is consistent with the values of obtained by Suzaku/HXD PIN with HXD nominal position23 and obtained by the large off-axis observations of the Crab by NuSTAR.36 On the other hand, the normalization is not consistent with either of these observations. Our best-fit normalization is just between by the HXD and by NuSTAR. This result does not immediately mean that the HXI effective area is inconsistent with the other observatories because the Crab flux can vary on a yearly timescale.37
Spatial Dependence of the Detector Response
In addition to the spatially integrated detector response used in the Crab spectral analysis, we also verified the reproducibility of positional difference of the detector response. In the standard analysis tool, the response is separated into the detection efficiency and response matrix. The matrix defines the relation between the incident photon energy and the output EPI values, which correspond to the spectral shape. The spatial dependence of the detection efficiency is implemented for pixel by pixel based on the Monte Carlo simulations of the detector response. On the other hand, the response matrix is integrated over all the detector area in order to reduce data size of response database files. Therefore, the spatial dependence of the detector response is taken into account only by the detection efficiency.
To demonstrate the accuracy of the spatial dependence of the detector response, noisy strips located close to the center of the field of view (FoV) of HXI2 provide a good example. The analysis thresholds of these noisy strips are 5.58 and 6.71 keV, which are much higher than those in the typical strips, 3.66 keV. Due to the higher analysis thresholds, detection efficiencies at low energy in these strips are significantly smaller than the other strips, resulting in a dark line at the center of the image as shown in Fig. 9.
The spectra extracted from a rectangular region covering the noisy strips in HXI2 are shown in the upper panel of Fig. 10. There is a clear difference of the low-energy spectra between HXI1 and HXI2. The count rate at 5 keV in HXI2 is smaller than that of HXI1 by a factor of 1.5 because of the higher analysis thresholds in the noisy strips in HXI2. As shown in the lower panel in Fig. 10, this large difference between HXI1 and HXI2 is reduced to better than level by applying the detector response, which is generated in a standard manner by assuming a point source located at the red cross in Fig. 9. In this analysis, a broken powerlaw model is assumed by following the NuSTAR observations of the Crab pulsar38 because the selected spectra are strongly affected by it. The best-fit parameters of powerlaw slopes and , break energy , normalization , and the constant parameter are obtained. This result demonstrates the accuracy of the spatial dependence of the HXI detector response.
Properties of DSSD Background
Before the screening of SAA2_HXI==0 described in Sec. 3, NXB in top-layer DSSD (layer 0) of the HXI is dominated by a hard powerlaw component as shown in Fig. 11. In prelaunch estimations, NXB in top-layer DSSD is expected to show a similar level with those in the other layers of DSSDs because their background is thought to be mainly caused by albedo neutrons, interacting via elastic scattering.39 In this sense, the fact that the middle layers (layers 1 to 3) of DSSDs show a similar level is as expected.
The powerlaw component in top-layer DSSD is due to the low-energy albedo electrons for the following two reasons. First, this component extends up to 100 keV as shown in Fig. 11, but lower layers do not show this component strongly. The difference cannot be explained if it is caused by . It means that this background component originates from particles with low penetrating power. The second reason is the distribution of the background rate. Figure 12 shows a trigger rate of the HXI2 top-layer DSSD as a function of latitude and longitude of the satellite. It extends larger than the SAA and has a hot region at above North America. The hot region does not simply depend on the geomagnetic cutoff rigidity because it is not seen in other regions. This distribution is consistent with the electron distribution () observed by DEMETER/IDP in orbit.40
NuSTAR might observe this electron background as well. According to a web page on the background filtering,41 a “tentacle”-like region of higher activity is found to be located around longitude and above about latitude. This region is consistent with the distribution of the low-energy electrons as shown in Figs. 12 and 14.
The albedo electrons cannot directly come into the HXI imager because baffles made of Pb/Sn shields are implemented to block the stray light and CXB out of the FoV. Thus, top-layer DSSD probably suffers from the electrons scattering on the HXI-S entrance window (two layers of -thick poly-carbonate sheets) or the extensible optical bench. To reduce this background, we should have baffled the entrance window from the low-energy albedo electrons.
Since the electron background strongly depends on the satellite position as shown in Fig. 12, a selection with the satellite position successfully reduces top-layer DSSD background down to 10% to 20% level as shown in Fig. 13. The definition of the regions discarded in this selection is shown in red lines in Fig. 14. This region is implemented in the standard screening procedure as SAA2_HXI. This selection reduces observation efficiency as well. Area of the selected regions (which excludes the SAA as well) is of the total area where the satellite orbits above, while those of Suzaku SAA definition is .
The electron background depends on the satellite position, even after excluding the high background regions using the satellite position selection. Figure 14 shows a distribution of count rate of HXI2 top-layer DSSD with corresponding to , which should have only background events. Clearly, the background level in top-layer DSSD is lower at and , while it is higher at . To investigate the position dependence of the top-layer DSSD background, spectra are extracted from six regions, which are divided at and . In Fig. 15, spectra from regions with highest and lowest background levels are plotted. Although there is a large difference at higher energies, it does not affect the scientific analysis because signals above 30 keV in top-layer DSSDs are discarded in the energy-dependent layer selection. At energies below 30 keV, the background level can change by a factor of 3 at maximum depending on the orbital phase.
Properties of CdTe-DSD Background
NXB in the CdTe-DSD (layer 4) of HXI is composed of many activation lines from radioactive isotopes induced by geomagnetically trapped protons in the SAA as shown in the cyan histogram in Fig. 13. Since orbit inclination angle of Hitomi is 31 deg, the HXI passes the SAA for 8 to 9 times a day. Low-energy protons trapped in the SAA generate radioactive isotopes inside the detectors via interactions between protons and heavy atoms, such as Cd, Te, and Bi, contained in the CdTe-DSDs and BGO shields. Gamma-ray photons and particles from these radioactive isotopes are the main cause of the CdTe-DSD background.
Since the radioactive isotopes are generated in the SAA passage, the CdTe-DSD background depends on the time after the SAA passage, which is defined as T_SAA. Figure 16 shows spectra sorted by T_SAA. In this figure, the SAA pass is defined as , where 5000 s roughly corresponds to one orbital period of Hitomi, and non-SAA pass is defined as . It is clear that a few lines at and are rapidly decaying after the SAA passage, but spectra at 30 to 80 keV, which are used for scientific analysis, show less variability .
The properties of activation background in the CdTe-DSD observed by the HXI give us essential information to understand the activation background and significantly improved the accuracy of the simulations of the activation background. Details of the simulation studies are described by Odaka.42 This result will be a great help for future hard x-ray missions, such as FORCE.43
Final Background Spectra
The final spectrum of the NXB and blank sky after all the processing and screening is shown Fig. 17. Compared with the NXB spectra before the energy-dependent layer selection (Fig. 13), the total background level is clearly reduced by ignoring the strong line at 20 to 30 keV in CdTe-DSD and the albedo electron component dominating the higher energy region in top-layer DSSD. The NXB level is as low as the preflight requirement of . The photon detection efficiency of the top-layer DSSD above 30 keV is ignorable, while that of the four layers of DSSDs (in total 2-mm thick) below 30 keV is . Therefore, the energy-dependent layer selection efficiently reduces the background with a small loss of detection efficiency. It demonstrates the effectiveness of the design of the stacked semiconductor detector for achieving better sensitivity.
Figure 17 shows one additional demonstration of the HXI performance. A clear spectral difference between the blank sky and NXB spectra can be seen below . It indicates that the HXI is able to detect the CXB below . Indeed, as shown in Fig. 18, the CXB spectra are significantly detected by the HXI. In this figure, a powerlaw model with a photon index of is overplotted as a historically measured spectral model of the CXB.44 Spectral fitting by this model with a fixed photon index provides the best-fit powerlaw normalization at 1 keV of . Here, the spectral fitting is performed in 5 to 10 keV because the 15- to 30-keV spectra deviate from the powerlaw model due to the variability of the albedo electron background as previously shown in Fig. 15. If the accuracy of background modeling is improved, it would be possible to detect the CXB spectra above 30 keV. A ratio of the best-fit normalization with the value reported by De Luca and Molendi44 is , which is roughly consistent with the CXB fluctuation, based on a relation of 45 with the HXI FoV and an assumed upper cutoff flux . Thus, this result shows that the HXI has a good sensitivity for extended sources, which enables detection of the CXB below 30 keV.
The HXI showed good performances and provided us important insights on the NXB in the hard x-ray band although it was lost after only two weeks of observations. The Crab spectra are well reproduced by the detector response constructed on the ground calibration data. The residual between the Crab spectra and the best-fit absorbed powerlaw model is less than at energies below 50 keV. The best-fit spectral parameters of the Crab are consistent with the historically reported values. This result indicates the correctness of the telescope effective area of the HXT and the detector response of the HXI. The NXB in top-layer DSSD is found to be dominated by the background due to low-energy albedo electrons. Utilizing its strong dependence on the latitude and longitude, it can be reduced to 10% to 20%. Even after this selection, the electron background in 20- to 30-keV varies by a factor of 3, depending on the orbital phase. The activation background in the CdTe-DSDs above 30 keV is more stable within . The final spectrum of the NXB after all the processing and screening satisfies the preflight requirement level of , and it enables to detect the CXB. The properties of the in-orbit background of the HXI would be useful for the future hard x-ray missions.
We acknowledge all the Hitomi team members, including many graduate students, for their great contributions to the HXI and the Hitomi project. We acknowledge support from the JSPS/MEXT KAKENHI Grant Nos. 24105007, 15H03639, 25287059, and 24244014 and the JSPS Core-to-Core Program. All U.S. members acknowledge support through the NASA Science Mission Directorate. Stanford and SLAC members acknowledge support via DoE contract to SLAC National Accelerator Laboratory DE-AC3-76SF00515 and NASA Grant No. NNX15AM19G. French members acknowledge support from the Centre National d’Etudes Spatiales.
Kouichi Hagino is an assistant professor at Tokyo University of Science. He received his BS and MS degrees in physics from the University of Tokyo in 2010 and 2012, respectively, and his PhD in physics from the University of Tokyo in 2015. He has been working on development of semiconductor detectors for high-energy astrophysics.