利用GEANT4和FLUKA研究质子诱发散裂反应出射中子双微分截面

在中高能质子诱发散裂反应相关核工程设计中,可靠的蒙特卡罗模拟程序结合核反应理论模型具有较好的理论指导意义。本工作中,利用GEANT4耦合INCL4和ABLA理论模型以及FLUKA耦合PEANUT模型模拟计算了几百MeV至几个GeV质子轰击Be,Al,Fe,W,U等靶后30°,60°,120°,150°出射角产生的散裂中子双微分截面,并与现有实验数据进行了比较。结果发现,FLUKA和GEANT4模拟计算较好地再现了Al,Fe,W,U等靶实验测量数据。然而,模拟结果明显低估了Be靶出射中子能量小于10 MeV能区的实验数据。For the design of nuclear engineering related with medium-high energy proton induced spallation reaction,the reliable Monte Carlo simulation codes coupled with nuclear reaction models have a good theoretical guidance.In this work,the production spallation neutron double differential cross sections at 30°,60°,120°,150°emission angle for Be,Al,Fe,W,U target materials at incident proton energies between several hundred MeV and GeV are theoretically calculated by using the GEANT4 coupled INCL4 and ABLA,and the FLUKA coupled PEANUT.The calculated results were compared with the available experimental data.It is found that the GEANT4 and FLUKA calculations well reproduced the experimental measurement of Al,Fe,W,U target materials.However,calculations obviously underestimated the emission neutrons of Be target for lower than 10 MeV energy range.     

400~1500 MeV质子轰击铅靶和钨靶的出射中子能谱的FLUKA和Geant4模拟研究

在ADS散裂靶系统的优化设计中,蒙特卡罗方法结合可靠的散裂反应理论模型进行中子学计算具有重要的作用。本工作利用Geant4程序中的INCLXX模型、BIC模型以及BERT模型和FLUKA程序分别模拟了597 MeV和1 500 MeV质子轰击薄铅靶不同出射角度的中子双微分截面,500,1 500 MeV质子轰击厚铅靶不同出射角度的中子双微分产额,以及400,600,800,1 000和1 200 MeV质子轰击厚钨靶在反角方向（175 °）的中子双微分产额,并与实验数据进行比较。研究表明,对于薄铅靶,Geant4程序的INCLXX模型和FLUKA程序模拟结果与实验符合得更好。能量在10～40 MeV范围内,BIC模型模拟结果明显高于实验数据,而BERT模型模拟结果略微低于实验数据。对于厚铅靶,在40 MeV左右所有的模拟结果都低于实验数据。对于厚钨靶,Geant4程序的BIC模型和FLUKA程序与实验数据符合得较好,INCLXX模型在能量高于60 MeV时模拟结果低于实验数据,BERT模型与实验数据差异较大。总体来看,Geant4程序的INCLXX模型和FLUKA程序进行ADS散裂靶相关的中子学的计算是合理和可靠的。The reliable Monte Carlo simulation codes coupled with nuclear reaction models play an important role in the neutronic calculation for the design and optimization of the ADS spallation target. In this work, the double differential cross sections at different angles produced from a thin lead target bombarded with 597 and 1 500 MeV protons, the neutron energy spectra at different angles produced from a thick lead target bombarded with 500 and 1 500 MeV protons, and the neutron energy spectra in the backward direction(175°) produced from a thick tungsten target bombarded with 400, 600, 800, 1 000 and 1 200 MeV protons are calculated with the Geant4 code coupled INCLXX, BIC and BERT models and the FLUKA code. The calculations are compared with the available experimental data. The results show that, for the thin lead target, the calculations with the Geant4 coupled INCLXX model and FLUKA code well reproduce the experimental results. In a energy range from 10 to 40 MeV, BIC model obviously overestimates the experimental results, and BERT model slightly underestimates the experimental results. For the thick lead target, all of the calculations underestimate the experimental results around 40MeV. For the thick tungsten target, the Geant4 coupled BIC model and FLUKA code well reproduce the experimental results. INCLXX model underestimates the experimental results above 60 MeV. BERT model bad reproduces the experimental results. Overall, the neutronic calculations with the Geant4 code coupled INCLXX model and FLUKA code for the ADS spallation target is reasonable and reliable.     

The fragmentation of 20Ne at 400 A MeV

The total charge-changing cross sections and partial cross sections for the production of projectile fragments are measured in the interactions of 400 A MeV ²⁰Ne with aluminum, carbon and polyethylene targets sandwiched with CR-39 plastic nuclear track detectors. The measured total charge-changing cross sections are compared with the predictions using the Bradt-Peters semi-empirical formula, and the NUCFRAG2 and PHITS models. It is shown that the measured experimental results are in good agreement with the theoretical model prediction and other experimental results, and it can be clearly seen that the partial cross sections for fragment production show obvious odd-even effects.     

Theoretical analysis on fully differential cross sections for C~(6+)impact ionization of helium

Fully differential cross sections （FDCS） are calculated within a four-body model for single ionization of helium by C6＋ impact at the incident energy of 100 MeV/a.u. （atomic unit）. The results are compared with experimental data and other theoretical predictions. It is shown that our results are in very good agreement with experiment for three small momentum transfers in the scattering plane; however, some significant discrepancies are still present at the largest momentum transfer in both the scattering plane and the perpendicular plane. In actuality, the problem has not been explained by the theory during the last decade. Accordingly, the contributions of different scattering amplitudes to FDCS are analyzed. It is found that for the largest momentum transfer the cross section arising from a destructive interference of the three amplitudes is much smaller than the experimental data. However, the cross section due to the constructive interference of two scattering amplitudes between projectile-ionized electron interaction and projectile-passive electron interaction almost approaches the experimental data.     

Emissions of fragments and mesons produced in heavy-ion collisions at CSR and related energies

Using a multisource ideal gas (MSIG) model, we reconstruct the transverse emission source in the momentum space for light fragments produced in reactions ⁸⁶Kr-¹²⁴Sn at 25 MeV/nucleon and b=7—10 fm based on the theoretical predictions of the Isospin-Dependent Quantum Molecular Dynamics (IDQMD) model. We show that the MSIG model can reasonably describe the IDQMD-predicted results for the azimuthal distribution and the transverse momentum dependence of elliptic flow v₂ and fourth-order anisotropic flow v₄ but can only qualitatively describe the transverse momentum spectra. The azimuthal distributions of nuclear fragments produced in heavy-ion collisions at intermediate energies are studied by the MSIG model. The calculated results are compared and found to be in agreement with the experimental data of Ca-Ca, Nb-Nb, and Au-Au collisions at 150—800 MeV/nucleon beam energies. Meanwhile, the angular distributions of pions and kaons produced in heavy-ion collisions at the low-energy end (1—2 GeV/nucleon) of high energies are investigated by the MSIG model, too. The calculated results are compared and found to be in agreement with the experimental data of the KaoS Collaboration.     

Calculation of the wide-angle neutron spectra from the ~9Be(d,xn) reaction in a thick beryllium target

The multi-layer computing model is developed to calculate wide-angle neutron spectra, in the range from0° to 180° with a 5° step, produced by bombarding a thick beryllium target with deuterons. The double-differential cross-sections(DDCSs) for the ~9 Be(d, xn) reaction are calculated using the TALYS-1.8 code. They are in agreement with the experimental data, and are much better than the PHITS-JQMD/GEM results at 15°, 30°, 45° and 60° neutron emission angles for deuteron energy of 10.0 MeV. In the TALYS-1.8 code, neutron contributions from direct reactions(break-up, stripping and knock-out reactions) are controlled by adjustable parameters, which describe the basic characteristics of typical direct reactions and control the relative intensity and the position of the ridgy hillock at the tail of DDCSs. It is found that the typical calculated wide-angle neutron spectra for different neutron emission angles and neutron angular distributions agree quite well with the experimental data for 13.5 MeV deuterons. The multi-layer computing model can reproduce the experimental data reasonably well by optimizing the adjustable parameters in the TALYS-1.8 code. Given the good agreement with the experimental data, the multi-layer computing model could provide better predictions of wide-angle neutron energy spectra, neutron angular distributions and neutron yields for the ~9 Be(d, xn) reaction neutron source.     

Study on the selection of BLM detector for C-ADS Injector Ⅱ

The injector of C-ADS (Chinese Accelerator Driven Sub-critical System) project is a high current, fully super-conducting proton accelerator. Meanwhile, a BLM system is indispensable for this facility, especially in low energy segments. This paper presents some basic simulations for 10 MeV proton by Monte Carlo program FLUKA, as well as the distributions on different secondary particles in three aspects: angular, energy spectrum, and current. These results are beneficial to selecting the detector type and its location and determining its dynamic range matching different requirements for both fast and slow beam loss. Furthermore, in this paper the major impact of the background is also analyzed, such as superconducting cavity X radiation and radiation caused by material activation. This work is meaningful in BLM system research.     

A modified equation of state for Xe at high pressures by molecular dynamics simulation

The exact equation of state （EOS） for the fission gas Xe is necessary for the accurate prediction of the fission gas behavior in uranium dioxide nuclear fuel, However, the comparison with the experimental data indicates that the applicable pressure ranges of existing EOS for Xe published in the literature cannot cover the overpressure of the rim fission gas bubble at the typical UO2 fuel pellet rim structure. Based on the interatomic potential of Xe, the pressure-volume-temperature data are calculated by the molecular dynamics （MD） simulation. The results indicate that the data of MD simulation with Ross and McMahan＇s potential [M. Ross and A. K. McMahan 1980 Phys. Rev. B 21 1658] are in good agreement with the experimental data. A preferable EOS for Xe is proposed based on the MD simulation. The comparison with the MD simulation data shows that the proposed EOS can be applied at pressures up to 550 MPa and 3 GPa and temperatures 900 K and 1373 K respectively. The applicable pressure range of this EOS is wider than those of the other existing EOS for Xe published in the literature.     

Comparison of Proton Shower Developments in the BGO Calorimeter of the Dark Matter Particle Explorer between GEANT4 and FLUKA Simulations

The DArk Matter Particle Explorer(DAMPE) is a satellite-borne detector for high-energy cosmic rays and y-rays. To fully understand the detector performance and obtain reliable physical results, extensive simulations of the detector are necessary. The simulations are particularly important for the data analysis of cosmic ray nuclei, which relies closely on the hadronic and nuclear interactions of particles in the detector material. Widely adopted simulation softwares include the GEANT4 and FLUKA,...     

Study on gamma response function of E J301 organic liquid scintillator with GEANT4 and FLUKA

The gamma response function is required for energy calibration of EJ301 （5 cm in diameter and 20 cm in height） organic liquid scintillator detector by means of gamma sources. The GEANT4 and FLUKA Monte Carlo simulation packages were used to simulate the response function of the detector for standard 22Na, 60Co, 137Cs gamma sources. The simulated results showed a good agreement with experimental data by incorporating the energy resolution function to simulation codes. The energy resolution and the position of the maximum Compton electron energy were obtained by comparing measured light output distribution with simulated one. The energy resolution of the detector varied from 21.2% to 12.4% for electrons in the energy region from 0.341 MeV to 1.12 MeV. The accurate position of the maximum Compton electron energy was determined at the position 81% of maximum height of Compton edges distribution. In addition, the relation of the electron energy calibration and the effective neutron detection thresholds were described in detail. The present results indicated that both packages were suited for studying the gamma response function of EJ301 detector.     

Simulation of residual activity in steel and copper targets induced by 950 MeV/nucleon uranium ions

Radioactive nuclides will be produced in structures and surrounding materials of the relativistic heavy ion collider NICA. Predicting radioactivity in the materials, air, cooling water, and ground is an important part of the technical project. Some Monte Carlo transport codes (FLUKA, GEANT4, SHIELD) can simulate the production of radionuclides induced by heavy ions. The preliminary verification of MC codes with experimental data is necessary before using them for radiation protection purposes. In the present work, comparisons of the experiment and GEANT4 and SHIELD simulations of induced activity in thick stainless steel and copper targets irradiated with 950 MeV/nucleon uranium ions are presented.     

利用GEANT4和TALYS研究中子与铁作用的次级中子双微分截面

铁的次级中子双微分截面对核装置的设计、运行与维护具有重要作用。相关实验数据缺乏,且评价数据不完善,需要使用可靠的核理论模型进行计算。本工作利用GEANT4程序结合不同的中子评价数据计算了8.17, 11.5, 14.1和18 MeV等入射能量下中子轰击薄铁靶不同出射角度的次级中子双微分截面;同时利用TALYS程序和GEANT4程序结合BIC、BERT和INCLXX模型计算了25.7, 65, 100和150 MeV等入射能量下中子轰击薄铁靶不同出射角度的次级中子双微分截面,并与实验数据进行对比。研究表明,在20 MeV以下能区,ENDF/B-VIII.0库的计算结果与实验数据符合较好,BROND-3.1、CENDL-3.1、JENDL-4.0u和JEFF-3.3库的计算结果与实验数据存在差异。在20~150 MeV能区,GEANT4程序的BERT模型和TALYS程序的计算结果与实验数据符合较好,INCLXX模型和BIC模型的计算结果与实验数据存在分歧。整体来看,需要对铁的中子评价数据和核反应理论模型做进一步研究。     

Ep≤63 GeV 质子入射铅、铋引起的中子产生双微分截面的模拟(英文)

加速器驱动次临界系统(ADS) 液态Pb-Bi 散裂靶的设计中,需要可靠的理论计算工具精确地预言几个GeV 能量范围的质子引起的散裂反应产生的各种粒子和核素。利用蒙特卡罗模拟软件包Geant4 计算研究了800 MeV至3 GeV 质子入射铅、铋材料引起的中子产生双微分截面。比较了Geant4 不同物理模型得到的模拟结果与现有的实验数据。其中,Geant4 的QGSP BERT和QGSP INCL ABLA 物理模型模拟结果很好地再现了实验数据。本工作证实了Geant4 蒙特卡罗模拟软件包适合用于能量高达3 GeV 的质子入射铅、铋引起的中子产生双微分截面的模拟计算。A detailed design of the liquid Pb-Bi spallation target of the Accelerator Driven Systems (ADS) requires powerful and reliable computational tools that can accurately predict particles and nuclides production by the proton induced spallation reactions in the energy range of a few GeV. In this paper, the neutron production double-differential cross sections for Pb and Bi target materials at incident proton kinetic energies between 800 MeV and 3 GeV are studied by calculations with Monte Carlo simulation package Geant4. The simulated results of Geant4 with several physics models are compared with available experimental data. The simulated results generated by QGSP BERT and QGSP INCL ABLA physics models of Geant4 well reproduce the available experimental data. The present results validated that Geant4 Monte Carlo simulation package is suitable for simulations of neutron production double-differential cross sections of proton induced reaction on Pb and Bi targets in the incident energy range up to 3 GeV.     

Study on the selection of BLM detector for C-ADS Injector Ⅱ

The injector of C-ADS (Chinese Accelerator Driven Sub-critical System) project is a high current, fully super-conducting proton accelerator. Meanwhile, a BLM system is indispensable for this facility, especially in low energy segments. This paper presents some basic simulations for 10 MeV proton by Monte Carlo program FLUKA, as well as the distributions on different secondary particles in three aspects: angular, energy spectrum, and current. These results are beneficial to selecting the detector type and its location and determining its dynamic range matching different requirements for both fast and slow beam loss. Furthermore, in this paper the major impact of the background is also analyzed, such as superconducting cavity X radiation and radiation caused by material activation. This work is meaningful in BLM system research.     

金刚石探测器用于C-ADS注入器Ⅱ束损探测的模拟研究

加速器驱动次临界系统C-ADS 注入器Ⅱ采用强流超导质子直线加速器,设计流强达到10 mA。强流质子束产生的束流损失有可能损伤超导腔,需要专用的束流损失监测系统进行监测,束流损失探测器(BLM) 需要在高能量沉积导致超导腔失超之前提供警报。通过MCNPX 模拟计算10 MeV 质子在半波谐振腔(HWR)不同位置损失产生的辐射场,比较选取超导腔管道进出口处4 个位置为推荐束损探测器放置的位置,结合HWR腔结构和束损探测器选择的影响因素,计算了次级辐射在金刚石探测器中的能量沉积以及1° ~ 5°不同质子入射角度对探测的影响。结果表明,根据不同位置处探测器的能量沉积关系可以推断出束损点;不同入射角度不会影响生成粒子的能量分布,只轻微影响生成粒子的数目。The Chinese Accelerator Driven Subcritical System (C-ADS) injector II consists of super-conduction accelerating section which is half wave resonator (HWR), the designed beam intensity is 10 mA. To avoid the damage to the resonator due to proton beam loss, special Beam Loss Monitor (BLM) system is essential. BLM system could provide alarm signal when high energy deposition occurs which may cause the resonator quenching. Radiation field of 10 MeV proton lost at different point of the HWR are simulated with MCNPX, BLM could be set at proper positions based on the simulation. Considering the structure of HWR and the BLM detector selecting influence factor, radiation energy deposition in the diamond detector are simulated with MCNPX when the proton incidence angle change from 1°  5°, Possible beam loss point can be deduced from the relationship of energy deposition in detectors at different locations. The results indicate that energy spectra of secondary particles are independent with incidence angle; the number of secondary particles may be influenced slightly.     

Energy spectra and temperature distributions of clusters produced during ion-beam sputtering of metals

A method for evaluating the energy spectra and temperature dependences of the yield of neutral and charged clusters that consist of N≥5 atoms and are produced by ion bombardment of metals is proposed. The results are presented in the form of simple formulas. Theoretical energy spectra of clusters emitted as a result of bombarding niobium, tantalum, and iron targets by atomic ions of gold or xenon and temperature dependences of the yield of silver clusters produced by bombarding the targets with xenon ions are compared with experimental data.

     

Influence of Bonner sphere response functions above 20 MeV on unfolded neutron spectra and doses

Monte Carlo (MC) codes for neutron transport calculations such as MCNP, MCNPX, FLUKA, PHITS, and GEANT4, crucially rely on cross sections that describe the interaction of neutrons with nuclei. For neutron energies below 20 MeV, evaluated cross sections are available that are validated against experimental data. In contrast, for high energies (above 20 MeV) experimental data are scarce and, for this reason, every neutron transport code is based on theoretical nuclear models to describe interactions of neutrons with nuclei in matter. Here we report on the calculation of a complete set of response functions for a Bonner spheres spectrometer (BSS), by means of GEANT4 using the Bertini and Binary Intranuclear Cascade (INC) models for energies above 20 MeV. The recent results were compared with those calculated by MCNP/LAHET and MCNP/HADRON MC codes. It turns out that, whatever code used, the response functions were rather similar for neutron energies below 20 MeV, for all 16 detector/moderator combinations of the considered BSS system. For higher energies, however, differences of more than a factor of 2 were observed, depending on neutron energy, detector/moderator combination, MC code, and nuclear model used. These differences are discussed in terms of neutron fluence rates measured at the environmental research station (UFS), “Schneefernerhaus”, (Zugspitze mountain, Germany, 2650 m a.s.l.) for energies below 0.4 eV (thermal neutrons), between 0.4 eV and 100 keV (epithermal neutrons), between 100 keV and 20 MeV (evaporation neutrons), and above 20 MeV (cascade neutrons). In terms of total neutron fluence rates, relative differences of up to 4% were obtained when compared to the standard MCNP/LAHET results, while in terms of total ambient dose equivalent, relative differences of up to 8% were obtained. Both the GEANT4 Binary INC and Bertini INC gave somewhat larger fluence and dose rates in the epithermal region. Most relevant for dose, however, those response functions calculated with the GEANT4 Bertini INC model provided about 18% less neutrons in the cascade region, which led to a roughly 13% smaller contribution of these neutrons to ambient dose equivalent. It is concluded that doses from secondary neutrons from cosmic radiation as deduced from BSS measurements are uncertain by about 10%, simply because of some differences in nuclear models used by various neutron transport codes.