第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 9048374

2 8187308

3 7408182

……

7 4965853

8 449329

……

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

很快，他便锁定了其中的第十八行：

18 1652989。

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

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

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。

见此情形。

徐云瞳孔顿时微微一缩。

这个结果的意思就是……

在一开始，y(xn＋1)－y(xn)/h≈f这个轨道上便存在有一颗粒子。

只是在撞击过程中它寿命终止或者跃迁失能了，所以最终没有被捕捉到。

想到这里。

徐云沉默片刻，走出图书馆。

拿出手机拨通了一个号码。

片刻过后。

手机接通，某个一听就知道很帅的声音从对头传了过来：

“喂，小徐？”

“嗯，是我，老师您这会儿有空吗？”

“刚出实验室，啥事儿？”

徐云组织了一番语言，说道：

“老师，我之前不是研究过一个Σ超子的课题吗？您还记得不？”

Σ超子是目前比较主流