Решение. Обозначим распадные константы:
λNi=t1/2(Ni)ln2,λCo=t1/2(Co)ln2. Если в начальный момент было N0 ядер никеля и не было кобальта, то
NNi(t)=N0e−λNit. Для дочернего кобальта в цепочке распадов:
NCo(t)=N0λNi−λCoλNi(e−λCot−e−λNit). Мощности нагрева:
PNi=λNiNNiENi,PCo=λCoNCoECo. Условие равенства мощностей:
λNiNNiENi=λCoNCoECo. После подстановки и сокращения N0 и λNi:
ENi=λNi−λCoλCoECo(e(λNi−λCo)t−1). Отсюда
e(λNi−λCo)t=1+λCoECoENi(λNi−λCo). Подставим числа в сутках (ln2≈0.693)
λNi=6.10.693≈0.114 суток−1,λCo=770.693≈9.00⋅10−3 суток−1. Тогда
λNi−λCo≈0.105 суток−1, 1+(9.0×10−3)⋅3.71.7⋅0.105≈1+5.36=6.36. Значит,
(λNi−λCo)t≈ln6.36≈1.85,t≈0.1051.85≈17.6 суток. Таким образом, мощность нагрева от кобальта сравнивается с никелевой примерно через 18 суток.
Теперь оценим массу никеля по светимости через 20 суток. Удобно перейти к секундам:
t=20⋅86400 с=1.728×106 с. Тогда
λNi=6.1⋅86400ln2≈1.31×10−6 с−1,λCo=77⋅86400ln2≈1.04×10−7 с−1. Получаем
NNi(20 суток)=N0e−λNit≈0.104N0, NCo(20 суток)=N0λNi−λCoλNi(e−λCot−e−λNit)≈0.821N0. Энергия одного распада в джоулях:
ENi=1.7×106⋅1.6×10−19≈2.7×10−13 Дж, ECo=3.7×106⋅1.6×10−19≈5.9×10−13 Дж. Полная мощность нагрева в момент наблюдения:
L=PNi+PCo=λNiNNiENi+λCoNCoECo. Подставляя найденные доли:
L=N0(λNi⋅0.104⋅ENi+λCo⋅0.821⋅ECo). Численно
λNi⋅0.104⋅ENi≈3.7×10−20 Вт и
λCo⋅0.821⋅ECo≈5.0×10−20 Вт на одно начальное ядро. Итак,
L≈8.7×10−20N0. При L=1.0⋅1036 Вт
N0≈8.7×10−201.0×1036≈1.15×1055. Обозначив 1 а.е.м. как u (если ее значение неизвестно, можно в качестве «замены» воспользоваться массой протона, с интересующей нас точностью результат от этого не изменится), получаем, что масса никеля
MNi=56uN0≈56⋅1.66×10−27⋅1.15×1055≈1.1×1030 кг≈0.54M⊙. Тем самым начальная масса 56Ni оказывается около 0.5 M⊙.