ВсОА 2023, Сбор 1 · анализ на данни — Задача 2. D12 Чёрный юмор

Q23S1D

Условие

В файле zT.dat даны значения температуры реликтового излучения и её неопределенности в зависимости от величины красного смещения (Gelo et al., 2022). Если предположить, что Вселенная расширяется адиабатически, количество фотонов сохраняется, а спектр реликтового излучения изначально чернотельный, то температура реликтового излучения будет зависеть от красного смещения как T(z)=T0(1+z)T(z) = T_0(1+z), где современное значение температуры фона T0=2.7255±0.0006T_0 = 2.7255 \pm 0.0006 K. Тем не менее, наблюдательные данные не идеально ложатся на такую модель. Gelo et al. (2022) предложили зависимость следующего вида:

T(z)T0=(1+z)3wr((m−3Ωm)+m(Ωm−1)(1+z)m−3(m−3)Ωm)wr.\frac{T(z)}{T_0} = (1+z)^{3w_r}\left(\frac{(m-3\Omega_m)+m(\Omega_m-1)(1+z)^{m-3}}{(m-3)\Omega_m}\right)^{w_r}.

Здесь wr,m,Ωmw_r, m, \Omega_m — параметры модели. При использовании стандартного уравнения состояния для излучения wr≡1/3w_r \equiv 1/3.

a. Минимизируя сумму квадратов невязок S2S^2, определите параметры mm и Ωm\Omega_m модели по имеющимся данным: постройте тепловую карту значений S2S^2 на плоскости параметров (m,Ωm)(m, \Omega_m) в диапазоне m∈[0.03:0.40]m \in [0.03:0.40], Ωm∈[0.20:0.40]\Omega_m \in [0.20:0.40] с шагом 0.01 по каждой переменной, определите точку с наименьшей S2S^2. Можно ли считать, что эта точка соответствует глобальному минимуму S2S^2?

b. Постройте зависимость T(z)T(z), проведите на графике оба варианта теоретической зависимости. Оцените S2S^2 для предлагаемой простой зависимости.

Примечание. В случае учёта неопределенностей параметров рассматривается взвешенная сумма квадратов невязок:

S2=∑i=1N1σi2(Ti−f(zi))2,S^2 = \sum_{i=1}^{N}\frac{1}{\sigma_i^2}(T_i-f(z_i))^2,

где σi\sigma_i — неопределённость измерения температуры TiT_i, f(Ti)f(T_i) — модельная зависимость температуры от красного смещения.

Решение

Внимание: Непълно решение

В архива няма официално решение на тази задача.


Оригинал в Архива: Q23S1D.pdf