发明名称 飞秒激光团簇尺寸的标定方法
摘要 一种用于飞秒激光团簇相互作用中飞秒激光团簇尺寸的标定方法,包括下列步骤:根据激光团簇相互作用实验中获得的离子飞行时间谱推导实验上的离子能谱;模拟团簇尺寸的对数正态分布;模拟团簇半径的分布;模拟团簇的离子能谱;确定最终的团簇半径分布和团簇尺寸分布;完成了团簇尺寸的标定。实验表明,本发明可以探测团簇的尺寸范围很大,能够对团簇平均尺寸实现绝对标定,并且可以得到团簇尺寸分布的准确信息,是一种非常有效的诊断团簇尺寸的方法。
申请公布号 CN103234491A 申请公布日期 2013.08.07
申请号 CN201310145187.X 申请日期 2013.04.24
申请人 中国科学院上海光学精密机械研究所 发明人 张辉;卢海洋;李松;李儒新
分类号 G01B15/00(2006.01)I 主分类号 G01B15/00(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯
主权项 1.一种飞秒激光团簇尺寸的标定方法,其特征在于该方法包括下列步骤:①根据激光团簇相互作用实验中获得的离子飞行时间谱fto(m,j)推导实验上的离子能谱fEo(m,j),包括下列步骤:&lt;1&gt;所述的离子飞行时间谱和离子能谱分别用矩阵fto(m,j)和fEo(m,j)表述,其中,m=1,2,…,M;j=1,2,矩阵fto(m,j)的第一列fto(m,1)表示离子飞行时间,矩阵fto(m,j)的第二列fto(m,2)表示对应的飞行时间分布,离子能谱矩阵fEo(m,j)的第一列fEo(m,1)表示离子能量数据,第二列fEo(m,2)表示对应的离子能量分布;&lt;2&gt;团簇分子组成为A<sub>x</sub>B<sub>y</sub>,A为重核原子,B为轻核原子,根据库伦爆炸原理得知飞行时间谱记录的是轻核离子的信号,元素B的质量为m<sub>B</sub>,飞行距离为L,飞行的时间为fto(m,1),则B离子的能量为:<maths num="0001"><![CDATA[<math><mrow><mi>fEo</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>m</mi><mi>B</mi></msub><mn>2</mn></mfrac><mo>*</mo><msup><mrow><mo>(</mo><mfrac><mi>L</mi><mrow><mi>fto</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>对应B离子的能量分布为:fEo(m,2)=-fto(m,2)*fto(m,1)<sup>3</sup>;②模拟团簇尺寸的对数正态分布fN(n,j),包括下列步骤;&lt;1&gt;团簇尺寸是对数正态分布的,设置团簇尺寸的最大值为Nmax,团簇尺寸峰值为Npeak和团簇对数尺寸的标准方差为σ,团簇对数尺寸的平均值μ满足:μ=ln(Npeak)+σ<sup>2</sup>;&lt;2&gt;所述的团簇尺寸的对数正态分布用矩阵fN(n,j)表述,其中n=1,2,…,Nmax,j=1,2,该矩阵的第一列fN(n,1)表示团簇尺寸,即团簇内分子的数目,第二列fN(n,2)表示对应的团簇尺寸分布,团簇尺寸为:fN(n,1)=n,对应的团簇尺寸的对数正态分布为:<maths num="0002"><![CDATA[<math><mrow><mi>fN</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt><mi>&sigma;</mi><mo>*</mo><mi>fN</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>ln</mi><mi>fN</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mn>2</mn></msup></mfrac><mo>]</mo><mo>;</mo></mrow></math>]]></maths>③模拟团簇半径的分布fR(i,j),包括下列步骤;&lt;1&gt;所述的团簇半径分布用矩阵fR(i,j)表述,其中i=1,2,…,I,j=1,2,该矩阵的第一列fR(i,1)表示团簇的半径,第二列fR(i,2)表示团簇的半径分布,团簇半径为:<img file="FDA00003095811900013.GIF" wi="439" he="124" />其中Rmax是团簇的最大半径;&lt;2&gt;建立一个长度为I的一维零矩阵N(i),对N(i)赋值为:<maths num="0003"><![CDATA[<math><mrow><mi>N</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mi>round</mi><mo>[</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mi>&pi;&rho;</mi><mo>*</mo><mi>fR</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mn>3</mn></msup><mo>]</mo><mo>,</mo></mrow></math>]]></maths>其中,ρ是团簇内团簇分子的平均密度,round表示数值取整,所述的N(i)即为团簇半径fR(i,1)对应的团簇尺寸;&lt;3&gt;根据团簇尺寸的对数正态分布fN(n,2)计算团簇半径的分布fR(i,2):对fN(n,1)进行扫描,当N(i-1)≤fN(n,1)≤N(i)时,团簇尺寸fN(n,1)对应的团簇尺寸分布fN(n,2)对团簇半径分布fR(i,2)有贡献,即团簇半径分布fR(i,2)是对团簇尺寸分布fN(n,2)在区间N(i-1)≤n≤N(i)求和:<maths num="0004"><![CDATA[<math><mrow><mi>fR</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mi>N</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>N</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></munderover><mi>fN</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>当i=1时,团簇半径分布fR(1,2)是对fN(n,2)在区间1≤n≤N(1)求和:<maths num="0005"><![CDATA[<math><mrow><mi>fR</mi><mrow><mo>(</mo><mn>1,2</mn><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></munderover><mi>fN</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>④模拟团簇的离子能谱fE(i,j),包括下列步骤;&lt;1&gt;所述的模拟的团簇离子能谱用矩阵fE(i,j)表述,其中i=1,2,…,I,j=1,2,该矩阵的第一列fE(i,1)是离子能量,矩阵的第二列fE(i,2)是对应的离子能量的分布;&lt;2&gt;团簇分子A<sub>x</sub>B<sub>y</sub>中,元素A的平均电离价态为+q<sub>A</sub>,元素B的平均电离价态为+q<sub>B</sub>,库伦爆炸时单个团簇分子产生的电荷数是q=(x*q<sub>A</sub>+y*q<sub>B</sub>)e,其中e是电子的电荷量;根据库伦爆炸理论,由团簇的半径fR(i,1)获得团簇爆炸后的B离子能量fE(i,1)为:<maths num="0006"><![CDATA[<math><mrow><mi>fE</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mrow><mn>3</mn><mi>&epsiv;</mi></mrow><mn>0</mn></msub></mfrac><msub><mi>&rho;qq</mi><mi>B</mi></msub><mi>e</mi><mo>*</mo><mi>fR</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>其中ε<sub>0</sub>是真空介电常数,ρ是团簇内团簇分子的平均密度;&lt;3&gt;对于给定的团簇半径fR(i,1),当fR(i′,1)≥fR(i,1)时,团簇半径的分布fR(i′,2)对离子能量分布fE(i,2)都有贡献,根据团簇的离子能量fE(i,1)和团簇半径fR(i,1)的关系得到,团簇半径的分布fR(i′,2)对离子能量分布fE(i,2)的贡献为:fR(i′,2)*fR(i,1)/fR(i′,1)<sup>3</sup>,即离子能量分布fE(i,2)是对所述的团簇半径分布fR(i′,2)在区间i≤i′≤I上求和:<maths num="0007"><![CDATA[<math><mrow><mi>fE</mi><mrow><mo>(</mo><mn>1,2</mn><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><msup><mi>i</mi><mo>&prime;</mo></msup><mo>=</mo><mi>i</mi></mrow><mi>I</mi></munderover><mi>fR</mi><mrow><mo>(</mo><msup><mi>i</mi><mo>&prime;</mo></msup><mo>,</mo><mn>2</mn><mo>)</mo></mrow><mo>*</mo><mi>fR</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mi>fR</mi><msup><mrow><mo>(</mo><msup><mi>i</mi><mo>&prime;</mo></msup><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mn>3</mn></msup><mo>;</mo></mrow></math>]]></maths>⑤确定最终的团簇半径分布fRo(i,j)和团簇尺寸分布fNo(n,j):将模拟得到的模拟离子能谱fE(i,j)与实验得到的实验离子能谱fEo(m,j)进行比较,通过不断的修改模拟的初始团簇参数:团簇尺寸的最大值Nmax、团簇尺寸峰值Npeak和团簇对数尺寸的标准方差σ,然后重复步骤②、③、④,直到所述的模拟离子能谱fE(i,j)与实验离子能谱fEo(m,j)重合为止,得到与所述的最终模拟离子能谱fE(i,j)相应的最终团簇半径分布fRo(i,j)和最终团簇尺寸分布fNo(n,j);⑥完成团簇尺寸的标定:利用所述的最终团簇半径分布fRo(i,j)和最终团簇尺寸分布fNo(n,j)按下列公式计算团簇尺寸和半径,团簇的平均半径是Rmean:<maths num="0008"><![CDATA[<math><mrow><mi>Rmean</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><mi>fRo</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mi>fRo</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><mi>fRo</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>团簇的平均尺寸是Nmean:<maths num="0009"><![CDATA[<math><mrow><mi>Nmean</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mi>max</mi></mrow></munderover><mi>fNo</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mi>fNo</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mi>max</mi></mrow></munderover><mi>fNo</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>完成团簇尺寸的标定。
地址 201800 上海市嘉定区800-211邮政信箱