Ionization for Hf- and W-L-shell by electron impact

Electron-induced Hf-, W-L-shell partial, total production cross sections, mean ionization cross sections and Hf-L_3-shell ionization cross sections (at two energies) have been measured as functions of electron energies (from near threshold to 36keV). The influence of electrons reflected from the backing of the thin targets on measured results was corrected using a model to relate to the electron transport process. The mean paths of electron multi-scattered in the target itself (including forward and backward scattering) were calculated by means of Monte Carlo program (EGS4) and they were used to correct measured results. A comparison with both theoretical predictions was given.     

Incomplete fusion reactions in 16O+165Ho

Excitation functions for evaporation residues of the system ¹⁶O + ¹⁶⁵Ho have been measured up to 100 MeV. Recoil range distribution of long lived reaction products were measured at ¹⁶O beam energy of 100 MeV. Detailed Monte Carlo simulation of recoil range distributions of products were performed with the help of PACE2 code, in order to extract the contributions of incomplete fusion in the individual channels. The results clearly show the incomplete fusion contributions in the tantalum and thulium products. This is confirmed by the predictions of breakup fusion model of the incomplete fusion.     

TREND code for reconstruction of neutron reaction cross sections and calculation of the Doppler broadening in the region of resolved resonances

The TREND code is developed for reconstruction of neutron cross sections in the region of resolved resonances from the data of the international evaluated nuclear data libraries using the Reich-Moore, Adler-Adler, and Breit-Wigner formalisms and for calculation of the Doppler broadening of the resulting cross sections within the classical approximation. The TREND code is incorporated into the MCU code package for Monte Carlo reactor calculations. The TREND and MCU codes have been used to evaluate the data of the ROSFOND data bank on uranium and plutonium isotopes in comparison with the data of the ICSBEP data bank.     

Measurements of total fast neutron cross sections for aluminum and barium by application of Si(Li)-detectors

Total neutron cross sections of²⁷Al and natural barium have been determined at bombarding energies in the range between 5.5 MeV and 9.6 MeV by transmission measurements. Si(Li)-detectors were used for neutron spectroscopy via charged particle reactions in the active volumes of the detectors. Inscattering corrections have been made by means of a Monte Carlo computer program. In most cases the measured cross sections have root-mean-square uncertainties of 1% to 3%. Cross sections of barium not yet measured have been calculated using optical model potential parameters obtained from an analysis of the experimental results.     

Model differential cross section of elastic scattering of electrons by atoms in the Monte Carlo simulation of the passage of electrons in a substance

A new model differential cross section is proposed for describing elastic scattering of electrons in simulating the passage of electrons through a substance by the Monte Carlo method. This differential cross section correctly describes the first and second transport scattering cross sections, but is characterized by the total elastic scattering cross section much smaller than the actual value of the total elastic scattering cross section. The application of this differential cross section makes it possible to considerably reduce the number of elastic collisions in the Monte Carlo simulation of passage of electrons in a substance and to model the passage of high-energy electrons using the individual collision model.     

Monte Carlo calculation of CR-39 efficiency for fast neutron detection using a combination of MCNP and SRIM codes,and comparison with experimental results

The elastic scattering of light nuclei, especially hydrogen, is widely used for detection of fast neutrons. Fast neutrons have a few interactions in solid state nuclear track detectors (SSNTDs), such as CR-39. It has been shown that to increase the efficiency of detector, a radiator layer, e.g. polyethylene may be used. In this work, the optimum thickness of the radiator layer has been determined. Some different thicknesses of aluminium layers have been used as discriminators for which the efficiency of detector to Am–Be neutron source has been calculated. A program FTRAC has been developed for the calculation of the detection efficiency. This software determines the specifications of recoil protons (using PTRAC card of MCNP code) and also the kinematics of interactions. Transport of recoil protons has been carried out using SRIM Monte Carlo code. A good agreement have been shown between measurement and simulation results.     

Six-fermion calculation of intermediate-mass Higgs boson production at future {e^+ e^-} colliders

The production of an intermediate-mass Higgs boson in processes of the kind fermions at the energies of future linear colliders is studied. The recently developed and fully automatic algorithm/code ALPHA is used to compute the tree-level scattering amplitudes for the reactions . The code has been interfaced with the Monte Carlo program HIGGSPV/WWGENPV, properly adapted to 6-fermion production, in order to provide realistic results, both in the form of cross sections and event samples at the partonic level. Phenomenological results, that incorporate the effects of initial-state radiation and beamstrahlung, are shown and commented, emphasizing the potentials of full six-fermion calculations for precise background evaluation as well as for detailed studies of the fundamental properties of the Higgs particle. Received: 22 May 1997 / Published online: 20 February 1998     

A Langevin description of the competition between fusion and deep-inelastic collisions close to the barrier

A Langevin Monte Carlo method based on the surface friction model is used to investigate the competition between fusion and deep-inelastic collisions close to the barrier. Data for the⁵⁸Ni+^112,124Sn systems are analyzed. The measured excitation functions for fusion are well reproduced, whereas the calculated deep-inelastic cross sections deviate somewhat from experiment. Predictions for the corresponding spin distributions are made. Contrary to what is usually assumed they are not well separated in 1-space.     

Extension of the Monte-Carlo code MOCADI to fusion-evaporation reactions

The Monte-Carlo code MOCADI has been extended, besides to fission and fragmentation reactions, also to fusion-evaporation reactions. The updated version includes the cross sections and the reaction kinematics for the production of evaporation residues calculated according to the statistical model code PACE2, while the atomic interaction with matter can be treated either with the code ATIMA or with the code TRIM. The present version of MOCADI can now be employed world-wide for design studies of new low-energy nuclear physics facilities and for the optimization of experiments involving fusion-evaporation reactions, e.g. mass measurements, laser spectroscopy, super-heavy synthesis.     

Effect of entrance channel parameters on the fusion of two heavy ions: Excitation functions of reaction products in16O +66Zn and37Cl +45Sc reactions

Excitation functions of reaction products formed in¹⁶O +⁶⁶Zn and³⁷Cl +⁴⁵Sc systems, leading to the same compound nucleus,⁸²Sr, were measured using recoil-catcher technique and off-line γ-ray spectrometry. The contribution of non-compound processes like transfer and incomplete fusion (ICF) reactions to the cross-sections of different evaporation residues were delineated by comparing the experimental data with the predictions of Monte Carlo simulation code PACE2. The results show that non-compound processes become a significant fraction of the total reaction cross-section in¹⁶O +⁶⁶ systems in the beam energy range studied, while³⁷Cl +⁴⁵Sc gives mainly compound nucleus products. The mass asymmetry dependence of the fusion and non-compound cross-sections have been analysed in terms of the static fusion model and sum rule model     

混合评价核数据库HENDL1.0/MG/MC研制

根据世界几个主要核评价数据库,如ENDF/B 6(美国)、JEF 2.2(欧盟)、JENDL 3.2(日本)、BROND 2.2(俄罗斯)、CENDL 2.1(中国)和FENDL 2(IAEA/NDS),兼顾聚变、裂变以及聚变 裂变次临界混合堆设计研究的多种需要,经过甄别、筛选,最后集成为包含213个核素的基本评价文件,名为HENDL1.0/E的核评价数据库.在此基础上,利用目前流行的群常数加工程序系统NJOY和输运截面制备程序TRANSX制作两套用于中子或/和光子输运计算的输运截面工作库:①参考Vitamin J能群结构制作了175群中子和42群光子、中子 光子耦合多群工作数据库HENDL1.0/MG,可用于离散纵标Sn法程序计算;②连续能群结构、紧凑ENDF(ACE)格式中子截面库HENDL1.0/MC,可用于蒙特卡罗方法输运计算,如MCNP.另外还制作了可用于燃耗(嬗变)计算的燃耗库BURNUP.LIB和响应函数库RESPONSE.LIB两个专用数据库.同时,也对HENDL1.0综合评价核数据库的有效性进行了抽样测试、基准检验和初步确认. A Hybrid Evaluated Nuclear Data Library(HENDL) named as HENDL1.0 has been developed by Fusion Design Study (FDS) team of Institute of Plasma Physics, Academia Sinica (ASIPP) to take into account the requirements in design and research relevant to fusion, fission and fusion-fission sub-critical hybrid reactor. HENDL1.0 contains one basic evaluated sub-library naming HENDL1.0/E and two processed working sub-libraries naming HENDL1.0/MG and HENDL1.0/MC, respectively. Through carefully comparing...     

考虑中子泄漏修正的快堆蒙特卡罗少群截面参数制作与验证方法

一种基于B1均匀化方程的泄漏修正模型在连续能量蒙特卡罗程序TRIPOLI4中得以实现并且用于制作少群截面参数。此蒙卡泄漏修正模型通过在连续能量的蒙卡模拟以及求解B1均匀化方程之间迭代,最终得到蒙卡模拟下的临界状态。通过此方法得到的少群截面参数较其他蒙卡以及确定论方法有两点显著优势:用于求解B1均匀化方程的少群常数是用通过临界状态的通量谱得到的;考虑了泄漏效应的蒙卡模拟可以更真实地反映组件计算时的能谱状态。为验证此泄漏修正模型,一个由连续能量的TRIPOLI4模拟而得到的数值临界实验被用于分析与比较。通过与其他蒙卡程序SERPENT以及确定论程序ECCO进行结果对比,可证明此B1泄漏修正方法能够给出更精确的用于堆芯计算的少群截面参数。     

Fast computation of MadGraph amplitudes on graphics processing unit (GPU)

Continuing our previous studies on QED and QCD processes, we use the graphics processing unit (GPU) for fast calculations of helicity amplitudes for general Standard Model (SM) processes. Additional HEGET codes to handle all SM interactions are introduced, as well as the program MG2CUDA that converts arbitrary MadGraph generated HELAS amplitudes (FORTRAN) into HEGET codes in CUDA. We test all the codes by comparing amplitudes and cross sections for multi-jet processes at the LHC associated with production of single and double weak bosons, a top-quark pair, Higgs boson plus a weak boson or a top-quark pair, and multiple Higgs bosons via weak-boson fusion, where all the heavy particles are allowed to decay into light quarks and leptons with full spin correlations. All the helicity amplitudes computed by HEGET are found to agree with those computed by HELAS within the expected numerical accuracy, and the cross sections obtained by gBASES, a GPU version of the Monte Carlo integration program, agree with those obtained by BASES (FORTRAN), as well as those obtained by MadGraph. The performance of GPU was over a factor of 10 faster than CPU for all processes except those with the highest number of jets.     

基于蒙特卡罗正算输运的全局减方差方法

蒙特卡罗方法是当前形势下辐射屏蔽计算的首选分析工具。小概率深穿透问题则是屏蔽计算的关键与亟待解决的核心问题,需要使用有效的减方差技巧。针对全局问题,利用蒙特卡罗正算输运得到的粒子通量或探测响应来构建权重窗参数,将现有的粒子位置偏移拓展到位置和能量偏倚。利用国际屏蔽基准题进行测试验证,通过使用该方法,粒子被引导到模型的所有位置。平均相对误差降低到10%以下,几乎所有网格区域都有粒子统计。结果表明,基于蒙特卡罗正算输运的输运偏倚参数构建方法能够实现全局减方差。     

Incomplete fusion reactions: Analysis of excitation functions and recoil range distributions in 16O + 51V

Incomplete fusion reactions were investigated by measuring the excitation functions of nine evaporation residues in ¹⁶O + ⁵¹V reaction in the beam energy 4-6 MeV/amu, using the well-known recoil catcher technique and gamma-ray spectrometry. The experimental data were compared with that obtained from Monte Carlo simulation calculations using the PACE2 code. The results indicate the presence of incomplete fusion process in the production of two alpha emission products. This was further confirmed by the measurement of recoil range distribution of these isotopes at 96 MeV beam energy. Calculations of the average angular momentum associated with these products revealed the peripheral nature of these ICF reactions. Received: 20 June 2001 / Accepted: 11 September 2001     

热辐射输运问题的高效蒙特卡罗模拟方法

惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升.     

Competition of fusion and quasifission in reactions leading to production of superheavy elements

The mechanism of fusion hindrance, an effect observed in the reactions of cold, warm, and hot fusion leading to production of superheavy elements, is investigated. A systematics of transfermium production cross sections is used to determine fusion probabilities. The mechanism of fusion hindrance is described as a competition of fusion and quasifission. Available evaporation residue cross sections in the superheavy region are reproduced satisfactorily. Analysis of the measured capture cross sections is performed and a sudden disappearance of the capture cross sections is observed at low fusion probabilities. A dependence of the fusion hindrance on the asymmetry of the projectile-target system is investigated using the available data. The most promising pathways for further experiments are suggested.     

多群输运-燃耗耦合蒙特卡罗计算

介绍了三维多群中子输运-燃耗耦合P3近似蒙特卡罗程序MCMG-BURN，该程序是在连续截面蒙特卡罗程序MCNP和反应堆栅元均匀化程序WIMS基础上发展而来的，使用多群截面为模拟临界实验堆和高通量工程实验堆，取得与MCNP和实验一致的结果，MCMG-BURN具有与MCNP相同的精度，但计算时间较MCNP要少得多。     

热辐射输运问题的隐式蒙特卡罗方法求解

热辐射与物质相互作用及辐射光子在物质中的传输是惯性约束聚变研究中的重要课题. 介绍了基于隐式蒙特卡罗方法的辐射输运方程,在该方程的积分-微分形式基础上,推导了利于蒙特卡罗方法模拟的等价的积分输运方程;基于积分方程设计数值模拟流程,编写三维蒙特卡罗数值模拟程序;针对热辐射输运典型问题及benchmark问题开展了数值实验,计算结果验证了方法的适应性及程序的正确性. 关键词： 热辐射 惯性约束聚变 输运方程 隐式蒙特卡罗