<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>October</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>SIMULATION OF SPECTRA OF CYLINDRICAL NEUTRON COUNTERS USING THE GEANT-4</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>A.V. Churakov</string-name>
          <email>churakov@nf.jinr.ru</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>A.K. Kurilkin</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>J. Saiz</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Andrey Churakov</institution>
          ,
          <addr-line>Alexei Kurilkin</addr-line>
          ,
          <country>Juan Saiz Lomas</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Dubna State University</institution>
          ,
          <addr-line>Dubna</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>Joint Institute for Nuclear Research</institution>
          ,
          <addr-line>Dubna</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff3">
          <label>3</label>
          <institution>University of York</institution>
          ,
          <addr-line>York, England</addr-line>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2019</year>
      </pub-date>
      <volume>4</volume>
      <issue>2019</issue>
      <fpage>408</fpage>
      <lpage>412</lpage>
      <abstract>
        <p>It is widely known that the amplitude spectrum of the helium proportional counter produced by the irradiation of thermal and cold neutrons has a peak of full absorption with energy of 768 keV and two small “shelves”, caused by boundary effects, i.e. the absorption of charged particles (proton or tritium nuclei) in the detector wall. The simulation of the amplitude spectra of cylindrical counters with different gas filling is presented in the paper. The possibility of the third peak, not coinciding with that of full absorption is shown, whilst the peak position depends on the ratio of the path length towards the counter diameter. The results obtained may be of interest for the development of low efficiency neutron detectors and neutron monitors.</p>
      </abstract>
      <kwd-group>
        <kwd>neutron</kwd>
        <kwd>detector</kwd>
        <kwd>proportional counter</kwd>
        <kwd>simulation</kwd>
        <kwd>GEANT4</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>The 3He-filled Cylindrical proportional counters are perhaps the most widely used neutron
detectors. They are simple, reliable and commercially available devices with high detection efficiency
of thermal neutrons and very low sensitivity to gamma rays. Proportional counters have a good enough
time resolution which allows their use on diffractometers operated in time-of-flight mode.</p>
      <p>An important characteristic of the proportional counters is the amplitude spectrum. It is known
that the amplitude spectrum of the helium proportional counter irradiated by thermal or cold neutrons
has a peak of full absorption with an energy of 768 KeV and two small "steps". These “wall effects”
arise from the hit of primary ionizing charged particles (proton or tritium nuclei) in the detector wall.
(fig. 1). When one of these particles products collides with the wall of the detector, its energy is
dissipated and does not contribute to the full energy peak.</p>
      <p>Gamma-rays</p>
    </sec>
    <sec id="sec-2">
      <title>2. Simulation assumption and tools</title>
      <p>
        In our work, we carried out a simulation of the amplitude spectrum of the proportional
counter. GEANT 4 [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ] was used as a modeling tool. It was designed to simulate the tracks of
elementary particles passing through matter by the Monte Carlo method. The version of GEANT4
used in our simulations is 10.03.03. In addition, the software packages ROOT [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ] on its version
6.14.06, for the storage and visualization of the data, and the program SRIM 2013 [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ], for calculations
involving particles of primary ionization, were also used.
      </p>
      <p>
        The simulated cylindrical counter was irradiated with a uniform flow of thermal neutrons. The
particle generator emulates the thermal neutron beam of energy (0.0253 eV). Neutrons pass through
the geometric model of the detector, taking into account their interaction with the walls of the detector.
The wall thickness was considered to be 0.2 mm, which corresponds to typical counter wall
thicknesses. The rate of neutrons was considered low enough to neglect double and triple events. In
the simulation, it was assumed that each interaction between a neutron and a 3He nucleus happening
in the sensitive volume of the detector, will be registered, and the induced charge proportional to the
energy losses of particles in the volume of the counter will be collected at the anode. When the
primary particles are completely absorbed in the volume of the counter, the fluctuation of the charge
induced on the anode has a normal distribution. The same assumption was believed to be true for all
other energy values of primary ionizing particles. The value of the energy resolution of the detector
determines the dispersion of the induced charge distribution. On the other hand, the modeling did not
take into account the part of the spectrum caused by gamma radiation, since the level of the gamma
background is determined mainly by the neutron source [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] and the proportional counters themselves
have a sufficiently low sensitivity to gamma radiation.
      </p>
    </sec>
    <sec id="sec-3">
      <title>3. Results and discussion</title>
      <p>
        We started by modelling the spectrum of a standard proportional counter. In many scientific
facilities, counters of the SNM-18 type and their numerous analogues have been used. For comparison
we chose the “Helium-18-200” counter, produced by JSC "Consensus" [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ], which is similar to the
SNM-18. On the left side of figure 2 the measured spectrum for these counters is presented. On the
right hand side, the simulation of the neutron component of the spectrum is shown. The simulation
results are in excellent agreement with the measured spectrum (fig. 2).
4500
4000
3500
3000
ts 2500
n
ou 2000
C
1500
1000
500
0
0 100 200 300 400 500 600 700 800 900 1000
      </p>
      <p>Energy (keV)</p>
      <p>It is necessary to say a few words about the efficiency of cylindrical detectors. Counters of
SNM-18 type are believed to absorb almost all thermal neutrons. However, this is not quite true. The
maximum efficiency of such counters is 86% in a narrow area strictly in the middle of the detector,
and the average efficiency of the detector, because of the cylindrical shape, will be 74 %. In
largediameter detectors, the bulk of neutrons are counted near the input window, so the effect of the
detector's cylindricity is smaller. However, the amount of gas required to fill the detector increases in
proportion to the square of the diameter, which causes a significant increase in the cost of the
experimental installation. It would be more rational to cover the area of interest with a double layer of
neutron counters, saving on the quite expensive gas.</p>
      <p>We simulated changes in amplitude spectra of SNM18 - type counters caused by the gas
mixtures. Helium is a fluid gas, so in some years, depending on the quality of the detectors, the gas
mixture may leak out of the sealed housing, and the efficiency and amplitude spectrum of the counter
will change. The results are shown at (fig. 3), where the pressure is expressed as a percentage of the
initial pressure value, 9 bar. The initial composition of the gas mixture is 8 bar 3He+ 0,95 bar Ar +0,05
bar CO2
100%
90%
80%
70%
60%
50%
40%
30%
20%
10%
5%
1%
100%
90%
80%
70%
60%
y
iicne 50%
c
ffE 40%
30%
20%
10%
0%</p>
      <p>At 100% to 30% of the initial pressure, the spectra look like the typical counter spectrum
(fig.1), and the decrease in pressure by half does not lead to any noticeable changes in the amplitude
spectrum. At 20% of the original pressure and below, the spectra are not similar to the classic, familiar
spectrum of the cylindrical counter. On the 20% spectrum we see the peak of full gathering, a strong
peak from the triton and a new peak arises, superimposed to it. On subsequent spectra, this new peak
almost completely suppresses the peak of full absorption.</p>
      <p>The explanation for this new peak is the following one; primary ionization energy is spent
mainly in the ionization of the gas mix and losses in the counter walls. The amount of primary
ionization charge collected at the detector anode depends on the number of electron-ion pairs
remaining in the detector volume. When the pressure decreases, the mileage of charged particles
increases. Calculations were carried out using the program SRIM 2013. At a pressure value of about
30% of the initial total, the particle range is of 16.97 mm, which is comparable to the diameter of the
detector. Further pressure reduction leads to an increase in the length of the track, and it no longer fits
completely into the volume of the detector. It follows from the simulations that this is the most likely
value of the track, that depends mainly on the ratio between the diameter of the detector and the
maximum length of the track.</p>
      <p>Amplitude spectra for different detector diameters are presented. The gas filling used was the
standard one. Up to about 5 mm, we did not observe the distortion of the spectrum and the appearance
of the third peak, since the mean free path of particles in the gas mixture 8 bar 3He + 0.95 bar Ar
+0.05 bar CO2 is about 5.1 mm, but at smaller diameters, the changes in the spectrum and the third
peak are clearly visible.
SNM-18
Max.Eff
1000
10
1
0
100
200
300
400
500
600
700
800
Energy (keV)</p>
    </sec>
    <sec id="sec-4">
      <title>4. Conclusion</title>
      <p>Amplitude spectra of 3He-filled proportional counters were simulated (Fig. 4). The appearance
of an additional peak (Landau peak) is shown. This peak is caused by the incomplete stacking of the
tracks of primary ionization particles in the counter volume. The position of the peak depends on the
ratio between the diameter of the counter and the range of the primary ionization particles. The results
of this work can be used in the development of some types of new neutron detectors, such as neutron
detectors with small gas pressure: monitors, low-background cold neutron detectors, neutron
calorimeters («neutron telescopes»). The appearance of the new peak under the low gas pressure must
be considered in the developments of thin position-sensitive proportional counter tubes with resistive
readout and during calibration of slow neutron counters.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>J.</given-names>
            <surname>Allison</surname>
          </string-name>
          , et all., // Geant4 - simulation toolkit // Nuclear Instruments and
          <source>Methods in Physics Research</source>
          , vol.
          <source>A 506</source>
          , pp
          <fpage>250</fpage>
          -
          <lpage>303</lpage>
          ,
          <year>2003</year>
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>Rene</given-names>
            <surname>Brun</surname>
          </string-name>
          and
          <article-title>Fons Rademakers, ROOT - An Object Oriented Data Analysis Framework</article-title>
          ,
          <source>Proceedings AIHENP'96 Workshop</source>
          , Lausanne, Sep.
          <year>1996</year>
          , Nucl.
          <source>Inst. &amp; Meth. in Phys. Res. A</source>
          <volume>389</volume>
          (
          <year>1997</year>
          )
          <fpage>81</fpage>
          -
          <lpage>86</lpage>
          . See also http://root.cern.ch/
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>James</surname>
            <given-names>F.</given-names>
          </string-name>
          <string-name>
            <surname>Ziegler</surname>
            , Jochen P. Biersack,
            <given-names>Matthias D. Ziegler SRIM</given-names>
          </string-name>
          .
          <source>The Stopping and Range of Ions in Matter</source>
          ,
          <year>2015</year>
          , Lulu Press, NY, USA. See also http://www.srim.org/
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>D.H.</given-names>
            <surname>Beddingfield</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.H.</given-names>
            <surname>Johnson</surname>
          </string-name>
          ,
          <string-name>
            <surname>H.O. Menlove.</surname>
          </string-name>
          <article-title>3He neutron proportional counter performance in high gamma-ray dose environments // Nucl</article-title>
          . Instr. and Meth.
          <source>in Phys. Res. A</source>
          ,
          <year>2000</year>
          , V.455 pp.
          <fpage>670</fpage>
          -
          <lpage>682</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5] htttp://consensus-group.
          <source>ru (accessed 28.09</source>
          .
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>