BeeTheory – Temeller – Teknik Not XI

Kayıp Parametrenin Belirlenmesi:
Adım 1 – Sistematik Korelasyon Analizi

Modeli değiştirmeden önce, bu not hangi gözlemlenebilir parametrenin artık hatayı en iyi tahmin ettiğini teşhis etmektedir. Not VIII’in 22 galaksilik kalibrasyon seti üzerinde çalışarak, mevcut modelin neyi ihmal ettiğini kesin olarak belirlemek için tahmin hatasının korelasyonunu fiziksel olarak anlamlı her değişkenle ve ardından her iki değişkenli kombinasyonla test ediyoruz.

1. İlk sonuç

Kayıp parametre merkezi yüzey yoğunluğudur

Merkezi baryonik yüzey yoğunluğu $\Sigma_d$ tahmin hatası ile önemsiz olmayan en güçlü korelasyona sahiptir: $r = +0,62$, $R^2 = 0,39$ tek başına.

İki değişkenli bir modelde $\Sigma_d$ ile $R_d$ disk boyutunun birleştirilmesi, tek başına $R_d$ ile $R^2 = 0,07$ olan artık varyansın $R^2 = 0,43$’ünü açıklamaktadır. RMS kalıntısı $19.5\%$’den $14.9\%$’a düşer.

Hem $R_d$ hem de $\Sigma_d$ absorbe edildikten sonra, hiçbir ek fiziksel gözlem kalıntı hakkında bilgi taşımaz.

2. Yöntem

22 galaksilik kalibrasyon seti üzerinde çalışarak (Not VIII), her galaksi için $text{err} = (V_text{tot} – V_f)/V_f$ tahmin hatasına ve ölçülebilir fiziksel parametrelerin bir listesine sahibiz. Hata ve her aday değişken arasındaki Pearson ve Spearman korelasyonlarını hesaplıyoruz, ardından formun iki değişkenli regresyonlarını test ediyoruz:

$$\text{err}(\%) \;=\; a \cdot R_d \;+\; b \cdot X \;+\; c$$

Burada $X$ her bir aday değişkendir. En iyi $X$, 22 galaksi üzerinde $R^2$ açıklanan varyansı en üst düzeye çıkaran değişkendir. Model çıktısından türetilen $V_\text{wave}$ veya $V_\text{tot}$ gibi kendi kendine referans veren değişkenler, hata ile korelasyonları totolojik olduğu için araştırmanın dışında bırakılmıştır.

3. Tek değişkenli korelasyonlar

Test edilen 24 aday değişken, hata ile mutlak Pearson korelasyonuna göre sıralanmıştır. Altın gölgeli satırlar modelin kendisinden türetilen değişkenlerdir (totolojik); kırmızı gölgeli satırlar ise $|r| > 0.5$ olan gerçek fiziksel gözlem değişkenleridir.

Değişken Açıklama Pearson $r$ p$-değeri Önem
Vw_over_VfVw / Vf oranı+0.9740.0000★★★
V_dynamicalV_dyn = √(GM_bar/Rd)+0.6320.0021★★★
log_Sigma_dlog₁₀(Σ_d)+0.6220.0026★★★
M_gasGaz kütlesi (M_sun)+0.6090.0034★★★
M_HIHI kütlesi (M_sun)+0.6090.0034★★★
THubble tipi-0.5850.0053★★
VbarBaryonik Vbar (km/s)+0.5820.0057★★
M_bar_over_Rd2M_bar / Rd²+0.5590.0084★★
VtotTahmini Vtot (km/s)+0.5550.0090★★
VwDalga Vw (km/s)+0.5500.0098★★
Vbar_over_VfVbar / Vf oranı+0.5190.0158★★
log_M_gaslog₁₀(M_gas)+0.5060.0193★★
log_M_barlog₁₀(M_bar)+0.5050.0196★★
M_barBaryonik kütle (M_sun)+0.4980.0214★★
log_M_starlog₁₀(M_star)+0.4490.0414★★
Sigma_dYüzey yoğunluğu (L/pc²)+0.4260.0544★★
M_star_over_Rd2M_star / Rd²+0.4260.0544★★
M_starYıldız kütlesi (M_sun)+0.3890.0815

Tablonun okunması

En yüksek tek korelasyon $V_\text{wave}/V_f = +0.974$’tür. Bu totolojiktir: yapı gereği, hata doğrudan $V_\text{wave}$ ile ölçeklenir, bu nedenle bu değişken harici bir fiziksel sürücüyü değil, sadece tahmin formülünün yapısını yansıtır.

Gerçek fiziksel gözlemlenebilirler arasında en yüksek korelasyonlar $\log(\Sigma_d) = +0.622$, $V_\text{dynamical} = +0.632$, $M_\text{gas} = +0.609$ ve Hubble tipi $T = -0.585$’dir. Bu dört sinyal fiziksel olarak bağlantılıdır: yoğun diskler daha büyük kütleli, daha erken tipte ve daha yüksek baryonik dinamik hıza sahip olma eğilimindedir. Asıl soru, hangisinin temel etken olduğudur.

4. Gereksiz değişkenlerin filtrelenmesi

En yüksek korelasyonlu değişkenlerin birçoğu, hatayı yönlendirdiği bilinen $R_d$ değişkeni ile güçlü bir korelasyona sahiptir. Asıl soru, hangisinin bağımsız bilgi taşıdığıdır.

Değişken R_d$ ile korelasyon Durum
$\log(M_\star)$$r = +0.88$R_d$ ile gereksiz
$\log(M_\text{bar})$$r = +0.87$R_d$ ile gereksiz
$\log(M_\text{gas})$$r = +0.86$R_d$ ile gereksiz
Hubble tipi $T$$r = -0.66$Kısmen gereksiz
$V_\text{dynamical}$$r = +0.50$Kısmen bağımsız
$M_\text{bar}/R_d^2$$r = -0.19$Bağımsız
$\log(\Sigma_d)$$r = +0.10$Bağımsız

Kütleler $R_d$ ile neredeyse mükemmel bir korelasyon gösterir: daha büyük bir disk basitçe daha fazla baryonik malzeme içerir. Dolayısıyla bu değişkenler esasen $R_d$’nin kendisiyle aynı bilgiyi taşır. Buna karşılık, $\Sigma_d$ (merkezi yüzey yoğunluğu) ve $M_\text{bar}/R_d^2$ (ortalama baryonik yüzey yoğunluğu) bu örnekte $R_d$ ile neredeyse ortogonaldir: “diskin ne kadar genişlediğinden” bağımsız olarak “maddenin ne kadar kompakt olduğu” yapısal özelliğini yakalarlar.

5. Yüzey yoğunluğuna karşı hata – görselleştirme

Hatanın sadece $\log_{10}(\Sigma_d)$ ‘ye karşı grafiği, Hubble tipine göre renklendirilmiştir:

Merkezi yüzey yoğunluğuna karşı tahmin hatası $Sigma_d$ – 22 kalibrasyon galaksisi 101001000 -50%-25%+0%+25%+50% Merkezi yüzey yoğunluğu $\Sigma_d$ ($L_\odot/$pc$^2$) – log ölçeği Tahmin hatası (%) err = 33,0-log(Σd) – 60Pearson r = 0,622, R² = 0,387 D631-7 Σd=115 err=-12%DDO064 Σd=120 err=+13%DDO154 Σd=45 err=+%4DDO161 Σd=35 err=+%11DDO168 Σd=100 err=-21%DDO170 Σd=25 err=+%15ESO116-G012 Σd=115 err=+14%ESO444-G084 Σd=60 err=+%6F561-1 Σd=30 err=-%15F563-1 Σd=20 err=-18%F563-V1 Σd=25 err=-36%F563-V2 Σd=30 err=-27%F565-V2 Σd=18 err=-39%F567-2 Σd=15 err=-23%F568-1 Σd=40 err=-%8F568-3 Σd=35 err=-12%F568-V1 Σd=20 err=-24%F571-8 Σd=50 err=+14%F574-1 Σd=30 err=-2%NGC2841 Σd=605 err=+18%NGC3198 Σd=153 err=+43% S0-SaSb-SbcSc-ScdSd-Im
Tek değişkenli uyum hatası$(\%) = 33\log(\Sigma_d) – 60$, Pearson $r = 0,62$, $R^2 = 0,39$.

Eğilim açık ve monotoniktir: daha yüksek merkezi yüzey yoğunluğuna sahip galaksiler BeeTheory tarafından sistematik olarak fazla tahmin edilirken, dağınık düşük yoğunluklu diskler az tahmin edilmektedir. On yıllık $\Sigma_d$ başına $33$ yüzde puanlık uyum eğimi, 15 ila 605 $L_\odot/\text{pc}^2$ aralığındaki tüm verilerle sağlam bir şekilde eşleşmektedir.

6. İki değişkenli modeller – karşılaştırma

Her bir aday değişkene $R_d$ eklendiğinde daha net bir sıralama elde edilir. Aşağıdaki tablo, $R_d$ her bir ikinci değişkenle eşleştirildiğinde açıklanan $R^2$ varyansını göstermektedir (totolojik kombinasyonlar hariç):

İki değişkenli model $R^2$ RMS kalıntı Notlar
$\text{err} = a R_d + c$ (tek değişkenli taban çizgisi)0.074$19.5\%$Referans, ikinci değişken yok
$\text{err} = a R_d + b f_\text{gas} + c$0.101$19.3\%$Önemsiz iyileşme
$\text{err} = a R_d + b \log M_\star + c$0.272$17.3\%$
$\text{err} = a R_d + b V_\text{bar} + c$0.345$16.4\%$
$\text{err} = a R_d + b \log M_\text{gas} + c$0.359$16.3\%$
$\text{err} = a R_d + b T + c$0.367$16.2\%$
$\text{err} = a R_d + b \log M_\text{bar} + c$0.373$16.1\%$
$\text{err} = a R_d + b\,V_\text{dynamical} + c$0.402$15.7\%$Güçlü
$\text{err} = a R_d + b \log\Sigma_d + c$0.430$15.3\%$R_d$’den bağımsız
$\text{err} = a R_d + b (M_\text{bar}/R_d^2) + c$0.459$14.9\%$En iyi totolojik olmayan model

En iyi iki değişkenli model

$$\text{err}(\%) \;=\; a\,R_d \;+\; b\,\frac{M_\text{bar}}{R_d^2} \;+\; c, \qquad R^2 = 0,46$$

M_\text{bar}/R_d^2$ değişkeni diskin ortalama baryonik yüzey yoğunluğudur, $\langle \Sigma_\text{bar} \rangle = M_\text{bar}/(\pi R_d^2)$. Diskin ne kadar büyük olduğundan bağımsız olarak, görünür maddenin ne kadar kompakt olduğu hakkında bilgi taşır. Bu, BeeTheory’nin şu anda hesaba katamadığı değişkendir.

7. Kapanış kontrolü – $R_d$ ve $\Sigma_d$ hesaplandıktan sonra geriye ne kalır

Eğer $R_d$ ve $\log \Sigma_d$ birlikte yapısal kusuru yakalıyorsa, iki değişkenli uyumun artığı her fiziksel gözlemlenebilir ile ilişkisiz olmalıdır. Bunun test edilmesi resmi kapanış kontrolüdür:

Değişken Kalıntı ile korelasyon Durum
$R_d$$+0.00$İnşaat tarafından
$\log \Sigma_d$$+0.00$İnşaat tarafından
$\log M_\star$$-0.05$Emilmiş
$\log M_\text{bar}$$+0.07$Emilmiş
$\log M_\text{gas}$$+0.14$Emilmiş
Hubble tipi $T$$-0.04$Emilmiş
$V_\text{dynamical}$$+0.08$Emilmiş
$V_\text{bar}$$+0.05$Emilmiş
$f_\text{gas}$$+0.28$Marjinal; anlamlılığın altında

R_d$ ve $\log \Sigma_d$ değişkenleri hesaba katıldıktan sonra, hiçbir fiziksel gözlem kalıntısı artık hata ile önemli bir korelasyon göstermez. Hatadaki yapısal bilgi bu iki değişken tarafından tamamen yakalanmıştır. Kalan %15$ RMS dağılımı, SPARC girdi parametrelerindeki gözlemsel belirsizlik ve bu toplam tanımlayıcıların hiçbiri tarafından yakalanamayan galaksiden galaksiye içsel değişkenlik ile tutarlıdır.

8. Fiziksel yorumlama

Mevcut BeeTheory modeli disk ölçek uzunluğunu $R_d$ iki yerde kullanır: baryonik dağılımın uzamsal ölçeği olarak (üstel profil $Sigma propto e^{-R/R_d}$) ve dalga çekirdeğinin tutarlılık uzunluğu olarak ($ell = c_text{disk},R_d$). Baryonik profilin genliği $\Sigma_0$ örtüktür ve entegre edildiğinde doğru yıldız kütlesini verecek şekilde ölçeklendirilir.

Yüzey yoğunluğu fiziksel olarak neyi temsil eder?

Ortalama baryonik yüzey yoğunluğu $langle Sigma_text{bar} rangle = M_text{bar}/(pi R_d^2)$ diskin birim alanı başına düşen kütledir. Aynı $R_d$ ancak farklı $\Sigma_d$ değerlerine sahip iki galaksi aynı geometrik genişliğe ancak farklı miktarlarda maddeye sahiptir. Mevcut model sadece geometrik kapsamı ($R_d$) dalga tutarlılığı uzunluğu ile ilgili olarak ele alır ve maddenin ne kadar yoğunlaştığını göz ardı eder. Bu tam da kalıntı analizinin eksik olarak tanımladığı parametredir.

Etkinin yönü

Korelasyon pozitiftir: hata yüzey yoğunluğu ile artar. Bu, sabit $R_d$ için, daha yoğun disklerin model tarafından aşırı tahmin edildiği anlamına gelir – dalga alanı dönme eğrisine göre çok güçlüdür. Tersine, belirli bir $R_d$ için, model dağınık düşük yoğunluklu diskleri eksik tahmin eder. Makul bir fiziksel yorum: dalga tutarlılık uzunluğu sadece kaynağın geometrik boyutuna değil, aynı zamanda yoğunluğuna da bağlı olmalıdır, daha yoğun madde daha yerel bir dalga tepkisi üretir. Bu durum doğal olarak yüksek$\Sigma$ disklerde dalga alanı genliğini bastırırken, düşük$\Sigma$ olanlarda genliği arttıracaktır.

9. Adım 1’in Özeti

1. 22 galaksilik kalibrasyon setinde tahmin hatası, gerçek fiziksel gözlemlenebilirler arasında en güçlü şekilde merkezi yüzey yoğunluğu $\Sigma_d$ ($r = +0.62$) ile ilişkilidir.

2. Başlangıçta güçlü korelasyon gösteren diğer değişkenlerin (yıldız kütlesi, gaz kütlesi, baryonik kütle) $R_d$ ile oldukça gereksiz olduğu ($R_d$ ile $\geq 0.86$ korelasyon) ve bu nedenle çok az yeni bilgi taşıdığı ortaya çıkmıştır.

3. En iyi totolojik olmayan iki değişkenli model $\text{err} = a\,R_d + b\,(M_\text{bar}/R_d^2) + c$, $R^2 = 0.46$ ve RMS kalıntısı $14.9\%$’dir. İkinci değişken diskin ortalama baryonik yüzey yoğunluğudur.

4. R_d$ ve $\Sigma_d$ hesaba katıldıktan sonra, başka hiçbir gözlemlenebilir kalıntı ile önemli bir korelasyona sahip değildir. Teşhis kapanmıştır.

5. Eksik parametre belirlenmiştir: mevcut BeeTheory modeli baryonik dağılımın geometrik boyutunu ($R_d$) hesaba katmakta ancak yüzey yoğunluğunu ($\Sigma_d$) hesaba katmamaktadır. Bir sonraki adım, $\Sigma_d$’yi dalga tutarlılığı uzunluğuna ikinci bir girdi olarak dahil etmek ve ardından modeli 22 galaksi seti üzerinde yeniden kurmaktır.


Referanslar. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: Spitzer Fotometrisi ve Doğru Dönme Eğrileri ile 175 Disk Galaksisi için Kütle Modelleri, AJ 152, 157 (2016). – Pearson, K. – Mathematical contributions to the theory of evolution III, Phil. Trans. R. Soc. A 187, 253 (1896). Korelasyon katsayısı. – Dutertre, X. – Bee Theory™: Yerçekiminin Dalga Tabanlı Modellemesi, v2, BeeTheory.com (2023).

BeeTheory.com – Dalga tabanlı kuantum yerçekimi – Teşhis adımı 1 – © Technoplane S.A.S. 2026