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.     

50MeV/u 12C离子实验靶区出射中子角分布实验测量

本文用氟和碳阈探测器对50MeV/u ¹²C离子引起的重离子反应实验靶区E_n>11MeV和E_n>20MeV中子的角分布进行了实验测量,结果显示出明显的前冲分布,较高能量的中子具有更强的前冲趋势.     

2.8MeV中子诱发238U裂变的瞬发中子谱高能段实验测量

在中国原子能科学研究院的高压倍加器上利用中子飞行时间方法测量了2.8 MeV中子引起238U裂变的瞬发中子能谱,通过增大样品量和设计合适的屏蔽体提高了效应/本底比,使得测量数据的不确定度满足预期目标（在5.5～14 MeV能区内,能量间隔为0.5 MeV条件下能谱的不确定度小于10%）。将实验几何、中子源能量分布及角分布、探测效率、束流的时间结构等实验条件输入到MCNP程序里,模拟了出射的中子谱,模拟结果与测量结果在不确定度范围内一致,验证了在入射中子能量较低时238>U的裂变瞬发中子谱评价数据的可靠性。The Prompt Fission Neutron Spectrum (PFNS) of 238U induced by 2.8 MeV neutron was measured using the Cockcroft-Walton accelerator in China Institute of Atomic Energy (CIAE). The signal/background ratio was improved by increasing the amount of sample mass and using an appropriate shielding system. The final uncertainty of neutron energy spectrum in 5.5～14 MeV region is less than 10% with a bin size of 0.5 MeV which has reached this project's anticipation. The experimental geometry, the angular distribution and energy distribution of neutron source, the detection efficiency and time structure of deuteron beam were inputted into the MCNP code to simulate the outgoing neutron spectrum. The simulated results agree with the experimental ones within the uncertainty. The result indicates that the evaluated PFNS for 238U at low neutron energy is reliable.     

25MeV/u 6He在9Be靶上的弹性转移反应研究

报告了25MeV/u ⁶He在⁹Be靶上的弹性散射和弹性转移的实验结果. 用光学模型和一阶DWBA的方法计算了弹性散射和³He集团转移的截面, 计算结果和实验测量基本一致. 计算表明³He在⁹Be中的谱幅度应该明显大于壳模型的理论计算值0.70.     

400 MeV/u碳离子打靶的屏蔽参数计算

为了给医用重离子加速器提供一种专用的快速计算屏蔽厚度或对蒙特卡罗模拟结果进行验证的方法.采用FLUKA程序完成了400 MeV/u碳离子打不同靶的屏蔽参数计算.首先研究了打厚靶(铁、水)产生的次级辐射场的角度分布及主要成分;接着给出了不同角度范围下周围剂量当量在屏蔽体中的衰减曲线,通过拟合数据进一步得到不同角度范围下的...     

80MeV/u 20Ne轰击9Be碎裂产物反应总截面的测量

介绍了用透射法测量中能区²⁰Ne打⁹Be靶碎裂产生的次级束与Si靶作用的核反应总截面的方法,以及实验的探测器布局、实验过程和实验结果.并对理论上预言有奇异结构的核¹²N,¹⁷Ne和¹⁷F的实验结果与其相邻核进行了比较.     

28MeV/u 18O与全阻止厚天然铅靶反应中丰中子汞同位素的截面测量

本次实验利用28MeV/u的¹⁸O束流和一套在线熔化铅靶装置以及4πΔEβ—γ符合技术,测量了该中能重离子反应中丰中子汞同位素^205—208Hg生成的相对截面比以及绝对截面的近似值.并外推出了²⁰⁹Hg的可能生成截面.发现由靶核²⁰⁸Pb削裂两个质子添加一个中子而生成²⁰⁷Hg的截面有一个突然的比较大的下降,这可以由反应的Q值的变化规律得到解释.     

晕核6He（25 MeV/u）的单、双中子转移截面测量

实验测量了25 MeV/u 的6He+9Be反应中的单、双中子转移微分截面. 利用反应耦合道模型初步分析了实验数据, 计算大体重现截面值. 转移反应对末态有相当的选择性, 并且双中子转移截面主要来自到基态的一步过程. 后角度截面的上升可能与连续态的耦合有关. Differential cross sections for quasi elastic scattering at forward angles and for 1n and 2n transfer reaction at backward angles induced by 6He at 25 MeV/u from 9Be target were measured. The experimental data were analyzed in the framework of the CRC model. The raise at backward angle of the 1n and 2n transfer cross section can not reproduced by the current calculation.     

贫铀球壳中D-T中子诱发的铀反应率的测量与分析

为校验次临界能源堆的概念设计,在R19.4/30.0 cm的贫铀球壳装置上采用活化法开展14 MeV中子学积分实验.布放6片贫铀活化片于球壳中与入射D离子束90°方向上的不同位置处活化,用HPGe探测器测量238U(n,γ)反应、238U(n,f)及235U(n,f)反应和238U(n,2n)各反应产物发射的特征γ射线,得到了相应的反应率.238U(n,γ)反应率的不确定度为3.6％-3.7％,238U(n,D和235U(n,f)反应率的不确定度为5.1％-5.9％,238U(n,2n)反应率的不确定为4.3％-4.7％.用MCNP5程序在ENDF66c数据库下进行模拟计算,238U(n,γ)反应率的计算值/实验值(C/E)为0.972-1.034,238U(n,f)和235U(n,f)反应率的C/E为0.983-1.058,238U(n,2n)反应率的C/E为0.979-1.019.     

120~430 MeV/u 碳离子的微剂量学特性研究(英文)

单粒子微剂量谱在放射治疗中是一个极其重要的参数,它可以用来评估辐射场的生物学效应。利用蒙特卡洛程序FLUKA模拟计算了由碳离子产生的混合辐射场能量沉积的微观模式。从已公开发表的文献中选取了实验测量300 MeV/u 碳离子的线能能谱,并与相同物理条件下模拟计算得到的线能能谱相比较,结果吻合得很好。此外,还计算了120~430 MeV/u 的碳离子的剂量平均线能能谱、频率平均线能和剂量平均线能。所得到的频率平均线能值为185~ 28.3 keV/m而剂量平均线能值则为272~ 64.1 keV/m。本文的结果对于制定碳离子放射治疗的治疗计划有着重要的意义.Microdosimetric single event spectrum is a significant parameter in radiotherapy, which can be used to evaluate the radiation biological effect. In this paper, microscopic patterns of energy deposition are simulated with Monte Carlo code FLUKA at mixed radiation fields during carbon ions therapy. The results are compared with experimental measured results at 300 MeV/u carbon ion and good agreement has been found. Meanwhile, dose-weighted lineal energy spectra, frequency averaged lineal energy values and dose averaged lineal energy values of carbon ion with energy from 120 to 430 MeV/u were calculated,too. The frequency averaged lineal energy values are from 185 to 28.3 keV/m while the dose averaged lineal energy values are from 272 to 64.1 keV/m. These studies are useful for treatment plan in carbon ion radiotherapy.     

PGNAA煤质在线测量中同时提高C、O和其它元素测量精度的慢化体模拟

提出两种慢化体设计方案,第一种方案直接在原有慢化体中引入孔道,第二种方案在慢化体中增加一层铅之后再引入孔道.利用蒙特卡罗方法对两种方案进行研究.计算结果表明,采用第一种方案可以有效地提升C、O元素的测量精度,但同时会降低H、Si等元素的测量精度;采用第二种方案可提高C元素和O元素的测量精度,同时也可以提高H、Si等其它元素的测量精度.     

Studies of K-shell x-ray energy shifts induced by MeV/u heavy ions

This paper reports that the K x-ray spectra of the thin target 47Ag, 48Cd, 49In and 50Sn were measured by an HPGe semi-conductor detector in collisions with 84.5 MeV 6C4+ ions. Our experiment revealed the Kα x-ray energy shifts were not obvious and the Kβ1 x-ray energy shifts were about 90～110 eV. The simple model of Burch et al has been previously used to calculate the K x-ray energy shifts due to an additional vacancy in 2p orbit. The present work extends the model of Burch to calculate the x-ray energy shifts of multiple ionized atoms induced by heavy ions with kinetic energy of MeV/u. In addition to our experimental results, many other experimental results are compared with the calculated values by using the model.     

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.     

超短超强激光与铜靶相互作用产生Kα源的蒙特卡罗模拟

超短超强激光打靶产生的超热电子，与固体靶相互作用时会产生Kα线辐射.由经典定标律给出了法线方向超热电子的温度.利用蒙特卡罗方法，对超热电子在固体靶中的传输进行了研究，模拟了不同靶厚度情况下Kα产额和角分布及不同电子温度下Kα光子的转化效率.计算结果与实验符合较好.结果表明：在一定电子温度下，随着靶厚度的增加Kα光子产额会达到饱和，并会使Kα光子发射的各向异性变得更加严重；存在最佳的电子温度，使Kα线转化效率最高. 关键词： 超短超强激光 超热电子 蒙特卡罗方法 Kα线')" href="#">Kα线     

Particle in Cell/Monte Carlo Collision Method for Simulation of RF Glow Discharges: Effect of Super Particle Weighting

A parallel Particle in Cell/Monte Carlo Collision (PIC/MCC) numerical code for glow discharge plasma simulations is developed and verified. This method is based on simultaneous solution of the Lorentz equations of motion of super particles, coupled with the Poisson's equation for electric field. Collisions between the particles are modelled by the Monte Carlo method. Proper choice of particle weighting is critically important in order to perform adequate and efficient PIC simulations of plasma. Herein, effects of particle weighting on the simulations of capacitive radio‐frequency argon plasma discharges are studied in details. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Isotopic production cross sections of fragmentation residues produced by 18O ions on a carbon target near 260 MeV/nucleon

The isotopic cross sections of residual nuclei produced in fragmentation reactions of ¹⁸O projectiles impinging on a carbon target at energies near 260 MeV/nucleon were measured at the HIRFL facility in Lanzhou (China). A full identification of atomic and mass numbers of fragments was achieved from the determination of their magnetic rigidity, energy loss, and time of flight. The production cross sections for a dozen of nitrogen, carbon, and boron isotopes were determined with uncertainties below 30% for most of the cases. The obtained cross sections for N and B isotopes show a rather good agreement with previous experimental data obtained with different projectile energies. The cross sections for some C isotopes seem to exhibit a dependence on the projectile energy. A comparison of the data and several theoretical model calculations are presented.     

Construction of modified embedded atom method potentials for Cu, Pt and Cu-Pt and modelling surface segregation in Cu3Pt alloys

In this work, surface segregation to Cu₃Pt surfaces is studied with the modified embedded atom method (MEAM). This work is triggered by the catalytic importance of Cu-Pt alloys, together with the contradictory experimental results for the surface segregation in Cu₃Pt(1 1 1) alloys based on low energy ion scattering (LEIS) [Y.G. Shen, D.J. O’Connor, K. Wandelt, R.J. MacDonald, Surf. Sci. 328 (1995) 21] and low energy electron diffraction (LEED) [Y. Gauthier, A. Senhaji, B. Legrand, G. Tréglia, C. Becker, K. Wandelt, Surf. Sci. 527 (2003) 71]. In order to accurately describe the segregation behaviour in the Cu₃Pt system, a reliable potential, that is also applicable to surface phenomena, is indispensable. Therefore, first, new MEAM parameters are derived, consistently based on ab initio density functional theory (DFT) calculations, according to a method that is a modification of previous work [P. van Beurden, G.J. Kramer, Phys. Rev. B 63 (2001) 165106]. Upon testing, these parameters prove to reproduce very well various surface properties of this system. Next, Monte Carlo (MC) simulations combined with the newly derived MEAM potentials are set up to investigate surface segregation to low index single crystal surfaces. For the Cu₃Pt(1 1 1) surface, our MC/MEAM simulations agree completely with the available LEIS evidence and contradict the unusual depth profile based on LEED. However, the slight Pt enrichment observed in the LEED experiments can be reproduced by assuming a slight Pt excess in the bulk of the sample. The simulated composition depth profile, on the other hand, does not agree with the LEED evidence. Also, for the Cu₃Pt(1 0 0) surface, the MC/MEAM results agree completely with LEIS experiments. For the Cu₃Pt(1 1 0) surface, finally, the MC/MEAM simulations show a somewhat deviating behaviour with respect to the experimental LEIS evidence. The possibility of a missing-row reconstruction is evaluated, but cannot explain the discrepancy for the Cu₃Pt(1 1 0) system. In order to further investigate the deviation from the experiments, additional DFT and MEAM calculations are performed in search of the preferred surface termination for Cu₃Pt(1 1 0). Both DFT and MEAM calculations agree on the pure Cu layer as the most stable surface termination. Although the experiment was extensively tested for reproducibility, it possibly reflects a metastable state. Finally, in view of the importance of small and less orderly particles in catalysis, the newly derived MEAM parameters are used in order to study the segregation to Cu₃Pt vicinal surfaces with {1 1 1} terraces, for which no experimental information is available yet.