CN105021292A - 一种获取同步辐射部分相干光的传播特性的方法 - Google Patents
一种获取同步辐射部分相干光的传播特性的方法 Download PDFInfo
- Publication number
- CN105021292A CN105021292A CN201510427324.8A CN201510427324A CN105021292A CN 105021292 A CN105021292 A CN 105021292A CN 201510427324 A CN201510427324 A CN 201510427324A CN 105021292 A CN105021292 A CN 105021292A
- Authority
- CN
- China
- Prior art keywords
- coherent light
- synchrotron radiation
- partially coherent
- plane
- mutual intensity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供一种获取同步辐射部分相干光的传播特性的方法,包括:1、将物面、像面、以及光学元件的入射面和出射面平均分为M*N个小面元;2、测量部分相干光在物面上的互强度分布,通过离散化模型获取部分相干光从物面传播到光学元件入射面上的互强度分布;3、获取光学元件的振幅透过率分布,并根据入射面互强度分布和振幅透过率分布获取部分相干光从光学元件入射面传播到光学元件出射面上的互强度分布;4、通过离散化模型,根据出射面互强度分布获取部分相干光从同步辐射光学元件出射面传播到像面上的互强度分布;5、根据像面互强度分布获取部分相干光的传播特性。本发明能够准确地实现部分相干光经过同步辐射光学元件的传播特性模拟。
Description
技术领域
本发明涉及同步辐射技术领域,尤其涉及一种获取同步辐射部分相干光的传播特性的方法。
背景技术
同步辐射技术在现代科学中的应用越来越广泛,越来越多的同步辐射装置相继建立起来,而且对同步辐射性能的要求也越来越高。经过三代同步辐射的发展,目前世界上的同步辐射装置有70余台,其中第三代同步辐射装置有25台,并且还在不断建设中。在光束线建设中,需要了解具有相干性的X光经过光束线传播之后相干性能的变化,从而优化具体的光束线设计。目前,技术人员已经开发设计了一些模型和软件来模拟分析X光在光束线中的传播。
例如,文献1【同步辐射光束线中部分相干光的传播研究,作者:孟祥雨;郭春磊;王勇;吴衍青;邰仁忠刊名:光学学报出版日期:2013期号:第7期页码:318-325ISSN:02532239】提出了一种模拟同步辐射部分相干光的传播特性的离散化模型,其特点是基于统计光学理论,建立适用于同步辐射软X射线的部分相干光的传播模型。该模型利用互强度描述部分相干光,并采取分割波前以及互强度叠加的方法以及合理的近似,以定量分析部分相干光传播一定距离之后相干性能的变化,以及传播距离、光源相干长度、光源尺寸等因素对相干性能的影响。然而,该模型只能模拟部分相干光在自由空间或者经过狭缝调制后的传播,并不能分析其经过聚焦反射镜、光栅、波带片等同步辐射光学元件对相干性的影响。
又如,ESRF(欧洲同步辐射光源)的SHADOW3软件能够分析部分相干光经过同步辐射光学元件的传播特性,其原理是利用互强度描述部分相干光,并将物面的互强度分解成N个复振幅叠加的形式,以计算每个复振幅经过光学元件的传播,然后再将传播后的复振幅叠加成互强度,最终得到像面上的互强度分布。然而,该软件将物面的互强度分解成N个复振幅叠加的形式为近似处理,并不能准确地描述部分相干光在物面上的互强度分布。
发明内容
针对上述现有技术的不足,本发明提供一种获取同步辐射部分相干光的传播特性的方法,以更加准确地实现部分相干光经过同步辐射光学元件的传播特性模拟。
为了实现上述目的,本发明采用如下技术方案:
一种获取同步辐射部分相干光的传播特性的方法,用于获取部分相干光经过同步辐射光学元件的传播特性,包括确定所述部分相干光的物面和像面,还包括以下步骤:
步骤1,将所述物面、所述像面、以及所述同步辐射光学元件的入射面和出射面平均分为M*N个小面元,其中,M表示行数,N表示列数;
步骤2,测量所述部分相干光在所述物面上的物面互强度分布,然后通过离散化模型,根据所述物面互强度分布获取所述部分相干光从所述物面传播到所述同步辐射光学元件的入射面上的入射面互强度分布;
步骤3,获取所述同步辐射光学元件的振幅透过率分布,并根据式(4)获取所述部分相干光从所述同步辐射光学元件的入射面传播到所述同步辐射光学元件的出射面上的出射面互强度分布:
J(Mij,Mkl)=t(Mij)t*(Mkl)J(Qij,Qkl) (4),
在式(4)中,J(Mij,Mkl)表示所述同步辐射光学元件的出射面上的第i行j列小面元的中心点Mij与第k行l列小面元的中心点Mkl的互强度,J(Qij,Qkl)表示所述同步辐射光学元件的入射面上的第i行j列小面元的中心点Qij与第k行l列小面元的中心点Qkl的互强度,t(Mij)表示所述点Mij的振幅透过率,t*(Mkl)表示所述点Mkl的振幅透过率的共轭,其中,1≤i≤M,1≤j≤N,1≤k≤M,1≤l≤N;
步骤4,通过所述离散化模型,根据所述出射面互强度分布获取所述部分相干光从所述同步辐射光学元件的出射面传播到所述像面上的像面互强度分布;以及,
步骤5,根据所述像面互强度分布获取所述部分相干光的传播特性。
进一步地,所述同步辐射光学元件为聚焦光学元件或光栅。
进一步地,所述聚焦光学元件为波荡器、平面镜、狭缝、柱面镜或超环面镜。
优选地,当所述同步辐射光学元件为聚焦光学元件时,将其作为薄透镜以计算其振幅透过率分布。
综上所述,本发明首先测量物面上的互强度分布,再计算物面到同步辐射光学元件入射面的互强度分布,然后计算同步辐射光学元件入射面到出射面的互强度分布,最终得到所述出射面到像面的互强度分布,根据所述像面上的互强度分布即可模拟部分相干光的传播特性。本发明由于采用了离散化的方式分割物面、光学元件的入射面和出射面、以及像面上的互强度分布,从而减少了对各面互强度的近似处理,能够更准确地模拟部分相干光经过光学元件后的传播。
附图说明
图1为部分相干光经过同步辐射光学元件的传播示意图。
具体实施方式
下面结合附图,给出本发明的较佳实施例,并予以详细描述。
本发明的获取同步辐射部分相干光的传播特性的方法可用于模拟部分相干光经过同步辐射光学元件的传播。如图1所示,具有相干性的X光(即部分相干光)从光源实物面1(简称物面1)经同步辐射光学元件2传播到光源成像面3(简称像面3)。在传播过程中,波传播到某一位置处等相位面组成的曲面称为波前,最前的波前即为物面1,最后的波前即为像面3。基于文献1记载的离散化方法,本发明将物面1、像面3、以及同步辐射光学元件2的入射面和出射面平均分为M*N(M表示行数,N表示列数)个小面元,在每个小面元内部,振幅和相位处处相等而且是全相干光,所以各面上的小面元内部任意一点与另一个小面元内部任意一点的互强度均相等,都等于两个小面元中心点的互强度。基于此,本发明通过如下步骤处理部分相干光经过同步辐射光学元件2的互强度传播:
首先,测量部分相干光在物面1上的互强度分布(为了与其它面上的互强度分布区分,可称为物面1互强度分布),其中,第i行j列(1≤i≤M,1≤j≤N)的小面元Pij与第k行l列(1≤k≤M,1≤l≤N)的小面元Pkl的互强度记为J(Pij,Pkl)。当光束从物面1传播到同步辐射光学元件2的入射面时,由于物面1到光学元件2入射面的传播属于自由空间传播,入射面上的互强度分布(为了与其它面上的互强度分布区分,可称为入射面互强度分布)可通过式(1),即离散化模型,计算得到:
在式(1)中,J(Qij,Qkl)表示光学元件2入射面上的小面元中心点Qij与Qkl的互强度,J(Pij,Pkl)表示物面1上的小面元中心点Pij与Pkl的互强度,非下标的i是虚数单位 表示物面1上的小面元中心点Pij与光学元件2入射面上的小面元中心点Qij之间的距离,表示物面1上的小面元中心点Pkl与光学元件2入射面上的小面元中心点Qkl之间的距离,χ(θ1)为X光从Pij到Qij的倾斜因子,χ(θ2)为X光从Pij到Qij的倾斜因子,λ为X光的波长,每个积分的范围都是小面元内部。
当X光到达同步辐射光学元件2后,自其入射面向出射面传播。此处的同步辐射光学元件2为聚焦光学元件(例如波荡器、平面镜、狭缝、柱面镜、超环面镜等)或光栅,当为聚焦光学元件2时,由于同步辐射聚焦光学元件2具有较好的聚焦能力以及良好的光学特性,可以将其近似看成理想的薄透镜,这在本申请中首次提出。
根据傅里叶光学可知,同步辐射光学元件2的坐标(i,j)点的振幅透过率t(i,j)为:
在式(2)中,a和b分别为光学元件2的长度和宽度,fx和fy为光学元件2在x方向和y方向的焦距。如果光学元件2为光栅,则t(i,j)=P(i,j),因为光栅没有聚焦的功能,振幅透过率没有变化,只受光栅长宽的影响。
由于入射到光学元件2上的坐标(i,j)点的光线以相同的坐标(u=i,v=j)从另一边出射,因此振幅散布函数为:
K(u,v,i,j)=δ(u-i,v-j)t(i,j) (3),
在式(3)中,δ(u-i,v-j)为狄拉克函数,该函数表示当u=i,v=j才有意义,否则K值为零。将振幅散布函数代入现有的光学元件互强度传播公式,即可得到:
J(Mij,Mkl)=∫∫∫∫J(Qij,Qkl)K(Mij,Qij)K*(Mkl,Qkl)dQijdQkl
=t(Mij)t*(Mkl)J(Qij,Qkl) (4),
在式(4)中,J(Mij,Mkl)表示光学元件2出射面上的小面元中心点Mij与Mkl的互强度,J(Qij,Qkl)表示光学元件2入射面上的小面元中心点Qij与Qkl的互强度,t(Mij)表示点Mij的振幅透过率,t*(Mkl)表示点Mkl的振幅透过率的共轭。
当X光从光学元件2出射后,由于光学元件2出射面到像面3的传播同样属于自由空间传播,像面3上的互强度分布(为了与其它面上的互强度分布区分,可称为像面互强度分布)同样可通过离散化模型计算得到。
当X光经过的光学元件2为聚焦光学元件时,其波长未发生变化,像面上的互强度分布如式(5)所示:
在式(5)中,J(Nij,Nkl)表示像面3上的小面元中心点Nij与Nkl的互强度,J(Mij,Mkl)表示光学元件2出射面上的小面元中心点Mij与Mkl的互强度,非下标的i是虚数单位 表示光学元件2出射面上的小面元中心点Mij与像面3上的小面元中心点Nij之间的距离,表示光学元件2出射面上的小面元中心点Mkl与像面3上的小面元中心点Nkl之间的距离,χ(θ1)为X光从Mij到Nij的倾斜因子,χ(θ2)为X光从Mij到Nij的倾斜因子,λ为X光的波长,每个积分的范围都是小面元内部。
当X光经过的光学元件2为光栅时,X光将发生色散形成多束波长不同的单色光,且像面上不同波长的光产生的互强度会叠加,即,像面上的互强度分布如式(6)表示:
其中,
在上式中,J(Nij,Nkl)表示像面3上的小面元中心点Nij与Nkl的互强度,J(Mij,Mkl)表示光学元件2出射面上的小面元中心点Mij与Mkl的互强度,非下标的i是虚数单位 表示光学元件2出射面上的小面元中心点Mij与像面3上的小面元中心点Nij之间的距离,表示光学元件2出射面上的小面元中心点Mkl与像面3上的小面元中心点Nkl之间的距离,χ(θ1)为X光从Mij到Nij的倾斜因子,χ(θ2)为X光从Mij到Nij的倾斜因子,λh为形成的相应单色光的波长,n表示分散成的单色光的数量,每个积分的范围都是小面元内部。
按照上述步骤获得像面3上的互强度分布后,即可定量分析具有相干性的X光的特性,即,其传播一定距离之后相干性能的变化,以及传播距离、光源相干长度、光源尺寸等因素对相干性能的影响。
以上所述的内容,仅为本发明的较佳实施例,并非用以限定本发明的范围,本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰,皆落入本发明专利的权利要求保护范围。
Claims (4)
1.一种获取同步辐射部分相干光的传播特性的方法,用于获取部分相干光经过同步辐射光学元件的传播特性,包括确定所述部分相干光的物面和像面,其特征在于,还包括以下步骤:
步骤1,将所述物面、所述像面、以及所述同步辐射光学元件的入射面和出射面平均分为M*N个小面元,其中,M表示行数,N表示列数;
步骤2,测量所述部分相干光在所述物面上的物面互强度分布,然后通过离散化模型,根据所述物面互强度分布获取所述部分相干光从所述物面传播到所述同步辐射光学元件的入射面上的入射面互强度分布;
步骤3,获取所述同步辐射光学元件的振幅透过率分布,并根据式(4)获取所述部分相干光从所述同步辐射光学元件的入射面传播到所述同步辐射光学元件的出射面上的出射面互强度分布:
J(Mij,Mkl)=t(Mij)t*(Mkl)J(Qij,Qkl) (4),
在式(4)中,J(Mij,Mkl)表示所述同步辐射光学元件的出射面上的第i行j列小面元的中心点Mij与第k行l列小面元的中心点Mkl的互强度,J(Qij,Qkl)表示所述同步辐射光学元件的入射面上的第i行j列小面元的中心点Qij与第k行l列小面元的中心点Qkl的互强度,t(Mij)表示所述点Mij的振幅透过率,t*(Mkl)表示所述点Mkl的振幅透过率的共轭,其中,1≤i≤M,1≤j≤N,1≤k≤M,1≤l≤N;
步骤4,通过所述离散化模型,根据所述出射面互强度分布获取所述部分相干光从所述同步辐射光学元件的出射面传播到所述像面上的像面互强度分布;以及,
步骤5,根据所述像面互强度分布获取所述部分相干光的传播特性。
2.根据权利要求1所述的获取同步辐射部分相干光的传播特性的方法,其特征在于,所述同步辐射光学元件为聚焦光学元件或光栅。
3.根据权利要求2所述的获取同步辐射部分相干光的传播特性的方法,其特征在于,所述聚焦光学元件为波荡器、平面镜、狭缝、柱面镜或超环面镜。
4.根据权利要求2或3所述的获取同步辐射部分相干光的传播特性的方法,其特征在于,当所述同步辐射光学元件为聚焦光学元件时,将其作为薄透镜以计算其振幅透过率分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510427324.8A CN105021292A (zh) | 2015-07-20 | 2015-07-20 | 一种获取同步辐射部分相干光的传播特性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510427324.8A CN105021292A (zh) | 2015-07-20 | 2015-07-20 | 一种获取同步辐射部分相干光的传播特性的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105021292A true CN105021292A (zh) | 2015-11-04 |
Family
ID=54411431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510427324.8A Pending CN105021292A (zh) | 2015-07-20 | 2015-07-20 | 一种获取同步辐射部分相干光的传播特性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105021292A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814384A (zh) * | 2015-11-27 | 2017-06-09 | 华北电力大学 | 核电厂点源辐射源强逆推方法及点源辐射源强逆推系统 |
CN115035035A (zh) * | 2022-05-09 | 2022-09-09 | 中国科学院上海高等研究院 | 一种x射线经过非均匀短程介质的传播的计算方法 |
CN116430581A (zh) * | 2023-03-03 | 2023-07-14 | 中国科学院上海高等研究院 | 部分相干光经过自由空间及光学元件的传播的计算方法 |
-
2015
- 2015-07-20 CN CN201510427324.8A patent/CN105021292A/zh active Pending
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814384A (zh) * | 2015-11-27 | 2017-06-09 | 华北电力大学 | 核电厂点源辐射源强逆推方法及点源辐射源强逆推系统 |
CN106814384B (zh) * | 2015-11-27 | 2019-03-15 | 华北电力大学 | 核电厂点源辐射源强逆推方法及点源辐射源强逆推系统 |
CN115035035A (zh) * | 2022-05-09 | 2022-09-09 | 中国科学院上海高等研究院 | 一种x射线经过非均匀短程介质的传播的计算方法 |
CN115035035B (zh) * | 2022-05-09 | 2024-05-10 | 中国科学院上海高等研究院 | 一种x射线经过非均匀短程介质的传播的计算方法 |
CN116430581A (zh) * | 2023-03-03 | 2023-07-14 | 中国科学院上海高等研究院 | 部分相干光经过自由空间及光学元件的传播的计算方法 |
CN116430581B (zh) * | 2023-03-03 | 2024-05-24 | 中国科学院上海高等研究院 | 部分相干光经过自由空间及光学元件的传播的计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schultheiss et al. | Hanbury Brown and Twiss measurements in curved space | |
CN101980056B (zh) | 基于m2因子和光闪烁指数确定大气湍流参数方法及装置 | |
Luquet et al. | Long range numerical simulation of acoustical shock waves in a 3D moving heterogeneous and absorbing medium | |
CN103884436B (zh) | 光束相位在线测量装置和测量方法 | |
Fu et al. | Fiber-optic curvature sensor with optimized sensitive zone | |
CN105021292A (zh) | 一种获取同步辐射部分相干光的传播特性的方法 | |
JP2014190705A5 (zh) | ||
Kulkarni et al. | Classical model of spontaneous parametric down-conversion | |
CN104090376A (zh) | 高数值孔径短焦距台阶相位型厚fzp的设计方法 | |
CN104634542A (zh) | 大口径光学元件二次曝光相位测量装置及测量方法 | |
CN106644104B (zh) | 一种基于谱反演法的离散雨滴介质的相位屏建模方法 | |
Banakh et al. | Wind profiling based on the optical beam intensity statistics in a turbulent atmosphere | |
CN103777263A (zh) | 亚波长矩形单周期光栅结构的制作方法及偶数分束器 | |
Bos et al. | Wave optics simulations of plane-wave propagation in non-Kolmogorov turbulence | |
Li et al. | Flexural wave control via the profile modulation of non-uniform Timoshenko beams | |
Cheng et al. | Propagation of scalar and vector vortex beams through turbulent atmosphere | |
CN106290094A (zh) | 应用于扬尘颗粒物在线监测的mie散射快速计算方法 | |
Toselli et al. | Scintillation analysis of LIDAR systems operating in weak-to-strong non-Kolmogorov turbulence: unresolved target case | |
CN105046010A (zh) | 一种卫星杂散光的分析方法 | |
CN102662237B (zh) | 可控小光程差干涉系统 | |
An et al. | Large survey camera collimation and alignment at the observation site by curvature sensing | |
He et al. | Imaging quality evaluation method of pixel coupled electro-optical imaging system | |
Fathnan et al. | Efficiency calculation of non-periodic metasurface based on modified near field to far field transformation | |
Li | 3D display based on complete digital optical phase conjugation | |
Fu et al. | Integrated simulation method of the scattering medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151104 |