Практика 1. Среда, единицы и моменты

Что считаем

Создаём lighthit.Medium для 450 нм и проверяем \(v_g=c/n_g\), длины поглощения, рассеяния и переноса. Затем применяем точные формулы лекции 1 к направленной вспышке. Пространственные моменты нормированы на выжившие фотоны.

import lighthit as lh
medium = lh.Medium(0.076, 0.019, 0.9, 1.3713237239062512)
print(medium.speed_m_per_ns)

Запуск

После установки и распаковки примеров:

python 01_medium.py

Полный скрипт · Общая среда.

Проверяемые величины

Величина	Значение на 450 нм
\(\ell_a\)	13.16 м
\(\ell_s\)	52.63 м
\(\ell_*\) при \(g=0.9\)	526.32 м
\(\mu_a/\mu_s'\)	40

summary.json содержит коэффициенты, скорость и контрольные числа. moments.csv содержит время, вероятность выживания и координату центра облака. В centre.svg центр выражен в транспортных единицах.

По лекции \(\langle z\rangle=\ell_*(1-e^{-v_gt/\ell_*})\). Формула описывает центр выжившего облака; число выживших при этом убывает как \(e^{-\mu_av_gt}\).

Самостоятельная проверка

Измените только \(g\). Какие длины изменятся? Измените только \(\mu_a\). Изменится ли нормированный центр при фиксированном времени?