-
氘核削裂反应是单核子转移反应的一种。图1是转移反应A(a, b)B的示意图,在这个过程中,粒子x从弹核a(a=b+x)转移到靶核A上,形成剩余核B(B=A+x)。通常我们把b和A称为a和B的核心,而把x称为价粒子(或价核子)。对氘核削裂反应A(d, p)B而言,x为被转移的中子,b为氘核中的质子。这个体系的哈密顿量可以根据不同的反应道表示为
$$ H = H_\alpha+T_\alpha+V_\alpha = H_\beta+T_\beta+V_\beta, $$ (1) 其中下标
$ \alpha $ 和$ \beta $ 分别表示入射道和出射道,$ H_{\alpha,\beta} $ 为内部坐标的哈密顿量,对应于原子核a、A、b、B的内部坐标波函数,$ T_{\alpha,\beta} $ 为动能算符,$ V_{\alpha,\beta} $ 为反应道相互作用($ V_\alpha = V_{\rm bA}+V_{\rm xA} $ ,$ V_\beta = V_{\rm bA}+V_{\rm bx} $ )。对于转移反应($ \alpha\neq\beta $ ),其跃迁振幅(transition amplitude)为[43]$$ T_{\beta\alpha} = \langle \varPsi_\beta^{(-)}|V_\alpha|\phi_\alpha\rangle = \langle \phi_\beta|V_\beta|\varPsi_\alpha^{(+)}\rangle, $$ (2) 这里前后两个表达式分别为跃迁振幅的前项(prior)和后项(post)形式,
$ \varPsi_{\alpha}^{(+)} $ 和$ \varPsi_\beta^{(-)} $ 分别是三体反应体系在初态({(b,x),A})和末态({b,(x,A)})的总波函数。它们的上下标表示它们的边界条件分别为出射球面波($ (+) $ )和入射球面波($ (-) $ )。$ \phi_\alpha $ 是如下方程的解:$$ (E-H_\alpha-T_\alpha)\phi_\alpha = 0, $$ (3) 这里的
$ E $ 为反应体系的质心系能量。$ \phi_\beta $ 对应的薛定谔方程可以同样写出。显然,$ \phi_\alpha $ 是平面波和入射道原子核的内部波函数$ \varPhi_\alpha({{\boldsymbol{A}}},{{\boldsymbol{a}}}) $ 的乘积:$$ \phi_\alpha = \exp({\rm{i}}{\boldsymbol{k}}_\alpha \boldsymbol\cdot {\boldsymbol{r}}_\alpha)\varPhi_\alpha({{\boldsymbol{A}}},{{\boldsymbol{a}}}), $$ (4) 其中
$ {\boldsymbol{k}}_\alpha = \sqrt{2\mu_\alpha E_{\rm{cm}}} $ 为入射粒子的动量,$ \mu_\alpha $ 为A-a系统的约化质量,$ E_{\rm{cm}} $ 为A-a的相对运动动能。转移反应A(a, b)B的a(= b+x)-A这个三体系统的波函数$ \varPsi_\alpha $ 满足薛定谔方程$$ (E-H_\alpha-T_\alpha)\varPsi_\alpha^{(+)} = V_\alpha\varPsi_\alpha^{(+)}{\text{。}} $$ (5) 形式上,这个方程的解为
$$ \begin{split} \varPsi_\alpha^{(+)} =& \phi_\alpha+ \frac{1}{E-H_\alpha-T_\alpha-V_\alpha+{\rm{i}}\varepsilon}V_\alpha\phi_\alpha \\ =& \left(1+\frac{1}{E-H+{\rm{i}}\varepsilon}\right)\phi_\alpha{\text{。}} \end{split} $$ (6) 把式(6)代入式(2),有(后项形式):
$$ T_{\beta\alpha} = \big\langle \phi_\beta|V_\beta+V_\beta\frac{1}{E-H+{\rm{i}}\varepsilon}V_\alpha|\phi_\alpha\big\rangle{\text{。}} $$ (7) 我们可以在方程(5)两边同时添加一个任意的辅助势
$ U_\alpha $ ,$$ (E-H_\alpha-T_\alpha-U_\alpha)\varPsi_\alpha^{(+)} = (V_\alpha-U_\alpha)\varPsi_\alpha^{(+)}, $$ (8) 则方程的解可以表示为
$$ \begin{split} \varPsi_\alpha^{(+)} =& \chi_\alpha^{(+)}+\frac{1}{E-H_\alpha-T_\alpha-U_\alpha+{\rm{i}}\varepsilon}(V_\alpha-U_\alpha)\varPsi_\alpha^{(+)}\\ =& \chi_\alpha^{(+)}+\frac{1}{E-H+{\rm{i}}\varepsilon}(V_\alpha-U_\alpha)\chi_\alpha^{(+)}{\text{。}} \end{split} $$ (9) 其中
$ \chi_\alpha^{(+)} $ 是如下方程的解:$$ (E-H_\alpha-T_\alpha-U_\alpha)\chi_\alpha^{(+)} = 0{\text{。}} $$ (10) 把式(9)代入式(2),我们有[43]:
$$ \begin{split} T_{\beta\alpha} =& \big\langle\phi_\beta|V_\beta-V_\alpha+U_\alpha|\chi_\alpha^{(+)}\big\rangle +\\& \big\langle \varPsi_\beta^{(-)}|V_\alpha-U_\alpha|\chi_\alpha^{(+)}\big\rangle{\text{。}} \end{split}$$ (11) 如果对出射道也引入一个任意的势
$ U_\beta $ ,跃迁振幅还可以表达如下:$$ \begin{split} T_{\beta\alpha} =& \big\langle\chi_\beta^{(-)}|V_\alpha-V_\beta+U_\beta|\phi_\alpha\big\rangle + \\&\big\langle \chi_\beta^{(-)}|V_\beta-U_\beta|\varPsi_\alpha^{(+)}\big\rangle{\text{。}}\end{split} $$ (12) -
原则上式(8)的辅助势
$ U_\alpha $ 可以作用于原子核的内部坐标,从而引起原子核内部状态波函数的改变。但是如果我们进一步假定$ U_\alpha $ 是$ \alpha $ 道(入射道)的光学势,即它只是a和A两个核质心之间距离的函数,那么$ \chi_\alpha^{(+)} $ 可以进一步分解为a和A的相对运动波函数$ \chi_\alpha^{(+)}({\boldsymbol{r}}_\alpha) $ 以及它们的内部运动波函数$ \varPhi_\alpha({{\boldsymbol{A}}},{{\boldsymbol{a}}}) $ 的乘积:$$ \chi_\alpha^{(+)}({{\boldsymbol{r}}_\alpha},{\boldsymbol{A}},{\boldsymbol{a}}) = \chi_{\alpha}^{(+)}({{\boldsymbol{r}}}_{\alpha})\varPhi_\alpha({\boldsymbol{A}},{\boldsymbol{a}}){\text{。}} $$ (13) 同样地,如果把
$ U_\beta $ 取作出射道的光学势,有:$$ \chi_\beta^{(-)} ({\boldsymbol{r}}_\beta, {\boldsymbol{B}},{\boldsymbol{b}}) = \chi_\beta^{(-)} ({\boldsymbol{r}}_\beta)\varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}}), $$ (14) 这里
$ \chi_\beta^{(-)} $ 是出射道核b和B的相对运动波函数,$ \varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}}) $ 为b和B的内部状态波函数。此外,在取
$ U_\alpha $ 和$ U_\beta $ 为入射道和出射道的光学势的情况下,式(11)和(12)右边的第一项是对角化的。把式(9, 13, 14)代入式(12),有:$$ \begin{split} T_{\beta\alpha} =& \big\langle\chi_\beta^{(-)}\varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}})|U_\alpha|\phi_\alpha\big\rangle \delta_{\alpha\beta}+ \\ &\big \langle \chi_\beta^{(-)}\varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}})|(V_\beta-U_\beta)+\\ &(V_\beta-U_\beta)\frac{1}{E-H+{\rm{i}}\varepsilon}(V_\alpha-U_\alpha)| \chi_\alpha^{(+)}{\varPhi_\alpha}({{\boldsymbol{A}}},{{\boldsymbol{a}}})\big\rangle{\text{。}} \end{split}$$ (15) 这就是两体反应跃迁振幅的一般形式。对于核子转移反应,由于
$ \beta\neq\alpha $ ,其跃迁振幅就只包含等号右边的第二项,即:$$ \begin{split} T_{\beta\alpha} = &\big\langle \chi_\beta^{(-)}\varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}})|(V_\beta-U_\beta)+\\ &(V_\beta-U_\beta)\frac{1}{E-H+{\rm{i}}\varepsilon}(V_\alpha-U_\alpha)| \chi_\alpha^{(+)}\varPhi_\alpha({{\boldsymbol{A}}},{{\boldsymbol{a}}})\big\rangle{\text{。}} \end{split}$$ (16) 注意,这里除了假定引入的辅助势
$ U_\alpha $ 和$ U_\beta $ 为入射道和出射道的光学势以外(即这些势只是a-A以及b-B之间的两体相互作用,只能引起a-A以及b-B体系的弹散或者非弹,而不引起核子交换等过程),没有引入其它近似。式(7)和式(16)中的跃迁振幅看起来简单,但是其严格求解却是复杂的,因为这些式子里面的
$ H $ 是三体体系的完整哈密顿量。人们通常需要引入各种近似来计算转移反应截面。 -
忽略式(7)中的第二项即可得到转移反应的PWBA,即
$$ \begin{split} T_{\beta\alpha}^{\rm{PWBA}} = &\big\langle \phi_\beta|V_\beta|\phi_\alpha\big\rangle = \big\langle \exp({\rm{i}}{\boldsymbol{k}}_\beta \cdot{\boldsymbol{r}}_\beta)\times \\&\varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}})|V_\beta| \exp({\rm{i}}{\boldsymbol{k}}_\alpha \cdot{\boldsymbol{r}}_\alpha)\varPhi_\alpha({{\boldsymbol{A}}},{{\boldsymbol{a}}})\big\rangle{\text{。}} \end{split} $$ (17) 式中
$ V_\beta = V_{\rm bA}+V_{\rm bx} $ ($ V_{\rm bA} $ 和$ V_{\rm bx} $ 分别是图1中b和A以及b和x之间的相互作用)。其中,$ V_{\rm bx} $ 对转移反应微分截面角分布形状的影响是主要的。$ V_{\rm bA} $ 对转移反应截面的贡献是平滑的本底,因此在PWBA计算中通常被忽略掉。在进一步分析PWBA之前,我们需要对转移反应跃迁振幅中原子核的内部波函数做一说明。在式(4, 13, 14)等式中,
$ \varPhi_\alpha $ 和$ \varPhi_\beta $ 分别是出射道和入射道原子核的内部坐标波函数,它们分别对应内部哈密顿量$ H_\alpha $ 和$ H_\beta $ 。在氘核削裂反应A(d, p)B中,$ \varPhi_\alpha $ 是A和氘核的内部坐标波函数($ \varPhi_{\rm A} $ 和$ \phi_{\rm d} $ )的乘积:$$ \varPhi_\alpha({{\boldsymbol{A}}},{{\boldsymbol{a}}}) = \varPhi_{J_A}^{M_A}({\boldsymbol{A}})\phi_d({\boldsymbol{r}}), $$ (18) 而
$ \varPhi_\beta $ 则是B的内部坐标波函数:$$ \varPhi_\beta({\boldsymbol{B}},{\boldsymbol{b}}) = \varPhi_{J_B}^{M_B}({\boldsymbol{A}},{\boldsymbol{r}}_n){\text{。}} $$ (19) 这里的
$ J_{A,B} $ 和$ M_{A,B} $ 分别是A和B的自旋及其投影。$ {\boldsymbol{r}} $ 和$ {\boldsymbol{r}}_n $ 分别对应图1中的$ {\boldsymbol{r}}_2 $ 和$ {\boldsymbol{r}}_1 $ 。$ \varPhi_{J_B}^{M_B}({\boldsymbol{A}},{\boldsymbol{r}}_n) $ 可以做如下展开:$$\begin{split} \varPhi_{J_B}^{M_B}({\boldsymbol{A}},{\boldsymbol{r}}_n) =& \sum\limits_{A'}\sum\limits_{\substack{nlj\\ M_j M_A}}S_{A';nlj}^{1/2}C_{jM_jJ_AM_A}^{J_BM_B} \times\\& \varPhi_{A',J_A M_A}({\boldsymbol{A}})\phi_{nlj M_j}({\boldsymbol{r}}_n){\text{。}} \end{split}$$ (20) 这里
$ A' $ 对应原子核A的不同的态,$ S_{A';nlj}^{1/2} $ 是谱幅度,它的平方即为谱因子$ S_{A';nlj} $ (注意谱因子总是正的,但谱幅度则可正可负),而$ \phi_{nlj M_j} $ 则是被转移粒子的单粒子波函数。谱因子$ S_{A';nlj} $ 体现了在多大程度上可以把B核(自旋及其投影为$ J_BM_B $ )看作是由A核(处于$ A' $ 态,自旋及其投影为$ J_AM_A $ )以及一个处于量子数为$ nlj $ 的单粒子轨道的中子n所组成的。为了简单起见,我们忽略了这些式子中粒子n的自旋波函数。在数值计算中需要对式(20)中的每一对
$\varPhi_{A',J_A M_A}({\boldsymbol{A}}) \phi_{nlj M_j}({\boldsymbol{r}}_n)$ 进行计算。下面我们取式 (20)等号右边的某一个成分,代入式(17),对氘核削裂反应A(d, p)B,有:$$ \begin{split} T_{\rm dp}^{\rm{PWBA}} =& S_{nlj}^{1/2}C_{jM_jJ_AM_A}^{J_BM_B} \big\langle \exp({\rm{i}}{\boldsymbol{k}}_{\rm p} \cdot{\boldsymbol{r}}_{\rm p})\phi_{nlj M_j}({\boldsymbol{r}}_n) \\& |V_{\rm np}({\boldsymbol{r}})| \exp({\rm{i}}{\boldsymbol{k}}_{\rm d} \cdot{\boldsymbol{r}}_{\rm d})\phi_{\rm d}({\boldsymbol{r}})\big\rangle{\text{。}} \end{split}$$ (21) 这里
${\boldsymbol{r}}_{\rm d}$ 和${\boldsymbol{r}}_{\rm p}$ 分别对应图1中的$ {\boldsymbol{r}}_\alpha $ 和$ {\boldsymbol{r}}_\beta $ 。注意,由于$V_{\rm np}$ 不作用在$ \varPhi_{J_A}^{M_A}({\boldsymbol{A}}) $ 的内部坐标上,可以利用A的内部波函数的正交性来避免对A的内部坐标的积分。对于(d, p)反应,利用图1中的坐标关系,有:
$$ \begin{split} {\boldsymbol{k}}_{\rm d}\cdot{\boldsymbol{r}}_{\rm d}-{\boldsymbol{k}}_{\rm p} \cdot{\boldsymbol{r}}_{\rm p} =& ({\boldsymbol{k}}_{\rm d}-\frac{A}{A+1}{\boldsymbol{k}}_{\rm p}){\boldsymbol{r}}_n +({\boldsymbol{k}}_{\rm p}-\frac{1}{2}{\boldsymbol{k}}_{\rm d}){\boldsymbol{r}} \\[1mm]\equiv&{\boldsymbol{q}} \cdot{\boldsymbol{r}}_n+{\boldsymbol{K}} \cdot{\boldsymbol{r}} , \\[-11pt] \end{split}$$ (22) 于是
$$ \begin{split}T_{\rm dp}^{\rm{PWBA}} \propto & \int {\rm e}^{-{\rm{i}}{\boldsymbol{k}}_{\rm p} \boldsymbol\cdot{\boldsymbol{r}}_{\rm p}}\phi_{nlj M_j}^\ast({\boldsymbol{r}}_n) V_{\rm np}({\boldsymbol{r}}) {\rm{e}}^{{\rm{i}}{\boldsymbol{k}}_{\rm d} \cdot{\boldsymbol{r}}_{\rm d}}\phi_{\rm d}({\boldsymbol{r}}){\rm{d}}{\boldsymbol{r}}_{\rm p} {\rm{d}}{\boldsymbol{r}}_{\rm d}= \\ &\int {\rm{e}}^{{\rm{i}}{\boldsymbol{K}} \cdot{\boldsymbol{r}}}V_{\rm np}({\boldsymbol{r}})\phi_{\rm d}({\boldsymbol{r}}){\rm{d}}{\boldsymbol{r}} \int {\rm{e}}^{{\rm{i}}i{\boldsymbol{q}} \boldsymbol\cdot{\boldsymbol{r}}_n}\phi_{nlj M_j}^\ast({\boldsymbol{r}}_n){\rm{d}}{\boldsymbol{r}}_n{\text{。}} \end{split}$$ (23) 可见在PWBA下,氘核削裂反应的跃迁振幅可写成两个傅立叶变换的乘积的形式。其中一个是对氘核波函数和n-p相互作用(
$V_{\rm np}$ )的乘积的变换,另一个是对被转移的中子在B中的单粒子波函数的变换。我们感兴趣的是第二个变换,因为无论从核中转移哪个中子,第一个变换都是相同的。我们定义第一个变换为$ G({\boldsymbol{K}}) $ 。第二个傅立叶变换中的平面波可以展开成如下形式:$$ {\rm{e}}^{{\rm{i}}{\boldsymbol{q}} \cdot{\boldsymbol{r}}_n} = 4\pi\sum\limits_{LM}i^L j_L(qr_n)Y_{LM}^*({\hat{\boldsymbol{q}}})Y_{LM}({\hat{\boldsymbol{r}}}_n) , $$ (24) 单粒子波函数可以写成径向波函数与球谐函数相乘的形式:
$$ \varphi_{nljM_j}({\boldsymbol{r}}_n) = \frac{1}{r_n}u_{nl}(r_n)Y_{lm}({\hat{\boldsymbol{r}}}_n) {\text{。}} $$ (25) 因此,跃迁振幅又可化为如下形式:
$$ T_{\rm dp}^{\rm{PWBA}} \propto G({\boldsymbol{K}})Y_{lm}^*({\hat{\boldsymbol{q}}})\int j_l(qr_n)u_{nl}(r_n)r_n {\rm{d}}r_n {\text{。}} $$ (26) 图2中展示了用PWBA计算的对应不同
$ l $ 值的反应微分截面(正比于$ |T|^2 $ )。可以看出,用PWBA计算的角分布对单粒子波函数的角动量$ l $ 有相当大的依赖性。PWBA假定a-A之间的相对运动由平面波描述,即忽略它们之间的相互作用在反应过程中的作用。在两个核距离很近($ r $ 很小)的情况下,这个图像显然是不恰当的。不过,直接核反应通常发生在核的表面。我们可以假设当$ r $ 小于某个半径$ R $ 时,原子核之间存在的强吸收作用会压低直接过程的贡献。因此可以忽略式(26)的积分中$ r<R $ 部分的贡献。同时,考虑到$ G({\boldsymbol{K}}) $ 比起$ j_l(qr) $ 随出射的质子角度变化较为缓慢,我们有$ T^{\rm{PWBA}}\propto\int_R^\infty j_l(qr)u_{nl}(r)r{\rm{d}}r $ 。此处的$ R $ 被称为截断半径(cutoff radius)。更进一步,若认为直接反应只发生在核的表面($ r\approx R $ ),则有$ T\propto j_l(qR) $ 。另外,随着$ l $ 的增加,$ j_l(qR) $ 的第一个也是最大的峰出现在较大的$ qR $ 处。可见,不同的角动量$ l $ 对应不同的角分布及其最大值对应的散射角度。因此,可以利用实验测得的微分截面角分布来确定单粒子波函数的轨道角动量,进而确定原子核的自旋宇称。图 2 (在线彩图)由PWBA计算的15O(d, p)16O微分截面角分布对被转移中子的角动量
$l$ 的依赖性,这里$E_{\rm d}=30$ MeV。红色实线、蓝色长虚线和黑色短虚线分别对应$l=0,1,2$ 时的结果PWBA中的截断半径
$ R $ 并没有很好的定义。它通常被当作自由参数,通过拟合转移反应截面角分布来确定。因此PWBA理论不具备良好的理论预言能力。此外,PWBA忽略了出射道和入射道的原子核之间的相互作用。这显然是不合适的。针对这些问题,在氘核削裂反应被PWBA理论证实是一个有效的研究原子核谱学的工具后,人们很快就提出了改进PWBA理论的需求[44-45]。这个过程经历了十年左右。直到20世纪60年代,扭曲波玻恩近似理论(Distorted Wave Born Approximation, DWBA)才被发展到比较完善的程度[46-49]。 -
DWBA下的转移反应跃迁振幅可以通过去掉式(16)中的第二项来得到。对于氘核削裂反应A(d, p)B,有:
$$ \begin{split} T_{\rm dp}^ {\rm{DWBA}} =& \big\langle \chi_{\rm p}^{(-)}\varPhi_{J_BM_B}({\boldsymbol{A}},{\boldsymbol{r}}_n)|V_{\rm pn}+V_{\rm pA}-U_{\rm pB}|\\ & \chi_{\rm d}^{(+)}\varPhi_{J_AM_A}({\boldsymbol{A}})\phi_{\rm d}({\boldsymbol{r}})\big\rangle{\text{。}} \end{split} $$ (27) 这里
$ \chi_{\rm d}^{(+)} $ 和$ \chi_{\rm p}^{(-)} $ 分别为入射道和出射道的扭曲波(distorted waves)。它们满足的方程为(为简单起见,假定A处于基态):$$ \begin{split} & (E-E_{\rm d}-T_{\rm dA}-U_{\rm dA})\chi_{\rm d}^{(+)} = 0, \\& (E-E_{\rm n}-T_{\rm pB}-U_{\rm pB})\chi_{\rm p}^{(-)} = 0{\text{。}} \end{split} $$ (28) 其中
$E_{\rm d}$ 和$E_{\rm n}$ 分别为中子在氘核以及在剩余核B中的束缚能。通过与式(3)中的平面波所满足的方程相比,可以看出,入射道和出射道的相对运动波函数$ \chi^{(+,-)} $ 由于这里引入的辅助势$ U_{\rm dA} $ 和$ U_{\rm pB} $ 的存在而被“扭曲”了(distorted)。原则上,$ U_{\rm dA} $ 和$ U_{\rm pB} $ 可以是任意的,但是选择它们为可以描述d-A和p-B体系的弹性散射数据的光学势是自然而方便的。式(27)中的$ V_{\rm pA}-U_{\rm pB} $ 被称作是剩余项(remnant term)。可见在PWBA中被忽略的$ V_{\rm pA} $ 项可在DWBA中被严格处理。但是,对于比较重的靶核($ {\rm B} = {\rm A+n}\approx {\rm A} $ ),通常可以认为$ V_{\rm pA}\approx U_{\rm pB} $ ,从而忽略剩余项。DWBA方法大大改善了对氘核削裂反应的描述能力。图3比较了12 MeV入射能条件下DWBA和PWBA理论对40Ca(d, p)41Ca反应的计算结果。可以看到,与PWBA相比,DWBA对反应截面的幅度及角分布形状的预言都有了很大的改善。同时,由于引入的光学势的内部吸收作用,DWBA理论计算中不需要像PWBA一样人为地引入截断半径。
图 3 PWBA和DWBA对40Ca(d, p)41Ca反应预测的比较,
$E_{\rm d}=12$ MeV。实线和虚线分别对应DWBA和PWBA的计算结果,实心圆点表示实验数据。此图来自文献[50]在前面的叙述中,我们假定反应前后原子核A、B的都处于基态。事实上,它们可以被激发到某个或某些激发态。如果这些激发态对应的反应道与基态的反应道之间的耦合效应比较强的话,则必须要考虑这些激发态之间的耦合效应。此时,转移反应可以用耦合道玻恩近似(Coupled Channel Born Approximation,CCBA)理论进行分析[43, 51]。此外,从图1可以看出,转移反应的入射道和出射道的原子核只是b, x, A 这三个集团的重新组合,它们属于同一个反应体系,因此,入射道和出射道之间也存在耦合。对于反应道之间的耦合,需要用耦合反应道(Coupled Reaction Channels,CRC)理论来进行处理[50, 52-57]。本文不对此做具体的讨论。
DWBA的数值计算比较简单,也能够很好地描述较大的入射能及靶核质量范围内的直接核反应,因此,直到今天,这个理论仍然被广泛地应用于对转移反应实验数据的理论分析中。但是,以氘核削裂反应理论为代表的直接核反应理论一直在不断发展。1968年,Philpott等[49]详细分析了(p, d)、(d, p)等核子转移反应的DWBA理论中的问题,并列出了如下需要进一步研究的内容:
$ (1) $ 光学势参数的不确定性,以及这些不确定性对核子转移反应截面的影响;$ (2) $ 对一些核来讲,入射道或出射道的光学势($ U_\alpha $ 和$ U_\beta $ )只能描述它们的弹性散射,这是不够的,需要考虑非弹散射道的耦合效应;$ (3) $ 光学势的非定域性及其对核子转移反应截面的影响;$ (4) $ 核心(B=A+n中的A核)激发的影响;$ (5) $ 弹核(d、3He等)中$ D $ 态波函数的贡献;$ (6) $ 用于抽取光学势的弹性散射数据通常只对光学势的表面敏感,在半径比较小的区域(比如原子核内部)的光学势应该是什么样子?$ (7) $ 弹核破裂(projectile breakup)的影响。在此后的半个多世纪的时间里,人们对氘核削裂反应的理论研究确实也集中在这些方面,并且都取得了不少的进展,简述如下:
$ (1) $ 众所周知,拟合少量的弹性散射数据得到的光学势参数具有比较大的不确定性。因此,在理论分析中应尽可能地采用系统学光学势[58-64]。系统学光学势通过同时拟合大量的弹性散射实验数据(覆盖比较大的入射能量的范围以及比较大的靶核质量的范围),可以在很大程度上减小光学势参数的不确定性,从而在直接核反应的系统学理论分析中得到比较一致的结果[65]。近年来,Nunes等[66]还分析了由光学势的不确定性引起的转移反应截面的不确定性;$ (2) $ 人们发展了耦合道玻恩近似理论,以有效地处理耦合道效应比较强的情况下的(如靶核具有较大的形变,从而容易引起非弹激发等)转移反应理论计算,并取得了很大的成功[51];$ (3) $ 从2013年以来,Tomofeyuk、Nunes等对光学势的非定域性及其在氘核削裂反应中的应用开展了深入的研究,详情可参考文献[67-71];$ (4) $ 2015年,Moro等[72]发展了包含核心激发的氘核削裂反应理论。他们发现一般情况下核心激发对转移反应的影响比较小。但是,在核心有比较大形变的情况下,核心激发的影响还是应该考虑的。这对与从氘核削裂反应来抽取单中子谱因子等信息具有一定的影响。$ (5) $ 由于数值计算能力的提高,在转移反应理论计算中考虑d、3He等弹核的$ D $ 态波函数的贡献是早已成为一个简单的事情。当前公开的计算机程序,如FRESCO[73]、TWOFNR[74]等,都允许考虑这些成分。$ (6) $ 通常的弹性散射实验数据确实对内部光学势不敏感,因此,通过拟合弹性散射数据得到的唯象光学势的内部具有相当大的不确定性。韵小艳等[75]通过系统分析从较低(约5 MeV/nucleon)到较高(约30 MeV/nucleon)入射能条件下在14C、36S和58Ni靶上的(d, p) 反应,发现采用微观光学势[64]得到的单粒子强度比用CH89[58]等唯象的光学势更接近于(e, e′p)反应的系统学行为。因此认为微观光学势更适合于转移反应的理论计算。由于微观光学势不是通过拟合弹性散射实验数据得到的,其内部或许具有更好的微观物理基础。当然,内部光学势的不确定性对于理论计算的转移反应截面的影响仍然是一个值得研究的问题。$ (7) $ 氘核的结合能很低,它在与原子核碰撞过程中很容易被激发到其连续能级,从而发生破裂反应。破裂反应道会与弹性散射、非弹散射、转移反应等其它反应道耦合,改变这些反应道的扭曲波波函数,从而影响它们的截面。作为一个最简单的三体核反应,氘核削裂反应一直是核反应机制研究的重要对象,其中对氘核破裂道的理论处理一直是其中的一个关键问题。我们接下来将重点讨论这一问题。 -
氘核削裂反应的跃迁振幅(式(12))中的
$ \varPsi_\alpha $ 是入射道d(p+n)+A三体体系的总波函数,它满足式(5)中的薛定谔方程。对于氘核削裂反应,在忽略$ {\rm A} $ 核的内部坐标自由度的情况下,式(5)的具体形式为(坐标关系见图1,令$ {\boldsymbol{R}} = {\boldsymbol{r}}_\alpha $ 、$ {\boldsymbol{r}} = {\boldsymbol{r}}_2 $ )$$ (T_{{\boldsymbol{R}}}+H_{{\boldsymbol{r}}}+V_{\rm pA}+V_{\rm nA})\varPsi({\boldsymbol{R}},{\boldsymbol{r}}) = E\varPsi({\boldsymbol{R}},{\boldsymbol{r}}), $$ (29) 其中
$H_{{\boldsymbol{r}}} = T_{{\boldsymbol{r}}}+V_{\rm np}$ 为n-p子系统的哈密顿量。对应的n-p子系统的薛定谔方程为$$ H_{{\boldsymbol{r}}}\phi({\boldsymbol{r}}) = (T_{{\boldsymbol{r}}}+V_{\rm np}(r))\phi({\boldsymbol{r}}) = \epsilon\phi({\boldsymbol{r}}){\text{。}} $$ (30) $ \phi({\boldsymbol{r}}) $ 既可以是n-p体系的束缚态(氘核),也可以是它的连续态。其满足的径向方程为$$ \left\{\frac{\hbar^2}{2\mu}\left[-\frac{{\rm{d}}^2}{{\rm{d}}r^2}+\frac{\ell(\ell+1)}{r^2}\right]+V_{\rm np}(r)\right\}\phi_\ell(r,\epsilon) = \epsilon\phi_\ell(r,\epsilon){\text{。}} $$ (31) 其中
$ \mu $ 是n-p体系的约化质量,$ \ell $ 为n-p之间的轨道角动量。氘核的束缚能很低,在与弹核的碰撞过程中很容易被激发到连续能级,因此,氘核的破裂反应道对氘核削裂反应的影响是一个不可忽略的因素。连续态离散化耦合道方法(Continuum Discretized Coupled Channels, CDCC)[76-77]是目前发展得最完善的处理氘核连续态耦合效应的方法。 -
方程(31)中氘核的连续态波函数之间的关系为
$$ \int \phi_\ell^{\ast}(r,k_1)\phi_\ell(r,k_2)r^2{\rm{d}}r = \delta(k_1-k_2){\text{。}} $$ (32) 显然它们不能构成一个完备正交基。连续态离散化耦合道方法通过构造如下的连续态分段函数(continuum bins)把氘核的连续态离散化[78]:
$$ \phi_{\ell,i}(r) = \sqrt{\frac{2}{\pi N_i}}\int\nolimits_{k_{i-1}}^{k_i} g_i(k)\phi_{\ell}(r,k){\rm{d}}k, $$ (33) 其中
$ g_i(k) $ 是积分的权重函数,通常取为1,$ N_i = \int_{k_{i-1}}^{k_i} $ $ |g_i(k)|^2{\rm{d}}k $ 是分段函数的归一化系数。n-p之间的相对动量$ k $ 和相对运动动能$ \epsilon $ 之间的关系为$ k = \sqrt{2\mu \epsilon/\hbar^2} $ 。可以证明,在能量范围
$ (\epsilon_{i-1},\epsilon_i) $ 和$ (\epsilon_{j-1},\epsilon_j) $ 不重叠的情况下,如此构造的连续态分段波函数满足如下正交关系:$$ \langle\phi_{\ell_1,i}|\phi_{\ell_2,j}\rangle = \int \phi^\ast_{\ell_1,i}\phi_{\ell_2,j}{\rm{d}}r = \delta_{\ell_1\ell_2}\delta_{ij}{\text{。}} $$ (34) 因此,这些分段波函数和氘核的基态波函数一起可以构成一个平方可积的基。三体波函数
$ \varPsi({\boldsymbol{R}},{\boldsymbol{r}}) $ 可以在这个在基下展开为$$ \varPsi({\boldsymbol{R}},{\boldsymbol{r}}) = \varPsi^{\rm{CDCC}}({\boldsymbol{R}},{\boldsymbol{r}}) = \sum\limits_{\ell = 0}^{\ell_{\max}}\sum\limits_{i = 0}^{N_\ell} \phi_{\ell,i}({\boldsymbol{r}})\chi_{\ell,i}({\boldsymbol{R}}) , $$ (35) 其中
$ N_\ell $ 是对每一个分波$ \ell $ 的连续态构造的分段函数的数目,它们由每一个分波$ \ell $ 的连续态的最高能量$ \epsilon_{\max,\ell} $ 以及分段波函数的能量间隔的大小决定$ N_\ell = N_\ell(\epsilon_{\max,\ell},\varDelta\epsilon) $ 。把式(35)代入式(29),即可得到扭曲波
$ \chi_{\ell,i} $ 所满足的耦合道方程组:$$ \begin{split} (E-T_R-\epsilon_i)\chi_{\ell,i}({\boldsymbol{R}}) =& \sum\limits_{\ell' = 0}^{\ell'_{\max}}\sum\limits_j^{N_{\ell'}} \big\langle\phi_{\ell,i}|V_{\rm nA}+\\& V_{\rm pA}|\phi_{\ell',j}\big\rangle \chi_{\ell',j}({\boldsymbol{R}}) {\text{。}}\\[-13pt] \end{split} $$ (36) 同样地,把式(35)代入式(12)即可在氘核削裂反应中包含氘核破裂反应道的贡献。对于氘核削裂反应,式(12)中的
$ \beta\neq\alpha $ ,因此有:$$ \begin{split} T_{\rm dp}^{\rm{CDCC}} =& \sum\limits_{\ell = 0}^{\ell_{\max}}\sum\limits_{i = 0}^{N_\ell} \big\langle \chi_{\rm p}^{(-)}\varPhi_{J_BM_B}({\boldsymbol{B}})|V_\beta-U_\beta|\times \\ & \phi_{\ell,i}({\boldsymbol{r}})\chi_{\ell,i}^{(+)}({\boldsymbol{R}}) \varPhi_{J_AM_A}({\boldsymbol{A}})\big\rangle{\text{。}} \end{split} $$ (37) 可以看到,此时n-p体系的基态和连续态都对转移反应截面有贡献。
CDCC方法已经被广泛地应用于氘核削裂反应理论分析中[11, 14, 79]。原则上CDCC计算中
$ \ell_{\max} $ 以及$ \epsilon_{\max,\ell} $ 的值要足够大,从而保证数值计算的收敛性(计算结果不再随着它们的增大而改变)。但是,目前CDCC方法在氘核削裂反应中的数值计算并不总能得到收敛的结果。对于重靶以及入射能比较高的情况下尤其如此[11]。此外,庞丹阳等[80]的研究表明,(d, p)反应的理论计算的收敛性要求$ \epsilon_{\max} $ 不小于20 MeV。这意味着对入射能小于20 MeV的条件下的氘核削裂反应的CDCC理论分析有必要引入闭合道(closed channel)的贡献。国内早在2008年,安海霞等[81-83]就对闭合道在氘核破裂反应中的作用进行了研究。因此,CDCC方法在氘核削裂反应上的应用仍然有待进一步改进。CDCC方法是处理氘核(弱束缚的入射核)连续态的有效方法。但是它的计算量比较大,并且在转移反应理论计算中不总能得到收敛的结果。因此,可近似处理连续态效应同时具有较高数值稳定性的方法,对于氘核削裂反应的分析是有利的。其中一个代表性的并且被广泛应用至今的方法是Johnson等[84]发展的绝热波近似理论 (adiabatic wave approximation, ADWA)。此外,人们还尝试用其它函数(如Weinberg state,或称Sturmian state)来展开转移反应的三体波函数[85-86]。庞丹阳等[80]发现,在氘核削裂反应中只需要取Weinberg state展开的第一项就可以得到和CDCC非常接近的结果。
-
当入射能在几十MeV以下的条件下,氘核的s-波在氘核削裂反应中起主要作用。如果只考虑n-p子系统的s-波相互作用以及s-波波函数,d(p+n)-A系统的三体波函数可以写成如下形式[84]:
$$ \varPsi^{(+)}({\boldsymbol{R}},{\boldsymbol{r}}) = \phi_0({\boldsymbol{r}})\chi_0({\boldsymbol{R}}) + \int {\rm{d}}{\boldsymbol{k}}\phi^{(+)}(\epsilon_k,{\boldsymbol{r}})\chi(\epsilon_k,{\boldsymbol{R}}), $$ (38) 其中
$ \phi_0 $ 和$ \phi^{(+)} $ 分别是n-p子系统的基态(氘核)和连续态波函数(连续态的能量为$ \epsilon_k $ ):$$ \begin{split} & H_r\phi_0 = -\epsilon_0\phi_0, \quad (\epsilon_0 = -2.226 \;{\rm{MeV}})\\ &H_r\phi^{(+)}(\epsilon_k) = \epsilon_k\phi^{(+)}(\epsilon_k), \quad \epsilon_k\geqslant 0{\text{。}} \end{split}$$ (39) $ \phi_0 $ 和$ \phi^{(+)}(\epsilon_k) $ 满足如下关系$$ \big\langle \phi_0|\phi^{(+)}(\epsilon_k)\big\rangle = 0, \quad \big\langle \phi^{(+)}(\epsilon_k)|\phi^{(+)}(\epsilon_{k'})\big\rangle = \delta({\boldsymbol{k}}-{\boldsymbol{k}}'), $$ (40) $$ |\phi_0\rangle\langle\phi_0|+ \int {\rm{d}}{\boldsymbol{k}}|\phi^{(+)}(\epsilon_k)\rangle\langle\phi^{(+)}(\epsilon_{k})| = 1{\text{。}} $$ (41) 这里,
${\boldsymbol{k}}^2 = 2\mu_{\rm np}\epsilon_{k}$ ,${\boldsymbol{k'}}^2 = 2\mu_{\rm np}\epsilon_{k'}$ 。把式(38)形式的三体波函数代入三体薛定谔方程,在
$ r\sim 0 $ 处(由于$V_{\rm np}$ 是短程相互作用,转移反应截面主要来自于$ r\sim 0 $ 附近的跃迁振幅贡献),有$$ \begin{split}& \big[E+\epsilon_0-T_R-\bar{V}(R)\big]\phi_0(0)\chi_0(R)+\\ &\int {\rm{d}}{\boldsymbol{k}}\big[E-\epsilon_k-T_R-\bar{V}(R)\big]\phi^{(+)}(\epsilon_k,0)\chi(\epsilon_k,R) = 0, \end{split}$$ (42) 其中
$\bar{V} = V_{\rm nA}+V_{\rm pA}$ 。注意到这个公式里两个中括号里面的差别在于$ \epsilon_k $ 和$ \epsilon_0 $ 。Johnson和Soper[84]提出的ADWA的要点在于把所有的连续态的能量$ \epsilon_k $ 用氘核的束缚态能量$ \epsilon_0 $ 来代替,此时方程变为$$ [E+\epsilon_0-T_R-\bar{V}]\bar{\chi}(R) = 0, $$ (43) 其中
$ \bar{\chi}(R) $ 是等效的氘核的扭曲波:$$ \bar{\chi}(R) = \chi_0(R)+\int {\rm{d}}{\boldsymbol{k}} \phi^{(+)}(\epsilon_k,0)\chi(\epsilon_k,R)/\phi_0(0), $$ (44) 而
$ \bar{V} $ 则是等效的d-A相互作用。从这里面可以看到,这个等效的扭曲波除了包含氘核基态所对应的扭曲波成分,还包含了氘核的连续态的贡献。用此等效氘核波函数代替式(27)中的$\chi_{\rm d}^{(+)}$ ,即得到绝热近似条件下的氘核削裂反应的跃迁振幅。$ \bar{\chi} $ 可以通过求解方程(43)得到,其中光学势$V_{\rm nA}$ 和$V_{\rm pA}$ 对应的中子和质子的入射能一般取作氘核入射能的一半($E_{\rm d}/2$ -rule)。Timofeyuk等[67]的研究表明,如果考虑了质子和中子光学势的非定域性,$E_{\rm d}/2$ -rule将需要被修正,详见文献[67]。值得指出的是,式(43)中的氘核的有效光学势$ \bar{V} $ 和式(28)中的$U_{\rm dA}$ 不同,前者只适用于氘核削裂反应的计算,理论上不能描述氘核的弹性散射,而后者是氘核的光学势,通常要求它可以准确地描述氘核与A核的弹性散射。Johnson-Soper提出的ADWA是零程的($ r\sim 0 $ ),并且只包含了氘核的s-波的贡献。有限程版本的ADWA见文献[85]。文献[87]发现,当氘核能量在20~30 MeV/nucleon以下时,有限程效应对ADWA计算的影响较小(
$ < 10\% $ ),但随着入射能的增加,有限程效应会越来越重要。不过,对于(d, p) 和(p, d)反应的分析,Johnson-Soper势的最简形式($ \bar{V} = V_{\rm nA}+V_{\rm pA} $ )迄今为止仍然得到最广泛的应用。在ADWA中,n-p体系的自由度被冻结(假定所有激发态的能量都等于氘核基态的结合能)。因此,ADWA要求束流能量应远高于对氘核削裂反应有贡献的氘核激发态的能量(式(35)中的$ \epsilon_{\max} $ )。文献[88-91]的计算表明,ADWA在$ E_{\rm d}\leqslant 20 $ MeV时对氘核削裂反应的n-p破裂贡献的描述不好(详见综述[76])。因此,人们期待ADWA在束流能量$ E_{\rm d}\gg $ $ 20 $ MeV时才足够精确。系统学研究[92-94]表明,ADWA在束流能量较高的(d, p)/(p, d)反应研究中优于DWBA。ADWA理论的数值计算和DWBA理论一样简单,同时又通过引入氘核与靶核的有效相互作用势在一定程度上包含了氘核破裂道的影响,因此它至今仍然被广泛地应用于(d, p)/(p, d)反应的理论分析中。 -
氘核削裂反应(A(d, p)B)是一个典型的三体反应。其中p, n, A三个集团共有三套雅可比坐标(见图4)。前面介绍的反应模型中的三体波函数只是一个雅可比坐标下的三体薛定谔方程的解。例如,式(12)中的
$ \varPsi_\alpha^{(+)} $ 只对应图4-(1)这一个雅可比坐标,因此,这个波函数在渐近区域只包含氘核的弹性散射和破裂反应道的成分。而完整的三体波函数应该包含所有的这三个雅可比坐标的成分:$$ \varPsi = \sum\limits_{n = 1}^3 \varPsi^{(n)}({\boldsymbol{r}}_n,{\boldsymbol{R}}_n) , $$ (45) 基于法捷耶夫方程组(Faddeev equations)的少体核反应理论允许我们在三体反应理论分析中同时包含这三个雅可比坐标的成分[95-98]。
法捷耶夫方程组由三个方程组成,它们与波函数的三个组分通过以下方式耦合(对于图4中的三体体系):
$$ \begin{split} (E-T_1-V_{\rm pn})\varPsi^{(1)} =& V_{\rm pn}(\varPsi^{(2)}+\varPsi^{(3)}) , \\ (E-T_2-V_{\rm pA})\varPsi^{(2)} =& V_{\rm pA}(\varPsi^{(3)}+\varPsi^{(1)}) , \\ (E-T_3-V_{\rm nA})\varPsi^{(3)} =& V_{\rm nA}(\varPsi^{(1)}+\varPsi^{(2)}){\text{。}} \end{split} $$ (46) 其中
$ T_1 $ 、$ T_2 $ 和$ T_3 $ 分别为三个雅可比坐标下的动能算符。可以看出,法捷耶夫方程能够在同一框架下对d(p+n)-A三体体系的所有反应道进行完备的描述.这些反应道包括氘核的弹性散射、破裂反应、(d, p)反应,以及(d, n)反应。2007年,Alt等[95]采用可分离化的光学势以及n-p相互作用,首次实现对d-12C体系的三体反应理论计算。目前,法捷耶夫方程在重核的氘核削裂反应中的应用面临的主要问题是对库仑相互作用的处理。通常的法捷耶夫方程组是在动量空间求解的,但包含长程库仑力的方程在动量空间无法直接求解。目前对库仑力的处理采用的是屏蔽势的方法[95, 97-98]。这种方法对于
$ Z\lesssim 20 $ 的靶核大致是适用的,但是在更重的靶核上难以给出收敛的结果[99]。2012年,Mukhamedzhanov[96]提出了一个推广的法捷耶夫方程方法,采用库仑波函数基求解法捷耶夫方程,从而克服带库仑势的求解问题。但是这个理论尚没有成功的数值实现。
-
摘要: 氘核削裂反应作为单核子转移反应的一种,在原子核的单粒子结构研究中一直起着重要的作用。同时,作为典型的三体核反应,它也是直接核反应理论研究的一个重要研究对象。本文简述了氘核削裂反应的研究历史,侧重介绍其反应理论的演进,以及平面波玻恩近似(Plane Wave Born Approximation,PWBA)、扭曲波玻恩近似(Distorted Wave Born Approximation,DWBA)、连续态离散化耦合道理论(Continuum Discretized Coupled Channels,CDCC)、绝热波近似理论(Adiabatic Wave Approximation,ADWA)和法捷耶夫方程组(Faddeev equations)方法之间的差别与联系。本文讨论了当前及未来氘核削裂反应理论研究的若干热点问题,包括中高能条件下的氘核削裂反应、转移到剩余核连续态的理论、色散光学势以及非定域光学势在氘核削裂反应中的应用等。本文展示了我们对中高能条件下的氘核削裂反应理论计算的一些初步成果。Abstract: As one type of single nucleon transfer reactions, deuteron stripping reaction has been playing an important role in the study of the single particle structure of atomic nuclei. Also, as a typical three-body nuclear reaction, deuteron stripping reaction is also an important subject in the development of direct nuclear reaction theories. We briefly review the history of the study of deuteron stripping reactions, focusing on the evolution of its theoretical descriptions, and the differences and the relations of the Plane Wave Born Approximation(PWBA), the distorted wave Born approximation(DWBA), the Continuum Discretized Coupled Channel method(CDCC), the Adiabatic Wave Approximation(ADWA), and the approach based on the Faddeev equations. Some important problems being currently addressed in deuteron stripping reaction theories and those should be addressed in near future are also discussed, including deuteron strippings at intermediate and high energies, theory for transfer to continuum states of the residue nuclei, applications of dispersive optical model potentials and non-local potentials, etc. Some preliminary results of our theoretical calculations for deuteron stripping at intermediate and high energies are presented.
-
图 3 PWBA和DWBA对40Ca(d, p)41Ca反应预测的比较,
$E_{\rm d}=12$ MeV。实线和虚线分别对应DWBA和PWBA的计算结果,实心圆点表示实验数据。此图来自文献[50]图 7 (在线彩图) 入射能量分别为45.34 MeV(a)、392 MeV(b)和800 MeV(c)时,16O(
${\rm p, d}$ )15O的正宇称激发态($\frac{5}{2}^{+}$ ,5.241 MeV)截面的理论计算与实验数据的比较(图标同图6)图 8 (在线彩图) 入射能量分别为45.34 MeV(a)、392 MeV(b)和800 MeV(c)时,16O(
${\rm p, d}$ )15O的正宇称激发态($\frac{3}{2}^{-}$ ,6.178 MeV)截面的理论计算与实验数据的比较(图标同图6) -
[1] LEE J, TOSTEVIN J A, BROWN B A, et al. Phys Rev C, 2006, 73: 044608. doi: 10.1103/PhysRevC.73.044608 [2] XU Y P, PANG D Y, YUN X Y, et al. Phys Rev C, 2018, 98: 044622. doi: 10.1103/PhysRevC.98.044622 [3] CHEN J, LOU J L, YE Y L, et al. Phys Lett B, 2018, 781: 412. doi: 10.1016/j.physletb.2018.04.016 [4] CHEN J, LOU J L, YE Y L, et al. Phys Rev C, 2018, 98: 014616. doi: 10.1103/PhysRevC.98.014616 [5] XU Y P, PANG D Y, YUN X Y, et al. Phys Lett B, 2019, 790: 308. [6] ATAR L, PASCHALIS S, BARBIERI C, et al. Phys Rev Lett, 2018, 120: 052501. doi: 10.1103/PhysRevLett.120.052501 [7] JIANG Y, LOU J L, YE Y L, et al. Chin Phys Lett, 2018, 35(8): 082501. doi: 10.1088/0256-307x/35/8/082501 [8] LIU W, LOU J L, YE Y L, et al. Nucl Sci Tech, 2020, 31: 20. doi: 10.1007/s41365-020-0731-y [9] CHEN J, WANG S M, FORTUNE H T, et al. Phys Rev C, 2021, 103: L031302. doi: 10.1103/PhysRevC.103.L031302 [10] MUKHAMEDZHANOV A M, BÉM P, BROWN B A, et al. Phys Rev C, 2003, 67: 065804. doi: 10.1103/PhysRevC.67.065804 [11] PANG D Y, MUKHAMEDZHANOV A M. Phys Rev C, 2014, 90: 044611. doi: 10.1103/PhysRevC.90.044611 [12] PANDHARIPANDE V R, SICK I, HUBERTS P K A D. Rev Mod Phys, 1997, 69: 981. doi: 10.1103/RevModPhys.69.981 [13] BARBIERI C. Phys Rev Lett, 2009, 103: 202502. doi: 10.1103/PhysRevLett.103.202502 [14] ONG H J, TANIHATA I, TAMⅡ A, et al. Phys Lett B, 2013, 725(4): 277. [15] TRIBBLE R E, BERTULANI C A, LA COGNATA M, et al. Reports on Progress in Physics, 2014, 77(10): 106901. doi: 10.1088/0034-4885/77/10/106901 [16] MUKHAMEDZHANOV A M, PANG D Y, KADYROV A S. Phys Rev C, 2019, 99: 064618. doi: 10.1103/PhysRevC.99.064618 [17] SHEN Y P, GUO B, LI Z H, et al. Phys Rev C, 2019, 99: 025805. doi: 10.1103/PhysRevC.99.025805 [18] JONES K L, ADEKOLA A S, BARDAYAN D W, et al. Nature, 2010, 465: 454. doi: 10.1038/nature09048 [19] SCHMITT K T, JONES K L, BEY A, et al. Phys Rev Lett, 2012, 108: 192701. doi: 10.1103/PhysRevLett.108.192701 [20] KOZUB R L, ARBANAS G, ADEKOLA A S, et al. Phys Rev Lett, 2012, 109: 172501. doi: 10.1103/PhysRevLett.109.172501 [21] SCHIFFER J P, HOFFMAN C R, KAY B P, et al. Phys Rev Lett, 2012, 108: 022501. doi: 10.1103/PhysRevLett.108.022501 [22] KAY B P, SCHIFFER J P, FREEMAN S J. Phys Rev Lett, 2013, 111: 042502. doi: 10.1103/PhysRevLett.111.042502 [23] BARDAYAN D W. Journal of Physics G: Nuclear and Particle Physics, 2016, 43(4): 043001. doi: 10.1088/0954-3899/43/4/043001 [24] 李根, 谭智威, 楼建玲, 等. 原子核物理评论, 2020, 37(3): 426. doi: 10.11804/NuclPhysRev.37.2019CNPC09 LI G, TAN Z W, LOU J L, et al. Nuclear Physics Review, 2020, 37(3): 426. (in Chinese) doi: 10.11804/NuclPhysRev.37.2019CNPC09 [25] LAWRENCE E O, LIVINGSTON M S, LEWIS G N. Phys Rev, 1933, 44: 56. doi: 10.1103/PhysRev.44.56 [26] OPPENHEIMER J R, PHILLIPS M. Phys Rev, 1935, 48: 500. doi: 10.1103/PhysRev.48.500 [27] SERBER R. Phys Rev, 1947, 72: 1008. doi: 10.1103/PhysRev.72.1008 [28] KEMPTON A E, BROWNE B C, MAASDORP R, et al. Proceedings of the Royal Society of London Series A - Mathematical and Physical Sciences, 1936, 157(891): 386. doi: 10.1098/rspa.1936.0202 [29] KONOPINSKI E J, TELLER E. Phys Rev, 1948, 73: 822. doi: 10.1103/PhysRev.73.822 [30] WOLFENSTEIN L. Phys Rev, 1949, 75: 1664. doi: 10.1103/PhysRev.75.1664 [31] BEIDUK F M, PRUETT J R, KONOPINSKI E J. Phys Rev, 1950, 77: 622. doi: 10.1103/PhysRev.77.622 [32] PRUETT J R, BEIDUK F M, KONOPINSKI E J. Phys Rev, 1950, 77: 628. doi: 10.1103/PhysRev.77.628 [33] MAYER M G. Phys Rev, 1949, 75: 1969. doi: 10.1103/PhysRev.75.1969 [34] HAXEL O, JENSEN J H D, SUESS H E. Phys Rev, 1949, 75: 1766. doi: 10.1103/PhysRev.75.1766.2 [35] BURROWS H B, GIBSON W M, ROTBLAT J. Phys Rev, 1950, 80: 1095. doi: 10.1103/PhysRev.80.1095 [36] BUTLER S T. Phys Rev, 1950, 80: 1095. doi: 10.1103/PhysRev.80.1095.2 [37] BUTLER S T, SYMONDS J L. Phys Rev, 1951, 83: 858. doi: 10.1103/PhysRev.83.858 [38] ROTBLAT J. Nature, 1951, 167: 1027. doi: 10.1038/1671027a0 [39] PRESTON R L, MARTIN H J, SAMPSON M B. Phys Rev, 1961, 121: 1741. doi: 10.1103/PhysRev.121.1741 [40] KULL L A. Phys Rev, 1967, 163: 1066. doi: 10.1103/PhysRev.163.1066 [41] KULL L A, KASHY E. Phys Rev, 1968, 167: 963. doi: 10.1103/PhysRev.167.963 [42] FICK D, WEISS U. Zeitschrift für Physik A Hadrons and nuclei, 1973, 265(1): 87. doi: 10.1007/BF01397730 [43] GLENDENNING N K. Direct Nuclear Reactions[M]. Singapore: World Scientific, 2004. [44] HOROWITZ J, MESSIAH A M L. Phys Rev, 1953, 92: 1326. doi: 10.1103/PhysRev.92.1326.2 [45] TOBOCMAN W. Phys Rev, 1954, 94: 1655. doi: 10.1103/PhysRev.94.1655 [46] BUCK B, HODGSON P E. The Philosophical Magazine: A Journal of Theoretical Experimental and Applied Physics, 1961, 6(71): 1371. doi: 10.1080/14786436108241232 [47] HINDS S, MIDDLETON R, PULLEN D J. Phys Lett B, 1961, 1: 12. doi: 10.1016/0031-9163(62)90261-5 [48] SATCHLER G R. Nucl Phys, 1964, 55: 1. doi: 10.1016/0029-5582(64)90124-5 [49] PHILPOTT R J, PINKSTON W T, SATCHLER G R. Nucl Phys A, 1968, 119: 241. doi: 10.1016/0375-9474(68)90300-X [50] SATCHLER G R. Direct Nuclear Reactions[M]. New York: Oxford University Press Inc., 1983. [51] TAMURA T. Physics Reports, 1974, 14: 59. Physics Reports, 1974, 14: 59. [52] STAMP A. Nuclear Physics, 1966, 83(1): 232.Nuclear Physics, 1966, 83(1): 232. [53] RAWITSCHER G H. Phys Rev, 1967, 163: 1223. doi: 10.1103/PhysRev.163.1223 [54] AUSTERN N. Phys Rev, 1969, 188: 1595. doi: 10.1103/PhysRev.188.1595 [55] ASCUITTO R J, GLENDENNING N K, SØRENSEN B. Nuclear Physics A, 1971, 170(1): 65. [56] IMANISHI B, ICHIMURA M, KAWAI M. . Phys Lett B, 1974, 52 (3): 267. [57] COTANCH S R, VINCENT C M. Phys Rev C, 1976, 14: 1739. doi: 10.1103/PhysRevC.14.1739 [58] VARNER R L, THOMPSON W J, MCABEE T L, et al. Physics Reports, 1991, 201(2): 57. [59] KONING A J, DELAROCHE J P. Nuclear Physics A, 2003, 713 (3): 231. [60] HAN Y, SHI Y, SHEN Q. Phys Rev C, 2006, 74: 044615. doi: 10.1103/PhysRevC.74.044615 [61] XU R, MA Z, VAN DALEN E N E, et al. Phys Rev C, 2012, 85: 034613. doi: 10.1103/PhysRevC.85.034613 [62] PANG D Y, DEAN W M, MUKHAMEDZHANOV A M. Phys Rev C, 2015, 91: 024611. doi: 10.1103/PhysRevC.91.024611 [63] ZHANG Y, PANG D Y, LOU J L. Phys Rev C, 2016, 94: 014619. doi: 10.1103/PhysRevC.94.014619 [64] XU R, MA Z, ZHANG Y, et al. Phys Rev C, 2016, 94: 034606. doi: 10.1103/PhysRevC.94.034606 [65] LIU X D, FAMIANO M A, LYNCH W G, et al. Phys Rev C, 2004, 69: 064313. doi: 10.1103/PhysRevC.69.064313 [66] LOVELL A E, NUNES F M, SARICH J, et al. Phys Rev C, 2017, 95: 024611. doi: 10.1103/PhysRevC.95.024611 [67] TIMOFEYUK N K, JOHNSON R C. Phys Rev Lett, 2013, 110: 112501. doi: 10.1103/PhysRevLett.110.112501 [68] TITUS L J, NUNES F M. Phys Rev C, 2014, 89: 034609. doi: 10.1103/PhysRevC.89.034609 [69] TITUS L J, NUNES F M, POTEL G. Phys Rev C, 2016, 93: 014604. doi: 10.1103/PhysRevC.93.014604 [70] GÓMEZ-RAMOS M, TIMOFEYUK N K. Phys Rev C, 2018, 98: 011601. doi: 10.1103/PhysRevC.98.011601 [71] TIMOFEYUK N K, JOHNSON R C.Progress in Particle and Nuclear Physics, 2020, 111: 103738. [72] GÓMEZ-RAMOS M, MORO A M, GÓMEZ-CAMACHO J, et al. Phys Rev C, 2015, 92: 014613. doi: 10.1103/PhysRevC.92.014613 [73] THOMPSON I J. Computer Physics Reports, 1988, 7(4): 167. [74] TOSTEVIN J A. University of Surrey Version of the Code Twofnr (of m. toyama, m. igarashi and n. kishida) and Code Front (Private Communication)[Z]. 1977. [75] YUN X Y, PANG D Y, XU Y P, et al. Sci China Phys Mech Astron, 2020, 63: 222011. doi: 10.1007/s11433-019-9389-6 [76] AUSTERN N, ISERI Y, KAMIMURA M, et al. Physics Reports, 1987, 154(3): 125. [77] YAHIRO M, OGATA K, MATSUMOTO T, et al. Progress of Theoretical and Experimental Physics, 2012, 2012(1): 01A206. doi: 10.1093/ptep/pts008 [78] TOSTEVIN J A, NUNES F M, THOMPSON I J. Phys Rev C, 2001, 63: 024617. doi: 10.1103/PhysRevC.63.024617 [79] MUKHAMEDZHANOV A M, PANG D Y, BERTULANI C A, et al. Phys Rev C, 2014, 90: 034604. doi: 10.1103/PhysRevC.90.034604 [80] PANG D Y, TIMOFEYUK N K, JOHNSON R C, et al. Phys Rev C, 2013, 87: 064613. doi: 10.1103/PhysRevC.87.064613 [81] AN H X, CAI C H. Communications in Theoretical Physics, 2008, 50(3): 719. doi: 10.1088/0253-6102/50/3/39 [82] 安海霞, 蔡崇海. 原子核物理评论, 2008, 25(3): 224. doi: 10.11804/NuclPhysRev.25.03.224 AN H X, CAI C H. Nuclear Physics Review, 2008, 25(3): 224. (in Chinese) doi: 10.11804/NuclPhysRev.25.03.224 [83] 安海霞, 蔡崇海. 原子核物理评论, 2010, 27(4): 411. doi: 10.11804/NuclPhysRev.27.04.411 AN H X, CAI C H. Nuclear Physics Review, 2010, 27(4): 411. (in Chinese) doi: 10.11804/NuclPhysRev.27.04.411 [84] JOHNSON R C, SOPER P J R. Phys Rev C, 1970, 1: 976. doi: 10.1103/PhysRevC.1.976 [85] JOHNSON R C, TANDY P C. Nuclear Physics A, 1974, 235(1): 56. doi: 10.1016/0375-9474(74)90178-X [86] LAID A, TOSTEVIN J A, JOHNSON R C. Phys Rev C, 1993, 48: 1307. doi: 10.1103/PhysRevC.48.1307 [87] NGUYEN N B, NUNES F M, JOHNSON R C. Phys Rev C, 2010, 82: 014611. doi: 10.1103/PhysRevC.82.014611 [88] RAWITSCHER G H. Phys Rev C, 1975, 11: 1152. doi: 10.1103/PhysRevC.11.1152 [89] ISERI Y, YAHIRO M, NAKANO M. Progress of Theoretical Physics, 1983, 69(3): 1038. doi: 10.1143/PTP.69.1038 [90] AMAKAWA H, AUSTERN N. Phys Rev C, 1983, 27: 922. doi: 10.1103/PhysRevC.27.922 [91] KAWAI M. Progress of Theoretical Physics Supplement, 1986, 89: 11. doi: 10.1143/PTPS.89.11 [92] HARVEY J D, JOHNSON R C. Phys Rev C, 1971, 3: 636. doi: 10.1103/PhysRevC.3.636 [93] SATCHLER G R. Phys Rev C, 1971, 4: 1485. doi: 10.1103/PhysRevC.4.1485 [94] WALES G L, JOHNSON R C. Nuclear Physics A, 1976, 274(1): 168. [95] ALT E O, BLOKHINTSEV L D, MUKHAMEDZHANOV A M, et al. Phys Rev C, 2007, 75: 054003. doi: 10.1103/PhysRevC.75.054003 [96] MUKHAMEDZHANOV A M, EREMENKO V, SATTAROV A I. Phys Rev C, 2012, 86: 034001. doi: 10.1103/PhysRevC.86.034001 [97] DELTUVA A. Phys Rev C, 2015, 91: 024607. doi: 10.1103/PhysRevC.91.024607 [98] DELTUVA A, ROSS A, NORVAIŠAS E, et al. Phys Rev C, 2016, 94: 044613. doi: 10.1103/PhysRevC.94.044613 [99] UPADHYAY N J, DELTUVA A, NUNES F M. Phys Rev C, 2012, 85: 054621. doi: 10.1103/PhysRevC.85.054621 [100] FRANEY M A, LOVE W G. Phys Rev C, 1985, 31: 488. doi: 10.1103/PhysRevC.31.488 [101] OGATA K. Private Communications. [102] VINCENT C M, FORTUNE H T. Phys Rev C, 1970, 2: 782. doi: 10.1103/PhysRevC.2.782 [103] KADYROV A S, BRAY I, MUKHAMEDZHANOV A M, et al. Annals of Physics, 2009, 324(7): 1516. doi: 10.1016/j.aop.2009.02.003 [104] MUKHAMEDZHANOV A M. Phys Rev C, 2011, 84: 044616. doi: 10.1103/PhysRevC.84.044616 [105] ESCHER J E, THOMPSON I J, ARBANAS G, et al. Phys Rev C, 2014, 89: 054605. doi: 10.1103/PhysRevC.89.054605 [106] CHARITY R J, SOBOTKA L G, DICKHOFF W H. Phys Rev Lett, 2006, 97: 162503. doi: 10.1103/PhysRevLett.97.162503 [107] MUELLER J M, CHARITY R J, SHANE R, et al. Phys Rev C, 2011, 83: 064605. doi: 10.1103/PhysRevC.83.064605 [108] LI R, SUN W L, SOUKHOVITSKIĨ E S, et al. Phys Rev C, 2013, 87: 054611. doi: 10.1103/PhysRevC.87.054611 [109] ZHAO X N, SUN W L, SOUKHOVITSKIĨ E S, et al. Journal of Physics G: Nuclear and Particle Physics, 2019, 46(5): 055103. doi: 10.1088/1361-6471/ab0010 [110] ZHAO X N, SUN W L, CAPOTE R, et al. Phys Rev C, 2020, 101: 064618. doi: 10.1103/PhysRevC.101.064618 [111] NGUYEN N B, WALDECKER S J, NUNES F M, et al. Phys Rev C, 2011, 84: 044611. doi: 10.1103/PhysRevC.84.044611 [112] PEREY F, BUCK B. Nuclear Physics, 1962, 32: 353. [113] GIANNINI M, RICCO G. Annals of Physics, 1976, 102(2): 458. [114] GIANNINI M, RICCO G, ZUCCHIATTI A. Annals of Physics, 1980, 124(1): 208. [115] PEREY F. Direct Interactions and Nuclear Reaction Mechanisms[M]. New York: Gordon and Breach Science Publishers, Inc, 1964. [116] MAHZOON M H, CHARITY R J, DICKHOFF W H, et al. Phys Rev Lett, 2014, 112: 162503. doi: 10.1103/PhysRevLett.112.162503 [117] TIAN Y, PANG D Y, MA Z Y. Int J Mod Phys E, 2015, 24: 1550006. doi: 10.1142/S0218301315500068 [118] LOVELL A E, BACQ P L, CAPEL P, et al. Phys Rev C, 2017, 96: 051601. doi: 10.1103/PhysRevC.96.051601 [119] ROSS A, TITUS L J, NUNES F M, et al. Phys Rev C, 2015, 92: 044607. doi: 10.1103/PhysRevC.92.044607 [120] TIAN Y, PANG D Y, MA Z Y. Phys Rev C, 2018, 97: 064615. doi: 10.1103/PhysRevC.97.064615 [121] JACKSON D F, JOHNSON R C. Phys Lett B, 1974, 49(3): 249. [122] JAGHOUB M I, LOVELL A E, NUNES F M. Phys Rev C, 2018, 98: 024609. doi: 10.1103/PhysRevC.98.024609