-
通过核子-核子相互作用重现有限核的性质是核物理的基本任务之一,基于非相对论或相对论框架构造的核力为有限核及核物质性质的计算提供了重要输入量。在非相对论框架下,常见的核子-核子相互作用包括零力程的Skyrme势和有限力程的Gogny势[35-36]。在相对论框架下,基于夸克-胶子自由度构造的核子-核子相互作用在实际应用中存在诸多困难,基于单玻色子交换模型在重子自由度下发展而来的相对论平均场理论以及基于手征对称性发展而来的手征有效核力在近年来得到了广泛关注与应用[37-38]。在量子强子动力学框架下,核子-核子相互作用可以通过交换$ \sigma $、$ \omega $、$ \delta $、$ \rho $、$ \pi $介子来实现。$ \sigma $是自旋标量,同位旋标量介子;$ \omega $是自旋矢量,同位旋标量介子;$ \delta $是自旋标量,同位旋矢量介子;$ \rho $是自旋和同位旋都是矢量的介子,它们共同组成了一套描述核力的交换粒子。在$ {\rm{NN}} $弹性散射研究中[18],相比于$ \rho $介子交换的贡献,$ \pi $介子交换对核子-核子弹性散射截面的贡献可以忽略。因此,在这里我们忽略$ \pi $介子的贡献,并将在后续的工作中进一步分析$ \pi $介子交换对核子-$ \Delta $弹性散射截面的贡献。我们之前的研究结果表明,核子-核子弹性散射截面敏感于核子-核子-介子的耦合常数,并且当$ g_{{\rm{N}} {\rm{N}}}^{\omega}>g_{{\rm{N}} {\rm{N}}}^{\sigma} $,$ g_{{\rm{N}} {\rm{N}}}^{\omega} $$ -g_{{\rm{N}} {\rm{N}}}^{\sigma} $ $ \lesssim $ 2时,理论计算结果可以描述$ \sigma^{\ast}_{{\rm{nn}}({\rm{pp}})} $和$ \sigma^{\ast}_{{\rm{np}}} $的能量依赖行为[18]。由于$ {\rm{N}} \Delta $散射截面与$ {\rm{NN}} $散射截面有着相同的性质,因此$ {\rm{N}} \Delta $散射截面的耦合常数的选择也需要满足$ {\rm{NN}} $散射截面中的限制条件。这里密度依赖的耦合常数可以通过拟合Dirac-Brueckner Hartree Fock(DBHF)自能来得到[39]。
系统的有效拉氏量表示为
$$ {\cal{L}} = {\cal{L}}_{\text{F}}+{\cal{L}}_{\text{I}}, $$ (1) 其中$ {\cal{L}}_{\text{F}} $是重子和介子的自由拉氏量,
$$ \begin{split} {\cal{L}}_{\text{F}} =\;& \bar{\varPsi}\left[{\rm{i}} \gamma_{\mu} \partial^{\mu}-m_{N}\right] \varPsi+\bar{\varPsi}_{\varDelta \nu}\left[{\rm{i}} \gamma_{\mu} \partial^{\mu}-m_{\varDelta}\right] \varPsi_{\varDelta}^{\nu}+ \\ & \dfrac{1}{2} \partial_{\mu} \sigma \partial^{\mu} \sigma+\dfrac{1}{2} \partial_{\mu} {\boldsymbol{\delta}} \partial^{\mu} {\boldsymbol{\delta}}-\dfrac{1}{4} F_{\mu \nu} \boldsymbol\cdot F^{\mu v}-\dfrac{1}{4} {\boldsymbol{L}}_{\mu\nu} \boldsymbol\cdot {\boldsymbol{L}}^{\mu \nu}- \\ & \dfrac{1}{2} m_{\sigma}^{2} \sigma^{2}-\dfrac{1}{2} m_{\delta}^{2} {\boldsymbol{\delta}}^{2} +\dfrac{1}{2} m_{\omega}^{2} \omega_{\mu} \omega^{\mu}+\dfrac{1}{2} m_{\rho}^{2} {\boldsymbol{\rho}}_{\mu} {\boldsymbol{\rho}}^{\mu}, \end{split} $$ (2) 这里,$F_{\mu \nu} \equiv \partial_{\mu} \omega_{v}-\partial_{v} \omega_{\mu} , {\boldsymbol{L}}_{\mu \nu} \equiv \partial_{\mu} {\boldsymbol{\rho}}_{v}-\partial_{v} {\boldsymbol{\rho}}_{\mu}$。相互作用部分$ {\cal{L}}_{\text{I}} $表示为
$$ \begin{split} {\cal{L}}_{\text{I}} =\; & g_{\rm{NN}}^{\sigma} \bar{\varPsi} \varPsi \sigma+g_{\rm{N N}}^{\delta} \bar{\varPsi}{\boldsymbol{\tau}}\boldsymbol\cdot \varPsi {\boldsymbol{\delta}} -g_{\rm{N N}}^{\omega} \bar{\varPsi} \gamma_{\mu} \varPsi \omega^{\mu}- \\ & g_{\rm{N N}}^{\rho} \bar{\varPsi} \gamma_{\mu} {\boldsymbol{\tau}} \boldsymbol\cdot \varPsi {\boldsymbol{\rho}}^{\mu}+g_{\varDelta \varDelta}^{\sigma} \bar{\varPsi}_{\varDelta} \varPsi_{\varDelta} \sigma +g_{\varDelta \varDelta}^{\delta} \bar{\varPsi}_{\varDelta} {\boldsymbol{\tau}} \boldsymbol\cdot \varPsi_{\varDelta} {\boldsymbol{\delta}}- \\ & g_{\varDelta \varDelta}^{\omega} \bar{\varPsi}_{\varDelta} \gamma_{\mu} \varPsi_{\varDelta} \omega^{\mu}-g_{\varDelta \varDelta}^{\rho} \bar{\varPsi}_{\varDelta} \gamma_{\mu} {\boldsymbol{\tau}} \boldsymbol\cdot \varPsi_{\varDelta} {\boldsymbol{\rho}}^{\mu}, \end{split} $$ (3) 这里有效拉氏密度采用密度依赖的耦合常数[39],其参数化形式为
$$ g^{i}_{{\rm{NN}}}(\rho_{b})=g^{i}_{\rm{NN}}(\rho_{0})p_{i}(\xi), \quad i=\sigma、\omega、\delta、\rho $$ (4) 其中 $\xi = {\rho_{b}}/{\rho_{0}}$,$ \rho_{0} $ 是对称核物质的饱和密度,耦合常数中密度依赖的修正部分$ p_{i}(\xi) $可以表达为
$$ p_{i}(\xi) = A_{i}\dfrac{1+B_{i}(\xi+D_{i})^{2}}{1+C_{i}(\xi+E_{i})^{2}}, $$ (5) 密度依赖的耦合常数参数由文献[17, 39]给出。这里,$ \psi $表示自旋为${1}/{2}$的Dirac旋量,$ \psi_{\varDelta} $表示自旋为${3}/{2}$的Rarita-Schwinger场旋量,$ \tau $为同位旋算符。核子质量$ m_{{\rm{N}}} $取939 MeV,$ \Delta $的质量$ m_{\Delta} $取1 232 MeV。关于$ \Delta $-$ \Delta $-介子耦合,由于重子有效质量的劈裂使得SU(6)对称性不能精确成立,因此重子八重态和十重态不能共用相同的耦合常数。为了能够重现饱和密度区截面的数值,耦合常数和重子截断质量需要同时考虑。对于$ \Delta $-$ \Delta $-介子耦合常数的取值,这里引入$ \Delta\Delta $-$ {\rm{NN}} $耦合常数比:
$$ f_{i} = \frac{g^{i}_{\varDelta\varDelta}}{g^{i}_{{\rm{NN}}}}, $$ (6) 其中$f_{\sigma} = 1.0,f_{\omega} = 0.8,f_{\rho} = 1.0$[40],$f_{\delta} = 1.0$。
-
接下来,我们推导$ \Delta $共振态的RBUU方程。对于非平衡问题的处理,闭合时间回路格林函数方法显得更加有效。在相互作用表象中,$ \Delta $格林函数定义为[34]
$$ {\rm{i}}G_{\alpha \beta}(1,2) = \left\langle T\left[\exp \left(-{\rm{i}} \oint {\rm{d}} x H_{I}(x)\right) \psi _{\varDelta \alpha}(1) \bar{\psi}_{\varDelta \beta }(2)\right]\right\rangle。 $$ (7) 在闭时回路框架下的四个$ \Delta $ 格林函数$ G^{--}_{\alpha \beta}(1,2) $、$ G^{-+}_{\alpha \beta}(1,2) $、$ G^{+-}_{\alpha \beta}(1,2) $、$ G^{++}_{\alpha \beta}(1,2) $表示为如下矩阵形式:
$$ G_{\alpha \beta}(1,2) = \left[\begin{array}{*{20}{ll}} G_{\alpha \beta}^{--}(1,2) & G_{\alpha \beta}^{-+}(1,2) \\ G_{\alpha \beta}^{+-}(1,2) & G_{\alpha \beta}^{++}(1,2) \end{array}\right]。 $$ (8) 将式(7)做微扰展开可得:
$$ \begin{split} {{\rm{i}} G(1,2)_{\alpha \beta}} =\;&{{\rm{i}} G_{\alpha \beta}^{0}(1,2)}+\\ &{\oint {\rm{d}} x_{3} \oint {\rm{d}} x_{4} G_{\alpha \nu}^{0}(1,4) \varSigma^{\nu \mu}(4,3) {\rm{i}} G_{\mu \beta}(3,2)}, \end{split} $$ (9) 其中$ \varSigma^{\nu \mu}(4,3) $是$ \Delta $的自能,在波恩近似下,有:
$$ \varSigma^{\nu \mu}(4,3) = \varSigma_{{\rm{HF}}}^{\nu \mu}(4,3)+\varSigma_{{\rm{Born}}}^{\nu \mu}(4,3); $$ (10) 这里$ \varSigma_{{\rm{HF}}}^{\nu \mu}(4,3) $包括Hartree项和Fock项,
$$ \varSigma_{{\rm{HF}}}^{\nu \mu}(4,3) = \varSigma_{{\rm{H}}}^{\nu \mu}(4,3)+\varSigma_{{\rm{F}}}^{\nu \mu}(4,3)。 $$ (11) 引入Rarita-Schwinger场算符,
$$ \begin{split} {\Lambda_{\alpha \beta}\left(\partial \right)} =\;&{\left({\rm{i}} \gamma\boldsymbol\cdot \partial -M_{\Delta}\right)g_{\alpha\beta}-{\rm{i}}\left( \gamma_{\alpha}\partial_{\beta}+\gamma_{\beta}\partial_{\alpha} \right)}+\\ &{{\rm{i}}\gamma_{\alpha}\gamma\boldsymbol\cdot \partial\gamma_{\beta}+M_{\varDelta}\gamma_{\alpha}\gamma_{\beta}}, \end{split} $$ (12) 将其作用在式(9)两端,引入辅助公式:
$$ \Lambda^{\lambda \alpha}\left(\partial_{1}\right)G_{\alpha \beta}^{0}(1,2) = \delta(1,2)g^{\lambda}_{\beta}, $$ (13) 化简可得:
$$ \Lambda^{\lambda \alpha}\left(\partial_{1}\right){\rm{i}}G_{\alpha \beta}(1,2) = {\rm{i}}\delta(1,2)g^{\lambda}_{\beta}+\oint {\rm{d}}x_{3}\varSigma^{\lambda\mu}(1,3){\rm{i}}G_{\mu\beta}(3,2)。 $$ (14) 对于$ G^{-+}_{\alpha\beta}(1,2) $,零阶$ \Delta $格林函数${\rm{i}}G^{0-+}_{\alpha\beta}\left(X,P\right)$通常表示为
$$ \begin{split} { {\rm{i}}G^{0-+}_{\alpha \beta}(x,P)} =\; &{P^{3/2}_{\alpha\beta}(P)\left[ -\left({\cancel{P}}+M_{\varDelta}\right) \frac{\pi}{E_{\varDelta}(P)}\right.}\times \\ &{ \delta\Big[P_{0}-E_{\varDelta}(P)\Big]f_{\varDelta}(X,P)\Biggr]}, \end{split} $$ (15) 这里,定义自旋投影算符$ P^{3/2}_{\alpha\beta}(P) $:
$$ P^{3/2}_{\alpha\beta}(P) = g_{\alpha\beta}-\frac{1}{3}\gamma_{\alpha}\gamma_{\beta}-\frac{1}{3P^{2}}\left({\cancel{P}}\gamma_{\alpha}P_{\beta}+P_{\alpha}\gamma_{\beta}{\cancel{P}} \right), $$ (16) 将式(14)整体做Winger变换,并将式(15)和式(16)带入,在半经典和准粒子近似下,可以得到$ \Delta $分布函数的RBUU方程:
$$ \begin{split} & \Bigg\{p_{\mu}\Big[\partial_{x}^{\mu}-\partial_{x}^{\mu} \varSigma_{\varDelta}^{\nu}(x) \partial_{\nu}^{p}+\partial_{x}^{\nu} \varSigma_{\varDelta}^{\mu}(x) \partial_{\nu}^{p}\Big]+m_{\varDelta}^{*} \partial_{x}^{\nu} \varSigma_{\varDelta}^{S}(x) \partial_{\nu}^{p}\Bigg\}\times \\ & \qquad \dfrac{f_{\varDelta}({\boldsymbol{x}}, {\boldsymbol{p}}, \tau)}{E_{\varDelta}^{*}(p)} = C^{\varDelta}(x, p), \end{split} $$ (17) $$ C^{\varDelta}(x, p) = \frac{1}{4} \int \frac{{\rm{d}} {\boldsymbol{p}}_{2}}{(2 \pi)^{3}} \sigma^{\varDelta}(s, t) v_{\varDelta}\left[F_{2}-F_{1}\right] {\rm{d}} \varOmega。 $$ (18) $ \sigma^{\Delta}(s,t) $ 为$ {\rm{N}} \Delta $弹性散射截面,$ v_{\Delta} $为Mϕller速度,$F_{2}$和$F_{1}$是泡利阻塞因子。$ {\rm{N}} \Delta \rightarrow {\rm{N}} \Delta $散射的跃迁概率表述为
$$ \begin{array}{*{20}{ll}} W^{\varDelta}\left(p, p_{2}, p_{3}, p_{4}\right) = G\left(p, p_{2}, p_{3}, p_{4}\right)+p_{3} \leftrightarrow p_{4} \\ \end{array}, $$ (19) 这里:
$$ G = \frac{g_{\Delta \Delta}^{A} g_{\Delta \Delta}^{B} g_{{\rm{N}} {\rm{N}}}^{A} g_{{\rm{N}} {\rm{N}}}^{B} T_{e} \varPhi_{e}}{16 E_{\varDelta}^{*}(p) E^{*}\left(p_{2}\right) E_{\varDelta}^{*}\left(p_{3}\right) E^{*}\left(p_{4}\right)}。 $$ (20) $ T_{e} $ 是同位旋矩阵, $ \varPhi_{e} $ 是自旋矩阵,它们是描述跃迁概率中最重要的两个输入量。散射截面和跃迁概率的关系表示为
$$ \begin{split} {\int v \frac{{\rm{d}}\sigma^{*}}{{\rm{d}}\varOmega } {\rm{d}} \varOmega} = &{\int \frac{{\rm{d}}^{3} p_{3}}{(2 \pi)^{3}} \int \frac{{\rm{d}}^{3} p_{4}}{(2 \pi)^{3}}(2 \pi)^{4} \delta^{4}(p+p_{2}-p_{3}-p_{4})}\times \\ & { W^{\varDelta}\left(p, p_{2}, p_{3}, p_{4}\right)}。\\[-11pt] \end{split} $$ (21) 考虑到核子的有限尺寸和短程关联效应,在计算反应截面中引入如下形状因子:
$$ F_{{\rm{NN}}i}(t) = \frac{\Lambda_{{\rm{NN}}i}^{2}}{\Lambda_{{\rm{NN}}i}^{2}-t}, $$ (22) 对于$ {\rm{N}} \Delta \rightarrow {\rm{N}} \Delta $反应道,核子截断质量取$\varLambda_{{\rm{NN}}\sigma}$= 1 100 ${\rm{MeV}}$、$\varLambda_{{\rm{NN}} \omega}$= 783 ${\rm{MeV}}$、$\varLambda_{{\rm{NN}} \delta}$= 983 ${\rm{MeV}}$、$\varLambda_{{\rm{NN}} \rho}$= 770 ${\rm{MeV}}$[17, 27],取${\varLambda_{\Delta \Delta i}}/{\varLambda_{{\rm{NN}} i }} = 0.4$,$ i $为$\sigma$、$ \omega $、$\delta$、$\rho$。
需要指出的是,本文重点关注的是密度、能量依赖的$ {\rm{N}} \Delta $弹性散射截面的同位旋劈裂效应,并且计算的结果将通过进一步的参数化引入至微观输运模型中,而在微观输运模型中,共振态粒子的衰变宽度已被细致考虑。此外,我们之前的研究表明$ \Delta $的衰变宽度的密度依赖性和质量劈裂效应都较弱[30, 41]。因此,在本文中将重子看作点粒子,并且忽略$ \Delta $的衰变宽度可能带来的影响。
-
从第1节的公式推导中可以看出,RBUU 输运理论中平均场和碰撞项分别来自自能展开的二阶项和四阶项,而且在数值计算中,通过平均场计算得到的重子有效质量将是介质修正的散射截面的输入量,可以认为RBUU 输运理论能够自洽地处理平均场和两体碰撞。其中$ \Delta $有效质量、核子质量和介子场的平均值有如下关系[42]:
$$ \begin{array}{*{20}{ll}} m_{{\rm{p}} / {\rm{n}}}^{*} = m_{{\rm{N}}}-g_{\sigma} \sigma \mp g_{\delta} \delta_{0},\\ m_{\varDelta^{++} / \varDelta^-}^{*} = m_{\varDelta}-g_{\sigma} \sigma \mp g_{\delta} \delta_{0}, \\ m_{\varDelta^+ / \varDelta^{0}}^{*} = m_{\varDelta}-g_{\sigma} \sigma \mp \frac{1}{3} g_{\delta} \delta_{0} 。 \end{array} $$ (23) 图1给出了同位旋不对称度 $\alpha = 0.6[\alpha = (\rho_{{\rm{n}}} -\rho_{{\rm{p}}})/ (\rho_{{\rm{n}}} +\rho_{{\rm{p}}})]$时,$ \Delta $有效质量随约化重子密度的变化关系,同时以同位旋不对称度 $ \alpha = 0 $时核子和$ \Delta $的有效质量随约化重子密度的变化作为参考,这里使用的饱和密度处核子的有效质量$m_{{\rm{N}}}^{\ast}/m_{{\rm{N}}}$为0.55。可以看到,当$ \alpha = 0 $时,重子的不同态并没有质量劈裂现象的出现(核子:红色实线;$ \Delta $:黑色实线)。当在有效拉氏量中引入 $ \delta $ 介子场之后,核子与$ \Delta $粒子的不同同位旋态的有效质量发生了明显的劈裂现象,并且$m_{\Delta^{++}}^{*} > m_{\varDelta^{+}}^{*}> m_{\varDelta^{0}}^{*}> m_{\Delta^{-}}^{*}$,$m_{{\rm{p}}}^{*} > m_{{\rm{n}}}^{*}$。相应地,$ \Delta $粒子的有效质量劈裂将会影响与不同同位旋态$ \Delta $粒子产生相关的反应道的散射截面,具体结果将在下文中展示和讨论。
图2展示了同位旋不对称度 $ \alpha = 0.6 $时,不同重子密度($ \rho = 0.5\rho_{0} $、$ \rho_{0} $、$ 2\rho_{0} $)下,同位旋无关(iso-indep.)和同位旋依赖(iso-dep.)的${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$的散射截面随质心系能量的变化。这里的同位旋无关指截面计算中仅考虑了同位旋标量介子$\sigma$和$\omega$的贡献,而同位旋依赖的截面包含$\sigma$,$ \omega $,$\delta$,$\rho$四种介子。这里密度依赖的耦合常数有效范围是$0.5 \sim 2\rho_{0}$。首先,从图中可以明显地看出,在低能区,散射截面随着质心系能量的增大迅速地降低;在高能区,散射截面随着质心系能量的增大而逐渐趋于平缓。其次,还可以看到相同颜色的两线条(iso-indep. 和iso-dep.)之间的差异随着密度的增大而逐渐减小并趋于靠近,表明同位旋矢量介子场对${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$的总散射截面的贡献随着密度的增大而逐渐减小。因为随着重子密度的增加,部分分反应道的散射截面以及总散射截面逐渐降低(如图3所示),因此,同位旋矢量介子场对它们的影响也逐渐变得不明显。需要指出的是,在较高能量下饱和密度处的散射截面与Dirac-Brueckner(DB)的计算结果[43]大致相同。
图 2 同位旋不对称度 $ \alpha = 0.6 $时,不同重子密度下同位旋无关和同位旋依赖的${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$的散射截面随质心系能量的变化关系(在线彩图)
图3展示了质心系能量为2.58 GeV时,${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$过程不同反应道的弹性散射截面随约化重子密度的变化情况。图中实线(虚线)分别表示质子(中子)与$ \Delta $不同同位旋态的弹性散射截面。在Born项中考虑不同介子交换时,${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$的各个分反应道对应的同位旋矩阵$ T_{e} $如表1所列。首先,可以看到不同反应道的弹性散射截面均存在着较强或较弱的密度依赖性,这是由于Born项中存在$ \sigma-\delta $、$ \sigma-\rho $、$\omega-\delta$、$\omega-\rho$交换,并且它们对截面有着较大的贡献。而在总的弹性散射截面中,这些交换项的作用相互抵消,因此,在总的弹性散射截面的计算中只存在纯的同位旋标量介子交换项和同位旋矢量介子交换项。其次,由于$\rm p\Delta^{++}$和$\rm n\Delta^{-}$的弹散过程有着相同的同位旋矩阵(如表1所列),同时,下面同组散射截面也有着与$ p\Delta^{++} $和$ n\Delta^{-} $相同的同位旋矩阵规律:$\rm n\Delta^{++}(\rm p\Delta^{-})$、$\rm p\Delta^{+}(\rm n\Delta^{0})$、$\rm n\Delta^{+}(\rm p\Delta^{0})$,因此,同一组弹散截面的跃迁概率只和自旋矩阵相关,且重子有效质量也是自旋矩阵的重要输入量之一,因此$\rm p\Delta^{++}$和$(\rm n\Delta^{-})$出现了不同的密度依赖关系,其他三组截面也有此类规律。此外,质子(实线)、中子(虚线)与$ \Delta $不同同位旋态的弹性散射截面也分别趋于靠近,当密度为0.5$ \rho_{0} $时,$\sigma_{{\rm{p}}\Delta^{++}}^{*}/\sigma_{{\rm{p}}\Delta^{+}}^{*}$= 2.88,$\sigma_{{\rm{n}}\Delta^{-}}^{*}/\sigma_{{\rm{n}}\Delta^{0}}^{*}$= 1.62,而在 1.5$ \rho_{0} $时,$\sigma_{{\rm{p}}\Delta^{++}}^{*}/ \sigma_{{\rm{p}}\Delta^{+}}^{*}$= 1.94,$\sigma_{{\rm{n}}\Delta^{-}}^{*}/\sigma_{{\rm{n}}\Delta^{0}}^{*}$= 1.25,表明$ \Delta $同位旋效应对分反应道的影响随密度增大逐渐减弱,并且可以推测在$ 2\rho_{0} $之后的高密区,$ {\rm{p}}\Delta^{*} $和$ {\rm{n}}\Delta^{*} $分别会可能趋于各自的固定数值。
图 3 同位旋不对称度 $ \alpha = 0.6 $,质心系能量为2.58 GeV时,${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$各分反应道弹性散射截面随重子数密度的变化情况(在线彩图)
表 1 $ {\rm{N}} \Delta \rightarrow {\rm{N}} \Delta $ 过程中考虑不同介子交换时各分反应道对应的同位旋矩阵$ T_{e} $
$\sigma-\sigma, \sigma-\omega\$
$\omega-\omega$$\delta-\delta, \delta-\rho,$
$\rho-\rho$$\sigma-\delta, \sigma-\rho,$
$\omega-\delta, \omega-\rho$${\rm{p}}\Delta^{++}({\rm{n}}\Delta^{-})$ 1 9/4 3/2 ${\rm{n}}\Delta^{++}({\rm{p}}\Delta^{-})$ 1 9/4 −3/2 ${\rm{p}}\Delta^{+}({\rm{n}}\Delta^{0})$ 1 1/4 1/2 ${\rm{n}}\Delta^{+}({\rm{p}}\Delta^{0})$ 1 1/4 −1/2
The Calculation of the In-medium Isospin-dependent NΔ → NΔ Cross Section Based on the Self-consistent RBUU Theory
-
摘要: 基于自洽的RBUU输运理论,研究了介质中同位旋相关的${\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$ 散射截面。发现同位旋效应对重子的有效质量以及不同重子密度区的总$ {\rm{N}} \Delta $弹性散射截面有着较为明显的影响。随着重子密度的增大,不同重子同位旋态之间的有效质量劈裂逐渐增大。在密度依赖的重子有效质量劈裂、耦合常数以及Born项中的$\sigma-\delta$、$\sigma-\rho$、$\omega-\delta$、$\omega-\rho$交换项的共同作用下,不同同位旋态分反应道的弹性散射截面呈现出了不同的密度依赖行为。总截面在低能区出现了明显的介质压低效应,高能区介质效应减弱。Abstract: Based on the self-consistent RBUU transport theory, the isospin-dependent in-medium $ {\rm{N}} \Delta \rightarrow {\rm{N}} \Delta$ cross section is investigated. It is found that the isospin effect has a relatively obvious influence on the effective mass of baryons and the total $ {\rm{N}} \Delta $ elastic cross section in different density regions. With the increase of baryon density, the effective mass splitting between different isospin states of baryons increases gradually. Under the joint effect of density-dependent baryon effective mass splitting, coupling constant, as well as Born terms such as $\sigma-\delta$, $\sigma-\rho$, $\omega-\delta$, $\omega-\rho$, the elastic cross-sections of sub-channels with different isospin states exhibit different density-dependent behaviors. The total cross-section has an obvious reduction effect of medium in the low-energy region, and the medium effect is weakened in the high-energy region.
-
Key words:
- Δ resonance /
- RBUU theory /
- isospin dependence /
- in-medium effect
-
表 1 $ {\rm{N}} \Delta \rightarrow {\rm{N}} \Delta $ 过程中考虑不同介子交换时各分反应道对应的同位旋矩阵$ T_{e} $
$\sigma-\sigma, \sigma-\omega\$
$\omega-\omega$$\delta-\delta, \delta-\rho,$
$\rho-\rho$$\sigma-\delta, \sigma-\rho,$
$\omega-\delta, \omega-\rho$${\rm{p}}\Delta^{++}({\rm{n}}\Delta^{-})$ 1 9/4 3/2 ${\rm{n}}\Delta^{++}({\rm{p}}\Delta^{-})$ 1 9/4 −3/2 ${\rm{p}}\Delta^{+}({\rm{n}}\Delta^{0})$ 1 1/4 1/2 ${\rm{n}}\Delta^{+}({\rm{p}}\Delta^{0})$ 1 1/4 −1/2 -
[1] HANAUSKE M, STEINHEIMER J, MOTORNENKO A, et al. Particles, 2019, 2(1): 44. doi: 10.3390/particles2010004 [2] OERTEL M, HEMPEL M, KLÄHN T, et al. Rev Mod Phys, 2017, 89(1): 015007. doi: 10.1103/RevModPhys.89.015007 [3] 肖志刚. 物理, 2020, 49: 137. doi: 10.7693/wl20200301 XIAO Z G. Physics, 2020, 49: 137. (in Chinese) doi: 10.7693/wl20200301 [4] YONG G, GUO Y. Nucl Phys Rev, 2020, 37(2): 136. doi: 10.11804/NuclPhysRev.37.2019068 [5] 徐骏. 中国科学: 物理学力学天文学, 2020, 49: 5. doi: 10.1360/SSPMA-2019-0068 XU J. Sci Sin Phys Mech Astro, 2020, 49: 5. (in Chinese) doi: 10.1360/SSPMA-2019-0068 [6] LIU Y, YONG G, ZUO W. Europhysics Letters, 2012, 99(4): 42001. doi: 10.1209/0295-5075/99/42001 [7] GUO W, YONG G, WANG Y, et al. Phys Lett B, 2013, 726(1-3): 211. doi: 10.1016/j.physletb.2013.07.056 [8] LI P, WANG Y, STEINHEIMER J, et al. Phys Lett B, 2021, 818: 136393. doi: 10.1016/j.physletb.2021.136393 [9] 刘恒金, 冯兆庆. 中国科学: 物理学力学天文学, 2023, 53: 50. doi: 10.1360/SSPMA-2022-0473 LIU H J, FENG Z Q. Sci Sin Phys Mech Astro, 2023, 53: 50. (in Chinese) doi: 10.1360/SSPMA-2022-0473 [10] ZOU L, LI M, GUO C, et al. Sci China Phys Mech Astron, 2016, 59(12): 122011. doi: 10.1007/s11433-016-0358-y [11] BOHNET A, OHTSUKA N, AICHELIN J, et al. Nucl Phys A, 1989, 494(2): 349. doi: 10.1016/0375-9474(89)90028-6 [12] HAN S, SHANG X, ZUO W, et al. Phys Rev C, 2022, 106(6): 064332. doi: 10.1103/PhysRevC.106.064332 [13] LI G Q, MACHLEIDT R. Phys Rev C, 1994, 49: 566. doi: 10.1103/PhysRevC.49.566 [14] PANDHARIPANDE V, PIEPER S C. Phys Rev C, 1992, 45(2): 791. doi: 10.1103/PhysRevC.45.791 [15] CUI Y, ZHANG Y, LI Z. Chin Phys C, 2020, 44(2): 024106. doi: 10.1088/1674-1137/44/2/024106 [16] HUBER S, AICHELIN J. Nucl Phys A, 1994, 573: 587. doi: 10.1016/0375-9474(94)90232-1 [17] LI Q, LI Z, ZHAO E. Phys Rev C, 2004, 69(1): 017601. doi: 10.1103/PhysRevC.69.017601 [18] LI Q, LI Z, MAO G. Phys Rev C, 2000, 62(1): 014606. doi: 10.1103/PhysRevC.62.014606 [19] CAI X Z, FENG J, SHEN W Q, et al. Phys Rev C, 1998, 58: 572. doi: 10.1103/PhysRevC.58.572 [20] LI P, WANG Y, LI Q, et al. Phys Lett B, 2022, 828: 137019. doi: 10.1016/j.physletb.2022.137019 [21] 李庆峰, 李祝霞, BLEICHER M, 等. 原子核物理评论, 2011, 28(2): 142. doi: 10.11804/NuclPhysRev.28.02.142 LI Q F, LI Z X, BLEICHER M, et al. Nucl Phys Rev, 2011, 28(2): 142. (in Chinese) doi: 10.11804/NuclPhysRev.28.02.142 [22] ZHANG Y, LI Z, DANIELEWICZ P. Phys Rev C, 2007, 75(3): 034615. doi: 10.1103/PhysRevC.75.034615 [23] LATTIMER J, PRAKASH M. Astrophys J, 2001, 550(1): 426. doi: 10.1086/319702 [24] HOFMANN M, MATTIELLO R, AMELIN N, et al. Nucl Phys A, 1994, 566: 15c. doi: 10.1016/0375-9474(94)90605-X [25] XIAO Z, YONG G, CHEN L, et al. Eur Phys J A, 2014, 50(2): 37. doi: 10.1140/epja/i2014-14037-6 [26] FERINI G, GAITANOS T, COLONNA M, et al. Phys Rev Lett, 2006, 97(20): 202301. doi: 10.1103/PhysRevLett.97.202301 [27] LI Q, ZHAO E. Mod Phys Lett A, 2003, 18(38): 2713. doi: 10.1142/S021773230301226X [28] LI Q, LI Z, SOFF S, et al. J Phys G, 2006, 32: 407. doi: 10.1088/0954-3899/32/4/001 [29] LI Q, LI Z. Phys Lett B, 2017, 773: 557. doi: 10.1016/j.physletb.2017.09.013 [30] 李庆峰, 李祝霞. 原子核物理评论, 2018, 35(4): 374. doi: 10.11804/NuclPhysRev.35.04.374 LI Q F, LI Z X. Nucl Phys Rev, 2018, 35(4): 374. (in Chinese) doi: 10.11804/NuclPhysRev.35.04.374 [31] LIU Y, WANG Y, CUI Y, et al. Phys Rev C, 2021, 103(1): 014616. doi: 10.1103/PhysRevC.103.014616 [32] LANG A, CASSING W, MOSEL U, et al. Nucl Phys A, 1992, 541(3): 507. doi: 10.1016/0375-9474(92)90189-Q [33] CUGNON J, MIZUTANI T, VANDERMEULEN J. Nucl Phys A, 1981, 352(3): 505. doi: 10.1016/0375-9474(81)90427-9 [34] MAO G, LI Z, ZHUO Y. Phys Rev C, 1996, 53(6): 2933. doi: 10.1103/PhysRevC.53.2933 [35] TONDEUR F, BRACK M, FARINE M, et al. Nucl Phys A, 1984, 420(2): 297. doi: 10.1016/0375-9474(84)90444-5 [36] ROBLEDO L, RODRÍGUEZ T, RODRÍGUEZ-GUZMÁN R. J Phys G, 2018, 46(1): 013001. doi: 10.1088/1361-6471/aadebd [37] LIU Z, XIA C, LU W, et al. Phys Rev C, 2018, 98(2): 024316. doi: 10.1103/PhysRevC.98.024316 [38] MACHLEIDT R, SAMMARRUCA F. Physica Scripta, 2016, 91(8): 083007. doi: 10.1088/0031-8949/91/8/083007 [39] HOFMANN F, KEIL C, LENSKE H. Phys Rev C, 2001, 64(3): 034314. doi: 10.1103/PhysRevC.64.034314 [40] RADUTA A R. Phys Lett B, 2021, 814: 136070. doi: 10.1016/j.physletb.2021.136070 [41] LI Q, LI Z. Sci China Phys Mech Astron, 2019, 62: 1. doi: 10.1007/s11425-017-9341-1 [42] SONG T, KO C M. Phys Rev C, 2015, 91(1): 014901. doi: 10.1103/PhysRevC.91.014901 [43] TER HAAR B, MALFLIET R. Phys Rev C, 1987, 36(4): 1611. doi: 10.1103/PhysRevC.36.1611