但我们已经提前知道了它的运动轨迹,那么完全可以事先就在那儿放一块干净的采样板。
然后双手离开现场,找个椅子做好,安静等它送上门来就行。
眼下有了Λ超子的信息,还有了公式模型,推导“落点”的环节也就非常简单了。
众所周知。
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。
见此情形。