Journal of Instrumentation       PAPER • OPEN ACCESS Demonstration of Gd-GEM detector design for neutron macromolecular crystallography applications To cite this article: D. Pfeiffer et al 2023 JINST 18 P04023   View the article online for updates and enhancements. You may also like The IFMIF-DONES Irradiation Modules Frederik Arbeiter, Urszula Wicek, Beatriz Brañas et al. - Radially excited pion: electromagnetic form factor and the box contribution to the muon's g-2 Ángel Salvador Miramontes López, Khepani Raya, Adnan Bashir et al. - DESI 2024 III: baryon acoustic oscillations from galaxies and quasars A.G. Adame, J. Aguilar, S. Ahlen et al. - This content was downloaded from IP address 194.27.118.1 on 30/07/2025 at 12:33 https://doi.org/10.1088/1748-0221/18/04/P04023 https://iopscience.iop.org/article/10.1088/1741-4326/add172 https://iopscience.iop.org/article/10.1088/1674-1137/add259 https://iopscience.iop.org/article/10.1088/1674-1137/add259 https://iopscience.iop.org/article/10.1088/1674-1137/add259 https://iopscience.iop.org/article/10.1088/1475-7516/2025/04/012 https://iopscience.iop.org/article/10.1088/1475-7516/2025/04/012 https://pagead2.googlesyndication.com/pcs/click?xai=AKAOjstKtd1z9mB9uyzcxdZqGJPkP_Myj5uVi2CDiwRM0ns_b0x67fwsfy3xXrb2P1Ob0-oxNmFgjpv-DpghuiZ6j3TkZbQZh_WAslFgDvPg-Ui2GFTPrqzD254fKBwyriKt8m-cauiLLF6AtWfQqWcwXVwCxyP1J6LuHXIUyG8naJZD8--0o2P8LsEhzPT4d5IDtkvTzpPnQzk6t5mEur63byIlSFHpyh9S94o4CfMxwai8v957DDxRp9ZrIwhabpINtx0zXTQ5zcqidC7pVVQTO9aoCyC4qQxq50hMwv2ueyUv1TT893r68mxkz4Bmn4xlkuGVIG8zsp1LqqUrTMWeBinmyFy1NCxvGItBXXqZISLhF0PAl6vMJvG1&sig=Cg0ArKJSzCUBZhDGerfM&fbs_aeid=%5Bgw_fbsaeid%5D&adurl=https://www.electrochem.org/248/registration%3Futm_source%3DIOP%26utm_medium%3Dbanner%26utm_campaign%3DIOP_248_Early_Reg%26utm_id%3DIOP%2B248%2BEarly%2BRegistration 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 Published by IOP Publishing for Sissa Medialab Received: November 16, 2022 Revised: March 16, 2023 Accepted: March 27, 2023 Published: April 18, 2023 Demonstration of Gd-GEM detector design for neutron macromolecular crystallography applications D. Pfeiffer,𝑎,𝑏,𝑐,∗ F. Brunbauer,𝑏 V. Cristiglio,𝑑 R. Hall-Wilton,𝑎,𝑐,𝑒 M. Lupberger, 𝑓 ,𝑔 M. Markó,ℎ,𝑖 H. Muller,𝑏, 𝑓 E. Oksanen,𝑎 E. Oliveri,𝑏 L. Ropelewski,𝑏 A. Rusu, 𝑗 J. Samarati,𝑎,𝑏 L. Scharenberg,𝑏,𝑒 M. van Stenis,𝑏 P. Thuiner𝑏,𝑘 and R. Veenhof𝑏,𝑙 𝑎European Spallation Source ESS ERIC (ESS), Box 176, SE-221 00 Lund, Sweden 𝑏European Organization for Nuclear Research (CERN), 1211 Geneva 23, Switzerland 𝑐Department of Physics, University of Milano-Bicocca, Piazza della Scienza 3, 20126 Milan, Italy 𝑑Institut Laue-Langevin, 71 avenue des Martyrs CS 20156, 38042 Grenoble, France 𝑒Sensors & Devices Centre, Fondazione Bruno Kessler, via Sommarive 18, 38123 Trento, Italy 𝑓 Physikalisches Institut, University of Bonn, Nußallee 12, 53115 Bonn, Germany 𝑔Helmholtz-Institut für Strahlen- und Kernphysik, University of Bonn, Nußallee 14–16, 53115 Bonn, Germany ℎNeutron Spectroscopy Department, Institute for Energy Security and Environmental Safety, Centre for Energy Research, 29–33 Konkoly-Thege Miklós út, Budapest, 1121, Hungary 𝑖Budapest Neutron Centre, 29–33 Konkoly-Thege Miklós út, Budapest, 1121, Hungary 𝑗SRS Technology, 30 Promenade des Artisans, 1217 Meyrin, Switzerland 𝑘Vienna University of Technology, 1040 Vienna, Austria 𝑙Bursa Uludağ University, Görükle Kampusu, 16059 Niüfer/Bursa, Turkey E-mail: dorothea.pfeiffer@ess.eu ∗Corresponding author. c© 2023 CERN. Published by IOP Publishing Ltd on behalf of Sissa Medialab. Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. https://doi.org/10.1088/1748-0221/18/04/P04023 mailto:dorothea.pfeiffer@ess.eu http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ https://doi.org/10.1088/1748-0221/18/04/P04023 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 Abstract: The European Spallation Source (ESS) in Lund, Sweden will become the world’s most powerful thermal neutron source. The Macromolecular Diffractometer (NMX) at the ESS requires three 51.2 × 51.2 cm2 detectors with reasonable detection efficiency, sub-mm spatial resolution, a narrow point-spread function (PSF), and good time resolution. This work presents measurements with the improved version of the NMX detector prototype consisting of a Triple-GEM (Gas Electron Multiplier) detector with a natural Gd converter and a low material budget readout. The detector was successfully tested at the neutron reactor of the Budapest Neutron Centre (BNC) and the D16 instrument at the Institut Laue-Langevin (ILL) in Grenoble. The measurements with Cadmium and Gadolinium masks in Budapest demonstrate that the point-spread function of the detector lacks long tails that could impede the measurement of diffraction spot intensities. On the D16 instrument at ILL, diffraction spots from Triose phosphate isomerase w/ 2-phosphoglycolate (PGA) inhibitor were measured both in the MILAND Helium-3 detector and the Gd-GEM. The comparison between the two detectors shows a similar point-spread function in both detectors, and the expected efficiency ratio compared to the Helium-3 detector. Both measurements together thus give good indications that the Gd-GEM detector fits the requirements for the NMX instrument at ESS. Keywords: Neutron detectors (cold, thermal, fast neutrons); Micropattern gaseous detectors (MSGC, GEM, THGEM, RETHGEM, MHSP, MICROPIC, MICROMEGAS, InGrid, etc) ArXiv ePrint: 2211.04335 https://arxiv.org/abs/2211.04335 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 Contents 1 Introduction 1 2 The NMX detector prototype 2 3 Test beam in Budapest 3 4 Test beam at ILL Grenoble 4 5 Conclusion 9 A Annex: detailed geometrical calculations 9 1 Introduction The European Spallation Source (ESS) [1] in Lund, Sweden is foreseen to start the scientific user programme in 2025, and the construction phase will be complete by 2027. It will become the world’s most powerful thermal neutron source with a significantly higher brilliance than at existing reactor sources such as the Institut Laue-Langevin (ILL) [2] or other spallation sources such as the Spallation Neutron Source (SNS) [3] and the Japan Proton Accelerator Research Complex (J-PARC) [4]. 15 neutron instruments constitute the initial instrument suite at ESS [5, 6]. For each of the instruments efficient thermal neutron detectors are a crucial component [7]. The Macromolecular Diffractometer NMX [6] requires three 51.2×51.2 cm2 detectors with reasonable detection efficiency (around 25% for cold neutrons in the center of the wavelength range from 1.8 to 10 Å) and sub-mm spatial resolution [8]. While spatial resolution is often used as a requirement for macromolecular crystallography detectors, it is worth noting that the shape and width of the point-spread function have an even larger effect on the ability to reliably measure the intensities of diffraction spots. Neutron image plates [9] that are used at reactor sources macromolecular crystallography instruments have a sufficient spatial resolution of 200 μm, but lack the time resolution needed at spallation source instruments. Scintillation-based detectors [10, 11] offer time resolution, but are currently limited to about 1 mm spatial resolution. For spallation source instruments the combination of solid neutron converters with Micro Pattern Gaseous Detectors (MPGDs) [12, 13] is a promising option to achieve a similar spatial resolution to the neutron image plate with sufficient time resolution, high-rate capabilities and a reasonable neutron detection efficiency. A time-resolved detector is necessary to be able to carry out a Time-of-Flight analysis of the neutrons in the diffraction spots, to disentangle the spatial overlap of reflections that occur at the same location, but at different wavelengths. The nanosecond time resolution of the detector and readout system [14] is far better than the required microsecond resolution. The rate capabilities of Triple-GEM detectors [15] and the readout system [16] together – 1 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 support the requirements for NMX, which are estimated to be less than one MHz per cm2. To obtain a good spatial resolution, the 𝜇TPC data analysis technique, based on the principle of the Time Projection Chamber (TPC), was developed. Charge produced close to the neutron converter has to drift to the amplification stage of the detector, whereas charge produced close to the first GEM is directly amplified. The 𝜇TPC method retains therefore the x and y-strip with the largest time as start of the ionisation track and thus as position where the neutron interacted in the neutron converter. The method has first been successfully applied to MPGDs with 10B4C converters [17]. A position resolution with a 𝜎 better than 200 μm was obtained by determining the starting point of the neutron-induced ionization track [18]. To improve the low detection efficiency (less than 5% at 1.8 Å), the 10B4C converter has been replaced with a more efficient gadolinium converter. A prototype of the NMX detector [19] consisting of a Triple-GEM detector [20] with Gd converter has first been tested in 2016 at the R2D2 beamline at the JEEP II research reactor at the Institute for Energy Technology (IFE) in Kjeller/Norway [21]. With the natural Gd converter, a position resolution with a 𝜎 better than 400 μm was measured. The properties of natural and isotope-enriched Gd (157Gd) neutron converters have been studied in detail in work package 4 of the European Horizon2020 grant BrightnESS [22]. This work presents measurement with the improved version of the NMX detector prototype at the neutron reactor of the Budapest Neutron Centre (BNC) [23] and at the D16 instrument at the Institut Laue-Langevin (ILL) [2] in Grenoble in the year 2018. The purpose of these measurements was to determine the spatial resolution of the detector and the point-spread function (PSF) with realistic samples. As it is a well-established practice in crystallography to approximate the PSF with a Gaussian [24], the paper subsequently uses Gaussian fits for this purpose. 2 The NMX detector prototype A schematic drawing of the NMX detector is displayed in figure 1(a). The NMX detector is a triple GEM detector [20] with a low material budget 𝑥/𝑦 strip readout. In contrast to the traditional 𝑥/𝑦 strip readout that contains an FR4 base [25], the low material budget readout uses a 100 μm thick Kapton foil to reduce the neutron scattering from about 20% to about 5%. The active area of the final NMX detector is 51.2 cm × 51.2 cm, with 1280 strips in 𝑥 and 𝑦. The strip pitch is 400 μm. In the present article though, a small prototype with 10 cm × 10 cm active area (256 strips in 𝑥 and 𝑦) has been tested. The detector was flushed with Ar/CO2 70/30 mixture at 5 l/h at room temperature and atmospheric pressure. The conversion volume or drift space was 10 mm long, and as cathode and neutron converter served a 25 μm Gadolinium foil. A drift field of 700 V/cm was chosen to avoid electron attachment to electronegative impurities and the subsequent loss of primary ionization electrons, while keeping the drift velocity smaller than 2.0 cm/μ at atmospheric pressure. The three GEM foils [26] in the detector are mounted at a distance of 2 mm from each other and powered with a resistor chain. The values of the resistors are stated in figure 1(a). The detector was operated at an effective gain of about 5000. This small prototype is read out with two RD51 VMM3a hybrids each in x and y directions (figure 1(b)). Two VMM3a ASICs [27, 28] with 64 channels per chip are mounted on a RD51 VMM3a hybrid [16, 29], resulting in 128 channels per hybrid. The Scalable Readout System (SRS) [30] serves as readout system to acquire the data from the RD51 [31] VMM3a hybrids. – 2 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 The detector is operated in backwards mode. Neutrons traverse the detector readout and the three GEM foils, and then impinge on the natural Gd cathode. Natural Gd contains 14.8% of 155Gd and 15.7% of 157Gd. Both of these isotopes have a very high neutron capture cross section [32]. After the neutron capture, gadolinium releases prompt gamma particles with an energy of up to 9 MeV and conversion electrons with energies ranging from 29 keV to 250 keV. The signals resulting from these conversion electrons stopped in the drift volume, are used to determine the position of the neutron impact on the cathode with the help of the uTPC method [18, 19]. 550 kΩ 500 kΩ 450 kΩ 1 MΩ 1 MΩ 1 MΩ 10 MΩ 2 mm 2 mm 2 mm …… 10 MΩ 10 MΩ 10 mm conversion electron neutron HV1: 4000 V HV2: 3300 V Gd sheet (25 um) (a) Schematic view of Triple-GEM detector (b) Gd-GEM 10 cm × 10 cm prototype Figure 1. Gd-GEM neutron detector. 3 Test beam in Budapest The detector prototype was tested in 2018 at the neutron reactor of the Budapest Neutron Centre (BNC). A monochromatic neutron beam at the ATHOS [33] triple-axis spectrometer with a wavelength of 3.5 Å was used for the measurements. Figure 2(a) shows a Cadmium mask with circular holes of 1.0 mm diameter. The holes of the mask have a minimum separation (centre to centre) of about 2 mm horizontally, 1.6 mm vertically, and 1.3 mm diagonally. To record this hole pattern with the detector, the beam was collimated to 1.1 cm in 𝑥 direction and 3.3 cm in 𝑦 direction. First data were acquired with the collimated beam directly hitting the detector without mask for 4 minutes. Then the mask was taped to the readout of the detector, and data were acquired for 84 minutes. To get rid of detector and electronics-specific inhomogeneities, the mask data were normalized with the direct beam data. Figure 2(b) shows the recorded image of the mask. The collimated beam was directed at the lowest three rows of the mask. The holes there can be clearly resolved. To further study the spatial resolution of the detector, the centre part of the readout of the detector was covered with a 25 μm thick Gd foil with a single 1 mm diameter hole. Data were acquired for 135 minutes. Figure 3(b) shows the measurement. As in the case of the Cadmium mask data, the single-hole data has been normalized with the direct beam data. The Gaussian (𝜎 600 μm) fit of the profile in 𝑥 direction is displayed in figure 3(b). The error bars indicate the statistical uncertainty on the counts. The Full-Width-Half-Maximum (FWHM) is 1.4 mm, and the full width at 5% of the height is 2.9 mm. – 3 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 (a) Cd mask, 1.0 mm holes 155 160 165 170 175 180 x strips [pitch 400 um] 180 185 190 195 200 205 210 215 220 y st rip s [p itc h 40 0 um ] 0 0.2 0.4 0.6 0.8 1 1.2 1.4 Entries : 71520337152033 (b) Measurement with Cd mask Figure 2. Recorded pattern of a Cd mask with circular holes of 1.0 mm diameter. Figure a) show a picture of the mask that was installed on the detector. The holes have a minimum separation (centre to centre) of about 2 mm horizontally, 1.6 mm vertically, and 1.3 mm diagonally. Figure b) shows the beam intensity recorded with the mask normalized to the beam intensity without mask. The beam was focused on the lowest three rows of the mask. The holes there can be clearly resolved. Figure 3(b) and the corresponding width suggests that the point-spread function of the detector lacks long tails that could impede the measurement of diffraction spot intensities. While a hole in a mask attached to the detector surface is not entirely comparable to a diffraction spot, it is a good measure of the point-spread function of such a spot due to detector instrumental effects. This result indicates that the detector is suitable for measuring the diffraction spot intensities in crowded diffraction patterns. 4 Test beam at ILL Grenoble Also in 2018 the detector has been tested against the MILAND detector (Helium-3 based wire chamber with a wire plane in 𝑥 and 𝑦 direction) of the D16 [34] instrument at the Institut Laue- Langevin (ILL) [2] in Grenoble. Monochromatic neutrons with a wavelength of 4.51 Å were used for the measurements. The beamline of the D16 instrument is displayed in figure 4(a). In the foreground of the picture, the sample and the Gd-GEM detector are visible. On the left-hand side the circular vessel of the D16 3He wire chamber can be seen in its retracted position. The neutron beam traverses a 3 mm natural B4C diaphragm and hits the sample (Triose phosphate isomerase w/ 2-phosphoglycolate (PGA) inhibitor) [35] (figure 4(b)). The Gd-GEM detector in its shielding box was mounted at a distance of 140 mm from the sample. The cathode of the detector, where the impinging neutrons are converted into conversion electrons, is 56 mm away from the beginning of the detector shielding box. The cathode of the detector was thus at a distance of 196 mm from the sample. – 4 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 re la tiv e in te n s ity 156 158 160 162 164 166 168 170 172 174 176 x coordinate [pitch 400 um] 164 166 168 170 172 174 176 178 180 182 184 y c o o rd in a te [ p itc h 4 0 0 u m ] Profile x (a) Measurement with Gd-foil mask with 1 mm hole 63 64 65 66 67 68 69 70 71 Position [mm] 0 20 40 60 80 100 120 Re la tiv e In te ns ity [% ] Fit parameter bkgd: 8.5 +/-0.5 height: 91 +/-10 mu: 66.6 +/-0.04 sigma: 0.6 +/-0.04 FWHM: 1.4 mm FW 5% max: 2.9 mm Profile Fit: bkgd + (height bkgd) * e (x )2/(2 * 2) (b) Profile with 1.4 mm FWHM Figure 3. Measurement with Gd-foil with 1 mm hole in front of detector readout. The Gaussian fitted to the profile has a FWHM of 1.4 mm. The binning corresponds to the 400 μm strip pitch of the detector readout. (a) D16 beamline with detector (b) B4C diaphragm and sample holder Figure 4. Beamline of D16 instrument at ILL in Grenoble. The Gd-GEM detector can be seen in the center of 4(a), whereas the MILAND detector is the circular chamber visible on the outer left. The neutron beam traverses a 3 mm 𝐵4𝐶 diaphragm and hits the PGA inhibitor sample. The Gd-GEM detector is mounted at a distance of 196 mm from the sample. The Gd-GEM and the MILAND detector were both fixed to a rotatable table seen at the bottom left corner. To take data with the MILAND detector, the Gd-GEM and its support have to be removed, since they block the neutrons coming from the sample. Subsequently, the MILAND detector is then moved forward towards the sample into its data acquisition position (332 mm from the sample). The – 5 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 angle 𝛾 indicates the rotation of the detectors with respect to the end of the beamline as indicated in figure 5(a). The sample itself is mounted on a goniometer, and can also be rotated around its axis. The angle 𝜔 thus indicates the rotation of the sample, as explained in figure 5(b). (a) Detector rotation 𝛾 (b) Sample orientation 𝜔 Figure 5. The detector and sample orientation (angles 𝛾 and 𝜔) can be adjusted from 0 to 360 degrees. The setup that was used at the D16 instrument was not optimized to study the unit cell of PGA inhibitor. The Gd-GEM detector could not be installed at the same location as the MILAND detector. Since the Gd-GEM with an active area of 10 cm x 10 cm is also nine times smaller than the MILAND detector, the Gd-GEM was installed closer to the sample. To summarize, it is not obvious which reflection measured on the MILAND detector belongs to which reflection on the Gd-GEM. Considering the different solid angles covered by the two detectors, geometrical calculations were needed to determine the matching reflection pairs. The details of the geometrical calculations can be found in the annex (section A). Figure 6 compares at the detector rotation 𝛾 = 60 ° the measurements in the MILAND detector with the Gd-GEM detector. The sample orientation was changed from 𝜔 = 5 ° in figure 6(a) to 𝜔 = 55 ° in figure 6(b). Whereas in the MILAND detector several diffraction spots are visible (8 spots at the right and 3 spots at the left), the Gd-GEM detector recorded only 1 spot in both measurements due to its smaller physical size. It is apparent that the centre of the Gd-GEM detector has been misaligned with respect to the centre of the MILAND detector. To investigate the shape and width of the diffraction spots, a Gaussian has been fitted to the profile in x-direction. Figure 7 displays the same diffraction spot (𝜔 = 55 °, 𝛾 = 60 °), on the left in the MILAND 3He detector and in the Gd-GEM on the right. The MILAN 3He detector acquired data for 300 s, whereas the Gd-GEM took data for 1620 s. Since the data of the Gd-GEM has been background subtracted, the signal to background ratio is slightly better in the Gd-GEM with 9.9 compared to 2.5 for the MILAND detector. The FWHM for the profile in the MILAND detector amounts to 3.5 mm, whereas in the Gd-GEM one gets 1.8 mm. These values cannot be compared – 6 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 102.4 102.4 0 0 Gd-GEM He3 detector ω 5° 𝜸 60° x 83.5 mm y 40.4 mm (a) 𝜔 = 5 ° 102.4 He3 detector ω 55° 𝜸 60° 102.4 0 0 Gd-GEM x 73.0 mm y 44.5 mm (b) 𝜔 = 55 ° Figure 6. Measurement of the diffraction spots at 𝜔 = 5 ° and 𝜔 = 55 ° in the MILAND detector. The detector rotation was 𝛾 = 60 °. The MILAND detector acquired in both measurements data for 300 s. Whereas in the MILAND detector several diffraction spots are visible (8 spots at the right and 3 spots at the left), the smaller Gd-GEM detector recorded only 1 spot in both measurements. The centre of the Gd-GEM detector has been misaligned with respect to the centre of the MILAND detector. Table 1. Measurements at 𝛾 = 60 ° detector rotation and a sample rotation of 𝜔 = 5 ° and 𝜔 = 55 °. The table shows the distance between spots in the MILAND and Gd-GEM detector, and the equivalent detector rotation. The equal equivalent detector rotation shows that the spots belong to the same reflection in both detectors. Measured/Calculated Detector Spot distance Detector rotation Gd-GEM 11.46(39) mm 3.35(12) ° MILAND 3He 20.63(17) mm 3.56(3) ° directly, since the diffraction spots are not completely round, and the 𝑥-axis of the Gd-GEM detector was rotated 45.55 ° with respect to the 𝑥-axis of the MILAND detector. The Gd-GEM detector produces neutron-induced secondary charged particles (electrons) of energies between 10 keV and 250 keV. Their different range in gas should not produce a deformation of the diffraction spots, since with the help of the uTPC method, the interaction point of the neutron on the cathode is reconstructed. For the MILAND 3He detector, there might be a small parallax effect, since the diffracted beam is not always perpendicular to the detector surface, and the position of the neutron is calculated from the charge with the centre of gravity method. The fit of a two-dimensional Gaussian to the same diffraction spots is shown in figure 8. The 𝜎𝑥 of the Gaussian fit amounts to 1.83(9) mm for the 3He detector and 1.0(1) mm for the Gd-GEM, whereas the values for 𝜎𝑦 are 1.15(5) mm and 0.85(9) mm respectively. The 2D Gaussian of the MILAND detector is rotated with respect to the 2D Gaussian in the Gd-GEM, the 𝜎𝑥 is thus not the 𝜎 along the x-axis but just one of the principal axes of the Gaussian. When comparing now the two – 7 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 235.0 237.5 240.0 242.5 245.0 247.5 250.0 252.5 255.0 Position [mm] 0 20 40 60 80 100 Co un ts Fit parameter bkgd: 29.8 +/-1.4 height: 103 +/-6 mu: 244.4 +/-0.1 sigma: 1.5 +/-0.1 FWHM: 3.5 mm FW 5% max: 7.3 mm Profile Fit: bkgd + (height bkgd) * e (x )2/(2 * 2) (a) MILAND detector 69 70 71 72 73 74 75 76 77 Position [mm] 75 50 25 0 25 50 75 100 125 Co un ts Fit parameter bkgd: 8.1 +/-6.2 height: 88 +/-17 mu: 73.0 +/-0.2 sigma: 0.8 +/-0.2 FWHM: 1.8 mm FW 5% max: 3.7 mm Profile Fit: bkgd + (height bkgd) * e (x )2/(2 * 2) (b) Gd-GEM detector Figure 7. Gaussian fit of profile in 𝑥-direction of the same diffraction spot (𝜔 = 55 °, 𝛾 = 60 °) in the MILAND detector and the Gd-GEM. The 3He detector acquired data for 300 s, whereas the Gd-GEM took data for 1620 s. The error bars represent the statistical uncertainty on the counts. larger 𝜎𝑥 of both detectors, the 𝜎𝑥 of the MILAND detector is 1.8 times larger than the 𝜎𝑥 of the Gd-GEM. For the smaller 𝜎𝑦 , the 𝜎𝑦 of the MILAND detector is 1.4 times larger than the 𝜎𝑦 of the Gd-GEM. The Gd-GEM was installed at 196 mm from the sample, whereas the 3He detector was positioned at 332 mm. Assuming an equal detector resolution and ignoring parallax effects, the size of the spot should be around 1.7 times smaller in the location of the Gd-GEM. The 𝜎 derived from the 2D fit indicates thus a similar point-spread function for both detectors. This indicates the absence of instrumental effects on the reconstruction of the diffraction spot in the Gd-GEM detector. x [mm] 235 240 245 250 255 y [mm] 225 230 235 240 245 z [ co un ts ] 40 50 60 70 80 90 100 amplitude = 104(4) background = 31.2(3) x = 244.46(8) y = 233.84(6) x = 1.83(9) y = 1.15(5) (a) MILAND detector x [mm] 68 70 72 74 76 y [mm] 40 42 44 46 48 z [ co un ts ] 10 20 30 40 50 60 70 amplitude = 68(8) background = 3(1) x = 72.9(1) y = 44.48(9) x = 1.0(1) y = 0.85(9) (b) Gd-GEM detector Figure 8. Fit of two-dimensional Gaussian to the same diffraction spot (𝜔 = 55 °, 𝛾 = 60 °) in the MILAND detector and the Gd-GEM. The 3He detector acquired data for 300 s, whereas the Gd-GEM took data for 1620 s. – 8 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 In the following paragraph, as a rough plausibility check, the total number of counts in the diffraction spots in the MILAN detector and the Gd-GEM are compared. The integral of the two-dimensional Gaussian fit amounts to 980.6 counts/𝑚𝑚2 for the MILAND detector, and 359.6 counts/(0.4 mm)2 for the Gd-GEM detector. Considering the different acquisition times of 300 s for the MILAND detector, and 1620 s for the Gd-GEM, one arrives at a rate of 3.3 Hz/m𝑚2 for the MILAND detector and 1.4 Hz/mm2 for the Gd-GEM. According to this data analysis, the MILAND detector is thus about 2.4 times more efficient than the Gd-GEM. The Gd-GEM with the traditional x/y readout has a measured efficiency of 11.8% for neutrons with a wavelength of 2 Å [19]. For a Gd-GEM detector with low material budget readout, as used in the measurements presented here, the detection efficiency has been simulated with Geant4 [36] to be around 25% with a threshold of 0 keV at 4.51 Å. The 3He wire chamber has a detection efficiency of about 70% at a wavelength at 2.5 Å [37]. Assuming a similar or higher efficiency of the 3He wire chamber at 4.51 Å, the efficiency calculated from the measurement of the diffraction spots slightly overestimates the efficiency of the Gd-GEM detector. 5 Conclusion The improved prototype of the NMX detector consisting of a Triple-GEM detector with Gd converter has been successfully tested at the neutron reactor of the Budapest Neutron Centre (BNC) [23] and at the D16 instrument at the Institut Laue-Langevin (ILL) [2] in Grenoble. The measurements with Cadmium and Gadolinium masks in Budapest demonstrate that the point-spread function of the detector lacks long tails that could impede the measurement of diffraction spot intensities. This result indicates that the detector is suitable for measuring the diffraction spot intensities in crowded diffraction patterns, as detailed in reference [8]. On the D16 instrument at ILL, diffraction spots from Triose phosphate isomerase w/ 2-phosphoglycolate (PGA) inhibitor were measured both in the D16 MILAND detector and the Gd-GEM. The comparison between the two detectors shows a similar point-spread function in both detectors, and the expected efficiency ratio compared to the MILAND detector. The high rate capabilities [15, 16], the excellent time resolution of GEM detectors [14] (not evaluated during the tests described here) and the lack of parallax effects are thus good indications that the Gd-GEM detector fits the requirements for the NMX instrument at ESS. A Annex: detailed geometrical calculations This annex contains the geometrical calculations that proof that the diffraction spots on the MILAND detector and the Gd-GEM detector (as compared in section 4) come from the same reflection. Figure 9(a) displays the schematic drawing of the Gd-GEM detector as it is seen from the sample. The neutrons pass first through the detector readout, then the three GEM foils and arrive finally at the cathode. If the detector rotation angle 𝛾 is increased (detector moves to the right), the diffraction spot will move to the left towards larger 𝑥 and 𝑦 values. If seen from the cathode side as in figure 9(b), an increase in 𝛾 means that the detector moves to the left, and the diffraction spot moves to the right. To see the effect of the detector movement both in the 𝑥 and in the 𝑦 coordinate, the detector was installed with an angle 𝛽 of approximately 45 ° between the beamline and the 𝑥-axis. The principles of the calculations below are explained in figure 10. A schematic drawing of the shift of a diffraction spot (due to detector rotation) on the detector plane is displayed in figure 10(a). – 9 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 𝜸 0° 𝜸 90° increase 𝜸 x/y strip readout Ar/CO2 70/30 outlet HV connections Ar/CO2 70/30 inlet increase x increase y (a) View from sample. 𝜸 90° 𝜸 0° Gd converter Ar/CO2 70/30 outlet HV connections Ar/CO2 70/30 inlet increase 𝜸 increase x increase y (b) View from back of detector. Figure 9. Schematic drawing of detector readout (𝑥/𝑦 strip readout), VMM hybrids, Gd converter and movement of diffraction spots on detector depending on rotational movement 𝛾 of detector. The neutrons coming from the sample traverse first the readout of the detector, as shown previously in figure 1(a). The shift between the two spots can be calculated using the theorem of Pythagoras. The angle between the beamline and the 𝑥-axis can be determined with the Arctan (inverse Tangens) function. Figure 10(b) shows the effect of the detector rotation 𝛾 on the position of the diffraction spots. If the distance d between the detector and the sample is known, 𝛾 can be calculated with the Arctan function. Figure 10(c) illustrates the effect of the sample rotation by the angle 𝜔. This sample rotation leads to a shift of the diffraction spot that is equivalent to a rotation of the detector by the hypothetical angle 𝛼. y-a xis 0 25 5 yx x-axis 0 255 (x2,y2)(x1,y1) shift s (a) Shift of diffraction spots diffraction spot sample ne ut ro ns detector position 1 position 2 detector di st an ce d distance d (b) Rotation of detector diffraction spots sample orientation changed by (c) Rotation of sample Figure 10. Schematic drawing of diffraction spot shift (due to detector rotation) on the detector plane(Figure 10(a)). Figure 10(b) shows the effect of the detector rotation 𝛾 on the position of the diffraction spots. Figure 10(c) displays the effect of the rotation of the sample by the angle 𝜔. This sample rotation leads to a diffraction spot shift that is equivalent to a rotation of the detector by the hypothetical angle 𝛼. Figure 11 shows the moving of a diffraction spot from PGA inhibitor. At a sample rotation of 𝜔 = 5 °, the increase in the detector rotation angle 𝛾 from 𝛾 = 57 ° to 𝛾 = 60 °, moves the spot from – 10 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 the position x=75.5 mm and y=33.3 mm to x=83.5 mm and y=40.4 mm. The positions of the spot on the detector have been determined by fitting a Gaussian distribution to the data. A detector rotation of 𝛾 = 3 °, resulted thus in a shift 𝑠57 °−60 ° of the spots. 𝑠57 °−60 ° = √︃ (𝑥60 − 𝑥57)2 + (𝑦60 − 𝑦57)2 = √︁ (83.67(26) mm − 75.49(31) mm)2 + (40.39(17) mm − 33.28(18) mm)2 = 10.84(47) mm (A.1) (a) 𝛾 = 57 ° (b) 𝛾 = 60 ° Figure 11. Measurement of diffraction spots from PGA inhibitor at a sample rotation of 𝜔 = 5 °. An increase in the detector rotation angle 𝛾 moves the spot to larger 𝑥 and 𝑦 values. Using a shift 𝑠57 °−60 ° = 10.84(47) mm and a detector sample distance 𝑑 = 196 mm, the change of the detector rotation angle 𝛾57 °−60 ° amounts to 𝛾57 °−60 ° = arctan 𝑠57 °−60 ° 𝑑 × 180 ° 𝜋 = arctan 10.84(47) mm 196 mm × 180 ° 𝜋 = 3.17(14)° (A.2) The actual angle 𝛽57 °−60 ° between the beamline and the 𝑥-axis of the detector be calculated with 𝛽57 °−60 ° = arctan 𝑦60 − 𝑦57 𝑥60 − 𝑥57 × 180 ° 𝜋 = arctan 40.39(17) mm − 33.28(18) mm 83.67(26) mm − 75.49(31) mm × 180 ° 𝜋 = arctan 7.11(25) mm 8.18(40) mm × 180 ° 𝜋 = 40.99(171)° (A.3) – 11 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 Measurements with a sample rotation of 𝜔 = 55 ° and a change of the detector rotation angle 𝛾 = 60 ° to 𝛾 = 62 °, resulted in a shift 𝑠60 °−62 ° = 6.67(42) 𝑚𝑚, a calculated detector rotation angle 𝛾60 °−62 ° = 1.95(12)° and an angle 𝛽60 °−62 ° = 45.55(258)° between the beamline and the 𝑥-axis of the detector. All rotational shifts align thus with the expectations, which proves that there is no systematic bias in position reconstruction. Table 2 summarizes these findings. Table 2. Comparison of applied detector rotation and detector rotation calculated from the distance between diffraction spots. Applied Measured/Calculated Sample rotation Detector rotation Spot distance Detector rotation 𝛽beamline 𝑥-axis 𝜔 = 5 ° 𝛾 = 57 to 60 ° 10.84(47) mm 3.17(14) ° 40.99(171) ° 𝜔 = 55 ° 𝛾 = 60 to 62 ° 6.67(42) mm 1.95(12) ° 45.55(258) ° For the comparison between the MILAND detector with the Gd-GEM detector, the sample orientation was changed from 𝜔 = 5 ° to 𝜔 = 55 °. In the Gd-GEM detector, the change in sample rotation from 𝜔 = 5 ° to 𝜔 = 55 ° shifted the diffraction spot by 𝑠𝐺𝑑−𝐺𝐸𝑀 = √︃ (𝑥55 − 𝑥5)2 + (𝑦55 − 𝑦5)2 = √︁ (83.67(26) mm − 72.97(16) mm)2 + (40.39(17) mm − 44.49(24) mm)2 = 11.46(39) mm (A.4) This shift of 𝑠Gd-GEM = 11.46(39) mm, caused by the rotation of the sample, is equivalent to a hypothetical rotation 𝛼Gd-GEM of the detector of 𝛼Gd-GEM = arctan 𝑠𝐺𝑑−𝐺𝐸𝑀 𝑑 × 180 ° 𝜋 = arctan 11.46(39) 𝑚𝑚 196𝑚𝑚 × 180 ° 𝜋 = 3.35(12)° (A.5) On the imaging plane of the MILAND 3He-detector, installed at 332 mm from the sample, the same change in sample orientation shifted the diffraction spot by 𝑠He3 = √︃ (𝑥55 − 𝑥5)2 + (𝑦55 − 𝑦5)2 = √︁ (244.38(13) mm − 244.89(20) mm)2 + (233.83(11) mm − 213.21(11) mm)2 = 20.63(17) mm (A.6) The hypothetical rotation 𝛼3He of the detector, equivalent to this shift of 𝑠3He = 20.63(17) mm, amounts to 𝛼He3 = arctan 𝑠Gd-GEM 𝑑 × 180 ° 𝜋 = arctan 20.63(17) mm 332 mm × 180 ° 𝜋 = 3.56(3)° (A.7) Since the two detector measurements were taken under the same conditions using the same detector and sample rotation, the angles 𝛼 are as expected identical within errors. – 12 – 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 Acknowledgments This work was partially funded by the EU Horizon 2020 framework, BrightnESS project 676548. The authors would like to acknowledge and thank BNC for providing the use of their facilities and beam time on the ATHOS instrument. Equally, they would like to thank the ILL for providing use of their facilities and beam time on the D16 instruments. References [1] European spallation source ESS ERIC, http://europeanspallationsource.se/. [2] Institut Laue-Langevin, http://www.ill.eu/. [3] Spallation Neutron Source, http://neutrons.ornl.gov/facilities/SNS/. [4] J-PARC, http://j-parc.jp/index-e.html. [5] S. Peggs et al., ESS Technical Design Report, Tech. Rep. ESS-DOC-274 (2013). [6] K. Andersen, D. Argyriou, A. Jackson, J. Houston, P. Henry, P. Deen et al., The instrument suite of the european spallation source, Nucl. Instrum. Meth. A 957 (2020) 163402. [7] O. Kirstein et al., Neutron Position Sensitive Detectors for the ESS, PoS Vertex2014 (2014) 029 [arXiv:1411.6194]. [8] M. Markó, G. Nagy, G. Aprigliano and E. Oksanen, Neutron macromolecular crystallography at the european spallation source, in Methods in Enzymology, Elsevier (2020), pp. 125–151. [9] N. Niimura, Y. Karasawa, I. Tanaka, J. Miyahara, K. Takahashi, H. Saito et al., An imaging plate neutron detector, Nucl. Instrum. Meth. A 349 (1994) 521. [10] T. Hosoya, T. Nakamura, M. Katagiri, A. Birumachi, M. Ebine and K. Soyama, Development of a new detector and DAQ systems for iBIX, Nucl. Instrum. Meth. A 600 (2009) 217. [11] L. Coates, A.D. Stoica, C. Hoffmann, J. Richards and R. Cooper, The macromolecular neutron diffractometer (MaNDi) at the Spallation Neutron Source, Oak Ridge: enhanced optics design, high-resolution neutron detectors and simulated diffraction, J. Appl. Crystallogr. 43 (2010) 570. [12] M. Titov and L. Ropelewski, Micro-pattern gaseous detector technologies and RD51 Collaboration, Mod. Phys. Lett. A 28 (2013) 1340022. [13] B. Guerard, R. Hall-Wilton and F. Murtas, Prospects in MPGDs development for neutron detection, Tech. Rep. RD51-NOTE-2014-003 (2014). [14] L. Scharenberg et al., Resolving soft X-ray absorption in energy, space and time in gaseous detectors using the VMM3a ASIC and the SRS, Nucl. Instrum. Meth. A 977 (2020) 164310. [15] P. Thuiner, Ion space-charge effects in multi-GEM detectors: challenges and possible solutions for future applications, Ph.D. thesis, Vienna University of Technology, Vienna, 2016, CERN-THESIS-2016-199. [16] D. Pfeiffer et al., Rate-capability of the VMM3a front-end in the RD51 Scalable Readout System, Nucl. Instrum. Meth. A 1031 (2022) 166548 [arXiv:2109.10287]. [17] C. Höglund, J. Birch, K. Andersen, T. Bigault, J.-C. Buffet, J. Correa et al., B4C thin films for neutron detection, J. Appl. Phys. 111 (2012) 104908. [18] D. Pfeiffer et al., The 𝜇TPC method: improving the position resolution of neutron detectors based on MPGDs, 2015 JINST 10 P04004 [arXiv:1501.05022]. – 13 – http://europeanspallationsource.se/ http://www.ill.eu/ http://neutrons.ornl.gov/facilities/SNS/ http://j-parc.jp/index-e.html https://doi.org/10.1016/j.nima.2020.163402 https://doi.org/10.22323/1.227.0029 https://arxiv.org/abs/1411.6194 https://doi.org/10.1016/bs.mie.2020.01.005 https://doi.org/10.1016/0168-9002(94)91220-3 https://doi.org/10.1016/j.nima.2008.11.123 https://doi.org/10.1107/s0021889810008587 https://doi.org/10.1142/S0217732313400221 https://doi.org/10.1016/j.nima.2020.164310 https://cds.cern.ch/record/2238855 https://doi.org/10.1016/j.nima.2022.166548 https://arxiv.org/abs/2109.10287 https://doi.org/10.1063/1.4718573 https://doi.org/10.1088/1748-0221/10/04/P04004 https://arxiv.org/abs/1501.05022 2 0 2 3 J I N S T 1 8 P 0 4 0 2 3 [19] D. Pfeiffer et al., First Measurements with New High-Resolution Gadolinium-GEM Neutron Detectors, 2016 JINST 11 P05011 [arXiv:1510.02365]. [20] C. Altunbas et al., Construction, test and commissioning of the triple-GEM tracking detector for COMPASS, Nucl. Instrum. Meth. A 490 (2002) 177. [21] IFE, http://ife.no/. [22] D. Pfeiffer and P. Thuiner, Brightness deliverable, work package 4: Innovation of key neutronic technologies: Detectors and moderators, Tech. Rep. European Spallation Source ERIC (2015), https://doi.org/https://doi.org/10.17199/BRIGHTNESS.D4.3. [23] Budapest neutron centre, https://www.bnc.hu/. [24] G. Winter et al., DIALS: implementation and evaluation of a new integration package, Acta Crystallogr. D 74 (2018) 85. [25] A. Bressan, R. de Oliveira, A. Gandi, J.C. Labbe, L. Ropelewski, F. Sauli et al., Two-dimensional readout of GEM detectors, Nucl. Instrum. Meth. A 425 (1999) 254. [26] F. Sauli, GEM: A new concept for electron amplification in gas detectors, Nucl. Instrum. Meth. A 386 (1997) 531. [27] G. De Geronimo, J. Fried, S. Li, J. Metcalfe, N. Nambiar, E. Vernon et al., VMM1 — An ASIC for micropattern detectors, in Proceedings pf the IEEE Nuclear Science Symposium and Medical Imaging Conference and 19th Workshop on Room-Temperature Semiconductor X-ray and Gamma-ray Detectors, Anaheim, CA, U.S.A., 27 October–3 November 2012, pp. 633–639. [28] G. Iakovidis, VMM3a, an ASIC for tracking detectors, J. Phys. Conf. Ser. 1498 (2020) 012051. [29] M. Lupberger et al., Implementation of the VMM ASIC in the Scalable Readout System, Nucl. Instrum. Meth. A 903 (2018) 91. [30] S. Martoiu, H. Muller, A. Tarazona and J. Toledo, Development of the scalable readout system for micro-pattern gas detectors and other applications, 2013 JINST 8 C03015. [31] RD51 collaboration, Development of micro-pattern gas detectors technologies, https://indico.cern.ch/event/765096/contributions/3295721/. [32] A.A. Harms and G. McCormack, Isotopic conversion in gadolinium-exposure neutron imaging, Nucl. Instrum. Meth. 118 (1974) 583. [33] G. Török, V.T. Lebedev, A. Nagy, G.P. Gordeev and G. Zsigmond, Inelastic neutron scattering facilities at the budapest neutron center, Phys. Status Solidi C 1 (2004) 3182. [34] D16 instrument layout at ILL, https://www.ill.eu/users/instruments/instruments-list/d16/description/ instrument-layout. [35] V. Kelpšas, B. Lafumat, M.P. Blakeley, N. Coquelle, E. Oksanen and C. von Wachenfeldt, Perdeuteration, large crystal growth and neutron data collection of leishmania mexicana triose-phosphate isomerase E65Q variant, Acta Crystallogr. F 75 (2019) 260. [36] GEANT4 collaboration, GEANT4 — a simulation toolkit, Nucl. Instrum. Meth. A 506 (2003) 250. [37] K. Zeitelhack, personal communication. – 14 – https://doi.org/10.1088/1748-0221/11/05/P05011 https://arxiv.org/abs/1510.02365 https://doi.org/10.1016/S0168-9002(02)00910-5 http://ife.no/ https://doi.org/https://doi.org/10.17199/BRIGHTNESS.D4.3 https://www.bnc.hu/ https://doi.org/10.1107/S2059798317017235 https://doi.org/10.1107/S2059798317017235 https://doi.org/10.1016/S0168-9002(98)01405-3 https://doi.org/10.1016/S0168-9002(96)01172-2 https://doi.org/10.1016/S0168-9002(96)01172-2 https://doi.org/10.1109/NSSMIC.2012.6551184 https://doi.org/10.1088/1742-6596/1498/1/012051 https://doi.org/10.1016/j.nima.2018.06.046 https://doi.org/10.1016/j.nima.2018.06.046 https://doi.org/10.1088/1748-0221/8/03/C03015 https://indico.cern.ch/event/765096/contributions/3295721/ https://doi.org/10.1016/0029-554X(74)90669-7 https://doi.org/10.1002/pssc.200405419 https://www.ill.eu/users/instruments/instruments-list/d16/description/instrument-layout https://www.ill.eu/users/instruments/instruments-list/d16/description/instrument-layout https://doi.org/10.1107/s2053230x19001882 https://doi.org/10.1016/S0168-9002(03)01368-8 Introduction The NMX detector prototype Test beam in Budapest Test beam at ILL Grenoble Conclusion Annex: detailed geometrical calculations