第276节(1 / 2)

加入书签

但我们已经提前知道了它的运动轨迹,那么完全可以事先就在那儿放一块干净的采样板。

然后双手离开现场,找个椅子做好,安静等它送上门来就行。

眼下有了Λ超子的信息,还有了公式模型,推导“落点”的环节也就非常简单了。

众所周知。

n及衰变的通解并不复杂。

比如存在衰变链a→b→c→d……,各种核素的衰变常数对应分别为λ1、λ2、λ3、λ4……

假设初始t0时刻只有a,则显然:n1=n1(0)exp(-λ1t)。

随后徐云又写下了另一个方程:

dn2/dt=λ1n1-λ2n2。

这是b原子核数的变化微分方程。

求解可得n2=λ1n1(0)[exp(-λ1t)-exp(-λ2t)]/(λ2-λ1)。

随后徐云边写边念:

“c原子核的变化微分方程是:dn3/dt=λ2n2-λ3n3,即dn3/dt+λ3n3=λ2n2……”

“代入上面的n2,所以就是n3=λ1λ2n1(0){exp(-λ1t)/[(λ2-λ1)(λ3-λ1)+exp(-λ2t)/[(λ1-λ2)(λ3-λ2)]+exp(-λ3t)/[(λ1-λ3)(λ2-λ3)]}……”

写完这些他顿了顿,简单验算了一遍。

确定没有问题后,继续写道:

“可以定义一个参数h,使得h1=λ1λ2/[(λ2-λ1)(λ3-λ1)],h2=λ1λ2/[(λ1-λ2)(λ3-λ2)],h3=λ1λ2/[(λ1-λ3)(λ2-λ3)]……”

“则n3可简作:n3=n1(0)[h1exp(-λ1t)+h2exp(-λ2t)+h3exp(-λ3t)]。”

写完这些。

徐云再次看向屏幕,将Λ超子的参数代入了进去:

“n=n1(0)[h1exp(-λ1t)+h2exp(-λ2t)+……hnexp(-λnt)],h的分子就是Πλi,i=1~n-1,即分子是λ1λ2λ3λ4……”

“Λ超子的衰变周期是17,所以h1的分母,就是除开Λ超子前一种衰变常数与Λ超子衰变常数λ1的差的积……”

半个小时后。

极光软件上现实出了一组数值。

a a 0 1000:

1 904.8374

2 818.7308

3 740.8182

……

7 496.58538 449.329

……

徐云没去看前面的数字,飞快的将鼠标下拉。

很快,他便锁定了其中的第十八行:

18 165.2989。

有了这一组数字,接下来的问题就非常简单了。

徐云将这种数字输入了极光模型,公式为:

f(t):=n(t)/n(0)=e^(-t/π)。

这里的“:=”是定义符号,它表示将右边的东西定义成左边的东西。

徐云现在为这个f(t)赋予了一个物理意义:

某个原子在时刻t依然存活(没有衰变)的概率。

n=n1(0)[h1exp(-λ1t)+h2exp(-λ2t)+……hnexp(-λnt)]这个公式描述了到时刻t还剩多少原子,徐云所作的是将剩下的原子数目比上最初的总原子数,这个量自然就是在那堆剩下的原子中能找到徐云想要的那个的概率。

非常简单,也非常好理解。

极光系统连接的是中科院的次级服务器,使用的是中科院超算“夜语”的部分算力。

因此只过了十多分钟。

他面前的屏幕上便显示出了一个结果:

t=0,f=1。

见此情形。

↑返回顶部↑

书页/目录