Simulation and photoelectron track reconstruction of Soft X-ray Polarimeter

Simulation and photoelectron track reconstruction of Soft X-ray Polarimeter 1 corresp first-author fogsn2012@foxmail.com fogsn2012@foxmail.com fogsn2012@foxmail.com https://orcid.org/0000-0001-5622-5864 1 corresp liuhb@gxu.edu.cn liuhb@gxu.edu.cn liuhb@gxu.edu.cn 1 1 1 1 1 1 2 1 3 corresp liuqian@ucas.ac.cn liuqian@ucas.ac.cn liuqian@ucas.ac.cn https://orcid.org/0000-0003-4658-6361 4 https://orcid.org/0000-0001-6046-6935 1 corresp lew@gxu.edu.cn lew@gxu.edu.cn lew@gxu.edu.cn Guangxi Key Laboratory for Relativistic Astrophysics, School of Physical Science and Technology, Guangxi University, Nanning 530004, China Guangxi Key Laboratory for Relativistic Astrophysics, School of Physical Science and Technology, Guangxi University, Nanning 530004, China PLAC, Key Laboratory of Quark & Lepton Physics (MOE), Central China Normal University, Wuhan 430079, China PLAC, Key Laboratory of Quark & Lepton Physics (MOE), Central China Normal University, Wuhan 430079, China School of Physical Science, University of Chinese Academy of Sciences, Beijing 100049, China School of Physical Science, University of Chinese Academy of Sciences, Beijing 100049, China Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China The soft X-ray polarimeter (SXP) is a detector with a wide energy range, large area, and large field of view. A SXP will be mounted on the Chinese Space Station and will mainly focus on detecting the polarization of transient soft X-ray (210) sources, especially gamma-ray bursts (GRBs). In this work, a polarimeter detector unit is taken as an example, and Geant4 and Garfield++ software are used to simulate the detection efficiency and track production. An improved track reconstruction algorithm is proposed and used to reconstruct two-dimensional images of the tracks. In this method, the initial emission angle of photoelectrons is reconstructed from the initial part of the track by shortening or extending the initial part of the track until the remaining track is straight, and the number of pixels is within an adjustable threshold. The modulation factor of the photoelectronic tracks after reconstruction reaches approximately 57% in the photon energy range of 7 keV to 10 keV. X-ray Polarimetry Track reconstruction algorithm 1 Introduction 1 X-ray polarization offers a new window in high-energy astrophysics. Sensitive X-ray polarimetry is a promising method for revealing crucial and unique information about the structures and physical processes of black holes, neutron stars, and all classes of X-ray sources in the universe. The measurement of X-ray polarization was first successfully achieved in an X-ray polarimeter aboard the OSO-8 satellite for Bragg peak measurement [ 1 1 ] although only two separate energies (2.6 keV and 5.2 keV)were measured. In 2001, Costa used a high-resolution gas detector to measure X-ray polarization in the laboratory[ 2 2 ]. Since then, X-ray polarization measurements have become more widespread. Various particle physics and astrophysics missions for detecting the X-ray polarization, including IXPE [ 3 3 ] and eXTP [ 4 4 ], are being prepared. Other missions, such as XIPE [ 5 5 ]and PRAXyS [ 6 6 ], have been planned but subsequently canceled. Recently, CubeSat performed new X-ray polarization measurements of the Crab Nebula[ 7 7 ] using gas detector technology. In these missions, the polarization of energetic astrophysical sources was mainly measured through a narrow field of view using a high-sensitivity X-ray polarization detector, which cannot respond rapidly to transient sources such as gamma-ray bursts (GRBs). Based on the radiation mechanism, it is known that the magnetic field configuration of the GRB is reflected in the radiation polarization. With this in mind, it is hoped that the magnetic field configuration of the jet can be studied by measuring the polarization of the GRB radiation to understand this extreme explosion phenomenon[ 8 8 ]. The SXP is currently in the conceptual design stage in which a detector with a large area and a field of view is being developed. The detector design directly affects the detector sensitivity and needs to be optimized through simulations. The obtained data must be processed by a track reconstruction algorithm to obtain the polarization of the X-ray source. A good track-reconstruction algorithm can improve the accuracy of the results. Algorithms for track reconstruction are also constantly being developed. Bellazzini et al. first proposed an algorithm for track reconstruction[ 9 9 , 10 10 ]. Analysis of the track moments allows some properties of their distribution, such as the barycenter, the major axis, and most importantly, the projection on the pixel plane of the point of absorption of the X-ray photons, to be analyzed. Differing from the traditional moment analysis method based on the track image as well as an algorithm based on the shortest path problem in graph theory [ 11 11 ] and adaptive cut methods [ 12 12 ], the algorithm proposed here is based on the detector condition and is called the iteration cut method. A brief introduction to SXP and the structure of the SXP detector unit is given in Sect. 2. The SXP Monte Carlo simulations are presented in Sect. 3. The algorithm for track reconstruction is proposed in Sect. 4, and the conclusions are presented in Sect. 5. 2 Detector structure 1 The SXP is composed of multiple detector arrays. They give an effective field of view of approximately 90° and a detection area of 300 cm 2 for detecting the polarization of transient X-ray transient sources. The SXP contains 96 detection modules. We use a single detector module to illustrate the algorithm in this study. A single detector module is shown in Fig. 1 Fig. 1 . The SXP will be mounted on the Chinese space station to detect X-ray polarization within the energy range of 2 keV to 10 keV. The photoelectric effect has the largest reaction cross-section in the soft X-ray energy range, and results in the generation of photoelectrons with lower energy. Sufficiently distinct tracks can be obtained using thick-gas electron multiplier(THGEM) [ 13 13 , 14 14 ] detection. THGEM has the advantages of high detection efficiency, simplicity, and stable operation. Track detection is the key to extracting the polarization of the X-ray source. The Topmetal-2 - chip is chosen for charge collection and reading out the electronic signals [ 15 15 ]. A field-programmable gate array (FPGA) is used to store information. The internal structure of the detector is divided into three parts, namely, the drift region, multiplier region, and induction region. As shown in Fig. 2 Fig. 2 , a photon incident in a direction perpendicular to the drift region of the detector is converted into a photoelectron through the photoelectric effect in the working gas. When the photoelectron passes through the working gas, it ionizes the gas, giving rise to electron-ion pairs. The electric field causes the electrons to drift toward the multiplier region, where they are amplified through avalanche generation. The photoelectron track is finally detected via two-dimensional imaging. The initial direction of the photoelectron is related to the polarization of the incident photons. The ejection probability is the highest when the initial direction has the same angle as the polarization. The polarization probability is given by the distribution of cos 2 θ [ 16 16 ], where θ is the polar angle with respect to the incident X-ray direction. The THGEMs are fabricated via standard PCB mechanical drilling and global etching processes, which allow for large-scale economic production of large-size detectors [ 17 17 ]. When two different negative voltages are applied to the extremities of the metal layer, electrons can drift through the etched hole and avalanche multiplication can occur, thus amplifying the charge for detection. The THGEM has several advantages, including discharge prevention, sub-millimeter position resolution, large gain, and ease of manufacture and maintenance. A thin THGEM [ 18 18 - 20 20 18-20 ] with 200 μm thickness, 200 μm hole diameter, 500 μm hole pitch, 30 μm rim width, and 10 μm thickness was produced. The Topmetal-2 - is a direct charge sensor that collects drift charges directly without post-processing and reads them out into a 2D projection image. The pixel array is patterned in a matrix of 72×72 pixels with an 83 μm pitch, resulting in a charge-sensitive area of 6×6 mm 2 . Each pixel contains a low-noise charge-sensitive preamplifier (CSA) to establish an analog signal and a discriminator with a tunable threshold to generate hits. The Topmetal-2 - is characterized by its small volume with the overall geometric size of 8 mm×9 mm×0.7 mm. The average equivalent noise charge measured is less than 15 e - . The low noise allows for improved trail detection. Furthermore, the sensor is resistant to radiation and remains functional even after 28 h of irradiation with a 60 Co source at the dose rate of 3.6 . 3 Monte Carlo simulations 1 The expected track was simulated using the Geant4 [ 21 21 ], Garfield++ [ 22 22 ], and COMSOL Multiphysics [ 23 23 ] software packages. Geant4 is a toolkit for simulating particle-matter interactions. The detector model was constructed in Geant4 and used to simulate the initial track. The structure of the SXP is shown in Fig. 3 Fig. 3 . Here, photons enter the detector through a collimated aperture and are detected by the photoelectric effect in the gas. The effect of the THGEM is modeled in the following simulations. The main criteria for selecting the gas are high detection efficiency and long track length. High detection efficiency is beneficial for the detection of weak sources, whereas long tracks are useful for reconstruction. As shown in Fig. 4 Fig. 4 , the detection efficiency and track length of different gases with the same thickness are simulated (panel(a) and panel (b)). It can be seen that the cross section of the photoelectric effect decreases with increasing energy of the incident photons, thus reducing the detection efficiency. The efficiency at 2 keV is lower than that at 3 keV because of the Be window, which acts as a shield against low-energy photons. Another factor worth considering is the gas diffusion coefficient. The diffusion coefficient of the gas should be as small as possible. Through experiments, it was found that the diffusion coefficients of the candidate gases were all within a certain range which did not affect the reconstruction algorithm. For the above reasons, a mix of Neon (80%) and dimethyl ether (DME 20%) was chosen as the detection gas. Ne-based gas mixtures have longer photoelectron tracks in the lower energy bands while retaining reasonable efficiency, while the DME absorbs UV photons produced by charge propagation in the gas with small transversal diffusion and relatively high gain [ 24 24 ]. After completing the track simulation in Geant4, Garfield++ and COMSOL Multiphysics were used to perform the track simulations in more detail. Garfield++ is a particle detector tool that is useful for detailed simulations in which a gas is used as a sensitive medium, while COMSOL Multiphysics is a multi-physical field simulation software based on finite element calculations. The gas components and the THGEM model of the detector were input into COMSOL Multiphysics to calculate the electric field distribution. The potential distribution calculated via COMSOL Multiphysics is shown in Fig. 5 Fig. 5 . The low-energy physics model in Geant4 was used to generate the secondary electrons and then reconstruct them into delta electrons in Garfield++ for simulations [ 25 25 ]. The Topmetal-2 - chip senses the charge signal through the exposed copper on its surface; it not only generates a signal at the incident location of each electron but also in the area before impact. To simulate the actual electronic signal collection, the charge value at each pixel position in the Garfield++ simulation was redistributed between the pixel and the surrounding pixels with a Gaussian distribution, and the appropriate σ value was selected for comparison with the experiment to obtain more accurate track simulations. Using Garfield++ simulations, which are intrinsically slow, a Gaussian distribution for the diffusion can be added to the slender electron track information produced by Geant4 to make the tracks thicker to simulate the signal sensed by the electrons on the Topmetal-2 - chip. The width of the Garfield++ track is first fitted with a Gaussian distribution to obtain the corresponding σ value. The fitted Gaussian distribution is then added to the Geant4 track to obtain a track width similar to that simulated with Garfield++. At the same time, different apertures in the THGEM directly affect the transverse diffusion of electrons, leading to different track thicknesses. The effect of different pore sizes on the track can be easily simulated by varying the σ value of the Gaussian distribution. Fig. 2 Fig. 2 shows an example of a track obtained through simulations. The size of each histogram bin is equivalent to the size of each pixel in the Topmetal-2 - . Panel (a) shows the track used in the reconstruction. The size of the square represents the charge value at the coordinate position. The label numbers at the horizontal and vertical axes represent the detector readout positions, which coincide with the Topmetal-2 - pixels. 4 Track reconstruction algorithm 1 4.1 Event selection 2 Several types of tracks can be difficult to reconstruct even under ideal simulation conditions. These include certain completely symmetrical tracks and certain elliptic or circular tracks. The experimental conditions are more complicated than the simulated ones. The presence of both electronic noise and multiple tracks in a frame may affect the reconstruction results. The isolated noise points can be removed through image erosion and dilation algorithm processing of the tracks. To estimate the track range, the coordinate values of the leftmost point ( x min ), lowest point ( y min ), rightmost point ( x max ), and highest point ( y max ) need to be considered. When the conditions x max - x min < 6 and y max - y min < 6 are satisfied, the track range is limited to six pixels, which may correspond to circular or elliptical tracks. By replacing the longer and shorter axes of x and y by L and S , respectively, a track that satisfies ( L max - L min )/( S max - S min ) < 1.5 is considered as an elliptical track and eliminated. The thresholds of 1.5 and 6 used here can be adjusted in the future. In addition, this process removes several difficult parts in the reconstruction process of symmetric charge distributions, such as finding the Bragg peak track, as explained in detail below. These tracks should be screened out in advance, while the remaining tracks are reconstructed in the following steps. 4.2 Event reconstruction 2 Electrons may be scattered multiple times as they pass through the gas, thus deviating from the initial ejection direction. When electrons are close to stopping, they release most of their energy in the form of a Bragg peak. Thus, only a modest charge is accumulated at the initial interaction point. In this case, monitoring only the start of the track is sufficient to achieve full reconstruction. The repeated iteration method is divided into the following four main steps, of which the first two are similar to those in the traditional method. 1) The barycenter and gravity line are calculated from the location and charge (as the weight) of the track in Fig. 2(a) Fig. 2(a) . This line roughly marks the direction of the track trend. 2) A line represented by the green line in Fig. 2(a) Fig. 2(a) is traced at the barycenter perpendicular to the gravity line. Because Bragg peaks are accumulated at the tail of the track, the head and tail can be separated using third-order moment calculations to determine the half-zone where the photoelectron interaction point is located. The third-order moment M d is calculated as <math display="block" id="M1"><mrow><msub><mi>M</mi><mtext>d</mtext></msub><mo>=</mo><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mtext>pixels</mtext></mrow></msub></mstyle><mo stretchy="false">(</mo><msup><mi>d</mi><mn>3</mn></msup><mo>×</mo><mi>q</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where d is the distance between the pixel and the line perpendicular to the gravity line ( d may be negative), q is the charge of the pixel, and ∑ pixels indicates summation over all the pixels. Thus, by comparing the d and M d of each pixel with similar or dissimilar signs, it can be determined whether the pixel is within the region where the photoelectron Bragg peak is located, or in the region where the photoelectrons are produced, which is indicated by the semicircle in the figure. Tracks with M_d = 0 are symmetrical in shape and difficult to reconstruct; thus, such tracks are rejected. 3) Multiple lines perpendicular to the gravity line are traced to re-divide the track, as shown by the multiple green lines in Fig. 2(a) Fig. 2(a) . The length of the track determines the direction of the subdivision, and for excessively long tracks, the track is subdivided along the direction that shortens the starting half of the track. The track length is calculated as <math id="M2"><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup><mo>/</mo><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>max</mtext></mrow></msubsup></mrow></math> . Each point on the track has the distance of <math id="M3"><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></math> from the gravity line and the distance of <math id="M4"><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>max</mtext></mrow></msubsup></mrow></math> from the line perpendicular to the gravity line. The sum of the squares of <math id="M5"><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></math> and <math id="M6"><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>max</mtext></mrow></msubsup></mrow></math> represents the square of the distance of the point from the barycenter M 2 . M 2 is given by <math display="block" id="M7"><mrow><msub><mi>M</mi><mn>2</mn></msub><mo>=</mo><mstyle displaystyle="true"><msubsup><mo>∑</mo><mi>n</mi><mrow><mi>n</mi><mo>=</mo><mtext>pixels</mtext></mrow></msubsup></mstyle><msup><mrow><mo stretchy="false">(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>−</mo><msub><mi>x</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>−</mo><msub><mi>y</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mn>,</mn></mrow></math> where x n and y n represent the coordinates of each pixel, and x 0 and y 0 are the coordinates of the barycenter. The track length value is calculated at the best applicable energy point of the traditional algorithm, which is approximately a Gaussian distribution, and the mean value of the Gaussian distribution is used as the initial threshold for determining if the track length is long or short. If the track length is smaller than this value, the track is a short track. Otherwise, it is a long track. At the same time, the mean value of the Gaussian distribution is used as one of the parameters that can be adjusted later. Furthermore, for tracks that are too short, a re-division needs to be performed along a direction that makes the half-zone of the starting track longer. The spacing between the green lines is set as a modifiable parameter. The degree of curvature of the track can be determined by calculating <math id="M8"><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mtext>pixels</mtext></mrow></msub></mstyle><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></math> . Unlike the third-order moment calculation, the influence of the charge value is removed in the second-order moment calculation and only the shape of the track is considered, thus preventing the uneven charge distribution of the track from affecting the results. The variation in the second-order moment is expressed by <math display="block" id="M9"><mrow><mrow><mo>|</mo><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mtext>out</mtext></mrow></munder><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup><mo>−</mo></mrow></mstyle><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mtext>last</mtext></mrow></munder><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></mstyle></mrow><mo>|</mo></mrow><mo>/</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mtext>pixels</mtext></mrow></munder><mrow><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup><mn>,</mn></mrow></mstyle></mrow></math> where <math id="M10"><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mtext>cut</mtext></mrow></msub></mstyle><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></math> is the summation of the second-order moments corresponding to the new green line after the move, and <math id="M11"><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mtext>last</mtext></mrow></msub></mstyle><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></math> is the second-order moment summation obtained after the last move. The reason for subtracting the two summations and taking the absolute value is because the result represents the amount of change in the second-order moment, and only the change in the trail needs to be detected. Each time the green line is moved, the point of charge at the other end of the green line is not involved in the calculation of the new second-order moment summation. The initial second-order moment sum of the first move is therefore calculated as ( <math id="M12"><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mtext>pixels</mtext></mrow></msub></mstyle><msubsup><mi>M</mi><mn>2</mn><mrow><mtext>min</mtext></mrow></msubsup></mrow></math> ). The variation is summed over the initial second-order moment to normalize the variation of the second-order moments for different tracks. The second-order moments of the green lines are computed in sequence, and the final split position chosen is the position with the most dramatic change. The above equation provides the maximum value at which a split will occur. To prevent the green line from dividing the outside of the track, a minimum track head length was set to preserve the first part of the track and the green line was forbidden from crossing this position. The minimum track head length is also a variable that can be adjusted. The actual parameters used were obtained after several rounds of empirical testing. 4) The barycenter and center of gravity line are calculated in the remaining half of the partition. The center of gravity is identified in the diagram by the interaction point (marked by an asterisk in Fig. 6(a) Fig. 6(a) )generated by the algorithm. The remaining track gravity line is illustrated in red in the diagram, while the arrow indicates the polarization direction obtained via the algorithm because the initial direction of photoelectron emission is also the polarization angle. Figure 6 Figure 6 shows the tracks obtained from the simulation and the reconstruction process. Some of these parameters can be modified at the beginning of the reconstruction, including the pixel blocks used for filtering circular tracks, and the value of the minimum preserved track head length. These parameters have a direct impact on all aspects of the reconstruction, including the selection of pre-reconstruction events. Monte Carlo simulations were used to generate a certain number of photoelectron tracks. The above algorithm was then used to analyze the tracks line by line. The initial emission angle of each track was obtained through parameter optimization. The overall angle distribution is shown in Fig. 7 Fig. 7 , where the horizontal coordinate represents the angle and the vertical coordinate represents the number of instances. 4.3 Event fitting analysis 2 The angular distribution of the K-shell photoelectron emission in the non-relativistic region is theoretically given by [ 26 26 ]: <math display="block" id="M13"><mrow><mfrac><mrow><mtext>d</mtext><msubsup><mi>σ</mi><mrow><mtext>ph</mtext></mrow><mi>K</mi></msubsup></mrow><mrow><mtext>d</mtext><mi>Ω</mi></mrow></mfrac><mo>∝</mo><mfrac><mrow><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>sin</mi></mrow></mstyle><mn>2</mn></msup><mi>ϕ</mi><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>cos</mi></mrow></mstyle><mn>2</mn></msup><mi>θ</mi></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>β</mi><mi>cos</mi><mi>ϕ</mi><mo stretchy="false">)</mo></mrow><mn>4</mn></msup></mrow></mfrac><mn>,</mn></mrow></math> where β is the ratio of the photoelectron velocity to the speed of light, θ is the azimuth angle of the X-ray electric vector, and φ is the polar angle with respect to the incident X-ray direction. The corresponding equation for the response amplitude after projecting the track onto a two-dimensional plane for processing is <math display="block" id="M14"><mrow><mi>M</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>A</mi><mo>×</mo><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>cos</mi></mrow></mstyle><mn>2</mn></msup><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mi>B</mi><mo stretchy="false">)</mo><mo>+</mo><mi>C</mi><mn>.</mn></mrow></math> The emission direction of polarized photoelectrons obeys such a distribution in theory. This equation was also used to fit the data shown in Fig. 7 Fig. 7 . In practice, modulation factors are often used to describe the polarization of the resulting photons because of slight deviations due to limited detection capability and the precision of the algorithm itself. The modulation factor is calculated as <math display="block" id="M15"><mrow><msub><mi>f</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>M</mi><mrow><mtext>max</mtext></mrow></msub><mo>−</mo><msub><mi>M</mi><mrow><mtext>min</mtext></mrow></msub></mrow><mrow><msub><mi>M</mi><mrow><mtext>max</mtext></mrow></msub><mo>+</mo><msub><mi>M</mi><mrow><mtext>min</mtext></mrow></msub></mrow></mfrac><mn>,</mn></mrow></math> where M max and M min correspond to the maximum and minimum values of the angle distribution in Fig. 7 Fig. 7 . The track reconstruction efficiency and the corresponding pre-screening efficiency of the SXP are shown in Fig. 8 Fig. 8 . The low-energy portion is difficult to reconstruct because the track lengths are too short. The middle energy portion is the most suitable area for the algorithm. At higher energies and longer track lengths, the modulation factor decreases because the parameters of the algorithm need to be adjusted. The circles show the efficiency of the pre-selection, and the ordinate the proportion of tracks retained for reconstruction. It can be observed that only a few cases have been screened out, and that the vast majority have been retained for reconstruction. Based on the conditions set in the simulations, a 4.5 keV X-ray polarization source was used in the experiment to reconstruct the measurements. Figure 9 Figure 9 shows the actual structure of the detector. The upper part of the structure is the SXP detector, while the center is the THGEM. The bottom part is the Topmetal-2 - . Both sides were equipped with gas pipes for ventilation and power supply connection wires for the application of high voltages. The signal collected by the Topmetal-2 - is shown in Fig. 10 Fig. 10 . The experimental tracks were reconstructed, and their angular distributions are shown in Fig. 11 Fig. 11 . The experimental tracks yielded the modulation factor of ∼ 28.60 ± 0.86%, while the event rejection rate was 46.68%. The modulation factor of the unpolarized X-ray source was 2:057%. The tracks obtained from the THGEM experiment were very rough and poorly reconstructed because of the large pore size. There were very few tracks that could be easily reconstructed, as shown in Fig. 10 Fig. 10 . In addition to the event selection rules described in Section 4.1 (the range for the circular track judgment was nine pixels), there were tracks distributed at the edge of the image in the experiment that were not involved in the reconstruction. For a track to be reconstructed, the track must have a sufficient number of pixels in the central area. Both the size of the central area and the number of pixels used as the threshold are adjustable input parameters. The main reason for the inconsistency between the experimental and simulation results is the influence of the experimental background; for example, the X-ray sources were not 100% polarized, and readout noise or dark current were present. Meanwhile, the spatial resolution was insufficient, resulting in possible inaccuracies in the reconstructed photoelectron emission direction. The spatial resolution will be optimized in future studies. 5 Discussion and Summary 1 The SXP adopts the most popular soft X-ray polarization detection method at present. The appropriate experimental gas and structure were selected through experiments and simulations. The iteration cut algorithm can provide reconstruction angles that are important for measuring the polarization and relevant for detection sources in the 210 energy range. The modulation factor reaches as high as 57 % above 7 . At low energies, the modulation factor is relatively low because of the short photoelectron tracks and the effect of electron diffusion in the THGEM avalanche. Experiments have shown that reducing the hole size of the THGEM structure can yield better results. The next step will be to optimize the hardware structure of the SXP and the algorithm, for example by using deep learning. It is envisaged that high-precision X-ray polarization detection will have a significant impact in revealing the physical mechanisms of astrophysical sources in the future. References H. Kestenbaum , G. Cohen , K. Long , et al. , The graphite crystal x-ray spectrometer on oso-8 . The Astrophysical Journal 210 , 805 - 809 ( 1976 ). doi: 10.1086/154889 http://doi.org/10.1086/154889 The graphite crystal x-ray spectrometer on oso-8 E. Costa , P. Soffitta , R. Bellazzini , et al. , An efficient photoelectric x-ray polarimeter for the study of black holes and neutron stars . Nature 411 , 662 - 665 ( 2001 ). doi: 10.1038/35079508 http://doi.org/10.1038/35079508 An efficient photoelectric x-ray polarimeter for the study of black holes and neutron stars M.C. Weisskopf , B. Ramsey , S.L. O’Dell , et al. , The imaging x-ray polarimetry explorer (ixpe) . Results in Physics 6 , 1179 - 1180 ( 2016 ). doi: 10.1117/12.2235240 http://doi.org/10.1117/12.2235240 The imaging x-ray polarimetry explorer (ixpe) S. Zhang , A. Santangelo , M. Feroci , et al. , The enhanced x-ray timing and polarimetry mission—extp . SCIENCE CHINA Physics, Mechanics & Astronomy 62 , 29502 ( 2019 ). doi: 10.1007/s11433-018-9309-2 http://doi.org/10.1007/s11433-018-9309-2 The enhanced x-ray timing and polarimetry mission—extp P. Soffitta , X. collaboration , et al. , Xipe, the x-ray imaging polarimetry explorer: Opening a new window in the x-ray sky . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 873 , 21 - 23 ( 2017 ). doi: 10.1016/j.nima.2017.02.025 http://doi.org/10.1016/j.nima.2017.02.025 Xipe, the x-ray imaging polarimetry explorer: Opening a new window in the x-ray sky W. Iwakiri , J. Black , R. Cole , et al. , Performance of the praxys x-ray polarimeter . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 838 , 89 - 95 ( 2016 ). doi: 10.1016/j.nima.2016.09.024 http://doi.org/10.1016/j.nima.2016.09.024 Performance of the praxys x-ray polarimeter H. Feng , H. Li , X. Long , et al. , Re-detection and a possible time variation of soft x-ray polarization from the crab . Nature Astronomy 4 , 511 - 516 ( 2020 ). doi: 10.1038/s41550-020-1088-1 http://doi.org/10.1038/s41550-020-1088-1 Re-detection and a possible time variation of soft x-ray polarization from the crab C. Ilie , Gamma-ray polarimetry: A new window for the nonthermal universe . Publications of the Astronomical Society of the Pacific 131 , 111001 ( 2019 ). doi: 10.1088/1538-3873/ab2a3a http://doi.org/10.1088/1538-3873/ab2a3a Gamma-ray polarimetry: A new window for the nonthermal universe R. Bellazzini , A. Brez , E. Costa , et al. , Photoelectric x-ray polarimetry with gas pixel detectors . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 720 , 173 - 177 ( 2013 ). doi: 10.1016/j.nima.2012.12.006 http://doi.org/10.1016/j.nima.2012.12.006 Photoelectric x-ray polarimetry with gas pixel detectors F. Muleri , P. Soffitta , L. Baldini , et al. , Spectral and polarimetric characterization of the gas pixel detector filled with dimethyl ether . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 620 , 285 - 293 ( 2010 ). doi: 10.1016/j.nima.2010.03.006 http://doi.org/10.1016/j.nima.2010.03.006 Spectral and polarimetric characterization of the gas pixel detector filled with dimethyl ether T. Li , M. Zeng , H. Feng , et al. , Electron track reconstruction and improved modulation for photoelectric x-ray polarimetry . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 858 , 62 - 68 ( 2017 ). doi: 10.1016/j.nima.2017.03.050 http://doi.org/10.1016/j.nima.2017.03.050 Electron track reconstruction and improved modulation for photoelectric x-ray polarimetry T. Kitaguchi , K. Black , T. Enoto , et al. , An optimized photoelectron track reconstruction method for photoelectric x-ray polarimeters . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 880 , 188 - 193 ( 2018 ). doi: 10.1016/j.nima.2017.10.070 http://doi.org/10.1016/j.nima.2017.10.070 An optimized photoelectron track reconstruction method for photoelectric x-ray polarimeters F. Sauli , Gem: A new concept for electron amplification in gas detectors . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 386 , 531 - 534 ( 1997 ). doi: 10.1016/S0168-9002(96)01172-2 http://doi.org/10.1016/S0168-9002(96)01172-2 Gem: A new concept for electron amplification in gas detectors R. Chechik , A. Breskin , C. Shalem , et al. , Thick gem-like hole multipliers: properties and possible applications . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 535 , 303 - 308 ( 2004 ). doi: 10.1016/j.nima.2004.07.138 http://doi.org/10.1016/j.nima.2004.07.138 Thick gem-like hole multipliers: properties and possible applications M. An , C. Chen , C. Gao , et al. , A low-noise cmos pixel direct charge sensor, topmetal-2- . Nuclear Instruments and Methods in Physics Research Section A-accelerators Spectrometers Detectors and Associated Equipment 810 , 144 - 150 ( 2016 ). doi: 10.1016/j.nima.2015.11.153 http://doi.org/10.1016/j.nima.2015.11.153 A low-noise cmos pixel direct charge sensor, topmetal-2- J. Cooper , Photoelectron-angular-distribution parameters for rare-gas subshells . Physical Review A 47 , 1841 ( 1993 ). doi: 10.1103/PhysRevA.47.1841 http://doi.org/10.1103/PhysRevA.47.1841 Photoelectron-angular-distribution parameters for rare-gas subshells Q. Liu , H. Liu , S. Chen , et al. , A successful application of thinner-thgems . Journal of Instrumentation 8 , C11008 ( 2013 ). doi: 10.1088/1748-0221/8/11/C11008 http://doi.org/10.1088/1748-0221/8/11/C11008 A successful application of thinner-thgems H.B. Liu , Y.H. Zheng , Y.G. Xie , et al. , The performance of thinner-thgem . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 659 , 237 - 241 ( 2011 ). doi: 10.1016/j.nima.2011.09.010 http://doi.org/10.1016/j.nima.2011.09.010 The performance of thinner-thgem H.B. Liu , Y.H. Zheng , Y.G. Xie , et al. , Study of the thgem detector with a reflective csi photocathode . Chinese Physics C 35 , 363 ( 2011 ). doi: 10.1088/1674-1137/35/4/008 http://doi.org/10.1088/1674-1137/35/4/008 Study of the thgem detector with a reflective csi photocathode H.B. Liu , Q. Liu , S. Chen , et al. , A study of thinner-thgem, with some applications . Journal of Instrumentation 7 , C06001 ( 2012 ). doi: 10.1088/1748-0221/7/06/C06001 http://doi.org/10.1088/1748-0221/7/06/C06001 A study of thinner-thgem, with some applications S. Agostinelli , J. Allison , K.a. Amako , et al. , Geant4—a simulation toolkit . Nuclear instruments and methods in physics research section A: Accelerators, Spectrometers, Detectors and Associated Equipment 506 , 250 - 303 ( 2003 ). doi: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 Geant4—a simulation toolkit R. Veenhof , Garfield, recent developments . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 419 , 726 - 730 ( 1998 ). doi: 10.1016/S0168-9002(98)00851-1 http://doi.org/10.1016/S0168-9002(98)00851-1 Garfield, recent developments C. Multiphysics , The platform for physics-based modeling and simulation . Avalible at: http://www.comsol.com/comsol-multiphysics. http://www.comsol.com/comsol-multiphysics. The platform for physics-based modeling and simulation R. Bellazzini , G. Spandre , Photoelectric polarimeters . X-ray Polarimetry: A New Window in Astrophysics 19 ( 2010 ). Photoelectric polarimeters D. Pfeiffer , L. De Keukeleere , C. Azevedo , et al. , Interfacing geant4, garfield++ and degrad for the simulation of gaseous detectors . Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 935 , 121 - 134 ( 2019 ). doi: 10.1016/j.nima.2019.04.110 http://doi.org/10.1016/j.nima.2019.04.110 Interfacing geant4, garfield++ and degrad for the simulation of gaseous detectors M. Gavrila , Relativistic k-shell photoeffect . Physical Review 113 , 514 ( 1959 ). doi: 10.1103/PhysRev.113.514 http://doi.org/10.1103/PhysRev.113.514 Relativistic k-shell photoeffect All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Xuefeng Huang, Jin Zhang, Huanbo Feng and Xichen Cai. The first draft of the manuscript was written by Xuefeng Huang and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. 10.1007/s41365-021-00903-0 This work was supported by the National Natural Science Foundation of China (Nos. U1731239, 12027803, 11851304, U1938201, 11575193, and U1732266) , the Guangxi Science Foundation (Nos. 2018GXNSFGA281007, 2017AD22006, 2018JJA110048) , and Key Research Program of Frontier Sciences, CAS (No. QYZDB-SSW368 SLH039) . 2020-10-09 2021-04-29 2021-05-07 2021-07-06 2021-07-01 2021-08-05 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2021 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-021-00903-0 http://dx.doi.org/10.1007/s41365-021-00903-0 . 2021 Nuclear Science and Techniques 07 32