生成超重核的熔合反应中的动力学形变效应

利用包含动力学势能面的双核模型对超重核的生成机制中的一些问题进行了研究。对双核系统的粒子交换势能面进行的计算结果表明,反应过程中原子核动力学形变对于粒子交换势能面的结构有显著的影响。进一步计算了生成超重核的熔合几率,结果显示,原子核的动力学形变导致内熔合位垒升高,进而明显降低了生成超重核的熔合几率。Some aspects in the fusion mechanism for the production of superheavy nuclei are investigated with the dinuclear system model with dynamical potential energy surface. The calculation results about the potential energy surface indicate that the inclusion of nuclear dynamical deformation affects the structure of potential energy surface significantly.The investigation on the fusion probability to synthesize superheavy nuclei indicates that the fusion probability decreases significantly due to the increase of the inner fusion barrier for the inclusion of the nuclear dynamical deformation.     

基于双核模型对超重元素合成机制的研究

在双核模型基础上，考虑了熔合与准裂变的竞争，通过数值法求解主方程，计算了⁵⁰Ti，⁵⁸Fe+²⁰⁸Pb，²⁰⁹Bi这4个反应系统通过冷熔合反应合成超重元素的激发函数，得到了与实验比较符合的结果.计算了不同入射能量时各角动量分波对熔合概率和超重核存活概率的影响以及对蒸发剩余截面的贡献.这些结果对进一步理解超重核的合成机制有重要意义. 关键词： 超重元素 双核模型 熔合反应 蒸发剩余截面     

弱束缚原子核引起的熔合反应机制研究

弱束缚原子核引起的熔合反应机制研究是近几十年中实验核物理研究的重要课题之一。相比于放射性核束,加速器提供的稳定弱束缚原子核束流的强度要高几个数量级,利用稳定弱束缚原子核作为弹核进行反应机制的研究,可以在保证统计性和准确性的基础上,深入研究原子核的破裂、转移等反应道对熔合过程的耦合作用。已有很多实验数据表明,在库仑位垒附近,弱束缚原子核引发的熔合反应有很多有趣的现象,例如完全熔合截面的“垒下增强”和“垒上压低”。本文主要回顾近年来弱束缚原子核“垒上压低”的研究结果,并探讨造成“垒上压低”的可能原因。完全熔合截面“垒上压低”的主要原因是弱束缚原子核在进入熔合位垒之前发生破裂,从而降低了完全熔合反应道的入射通量。同时,实验研究表明完全熔合截面压低的程度可能与靶核质量数以及靶核结构相关。目前,在实验上对弱束缚原子核引起的熔合反应研究主要有3种测量方法,分别为$\gamma$射线测量方法、带电粒子测量方法以及带电粒子-$\gamma$射线符合测量的方法。其中,带电粒子-$\gamma$射线符合测量的方法在反应道鉴别方面具有明显的优势。本文对这3种测量方法进行了概要介绍,并就国内外对运用这3种方法开展的研究进行了介绍,包括本研究组在此方面的研究工作。此外,对弱束缚原子核引起的熔合反应近期在理论方面的研究工作也做了些介绍。     

热熔合反应合成超重核产生截面的同位素依赖性

通过在形成超重核的重离子俘获和熔合过程中引入位垒分布函数的方法对双核模型做了进一步发展. 超重核形成过程中的俘获、熔合和蒸发3个阶段分别采用了半经验的耦合道模型、数值求解主方程和统计蒸发模型的方法来描述. 计算了近年来Dubna小组利用热熔合反应⁴⁸Ca(²⁴³Am, 3n—5n)^288—286115和⁴⁸Ca(²⁴⁸Cm, 3n—5n)^293—291116合成超重新核素的蒸发余核激发函数. 系统分析了⁴⁸Ca轰击锕系元素U,Np,Pu,Am,Cm合成超重核Z=112—116产生截面的同位素依赖性. 给出了合成超重新核素最佳的弹靶组合和入射能量, 即有最大的超重核产生截面. 计算说明, 壳修正能和中子分离能是影响超重核生成截面产生同位素依赖性的主要因素.     

重离子碰撞中熔合位垒的研究

熔合位垒的研究对熔合反应以及超重核合成有重要的意义。在改进的同位旋相关的量子分子动力学(ImIQMD)模型框架下,提取了熔合反应体系40Ca+40Ca,48Ca+208Pb,48Ca+204Pb和16O+154Sm的熔合位垒。研究了壳修正能对熔合位垒的影响、动力学位垒的能量依赖性、同位旋效应以及形变核的方向效应。计算发现壳修正能降低了熔合反应的位垒。在研究动力学位垒的能量依赖性时,发现位垒高度和位垒半径表现出相反的能量依赖行为。在动力学反应中,当两个核距离接近时,缺中子体系的库仑势同样表现出一定的能量依赖性。对于丰中子体系,由于中子屏蔽作用,库仑势基本没有能量依赖性。     

重离子位垒下熔合——多维位垒透射及动力学过程的研究

前言入射能量接近和低于库仑位垒能区的重离子熔合反应是当前重离子反应研究中一个十分活跃的前沿领域。实验和理论的探索已经证实并描述了高于位垒的低能区重离子(轻核到中重核)反应主要经历的全熔合反应机制。两个硬球在一维(径向距离)库仑+核势相互作用下碰撞、并越过熔合位垒熔合的图像给出了对于垒上低能区全熔合截面的相当好的拟合。然而,相当简化的一维位垒模型的成功是由于垒上低能区全熔合截面一般说来并不灵敏于其它自由度。     

双核系统核子交换势能面

实验证实,用重离子碰撞合成超重原子核时,准裂变对熔合的抑制是非常重要的.碰撞中双核系统间核子转移所形成的势能面称之为驱动势,它制约核子转移,因而决定熔合与准裂变的竞争.双核系统势能面还提供重离子碰撞合成超重原子核的最佳激发能和最佳弹靶组合的信息.     

超重原子核与新元素研究

当前,原子核物理研究的一个重要前沿是探索原子核的电荷与质量极限,研究超重原子核与超重元素的性质,以及合成超重原子核。20世纪60年代,基于量子壳效应,理论预言质子数为114、中子数为184的原子核及其相邻核具有较长的寿命,甚至可能是稳定的,形成一个超重稳定岛。这个理论预言促进了重离子加速器及相关探测设备的建造,推动了重离子物理的发展。到目前,已经合成到了118号元素,填满了元素周期表的第7行。然而,合成更重的超重元素或包含更多中子的超重原子核面临着很多挑战,需要理论与实验密切结合,探索超重原子核的性质与合成机制,以登上超重稳定岛。文章概要评述超重原子核与新元素研究。首先介绍超重原子核与超重元素研究的背景及理论预言,包括超重核存在的根源、理论预言的概况等。之后简要给出实验合成超重核取得的主要进展和新元素命名情况。关于合成更重的超重元素面临的挑战,文章将针对利用重离子熔合蒸发反应合成超重核的截面低、所合成的超重核缺中子等情况展开讨论。最后评述近年来超重原子核结构性质、衰变、裂变与合成机制等方面的理论研究进展,包括超重核区的幻数和超重岛的位置,超重核的稳定性,利用重离子熔合蒸发反应合成超重核的三步过程及其复杂性,利用多核子转移合成超重核的探索,等等。The exploration of charge and mass limits of atomic nuclei and the synthesis of long-lived or stable superheavy nuclei (SHN) are at the frontier of modern nuclear physics. In the 1960s, based on the stability originating from quantum shell effects, the possible existence of an island of stability around 298114 was predicted. This prediction advanced the construction of heavy ion accelerators and detectors and the development of heavy ion physics. So far, superheavy elements (SHE) with Z up to 118 have been synthesized via heavy ion fusion reactions in laboratories. Recently the IUPAC/IUPAP Joint Working Party (JWP) concluded that criteria for the discovery of new elements have been met for those with Z=113, 115, 117 and 118. Therefore the seventh period of the periodic table of elements is completed. To synthesize even heavier elements or more neutron-rich SHN by using heavy ion fusion reactions, one confronts many challenges. More efforts should be made to study the properties of SHN both experimentally and theoretically. In this short review on the study on SHN and SHE, we will first introduce the background and theoretical predictions of SHN, including the origin of the possible existence of SHN and the predicted island of stability of SHN, etc. Then we will present progresses made up to now concerning the synthesis of SHN and the naming of the four new elements. As for the challenges nuclear physicists confront in synthesizing even heavier SHEs, we will detail those connected with heavy ion fusion-evaporation reactions, namely, the tiny cross sections to produce SHN and the fact that only neutron-deficient SHNs can be synthesized. Finally we will discuss some theoretical progresses on the study of SHN, including the structure of SHN and proton and neutron magic numbers after 208Pb, the stability and the synthesis mechanism of SHN as well as what we should focus on in the future.     

用QMD模型研究晕核 11 Be+208Pb的熔合机制

通过用QMD模型研究晕核^11Be ^208Pb的近垒熔合反应,发现晕核引起的熔合反应中,并存着两种相互竞争的机制：一方面当入射晕核^11Be靠近靶核时,由于^11Be是弱束缚体系,与靶核的相互作用可使其很容易破裂或少数核子被靶核俘获形成核子转移反应,从而对于熔合表现出压制;另一方面当^11Be的少数中子进入靶核并与靶核相互作用而使得靶核有些激发,而使局部半径增大,导致熔合势垒降低,熔合截面增强。用QMD模型计算出的熔合截面与实验值基本符合,垒附近表现出增强效应。     

五参数哑铃模型中核熔合反应的新奇不稳定度

在原子核熔合反应或裂变反应的动力学过程中,颈部变量都起着非常重要的作用。通常情况下,模型中变量的定义根据各种描述宏观形变的参数的不同而不同。为了进一步研究颈部增长在核熔合反应中的作用,通过引入弹核和靶核的形变参量来拓展哑铃模型。在此框架下,计算了从两个接触核到生成一个复合核过程的核熔合反应系统在不同形变下的势能曲面,同时又对核熔合反应体系的惯性张量和粘滞张量进行了计算,这些物理量的计算都为以后的朗之万动力学研究打下了基础。计算结果表明,在改进后的五参数哑铃模型之中,颈部变量是不稳定度并在核反应中起重要作用。     

Systematic study of survival probability of excited superheavy nuclei

The stability of excited superheavy nuclei (SHN) with 100 Z 134 against neutron emission and fission is investigated by using a statistical model. In particular, a systematic study of the survival probability against fission in the 1n-channel of these SHN is made. The present calculations consistently take the neutron separation energies and shell correction energies from the calculated results of the finite range droplet model which predicts an island of stability of SHN around Z = 115 and N = 179. It turn...     

Systematic study of survival probability of excited superheavy nuclei

The stability of excited superheavy nuclei (SHN) with 100 ⩽ Z ⩽ 134 against neutron emission and fission is investigated by using a statistical model. In particular, a systematic study of the survival probability against fission in the 1n-channel of these SHN is made. The present calculations consistently take the neutron separation energies and shell correction energies from the calculated results of the finite range droplet model which predicts an island of stability of SHN around Z = 115 and N = 179. It turns out that this island of stability persists for excited SHN in the sense that the calculated survival probabilities in the 1n-channel of excited SHN at the optimal excitation energy are maximized around Z = 115 and N = 179. This indicates that the survival probability in the 1n-channel is mainly determined by the nuclear shell effects.

     

理论预测超重核274-291Cn和266-287Ds的衰变模式

自发裂变和α衰变是影响超重核稳定性的两个主要因素。为了探索270Ds附近的长寿命的超重核,系统地计算了电荷数在104 ≤ Z ≤ 112范围内的α衰变与自发裂变之间的竞争。采用推广的液滴模型和唯象的解析公式计算了α衰变半衰期。基于包括壳效应和同位旋效应的WKB近似方法估算了相同超重核的自发裂变半衰期,进而预测了未知超重核274-276,279Cn与267-269Ds的衰变模式。The stability of superheavy nuclei (SHN) is controlled mainly by spontaneous fission and α decay processes. To investigate whether long lived SHN could really exist around 270Ds, the competition between α decay and spontaneous fission in the region 104 ≤ Z ≤ 112 are studied systematically. The α decay half-lives are investigated by employing a generalized liquid drop model (GLDM) and phenomenological analytical formula. Calculations of spontaneous fission half-lives for the same SHN are carried out based on the Wenzel-Kramers-Brillouin(WKB) approximation with both the shell effect and the isospin effect included. Decay modes are predicted for the unknown nuclei 274-276,279Cn and 267-269Ds.     

Fission fragment mass yields of Th to Rf even-even nuclei

Fission properties of the actinide nuclei are deduced from theoretical analysis. We investigate potential energy surfaces and fission barriers and predict the fission fragment mass yields of actinide isotopes. The results are compared with experimental data where available. The calculations were performed in the macroscopic-microscopic approximation with the Lublin-Strasbourg Drop (LSD) for the macroscopic part, and the microscopic energy corrections were evaluated in the Yukawa-folded potential. The Fourier nuclear shape parametrization is used to describe the nuclear shape, including the non-axial degree of freedom. The fission fragment mass yields of the nuclei considered are evaluated within a 3D collective model using the Born-Oppenheimer approximation.     

On the multiple-humped fission barriers and half-lives of actinides

The energy of actinide nuclei has been determined within a generalized liquid drop model taking into account the proximity energy, the mass and charge asymmetry, an accurate nuclear radius in adding the shell and pairing energies. Double and triple-humped potential barriers appear. The second maximum corresponds to the transition from compact and creviced one-body shapes to two touching ellipsoids. A third minimum and third peak appear in special asymmetric exit channels where one fragment is almost a magic nucleus with a quasi-spherical shape while the other one evolves from oblate to prolate shapes. The heights of the double and triple-humped fission barriers agree precisely with the experimental results in all the actinide region. The predicted half-lives follow the experimental data trend.     

Spontaneous Fission Barriers Based on a Generalized Liquid Drop Model

The barrier against the spontaneous fission has been determined within the Generalized Liquid Drop Model （GLDM） including the mass and charge asymmetry, and the proximity energy. The shell correction of the spherical parent nucleus is calculated by using the Strutinsky method, and the empirical shape-dependent shell correction is 6mp10yed during the deformation process. A quasi-molecular shape sequence has been defined to describe the whole process from one-body shape to two-body shape system, and a two-touching-ellipsoid is adopted when the superdeformed one-body system reaches the rupture point. On these bases the spontaneous fission barriers are systematically studied for nuclei from 2a~Th to 249 Cm for different possible exiting channels with the different mass and charge asymmetries. The double, and triple bumps are found in the fission potential energy in this region, which roughly agree with the experimental results. It is found that at around Sn-like fragment the outer fission barriers are lower, while the partner of the Sn-like fragment is in the range near l~SRu where the ground-state mass is lowered by allowing axially symmetric shapes. The preferable fission channels are distinctly pronounced, which should be corresponding to the fragment mass distributions.     

Asymmetry of mass and charge division in spontaneous fission

The order-disorder model has been used to explain asymmetry of mass and charge division and related phenomena in fission. According to this model the fission process involves two steps consisting of charge polarisation into two ‘impending fragment clusters’ with beta stable neutron numbers and subsequent distribution of the balance neutrons between the two. Mode of elemental division of the fissioning nuclei is attributed to the charge polarisation in the first step. Theory of reaction rate has been applied to the system. The frequency term is obtained by applying the conditional stochastic process under charge polarisation constraint and the energy-dependent term is given by the condition of minimum in free energy of the system. Using this, the relative probability of polarisation into given charge pair is arrived at. The model uses stable neutron numbers for the charges as the only input. No explicit assumption or quantification on the preference of formation of shell closure species in fission is necessary. The statistics developed on the principle of equala priori probability of all charge polarisation is used. The shell effects come into play only in deciding a stable neutron number for the charges. The total isotopic yield distribution for a number of fission reactions shows asymmetry in the actinide region which reduces with increasing mass/charge of the fissioning nuclide and bunching of the higherz peaks. Although the mass yields obtained therefrom for a number of fission reactions agree with experimental results, the peaks of the distributions are slightly shifted away from the symmetric point and the distributions are somewhat narrower. Charge distribution parameters obtained from these results are also presented.     

Potential energy surfaces and fission fragment mass yields of even-even superheavy nuclei

Potential energy surfaces and fission barriers of superheavy nuclei are analyzed in a macroscopic-microscopic model. The Lublin-Strasbourg Drop (LSD) model is used to obtain the macroscopic part of the energy, whereas the shell and pairing energy corrections are evaluated using the Yukawa-folded potential; a standard flooding technique is utilized to determine barrier heights. A Fourier shape parametrization containing only three deformation parameters is shown to effectively reproduce the nuclear shapes of nuclei approaching fission. In addition, a non-axial degree of freedom is taken into account to better describe the structure of nuclei around the ground state and in the saddle region. In addition to the symmetric fission valley, a new highly asymmetric fission mode is predicted in most superheavy nuclei. The fission fragment mass distributions of the considered nuclei are obtained by solving 3D Langevin equations.