Сформулируйте способ оценить массу галактического скопления по наблюдаемому гравитационному линзированию «кользящим» образом, включая учет распределения тёмной материи и неоднородностей
Коротко — пошаговый «кользящий» (скользящий апертурный + итеративный) алгоритм для оценки массы скопления по гравитационному линзированию с учётом тёмной материи и неоднородностей: 1) Подготовка данных: - получить оценки компонент сдвига галактик \(\gamma=(\gamma_1,\gamma_2)\), их координаты и распределение/фото-z \(n(z)\); исправить PSF, калибровку формы, отсеять звёзды. 2) Критическая поверхностная плотность: - для каждой пары линз–источник вычислить \(\Sigma_{\rm crit}=\frac{c^2}{4\pi G}\frac{D_s}{D_d D_{ds}}\), где \(D_d,D_s,D_{ds}\) — угловые диаметральные расстояния. 3) Скользящее (апертурное) картирование массы: - в каждой позиции центра апертуры \(\theta_0\) вычислять апертурную массу через тангенциальную компоненту сдвига: \[ M_{\rm ap}(\theta_0)=\int d^2\theta\;Q(|\theta-\theta_0|)\;\gamma_t(\theta;\theta_0), \] где \(Q\) — компенсированный фильтр (например, полиномиальный или Шримера) с радиусом апертуры \(R\). Путём скольжения апертуры по полю и варьирования \(R\) строить многошкальную карту плотности и находить пики (кандидаты-гало). 4) Непараметрическая реконструкция (для глобальной карты): - выполнить обратную задачу Kaiser–Squires (в частотной форме): \[ \widehat\kappa(\boldsymbol\ell)=\frac{\ell_1^2-\ell_2^2 + 2i\ell_1\ell_2}{\ell^2}\,\widehat\gamma(\boldsymbol\ell), \] и получить \(\kappa(\theta)\). Гладить/оставлять многомасштабные компоненты (wavelet) для выделения субструктур. 5) Параметрическая подгонка пиков (учёт DM-профиля и субгало): - модель основной массы как NFW: \[ \rho(r)=\frac{\rho_s}{(r/r_s)(1+r/r_s)^2}, \] проектированная поверхностная плотность \(\Sigma_{\rm NFW}(R)\) даёт предсказанную \(\kappa(R)=\Sigma_{\rm NFW}(R)/\Sigma_{\rm crit}\) и сдвиг \(\gamma^{\rm model}(R)\). Для каждого пика подгонять параметры \((M_{200},c,\text{центр})\); включать дополнительные маломасштабные гало для субструктуры. 6) Итеративное «скольжение/очищение» (peeling): - для обнаруженного пика: подогнать парам. гало, из предсказанного \(\gamma^{\rm model}\) вычесть вклад из наблюдаемого поля, затем повторить поиск пиков на остатке (чтобы учесть перекрытия гало и разложить неоднородности). 7) Объединение сильного и слабого линзирования, томография: - использовать сильные линзовые конфигурации для точной локальной нормировки (ломает mass-sheet), и многозонную (tomographic) информацию по \(n(z)\) для ограничения распределения вдоль луча. 8) Байесовская оценка и учёт шумов: - строить правдоподобие по компонентам сдвига: \[ \chi^2=(\gamma^{\rm obs}-\gamma^{\rm model})^T C^{-1}(\gamma^{\rm obs}-\gamma^{\rm model}), \] и получать постериоры по MCMC/nested-sampling для параметров гало. В ковариацию \(C\) включить шум форм (shape noise), ковариацию фонов из крупномасштабной структуры (LSS), погрешности по z и систематики калибровки. 9) Масса в физическом радиусе: - итоговую массу внутри углового радиуса \(\theta_R\) переводят в физический радиус \(R=D_d\theta_R\) и вычисляют \[ M(<R)=\Sigma_{\rm crit}\int_{|\theta|<\theta_R} d^2\theta\;\kappa(\theta) =2\pi\int_0^{R}\Sigma(R')\,R'\,dR'. \] 10) Контрольные проверки и систематика: - оценить биас/разброс через имитации и мок-данные (включая реалистичные субструктуры и LSS); - проверять влияние mass-sheet degeneracy (решать через magnification: измерения числа/яркости/размеров источников) и чувствительность к выбору фильтра/масштаба апертуры. Краткое резюме: сочетание скользящего апертурного картирования \(M_{\rm ap}\) (для обнаружения и локальной оценки), глобальной непараметрической реконструкции (\(\kappa\) через KS) и параметрической подгонки NFW (с итеративным вычитанием субгало), плюс томография и объединение со strong-lensing/magnification, даёт надёжную оценку массы скопления с учётом тёмной материи и неоднородностей.
1) Подготовка данных:
- получить оценки компонент сдвига галактик \(\gamma=(\gamma_1,\gamma_2)\), их координаты и распределение/фото-z \(n(z)\); исправить PSF, калибровку формы, отсеять звёзды.
2) Критическая поверхностная плотность:
- для каждой пары линз–источник вычислить \(\Sigma_{\rm crit}=\frac{c^2}{4\pi G}\frac{D_s}{D_d D_{ds}}\), где \(D_d,D_s,D_{ds}\) — угловые диаметральные расстояния.
3) Скользящее (апертурное) картирование массы:
- в каждой позиции центра апертуры \(\theta_0\) вычислять апертурную массу через тангенциальную компоненту сдвига:
\[
M_{\rm ap}(\theta_0)=\int d^2\theta\;Q(|\theta-\theta_0|)\;\gamma_t(\theta;\theta_0),
\]
где \(Q\) — компенсированный фильтр (например, полиномиальный или Шримера) с радиусом апертуры \(R\). Путём скольжения апертуры по полю и варьирования \(R\) строить многошкальную карту плотности и находить пики (кандидаты-гало).
4) Непараметрическая реконструкция (для глобальной карты):
- выполнить обратную задачу Kaiser–Squires (в частотной форме):
\[
\widehat\kappa(\boldsymbol\ell)=\frac{\ell_1^2-\ell_2^2 + 2i\ell_1\ell_2}{\ell^2}\,\widehat\gamma(\boldsymbol\ell),
\]
и получить \(\kappa(\theta)\). Гладить/оставлять многомасштабные компоненты (wavelet) для выделения субструктур.
5) Параметрическая подгонка пиков (учёт DM-профиля и субгало):
- модель основной массы как NFW:
\[
\rho(r)=\frac{\rho_s}{(r/r_s)(1+r/r_s)^2},
\]
проектированная поверхностная плотность \(\Sigma_{\rm NFW}(R)\) даёт предсказанную \(\kappa(R)=\Sigma_{\rm NFW}(R)/\Sigma_{\rm crit}\) и сдвиг \(\gamma^{\rm model}(R)\). Для каждого пика подгонять параметры \((M_{200},c,\text{центр})\); включать дополнительные маломасштабные гало для субструктуры.
6) Итеративное «скольжение/очищение» (peeling):
- для обнаруженного пика: подогнать парам. гало, из предсказанного \(\gamma^{\rm model}\) вычесть вклад из наблюдаемого поля, затем повторить поиск пиков на остатке (чтобы учесть перекрытия гало и разложить неоднородности).
7) Объединение сильного и слабого линзирования, томография:
- использовать сильные линзовые конфигурации для точной локальной нормировки (ломает mass-sheet), и многозонную (tomographic) информацию по \(n(z)\) для ограничения распределения вдоль луча.
8) Байесовская оценка и учёт шумов:
- строить правдоподобие по компонентам сдвига:
\[
\chi^2=(\gamma^{\rm obs}-\gamma^{\rm model})^T C^{-1}(\gamma^{\rm obs}-\gamma^{\rm model}),
\]
и получать постериоры по MCMC/nested-sampling для параметров гало. В ковариацию \(C\) включить шум форм (shape noise), ковариацию фонов из крупномасштабной структуры (LSS), погрешности по z и систематики калибровки.
9) Масса в физическом радиусе:
- итоговую массу внутри углового радиуса \(\theta_R\) переводят в физический радиус \(R=D_d\theta_R\) и вычисляют
\[
M(<R)=\Sigma_{\rm crit}\int_{|\theta|<\theta_R} d^2\theta\;\kappa(\theta)
=2\pi\int_0^{R}\Sigma(R')\,R'\,dR'.
\]
10) Контрольные проверки и систематика:
- оценить биас/разброс через имитации и мок-данные (включая реалистичные субструктуры и LSS);
- проверять влияние mass-sheet degeneracy (решать через magnification: измерения числа/яркости/размеров источников) и чувствительность к выбору фильтра/масштаба апертуры.
Краткое резюме: сочетание скользящего апертурного картирования \(M_{\rm ap}\) (для обнаружения и локальной оценки), глобальной непараметрической реконструкции (\(\kappa\) через KS) и параметрической подгонки NFW (с итеративным вычитанием субгало), плюс томография и объединение со strong-lensing/magnification, даёт надёжную оценку массы скопления с учётом тёмной материи и неоднородностей.