Değişken-faz yöntemi (variable-phase method), nonrelativistik iki-cisim saçılmasındaki ikinci mertebeden radyal Schrödinger denklemini, çözümün serbest Riccati–Bessel ve Riccati–Neumann bileşenlerinin katsayılarını izleyen birinci mertebeden diferansiyel denklemler sistemine dönüştüren bir yöntemdir. Bu çalışma, söz konusu yöntemi hem regular hem de irregular radyal çözümleri aynı denklem sistemi içinde hesaplayacak biçimde geliştiriyor ve özellikle Yukawa potansiyelindeki uzun menzilli karanlık madde saçılması için kapsamlı nümerik veri üretiyor. Regular çözüm faz kaymalarını ve Sommerfeld faktörlerini verirken, irregular çözüm temas-benzeri inelastik süreçlerin öz-tutarlı unitarizasyonunda gerekli olan sonlu \(s\)-dalga Green fonksiyonu bilgisinin çıkarılmasını sağlıyor.

Çalışmanın fiziksel motivasyonu, etkileşim menzilinin parçacıkların de Broglie dalga boyuyla karşılaştırılabilir veya ondan büyük olduğu durumlarda pertürbasyon teorisinin yetersiz kalabilmesidir. Özellikle karanlık madde fenomenolojisinde uzun menzilli kuvvetler termal kalıntı bolluğunu, dolaylı tespit için yok olma hızlarını ve kendiliğinden etkileşen karanlık madde senaryolarındaki galaktik yapı oluşumunu etkileyebilir. Böyle bir rejimde tekrarlanan iki-parçacık etkileşimlerinin bütün mertebelerde hesaba katılması gerekir ve bunun için Schrödinger problemi nonpertürbatif olarak çözülür.

Yazarlar yöntemi Yukawa potansiyeline uygulayarak çekici ve itici etkileşimleri; Born, rezonans, klasik/yarı-klasik ve Coulomb-benzeri rejimleri; faz kaymalarını, Sommerfeld faktörlerini, elastik/transfer/viskozite kesitlerini ve kısmi-dalga yakınsamasını birlikte inceliyor. Nümerik sonuçlar uygun olduğu bölgelerde Born, WKB ve Coulomb analitik limitleriyle karşılaştırılıyor. Ayrıca \(\ell=0\) irregular çözümünden elde edilen \(T(k;k_\diamond)\) niceliği hesaplanıyor.

Yayımlanan veri tabloları \(\ell=0\) ile \(\ell=100\) arasındaki faz kaymalarını ve Sommerfeld faktörlerini içeriyor. Elastik, transfer ve viskozite kesitleri için kullanılan kısmi dalga toplamları sırasıyla \(\ell=100\), \(99\) ve \(98\)'e kadar uzatılıyor. Çalışmanın kendi doğrulamasına göre tablolanmış parametre uzayında bu seçimler kesme hatalarını %1'in altında tutuyor. Bununla birlikte çalışma gerçek bir karanlık madde deneyinde etkileşim ölçmemekte ve belirli bir karanlık madde parçacık modelinin doğada doğru olduğunu göstermemektedir; sunduğu sonuçlar belirli saçılma potansiyellerinin Schrödinger düzeyindeki teorik ve nümerik çözümüdür.

Neden uzun menzilli saçılma ayrı bir problem?

İki parçacığın etkileşim menzili onların de Broglie dalga boyundan çok küçükse etkileşim birçok durumda kısa menzilli bir düzeltme gibi ele alınabilir. Ancak kuvvetin menzili dalga boyuyla karşılaştırılabilir olduğunda parçacıklar birbirlerini uzun bir uzaysal bölge boyunca hisseder. Dalga fonksiyonunun bu bölgedeki tekrar tekrar etkileşimi, düşük mertebeli pertürbatif açılımların kaçırabileceği güçlü etkiler yaratabilir.

Çalışmanın başlangıç noktası iki nonrelativistik parçacık için zamandan bağımsız Schrödinger denklemidir:

\[ -\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} \]

Burada \(\mu\) iki parçacığın indirgenmiş kütlesidir, \(V(r)\) gerçek, yerel ve merkezî potansiyeldir, \(\mathbf{k}\) kütle merkezi sistemindeki momentumdur ve

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

ilişkileri kullanılır. \(v_{\rm rel}\), iki parçacığın göreli hızıdır. Doğal birimlerde momentum, kütle ve ters uzunluk birimleri birbirleriyle ilişkilidir.

Kısmi-Dalga Açılımı Neden Kullanılıyor?

Kısmi-dalga açılımı, merkezî bir potansiyelde üç boyutlu saçılma problemini farklı açısal momentum kuantum sayıları \(\ell\) için bağımsız bir-boyutlu radyal problemlere ayırır. Böylece tam dalga fonksiyonu, her biri belirli bir açısal momentum kanalını temsil eden bileşenlerin toplamı olarak çözülebilir ve toplam saçılma kesiti bu kanalların kontrollü toplamından elde edilir.

Tam dalga fonksiyonu:

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

şeklinde ayrılır. \(P_\ell\), \(\ell\)'inci Legendre polinomudur. Radyal fonksiyon için

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

tanımı yapıldığında problem:

\[ \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} \]

biçimindeki radyal Schrödinger denklemine dönüşür.

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

terimi merkezkaç bariyeridir. \(\ell=0\) için bu terim yoktur; \(\ell\) büyüdükçe özellikle küçük \(r\)'de giderek güçlü bir engel oluşturur. Bu nedenle düşük ve yüksek kısmi dalgaların rezonans davranışları aynı değildir.

Çalışmanın potansiyel sınıfı

Genel yöntem şu koşulu sağlayan potansiyellere uygulanır:

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

Bu koşul, potansiyelin orijinde \(1/r^2\)'den daha güçlü tekillik göstermemesini ve büyük uzaklıkta yeterince hızlı azalmasını güvence altına alır. Bu varsayım altında radyal denklemin iki lineer bağımsız ailesi vardır:

\[ \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}\) regular çözüm, \(\mathcal G_{k,\ell}\) ise irregular çözüm ailesidir.

Regular ve irregular ne anlama geliyor?

Fiziksel saçılma durumu orijinde sonlu olmalıdır. Tam radyal dalga fonksiyonu \(u/(kr)\) olduğundan, orijinde diverge eden irregular bileşen fiziksel gelen saçılma durumunda kullanılmaz. Buna karşılık irregular çözüm matematiksel olarak gereksiz değildir; Green fonksiyonunun ve temas-etkileşimlerinin yeniden normalize edilmiş/ünitarize edilmiş yapısının hesaplanmasında gerekli bilgi taşır.

Faz kayması

Gerçek bir potansiyelde elastik saçılmada her kısmi dalga için olasılık korunur. Kısmi-dalga saçılma genliği:

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

ilişkisiyle tanımlanan \(\theta_\ell(k)\) faz kaymasına bağlanır. Burada \(\mathscr C_\ell\), ayırt edilebilir ve özdeş parçacıklarda kullanılan simetri faktörüdür.

Bundan:

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

elde edilir. Faz kayması, etkileşimin serbest dalgaya göre oluşturduğu faz değişimini nicelleştirir.

Her kısmi dalganın elastik kesiti:

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

ve kısmi-dalga unitarite sınırı:

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

olarak tanımlanır. \(\sin^2\theta_\ell=1\) olduğunda o kısmi dalga kendi unitarite sınırına ulaşır.

Sommerfeld Faktörü Neyi Ölçüyor?

Sommerfeld faktörü \(S_\ell(k)\), uzun menzilli potansiyelin kısa mesafedeki dalga fonksiyonu genliğini serbest duruma göre ne kadar değiştirdiğini ölçer; bu nedenle yok olma gibi temas-benzeri kısa-mesafeli inelastik süreçlerin hızlarına uzun menzilli etkileşimin yaptığı çarpan düzeyindeki düzeltmeyi belirler.

Tanımın özü:

\[ 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\) kısa mesafedeki dalga fonksiyonu genliğinin serbest duruma göre büyüdüğünü, \(S_\ell<1\) ise bastırıldığını ifade eder. Bu çalışma Sommerfeld faktörünü yeni bir fiziksel kuvvet olarak tanımlamaz; uzun menzilli etkileşimin mevcut kısa-mesafe sürecine yaptığı dalga-fonksiyonu düzeltmesini hesaplar.

Irregular çözüm ve \(T(k;k_\diamond)\)

\(s\)-dalga için, orijin yakınındaki potansiyel:

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

biçimindeyse irregular çözüm Frobenius açılımıyla:

\[ \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} \]

şeklinde yazılır.

Temas-benzeri \(s\)-dalga etkileşimlerinde Green fonksiyonunun \(r=r'=0\) izdüşümü ultraviolet olarak diverge eder. Ancak çalışmanın ele aldığı yeniden-normalize edilebilir durumda bu diverjansın momentumdan bağımsız bölümü bir referans momentum \(k_\diamond\) kullanılarak çıkarılabilir. Boyutsuz sonlu büyüklük:

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

olarak tanımlanır. \(T(k;k_\diamond)\), makalenin yayımladığı veri ürünlerinden biridir ve temas-benzeri inelastik süreçlerin öz-tutarlı unitarizasyonu için kullanılır.

Değişken-Faz Yöntemi Schrödinger Denklemini Nasıl Dönüştürüyor?

Değişken-faz yaklaşımı, radyal dalga fonksiyonunu serbest problemin iki bağımsız Riccati–Bessel/Riccati–Neumann çözümünün konuma bağlı katsayılarla birleşimi olarak yazar; böylece ikinci-mertebe Schrödinger denklemi, bu iki katsayı için birbiriyle bağlı iki birinci-mertebe diferansiyel denkleme dönüşür.

Boyutsuz değişken:

\[ x\equiv kr \]

seçilir ve denklem:

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

biçimine getirilir. Burada:

\[ \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} \]

Serbest çözümler:

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

Riccati–Bessel ve Riccati–Neumann fonksiyonlarıdır. Wronskian bağıntısı:

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

olduğu için iki fonksiyon lineer bağımsızdır.

Yeni çözüm:

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

ve türevi:

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

olarak seçilir. Bunun sonucunda çalışmanın temel hesaplama sistemi elde edilir:

\[ 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} \]

Bu dönüşüm önemlidir: artık doğrudan hızla salınabilen ikinci-mertebe radyal denklemi çözmek yerine iki katsayının birinci-mertebe evrimi izlenir.

Regular çözümün sınır koşulları

Regular aile için orijinde:

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

alınır. Nümerik entegrasyon \(x=0\)'da başlatılamadığı için çalışma küçük fakat sonlu \(x_{\min}\)'de analitik küçük-\(x\) açılımlarını başlangıç koşulu olarak kullanır.

Faz kayması nasıl çıkarılıyor?

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

tanımı altında:

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

Kaynak ayrıca doğrudan değişken fazı izleyen standart denklemi türetir:

\[ \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} \]

ve:

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

Ancak yazarlar veri üretiminde doğrudan Eşitlik (3.17)'yi değil, \(c\) ve \(s\) için Eşitlik (3.8) sistemini kullanıyor; çünkü pratikte bunun nümerik olarak daha kararlı olduğunu bildiriyorlar.

Sommerfeld faktörünün değişken-faz biçimi

Regular çözümün büyük-\(x\) katsayıları kullanıldığında:

\[ 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} \]

Bu formül yöntemin önemli hesaplama avantajlarından biridir: aynı \(c,s\) entegrasyonu hem faz kaymasını hem Sommerfeld faktörünü verir.

Irregular çözümün ters yönde entegrasyonu

Büyük uzaklıktaki irregular sınır koşulu regular çözümün asymptotik katsayılarından oluşturulur:

\[ 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} \]

Nümerik uygulamada regular çözüm önce \(x_{\min}\)'den \(x_{\max}\)'a doğru entegre edilir. Asimptotik katsayılar belirlendikten sonra irregular çözüm \(x_{\max}\)'dan tekrar küçük \(x\)'e doğru entegre edilir.

\(s\)-dalga için \(T(k;k_\diamond)\) doğrudan:

\[ 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} \]

biçiminde ifade edilebilir. İki terim aynı fiziksel \(r\) değerinde değerlendirildiğinde momentumdan bağımsız diverjans birbirini götürür. Yazarlar nümerik uygulamada diverjan parçayı analitik olarak çıkarıp kalan lineer katsayıyı ayrı ayrı hesaplamayı tercih etmektedir.

Yukawa potansiyeli

Uygulama örneği:

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

Burada \(\alpha\) bağlaşım gücüdür ve pozitif ya da negatif olabilir; \(m_{\rm med}\) etkileşimi taşıyan aracının kütlesidir. Sonlu \(m_{\rm med}\), Coulomb potansiyeline kıyasla etkileşimi büyük uzaklıkta üstel olarak perdeler.

\(\zeta\) ve \(\xi\) boyutsuz parametreleri

Çalışma Yukawa saçılmasını iki boyutsuz parametreyle düzenler:

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

\(\zeta\), bağlaşım gücünü göreli hıza göre ölçer. Mutlak değeri düşük hızlarda büyüyebilir. \(\xi\) ise etkileşimin kütle ve mediator ölçekleriyle ilişkili menzil/güç parametresidir. Boyutsuzlaştırılmış Yukawa potansiyeli:

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

şeklini alır.

Nümerik entegrasyon neden sabit bir aralıkta yapılamıyor?

Parametre uzayında potansiyel, merkezkaç bariyeri ve parçacık enerjisi arasındaki karakteristik ölçekler birçok mertebe değişebilir. Bu nedenle her parametre noktasında aynı \(x_{\min}\) ve \(x_{\max}\) değerlerinin kullanılması bazı bölgelerde ya yanlış başlangıç koşulu ya da asimptotik bölgeye ulaşmadan biten entegrasyon üretir.

Çalışma şu ölçekleri tanımlar:

\[ 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} \]

ve \(\ell\geq1\) için merkezkaç ile potansiyel terimlerinin ilk kesişimini kapsayan:

\[ 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} \]

Ayrıca enerji–merkezkaç dengesi:

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

ile karakterize edilir. \(W_0\), Lambert \(W\) fonksiyonunun ana dalıdır.

Alt entegrasyon sınırı:

\[ 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} \]

olarak belirlenir. Kaynağın beraberindeki kod:

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

kullanır.

Üst sınır:

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

ve kodda:

\[ \Lambda=10^4 \]

olarak seçilir. Böylece entegrasyon alanı sabit değil, her fiziksel parametre noktasına göre ölçeklenir.

Saçılma rejimlerini belirleyen yaklaşık kısmi dalga ölçeği

Kuvvet menziline karşılık gelen klasik açısal momentum ölçeği:

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

Bu nicelik büyük olduğunda çok sayıda kısmi dalga saçılmaya katkıda bulunur.

Çekici Yukawa potansiyelindeki rezonanslar

\(s\)-dalga için ilk bağlı-durum eşiği yaklaşık:

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

olarak verilir. Bağlı durum sıfır enerjiye yaklaştığında \(\theta_0\) yaklaşık \(\pi/2\) modulo \(\pi\)'ye yaklaşır ve \(s\)-dalga elastik kesiti kendi unitarite sınırına yaklaşabilir.

Yüksek \(\ell\)'lerde merkezkaç bariyeri nedeniyle yapı farklıdır. Eşik altındaki pozitif enerjili yarı-bağlı durumlar bariyer arkasında geçici olarak tutulabilir ve dar rezonanslara yol açabilir. Kaynak bunun yüksek kısmi dalgalarda \(s\)-dalga davranışının basit bir kopyası olmadığını özellikle açıklar.

Bağlı-durum eşikleri arasında Ramsauer–Townsend anti-rezonansları oluşabilir ve düşük enerjili elastik saçılmayı güçlü biçimde bastırabilir. İtici Yukawa potansiyeli ise bağlı-durum rezonanslarını desteklemez.

Sommerfeld faktörünün rezonanstaki davranışı

Çekici Yukawa potansiyelinde bağlı-durum eşiklerinde kaynak:

\[ S_0\propto\zeta^2 \]

ölçeklenmesini raporlar. Rezonanstan uzakta sonlu Yukawa menzili nedeniyle büyük \(\zeta\)'da Sommerfeld artışı doygunlaşabilir. Buna karşılık sıkı Coulomb limitinde:

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

davranışı vardır.

Irregular çözümden sonlu katsayının çıkarılması

Yukawa potansiyelinde küçük \(x\)'te irregular çözümün ilgili katsayısı analitik logaritmik tekillik çıkarıldıktan sonra hesaplanır. Çalışma:

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

tanımını kullanır ve sonuç:

\[ 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} \]

şeklinde değerlendirilir. Tanımdaki momentumdan bağımsız kayma iki momentum arasındaki farkta yok olur.

Elastik, transfer ve viskozite kesitleri neden farklı?

Toplam elastik kesit yalnız toplam saçılma olasılığını ölçerken karanlık madde halolarındaki momentum ve enerji taşınımı saçılma açısına da bağlıdır. Bu nedenle çalışma üç ayrı kesit kullanır:

\[ \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\) ileri saçılmayı daha az ağırlıklandırarak net momentum aktarımına odaklanır. \(\sigma_V\) ise hem tam ileri hem tam geri saçılmayı azaltan \(\sin^2\vartheta\) ağırlığını kullanır.

Ayırt edilebilir parçacıklar için kısmi-dalga biçimleri:

\[ \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} \]

Transfer ve viskozite kesitlerinde komşu kısmi dalgaların faz farkları yer aldığı için yüksek-\(\ell\) kuyruğu elastik kesite göre ek bastırma gösterebilir. Bu özellik daha sonra truncation analizinde önemli hale gelir.

Born yaklaşımı

Birinci Born yaklaşımı tam saçılma durumunu serbest gelen dalgayla değiştirir. Küresel simetrik potansiyel için:

\[ 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} \]

Yukawa potansiyeli için diferansiyel kesit:

\[ \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} \]

Born faz kayması:

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

burada \(Q_\ell\), ikinci tür Legendre fonksiyonudur.

Yazarlar:

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

parametresini tanımlayarak büyük-\(\ell\) kuyruğunun üstel bastırılmasını açık hale getirir:

\[ \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} \]

Bu ifade kısmi-dalga toplamının neden yakınsadığını ve \(\eta\) küçüldüğünde yakınsamanın neden yavaşladığını anlamak için merkezi önemdedir.

WKB yaklaşımı

Zayıf-potansiyel WKB limitinde Yukawa faz kayması:

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

olarak elde edilir. \(K_0\), ikinci tür değiştirilmiş Bessel fonksiyonudur. Bu yaklaşım özellikle çok sayıda kısmi dalganın katkıda bulunduğu yarı-klasik bölgenin ölçeklenmesini anlamak için kullanılır.

Born ve WKB neyin yerine geçmiyor?

Makale Born veya WKB formüllerini bütün parametre uzayında exact sonuç olarak kullanmıyor. Tam nümerik variable-phase çözümü temel hesaplamadır. Born zayıf-etkileşim ve yüksek-\(\ell\) kuyruğunu, WKB ise yarı-klasik davranışı anlamak ve nümerik sonuçları kontrol etmek için analitik karşılaştırma araçları olarak kullanılır.

Çalışmanın Yöntemi ve Bulguları

Hesaplama süreci

1. İki-cisim Schrödinger problemi kısmi dalgalara ayrılır.
2. \(x=kr\) ile problem boyutsuzlaştırılır.
3. Radyal çözüm Riccati–Bessel ve Riccati–Neumann fonksiyonlarının konuma bağlı katsayılarıyla ifade edilir.
4. Eşitlik (3.8)'deki iki birinci-mertebe diferansiyel denklem çözülür.
5. Regular çözüm \(x_{\min}\)'den \(x_{\max}\)'a entegre edilerek faz kayması ve Sommerfeld faktörü çıkarılır.
6. Aynı asymptotik katsayılar irregular çözümün \(x_{\max}\)'daki sınır koşullarını oluşturur.
7. Irregular çözüm ters yönde küçük \(x\)'e entegre edilir.
8. Bilinen orijin tekilliği analitik olarak çıkarılarak \(T(k;k_\diamond)\) hesaplanır.
9. Faz kaymalarından elastik, transfer ve viskozite kesitleri oluşturulur.
10. Kısmi-dalga toplamının kesme hatası analitik ve doğrudan nümerik testlerle denetlenir.

Kaynakta yayımlanan veri ürünleri

Kısmi-dalga truncation problemi

Sonsuz toplam nümerik olarak sonlu \(\ell_{\max}\)'ta kesilmek zorundadır:

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

Gerçek eksik kuyruk:

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

Bir sonraki tek kısmi-dalga katkısını ölçen artımlı göreli hata:

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

Ancak kaynak önemli bir uyarı yapar: bu nicelik yalnız bir sonraki terimi ölçtüğünden bütün kalan kuyruğun hatasını küçümseyebilir.

Kısmi-Dalga Toplamının Yakınsadığı Nasıl Kontrol Ediliyor?

Çalışma, yüksek-\(\ell\) Born faz kaymasının \(\exp[-(\ell+1/2)\eta]\) biçimindeki düşüşünden yararlanarak kalan kısmi-dalga kuyruğunu yaklaşık geometrik seri olarak tahmin eder; bu sayede yalnız son eklenen terime bakmak yerine toplam eksik kuyruğun göreli büyüklüğü hesaplanabilir ve \(\ell_{\max}\)'ın yeterli olup olmadığı nicel olarak denetlenebilir.

Elastik kesit için analitik kuyruk tahmini:

\[ \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} \]

Transfer kesiti için:

\[ \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} \]

Viskozite kesiti için:

\[ \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} \]

Üç formülde de \(\ell_{\max}\) büyüdükçe hata üstel olarak azalır. Kontrol parametresi \(\eta\)'dır. Eğer:

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

ise:

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

küçülür ve üstel bastırma zayıflar. Dolayısıyla yarı-klasik/Coulomb-benzeri bölgeye gidildikçe daha fazla kısmi dalga gerekir.

Geometrik kuyruk yaklaşımında tam göreli hata ile bir sonraki terim arasındaki ilişki:

\[ \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} \]

şeklindedir. Kaynak, \(\eta\simeq0.1\) için:

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

olduğunu belirtir. Dolayısıyla yalnız “son terim ne kadar küçük?” sorusuna bakmak toplam hata için fazla iyimser olabilir.

En zor yakınsama testi

Kaynak, tablolanan alanın en zor yakınsama noktası olarak:

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

noktasını ayrıca test eder. Doğrudan kontrol için nümerik faz kaymaları \(\ell=164\)'e, Born kuyruğu ise \(\ell=1000\)'e kadar götürülür.

Üç değer de %1, yani \(10^{-2}\) eşiğinin altındadır. Bu sonuç “100 kısmi dalga her olası Yukawa problemi için daima yeterlidir” anlamına gelmez; yalnız çalışmanın yayımladığı parametre alanındaki doğrulanmış truncation kontrolüdür.

Coulomb limitiyle doğrulama

Sıkı Coulomb potansiyeli:

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

için analitik çözümler mevcuttur. Çalışmanın kullandığı Coulomb faz kayması:

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

Sommerfeld faktörü:

\[ 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} \]

ve sonlu Green-fonksiyonu bilgisi:

\[ 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} \]

olarak verilir.

Burada önemli bir ayrıntı vardır: sonlu mediator kütleli Yukawa faz kayması \(\xi\rightarrow\infty\) olduğunda kısmi-dalga bazında Coulomb faz kaymasına doğrudan noktasal olarak yakınsamaz. Ortak, logaritmik büyüyen bir faz birikir. Buna karşılık Sommerfeld faktörleri ve \(T(k;k_\diamond)\) iyi tanımlı Coulomb limitlerine sahiptir. Dolayısıyla “Yukawa potansiyeli Coulomb'a yaklaşıyor” ifadesi her hesaplanan niceliğin aynı şekilde noktasal yakınsadığı anlamına gelmez.

Born, WKB ve nümerik çözümün rolleri

Çalışmanın desteklediği sonuçlar

• Regular ve irregular Schrödinger çözümleri aynı variable-phase birinci-mertebe denklem sistemiyle nümerik olarak elde edilebilir.
• Regular çözümün asymptotik katsayıları aynı hesaplamadan hem faz kaymasını hem Sommerfeld faktörünü verir.
• Irregular \(s\)-dalga çözümü, temas-benzeri inelastik süreçlerin unitarizasyonunda kullanılan \(T(k;k_\diamond)\) bilgisini sağlayabilir.
• Yukawa potansiyeli için fiziksel ölçeklerden türetilen dinamik \(x_{\min}\) ve \(x_{\max}\) seçimi, geniş parametre taramalarında sabit entegrasyon sınırlarından daha güvenilir bir prosedür sağlar.
• Born, WKB ve Coulomb limitleri nümerik çözümün farklı parametre bölgelerinde tutarlılık kontrollerini sağlar.
• Kısmi-dalga truncation hatası yüksek-\(\ell\) Born kuyruğu kullanılarak analitik olarak tahmin edilebilir.
• Çalışmanın yayımladığı parametre alanında \(\ell=100,99,98\) kesmeleri elastik, transfer ve viskozite kesitlerinin truncation hatalarını %1'in altında tutar.
• Kod farklı gerçek, yerel ve merkezî potansiyellere uyarlanabilir; ancak \(T(k;k_\diamond)\)'nin çıkarılması ayrıca potansiyelin gerekli kısa-mesafe davranışının bilinmesini gerektirir.

Çalışmanın desteklemediği genellemeler

• Çalışma karanlık maddenin varlığını veya belirli bir karanlık madde modelini deneysel olarak doğrulamamaktadır.
• Yukawa potansiyelinin bütün karanlık madde etkileşimlerini temsil ettiği iddia edilmemektedir.
• Yayımlanan \(\ell_{\max}\) değerlerinin tüm olası parametre uzaylarında yeterli olduğu gösterilmemiştir.
• Born ve WKB yaklaşımları bütün parametrelerde exact nümerik çözüm yerine kullanılamaz.
• İtici Yukawa potansiyelinde çekici potansiyeldeki bağlı-durum rezonanslarının bulunması beklenmez ve kaynak böyle bir sonuç raporlamaz.
• Tablolanmayan bound-state formation ve bremsstrahlung overlap integralleri veri setinde hazır sonuç olarak verilmemiştir; kaynak, gerekli tam dalga fonksiyonlarının yayımlanan kodla üretilebileceğini belirtir.

Yeniden üretilebilirlik

Çalışma yalnız denklemleri yayımlamakla kalmıyor. Kullanılan Mathematica kodu açık olarak yayımlanmış; Yukawa potansiyeline ait tablolanmış sonuçlar ayrı bir Zenodo veri setinde arşivlenmiştir. Kaynak, depoda parametre ızgarasını, veri biçimini ve interpolasyon fonksiyonlarını açıklayan kullanıcı kılavuzu bulunduğunu belirtmektedir.

Kaynak ve Yöntem Notu

Özgün çalışma: Variable-phase methods and data for long-range scattering

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

Yazar sırası: Marcos M. Flores → Kalliopi Petraki → Thomas Pignard.

Kurumlar: Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, Paris, France. Marcos M. Flores ayrıca Department of Physics, University of Oslo, Norway afiliyasyonuna sahiptir.

Yazar iletişimleri: Kaynak başlık sayfasında üç yazarın da e-posta adresi verilmiştir. Belge tek bir kişiyi açıkça “corresponding author” olarak etiketlemediği için Verianla tarafından tek bir sorumlu yazar atanmadı.

Eş katkı: Kaynakta eş katkı beyanı tespit edilmemiştir.

Platform: arXiv.

arXiv: 2610.30261v1 [hep-ph].

arXiv DOI: 10.48550/arXiv.2610.30261. Bu DOI arXiv kaydının kalıcı kimliğidir; hakemli dergi DOI'si olarak yorumlanmamalıdır.

arXiv gönderim tarihi: 21 Ekim 2026.

PDF üzerinde yazılı tarih: 23 Ekim 2026.

Tarih notu: PDF başlık sayfasındaki 23 Ekim 2026 tarihi ile arXiv v1 gönderim kaydındaki 21 Ekim 2026 aynı değildir. Her iki bilgi kendi bağlamıyla korunmuştur.

Kaynak türü: Araştırma preprinti; teorik ve nümerik yüksek enerji fiziği / astrofiziksel fenomenoloji çalışması.

Hakemlik durumu: 4 Kasım 2026 itibarıyla resmî arXiv kaydında hakemli bir dergi yayınına ait journal reference gösterilmemektedir. Bu nedenle kaynak hakemli makale gibi sunulmamıştır.

Ana arXiv kategorisi: High Energy Physics – Phenomenology (hep-ph).

Çapraz kategoriler: Cosmology and Nongalactic Astrophysics (astro-ph.CO) ve High Energy Astrophysical Phenomena (astro-ph.HE).

Çalışma uzunluğu: 48 sayfa, 13 şekil.

Temel yöntem: Kısmi-dalga ayrışımı, radyal Schrödinger denklemi, variable-phase birinci-mertebe sistemi, outward regular çözüm, inward irregular çözüm, Sommerfeld faktörleri, faz kaymaları, Green-fonksiyonu subtraction işlemi, Yukawa potansiyeli, Born/WKB/Coulomb analitik kontrolleri ve kısmi-dalga truncation analizi.

Kod: YukawaVPS adlı Mathematica kod deposu çalışma ile birlikte açık olarak yayımlanmıştır.

Veri: “YukawaVPS data: phase shifts, Sommerfeld factors, and cross sections for Yukawa scattering” veri seti, Zenodo DOI 10.5281/zenodo.23475703 altında arşivlenmiştir.

Finansman: Çalışma Avrupa Birliği Horizon 2020 araştırma ve yenilik programı, Grant Agreement No. 101002846, ERC CoG CosmoChart tarafından desteklenmiştir. Marcos M. Flores ayrıca Norwegian Research Council FRIPRO Grant, Project ID 353561 “DarkTurns” desteğini belirtmektedir.

Çıkar çatışması: Kaynak PDF'de ayrı bir çıkar çatışması beyanı tespit edilmemiştir; Verianla kaynakta bulunmayan bir beyan üretmemiştir.

CRediT / yazar katkıları: Kaynak PDF'de ayrı CRediT katkı tablosu veya ayrıntılı yazar-katkı beyanı tespit edilmemiştir.

Lisans: Kaynak PDF'nin yeniden kullanımına ilişkin belirli bir Creative Commons lisansı güvenilir biçimde doğrulanmadığından açık erişim, açık lisans olarak yorumlanmamıştır. Kaynak grafikler birebir yeniden kullanılmamıştır.

Görsel yöntemi: Verianla içindeki yöntem şeması kaynakta açıklanan matematiksel ilişkiler temelinde bağımsız olarak hazırlanmış ölçeksiz SVG'dir; kaynak çalışmanın 13 şeklinin grafik tasarımının kopyası değildir. Kaynak grafiklerindeki ham veriler PDF'den tahmin edilerek yeni veri setine dönüştürülmemiştir.

Temel bilimsel sınır: Çalışma nonrelativistik Schrödinger saçılma problemini çözmektedir. Sonuçlar belirli potansiyeller ve parametreler için teorik/nümerik hesaplamalardır; doğrudan astrofiziksel gözlem, parçacık detektörü sonucu veya karanlık maddenin mikrofiziğine ilişkin deneysel doğrulama değildir.

Genellenebilirlik sınırı: Variable-phase denklemleri Eşitlik (2.7)'deki gerçek, yerel, merkezî potansiyel sınıfına uygulanabilir. \(T(k;k_\diamond)\)'nin extraction işlemi ayrıca Eşitlik (2.22)'deki uygun kısa-mesafe davranışını gerektirir.