1 Introduction
The X-rays emitted by galaxy clusters, black holes and neutron stars are important for studying the extraordinary gravitational [1], electromagnetic, and nuclear-physics environments inside neutron stars, and exploring the exotic states of matter, of which the density and pressure are higher than in atomic nuclei [2]. X-rays of pulsars are attractive in autonomous navigation due to their properties of stable, periodic and predictable signatures [3, 4].
In recent years, several projects of celestial soft X-ray detection in space have been carried out or planned, such as the X-ray Imaging Spectrometer of Suzaku [5], Neutron Star Interior Composition Explorer (NICER) [2] and Large Observatory for X-ray Timing (LOFT) [6]. The Suzaku satellite used charge-coupled devices to observe celestial X-ray sources of 0.2–12 keV. The NICER, scheduled to launch in early 2017, adopts silicon drift detector (SDD) for X-ray timing and spectroscopy of neutron stars in the soft X-ray band. The LOFT, proposed to launch around 2025, will adopt SDD for exploiting and diagnostics of rapid X-ray flux and spectra in 2–50 keV from neutron stars and black holes.
Being able to work at close to room temperature with a Peltier cooler, SDD is considered as an appropriate choice for soft X-ray detection, with high energy resolution and count rate. In SDD, an incoming photon generates a number of electrons and holes. The holes drift towards the back side of the detector, while the electrons drift towards the anode electrode with a proper bias of a set of cathodes [7]. The accumulated charge at the anode is connected to the gate of a field effect transistor, which forms the first stage of a charge sensitive amplifier.
A prototype high resolution soft X-ray spectrometer is under development at USTC, aimed at conducting astrophysics studies, and pushing forward X-ray based navigation and X-ray communication in space [8]. The detector, SDD H30 [9] made by KETEK GmbH, Germany, has a maximum input count rate of 1000 kcps, and excellent energy resolution. With an active area of 30 mm2, it is of about 80% at 1 keV and 98% at 10 keV in efficiency. In this paper, the readout electronics and test results are described in detail.
2 The readout electronics
2.1. Architecture
High-resolution time and charge measurement are widely used in nuclear physics [10-12]. Fig. 1 shows a block diagram of the soft X-ray detection (SXD) readout electronics for time and energy measurement. It is mainly composed of two sections: analog and a digital electronics. In the analog electronics section, a pre-amplifier module (PA), a high voltage module(HV) and a temperature control module (TC) are specially designed by the Original Equipment Manufacturer of SDD. All the other modules after the PA are developed by our group. The PA output goes to two channels: the one marked as 'Q_out’ from a slow shaper is for energy measurement, and the other marked as 'T_out’ from a fast shaper is for time measurement.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F001.jpg)
In the digital electronics section, the 'Q_out’ is fed into an analog peak hold module, and then sampled by an analog-to-digital-converter (ADC), while the 'T_out’ is converted into a hit signal by a constant fraction timing circuit (CFT), and then sent to a time-to-digital-converter (TDC) for digitization. The ADC and TDC signals are buffered and fed into a data processor (DP). The system can be monitored and controlled by the DP according to the standard Controller Area Network (CAN) protocol. The external power supplies of the whole electronics are ±6 V, +5 V and +12 V.
2.2. Pre-amplifier
As shown in Fig. 2, a matched PA [13] for SDD is used to amplify the output current pulse of the detector, which offers an ultra-low noise, and a ramped reset type CSA with high gain (4.5 mV/keV ±15%). The amplitude of the ramped output signal is between ± 2 V with a reset duration of <5 µs, and the period of the reset ranges from a few milliseconds to a few hundred milliseconds, which depends on the frequency of the incoming photons and the temperature of the detector.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F002.jpg)
The output of PA is a voltage step when the signal produced by an X-ray photon appears. The rise time (from 10 ns to 100 ns [14]) mainly depends on the location where X-rays interact with the SDD, while the voltage amplitude depends on the energy deposition of the incoming photon.
2. 3.Energy measurement circuit
For energy measurement, the signal amplitude is detected efficiently, achieved via the analog peak detection and ADC method. As shown in Fig. 3, the PA is followed by a slow shaping circuit consisting of a CR-RC2 filter, which is employed to form a quasi-Gaussian signal. The peak value of the quasi-Gaussian signal is held by a peak hold module (PH300 from Amp-Tek) and then digitalized by a 14-bit ADC (AD9243) with 3 MHz sampling clock. The pole zero cancellation circuit is not necessary in case that being the signal form PA a step voltage [15].
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F003.jpg)
Noise contribution of the SDD and readout electronics is analyzed, to find the performance-limiting factor of the system.
Energy resolution of the detector, staged in the datasheet, is 126.9 eV @ 5.9 keV [16] (at below −60°C using a Peltier cooler, peaking time of 16 µs). For the slow shaping circuit, a Pspice simulation shows that the noise RMS is 3.4 mV, which is equivalent to 44 eV @ 5.9 keV. A 14-bit ADC with a quantization error of 0.12 mV (1.5 eV @ 5.9 keV) is selected, which shall be a minor role in affecting the total energy resolution.
In addition, the rise time of the real signal will cause extra loss of the shaped signal height, which is termed as the Ballistic Deficit. Although increasing the peaking time (i.e. the time constant of CR and RC) can decrease the Ballistic Deficit, it will increase the dead time of the energy measurement. A peaking time of 2 µs is determined finally to guarantee a proper counting rate and the energy resolution. Test results with a signal generator show that the actual electronic noise (FWHMelectronic) of the energy measurement circuit is equivalent to 68 eV @ 5.9 keV as a constant in the dynamic range with a peaking time of 2 µs.
2.4. Time measurement circuit
Compared with leading edge timing, a constant fraction timing circuit essentially eliminates the amplitude-dependent time walk for signals and has better timing resolution [17]. Fig. 4 explains the CFT principle briefly, and a delay time '∆t’ of 8 ns as an attenuation coefficient 'q’ of 0.5 is chosen, leading to a fixed delay 't’ of 16 ns.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F004.jpg)
As shown in Fig. 5, the PA is followed by a fast shaping circuit, which consists of an AC coupled circuit (CR with a time constant of 0.5 µs) and an amplifier, to condition the signals for time measurement. The fast shaper also drives three branches of the CFT circuits: the first branch is delayed by a time of 'td’ (8 ns), the second one is attenuated by a coefficient of 'q’ (0.5) and the third one is compared with an adjustable threshold for noise canceling. The output of both comparators is fed into an 'AND’ gate for coincidence and then digitized by an FPGA based TDC.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F005.jpg)
Since accuracy of time measurement is highly sensitive to CFT circuit, stringent hardware design is considered. A 100 Ω impedance fixed delay line (1919-20B, Data Delay Device) is used to achieve high precision delay and impedance matching. An ultra-high-speed comparator (MAX9601, MAXIM Inc.) is selected, with 30 ps propagation delay dispersion and differential PECL outputs. A programmable delay line (NB6L295 from ON semiconductor) configured at 5 ns, is used to delay the CFT branch to guarantee the timing constrain. The AND gate (MC100LVEL05,ON Semiconductor) has two differential ECL inputs with a propagation delay of 340 ps.
2.5. TDC integrated in FPGA
To digitize the output signals from the CFT module, a TDC is implemented in a XILINX FPGA (XC5VLX110), which reduces the system complexity and provides good flexibility [18, 19]. Fig. 6 shows its block diagram. The arrival time of the 'hit in’ signal is digitized by a 3-bit 'fine time’ and a 29-bit 'coarse counter’ modules, and then encoded to a final 32-bit time measurement result, which is finally read out by a FIFO.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F006.jpg)
A 50 MHz external clock is fed into the FPGA, and then used by an internal phase lock loop (PLL) to generate four 100 MHz clock outputs with 90° phase interval (0°, 90°, 180°, 270°), which is equivalent to 2.5 ns least significant bit (LSB) and much less than the transit time spread of SDD.
3 Test results
The test platform to evaluate performance of the system mainly consists of a power supply (KEITHLEY 2230-30-1), a signal source (Tektronix AFG3252), an oscilloscope (Tektronix DPO5104), a 55Fe X-ray source, and the aforementioned analog electronics module and digital electronics module. A series of tests are conducted to evaluate performance of the SXD readout electronics.
3.1. Energy measurement results
On the SXD system, using the 55Fe source, the energy resolution is measured at 153.4 eV @ 5.9 keV FWHM (Fig. 7a). Statistical distribution of time intervals of two photons is shown in Fig. 7(b), where the counts have the characteristic of Poisson distribution and behave in accordance with the theoretical property of radioisotope.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F007.jpg)
The energy resolution and linearity (properties of the detector and electronics) X-rays of different energies are listed in Table 1. The X-rays are electron beam-induced fluorescence from samples of Cu, Fe, Cr and Ti [20]. The linearity of one representative SDD indicates that the integral non-linearity (INL) of less than 0.1% can be achieved. It is worth noting that the minor INL can be corrected by off-line analysis.
X-ray | Cu-Kβ | Cu-Kα | Fe-Kβ | Fe-Kα | Cr-Kβ | Cr-Kα | Ti-Kβ | Ti-Kα |
---|---|---|---|---|---|---|---|---|
Energy (keV) | 8.91 | 8.04 | 7.06 | 6.40 | 5.95 | 5.41 | 4.93 | 4.51 |
Meana (code) | 5583 | 5072 | 4479 | 4091 | 3821 | 3503 | 3212 | 2960 |
INL (%) | −0.05 | 0.07 | 0.06 | 0.03 | 0.01 | 0.07 | −0.02 | −0.04 |
FWHMall (eV) | 163.2 | 162.1 | 149.1 | 148.0 | 138.8 | 137.6 | 127.7 | 126.8 |
FWHMFanob (eV) | 148.4 | 147.2 | 132.7 | 131.5 | 121.0 | 119.6 | 108.1 | 107.0 |
3.2. Time measurement results
The RMS of time measurement is basically contributed by two factors: the noise of circuits and the quantization error. A 10 kHz pulse (30 ns rise time, 50 mV amplitude, generated by a Tektronix AFG3252) is utilized to drive the fast shaper module directly, to evaluate the performance of time measurement. Test results of the SXD (including contribution of both analog and digital electronics) are shown in Fig. 8. It can be seen that the standard deviation (STD) is 3.7 ns, which is much less than the transit time spread of SDD.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F008.jpg)
As mentioned in Section 3.3, CFT essentially eliminates the amplitude-dependent time walk for signals. Fig. 9(a) shows the tests to observe performance of the CFT circuit. The first channel (CH1) is the reference signal with a frequency of 10 kHz and 2.5 ns rise time; while CH2 is created based on waveform of the fast shaper, which has different amplitude (0.4–4 V) with 10 ns rise time. Fig. 9(b) shows the performance of CFT circuit. The propagation delay of CFT circuit is 20 ns with the time walk less than 3 ns.
-201701/1001-8042-28-01-014/alternativeImage/1001-8042-28-01-014-F009.jpg)
4 Conclusion
The readout electronics system for a prototype soft X-ray spectrometer with SDD is successfully implemented, and good time and energy measurement performance are achieved. A time measurement resolution below 5 ns is obtained. The energy resolution is 153.4 eV FWHM @ 5.9 keV with analog peak detection and ADC method. XRF tests indicate that the INL is less than 0.1% with a dynamic range of 1 keV to 10 keV. All these show that the SXD system has broad potential for X-ray timing and observation missions in space.
Gravitational waves from the axial oscillation of neutron star considering non-Newtonian gravity
. Nucl Sci Tech, 24 (2013) 050508. doi: 10.13538/j.1001-8042/nst.2013.05.008The Neutron star Interior Composition ExploreR (NICER): an Explorer mission of opportunity for soft X-ray timing spectroscopy
. SPIE, Space Telescopes and Instrumentation 2012. doi: 10.1117/12.926396Recursive estimation of spacecraft position using X-ray pulsar time of arrival measurements
.Spacecraft Navigation Using X-ray Pulsars
. NRL Review. 2006.The X-ray Observatory Suzaku
. Astronomical Society of Japan, 2007. doi: 10.1093/pasj/59.sp1.S1LOFT-the Large Observatory For X-ray Timing
.X-ray Silicon Drift Detector-CMOS Front-End System with High Energy Resolution at Room Temperature
. IEEE T Nucl Sci. doi: 10.1109/TNS.2015.2513602The Role of X-rays in Future Space Navigation and Communication
.A 16-Channel high-resolution time and charge measurement module for the external target experiment in the CSR of HIRFL
. Nucl Sci Tech, 25, 010401(2014). doi: 10.13538/j.1001-8042/nst.25.010401A digitalizing board for the prototype array of LHAASO WCDA
. Nucl Sci Tech, 22(2011) 178-148. doi: 10.13538/j.1001-8042/nst.22.178-184Design and Implementation of a High Resolution DAQ System for an (e, 2e+ion) Electron Momentum Spectrometer
. IEEE T Nucl Sci, 2015, Vol. 62, No. 6. doi: 10.1109/TNS.2015.2497282Characterization of the silicon drift detector for NICER instrument
.Study of recursive model for pole-zero cancellation circuit
. Nucl Sci Tech, 25, 010403(2014). doi: 10.13538/j.1001-8042/nst.25.010403A high precision time-to-digital converter based on multi-phase clock implemented within FPGA
. Nucl Sci Tech, 21(2010)123-128. doi: 10.13538/j.1001-8042/nst.21.123-128A field-programmable-gate-array based time digitizer for the time-of-flight mass spectrometry
. Rev Sci Instrum, 85, 045115(2014). doi: 10.1063/1.4870922