Mbinu ya awamu inayobadilika (variable-phase method) ni mbinu inayobadilisha mlinganyo wa radial Schrödinger wa daraja la pili katika mtawanyiko usio wa kirelativisti wa miili miwili kuwa mfumo wa milinganyo tofauti ya daraja la kwanza unaofuatilia viwianishi vya vipengele huru vya Riccati–Bessel na Riccati–Neumann vya suluhisho. Utafiti huu unaendeleza mbinu hiyo ili kukokotoa suluhisho za radial regular na irregular ndani ya mfumo mmoja wa milinganyo, na hasa unazalisha data nyingi za kihesabu kwa mtawanyiko wa masafa marefu wa mada nyeusi katika potenshali ya Yukawa. Suluhisho regular hutoa mabadiliko ya awamu na vipengele vya Sommerfeld, ilhali suluhisho irregular huruhusu kutolewa kwa taarifa ya mwisho ya Green function ya wimbi la \(s\), inayohitajika katika unitarization inayojitosheleza ya michakato isiyo elastiki inayofanana na contact.

Msukumo wa kifizikia wa utafiti ni kwamba nadharia ya perturbation inaweza kutotosha katika hali ambapo masafa ya mwingiliano yanaweza kulinganishwa na au kuwa makubwa kuliko urefu wa wimbi wa de Broglie wa chembe. Hasa katika fenomenolojia ya mada nyeusi, nguvu za masafa marefu zinaweza kuathiri wingi wa mabaki ya joto, viwango vya annihilation kwa ugunduzi usio wa moja kwa moja, na uundaji wa miundo ya galaksi katika matukio ya mada nyeusi inayojingiliana yenyewe. Katika rejimu kama hiyo, mwingiliano unaojirudia wa chembe mbili lazima uzingatiwe katika daraja zote, na kwa hiyo tatizo la Schrödinger hutatuliwa kwa njia isiyo ya perturbative.

Waandishi wanatumia mbinu hiyo kwenye potenshali ya Yukawa na kuchunguza kwa pamoja mwingiliano wa kuvutana na kusukumana; rejimu za Born, resonance, classical/semi-classical na zinazofanana na Coulomb; mabadiliko ya awamu, vipengele vya Sommerfeld, cross section za elastic/transfer/viscosity na convergence ya partial-wave. Matokeo ya kihesabu yanalinganishwa, inapofaa, na mipaka ya kiuchanganuzi ya Born, WKB na Coulomb. Aidha, kutoka kwa suluhisho irregular la \(\ell=0\), kiasi \(T(k;k_\diamond)\) kinakokotolewa.

Jedwali za data zilizochapishwa zinajumuisha mabadiliko ya awamu na vipengele vya Sommerfeld kutoka \(\ell=0\) hadi \(\ell=100\). Jumla za partial-wave zinazotumiwa kwa cross section za elastic, transfer na viscosity zinaendelezwa hadi \(\ell=100\), \(99\) na \(98\), mtawalia. Kulingana na uthibitishaji wa utafiti wenyewe, katika nafasi ya vigezo iliyowekwa kwenye jedwali, uchaguzi huu huweka makosa ya kukata chini ya %1. Hata hivyo, utafiti haupimi mwingiliano katika jaribio halisi la mada nyeusi wala hauonyeshi kwamba modeli fulani ya chembe ya mada nyeusi ni sahihi kiasili; matokeo yake ni suluhisho za kinadharia na kihesabu katika kiwango cha Schrödinger kwa potenshali fulani za mtawanyiko.

Kwa nini mtawanyiko wa masafa marefu ni tatizo tofauti?

Ikiwa masafa ya mwingiliano wa chembe mbili ni madogo sana kuliko urefu wao wa wimbi wa de Broglie, mara nyingi mwingiliano unaweza kutazamwa kama marekebisho ya masafa mafupi. Lakini masafa ya nguvu yanapolingana na urefu wa wimbi, chembe huhisiana katika eneo kubwa la anga. Mwingiliano unaojirudia wa wave function katika eneo hilo unaweza kuzalisha athari kubwa ambazo upanuzi wa perturbative wa daraja la chini unaweza kukosa.

Sehemu ya kuanzia ya utafiti ni mlinganyo wa Schrödinger usiotegemea muda kwa chembe mbili zisizo za kirelativisti:

\[ -\frac{\nabla^2}{2\mu}\psi_{\mathbf{k}}(\mathbf r) +V(r)\psi_{\mathbf{k}}(\mathbf r) = E_{\mathbf{k}}\psi_{\mathbf{k}}(\mathbf r). \tag{2.1} \]

Hapa \(\mu\) ni reduced mass ya chembe mbili, \(V(r)\) ni potenshali halisi, ya eneo maalum na ya kati, \(\mathbf{k}\) ni momentum katika mfumo wa kituo cha molekuli, na

\[ E_{\mathbf{k}}=\frac{k^2}{2\mu}, \qquad k=\mu v_{\rm rel} \]

mahusiano haya hutumika. \(v_{\rm rel}\) ni kasi ya jamaa ya chembe mbili. Katika natural units, vitengo vya momentum, mass na inverse length vinahusiana.

Kwa Nini Upanuzi wa Partial-Wave Hutumika?

Upanuzi wa partial-wave hugawa tatizo la mtawanyiko wa pande tatu katika potenshali ya kati kuwa matatizo huru ya radial ya mwelekeo mmoja kwa namba tofauti za quantum za angular momentum \(\ell\). Kwa njia hii, wave function kamili inaweza kutatuliwa kama jumla ya vipengele ambavyo kila kimoja kinawakilisha njia fulani ya angular momentum, na cross section ya jumla ya mtawanyiko hupatikana kutoka kwa jumla inayodhibitiwa ya njia hizi.

Wave function kamili:

\[ \psi_{\mathbf{k}}(\mathbf r) = \sum_{\ell=0}^{\infty} (2\ell+1) P_\ell(\cos\vartheta)\, \psi_{k,\ell}(r) \tag{2.4} \]

hugawanywa kwa namna hii. \(P_\ell\) ni polynomia ya Legendre ya \(\ell\)-th. Kwa radial function, tukifafanua

\[ u_{k,\ell}(r) \equiv kr\,\psi_{k,\ell}(r) \tag{2.5} \]

tatizo hugeuka kuwa:

\[ \left[ -\frac{1}{2\mu}\frac{d^2}{dr^2} + \frac{\ell(\ell+1)}{2\mu r^2} + V(r) \right] u_{k,\ell}(r) = E_{\mathbf{k}}u_{k,\ell}(r) \tag{2.6} \]

mlinganyo wa radial Schrödinger.

\[ \frac{\ell(\ell+1)}{2\mu r^2} \]

ni centrifugal barrier. Kwa \(\ell=0\) term hii haipo; kadiri \(\ell\) inavyoongezeka, hasa kwa \(r\) ndogo, inakuwa kizuizi chenye nguvu zaidi. Kwa hiyo tabia za resonance za partial-wave za chini na za juu si sawa.

Darasa la potenshali katika utafiti

Mbinu ya jumla hutumika kwa potenshali zinazotimiza sharti:

\[ \int_0^\infty dr\,r\,|V(r)|<\infty. \tag{2.7} \]

Sharti hili huhakikisha kwamba potenshali haina singularity yenye nguvu kuliko \(1/r^2\) kwenye origin na kwamba inapungua kwa kasi ya kutosha katika umbali mkubwa. Chini ya dhana hii, mlinganyo wa radial una familia mbili zinazojitegemea kwa mstari:

\[ \mathcal F_{k,\ell}(r)\propto r^{\ell+1}, \qquad \mathcal G_{k,\ell}(r)\propto r^{-\ell}, \qquad r\rightarrow0. \tag{2.8} \]

\(\mathcal F_{k,\ell}\) ni familia ya suluhisho regular, na \(\mathcal G_{k,\ell}\) ni familia ya suluhisho irregular.

Regular na irregular zinamaanisha nini?

Hali ya kimwili ya mtawanyiko lazima iwe finite kwenye origin. Kwa kuwa full radial wave function ni \(u/(kr)\), sehemu irregular inayodiverge kwenye origin haitumiki katika hali ya kimwili ya wimbi linaloingia. Hata hivyo, suluhisho irregular si lisilo na umuhimu wa kihisabati; hubeba taarifa inayohitajika kukokotoa Green function na muundo uliorenormalize/kuunitarize wa contact interactions.

Mabadiliko ya awamu

Katika potenshali halisi, uwezekano huhifadhiwa kwa kila partial-wave katika mtawanyiko elastic. Amplitude ya mtawanyiko wa partial-wave:

\[ 1+ \frac{2ik}{\sqrt{\mathscr C_\ell}} f_{k,\ell} = e^{2i\theta_\ell(k)} \tag{2.12} \]

inahusishwa na mabadiliko ya awamu \(\theta_\ell(k)\) yanayofafanuliwa na uhusiano huu. Hapa \(\mathscr C_\ell\) ni symmetry factor inayotumika kwa chembe zinazoweza kutofautishwa na zinazofanana.

Kutokana na hili:

\[ f_{k,\ell} = \sqrt{\mathscr C_\ell}\, \frac{1}{k} e^{i\theta_\ell(k)} \sin\theta_\ell(k) \tag{2.13} \]

hupatikana. Mabadiliko ya awamu huweka kwa kiasi mabadiliko ya phase yanayosababishwa na mwingiliano ikilinganishwa na free wave.

Cross section elastic ya kila partial-wave:

\[ \sigma_\ell^{\rm elas}(k) = \sigma_\ell^{\rm U}(k) \sin^2\theta_\ell(k), \tag{2.16} \]

na kikomo cha partial-wave unitarity:

\[ \sigma_\ell^{\rm U}(k) = \mathscr C_\ell \frac{4\pi}{k^2}(2\ell+1) \tag{2.17} \]

hufafanuliwa. Wakati \(\sin^2\theta_\ell=1\), partial-wave hiyo hufikia kikomo chake cha unitarity.

Faktori ya Sommerfeld Hupima Nini?

Faktori ya Sommerfeld \(S_\ell(k)\) hupima ni kwa kiasi gani potenshali ya masafa marefu hubadilisha amplitude ya wave function katika umbali mfupi ikilinganishwa na hali huru; hivyo huamua marekebisho ya kiwango cha multiplicative ambayo mwingiliano wa masafa marefu huleta kwenye viwango vya michakato isiyo elastiki ya umbali mfupi inayofanana na contact, kama annihilation.

Msingi wa ufafanuzi:

\[ S_\ell(k) = \lim_{r\rightarrow0} \left| \frac{ \mathcal F_{k,\ell}(r) }{ \mathcal F_{k,\ell}^{\rm free}(r) } \right|^2. \tag{2.20} \]

\(S_\ell>1\) ina maana amplitude ya wave function katika umbali mfupi imeongezeka ikilinganishwa na hali huru, na \(S_\ell<1\) ina maana imebanwa. Utafiti huu haufafanui faktori ya Sommerfeld kama nguvu mpya ya kimwili; unakokotoa marekebisho ya wave function ambayo mwingiliano wa masafa marefu huleta kwenye mchakato uliopo wa umbali mfupi.

Suluhisho irregular na \(T(k;k_\diamond)\)

Kwa wimbi la \(s\), ikiwa potenshali karibu na origin ni:

\[ V(r) = -\frac{\alpha_{-1}}{r} + O(r^0) \tag{2.22} \]

basi suluhisho irregular linaweza kuandikwa kwa upanuzi wa Frobenius kama:

\[ \mathcal G_{k,0}(r) = \mathcal G_{k,0}(0) \left[ 1+ g_0^{(1)}(k)r - 2\mu\alpha_{-1}r\ln(\mu r) + O(r^2,r^2\ln r) \right] \tag{2.23} \]

.

Katika contact-like \(s\)-wave interactions, projection ya Green function katika \(r=r'=0\) hudiverge kwenye ultraviolet. Lakini katika hali inayoweza kurenormalize inayochunguzwa na utafiti, sehemu ya divergence isiyotegemea momentum inaweza kuondolewa kwa kutumia reference momentum \(k_\diamond\). Kiasi finite kisicho na dimension:

\[ T(k;k_\diamond) = \frac{1}{k} \left[ g_0^{(1)}(k) - g_0^{(1)}(k_\diamond) \right] \tag{2.24} \]

hufafanuliwa. \(T(k;k_\diamond)\) ni mojawapo ya bidhaa za data zinazochapishwa na makala na hutumika katika unitarization inayojitosheleza ya michakato isiyo elastiki inayofanana na contact.

Mbinu ya Awamu Inayobadilika Hubadilishaje Mlinganyo wa Schrödinger?

Njia ya awamu inayobadilika huandika radial wave function kama mchanganyiko wa suluhisho mbili huru za Riccati–Bessel/Riccati–Neumann za tatizo huru zenye coefficients zinazotegemea nafasi; hivyo mlinganyo wa Schrödinger wa daraja la pili hubadilika kuwa milinganyo miwili ya tofauti ya daraja la kwanza iliyounganishwa kwa coefficients hizi mbili.

Kigezo kisicho na dimension:

\[ x\equiv kr \]

huchaguliwa na mlinganyo huandikwa:

\[ \chi_{k,\ell}''(x) + b_\ell(x)\chi_{k,\ell}(x) = \widetilde V(x)\chi_{k,\ell}(x) \tag{3.1} \]

ambapo:

\[ \widetilde V(x) = \frac{2\mu}{k^2}V(x/k), \qquad b_\ell(x) = 1-\frac{\ell(\ell+1)}{x^2}. \tag{3.2} \]

Suluhisho huru ni:

\[ \hat j_\ell(x)=xj_\ell(x), \qquad \hat n_\ell(x)=xn_\ell(x) \tag{3.3} \]

ambazo ni functions za Riccati–Bessel na Riccati–Neumann. Uhusiano wa Wronskian:

\[ \hat j_\ell(x)\hat n_\ell'(x) - \hat j_\ell'(x)\hat n_\ell(x) = 1 \tag{3.4} \]

unamaanisha functions hizi mbili zinajitegemea kwa mstari.

Suluhisho jipya huchaguliwa kama:

\[ \chi_{k,\ell}(x) = c_{k,\ell}(x)\hat j_\ell(x) - s_{k,\ell}(x)\hat n_\ell(x) \tag{3.6a} \]

na derivative yake:

\[ \chi_{k,\ell}'(x) = c_{k,\ell}(x)\hat j_\ell'(x) - s_{k,\ell}(x)\hat n_\ell'(x) \tag{3.6b} \]

. Matokeo yake ni mfumo mkuu wa ukokotoaji wa utafiti:

\[ c_{k,\ell}'(x) = -\widetilde V(x)\hat n_\ell(x) \left[ c_{k,\ell}(x)\hat j_\ell(x) - s_{k,\ell}(x)\hat n_\ell(x) \right], \tag{3.8a} \]

\[ s_{k,\ell}'(x) = -\widetilde V(x)\hat j_\ell(x) \left[ c_{k,\ell}(x)\hat j_\ell(x) - s_{k,\ell}(x)\hat n_\ell(x) \right]. \tag{3.8b} \]

Mabadiliko haya ni muhimu: badala ya kutatua moja kwa moja mlinganyo wa radial wa daraja la pili ambao unaweza kuoscillate haraka, evolution ya daraja la kwanza ya coefficients mbili hufuatiliwa.

Masharti ya mpaka ya suluhisho regular

Kwa familia regular, kwenye origin:

\[ \lim_{x\rightarrow0} c_{k,\ell}^{\rm rgl}(x)=1, \qquad \lim_{x\rightarrow0} s_{k,\ell}^{\rm rgl}(x)=0 \tag{3.10} \]

hutumiwa. Kwa kuwa integration ya kihesabu haiwezi kuanzishwa kwenye \(x=0\), utafiti hutumia kwenye \(x_{\min}\) ndogo lakini finite upanuzi wa kiuchanganuzi wa small-\(x\) kama masharti ya kuanzia.

Mabadiliko ya awamu yanatolewaje?

\[ t_{k,\ell}(x) = \frac{ s_{k,\ell}^{\rm rgl}(x) }{ c_{k,\ell}^{\rm rgl}(x) } \tag{3.11} \]

chini ya ufafanuzi huu:

\[ \tan\theta_\ell(k) = \lim_{x\rightarrow\infty} \frac{ s_{k,\ell}^{\rm rgl}(x) }{ c_{k,\ell}^{\rm rgl}(x) }. \tag{3.14} \]

Chanzo pia hutolea mlinganyo wa kawaida unaofuatilia moja kwa moja awamu inayobadilika:

\[ \bar\theta_{k,\ell}'(x) = -\widetilde V(x) \left[ \cos\bar\theta_{k,\ell}(x)\hat j_\ell(x) - \sin\bar\theta_{k,\ell}(x)\hat n_\ell(x) \right]^2. \tag{3.17} \]

na:

\[ \theta_\ell(k) = \lim_{x\rightarrow\infty} \bar\theta_{k,\ell}(x). \tag{3.18} \]

Hata hivyo, waandishi katika uzalishaji wa data hawatumii moja kwa moja Mlinganyo (3.17), bali mfumo wa Mlinganyo (3.8) kwa \(c\) na \(s\); kwa sababu wanaripoti kwamba kwa vitendo ni thabiti zaidi kihesabu.

Muundo wa variable-phase wa faktori ya Sommerfeld

Coefficients za large-\(x\) za suluhisho regular zinapotumika:

\[ S_\ell(k) = \left[ \left(c_{k,\ell;\infty}^{\rm rgl}\right)^2 + \left(s_{k,\ell;\infty}^{\rm rgl}\right)^2 \right]^{-1}. \tag{3.23} \]

Formula hii ni mojawapo ya faida muhimu za ukokotoaji za mbinu: integration hiyo hiyo ya \(c,s\) hutoa mabadiliko ya awamu na faktori ya Sommerfeld.

Integration ya suluhisho irregular kwa mwelekeo wa kinyume

Sharti la mpaka la irregular katika umbali mkubwa hujengwa kutoka kwa coefficients za asymptotic za suluhisho regular:

\[ c_{k,\ell;\infty}^{\rm irr} = s_{k,\ell;\infty}^{\rm rgl}, \qquad s_{k,\ell;\infty}^{\rm irr} = -c_{k,\ell;\infty}^{\rm rgl}. \tag{3.24} \]

Katika utekelezaji wa kihesabu, suluhisho regular huintegratiwa kwanza kutoka \(x_{\min}\) hadi \(x_{\max}\). Baada ya coefficients za asymptotic kubainishwa, suluhisho irregular huintegratiwa tena kutoka \(x_{\max}\) kuelekea \(x\) ndogo.

Kwa wimbi la \(s\), \(T(k;k_\diamond)\) inaweza kuandikwa moja kwa moja kama:

\[ T(k;k_\diamond) = \lim_{r\rightarrow0} \left[ \frac{ \chi_{k,0}^{{\rm irr}\,'}(kr) }{ \chi_{k,0}^{\rm irr}(kr) } - \frac{k_\diamond}{k} \frac{ \chi_{k_\diamond,0}^{{\rm irr}\,'}(k_\diamond r) }{ \chi_{k_\diamond,0}^{\rm irr}(k_\diamond r) } \right] \tag{3.30} \]

. Terms hizi mbili zinapotathminiwa kwenye thamani ileile ya kimwili ya \(r\), divergence isiyotegemea momentum hufutiana. Waandishi wanapendelea katika utekelezaji wa kihesabu kuondoa sehemu inayodiverge kwa uchanganuzi na kukokotoa coefficient ya mstari iliyobaki kando.

Potenshali ya Yukawa

Mfano wa matumizi:

\[ V_Y(r) = -\frac{\alpha}{r} e^{-m_{\rm med}r}. \tag{4.1} \]

Hapa \(\alpha\) ni nguvu ya coupling na inaweza kuwa chanya au hasi; \(m_{\rm med}\) ni mass ya mediator anayebeba mwingiliano. \(m_{\rm med}\) finite huscreen mwingiliano kwa exponential katika umbali mkubwa ikilinganishwa na potenshali ya Coulomb.

Vigezo visivyo na dimension \(\zeta\) na \(\xi\)

Utafiti hupanga mtawanyiko wa Yukawa kwa vigezo viwili visivyo na dimension:

\[ \zeta \equiv \frac{\mu\alpha}{k} = \frac{\alpha}{v_{\rm rel}}, \qquad \xi \equiv \frac{\mu|\alpha|}{m_{\rm med}}. \tag{4.2} \]

\(\zeta\) hupima nguvu ya coupling dhidi ya relative velocity. Absolute value yake inaweza kuwa kubwa katika kasi ndogo. \(\xi\) ni kigezo cha range/nguvu kinachohusiana na scale za mass na mediator za mwingiliano. Potenshali ya Yukawa iliyofanywa dimensionless:

\[ \widetilde V_Y(x) = -\frac{2\zeta}{x} e^{-x|\zeta|/\xi}. \tag{4.3} \]

huchukua umbo hili.

Kwa nini integration ya kihesabu haiwezi kufanywa kwenye interval ileile?

Katika nafasi ya vigezo, scale bainifu kati ya potenshali, centrifugal barrier na nishati ya chembe zinaweza kubadilika kwa orders nyingi. Kwa hiyo kutumia values zilezile za \(x_{\min}\) na \(x_{\max}\) katika kila parameter point kunaweza katika baadhi ya maeneo kutoa ama sharti la mwanzo lisilo sahihi au integration inayomalizika kabla ya kufikia eneo la asymptotic.

Utafiti hufafanua scale zifuatazo:

\[ x_1(\zeta,\xi)=\frac{\xi}{|\zeta|}, \tag{4.5b} \]

\[ x_2(\ell,\zeta) = \frac{\ell+1}{|\zeta|}, \tag{4.5c} \]

\[ x_3(\zeta,\xi) = \frac{\xi}{|\zeta|} W_0\left( \frac{2\zeta^2}{\xi} \right), \tag{4.5d} \]

na kwa \(\ell\geq1\), inayojumuisha makutano ya kwanza ya terms za centrifugal na potenshali:

\[ x_4(\ell,\zeta,\xi) = \begin{cases} \dfrac{\xi}{|\zeta|} \left[ -W_0\left( -\dfrac{\ell(\ell+1)}{2\xi} \right) \right], & \xi\geq\dfrac{e}{2}\ell(\ell+1), \\[8pt] +\infty, & \xi<\dfrac{e}{2}\ell(\ell+1). \end{cases} \tag{4.5e} \]

Pia balance ya nishati–centrifugal:

\[ x_5(\ell) = \sqrt{\ell(\ell+1)} \tag{4.5f} \]

hutambulishwa. \(W_0\) ni tawi kuu la Lambert \(W\) function.

Kikomo cha chini cha integration:

\[ x_{\min}(\ell=0,\zeta,\xi) = \varepsilon \min\{1,x_1,x_2,x_3\}, \tag{4.6a} \]

\[ x_{\min}(\ell\geq1,\zeta,\xi) = \varepsilon \min\{1,x_1,x_2,x_3,x_4\} \tag{4.6b} \]

hufafanuliwa. Code inayoambatana na chanzo hutumia:

\[ \varepsilon=10^{-7} \]

.

Kikomo cha juu:

\[ x_{\max}(\ell,\zeta,\xi) = \Lambda \max\{1,x_1,x_3,x_5\}, \tag{4.7} \]

na kwenye code huchaguliwa:

\[ \Lambda=10^4 \]

. Kwa hiyo eneo la integration si fixed, bali huscale kulingana na kila parameter point ya kimwili.

Scale ya takriban ya partial-wave inayobainisha rejimu za mtawanyiko

Scale ya classical angular momentum inayolingana na range ya nguvu ni:

\[ \ell_{\rm clas} \sim r_{\rm force}k \sim \frac{\mu v_{\rm rel}}{m_{\rm med}} = \frac{\xi}{|\zeta|}. \tag{4.8} \]

Kiasi hiki kinapokuwa kikubwa, partial-wave nyingi huchangia mtawanyiko.

Resonance katika potenshali ya Yukawa ya kuvutana

Kwa wimbi la \(s\), threshold ya kwanza ya bound-state hutolewa takriban kama:

\[ \xi_{1,0}\simeq0.84 \]

. Bound state inapokaribia nishati sifuri, \(\theta_0\) hukaribia takriban \(\pi/2\) modulo \(\pi\), na cross section elastic ya wimbi la \(s\) inaweza kukaribia kikomo chake cha unitarity.

Kwa \(\ell\) kubwa, muundo ni tofauti kutokana na centrifugal barrier. Hali quasi-bound zenye nishati chanya chini ya threshold zinaweza kunaswa kwa muda nyuma ya barrier na kuleta resonance nyembamba. Chanzo kinaeleza wazi kwamba hii si nakala rahisi ya tabia ya wimbi la \(s\) katika partial-wave za juu.

Kati ya thresholds za bound-state kunaweza kutokea Ramsauer–Townsend anti-resonances na kubana sana mtawanyiko elastic wa nishati ndogo. Potenshali ya Yukawa ya kusukumana haiungi mkono resonance za bound-state.

Tabia ya faktori ya Sommerfeld kwenye resonance

Katika potenshali ya Yukawa ya kuvutana, kwenye thresholds za bound-state chanzo kinaripoti scaling:

\[ S_0\propto\zeta^2 \]

. Mbali na resonance, kwa sababu ya range finite ya Yukawa, enhancement ya Sommerfeld inaweza kusaturate kwa \(\zeta\) kubwa. Kwa upande mwingine, katika strict Coulomb limit:

\[ S_0^{\rm C}\propto\zeta \qquad (\zeta\gg1) \]

ndiyo tabia.

Kutoa coefficient finite kutoka suluhisho irregular

Katika potenshali ya Yukawa, coefficient husika ya suluhisho irregular kwa \(x\) ndogo hukokotolewa baada ya kuondoa analytic logarithmic singularity. Utafiti hutumia:

\[ \widetilde g_0^{(1)}(k) = g_0^{(1)}(k) + 2\alpha\mu\ln|\alpha| \tag{4.14} \]

na matokeo hutathminiwa kama:

\[ T(k;k_\diamond) = \frac{\widetilde g_0^{(1)}(k)}{k} - \frac{k_\diamond}{k} \frac{\widetilde g_0^{(1)}(k_\diamond)}{k_\diamond} \tag{4.15} \]

. Shift isiyotegemea momentum katika ufafanuzi hutoweka katika tofauti kati ya momenta mbili.

Kwa nini cross section za elastic, transfer na viscosity zinatofautiana?

Cross section ya jumla ya elastic hupima tu uwezekano wa jumla wa mtawanyiko, ilhali uhamishaji wa momentum na nishati katika halo za mada nyeusi hutegemea pia pembe ya mtawanyiko. Kwa hiyo utafiti hutumia cross section tatu tofauti:

\[ \sigma_{\rm elas} = \int d\Omega\, \frac{d\sigma}{d\Omega}, \tag{5.1a} \]

\[ \sigma_T = \int d\Omega\, (1-\cos\vartheta) \frac{d\sigma}{d\Omega}, \tag{5.1b} \]

\[ \sigma_V = \int d\Omega\, \sin^2\vartheta \frac{d\sigma}{d\Omega}. \tag{5.1c} \]

\(\sigma_T\) hupunguza uzito wa forward scattering na kulenga uhamishaji halisi wa momentum. \(\sigma_V\) hutumia uzito wa \(\sin^2\vartheta\), unaopunguza scattering ya moja kwa moja mbele na nyuma.

Kwa chembe zinazoweza kutofautishwa, miundo ya partial-wave ni:

\[ \sigma_\ell^{\rm elas} = \frac{4\pi}{k^2} (2\ell+1)\sin^2\theta_\ell, \tag{5.3a} \]

\[ \sigma_\ell^T = \frac{4\pi}{k^2} (\ell+1) \sin^2(\theta_{\ell+1}-\theta_\ell), \tag{5.3b} \]

\[ \sigma_\ell^V = \frac{4\pi}{k^2} \frac{(\ell+1)(\ell+2)}{2\ell+3} \sin^2(\theta_{\ell+2}-\theta_\ell). \tag{5.3c} \]

Kwa kuwa cross section za transfer na viscosity zina tofauti za awamu za partial-wave zilizo jirani, mkia wa high-\(\ell\) unaweza kuonyesha suppression ya ziada ikilinganishwa na cross section ya elastic. Sifa hii huwa muhimu baadaye katika uchanganuzi wa truncation.

Makadirio ya Born

Makadirio ya kwanza ya Born hubadilisha hali kamili ya mtawanyiko kwa free incoming wave. Kwa potenshali yenye spherical symmetry:

\[ f_k^{\rm Born}(\vartheta) = -\frac{2\mu}{q} \int_0^\infty dr\,r\,V(r)\sin(qr), \qquad q=2k\sin(\vartheta/2). \tag{5.4} \]

Kwa potenshali ya Yukawa, differential cross section ni:

\[ \frac{d\sigma^{\rm Born}}{d\Omega} = \frac{ 4\mu^2\alpha^2 }{ \left[ m_{\rm med}^2+ 2k^2(1-\cos\vartheta) \right]^2 }. \tag{5.5} \]

Mabadiliko ya awamu ya Born:

\[ \theta_\ell^{\rm Born,Y}(k) = \zeta Q_\ell \left[ 1+ \frac{1}{2} \left( \frac{\zeta}{\xi} \right)^2 \right], \tag{5.8} \]

ambapo \(Q_\ell\) ni Legendre function ya aina ya pili.

Waandishi hufafanua:

\[ \cosh\eta = 1+ \frac{1}{2} \left( \frac{\zeta}{\xi} \right)^2, \qquad \sinh(\eta/2) = \frac{|\zeta|}{2\xi} \tag{5.9} \]

na kufanya suppression ya exponential ya mkia wa high-\(\ell\) ionekane wazi:

\[ \theta_{\ell\rightarrow\infty}^{\rm Born,Y} \approx \zeta e^{-(\ell+1/2)\eta} \sqrt{ \frac{\pi}{ 2(\ell+1/2)\sinh\eta } }. \tag{5.10} \]

Usemi huu ni muhimu sana kuelewa kwa nini jumla ya partial-wave huconverge na kwa nini convergence hupungua kasi wakati \(\eta\) inapungua.

Makadirio ya WKB

Katika weak-potential WKB limit, mabadiliko ya awamu ya Yukawa hupatikana kama:

\[ \theta_\ell^{\rm WKB,Y}(k) \approx \zeta K_0 \left[ \left(\ell+\frac{1}{2}\right) \frac{|\zeta|}{\xi} \right], \tag{5.11} \]

. \(K_0\) ni modified Bessel function ya aina ya pili. Makadirio haya hutumika hasa kuelewa scaling ya eneo la semi-classical ambako partial-wave nyingi huchangia.

Born na WKB hazibadilishi nini?

Makala haitumii fomula za Born au WKB kama matokeo exact katika nafasi yote ya vigezo. Suluhisho kamili la kihesabu la variable-phase ndilo ukokotoaji mkuu. Born hutumika kama chombo cha kiuchanganuzi cha kulinganisha kwa weak-interaction regime na mkia wa high-\(\ell\), huku WKB ikitumika kuelewa tabia ya semi-classical na kukagua matokeo ya kihesabu.

Mbinu na Matokeo ya Utafiti

Mchakato wa ukokotoaji

1. Tatizo la Schrödinger la miili miwili hugawanywa katika partial-wave.
2. Tatizo hufanywa dimensionless kwa \(x=kr\).
3. Suluhisho radial huandikwa kwa coefficients zinazotegemea nafasi za functions za Riccati–Bessel na Riccati–Neumann.
4. Milinganyo miwili ya tofauti ya daraja la kwanza katika Mlinganyo (3.8) hutatuliwa.
5. Suluhisho regular huintegratiwa kutoka \(x_{\min}\) hadi \(x_{\max}\), na mabadiliko ya awamu pamoja na faktori ya Sommerfeld hutolewa.
6. Coefficients zilezile za asymptotic huunda masharti ya mpaka ya suluhisho irregular kwenye \(x_{\max}\).
7. Suluhisho irregular huintegratiwa kwa mwelekeo wa kinyume hadi \(x\) ndogo.
8. Singularity inayojulikana ya origin huondolewa kwa uchanganuzi na \(T(k;k_\diamond)\) hukokotolewa.
9. Cross section za elastic, transfer na viscosity hujengwa kutoka mabadiliko ya awamu.
10. Kosa la truncation la jumla ya partial-wave hukaguliwa kwa majaribio ya kiuchanganuzi na ya moja kwa moja ya kihesabu.

Bidhaa za data zilizochapishwa katika chanzo

Tatizo la partial-wave truncation

Jumla isiyo na mwisho lazima ikatwe kihesabu kwenye \(\ell_{\max}\) finite:

\[ \sigma_\star^{(\ell_{\max})} = \sum_{\ell=0}^{\ell_{\max}} \sigma_{\ell,\star}, \qquad \star\in\{\mathrm{elas},T,V\}. \tag{5.14} \]

Mkia halisi unaokosekana:

\[ \Delta\sigma_\star^{(\ell_{\max})} = \sigma_\star - \sigma_\star^{(\ell_{\max})} = \sum_{\ell=\ell_{\max}+1}^{\infty} \sigma_{\ell,\star}. \tag{5.15} \]

Kosa la jamaa la incremental linalopima mchango mmoja wa partial-wave unaofuata:

\[ \varepsilon_{\ell_{\max}}^\star = \frac{ \left| \sigma_\star^{(\ell_{\max}+1)} - \sigma_\star^{(\ell_{\max})} \right| }{ \sigma_\star^{(\ell_{\max})} }. \tag{5.16} \]

Hata hivyo, chanzo kinatoa tahadhari muhimu: kwa kuwa kiasi hiki hupima term moja tu inayofuata, kinaweza kudharau kosa la mkia wote uliobaki.

Convergence ya Jumla ya Partial-Wave Hukaguliwaje?

Utafiti hutumia kupungua kwa mabadiliko ya awamu ya Born ya high-\(\ell\) katika umbo la \(\exp[-(\ell+1/2)\eta]\) kukadiria mkia uliobaki wa partial-wave kama mfululizo wa takriban wa kijiometria; hivyo badala ya kuangalia term ya mwisho iliyoongezwa pekee, ukubwa wa jamaa wa mkia wote unaokosekana unaweza kukokotolewa na kutathmini kwa kiasi ikiwa \(\ell_{\max}\) inatosha.

Makadirio ya kiuchanganuzi ya mkia kwa cross section elastic:

\[ \frac{ \Delta\sigma_{\rm elas}^{(\ell_{\max})} }{ \sigma_{\ell=0}^{\rm U} } \approx \frac{ 2\pi\xi^4 }{ 4\xi^2+\zeta^2 } e^{-2(\ell_{\max}+1)\eta}. \tag{5.18} \]

Kwa cross section transfer:

\[ \frac{ \Delta\sigma_T^{(\ell_{\max})} }{ \sigma_{\ell=0}^{\rm U} } \approx \frac{\pi\zeta^2}{2\sinh\eta} \tanh(\eta/2) e^{-(2\ell_{\max}+3)\eta}. \tag{5.21} \]

Kwa cross section viscosity:

\[ \frac{ \Delta\sigma_V^{(\ell_{\max})} }{ \sigma_{\ell=0}^{\rm U} } \approx \frac{\pi\zeta^2}{2} e^{-2(\ell_{\max}+2)\eta}. \tag{5.24} \]

Katika fomula zote tatu, kosa hupungua kwa exponential kadiri \(\ell_{\max}\) inavyoongezeka. Kigezo cha kudhibiti ni \(\eta\). Ikiwa:

\[ \frac{|\zeta|}{\xi}\ll1, \]

basi:

\[ \eta\simeq\frac{|\zeta|}{\xi} \]

hupungua na suppression ya exponential hudhoofika. Kwa hiyo kadiri tunavyoingia katika eneo la semi-classical/linalofanana na Coulomb, partial-wave zaidi zinahitajika.

Katika approximation ya mkia wa kijiometria, uhusiano kati ya kosa kamili la jamaa na term inayofuata ni:

\[ \frac{ \Delta\sigma_\star^{(\ell_{\max})} }{ \sigma_\star^{(\ell_{\max})} } \approx \frac{ \varepsilon_{\ell_{\max}}^\star }{ 1-e^{-2\eta} }, \qquad \ell_{\max}\gtrsim\ell_{\rm elas}^{\rm WKB}. \tag{5.25} \]

Chanzo kinaeleza kwamba kwa \(\eta\simeq0.1\):

\[ (1-e^{-2\eta})^{-1}\simeq5.5 \]

. Kwa hiyo kuangalia tu swali “term ya mwisho ni ndogo kiasi gani?” kunaweza kuwa na matumaini kupita kiasi kwa kosa la jumla.

Jaribio gumu zaidi la convergence

Chanzo kinajaribu kando kama point ngumu zaidi ya convergence katika eneo lililowekwa kwenye jedwali:

\[ (\zeta,\xi) = (10^4,10^5) \]

. Kwa ukaguzi wa moja kwa moja, mabadiliko ya awamu ya kihesabu yanaendelezwa hadi \(\ell=164\), na mkia wa Born hadi \(\ell=1000\).

Thamani zote tatu ziko chini ya %1, yaani chini ya kizingiti cha \(10^{-2}\). Matokeo haya hayamaanishi kwamba “partial-wave 100 zinatosha kila wakati kwa tatizo lolote la Yukawa”; ni udhibiti wa truncation uliothibitishwa tu katika eneo la vigezo lililochapishwa na utafiti.

Uthibitishaji kwa kikomo cha Coulomb

Kwa potenshali kali ya Coulomb:

\[ V_C(r) = -\frac{\alpha}{r} \]

kuna suluhisho za kiuchanganuzi. Mabadiliko ya awamu ya Coulomb yanayotumiwa na utafiti ni:

\[ \theta_\ell^C(k) = \arg\Gamma(1+\ell-i\zeta), \tag{B.2a} \]

faktori ya Sommerfeld:

\[ S_\ell^C(k) = \frac{2\pi\zeta}{1-e^{-2\pi\zeta}} \prod_{j=1}^{\ell} \left( 1+\frac{\zeta^2}{j^2} \right) \tag{B.2b} \]

na taarifa finite ya Green function:

\[ T^C(k;k_\diamond) = 2\zeta \left[ \ln\left( \frac{\zeta}{\zeta_\diamond} \right) + {\rm Re}\,H_{i\zeta_\diamond} - {\rm Re}\,H_{i\zeta} \right] \tag{B.2c} \]

hutolewa.

Kuna jambo muhimu hapa: mabadiliko ya awamu ya Yukawa yenye mass finite ya mediator hayakaribii moja kwa moja mabadiliko ya awamu ya Coulomb point-by-point katika basis ya partial-wave wakati \(\xi\rightarrow\infty\). Awamu ya pamoja inayokua logarithmically hukusanyika. Kwa upande mwingine, vipengele vya Sommerfeld na \(T(k;k_\diamond)\) vina mipaka ya Coulomb iliyofafanuliwa vizuri. Kwa hiyo kauli kwamba “potenshali ya Yukawa inakaribia Coulomb” haimaanishi kila kiasi kinachokokotolewa kinaconverge point-by-point kwa njia ileile.

Majukumu ya Born, WKB na suluhisho la kihesabu

Matokeo yanayoungwa mkono na utafiti

• Suluhisho za Schrödinger regular na irregular zinaweza kupatikana kihesabu kwa mfumo uleule wa milinganyo ya daraja la kwanza wa variable-phase.
• Coefficients za asymptotic za suluhisho regular hutoa mabadiliko ya awamu na faktori ya Sommerfeld kutoka ukokotoaji uleule.
• Suluhisho irregular la wimbi la \(s\) linaweza kutoa taarifa ya \(T(k;k_\diamond)\) inayotumika katika unitarization ya michakato isiyo elastiki inayofanana na contact.
• Kwa potenshali ya Yukawa, uchaguzi dynamic wa \(x_{\min}\) na \(x_{\max}\) unaotokana na scale za kimwili hutoa utaratibu unaotegemeka zaidi kuliko mipaka fixed ya integration katika scans pana za vigezo.
• Mipaka ya Born, WKB na Coulomb hutoa ukaguzi wa consistency wa suluhisho la kihesabu katika maeneo tofauti ya vigezo.
• Kosa la partial-wave truncation linaweza kukadiriwa kiuchanganuzi kwa kutumia mkia wa Born wa high-\(\ell\).
• Katika eneo la vigezo lililochapishwa na utafiti, truncations za \(\ell=100,99,98\) huweka makosa ya truncation ya cross section za elastic, transfer na viscosity chini ya %1.
• Code inaweza kubadilishwa kwa potenshali mbalimbali halisi, za eneo maalum na za kati; lakini extraction ya \(T(k;k_\diamond)\) pia inahitaji kujua tabia inayohitajika ya umbali mfupi ya potenshali.

Ujumlishaji ambao utafiti hauungi mkono

• Utafiti hauthibitishi kwa majaribio uwepo wa mada nyeusi au modeli fulani ya mada nyeusi.
• Haidaiwi kwamba potenshali ya Yukawa inawakilisha mwingiliano wote wa mada nyeusi.
• Haijaonyeshwa kwamba values zilizochapishwa za \(\ell_{\max}\) zinatosha katika nafasi zote zinazowezekana za vigezo.
• Makadirio ya Born na WKB hayawezi kutumika badala ya suluhisho exact la kihesabu katika vigezo vyote.
• Katika potenshali ya Yukawa ya kusukumana, resonance za bound-state za potenshali ya kuvutana hazitarajiwi na chanzo hakiripoti matokeo kama hayo.
• Bound-state formation na bremsstrahlung overlap integrals ambazo hazijawekwa kwenye jedwali hazitolewi kama matokeo tayari katika data set; chanzo kinasema full wave functions zinazohitajika zinaweza kuzalishwa kwa code iliyochapishwa.

Uzalishaji upya

Utafiti hauishii tu kuchapisha milinganyo. Code ya Mathematica iliyotumiwa imechapishwa wazi; matokeo ya jedwali ya potenshali ya Yukawa yamehifadhiwa katika data set tofauti ya Zenodo. Chanzo kinasema repo ina user guide inayoeleza parameter grid, data format na interpolation functions.

Maelezo ya Chanzo na Mbinu

Utafiti asilia: Variable-phase methods and data for long-range scattering

Waandishi: Marcos M. Flores, Kalliopi Petraki, Thomas Pignard.

Mpangilio wa waandishi: Marcos M. Flores → Kalliopi Petraki → Thomas Pignard.

Taasisi: Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, Paris, France. Marcos M. Flores pia ana afiliasheni ya Department of Physics, University of Oslo, Norway.

Mawasiliano ya waandishi: Ukurasa wa kichwa wa chanzo una anwani za barua pepe za waandishi wote watatu. Kwa kuwa hati haimtaji mtu mmoja wazi kama “corresponding author”, Verianla haijamteua mwandishi mmoja wa mawasiliano.

Mchango sawa: Hakuna tamko la mchango sawa lililotambuliwa katika chanzo.

Jukwaa: arXiv.

arXiv: 2610.30261v1 [hep-ph].

arXiv DOI: 10.48550/arXiv.2610.30261. DOI hii ni kitambulisho cha kudumu cha rekodi ya arXiv; haipaswi kutafsiriwa kama DOI ya jarida lililopitiwa na wahakiki.

Tarehe ya kuwasilisha arXiv: 21 Oktoba 2026.

Tarehe iliyoandikwa kwenye PDF: 23 Oktoba 2026.

Maelezo ya tarehe: Tarehe 23 Oktoba 2026 kwenye ukurasa wa kichwa wa PDF si sawa na tarehe 21 Oktoba 2026 katika rekodi ya kuwasilisha arXiv v1. Taarifa zote mbili zimehifadhiwa katika muktadha wake.

Aina ya chanzo: Preprint ya utafiti; utafiti wa kinadharia na kihesabu wa high energy physics / astrophysical phenomenology.

Hali ya peer review: Kufikia 4 Novemba 2026, rekodi rasmi ya arXiv haionyeshi journal reference ya uchapishaji wa jarida uliopitiwa na wahakiki. Kwa hiyo chanzo hakijawasilishwa kama makala iliyopitiwa na wahakiki.

Kategoria kuu ya arXiv: High Energy Physics – Phenomenology (hep-ph).

Kategoria mtambuka: Cosmology and Nongalactic Astrophysics (astro-ph.CO) na High Energy Astrophysical Phenomena (astro-ph.HE).

Urefu wa utafiti: 48 kurasa, 13 michoro.

Mbinu kuu: Partial-wave decomposition, radial Schrödinger equation, variable-phase first-order system, outward regular solution, inward irregular solution, Sommerfeld factors, phase shifts, Green-function subtraction, Yukawa potential, Born/WKB/Coulomb analytical checks na partial-wave truncation analysis.

Code: Repo ya code ya Mathematica iitwayo YukawaVPS imechapishwa wazi pamoja na utafiti.

Data: Data set “YukawaVPS data: phase shifts, Sommerfeld factors, and cross sections for Yukawa scattering” imehifadhiwa chini ya Zenodo DOI 10.5281/zenodo.23475703.

Ufadhili: Utafiti uliungwa mkono na programu ya utafiti na ubunifu ya Umoja wa Ulaya Horizon 2020, Grant Agreement No. 101002846, ERC CoG CosmoChart. Marcos M. Flores pia anataja msaada wa Norwegian Research Council FRIPRO Grant, Project ID 353561 “DarkTurns”.

Mgongano wa maslahi: Hakuna tamko tofauti la mgongano wa maslahi lililotambuliwa katika PDF ya chanzo; Verianla haijazalisha tamko ambalo halipo katika chanzo.

CRediT / michango ya waandishi: Hakuna jedwali tofauti la michango ya CRediT au tamko la kina la michango ya waandishi lililotambuliwa katika PDF ya chanzo.

Leseni: Kwa kuwa leseni maalum ya Creative Commons kwa matumizi upya ya PDF ya chanzo haikuweza kuthibitishwa kwa uaminifu, open access haijatafsiriwa kama open license. Grafu za chanzo hazijatumika tena moja kwa moja.

Mbinu ya kuona: Mchoro wa mbinu ndani ya Verianla ni SVG huru isiyo na scale iliyotayarishwa kwa msingi wa mahusiano ya kihisabati yaliyoelezwa katika chanzo; si nakala ya muundo wa kuona wa michoro 13 ya utafiti chanzo. Data ghafi kwenye grafu za chanzo haijakadiriwa kutoka PDF na kubadilishwa kuwa data set mpya.

Kikomo kikuu cha kisayansi: Utafiti hutatua tatizo la mtawanyiko wa Schrödinger lisilo la kirelativisti. Matokeo ni hesabu za kinadharia/kihesabu kwa potenshali na vigezo fulani; si uchunguzi wa moja kwa moja wa astrophysics, matokeo ya detector ya chembe, au uthibitisho wa majaribio wa microphysics ya mada nyeusi.

Kikomo cha ujumlishaji: Milinganyo ya variable-phase inaweza kutumika kwa darasa la potenshali halisi, za eneo maalum na za kati katika Mlinganyo (2.7). Extraction ya \(T(k;k_\diamond)\) pia inahitaji tabia inayofaa ya umbali mfupi ya Mlinganyo (2.22).