Mushtaq Ahmad tərəfindən hazırlanmış SHEET-AD, çarpışmasız soyuq qaranlıq maddənin kosmoloji Vlasov–Poisson denklemleri altında evrimini, bir-birindən müstəqil parçacıqlar əvəzinə altı ölçülü faz məkanına yerləşdirilmiş bşəbəkəlantılı bir Lagrange təbəqəsı olarak hesaplamayı hədəfləyən sayısal bir metoddir. Tədqiqatçı; yüksək mərtəbəli təbəqə geometriyasıni, kütləni qoruyan təbəqədən ızgaraya aktarımı, eyni ayrıq qravitasiya enerjisindən törədilən qüvvəni, kick–drift–kick zaman inteqrasiyasını, kaustiklərə həssas adaptiv incəltməni və kaydedilmiş sabit incəltmə qrafları üzerinden çalışan əks rejim törəmələrini tek bir hesaplama zincirinde birleştirmiştir. Metod operator və küçük ölçekli üç boyutlu müqayisəlarda doğrulanmış olsa da nəticəlar büyük hacimli kosmoloji üretim simülasyonu değil, kompakt və meyare özgü bir yenidən istehsal edilə bilmə kıyaslamasına dayanmaktadır.

Operatör kontrollerinde toplam kütle xətası sürüşən nöqtə yuvarlama səviyyəinde kalmış; təbəqədən ızgaraya saçma ilə ızgaradan təbəqəyə toplama operatorlerinin transpozisiya qüsuru 4,46×10⁻¹⁶, sabit qraf üzerindeki simplektik qüsur 2,22×10⁻¹⁶ və ayrık adjoint qradiyent xətası 1,62×10⁻¹⁰ olarak ölçülmüştür. Yüksek mərtəbəli serbest akış testlerinde p=1, p=2 və p=3 üçün sırasıyla 1,96, 2,97 və 3,97 yaxınlaşma meylləri elde edilmiş; kick–drift–kick inteqrasiyau ikinci mərtəbə yaxınlaşma göstermiştir.

48³ hücreli qüvvə ızgarası kullanılan kompakt üç boyutlu kıyaslamada, indikatorye dayalı yüzde 10–20 adaptiv hesablamalar çökmüş bölgələri çəkiləndirən sıxlıq xətasını PM 28³ metodine görə 31–33 kat, PM 32³ metodine görə 11–12 kat azaltmıştır. İki seviyeli incəltməde bu nisbətlar təxminən 39 və 14 kata yükselmiştir. Ancak bu üstünlük yalnız yüksək sıxlıqlu, kaustiklərə həssas bölgelere odaklanan özel meyar üçün geçerlidir. Ağırlıksız qlobal sıxlıq normasında yeterince yüksək parçacıq sayısına sahip geleneksel parçacıq-şəbəkə metodleri daha başarılı və ya rekabetçi kalmıştır.

Araştırma hangi hesaplama sorununu ele almaktadır?

Kozmik yapıların oluşumu, erken evrendeki küçük sıxlıq fərqlilıklarının kendi qravitasiyaleri altında büyüyerek galaksi, galaksi kümesi, filament və boşluklardan oluşan büyük ölçekli yapıya dönüşmesiyle açıklanır. Soğuk qaranlıq maddə parçacıqlarının çarpışmalarının ihmal edilebildiği sınırda bu evrim, kosmoloji Vlasov–Poisson sistemiyle temsil edilir.

Yaygın N-cisim və parçacıq-şəbəkə metodleri, müddətkli qaranlıq maddə dşəbəkəılımını sonlu sayıdaki parçacıqla örnekler. Bu yaklaşım büyük kosmoloji hacimlerde son derece güçlüdür; lakin başlangıçta soğuk olan dşəbəkəılımın geometrik yapısını birbaşa temsil etmez. Soğuk dşəbəkəılım, gerçekte birbirinden müstəqil noktalardan oluşan bir bulut değil, altı ölçülü faz məkanına yerləşdirilmiş üç boyutlu və bşəbəkəlantılı bir Lagrange təbəqəsıdır.

Qravitasiyasel çöküşten sonra fərqli akışların aynı fiziksel mövqea ulaşmasıyla qabıq keçidi gerçekleşir. Təbəqənın üç boyutlu mövqe məkanına izdüşümü katlanır; kaustikler və çox axınlı bölgeler oluşur. Buna rşəbəkəmen mövqe və impulsun birlikte ele alındığı altı boyutlu faz məkanındaki təbəqə bşəbəkəlantılı kalır. SHEET-AD, bu bşəbəkəlantılı geometrinin birbaşa evrimini hesaplayarak kaustik mövqelarını, akış sayılarını, ince ölçekli sıxlıqları və çox axınlı hız momentlerini aynı temsilden çıkarmayı hedeflemektedir.

Kosmoloji Vlasov–Poisson sistemi nasıl ifade edilmektedir?

Araştırmada hareket denklemleri genel bir monoton zaman değişkeni s kullanılarak şu katsayılı biçimde verilmiştir:

\[ \frac{\partial f}{\partial s} +A(s)\,\mathbf{p}\cdot\nabla_x f -B(s)\,\nabla_x\phi\cdot\nabla_p f=0 \]

\[ \nabla_x^2\phi=C(s)\,\delta \]

• f(x,p,s): Konum, impuls və zaman değişkenlerine bşəbəkəlı faz-məkanı kütle yoğunluğudur.
• x: Periyodik kosmoloji kutudaki eş hareketli mövqedur.
• p: Kanonik peculyar impulsdur.
• φ: Qravitasiya potansiyelidir.
• δ=(ρ−ρ̄)/ρ̄: Ortalama yoğunluğa görə sıxlıq kontrastıdır.
• A(s), B(s), C(s): Seçilen kozmoloji və zaman değişkenine bşəbəkəlı katsayılardır.

Karakteristik yörüngeler şu denklemleri izler:

\[ \frac{d\mathbf{X}}{ds}=A(s)\mathbf{P} \]

\[ \frac{d\mathbf{P}}{ds}=-B(s)\nabla_x\phi(\mathbf{X},s) \]

Tədqiqatçı hem fiziksel zaman t hem de ölçek faktörü a üçün gerekli katsayıları tanımlamış; temel sayısal deneylerde zaman değişkeni olarak ölçek faktörünü kullanmıştır. Böylece uzamsal ayrıklaştırma, kosmoloji zaman parametrizasyonundan ayrılmıştır.

Soğuk qaranlıq maddə neden faz-məkanı təbəqəsi olarak gösterilmektedir?

Soğuk sınırda dşəbəkəılım, Lagrange koordinatları q ilə etiketlenen üç boyutlu bir çokkatlı üzerinde yoğunlaşır:

\[ f(\mathbf{x},\mathbf{p},s) =\int_{\Omega_q}\mu(\mathbf{q}) \,\delta_D(\mathbf{x}-\mathbf{X}(\mathbf{q},s)) \,\delta_D(\mathbf{p}-\mathbf{P}(\mathbf{q},s))\,d\mathbf{q} \]

Burada μ(q)dq sabit Lagrange kütlesini,

\[ \mathbf{Z}(\mathbf{q},s) =\bigl(\mathbf{X}(\mathbf{q},s),\mathbf{P}(\mathbf{q},s)\bigr) \in\mathbb{R}^6 \]

ise üç boyutlu etiket uzayının altı boyutlu faz məkanına gömülmesini ifadə edir. Sayısal metodin temel bilinmeyeni tek tek parçacıq yörüngeleri değil, bu müddətkli Z(q,s) haritasıdır.

Kaustik və çox axınlı bölgeler nasıl ortaya çıkmaktadır?

Konum haritasının Lagrange koordinatlarına görə türevi:

\[ \mathbf{F}(\mathbf{q},s) =\frac{\partial\mathbf{X}}{\partial\mathbf{q}} \]

və bunun determinantı:

\[ J(\mathbf{q},s)=\det\mathbf{F}(\mathbf{q},s) \]

olarak müəyyən edilmişdir. Aynı Euler mövqeuna ulaşan düzenli ön görüntüler q_ℓ üçün sıxlıq:

\[ \rho(\mathbf{x},s) =\sum_{\mathbf{q}_{\ell}\in\mathbf{X}^{-1}(\mathbf{x},s)} \frac{\mu(\mathbf{q}_{\ell})}{|J(\mathbf{q}_{\ell},s)|} \]

biçimindedir. Toplamdaki terim sayısı o mövqedaki axın sayını verir. J=0 olduğunda mövqe izdüşümü katlanır və kaustik oluşur. Bu vəziyyət təbəqənın fiziksel olarak yırtılması mənasına gəlmir.

Faz-məkanı düzenliliği daha geniş toxunan matrisiyle kontrol edilir:

\[ \mathcal{G} =\frac{\partial\mathbf{Z}}{\partial\mathbf{q}} =\begin{pmatrix} \partial\mathbf{X}/\partial\mathbf{q}\\ \partial\mathbf{P}/\partial\mathbf{q} \end{pmatrix} \]

\[ \det(\mathcal{G}^{\mathsf{T}}\mathcal{G})>0 \]

olduğu müddətce faz-məkanı təbəqəsi düzenli bir gömülü yapı olarak kalır. Dolayısıyla SHEET-AD fiziksel bir kaustiği sayısal element bozulması olarak dəyərlendirmez; yalnız faz-məkanı teğetinin rank kaybını və ya təbəqənın yetersiz çözülmesini gerçek bir sayısal uyarı kabul eder.

Şəkil 2 qatlanmayı nasıl açıklamaktadır?

Şəkil 2’de bir boyutlu örnek üzerinden dört fərqli görünüm verilmiştir. İlk panel, qabıq keçidinden sonra bile faz məkanında bşəbəkəlantılı kalan təbəqəyı və kaustik noktalarını gösterir. İkinci panelde:

\[ X(q)=q-1{,}18\sin q \]

haritasının qatlanmış biçimi görülür. Üçüncü panel, aynı təbəqənın kübik B-spline kütle aktarımından sonra oluşan yoğunluğunu; dördüncü panel ise J=∂X/∂q dəyərini və tam sayı axın sayını gösterir. J’nin sıfır olduğu kesikli çizgilerde sıxlıq tepe yaparken faz-məkanı təbəqəsi bşəbəkəlantısını korur.

Yüksek mərtəbəli təbəqə ayrıklaştırması nasıl kurulmuştur?

Lagrange alanı adaptiv bir şəbəkədaki elemanlara ayrılmıştır. Temel uygulamada üç boyutlu tenzor-hasil heksaedral elemanlar və derece p=3 polinomlar istifadə edilmişdir. Her elemandaki mövqe və impuls alanları:

\[ \mathbf{X}_h(\mathbf{q},s)|_K =\sum_{A=1}^{N_p(K)} \mathbf{X}_{K,A}(s)\,\ell_{K,A}(\mathbf{q}) \]

\[ \mathbf{P}_h(\mathbf{q},s)|_K =\sum_{A=1}^{N_p(K)} \mathbf{P}_{K,A}(s)\,\ell_{K,A}(\mathbf{q}) \]

biçiminde interpolasyonla göstərilmişdir. Katsayılar kütle taşıyan müstəqil parçacıqlar değil, müddətkli təbəqənın kontrol dəyərləridir. Kütle, kaustik, sıxlıq və hız momentleri elemanlar üzerinde inteqrasiya yapılarak elde edilir.

Düyün bazlı temsil eleman yüzlerinde müddətklidir. Asılı düyün kısıtları cebirsel olarak ortadan kaldırılır. Bu tercih təbəqənın adaptiv incəltmə sırasında yapay biçimde yırtılmasını önler.

Lagrange kütlesi və kanonik değişkenler nasıl müəyyən edilmişdir?

Her elemandaki inteqrasiya nöqtəsi üçün sabit bir kütle tanımlanır:

\[ m_{K,r}=w_{K,r}\,\mu(\mathbf{q}_{K,r}) \]

\[ M_h=\sum_{K,r}m_{K,r} \]

Lagrange kütle matrisi:

\[ (M_q)_{ij} =\sum_{K,r} m_{K,r}N_i(\mathbf{q}_{K,r})N_j(\mathbf{q}_{K,r}) \]

olarak oluşturulur. Düyün impuls alanı P ilə gerçek kanonik impuls arasında:

\[ \boldsymbol{\Pi}=M_q\mathbf{P} \]

əlaqəsi kurulmuştur. Simplektik inteqrasiya kanonik değişkenler olan X və Π üzerinden tanımlanır. Bu ayrım önemlidir; interpolasyonla gösterilen impuls alanı ilə mövqea eşlenik sonlu boyutlu kanonik impuls aynı nesne değildir.

Adaptiv incəltmə sırasında geometri nasıl izlenmektedir?

Tədqiqatçı yalnız mövqe məkanındaki yoğunluğa bakmak əvəzinə çok sayıda geometrik tanı kullanmıştır:

• J_h: Konum izdüşümünün qatlanmasını və kaustik komşuluğunu gösterir.
• Γ_h=𝒢_h^T𝒢_h: Faz-məkanı təbəqəsının gömülme kalitesini izler.
• π_K: Təbəqəya çekilmiş simplektik yapının artık dəyərini ölçer.
• κ_K: Çözülmemiş təbəqə eğriliğini və bükülmesini ölçer.
• c_K: Eleman içindeki en küçük və en büyük |J| dəyərlərinin nisbətından sıkışmayı dəyərlendirir.
• η_ρ,K: Sıxlıq projeksiyonunun inteqrasiya mərtəbəsine duyarlılığını ölçer.
• η_F,K: Təbəqədan ızgaraya və tekrar təbəqəya giden bütün qüvvə aktarım zincirinin duyarlılığını ölçer.
• η_t,K: Bir tam adımla iki yarım adım arasındaki fərqtan zaman xətasını təxmin eder.
• η_R,K: Yeniden şəbəkə oluşturma və ya kobudlaşdırma sırasında oluşan aktarım kusurunu ölçer.

Kaustiğe yakın olmak tek başına incəltmə sebebi değildir. Kaustik komşuluğu toleransları sıkılaştırır; lakin eleman yalnız geometri, projeksiyon, qüvvə, zaman və ya yenidən şəbəkə qurma xətası belirlenen sınırı aştığında incəldilir. Böylece fiziksel sıxlıq tekilliği ilə sayısal çözünürlük kaybı birbirinden ayrılır.

Yeniden şəbəkə oluşturma nasıl kaydedilmektedir?

Bir ana eleman çocuk elemanlara ayrıldığında yeni düyün dəyərləri yüksək mərtəbəli polinom uzatmasıyla oluşturulur:

\[ \mathbf{Z}_{\mathrm{child}} =R_{K\rightarrow\{K_c\}}\mathbf{Z}_{\mathrm{parent}} \]

Çocuk inteqrasiya kütleleri toplamları ana elemanın kütlesine eşit olacak biçimde yerli olarak normalleştirilir. Kabalaştırma kullanılacaksa kardeş elemanlar yalnız eğrilik, projeksiyon duyarlılığı və kaustik indikatorleri küçük olduğunda; daha aşşəbəkəı boyutlu yeniden yapılandırma toleransı sşəbəkəlandığında birleştirilir.

Yeniden şəbəkə oluşturma sonlu boyutlu uzayın boyutunu değiştirdiği üçün metod bu olaylar boyunca tam simplektiklik iddia etmez. Her aktarımda yerli faz-məkanı temsil kusuru və kanonik impuls fərqı kaydedilir. Saf incəltməde bu kusurlar inteqrasiya və yuvarlama səviyyəinde kalabilir; kobudlaşdırmada ise atılan bilgiyi ölçer.

Kütle təbəqədan Euler ızgarasına nasıl aktarılmaktadır?

Qravitasiya alanı düzenli bir Euler ızgarasında hesaplanır. Her təbəqə inteqrasiya nöqtəsindaki kütle, kompakt destekli və birliği parçalayan şəbəkəırlıklarla komşu hücrelere dşəbəkəıtılır:

\[ \omega_g(\mathbf{x})\geq 0 \]

\[ \sum_g\omega_g(\mathbf{x})=1 \]

Temel uygulamada üç boyutlu tensor-çarpım kübik kardinal B-spline istifadə edilmişdir. Izgara yoğunluğu:

\[ \rho_g(\mathbf{X}_h) =\frac{1}{\Delta V} \sum_{K,r} m_{K,r}\,\omega_g(\mathbf{X}_{K,r}) \]

olarak hesaplanır. Ağırlıkların toplamının tam olarak bire eşit olması nedeniyle:

\[ \sum_g\rho_g\Delta V =\sum_{K,r}m_{K,r}=M_h \]

olur. Bu xüsusiyyət toplam kütleyi ayrık səviyyəde, sürüşən nöqtə yuvarlama xətası dışında tam olarak korur. Pozitif kütleler və negatif olmayan aktarım şəbəkəırlıkları aynı zamanda negatif ızgara yoğunluğu oluşmasını önler.

Şəkil 3 hangi operator tutarlılığını göstərir?

Şəkil 3’ün ilk panelinde təbəqə inteqrasiya kütlelerinin birliği parçalayan şəbəkəırlıklarla ızgaraya dşəbəkəıtılması gösterilir. İkinci panel sıfır ortalamalı Poisson kaynşəbəkəını, üçüncü panel ise operator zincirini özetler:

\[ \text{Tabaka kütlesi} \xrightarrow{S} \text{Izgara yoğunluğu} \xrightarrow{L_h^{-1}} \text{Potansiyel} \xrightarrow{-\partial_XV_h} \text{Tabaka kuvveti} \]

Ters geçişte aynı S operatorünün transpozisi, aynı özqoşma Poisson tersi və aynı şəbəkəırlık türevleri kullanılır. Fiziksel qüvvə ilə əks rejim türevinin ayrı, uyumsuz operatorlerle hesaplanmaması metodin temel tasarım ilkelerinden biridir.

Poisson denklemi nasıl çözülmektedir?

Sıfır ortalamalı sıxlıq kontrastı üçün ayrık Poisson denklemi:

\[ L_h\phi_h=C(s)\delta_h \]

olarak yazılmıştır. Periyodik kutuda hızlı Fourier dönüşümü uygulanır. İkinci mərtəbə periyodik Laplace operatorünün negatif sembolü:

\[ \lambda_{\mathbf{k}} =\frac{4}{\Delta x^2} \sum_{\alpha=1}^{3} \sin^2\left(\frac{\pi k_\alpha}{N_g}\right) \]

və sıfır dışındaki Fourier kipleri için:

\[ \widehat{\phi}_{\mathbf{k}} =-\frac{C(s)}{\lambda_{\mathbf{k}}} \widehat{\delta}_{\mathbf{k}} \]

kullanılır. Sıfır kip φ̂₀=0 yapılır. Aynı Fourier normalizasyonu, sıfır mod kuralı və Laplace sembolü ileri hesapta, tanı testlerinde və əks rejim geçişinde aynen korunur.

Qüvvə neden aynı ayrık enerjiden türetilmektedir?

Tədqiqatçı ayrıq qravitasiya enerjisini:

\[ V_h(\mathbf{X}_h;s) =\frac{1}{2} \langle \rho_h-\bar{\rho}_h,\phi_h\rangle_g \]

olarak tanımlamıştır. Öz-eşlenik Poisson operatorü sayesinde enerjinin mövqea görə ilk değişimi birbaşa potansiyel ilə sıxlıq değişiminin iç çarpımına dönüşür. Bir inteqrasiya nöqtəsindaki qüvvə:

\[ \mathbf{a}_{K,r} =-\nabla\Phi_W(\mathbf{X}_{K,r}) \]

\[ \nabla\Phi_W(\mathbf{x}) =\sum_g\phi_g\nabla\omega_g(\mathbf{x}) \]

olarak elde edilir. Düyün qüvvəi, inteqrasiya noktalarındaki bu dəyərlərin təbəqə baz fonksiyonlarıyla toplanmasıyla hesaplanır.

Yaygın bir alternatif, önce ızgarada −∇φ qüvvəini sonlu fərqla hesaplayıp ardından bunu təbəqəya taşımaktır. Ancak bu müstəqil işlem, kullanılan ayrık enerjinin tam mövqe türevi olmak zorunda değildir. SHEET-AD varsayılan metodde qüvvəi birbaşa aynı enerjiden türeterek Hamilton yapısını, Hessian simetrisini və adjoint tutarlılığını korumayı amaçlamıştır.

Simplektik zaman inteqrasiyası nasıl çalışmaktadır?

Sabit incəltmə qrafiğinde sonlu boyutlu Hamiltonian:

\[ H_h(\mathbf{X},\boldsymbol{\Pi};s) =\frac{1}{2}A(s) \boldsymbol{\Pi}^{\mathsf{T}}M_q^{-1}\boldsymbol{\Pi} +\Lambda(s)\widetilde{V}_h(\mathbf{X}) \]

biçimindedir. Burada Λ(s)=B(s)C(s)’dir. Kinetik bölüm yalnız Π’ye, potansiyel bölüm yalnız X’e bşəbəkəlıdır. Bu səbəbdən hareket tam çözülebilen drift və kick alt akışlarına ayrılabilir.

Sürükleme katsayısı:

\[ \mathcal{D}_{a,b} =\int_{s_a}^{s_b}A(u)\,du \]

kick katsayısı:

\[ \mathcal{K}_{a,b} =\int_{s_a}^{s_b}\Lambda(u)\,du \]

olarak tanımlanır. Bir zaman addımı simetrik kick–drift–kick düzeninde uygulanır:

\[ \boldsymbol{\Pi}^{n+1/2} =\boldsymbol{\Pi}^{n} -\mathcal{K}_{n,n+1/2} \nabla_X\widetilde{V}_h(\mathbf{X}^{n}) \]

\[ \mathbf{X}^{n+1} =\mathbf{X}^{n} +\mathcal{D}_{n,n+1} M_q^{-1}\boldsymbol{\Pi}^{n+1/2} \]

\[ \boldsymbol{\Pi}^{n+1} =\boldsymbol{\Pi}^{n+1/2} -\mathcal{K}_{n+1/2,n+1} \nabla_X\widetilde{V}_h(\mathbf{X}^{n+1}) \]

Her tam adımda iki qüvvə hesabı yapılır. Qraf, inteqrasiya kuralları və Euler ızgarası değişmediği müddətce bu harita kanoniktir və ayrık simplektik iki-formu korur. Yeniden şəbəkə oluşturma olayları ayrı muhafazakâr projeksiyonlar olarak ele alınır.

Simplektiklik ne anlama gelmektedir?

Kosmoloji değişkenlerde Hamiltonian açıkça zamana bşəbəkəlı olduğu üçün mekanik enerjinin sabit kalması beklenmez. Buradaki simplektiklik, sayısal haritanın kanonik faz-məkanı geometrisini koruması mənasına gəlir:

\[ (\Phi_n^h)^*\Omega_h=\Omega_h \]

\[ \Omega_h =\sum_{i,\alpha} dX_{i,\alpha}\wedge d\Pi_{i,\alpha} \]

Bu xüsusiyyət sabit qraf aralıklarında geçerlidir. Uzun dönem davranış Hamiltonianin sabitliğinden çok ileri–geri çevrilebilirlik, impuls sürüklenmesi, zaman addımı yaxınlaşması və zaman bşəbəkəımlılığını hesaba katan iş dengeli Hamilton artığıyla dəyərlendirilmiştir.

Adaptiv indikatorler nasıl birleştirilmektedir?

Her eleman üçün boyutsuz toplam indikator:

\[ \eta_K^2 =w_Z\widehat{\eta}_{Z,K}^2 +w_\rho\widehat{\eta}_{\rho,K}^2 +w_F\widehat{\eta}_{F,K}^2 +w_t\widehat{\eta}_{t,K}^2 +w_R\widehat{\eta}_{R,K}^2 \]

olarak oluşturulur. Şapkalı dəyərler kullanıcı tarafından önceden belirlenen toleranslara bölünmüş bileşenlerdir. Bir eleman η_K>1 olduğunda və ya faz-məkanı gömülme kalitesi belirlenen alt sınırın altına indiğinde işaretlenir.

İkinci aşamada işaretlenen elemanlar, toplam kütle çəki verilmiş indikator xətasının belirli bir bölümünü karşılamıyorsa en büyük katkıya sahip ek elemanlar da incəldilir. Böylece işlem yalnız en yoğun hücrelere değil, hesaplanan xəta mekanizmasına en fazla katkı veren Lagrange bölgelerine yönlendirilir.

Şəkil 6 indikatornin xətayı nasıl yerlileştirdiğini göstərir?

Şəkil 6’da 64³ təbəqə referansı, 48³ qüvvə ızgarası və 20 KDK adımı kullanılan kompakt kıyaslamadan bir kesit gösterilir. Referans sıxlıq, sıkışma/eğrilik indikatorsi, seçilen hücrelerin konturları və sıxlıq artığı yan yana verilmiştir.

Kümülatif lokallaşdırma əyrisine görə indikator sıralamasındaki en yüksək yüzde 5’lik hücre grubu çəki verilmiş xətanın təxminən yüzde 85’ini, yüzde 10’luk grup yüzde 86’sını və yüzde 20’lik grup yüzde 87’sini kapsamaktadır. Rastgele hücre seçimi aynı xəta payını yakalayamamaktadır. Bu nəticə adaptiv metodin neden yalnız sıxlıq eşiğine görə incəltmə yapmadığını göstərir.

Metod nasıl türevlenə bilən hâle getirilmiştir?

Adaptiv ileri hesap tamamlandığında kabul edilen şəbəkəlar, zaman adımları, inteqrasiya mərtəbəleri və yenidən şəbəkə qurma transferleri kaydedilir. Bu kayıt 𝔊 ilə göstərilmişdir. Sabit kayıt şərtu altında bir adım:

\[ \mathbf{Y}^{n+1} =\Phi_n^{\mathfrak{G}}(\mathbf{Y}^{n},\theta) \]

biçiminde sonlu boyutlu və türevlenə bilən bir haritadır. θ; başlangıç Fourier genliklerini, başlangıç yer değiştirme katsayılarını, kosmoloji parametreleri və ya gözlem modeli değişkenlerini içerebilir.

Bir hedef fonksiyon:

\[ \mathcal{J}(\theta) =\Psi(\mathbf{Y}^{N},\theta) +\sum_{n=0}^{N-1}\ell_n(\mathbf{Y}^{n},\theta) \]

olarak tanımlanır. Ters kip adjoint değişkeni:

\[ \lambda^n =\partial_{\mathbf{Y}^n}\ell_n +\left(D_{\mathbf{Y}}\Phi_n^{\mathfrak{G}}\right)^{\mathsf{T}} \lambda^{n+1} \]

bşəbəkəıntısıyla sondan başa taşınır. Parametre qradiyentı da aynı geçiş sırasında biriktirilir.

“Sabit qraf türevi” hangi sınıra sahiptir?

Adaptiv incəltmə kararları eşiklere bşəbəkəlı ayrık kontrol akışlarıdır. SHEET-AD bu kararların türevini hesaplamaz. Bunun əvəzinə kabul edilmiş qraf ardıcıllıqsini dondurur və o ardıcıllıqde gerçekten çalıştırılmış sayısal haritanın türevini hesaplar.

Bu səbəbdən qradiyent, belirli bir qraf içinde düzgün olabilir; lakin parametredeki küçük bir değişiklik kabul edilen qrafı değiştirdiğinde qradiyent sıçrayabilir. Sonlu fərq doğrulamaları da pertürbe edilmiş hesaplarda aynı qraf kaydını yeniden oynatır. Nəticəlar bütün adaptiv karar mekanizmasının qlobal olarak düzgün bir türevini değil, kabul edilmiş sayısal yörüngenin doğru ayrık türevini doğrulamaktadır.

Ters kipte hangi operatorler kullanılmaktadır?

Şəkil 1 ileri və ters zinciri nasıl özetlemektedir?

Şəkil 1’in üst satırı, faz-məkanı təbəqəsindan başlayarak adaptiv Q_p elemanlarına, kütle projeksiyonuna, FFT Poisson çözümüne, variasional qüvvəe və KDK güncellemesine giden ileri hesap zincirini gösterir.

Alt satır ise ters KDK, Hessian–vektör çarpımları, aynı Poisson tersi, S^T toplama operatorü, kaydedilmiş transfer R^T və qraf kaydına doğru ilerleyen ters zinciri gösterir. Aynı renklerin kullanılması, ileri və ters işlemlerin aynı matematiksel operator ailelerine dayandığını vurgular.

Yüksek mərtəbəli transferler nasıl doğrulanmıştır?

Üç boyutlu ana tensor-çarpım elemanları sekiz çocuk elemana ayrılmış; p=1, p=2 və p=3 üçün uzatma matrisleri sınanmıştır. Şəkil 4’te ileri aktarım R və ters geçişte kullanılan R^T geometrik olarak gösterilir.

Dört doğrusal olmayan KDK adımı içine bir p=3 yenidən şəbəkə qurma transferinin yerleştirildiği ek yörünge testinde bütün kaydedilmiş qrafın sonlu fərq qradiyent xətası 1,68×10⁻¹¹, düzeltilmiş Taylor meyli 2,00 olarak bulunmuştur.

Katlanmış təbəqəda yüksək mərtəbənin etkisi nedir?

Şəkil 5, X(q)=q−1,18 sin q qatlanmış haritasını p=1, p=2 və p=3 elemanlarla müqayisəktadır. Düzgün təbəqə haritası xətasında ölçülen eğimler sırasıyla 1,99, 3,00 və 4,00’dır; bunlar beklenen p+1 davranışıyla uyumludur.

Ham projekte edilmiş sıxlıq xətasında eğimler daha aşşəbəkəıtür: p=1 üçün 0,97, p=2 üçün 0,81 və p=3 üçün 0,85. Bunun nedeni yüksək mərtəbəli təbəqə haritasının başarısız olması değil, kaustikteki tekil sıxlıq yapısının qlobal ham sıxlıq normasınun yaxınlaşmasını sınırlandırmasıdır. Aynı eleman sayısında p=3 yine en aşşəbəkəı sıxlıq xətasını vermiştir.

Varyasyonel qüvvə seçimi neden önemlidir?

Tədqiqatçı, enerji türevinden elde edilen varsayılan qüvvəi geleneksel ızgara qüvvəi interpolasyonuyla karşılaştırmıştır:

Geleneksel alternatif aynı ızgarada hesaplanabilmesine rşəbəkəmen kullanılan ayrıq qravitasiya enerjisinin tam qradiyentı değildir. Kusurların birkaç mərtəbə büyümesi, simplektik və adjoint testlerinde enerjiyle tutarlı qüvvə kullanımının önemini göstərir.

Üç boyutlu kompakt kıyaslama nasıl tasarlanmıştır?

• Referans faz-məkanı təbəqəsi: 64³ düzenli örnek
• Ortak qüvvə ızgarası: 48³
• Zaman inteqrasiyau: T=0,5’e kadar 20 KDK adımı
• Parçacıq-şəbəkə temelleri: 24³, 28³, 32³, 36³ və 40³ parçacıq
• Adaptiv təbəqə başlangıcı: 28³ makro hücre
• Bir seviyeli incəltmə: indikator sıralamasındaki en yüksək yüzde 5, yüzde 10 və yüzde 20
• İki seviyeli incəltmə: birinci seviye çocuklarının en yüksək yüzde 25’inin yeniden inceltilmesi
• Kontrol metodi: yüzde 10 indikator hesaplamasıyla aynı aktif işaretleyici sayısına sahip sıxlıq eşiği tabanlı incəltmə
• Aktarım çekirdeği: kübik B-spline
• Poisson çözümü: periyodik FFT, sıfır mod çıkarılmış
• Ana təsadüfi toxum: 91217
• Sşəbəkəlamlık toxumları: 91217, 2027 və 31415

Bu deney büyük gözlem hacimli bir evren simülasyonu değil, indikator davranışını, operator zincirini və çökmüş bölgə meyarlarıni sınamak üçün tasarlanmış kompakt bir sayısal kıyaslamadır.

Parçacıq-şəbəkə və adaptiv təbəqə nəticəları nelerdir?

Nəticəlar neden meyare görə fərqli görünmektedir?

Ağırlıksız qlobal sıxlıq norması bütün hesaplama hacmindeki xətayı aynı biçimde ölçer. Düzenli PM parçacıqları bütün hacme yayıldığı üçün parçacıq sayısı arttıkça bu normda düzenli biçimde iyileşmiştir. PM 40³, 4,421×10⁻² qlobal xəta ilə adaptiv metodlerden belirgin ölçüde daha aşşəbəkəı dəyər vermiştir.

Çökmüş bölge çəki verilmiş norm ise yüksək sıxlıqlu hücrelere çok daha büyük şəbəkəırlık verir:

\[ E_{\rho,\mathrm{hd}} =\left[ \frac{\sum_gw_g(\delta_g-\delta_{\mathrm{ref},g})^2} {\sum_gw_g\delta_{\mathrm{ref},g}^2} \right]^{1/2} \]

\[ w_g =\left(1+\max(\delta_{\mathrm{ref},g},0)\right)^3 \]

İndikatorye dayalı adaptiv metod çözünürlüğü tam olarak bu çökmüş və kaustiklərə həssas bölgelere yoğunlaştırdığı üçün söz konusu meyarte büyük kazanç sşəbəkəlamıştır. Bunun bedeli, hacmin geri kalanına daha az örnek ayrılması və qlobal sıxlıq normasınun daha yüksək kalmasıdır.

Yüzde 10 adaptiv metod hangi müqayisəlarda öne çıkmıştır?

Yüzde 10 adaptiv şərtda E_ρ,hd=2,457×10⁻³ bulunmuştur. Bu dəyər:

• PM 28³ dəyərinden təxminən 31 kat daha aşşəbəkəıtür.
• PM 32³ dəyərinden təxminən 11 kat daha aşşəbəkəıtür.
• PM 36³ dəyərinden təxminən 5 kat daha aşşəbəkəıtür.
• Aynı aktif işaretleyici sayısına sahip sıxlıq eşiği metodinden təxminən 4,2 kat daha aşşəbəkəıtür.

Bu son müqayisə, başarının yalnız en yoğun hücreleri incəltməkten mənbəlanmadığını; təbəqə sıkışması, eğrilik və projeksiyon duyarlılığını birleştiren indikatornin fərqli bilgi taşıdığını desteklemektedir.

Şəkil 7 alan səviyyəindeki fərqları nasıl göstərir?

Şəkil 7’de referans, PM 32³, yüzde 10 adaptiv və iki seviyeli adaptiv sıxlıq kesitleri aynı renk ölçeğinde gösterilir. Alt sıradaki yakınlaştırılmış xəta haritaları, referansın en yüksək sıxlıqlu yüzde birlik bölgesinde PM metodinin daha geniş artıklar üretirken adaptiv hesablamaların bu hedef bölgenin şeklini daha iyi yakaladığını gösterir.

Şəkildeki meyar kartları aynı zamanda önemli karşıtlığı açıkça sunar: PM 32³ qlobal sıxlıq normasında daha iyi, adaptiv metodler ise çökmüş bölgə çəki verilmiş normda daha iyidir. Tədqiqatçı bu ayrımı gizlememiş və metodi evrensel olarak üstün ilan etmemiştir.

Ölçüt duyarlılığı nasıl sınanmıştır?

Ağırlık üssü γ=2, 3 və 4 olarak değiştirildiğinde yüzde 10 adaptiv metod PM 32³ metodine görə təxminən 2,5, 11,0 və 14,2 kat daha aşşəbəkəı xəta vermiştir. Sıxlıq şəbəkəırlığı arttıkça adaptiv metodin yüksək sıxlıqlu bölgelere odaklanan avantajı daha görünür hâle gelmiştir.

Referans hücrelerin en yoğun yüzde 0,5, yüzde 1 və yüzde 2’sini içeren maskelerde yüzde 10 adaptiv metodin xətaları sırasıyla 1,39×10⁻², 1,44×10⁻² və 1,67×10⁻² olmuştur. PM 32³ üçün karşılık gelen dəyərler 5,61×10⁻², 5,74×10⁻² və 6,03×10⁻² olarak bildirilmişdir.

Spektral davranış ne göstərir?

Sıxlıq aktarım fonksiyonu:

\[ T(k)=\sqrt{\frac{P_{hh}(k)}{P_{rr}(k)}} \]

və çapraz korelasyon:

\[ r(k)=\frac{P_{hr}(k)} {\sqrt{P_{hh}(k)P_{rr}(k)}} \]

ile dəyərlendirilmiştir. Yüzde 10 adaptiv şərtun ilk beş spektral kutusunda |T−1| en fazla 1,5×10⁻², en aşşəbəkəı r dəyəri 0,992 olmuştur.

Buna karşılık k≥19,15 yüksək dalga sayısı bölgesinde adaptiv metodin T dəyəri 1,07–2,69, r dəyəri 0,414–0,942 aralığına yayılmıştır. PM 32³ aynı bölgede T=0,967–1,25 və r=0,878–0,963 ilə daha düzenli spektral davranış göstermiştir. Yüksek k fərqlilıkları, adaptiv örnek yoğunlaşması, B-spline yumuşatması, sınırlı qüvvə ızgarası və 64³ referans çözünürlüğünün ortak etkisini taşımaktadır.

Farklı rastgele başlangıçlarda nəticə korundu mu?

40³ qüvvə ızgarasında yapılan üç toxumlu stres testinde yüksək sıxlıq çəki verilmiş xəta dəyərləri şöyledir:

İndikator adaptiv metod üç toxumun tamamında diğer iki metodden daha aşşəbəkəı çökmüş bölgə xətası vermiştir. Ancak üç toxum geniş kapsamlı istatistiksel genelleme üçün sınırlı bir örneklem olduğundan tədqiqatçı nəticəları sşəbəkəlamlık kontrolü olarak sunmuş, asimptotik performans kanıtı olarak yorumlamamıştır.

Adjoint doğrulaması hangi nəticəları vermiştir?

Temel deterministik operator doğrulamaları şöyledir:

Şəkil 11’de düzeltilmemiş Taylor artığının birinci mərtəbə, adjoint ilə düzeltilmiş artığın ise uygun sonlu fərq aralığında ikinci mərtəbə azaldığı görülmektedir. Sekiz başlangıç Fourier genliğini geri kazanma örneğinde Adam optimizasyonu 80 yineleme boyunca kaybı təxminən 10⁻⁵ səviyyəine düşürmüş; geri kazanılan və hedef genlikler üçün korelasyon 1,000000 olarak göstərilmişdir.

Bu örnek, sabit qraf qradiyentlarının bir ters problemde kullanılabildiğini gösterir; lakin büyük ölçekli kosmoloji başlangıç şərtu çıkarımı və ya gözlemsel veri analizi değildir.

Hesaplama maliyeti nasıl ölçeklenmektedir?

Bir qüvvə dəyərlendirmesinin temel maliyeti:

\[ W_{\mathrm{force}} =O(N_qn_\omega) +O(N_g\log N_g) +O\bigl(N_q(n_\omega+N_p)\bigr) \]

şeklinde verilmiştir. Burada N_q təbəqə projeksiyon inteqrasiya nöqtəsi sayısı, n_ω aktarım çekirdeğinin hücre desteği, N_g qüvvə ızgarası hücre sayısı və N_p eleman baz fonksiyonu sayısıdır.

Üç boyutlu kübik B-spline 4×4×4=64 hücrelik kompakt destek kullanır. Bir KDK adımı iki qüvvə dəyərlendirmesine ek olarak düyün sayısıyla doğrusal drift və kick işlemleri, indikator hesapları və yalnız işaretlenen elemanlara uygulanan yerli yenidən şəbəkə qurma maliyeti içerir.

Şəkil 12’de kütle saçma müddətsinin aktif inteqrasiya işaretleyicisi sayısıyla təxminən doğrusal, FFT Poisson müddətsinin ise N_glogN_g eğilimiyle uyumlu olduğu göstərilmişdir. Bildirilen müddətler tek düyünlü deterministik referans uygulamasına aittir; üretim səviyyəinde paralel ölçeklenme sonucu değildir.

Metodin güclü tərəfləri nelerdir?

• Soğuk qaranlıq maddəyi bşəbəkəlantısız parçacıqlar əvəzinə bşəbəkəlantılı faz-məkanı təbəqəsi olarak temsil eder.
• Kaustikleri fiziksel təbəqə qatlanmaları olarak korur və eleman bozulmasıyla karıştırmaz.
• Kütle aktarımı birliği parçalayan şəbəkəırlıklarla yuvarlama səviyyəinde korunur.
• Qüvvə, ileri hesapta kullanılan aynı ayrıq qravitasiya enerjisinden türetilir.
• Sabit qraf üzerindeki zaman inteqrasiyası kanonik və ikinci mərtəbədir.
• İleri və ters işlemler aynı ayrık operatorlerin devriklerine dayanır.
• p=3’e kadar yüksək mərtəbəli yenidən şəbəkə qurma transferleri birbaşa doğrulanmıştır.
• Adaptiv kararlar sıxlıq yüksəkliğine değil geometri, sıkışma, projeksiyon və qüvvə duyarlılığına bşəbəkəlanmıştır.
• Yeniden üretilebilirlik paketi kod, veri, yapılandırma, toxum və yoxlama cəmlərinı içerir.
• Nəticəlar qlobal norm və çökmüş bölgə normu ayrımı açıkça korunarak sunulmuştur.

Çalışmanın temel məhdudiyyətləri nelerdir?

• Hakemlik: Araştırma henüz hakem dəyərlendirmesinden geçmemiş bir preprinttir.
• Kompakt kıyaslama: Ana üç boyutlu deney 48³ qüvvə ızgarası və 64³ referansla sınırlıdır.
• Büyük hacim gösterimi yoktur: Metod gözlemsel anket ölçeğinde və ya büyük kosmoloji kutularda sınanmamıştır.
• Doğrudan ColDICE müqayisəsı yoktur: Benzer faz-təbəqəsı çözücüleriyle ortak kod səviyyəinde performans kıyaslaması yapılmamıştır.
• Qraf kararları türevlenmemektedir: Adjoint yalnız kabul edilmiş və dondurulmuş qrafın türevidir.
• Yeniden şəbəkə oluşturmada tam simplektiklik yoktur: Simplektiklik sabit qraf aralıkları üçün geçerlidir.
• Küresel sıxlıq üstünlüğü yoktur: Yeterince çözülmüş düzenli PM metodleri qlobal sıxlıq normasında daha iyidir.
• Yüksek dalga sayısı davranışı: Adaptiv metod yüksək k bölgelerinde önemli aktarım və faz sapmaları göstermiştir.
• Sınırlı toxum sayısı: Çoklu başlangıç sşəbəkəlamlığı yalnız üç toxumla dəyərlendirilmiştir.
• Performans ölçeği: Zamanlamalar üretim tipi çok düyünlü CPU/GPU ölçeklenmesini göstermemektedir.
• İndikatör sertifikalı değildir: Toplam adaptiv indikator doğrulanmış katı bir qlobal xəta üst sınırı değildir.
• Ters problem örneği küçüktür: Sekiz parametreli genlik geri kazanımı gerçek gözlemsel kozmoloji çıkarımı əvəzinə qradiyent doğrulamasıdır.

Araştırmanın desteklediği nəticəlar

• SHEET-AD’ın kütle projeksiyonu gerçekleştirilen testlerde toplam kütleyi yuvarlama səviyyəinde korumuştur.
• Varyasyonel qüvvə, seçilen ayrıq qravitasiya enerjisinin mövqea görə türeviyle sonlu fərq doğruluğunda uyuşmuştur.
• KDK güncellemesi sabit qraf testlerinde ikinci mərtəbə və simplektik davranış göstermiştir.
• Kaydedilmiş qraf üzerindeki əks rejim operatorleri ileri hesap haritasının qradiyentını doğru biçimde üretmiştir.
• Yüksek mərtəbəli təbəqə geometrisi düzgün haritalarda beklenen p+1 yaxınlaşmasına ulaşmıştır.
• İndikatorye dayalı incəltmə kompakt kıyaslamadaki çökmüş bölgə xətasını PM və sıxlıq eşiği temellerine görə azaltmıştır.

Araştırmanın kanıtlamadığı nəticəlar

• SHEET-AD bütün kosmoloji gözlenebilirlerde parçacıq-şəbəkə metodlerinden üstün değildir.
• Metodin büyük gözlem hacimlerinde daha hızlı və ya daha aşşəbəkəı bellekli olduğu gösterilmemiştir.
• Adaptiv qraf kararlarının kendisinin düzgün biçimde türevlenə bilən olduğu gösterilmemiştir.
• Yeniden şəbəkə oluşturma olayları boyunca tam simplektiklik sübut edilməmişdir.
• ColDICE və ya başka gelişmiş faz-təbəqəsı kodlarına karşı birbaşa üstünlük gösterilmemiştir.
• Yüksek k güç spektrumu və faz doğruluğunda evrensel üstünlük elde edilmemiştir.
• Gerçek galaksi gözlemleri və ya kosmoloji anket verileri üzerinde parametre çıkarımı yapılmamıştır.

Bilimsel və teknolojik anlamı nedir?

SHEET-AD’ın temel dəyəri, kosmoloji ileri modelin fiziksel və sayısal bileşenlerini tek bir operator-tutarlı yapı altında birleştirmesidir. Faz-məkanı geometrisi, kütle aktarımı, Poisson çözümü, qüvvə, zaman inteqrasiyası və əks rejim qradiyentları birbirinden müstəqil yazılım modülleri olarak değil, aynı ayrık modelin matematiksel parçaları olarak ele alınmıştır.

Bu yaklaşım xüsusiyyətle kaustik mövqeu, çok axın sayı, yüksək sıxlıq kuyrukları, çökmüş bölgəlerdeki hız alanları və ya başlangıç şərtlərına görə türev gerektiren ters problemler açısından önem taşıyabilir. Buna karşılık bütün hacmin güç spektrumu ya da qlobal sıxlıq norması gibi hedeflerde düzenli parçacıq-şəbəkə metodleri daha verimli olabilir.

Araştırma bu nedenle SHEET-AD’ı standart büyük hacimli N-cisim kodlarının əvəzinə geçen evrensel bir çözüm olarak değil; faz-məkanı sadakati, kaustik duyarlılığı, denetlenebilir xəta bütçesi və doğrulanmış sabit qraf qradiyentlarının önemli olduğu görəvleri tamamlayıcı bir araç olarak mövqelandırmaktadır.

Çalışmanın Metodi və Tapıntıları

Metodsel yapı

Doğrulama hiyerarşisi

Ana tapıntılar

• Kütle aktarımı bir, iki və üç boyutlu testlerde 10⁻¹⁶ mərtəbəsinde kütle korumuştur.
• Saçma və toplama operatorleri 4,46×10⁻¹⁶ transpozisiya qüsuruyla uyuşmuştur.
• Enerjiden türetilen qüvvə, sonlu fərq enerji qradiyentıyla 1,48×10⁻⁹ səviyyəinde uyuşmuştur.
• p=1, 2 və 3 təbəqə haritaları düzgün problemlerde təxminən ikinci, üçüncü və dördüncü mərtəbə yakınsamıştır.
• Sabit qraf KDK inteqrasiyau ikinci mərtəbə və yuvarlama səviyyəinde simplektik davranış göstermiştir.
• Sabit qraf adjoint qradiyentı merkezî sonlu fərqla 1,62×10⁻¹⁰ səviyyəinde uyuşmuştur.
• p=3 dâhil kaydedilmiş yenidən şəbəkə qurma transferleri əks rejim transpozisiya və sonlu fərq testlerini geçmiştir.
• İndikator adaptiv hesablamalar çökmüş bölgə meyaründe benzer maliyetli PM temellerinden daha aşşəbəkəı xəta vermiştir.
• Düzenli PM hesaplamaları qlobal sıxlıq normasında daha başarılı kalmıştır.
• Yüksek k spektral davranışı adaptiv metodin evrensel alan doğruluğu sşəbəkəlamadığını göstermiştir.

Şəkillerin temel mesajları

Yeniden üretilebilirlik

Tədqiqatçı, sayısal nəticəların yeniden oluşturulabilmesi üçün tek sürücülü bir iş akışı, bşəbəkəımlılık kilitleri, kapsayıcı tarifleri, deterministik təsadüfi toxumlar, yapılandırma dosyaları, doğrulama testleri, ham sayısal çıktılar, şəkil üretim betikleri və SHA-256 yoxlama cəmləri hazırlamıştır.

make reproduce komutu referans doğrulama çıktılarını, şəkilleri və yoxlama cəmlərinı yeniden üretmek üzere müəyyən edilmişdir. Ana doğrulama və adaptiv kıyaslama kampanyaları ayrı Python betikleriyle, yayın şəkilleri ise ayrı görselleştirme betiğiyle oluşturulmaktadır.

Zenodo arşivi, makaleyle əlaqəli yenidən istehsal edilə bilmə paketini veri dəsti olarak sunmaktadır. Makale DOI’si ilə veri/kod DOI’si birbirinden ayrılmalıdır:

• Makale preprint DOI’si: 10.2139/ssrn.7029053
• Yeniden üretilebilirlik paketi DOI’si: 10.5281/zenodo.20794002

Mənbə və Metod Notu

• Çalışmanın tam özgün adı: SHEET-AD: A Fixed-Graph Symplectic and Differentiable Phase-Sheet Discretization for the Cosmological Vlasov–Poisson System
• Yazar: Mushtaq Ahmad
• Yazar sırası: Tek yazarlı çalışmadır.
• Eş katkı və ya eş birinci yazarlık: Uygulanabilir değildir; eş katkı beyanı mövcud deyil.
• İletişim yazarı: Mushtaq Ahmad
• İletişim adresi: mushtaq.sial@nu.edu.pk
• Kurum 1: National University of Computer and Emerging Sciences, Islamabad, Pakistan
• Kurum 2: Center for Theoretical Physics, Khazar University, 41 Mehseti Street, Bakü AZ1096, Azerbaycan
• Kurum kaydı uyarısı: Zenodo kaydında ilk kurum “Lahore Campus” olarak gösterilirken çalışma metninde “Islamabad” yer almaktadır. Kurumsal mövqe fərqlilığı doğrulama gerektirmektedir.
• Mənbə türü: Hesaplamalı kozmoloji və sayısal metod alanında preprint araştırma makalesi
• Hakemlik vəziyyətu: Hakem dəyərlendirmesinden geçmemiştir.
• Preprint platformu: SSRN
• SSRN kayıt numarası: 7029053
• Makale DOI’si: 10.2139/ssrn.7029053
• Yıl: 2026
• Dergi: Belirli bir hakemli dergi adı doğrulanamamıştır.
• Yayınevi: Metinde “Preprint submitted to Elsevier” ifadesi bulunsa da özgün hakemli yayınevi və ya kabul edilmiş dergi sürümü doğrulanamamıştır. SSRN, preprint platformudur.
• Resmî preprint bşəbəkəlantısı:SSRN çalışma kaydı
• Yeniden üretilebilirlik paketi: Zenodo veri və kod arşivi
• Yeniden üretilebilirlik DOI’si: 10.5281/zenodo.20794002
• Resmî arşiv bşəbəkəlantısı:Zenodo yenidən istehsal edilə bilmə kaydı

Bu Verianla makalesi, araştırmanın tam metnindeki fiziksel model, sayısal ayrıklaştırma, kütle projeksiyonu, Poisson çözümü, enerji tabanlı qüvvə, simplektik zaman inteqrasiyası, adaptiv indikatorler, ayrık adjoint, yüksək mərtəbəli transferler, operator doğrulamaları və üç boyutlu müqayisəlar temel alınarak hazırlanmıştır. Araştırmada bulunmayan bir kosmoloji nəticə, büyük hacimli performans iddiası və ya gözlemsel veri başarısı eklenmemiştir.

Ana nəticə, metodin bütün sıxlıq alanlarında parçacıq-şəbəkə metodlerinden üstün olduğu biçiminde yorumlanmamalıdır. Gösterilen avantaj, kompakt kıyaslamadaki yüksək sıxlıqlu və kaustiklərə həssas bölgelere özel çəki verilmiş xəta meyarüne aittir. Ağırlıksız qlobal sıxlıq normasında yüksək çözünürlüklü düzenli PM metodleri daha aşşəbəkəı xəta vermiştir. Sabit qraf qradiyentları kabul edilmiş sayısal qrafın türevidir; adaptiv eşik kararlarının türevi değildir. Metodin büyük hacimli kosmoloji hesaplarda, gerçek gözlemsel ters problemlerde və gelişmiş faz-təbəqəsı kodlarına karşı birbaşa kıyaslamada ayrıca sınanması gerekmektedir.