CN107664648B - A kind of X-ray differential phase contrast microscopic system and its two-dimensional imaging method - Google Patents
A kind of X-ray differential phase contrast microscopic system and its two-dimensional imaging method Download PDFInfo
- Publication number
- CN107664648B CN107664648B CN201610617865.1A CN201610617865A CN107664648B CN 107664648 B CN107664648 B CN 107664648B CN 201610617865 A CN201610617865 A CN 201610617865A CN 107664648 B CN107664648 B CN 107664648B
- Authority
- CN
- China
- Prior art keywords
- image
- grating
- sample
- annular
- angle signal
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 106
- 238000004458 analytical method Methods 0.000 claims abstract description 80
- 230000004044 response Effects 0.000 claims description 99
- 238000010521 absorption reaction Methods 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 26
- 238000006073 displacement reaction Methods 0.000 claims description 25
- 230000014509 gene expression Effects 0.000 claims description 22
- 230000003287 optical effect Effects 0.000 claims description 21
- 238000005316 response function Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 18
- 230000010363 phase shift Effects 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 14
- 230000009131 signaling function Effects 0.000 claims description 12
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 3
- 239000010931 gold Substances 0.000 claims description 3
- 229910052737 gold Inorganic materials 0.000 claims description 3
- 229910001385 heavy metal Inorganic materials 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 230000007423 decrease Effects 0.000 description 15
- 230000001427 coherent effect Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 239000002245 particle Substances 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- -1 bubbles Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000005338 frosted glass Substances 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000013081 microcrystal Substances 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000002135 phase contrast microscopy Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000005469 synchrotron radiation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/083—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/045—Investigating materials by wave or particle radiation combination of at least 2 measurements (transmission and scatter)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/605—Specific applications or type of materials phases
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Toxicology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种X射线微分相位衬度显微镜系统及其二维成像方法,所述X射线微分相位衬度显微镜系统包括:用于产生X射线的光源;以及沿X射线传播方向依次设置的聚光镜、中心光阑、分束光栅、针孔、样品台、物镜、环形分析光栅和成像探测器。本发明的有益效果是:该X射线微分相位衬度显微镜系统,仅在传统的X射线显微镜中增加分束光栅和环形分析光栅,就能实现相位衬度定量成像,具有结构简单、易于推广的优点。另外,可以将X射线光源、聚光镜、分束光栅集成为一个X射线环形栅源元件,则整个X射线微分相位衬度显微镜系统长度可以进一步缩短,不仅可以降低X射线显微镜系统的制造成本,而且光的利用效率也能进一步提高。
The invention discloses an X-ray differential phase contrast microscope system and a two-dimensional imaging method thereof. The X-ray differential phase contrast microscope system includes: a light source for generating X-rays; Condenser, central diaphragm, beamsplitter grating, pinhole, sample stage, objective, annular analysis grating and imaging detector. The beneficial effects of the present invention are: the X-ray differential phase contrast microscope system can realize phase contrast quantitative imaging only by adding a beam splitting grating and a ring analysis grating to the traditional X-ray microscope, and has a simple structure and is easy to popularize advantage. In addition, the X-ray light source, condenser lens, and beam splitting grating can be integrated into an X-ray ring grating source element, and the length of the entire X-ray differential phase contrast microscope system can be further shortened, which not only reduces the manufacturing cost of the X-ray microscope system, but also The utilization efficiency of light can also be further improved.
Description
技术领域technical field
本发明涉及纳米分辨X射线显微镜成像技术领域,尤其涉及X射线微分相位衬度显微镜系统及其二维成像方法。The invention relates to the technical field of nano-resolution X-ray microscope imaging, in particular to an X-ray differential phase contrast microscope system and a two-dimensional imaging method thereof.
背景技术Background technique
物质对X射线的作用可以用折射率代表,n=1-δ+iβ,其中δ是折射率实部衰减率,β是吸收项。δ的物理意义是,经过单位长度物质的X射线相对于经过单位长度真空的X射线产生的波面移动;β的物理意义是经过单位长度物质的X射线产生的波幅下降。根据δ和β的物理意义,可得X射线经过一个样品后,其复振幅可以表达为其中The effect of matter on X-rays can be represented by the refractive index, n=1-δ+iβ, where δ is the attenuation rate of the real part of the refractive index, and β is the absorption term. The physical meaning of δ is that the X-rays passing through a unit length of matter move relative to the wavefront generated by the X-rays passing through a unit length of vacuum; the physical meaning of β is the decrease in the amplitude of X-rays passing through a unit length of matter. According to the physical meaning of δ and β, after the X-ray passes through a sample, its complex amplitude can be expressed as in
为经过样品的X射线相位相对经过相同长度真空的X射线相位的变化量,简称相移,l为X射线经过样品的路径;因为在样品外的积分值为零,所以积分的上下限可以扩展到无穷;It is the change amount of the X-ray phase passing through the sample relative to the X-ray phase passing through the vacuum of the same length, referred to as the phase shift, and l is the path of the X-ray passing through the sample; because the integral value outside the sample is zero, the upper and lower limits of the integral can be expanded to infinity;
为吸收,μ为线性吸收系数。因为在硬X射线波段轻元素的δ比β大三个量级以上,所以相位变化有可能产生比吸收变化大得多的光强信号。荷兰科学家泽尼克(Zernik)是相位衬度成像的第一人。早在1935年,他就在可见光波段提出了相位衬度显微镜的理论和方法,为此他荣获了1953年的诺贝尔物理学奖。目前人们已经成功地将泽尼克相位衬度显微镜方法推广到以波带片作为物镜的X射线显微镜,研制成功利用相移环获得相位衬度的X射线相位衬度显微镜,如图1所示,该显微镜按照X射线传播方向依次包括X射线光源10、聚光镜1、中心光阑2、针孔3、样品台4、物镜5、相移环6和成像探测器7。然而,泽尼克发明的相位衬度成像存在明显的不足:其应用范围局限于吸收可以忽略的弱相位样品。当吸收不能忽略时,它不能区分吸收产生的光强响应和相移引起的光强响应;当样品不是弱相位物体,特别样品的相移超过π/4时,相移和光强响应的关系不再满足线性关系,光强响应将随着相位周期产生振荡。这些不足为定量研究样品密度结构带来困难。is the absorption, and μ is the linear absorption coefficient. Because δ is more than three orders of magnitude larger than β for light elements in the hard X-ray band, phase changes have the potential to generate light intensity signals that are much larger than absorption changes. Dutch scientist Zernik (Zernik) is the first person in phase contrast imaging. As early as 1935, he proposed the theory and method of phase contrast microscopy in the visible light band, for which he won the 1953 Nobel Prize in Physics. At present, people have successfully extended the Zernike phase contrast microscope method to the X-ray microscope with a zone plate as the objective lens, and successfully developed an X-ray phase contrast microscope that uses a phase shift ring to obtain phase contrast, as shown in Figure 1. The microscope includes an X-ray light source 10, a condenser lens 1, a central diaphragm 2, a pinhole 3, a sample stage 4, an objective lens 5, a phase shift ring 6 and an imaging detector 7 in sequence according to the X-ray propagation direction. However, the phase-contrast imaging invented by Zernike has obvious shortcomings: its application range is limited to weak-phase samples with negligible absorption. When the absorption cannot be ignored, it cannot distinguish between the light intensity response caused by absorption and the light intensity response caused by phase shift; when the sample is not a weak phase object, especially when the phase shift of the sample exceeds π/4, the relationship between phase shift and light intensity response The linear relationship is no longer satisfied, and the light intensity response will oscillate with the phase period. These deficiencies make it difficult to quantitatively study the density structure of samples.
发明内容Contents of the invention
本发明的目的之一,提供一种能够定量研究样品密度结构的灵敏度高于传统吸收X射线显微镜的X射线微分相位衬度显微镜系统。One of the objectives of the present invention is to provide an X-ray differential phase contrast microscope system capable of quantitatively studying the density structure of a sample with higher sensitivity than traditional absorption X-ray microscopes.
本发明的目的之二,提供一种用于所述高灵敏度的X射线微分相位衬度显微镜系统的二维成像方法。The second object of the present invention is to provide a two-dimensional imaging method for the high-sensitivity X-ray differential phase contrast microscope system.
根据本发明的一方面,公开了一种X射线微分相位衬度显微镜系统,其可以包括:用于产生X射线的光源;以及沿X射线传播方向依次设置的聚光镜、中心光阑、分束光栅、针孔、样品台、物镜、环形分析光栅和成像探测器;其中,所述产生X射线的光源为单色X射线光源;所述聚光镜为椭球毛细管、锥形毛细管、波带片或者其它具有聚焦功能的X射线光学元件,用于产生照明样品的聚焦光束;所述中心光阑位于聚光镜出口中心位置,用于阻挡来自X射线光源的直通光,形成照明样品的空心锥光束;所述分束光栅为周期为微米量级的吸收光栅,位于聚光镜出口的中心光阑后,用于对来自聚光镜和中心光阑的空心锥光束进行分束,使其成为具有空间周期结构的空心锥光束;所述针孔,位于分束光栅后,用于阻挡来自X射线光源的直通光和杂散光;所述样品台,位于物面,用于承载样品,并能平移和旋转样品;所述物镜为波带片或者其它具有X射线透镜成像功能的光学元件,用于把位于物面的具有纳米量级的样品结构放大为像面上具有微米量级结构的像,并使分束光栅被空心锥光束照明的环形部分在物镜后焦面附近形成环形光栅像及环形光栅像光束;所述环形分析光栅为周期为微米或者亚微米量级的吸收光栅,位于物镜后焦面附近,其形状和尺寸与所述分束光栅在物镜后焦面附近的环形光栅像相同,用于对所述物镜后焦面附件的环形光栅像光束进行滤波;所述成像探测器位于像面,用于拍摄样品的二维放大像。According to one aspect of the present invention, an X-ray differential phase contrast microscope system is disclosed, which may include: a light source for generating X-rays; and a condenser lens, a central diaphragm, and a beam splitting grating arranged in sequence along the X-ray propagation direction , pinhole, sample stage, objective lens, annular analysis grating and imaging detector; Wherein, the light source that described X-ray generation is monochromatic X-ray light source; Described condenser is ellipsoid capillary, tapered capillary, zone plate or other The X-ray optical element with focusing function is used to generate a focused beam for illuminating the sample; the central diaphragm is located at the center of the condenser exit and is used to block the through light from the X-ray light source to form a hollow cone beam for illuminating the sample; The beam splitting grating is an absorption grating with a period of micron scale, located behind the central diaphragm at the exit of the condenser, and is used to split the hollow cone beam from the condenser and the central diaphragm, making it a hollow cone beam with a spatially periodic structure ; The pinhole, located behind the beam splitting grating, is used to block the through light and stray light from the X-ray light source; the sample stage, located on the object plane, is used to carry the sample, and can translate and rotate the sample; the objective lens It is a zone plate or other optical element with X-ray lens imaging function, which is used to magnify the nanoscale sample structure on the object plane into an image with micron scale structure on the image plane, and make the beam splitting grating hollow The annular part of the cone beam illumination forms an annular grating image and an annular grating image beam near the rear focal plane of the objective lens; the annular analysis grating is an absorption grating with a period of micron or submicron order, located near the rear focal plane of the objective lens, and its shape and The size is the same as the annular grating image near the rear focal plane of the objective lens, and is used to filter the annular grating image beam attached to the rear focal plane of the objective lens; the imaging detector is located on the image plane, and is used to photograph samples 2D magnified image of .
可选地,所述分束光栅和所述环形分析光栅的材质为金或其它重金属。Optionally, the beam splitting grating and the ring analysis grating are made of gold or other heavy metals.
可选地,在所述分束光栅和所述环形分析光栅厚度不大的情况下,分束光栅的厚度为π相移的奇数倍,环形分析光栅的厚度为π相移的偶数倍。Optionally, when the thicknesses of the beam splitting grating and the annular analysis grating are not large, the thickness of the beam splitting grating is an odd multiple of the π phase shift, and the thickness of the annular analysis grating is an even multiple of the π phase shift.
根据本发明的另一方面,公开了一种用于所述X射线微分相位衬度显微镜系统的二维成像方法,包括如下步骤:According to another aspect of the present invention, a two-dimensional imaging method for the X-ray differential phase contrast microscope system is disclosed, comprising the following steps:
S1,开启并调整光源:使所述光源产生的X射线光束对准聚光镜;S1, turn on and adjust the light source: aim the X-ray beam generated by the light source at the condenser;
S2,调整聚光镜、中心光阑和针孔:使聚光镜和中心光阑产生的空心锥光束对准所述样品台上承载样品位置,使中心光阑和针孔阻挡来自所述X射线光源的直通光和杂散光;S2, adjust the condenser, central diaphragm and pinhole: make the hollow cone beam generated by the condenser and the central diaphragm align with the position of the sample on the sample stage, so that the central diaphragm and the pinhole block the direct passage from the X-ray light source light and stray light;
S3,调整物镜:使物镜聚焦形成的成像光束对准位于像面的成像探测器;S3, adjusting the objective lens: aligning the imaging light beam formed by the objective lens with the imaging detector located on the image plane;
S4,馈入分束光栅:在聚光镜出口的中心光阑后,馈入分束光栅,使空心锥光束照明分束光栅的环形部分在物镜后焦面附近形成环形光栅像及环形光栅像光束;S4, feeding into the beam splitting grating: after the central diaphragm at the exit of the condenser, feed into the beam splitting grating, so that the hollow cone beam illuminates the annular part of the beam splitting grating to form a ring grating image and a ring grating image beam near the rear focal plane of the objective lens;
S5,馈入并调整环形分析光栅:使环形分析光栅对准分束光栅在物镜后焦面附近形成的环形光栅像,调整环形分析光栅栅条方向,使环形分析光栅栅条平行于分束光栅栅条;S5, feed and adjust the ring analysis grating: make the ring analysis grating align with the ring grating image formed by the beam splitting grating near the rear focal plane of the objective lens, adjust the direction of the ring analysis grating bars so that the ring analysis grating bars are parallel to the beam splitting grating grid;
S6,测角度信号响应曲线:沿着垂直于光轴和光栅栅条方向逐步移动分束光栅,使分束光栅的环形光栅像相对环形分析光栅发生剪切位移,用成像探测器在像面测得亮度随环形光栅像位移变化的角度信号响应曲线;S6, angle measurement signal response curve: move the beam splitting grating step by step along the direction perpendicular to the optical axis and the grating bars, so that the annular grating image of the beam splitting grating has a shear displacement relative to the annular analysis grating, and the imaging detector is used to measure the angle on the image plane. Obtain the angle signal response curve of the brightness changing with the displacement of the annular grating image;
S7,用余弦曲线拟合角度信号响应曲线:角度信号响应曲线类似于余弦曲线,用余弦曲线拟合测得的角度信号响应曲线,使其获得余弦曲线的解析表达;S7, fitting the angle signal response curve with a cosine curve: the angle signal response curve is similar to a cosine curve, and the measured angle signal response curve is fitted with a cosine curve to obtain an analytical expression of a cosine curve;
S8,拍摄样品二维放大像:把分束光栅的环形光栅像分别固定在所述角度信号响应曲线的谷位、上坡位、峰位、下坡位,在样品台上放置样品,拍摄样品的谷位放大像、上坡放大像、峰位放大像、下坡放大像;S8, take a two-dimensional enlarged image of the sample: fix the annular grating image of the beam splitting grating at the valley position, uphill position, peak position, and downhill position of the angle signal response curve respectively, place the sample on the sample stage, and photograph the sample The valley zoom image, uphill zoom image, peak zoom image, downhill zoom image;
S9,提取样品在物空间的定量二维像:从样品的谷位放大像、上坡放大像、峰位放大像、下坡放大像中提取样品在物空间的吸收像、折射像和散射像;从四种放大像中的任意三种放大像中提取样品在物空间的吸收像、折射像和散射像。S9, extract the quantitative two-dimensional image of the sample in the object space: extract the absorption image, refraction image and scattering image of the sample in the object space from the valley magnification image, uphill magnification image, peak magnification image, and downhill magnification image of the sample ; Extract the absorption, refraction and scattering images of the sample in the object space from any three of the four magnification images.
在所述步骤S6中,考虑到样品对焦深细光束的的折射和散射影响,即考虑到样品产生的角度信号对成像过程的影响,故样品用角度信号函数f(xo,yo,ψx)表达为In the step S6, considering the influence of the refraction and scattering of the sample's focal depth thin beam, that is, considering the influence of the angle signal generated by the sample on the imaging process, the sample uses the angle signal function f(x o , y o , ψ x ) is expressed as
其中(xo,yo,zo)为样品处空间坐标,M(xo,yo)、θx(xo,yo)和分别为样品的吸收、折射角和散射角方差,也是X射线微分相位衬度显微镜二维定量成像中要求的样品在物空间的吸收像、折射像和散射像,它们的表达式为Where (x o , y o , z o ) are the spatial coordinates of the sample, M(x o , y o ), θ x (x o , y o ) and are the variances of the absorption, refraction and scattering angles of the sample, respectively, and are also the absorption, refraction and scattering images of the sample in the object space required in the two-dimensional quantitative imaging of the X-ray differential phase contrast microscope, and their expressions are
其中(xo,yo,zo)为物面处空间坐标,μ为线性吸收系数,δ是折射率实部衰减率,ωx是垂直于光栅栅条方向的线性散射系数。Where (x o , y o , z o ) are the spatial coordinates at the object plane, μ is the linear absorption coefficient, δ is the attenuation rate of the real part of the refractive index, and ω x is the linear scattering coefficient perpendicular to the direction of the grating bars.
在所述步骤S6中测得的是角度信号响应曲线,其用拟合余弦曲线解析表达为What measure in described step S6 is the angle signal response curve, and it is analytically expressed with fitting cosine curve as
或者表达为or expressed as
其中η为物镜衍射效率,B0为物面处入射光束亮度,xp为环形分析光栅处位置坐标,为环形光栅像位移相对物面的角位移,do为样品相对物镜的物距,p为环形分析光栅的周期,为角度信号响应曲线平均值,Rmax和Rmin分别为角度信号响应曲线的最大值和最小值,为角度信号响应曲线的可见度,n为调制参数,n=0,1,2,3分别对应谷位响应曲线、上坡响应曲线、峰位响应曲线和下坡响应曲线。Wherein η is the diffraction efficiency of the objective lens, B 0 is the brightness of the incident beam at the object plane, and xp is the position coordinate at the annular analysis grating, is the angular displacement of the annular grating image displacement relative to the object plane, d o is the object distance of the sample relative to the objective lens, p is the period of the annular analysis grating, is the average value of the angle signal response curve, R max and R min are the maximum and minimum values of the angle signal response curve, respectively, is the visibility of the angle signal response curve, n is the modulation parameter, and n=0, 1, 2, 3 correspond to the valley response curve, upslope response curve, peak response curve and downslope response curve respectively.
在所述步骤S8中,成像探测器拍摄的二维放大像用角度信号成像方程描述,角度信号成像方程是角度信号函数和角度信号响应函数的卷积,其表达式为:In said step S8, the two-dimensional enlarged image taken by the imaging detector is described by the angle signal imaging equation, and the angle signal imaging equation is the convolution of the angle signal function and the angle signal response function, and its expression is:
其中(xo,yo)为物面坐标,(xi,yi)为像面坐标,η为物镜衍射效率,B0为物面处入射光束亮度,n=0,1,2,3分别对应谷位放大成像方程、上坡放大成像方程、峰位放大成像方程和下坡放大成像方程。Where (x o , y o ) is the coordinates of the object plane, ( xi , y i ) is the coordinates of the image plane, η is the diffraction efficiency of the objective lens, B 0 is the brightness of the incident beam at the object plane, n=0,1,2,3 Corresponding to the valley magnification imaging equation, uphill magnification imaging equation, peak magnification imaging equation and downhill magnification imaging equation.
在所述步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线谷位,在样品台上放置样品,成像探测器拍摄到样品的谷位放大像,其表达式为In the step S8, the annular grating image of the beam splitting grating is fixed at the valley of the angle signal response curve, the sample is placed on the sample stage, and the imaging detector captures the enlarged image of the valley of the sample, and its expression is
在所述步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线上坡位,在样品台上放置样品,成像探测器拍摄到样品的上坡放大像,其表达式为In the step S8, the annular grating image of the beam splitting grating is fixed on the slope position of the angle signal response curve, the sample is placed on the sample stage, and the imaging detector captures the uphill enlarged image of the sample, its expression for
在所述步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线峰位,在样品台上放置样品,成像探测器拍摄到样品的峰位放大像,其表达式为In the step S8, the annular grating image of the beam-splitting grating is fixed at the peak position of the angle signal response curve, the sample is placed on the sample stage, and the imaging detector captures an enlarged image of the peak position of the sample, and its expression is
在所述步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线下坡位,在样品台上放置样品,成像探测器拍摄到样品的下坡放大像,其表达式为In the step S8, the annular grating image of the beam splitting grating is fixed at the downhill position of the angle signal response curve, the sample is placed on the sample stage, and the imaging detector captures the downhill enlarged image of the sample, its expression for
由于because
BV(xi,yi)+BP(xi,yi)=BU(xi,yi)+BD(xi,yi)B V ( xi ,y i )+B P ( xi ,y i )= BU ( xi ,y i )+B D ( xi ,y i )
所以谷位放大像、上坡放大像、峰位放大像、下坡放大像中三幅像是独立的,其中任一幅像都可以用其余三幅像表达出来;以上各式中,(xi,yi)为像面坐标,η为物镜的衍射效率,B0为照明物面光束的亮度,为角度信号响应曲线的平均值,Rmax和Rmin分别为角度信号响应曲线的最大值和最小值,为角度信号响应曲线的可见度,do为样品相对物镜的物距,p为环形分析光栅周期。Therefore, the three images in the valley magnified image, uphill magnified image, peak magnified image, and downhill magnified image are independent, and any one of them can be expressed by the other three images; among the above formulas, (x i , y i ) are the coordinates of the image plane, η is the diffraction efficiency of the objective lens, B 0 is the brightness of the light beam illuminating the object plane, is the average value of the angle signal response curve, R max and R min are the maximum and minimum values of the angle signal response curve respectively, is the visibility of the angle signal response curve, d o is the object distance of the sample relative to the objective lens, and p is the ring analysis grating period.
在所述步骤S9中,所述的从样品的谷位放大像、上坡放大像、峰位放大像、下坡放大像中提取样品在物空间的吸收像、折射像和散射像方法如下:In the step S9, the method of extracting the absorption image, refraction image and scattering image of the sample in the object space from the valley magnification image, uphill magnification image, peak position magnification image and downhill magnification image of the sample is as follows:
提取样品在物空间吸收像的公式为The formula for extracting the absorption image of the sample in the object space is
提取样品在物空间折射像的公式为The formula for extracting the refraction image of the sample in object space is
在散射可以忽略的条件下,提取样品在物空间折射像的公式简化为Under the condition that the scattering can be ignored, the formula for extracting the refraction image of the sample in the object space is simplified as
提取样品在物空间散射像的公式为The formula for extracting the scattering image of the sample in the object space is
在折射可以忽略的条件下,提取样品在物空间散射像的公式简化为Under the condition that the refraction can be ignored, the formula for extracting the scattering image of the sample in the object space is simplified as
以上各式中,(xo,yo)为物面坐标,(xi,yi)为像面坐标,do为样品相对物镜的物距, In the above formulas, (x o , y o ) are object plane coordinates, ( xi , y i ) are image plane coordinates, d o is the object distance of the sample relative to the objective lens,
作为一种可选的实施方式,所述X射线光源、聚光镜和分束光栅可以集成为一个周期为微米量级的X射线环形栅源的元件,此时在所述X射线环形栅源的出口处,沿光轴依次放置两个环形通光孔,使所述环形栅源发出的光束经过两个环形通光孔后,形成照明样品的空心锥光束,并在物镜后焦面附近形成环形栅源像及环形栅源像光束,所述环形分析光栅的形状和尺寸与所述环形栅源像相同。As an optional implementation, the X-ray light source, condenser lens and beam splitting grating can be integrated into an element of an X-ray ring grating source whose period is on the order of microns. At this time, at the exit of the X-ray ring grating source Place two ring-shaped light holes in sequence along the optical axis, so that the light beam emitted by the ring-shaped grating source passes through the two ring-shaped light holes to form a hollow cone beam that illuminates the sample, and forms a ring grating near the rear focal plane of the objective lens. The source image and the annular grating source image light beam, the shape and size of the annular analysis grating are the same as the annular grating source image.
本发明公开的X射线微分相位衬度显微镜系统,仅在传统的X射线显微镜中增加分束光栅和环形分析光栅,就能实现相位衬度定量成像,具有结构简单、易于推广的优点。另外,可以将X射线光源、聚光镜、分束光栅集成为一个X射线环形栅源元件,并在环形栅源出口沿光轴方向依次安置两个环形通光孔,则整个X射线微分相位衬度显微镜系统长度可以进一步缩短,不仅可以降低X射线显微镜系统的制造成本,而且光的利用效率也能进一步提高,具有广阔的应用前景。The X-ray differential phase contrast microscope system disclosed by the invention can realize phase contrast quantitative imaging only by adding a beam splitting grating and a ring analysis grating to the traditional X-ray microscope, and has the advantages of simple structure and easy popularization. In addition, the X-ray light source, condenser lens, and beam splitting grating can be integrated into an X-ray annular grating source element, and two annular apertures are arranged in sequence along the optical axis at the exit of the annular grating source, and the entire X-ray differential phase contrast The length of the microscope system can be further shortened, not only can reduce the manufacturing cost of the X-ray microscope system, but also can further improve the light utilization efficiency, which has broad application prospects.
附图说明Description of drawings
应当理解,下述所有附图中,相同的附图标记表示相同或相应的部件和特征。It should be understood that in all the following drawings, the same reference numerals indicate the same or corresponding components and features.
图1为现有技术中的X射线相位衬度显微镜系统的光路示意图。Fig. 1 is a schematic diagram of an optical path of an X-ray phase contrast microscope system in the prior art.
图2为根据本发明的一实施例的X射线微分相位衬度显微镜系统的示意图。FIG. 2 is a schematic diagram of an X-ray differential phase contrast microscope system according to an embodiment of the present invention.
图3为根据本发明的一实施例的X射线微分相位衬度显微镜系统的二维成像方法的流程框图。Fig. 3 is a flowchart of a two-dimensional imaging method of an X-ray differential phase contrast microscope system according to an embodiment of the present invention.
图4为根据本发明的实施例的X射线微分相位衬度显微镜系统的折射角信号成像原理图。Fig. 4 is a schematic diagram of refraction angle signal imaging of an X-ray differential phase contrast microscope system according to an embodiment of the present invention.
图5为根据本发明的实施例的X射线微分相位衬度显微镜系统的散射角信号成像原理图。Fig. 5 is a schematic diagram of scattering angle signal imaging of an X-ray differential phase contrast microscope system according to an embodiment of the present invention.
图6为X射线微分相位衬度显微镜系统的四种角度信号响应函数曲线图。Fig. 6 is a curve diagram of four angle signal response functions of the X-ray differential phase contrast microscope system.
图7为根据本发明的另一实施例的基于环形栅源的X射线微分相衬显微镜系统的示意图。Fig. 7 is a schematic diagram of an X-ray differential phase contrast microscope system based on a ring grid source according to another embodiment of the present invention.
附图中各附图标记的简单说明如下:A brief description of each reference sign in the accompanying drawings is as follows:
1:聚光镜1: Condenser
2:中心光阑2: Central aperture
3:针孔3: pinhole
4:样品台4: Sample stage
5:物镜5: objective lens
6:相移环6: Phase shift ring
7:成像探测器7: Imaging detector
3’:分束光栅3': beam splitting grating
6’:环形分析光栅6': Ring analysis grating
10:X射线光源10: X-ray light source
10’:X射线环形栅源10': X-ray ring grid source
20:第一环形通光孔20: The first annular light hole
30:第二环形通光孔。30: the second annular light hole.
具体实施方式Detailed ways
下面结合附图对本发明公开的实施例的X射线微分相位衬度显微镜系统及用于该系统的二维成像方法做详细的说明。The X-ray differential phase contrast microscope system of the embodiment disclosed in the present invention and the two-dimensional imaging method used in the system will be described in detail below with reference to the accompanying drawings.
图2为根据本发明的一实施例的X射线微分相位衬度显微镜系统的示意图。如图2所示,X射线微分相位衬度显微镜系统100按照X射线传播方向依次包括X射线光源10(图中未画出)、聚光镜1、中心光阑2、分束光栅3’、针孔3、样品台4、物镜5、环形分析光栅6’和成像探测器7构成。各元件的性质、结构和功能分别叙述如下:FIG. 2 is a schematic diagram of an X-ray differential phase contrast microscope system according to an embodiment of the present invention. As shown in Figure 2, the X-ray differential phase contrast microscope system 100 includes an X-ray light source 10 (not shown in the figure), a condenser lens 1, a central diaphragm 2, a beam splitting grating 3', and a pinhole in sequence according to the X-ray propagation direction. 3. A sample stage 4, an objective lens 5, an annular analysis grating 6' and an imaging detector 7 are formed. The nature, structure and function of each component are described as follows:
X射线光源10:X射线物镜一般是具有透镜成像功能的衍射光学元件,例如波带片,所以X射线微分相位衬度显微镜系统所用的X射线光源10可以为单色X射线光源,如电子撞击金属靶产生的特征谱线X射线光源、激光等离子X射线光源、同步辐射单色X射线光源。X-ray source 10: The X-ray objective lens is generally a diffractive optical element with lens imaging function, such as a zone plate, so the X-ray source 10 used in the X-ray differential phase contrast microscope system can be a monochromatic X-ray source, such as electron impact X-ray sources with characteristic spectral lines produced by metal targets, laser plasma X-ray sources, and synchrotron radiation monochromatic X-ray sources.
聚光镜1:可以为椭球毛细管、锥形毛细管,波带片或者其它具有聚焦功能的X射线光学元件,其作用是为样品提供聚焦照明光束。Condenser 1: It can be an ellipsoidal capillary, a tapered capillary, a zone plate or other X-ray optical elements with a focusing function, and its function is to provide a focused illumination beam for the sample.
中心光阑2:位于聚光镜1出口中心部分,用于阻挡来自所述X射线光源的过强的直通光照明样品和探测器,使照明样品的聚焦光束成为空心锥光束。Central diaphragm 2: located at the center of the exit of condenser 1, it is used to block the excessively strong through light from the X-ray source from illuminating the sample and detector, so that the focused beam illuminating the sample becomes a hollow cone beam.
分束光栅3’:可以是周期为微米量级的吸收光栅,位于聚光镜1出口的中心光阑2后,用于对来自聚光镜1和中心光阑2的空心锥光束进行分束,使其成为具有空间周期结构的空心锥光束。Beam splitting grating 3': It can be an absorption grating with a period of micron order, located behind the central diaphragm 2 at the exit of the condenser 1, and is used to split the hollow cone beam from the condenser 1 and the central diaphragm 2, making it a Hollow cone beam with spatially periodic structure.
针孔3:位于分束光栅3’后样品台前,用于阻挡来自所述X射线光源的过强的直通光以及杂散光照明样品和探测器。Pinhole 3: Located in front of the sample stage behind the beam splitting grating 3', it is used to block the excessively strong through light and stray light from the X-ray light source to illuminate the sample and detector.
样品台4:位于针孔3和物镜5之间,用于承载样品、平移和旋转样品,其所在位置垂直于光轴的平面称为物面。Sample stage 4: Located between the pinhole 3 and the objective lens 5, it is used to carry the sample, translate and rotate the sample, and the plane perpendicular to the optical axis is called the object plane.
物镜5:可以为具有X射线透镜成像功能的光学元件,例如波带片,是一块周期随半径增大逐渐变小的圆形光栅,对单色X射线具有透镜功能,因此被称为X射线透镜;其作用有二,第一个作用是把位于物面具有纳米量级的样品结构放大成像于像面上,形成具有微米量级结构的像,第二个作用是把分束光栅上被空心锥光束照明的环形部分成像于物镜5后焦面附近形成环形光栅像及环形光栅像光束。Objective lens 5: It can be an optical element with X-ray lens imaging function, such as a zone plate, which is a circular grating whose period gradually becomes smaller as the radius increases. It has a lens function for monochromatic X-rays, so it is called X-ray Lens; it has two functions. The first function is to magnify and image the nanometer-scale sample structure on the object surface on the image plane to form an image with a micron-scale structure. The annular portion illuminated by the hollow cone beam is imaged near the rear focal plane of the objective lens 5 to form an annular grating image and an annular grating image beam.
环形分析光栅6’:可以为周期为微米量级或者亚微米量级的吸收光栅,位于物镜5后焦面附近,其形状和尺寸与无样品时分束光栅3’被空心锥光束照明的环形部分在物镜5后焦面附近的环形光栅像相同,用于对所述物镜5后焦面附近的环形光栅像光束进行滤波。Annular analysis grating 6': It can be an absorption grating with a period of micron or submicron level, located near the rear focal plane of the objective lens 5, and its shape and size are the same as the annular part of the beam splitting grating 3' illuminated by the hollow cone beam when there is no sample The annular grating image near the rear focal plane of the objective lens 5 is the same, and is used to filter the annular grating image beam near the rear focal plane of the objective lens 5 .
成像探测器7:由二维面阵像素排列而成(例如X射线CCD),每个像素都具有独立探测光强的功能,用于测量角度信号响应曲线、探测光强空间位置变化、拍摄样品的放大像,其所在位置垂直于光轴的平面称为像面。Imaging detector 7: Arranged by two-dimensional area array pixels (such as X-ray CCD), each pixel has the function of independently detecting light intensity, which is used to measure the angle signal response curve, detect the spatial position change of light intensity, and photograph samples The magnified image of , the plane where its position is perpendicular to the optical axis is called the image plane.
进一步地,分束光栅3’和环形分析光栅6’虽然都是吸收光栅,但是在光栅厚度不大的情况下,会有部分X射线穿过光栅条,使穿过的X射线产生相移,为了避免相移带来的负面效应,分束光栅3’的厚度为π相移的奇数倍,环形分析光栅6’的厚度为π相移的偶数倍。Further, although the beam splitting grating 3' and the annular analysis grating 6' are both absorption gratings, when the thickness of the grating is not large, some X-rays will pass through the grating strips, causing the passing X-rays to generate a phase shift, In order to avoid the negative effect brought by the phase shift, the thickness of the beam splitting grating 3' is an odd multiple of the π phase shift, and the thickness of the annular analysis grating 6' is an even multiple of the π phase shift.
另外,分束光栅3’和环形分析光栅6’的材质为金或其它重金属。In addition, the beam splitting grating 3' and the ring analysis grating 6' are made of gold or other heavy metals.
本发明的实施例的X射线微分相位衬度显微镜系统是传统的X射线显微镜与分束光栅3’和环形分析光栅6’结合而成,在传统的X射线显微镜系统的成像过程基础上,增加了两个方面的物理内容:第一,考虑样品对焦深细光束的折射和散射影响,即考虑样品产生的角度信号对成像过程的影响;第二,考虑样品的角度信号调控分束光栅的环形光栅像在环形分析光栅上的剪切位移,使角度信号获得光强响应。The X-ray differential phase contrast microscope system of the embodiment of the present invention is a traditional X-ray microscope combined with a beam splitting grating 3' and an annular analysis grating 6'. On the basis of the imaging process of the traditional X-ray microscope system, an additional Two aspects of physical content are discussed: first, consider the influence of refraction and scattering of the narrow beam at the focal depth of the sample, that is, consider the influence of the angle signal generated by the sample on the imaging process; second, consider the angle signal of the sample to regulate the annular The shear displacement of the grating image on the annular analysis grating makes the angle signal obtain the light intensity response.
图3为根据本发明的一实施例的X射线微分相位衬度显微镜系统的二维成像方法的流程框图,如图3所示,用于所述的X射线微分相位衬度显微镜系统的二维成像方法,包括如下步骤:Fig. 3 is a flow chart of a two-dimensional imaging method of an X-ray differential phase contrast microscope system according to an embodiment of the present invention. As shown in Fig. 3, the two-dimensional An imaging method, comprising the steps of:
S1,开启并调整光源:使所述光源产生的X射线光束对准聚光镜;S1, turn on and adjust the light source: aim the X-ray beam generated by the light source at the condenser;
S2,调整聚光镜、中心光阑和针孔:使聚光镜和中心光阑产生的空心锥光束对准所述样品台上承载样品位置,使中心光阑和针孔阻挡来自所述X射线光源的直通光和杂散光;S2, adjust the condenser, central diaphragm and pinhole: make the hollow cone beam generated by the condenser and the central diaphragm align with the position of the sample on the sample stage, so that the central diaphragm and the pinhole block the direct passage from the X-ray light source light and stray light;
S3,调整物镜:使物镜聚焦形成的成像光束对准位于像面的成像探测器;S3, adjusting the objective lens: aligning the imaging light beam formed by the objective lens with the imaging detector located on the image plane;
S4,馈入分束光栅:在聚光镜出口的中心光阑后,馈入分束光栅,使空心锥光束照明分束光栅的环形部分在物镜后焦面附近形成环形光栅像及环形光栅像光束;S4, feeding into the beam splitting grating: after the central diaphragm at the exit of the condenser, feed into the beam splitting grating, so that the hollow cone beam illuminates the annular part of the beam splitting grating to form a ring grating image and a ring grating image beam near the rear focal plane of the objective lens;
S5,馈入并调整环形分析光栅:使环形分析光栅对准分束光栅在物镜后焦面附近形成的环形光栅像,调整环形分析光栅栅条方向,使环形分析光栅栅条平行于分束光栅栅条;S5, feed and adjust the ring analysis grating: make the ring analysis grating align with the ring grating image formed by the beam splitting grating near the rear focal plane of the objective lens, adjust the direction of the ring analysis grating bars so that the ring analysis grating bars are parallel to the beam splitting grating grid;
S6,测角度信号响应曲线:沿着垂直于光轴和光栅栅条方向逐步移动分束光栅,使分束光栅的环形光栅像相对环形分析光栅发生剪切位移,用成像探测器在像面测得亮度随环形光栅像位移变化的角度信号响应曲线;其中,成像探测器各个像素上测得的角度信号响应曲线基本相同,测量角度信号响应曲线时,探测器获得的光强信号越强,分束光栅移动步长越密集,测得的数据信噪比越高;当分束光栅的环形光栅像和环形分析光栅完全重合时,环形光栅像光束通过率最低,整个视场成为暗场,与角度信号响应曲线的最低位置—谷位相对应;当分束光栅的环形光栅像和环形分析光栅错开四分之一周期时,环形光栅像光束一半通过,另一半被阻断,整个视场成为半亮场,与角度信号响应曲线的上坡位相对应;当分束光栅的环形光栅像和环形分析光栅错开二分之一周期时,环形光栅像光束通过率最高,整个视场成为亮场,与角度信号响应曲线的最高位置—峰位相对应;当分束光栅的环形光栅像和环形分析光栅错开四分之三周期时,环形光栅像光束一半被阻断,另一半通过,整个视场成为半亮场,与角度信号响应移曲线的下坡位相对应;S6, angle measurement signal response curve: move the beam splitting grating step by step along the direction perpendicular to the optical axis and the grating bars, so that the annular grating image of the beam splitting grating has a shear displacement relative to the annular analysis grating, and the imaging detector is used to measure the angle on the image plane. The angle signal response curve of the brightness changing with the image displacement of the ring grating is obtained; the angle signal response curves measured on each pixel of the imaging detector are basically the same, and when the angle signal response curve is measured, the stronger the light intensity signal obtained by the detector, the more The denser the moving step of the beam grating, the higher the signal-to-noise ratio of the measured data; when the annular grating image of the beam splitting grating and the annular analysis grating completely overlap, the beam passing rate of the annular grating image is the lowest, and the entire field of view becomes a dark field. The lowest position of the signal response curve corresponds to the valley position; when the annular grating image of the beam splitting grating and the annular analysis grating are staggered by a quarter period, half of the annular grating image beam passes through, and the other half is blocked, and the entire field of view becomes a semi-bright field , corresponding to the uphill position of the angle signal response curve; when the annular grating image of the beam splitting grating and the annular analysis grating are staggered by one-half period, the beam passing rate of the annular grating image is the highest, and the entire field of view becomes a bright field, which is consistent with the angular signal response The highest position of the curve corresponds to the peak position; when the annular grating image of the beam splitting grating and the annular analysis grating are staggered by three quarters of the period, half of the beam of the annular grating image is blocked, and the other half passes through, and the entire field of view becomes a semi-bright field, which is the same as The angle signal corresponds to the downhill position of the shift curve;
S7,用余弦曲线拟合角度信号响应曲线:角度信号响应曲线类似于余弦曲线,用余弦曲线拟合测得的角度信号响应曲线,使其获得余弦曲线的解析表达;其中,拟合角度信号响应曲线时,探测器获得的光强信号越强,分束光栅移动步长越密集,测得的数据信噪比越高,拟合的余弦曲线越精确;S7, fitting the angle signal response curve with a cosine curve: the angle signal response curve is similar to a cosine curve, and the measured angle signal response curve is fitted with a cosine curve to obtain an analytical expression of a cosine curve; wherein, the fitting angle signal response When the curve is drawn, the stronger the light intensity signal obtained by the detector, the denser the moving step of the beam splitting grating, the higher the signal-to-noise ratio of the measured data, and the more accurate the fitted cosine curve;
S8,拍摄样品二维放大像:把分束光栅的环形光栅像分别固定在所述角度信号响应曲线的谷位、上坡位、峰位、下坡位,在样品台上放置样品,拍摄样品的谷位放大像、上坡放大像、峰位放大像、下坡放大像;S8, take a two-dimensional enlarged image of the sample: fix the annular grating image of the beam splitting grating at the valley position, uphill position, peak position, and downhill position of the angle signal response curve respectively, place the sample on the sample stage, and photograph the sample The valley zoom image, uphill zoom image, peak zoom image, downhill zoom image;
S9,提取样品在物空间的定量二维像:从样品的谷位放大像、上坡放大像、峰位放大像、下坡放大像中提取样品在物空间的吸收像、折射像和散射像;因为谷位放大像与峰位放大像之和等于上坡放大像与下坡放大像之和,所以上述四种放大像中只有三种放大像是独立的,可以从四种放大像中的任意三种放大像中提取样品在物空间的吸收像、折射像和散射像。S9, extract the quantitative two-dimensional image of the sample in the object space: extract the absorption image, refraction image and scattering image of the sample in the object space from the valley magnification image, uphill magnification image, peak magnification image, and downhill magnification image of the sample ; Because the sum of the valley magnification and the peak magnification is equal to the sum of the uphill magnification and the downhill magnification, only three of the above four magnifications are independent, and can be obtained from the four magnifications. Extract the absorption image, refraction image and scattering image of the sample in the object space from any three magnified images.
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S6中,考虑到样品对焦深细光束的的折射和散射影响,即考虑到样品产生的角度信号对成像过程的影响,因此,需要建立样品对X射线作用的数学模型,求出样品用角度信号函数的数学表达式,具体方法如下所述:In the step S6 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the influence of refraction and scattering of the narrow beam in the focal depth of the sample is considered, that is, the influence of the angle signal generated by the sample on the imaging process is considered, so , it is necessary to establish a mathematical model of the sample’s action on X-rays, and obtain the mathematical expression of the sample’s angle signal function. The specific method is as follows:
定义样品的直角坐标系(xo,yo,zo),对样品中一点进行定义,在二维成像中,样品所在物面上一点(xo,yo),不是一个二维几何点,而是一个以(xo,yo)为中心的物面积元ΔxΔy,Δx和Δy的大小由物镜的数值孔径和探测器分辨率决定。在此特别说明,下面提到的物点,其含义均为物面积元。Define the Cartesian coordinate system (x o , y o , z o ) of the sample, and define a point in the sample. In two-dimensional imaging, a point (x o , y o ) on the object surface where the sample is located is not a two-dimensional geometric point , but an object area element ΔxΔy centered on (x o , y o ), the size of Δx and Δy is determined by the numerical aperture of the objective lens and the resolution of the detector. It is specially stated here that the object point mentioned below means the object area element.
吸收(包括非弹性散射)是一个X射线能量在样品中转化为热能的耗散过程,样品中一点(xo,yo)对焦深细光束吸收,引起焦深细光束亮度衰减。样品中一点(xo,yo)对通过该点的焦深细光束的吸收可以表达为:Absorption (including inelastic scattering) is a dissipation process in which X-ray energy is converted into heat energy in the sample. A point (x o , y o ) in the sample absorbs the focal depth thin beam, causing the brightness of the focal depth thin beam to attenuate. The absorption of a point (x o , y o ) in the sample to the focal depth thin beam passing through this point can be expressed as:
其中表示狄拉克函数(为避免和折射率实部衰减率δ混淆,用顶部带标记的表示狄拉克函数),为入射在样品上的焦深细光束的偏转角矢量,in Represents the Dirac function (to avoid confusion with the real part of the refractive index decay rate δ, marked with the top band represents the Dirac function), is the deflection angle vector of the depth-of-focus thin beam incident on the sample,
其中μ(xo,yo,zo)为样品的线性吸收系数。式(1)的物理意义为,吸收引起焦深细光束能量损失,导致焦深细光束亮度降低,但不改变焦深细光束传播方向,是零角度信号。式(1)还可以表示为分量形式:Where μ(x o , y o , z o ) is the linear absorption coefficient of the sample. The physical meaning of formula (1) is that the energy loss of the depth-of-focus thin beam is caused by absorption, which leads to the decrease of the brightness of the depth-of-focus thin beam, but does not change the propagation direction of the depth-of-focus thin beam, which is a zero-angle signal. Equation (1) can also be expressed in component form:
其中ψx和ψy分别为沿xo方向和yo方向的分量。where ψ x and ψ y are respectively Components along the x o direction and y o direction.
图4为根据本发明的实施例的X射线微分相位衬度显微镜系统的折射角信号成像原理图,局部放大图中的箭头描绘焦深细光束的折射方向。如图4所示,来自分束光栅的聚焦于样品的空心锥光束,在物面可以分割为多个焦深细光束,每一个焦深细光束照明一个物点(xo,yo),物点(xo,yo)对经过该点的焦深细光束产生折射,引起物镜后焦面附近的环形光栅像相对环形分析光栅发生位移,使通过环形分析光栅到达像面像点(xi,yi)的光子数发生增加或者减少,(xo,yo)和(xi,yi)之间是共轭成像关系,相互唯一对应。Fig. 4 is a schematic diagram of refraction angle signal imaging of the X-ray differential phase contrast microscope system according to an embodiment of the present invention, and the arrows in the partially enlarged diagram depict the refraction direction of the depth-of-focus thin beam. As shown in Figure 4, the hollow cone beam focused on the sample from the beam splitting grating can be split into multiple depth-of-focus beamlets on the object plane, and each depth-of-focus beamlet illuminates an object point (x o , y o ), The object point (x o , y o ) refracts the focal depth thin beam passing through this point, causing the annular grating image near the rear focal plane of the objective lens to shift relative to the annular analysis grating, so that the image point (x The number of photons of i , y i ) increases or decreases, and the relationship between (x o , y o ) and ( xi , y i ) is a conjugate imaging relationship, which uniquely corresponds to each other.
折射是一个能量守恒的过程,根据图4,样品中一点(xo,yo)对焦深细光束的折射可以表达为Refraction is a process of energy conservation. According to Fig. 4, the refraction of a point (x o , y o ) in the sample at the focal depth of the thin beam can be expressed as
其中为折射角矢量,式(4)的物理意义为,折射改变焦深细光束的传播方向,但不损失焦深细光束能量,不导致焦深细光束亮度降低。式(4)还可以写为分量形式,in is the refraction angle vector, and the physical meaning of formula (4) is that refraction changes the propagation direction of the depth-of-focus thin beam, but does not lose the energy of the depth-of-focus thin beam, and does not cause the brightness of the depth-of-focus thin beam to decrease. Equation (4) can also be written in component form,
其中θx(xo,yo)和θy(xo,yo)分别为沿x方向和y方向的分量。令δ(xo,yo,zo)代表样品折射率实部衰减率,则在(xo,yo)直角坐标系的分量表达式为Where θ x (x o , y o ) and θ y (x o , y o ) are respectively Components along the x-direction and y-direction. Let δ(x o ,y o ,z o ) represent the attenuation rate of the real part of the refractive index of the sample, then The component expression in (x o , y o ) Cartesian coordinate system is
图5为根据本发明的实施例的X射线微分相位衬度显微镜系统的散射角信号成像原理图,局部放大图中箭头描绘的是散射引起焦深细光束从一个传播方向扩展为多个方向。如图5所示,来自分束光栅的聚焦于样品的空心锥光束,在物面可以分割为多个焦深细光束,每一个焦深细光束照明一个物点(xo,yo),物点(xo,yo)对经过该点的焦深细光束产生散射,引起物镜后焦面附近的环形光栅像的条纹模糊,使通过环形分析光栅到达像面像点(xi,yi)的光子数发生增加或者减少,(xo,yo)和(xi,yi)之间是共轭成像关系,相互唯一对应。FIG. 5 is a schematic diagram of scattering angle signal imaging of an X-ray differential phase contrast microscope system according to an embodiment of the present invention. The arrows in the partially enlarged figure depict that scattering causes the focal depth thin beam to expand from one propagation direction to multiple directions. As shown in Figure 5, the hollow cone beam focused on the sample from the beam splitting grating can be split into multiple depth-of-focus beamlets on the object plane, and each depth-of-focus beamlet illuminates an object point (x o , y o ), The object point (x o , y o ) scatters the focal depth thin beam passing through this point, which causes the fringes of the annular grating image near the back focal plane of the objective lens to be blurred, so that the image point (x i , y The number of photons in i ) increases or decreases, and the relationship between (x o , y o ) and ( xi , y i ) is a conjugate imaging relationship, which uniquely corresponds to each other.
散射(在此指面积元内部小颗粒引起的折射)是一个能量守恒的过程,散射和折射的不同之处在于,折射把样品物面上一个面积元作为一个整体来研究,即把样品物面上一个面积元作为一块棱镜,而散射则把这个面积元作为一块磨砂玻璃,研究其内部的不均匀性质,例如面积元内部的颗粒、气泡、微晶、杂质等。因此,对于每个面积元,出射的焦深细光束只有一个折射方向,却有多个散射方向。换言之,散射是一个焦深细光束分散的过程。因为样品有一定厚度,在面积元内部沿着焦深细光束传播方向,各小颗粒分布是随机的,前后两个小颗粒产生的折射是相互独立的,小颗粒每次折射偏离焦深细光束入射方向的角度是随机的,所以根据中心极限定理,散射角是以入射角为中心的正态统计分布,可以用方差来描述散射角分布范围。Scattering (referring to the refraction caused by small particles inside the area element) is a process of energy conservation. The difference between scattering and refraction is that refraction studies an area element on the sample surface as a whole, that is, the sample surface The last area element is used as a prism, and the scattering uses this area element as a piece of frosted glass to study the inhomogeneous properties inside it, such as particles, bubbles, microcrystals, impurities, etc. inside the area element. Therefore, for each area element, the exiting depth-of-focus beamlet has only one refraction direction, but multiple scattering directions. In other words, scattering is a process in which the focal depth of a thin beam is dispersed. Because the sample has a certain thickness, along the propagation direction of the focal depth thin beam inside the area element, the distribution of each small particle is random, and the refraction produced by the two small particles before and after is independent of each other, and each refraction of the small particle deviates from the depth of focus thin beam The angle of the incident direction is random, so according to the central limit theorem, the scattering angle is a normal statistical distribution centered on the incident angle, and the variance can be used to describe the distribution range of the scattering angle.
根据图5,焦深细光束射入样品时,随着焦深细光束在样品中穿行,面积元内部小颗粒折射事件的不断发生,散射角方差逐渐加宽。因为利用光栅采集样品的角度信号,所以有必要把散射角信号进行分解。样品一点(xo,yo)对焦深细光束的散射在(xo,yo)直角坐标系的分量表达式为According to Figure 5, when the depth-of-focus thin beam enters the sample, as the depth-of-focus thin beam travels through the sample, the refraction events of small particles inside the area element continue to occur, and the variance of the scattering angle gradually widens. Since the angle signal of the sample is collected by a grating, it is necessary to decompose the scattering angle signal. The component expression of the scattering of the focused deep thin beam at the sample point (x o , y o ) in the (x o , y o ) Cartesian coordinate system is
其中和为(xo,yo)点处样品整体厚度分别在xo方向和yo方向产生的散射角方差。根据式(7)和式(8),虽然散射不损失光束能量,但是散射引起焦深细光束扩散角增大。令代表样品整体厚度的散射角方差,下标κ可以是x或者y,则是X射线经过路径上一系列Δzo薄片的微分散射角方差之和,所以样品整体厚度的散射角方差可以表示为微分散射角方差的积分,即:in and is the variance of the scattering angle generated by the overall thickness of the sample at point (x o , y o ) in the x o direction and y o direction, respectively. According to Equation (7) and Equation (8), although scattering does not lose beam energy, scattering causes an increase in the spread angle of the narrow beam at the depth of focus. make represents the scattering angle variance of the overall thickness of the sample, the subscript κ can be x or y, then is the differential scattering angle variance of a series of Δz o slices along the X-ray path The sum of , so the variance of the scattering angle across the thickness of the sample can be expressed as the integral of the variance of the differential scattering angle, namely:
其中ωκ(xo,yo,zo)为垂直于光束沿κ方向的线性散射系数。where ω κ (x o , y o , z o ) is the linear scattering coefficient perpendicular to the beam along the κ direction.
综合考虑上述三种作用,样品中一点(xo,yo)对通过该点的焦深细光束的作用可以用样品的角度信号函数f(xo,yo,ψκ)表达为:Considering the above three effects comprehensively, the effect of a point (x o , y o ) in the sample on the focal depth thin beam passing through this point can be expressed as :
其中下标κ可以是x或者y。where the subscript κ can be x or y.
根据式(10),可知在垂直于光束沿着κ方向,出射X射线携带了样品的三种角度信号,零角度信号:M(xo,yo)、折射角信号:θκ(xo,yo)、散射角信号:它们均可以表示为线积分,这就为计算机断层成像利用投影数据获得样品三维结构奠定了数学基础。According to formula (10), it can be seen that along the κ direction perpendicular to the beam, the outgoing X-rays carry three angle signals of the sample, zero angle signal: M(x o , y o ), refraction angle signal: θ κ (x o ,y o ), scattering angle signal: All of them can be expressed as line integrals, which lays a mathematical foundation for computer tomography to use projection data to obtain the three-dimensional structure of samples.
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S6中,考虑到样品的角度信号调控分束光栅的环形光栅像在环形分析光栅上的剪切位移,使角度信号获得光强响应。因此,需要建立X射线微分相位衬度显微镜系统的角度信号响应函数,具体方法如下所述:In the step S6 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, considering the angle signal of the sample to regulate the shear displacement of the ring grating image of the beam splitting grating on the ring analysis grating, the angle signal is obtained by optical Strong response. Therefore, it is necessary to establish the angle signal response function of the X-ray differential phase contrast microscope system, and the specific method is as follows:
在传统的X射线显微镜中,无样品时的光束传播过程可以描述如下。来自X射线光源的扩散光束经过聚光镜会聚和中心光阑阻挡,形成聚焦于样品的空心锥光束,可以在物面按照空间相干面积,把聚焦的空心锥光束分割为多个焦深细光束,每一个焦深细光束照明一个物点位置。在焦深细光束内部是相干光,而两相邻焦深细光束之间互不相干。经过物面后,各焦深细光束又变成一个个扩散的相干空心锥光束,经过物镜的聚焦作用,各扩散的相干空心锥光束又变成一个个聚焦的相干空心锥光束,在像面上形成一个个像点。在传播过程中,在相干空心锥光束内部是相干的,而两相干空心锥光束之间是不相干的。在聚焦时,各相干空心锥光束都形成焦深细光束、在物面和像面上互相分离,可是在聚焦前后,例如在聚光镜、中心光阑、物镜、物镜后焦面等位置,这些相干空心锥光束在空间上是互相重合的。在传统的X射线显微镜成像理论中,只考虑样品的吸收和相移,不考虑样品的折射和散射,或者即使样品存在折射和散射,可是不存在把折射和散射转变为光强的响应机制,传统的X射线显微镜成像理论认为,放入样品后,除了样品吸收引起光束亮度下降和相移外,光束传播过程和无样品时没有不同。In a conventional X-ray microscope, the beam propagation process without a sample can be described as follows. The diffuse beam from the X-ray source is converged by the condenser and blocked by the central diaphragm to form a hollow cone beam focused on the sample. The focused hollow cone beam can be divided into multiple focal depth thin beams on the object surface according to the spatial coherence area. A depth-of-focus thin beam illuminates an object point position. Inside the depth-of-focus thin beams is coherent light, while two adjacent depth-of-focus thin beams are incoherent. After passing through the object plane, each depth-of-focus thin beam becomes a diffused coherent hollow cone beam. After the focusing effect of the objective lens, each diffused coherent hollow cone beam becomes a focused coherent hollow cone beam. Each image point is formed. In the process of propagation, there is coherence inside the coherent hollow cone beams, but incoherence between two coherent hollow cone beams. When focusing, each coherent hollow cone beam forms a thin beam of focal depth and is separated from each other on the object plane and image plane, but before and after focusing, such as in the position of the condenser lens, the central diaphragm, the objective lens, and the back focal plane of the objective lens, these coherent Hollow cone beams are mutually coincident in space. In the traditional X-ray microscope imaging theory, only the absorption and phase shift of the sample are considered, and the refraction and scattering of the sample are not considered, or even if there is refraction and scattering in the sample, there is no response mechanism that converts refraction and scattering into light intensity. The traditional X-ray microscope imaging theory believes that after the sample is placed, the beam propagation process is the same as when there is no sample, except that the beam brightness decreases and the phase shift is caused by the absorption of the sample.
在传统的X射线显微镜系统的聚光镜出口、中心光阑后插入一块分束光栅,在物镜后焦面附近插入一块环形分析光栅,就形成了本发明的一实施例中的X射线微分相位衬度显微镜。分束光栅的作用是,把空心锥光束分割为具有空间结构的空心锥光束,并使分束光栅被空心锥光束照明的环形部分在物镜后焦面附近形成环形光栅像及环形光栅像光束。环形分析光栅的形状和尺寸与分束光栅被空心锥光束照明的环形部分在物镜后焦面附近形成环形光栅像完全相同,其作用是对所述物镜后焦面附近的环形光栅像光束进行滤波。由此可知,在X射线微分相位衬度显微镜中,无样品时的光束传播过程可以描述如下:来自聚光镜和中心光阑的空心锥光束,经过分束光栅分束,形成具有空间结构的聚焦于样品的空心锥光束。在物面按照空间相干面积,把聚焦的空心锥光束分割为多个焦深细光束,每一个焦深细光束照明一个物点。在焦深细光束内部是相干光,而两相邻焦深细光束之间互不相干。经过物面后,各焦深细光束又变成一个个扩散的相干空心锥光束,经过物镜的聚焦作用,各扩散的相干空心锥光束又变成一个个聚焦的相干空心锥光束,不仅在物镜后焦面附近形成各自的环形光栅像光束,而且在像面上形成一个个像点。虽然各相干空心锥光束在物镜后焦面附近形成的环形光栅像光束互相重合,但是它们各自独立传播,互不相干。A beam splitting grating is inserted behind the exit of the condenser lens and the central diaphragm of the traditional X-ray microscope system, and an annular analysis grating is inserted near the rear focal plane of the objective lens, forming the X-ray differential phase contrast in an embodiment of the present invention microscope. The function of the beam splitting grating is to split the hollow cone beam into hollow cone beams with spatial structure, and make the annular part of the beam splitting grating illuminated by the hollow cone beam form a ring grating image and a ring grating image beam near the rear focal plane of the objective lens. The shape and size of the annular analysis grating are exactly the same as that of the annular part of the beam splitting grating illuminated by the hollow cone beam to form an annular grating image near the rear focal plane of the objective lens, and its function is to filter the annular grating image beam near the rear focal plane of the objective lens . It can be seen that in the X-ray differential phase contrast microscope, the beam propagation process when there is no sample can be described as follows: the hollow cone beam from the condenser and the central diaphragm is split by the beam splitting grating to form a beam with a spatial structure focused on Hollow cone beam of the sample. According to the spatial coherence area on the object plane, the focused hollow cone beam is divided into multiple depth-of-focus thin beams, and each depth-of-focus thin beam illuminates an object point. Inside the depth-of-focus thin beams is coherent light, while two adjacent depth-of-focus thin beams are incoherent. After passing through the object plane, each depth-of-focus thin beam becomes a diffused coherent hollow cone beam. After the focusing effect of the objective lens, each diffused coherent hollow cone beam becomes a focused coherent hollow cone beam, not only in the objective lens Respective annular grating image beams are formed near the back focal plane, and image points are formed on the image plane. Although the annular grating image beams formed by the coherent hollow cone beams near the back focal plane of the objective lens coincide with each other, they propagate independently and are independent of each other.
若在X射线微分相位衬度显微镜的物面上放置样品,则样品中一点对通过该点的焦深细光束产生折射和散射作用,使该焦深细光束产生偏转和发散的两种角度信号,这两种角度信号必然会引起该焦深细光束的环形光栅像在环形分析光栅上发生位置移动和条纹模糊,引起通过环形分析光栅的光子数增加或者减少,样品中一点的角度信号就会在像面的像点上获得光强响应。因此,环形分析光栅对分束光栅的环形光栅像的滤波作用,为X射线微分相位衬度显微镜提供了角度信号响应成像机制。因为样品折射引起分束光栅的环形光栅像的位置移动,和无样品时人为移动分束光栅等价,所以在无样品时,可以人为移动分束光栅,利用成像探测器一个像素测得环形光栅像和环形分析光栅之间剪切位移产生的光强响应,从而测得X射线微分相位衬度显微镜对物面上一点的角度信号响应函数。至于样品散射引起分束光栅的环形光栅像产生条纹模糊,可以看作环形光栅像各局部相对环形分析光栅做无规律的剪切位移,由此产生的光强响应也可以用测得的角度信号响应函数进行解释。因为各物点的环形光栅像在物镜后焦面附近是重合的,所以可以用一块分束光栅和一块环形分析光栅就能把物面上各物点的角度信号转变为像面上各像点的光强响应。换言之,样品各点的折射角信号和散射角信号转换为响应光强的机制均相同,具有统一的角度信号响应函数。因此在无样品时,可以人为步进移动分束光栅一次,利用成像探测器各个像素可以同时并行测得各自的角度信号响应函数。If the sample is placed on the object surface of the X-ray differential phase contrast microscope, a point in the sample will refract and scatter the focal depth thin beam passing through this point, causing the focal depth thin beam to produce two angle signals of deflection and divergence , these two kinds of angle signals will inevitably cause the position shift and fringe blurring of the annular grating image of the depth-of-focus thin beam on the annular analysis grating, which will cause the number of photons passing through the annular analysis grating to increase or decrease, and the angle signal of a point in the sample will be The light intensity response is obtained on the image points of the image plane. Therefore, the filtering effect of the ring analysis grating on the ring grating image of the beam splitting grating provides an angle signal response imaging mechanism for the X-ray differential phase contrast microscope. Because the refraction of the sample causes the position of the annular grating image of the beam-splitting grating to move, which is equivalent to moving the beam-splitting grating artificially when there is no sample, so when there is no sample, the beam-splitting grating can be artificially moved, and the annular grating can be measured with one pixel of the imaging detector The light intensity response generated by the shear displacement between the image and the annular analysis grating is used to measure the angular signal response function of the X-ray differential phase contrast microscope to a point on the object surface. As for the fringe blurring of the ring grating image of the beam splitting grating caused by sample scattering, it can be regarded as the irregular shear displacement of each part of the ring grating image relative to the ring analysis grating, and the resulting light intensity response can also be measured by the angle signal The response function is explained. Because the ring grating images of each object point are coincident near the rear focal plane of the objective lens, a beam splitting grating and a ring analysis grating can be used to convert the angle signal of each object point on the object plane into each image point on the image plane light intensity response. In other words, the mechanism of converting the refraction angle signal and scattering angle signal at each point of the sample into the response light intensity is the same, and has a unified angle signal response function. Therefore, when there is no sample, the beam-splitting grating can be manually stepped once, and each pixel of the imaging detector can be used to measure the respective angular signal response functions in parallel.
由于分束光栅和环形分析光栅仅能对垂直光栅栅条方向角度信号产生光强响应,而对平行于光栅栅条方向的角度信号没有响应,所以只须考虑垂直于光栅条方向的角度信号函数。根据分束光栅的环形光栅像和环形分析光栅的剪切位移,可知测得的角度信号响应曲线为一周期振荡曲线。由于折射和散射都是小角度信号,所以不是整条角度信号响应曲线在起作用,而是角度信号响应曲线局部在起作用。因此,根据分束光栅步进扫描起点不同,可以把角度信号响应函数曲线可以分为四种类型:Since the beam splitting grating and the annular analysis grating can only produce light intensity responses to the angle signal perpendicular to the direction of the grating bars, but have no response to the angle signal parallel to the direction of the grating bars, so only the angle signal function perpendicular to the direction of the grating bars needs to be considered . According to the ring grating image of the beam splitting grating and the shear displacement of the ring analysis grating, it can be known that the measured angle signal response curve is a periodic oscillation curve. Since both refraction and scattering are small-angle signals, it is not the entire angle signal response curve that is at work, but the angle signal response curve is locally at work. Therefore, according to the different start points of the beam-splitting grating step scan, the angle signal response function curves can be divided into four types:
(i)以分束光栅的环形光栅像和环形分析光栅完全重合为零点,随着分束光栅步进扫描,分束光栅的环形光栅像相对环形分析光栅发生剪切位移,测得亮度逐步上升的曲线,称为谷位响应曲线;(i) The ring-shaped grating image of the beam-splitting grating and the ring-shaped analysis grating completely coincide as the zero point. With the step-by-step scanning of the beam-splitting grating, the ring-shaped grating image of the beam-splitting grating undergoes a shear displacement relative to the ring analysis grating, and the measured brightness gradually increases The curve is called the valley response curve;
(ii)以分束光栅的环形光栅像和环形分析光栅错开四分之一周期为零点,随着分束光栅步进扫描,分束光栅的环形光栅像相对环形分析光栅发生剪切位移,沿着剪切位移增大方向,测得亮度逐步线性增加曲线,沿着剪切位移减小的方向,测得亮度逐步线性下降曲线,称为上坡响应曲线;(ii) Taking the staggered quarter period between the annular grating image of the beam splitting grating and the annular analysis grating as the zero point, along with the step scanning of the beam splitting grating, the shear displacement of the annular grating image of the beam splitting grating relative to the annular analysis grating occurs along the Along the direction of increasing shear displacement, the measured brightness gradually increases linearly, and along the direction of shear displacement decreases, the measured brightness gradually decreases linearly, which is called the uphill response curve;
(iii)以分束光栅的环形光栅像和环形分析光栅错开二分之一周期为零点,随着分束光栅步进扫描,分束光栅的环形光栅像相对环形分析光栅发生剪切位移,测得亮度逐步下降的曲线,称为峰位响应曲线;(iii) The ring grating image of the beam splitting grating and the ring analysis grating are staggered by half of the period as the zero point. With the step scanning of the beam splitting grating, the ring grating image of the beam splitting grating has a shear displacement relative to the ring analysis grating. The curve that the brightness gradually decreases is called the peak response curve;
(iv)以分束光栅的环形光栅像和环形分析光栅错开四分之三周期为零点,随着分束光栅步进扫描,分束光栅的环形光栅像相对环形分析光栅发生剪切位移,沿着剪切位移增大方向,测得亮度逐步线性下降曲线,沿着剪切位移减小的方向,测得亮度逐步线性上升曲线,称为下坡响应曲线。(iv) Taking the staggered three-quarters cycle of the annular grating image of the beam splitting grating and the annular analysis grating as the zero point, along with the step scanning of the beam splitting grating, the shear displacement of the annular grating image of the beam splitting grating relative to the annular analysis grating occurs along the Along the direction of shear displacement increase, the measured brightness gradually decreases linearly, and along the direction of shear displacement decreases, the measured brightness gradually increases linearly, which is called the downhill response curve.
图6为X射线微分相位衬度显微镜系统的四种角度信号响应函数曲线图。根据图6,可知角度信号响应曲线类似于余弦曲线,可用余弦曲线拟合近似,使角度信号响应曲线可以用拟合的余弦曲线进行解析表达,有Fig. 6 is a curve diagram of four angle signal response functions of the X-ray differential phase contrast microscope system. According to Figure 6, it can be seen that the angle signal response curve is similar to a cosine curve, which can be approximated by fitting a cosine curve, so that the angle signal response curve can be analytically expressed by a fitted cosine curve, and
或者表达为or expressed as
其中xp为环形分析光栅处位置坐标,为分束光栅的环形光栅像位移相对物面的角位移,do为样品相对物镜的物距,p为环形分析光栅的周期,为角度信号响应曲线平均值,Rmax和Rmin分别为角度信号响应曲线的最大值和最小值,为角度信号响应曲线的可见度,n为调制参数,n=0,1,2,3分别对应于图6中的(a)谷位响应曲线、(b)上坡响应曲线、(c)峰位响应曲线和(d)下坡响应曲线。Where x p is the position coordinates of the annular analysis grating, is the angular displacement of the annular grating image displacement of the beam splitting grating relative to the object plane, d o is the object distance of the sample relative to the objective lens, p is the period of the annular analysis grating, is the average value of the angle signal response curve, R max and R min are the maximum and minimum values of the angle signal response curve, respectively, is the visibility of the angle signal response curve, n is the modulation parameter, and n=0, 1, 2, 3 respectively correspond to (a) valley response curve, (b) uphill response curve, (c) peak position in Figure 6 Response curves and (d) downhill response curves.
根据角度信号响应函数的推导过程,可知角度信号响应函数描述一个物点的角度信号调控该物点的成像光子数在经过环形分析光栅时的通过率。According to the derivation process of the angle signal response function, it can be seen that the angle signal response function describes the angle signal of an object point and regulates the pass rate of the imaging photon number of the object point when passing through the annular analysis grating.
根据亮度不变定理,在无能量衰减的理想光学系统中,亮度是一个守恒量。虽然真实光学系统的传输效率不可能达到100%,但是可以把真实光学系统在像面的亮度可以看作物面亮度和一个传输效率的乘积。在X射线微分相位衬度显微镜系统中,从物面到像面,除了样品吸收引起亮度下降外,引起亮度下降的主要贡献来自物镜的衍射效率和环形分析光栅的通过率,因而像面亮度可以看作物面的亮度和物镜衍射效率和环形分析光栅通过率的乘积。无样品时,设物面处入射光束亮度为B0,物镜的衍射效率为η,则在物镜后焦面附近的分束光栅环形光栅像的亮度为ηB0。根据式(12),无样品时,环形分析光栅的通过率就是角度信号响应函数,因而环形分析光栅后的亮度,或者像面亮度为According to the constant brightness theorem, in an ideal optical system without energy attenuation, brightness is a conserved quantity. Although the transmission efficiency of a real optical system cannot reach 100%, the brightness of a real optical system on the image plane can be regarded as the product of the plane brightness and a transmission efficiency. In the X-ray differential phase contrast microscope system, from the object plane to the image plane, in addition to the decrease in brightness caused by the absorption of the sample, the main contribution to the decrease in brightness comes from the diffraction efficiency of the objective lens and the pass rate of the annular analysis grating, so the brightness of the image plane can be Look at the brightness of the crop surface and the product of the diffraction efficiency of the objective lens and the pass rate of the annular analysis grating. When there is no sample, assuming that the brightness of the incident beam at the object plane is B 0 , and the diffraction efficiency of the objective lens is η, then the brightness of the beam-splitting ring grating image near the rear focal plane of the objective lens is ηB 0 . According to formula (12), when there is no sample, the pass rate of the annular analysis grating is the angle signal response function, so the brightness behind the annular analysis grating, or the brightness of the image plane is
式(13)就是考虑了物镜衍射效率的角度信号响应函数。Equation (13) is the angular signal response function considering the diffraction efficiency of the objective lens.
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S8中,成像探测器拍摄的二维放大像用角度信号成像方程描述。因此,需要建立角度信号成像方程,具体方法如下所述:In step S8 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the two-dimensional enlarged image captured by the imaging detector is described by an angle signal imaging equation. Therefore, it is necessary to establish an angle signal imaging equation, and the specific method is as follows:
对于样品放大成像而言,探测器各个像素对物面各点的成像是并行且相互独立的,因此,只须讨论探测器一个像素对样品中一点的成像,就能建立X射线微分相位衬度显微镜成像方程。因为样品各点折射和散射转换为光强的机制均相同,在物面上具有平移不变性,所以可以根据上述角度信号函数和角度信号响应函数的卷积推导出X射线微分相位衬度显微镜的角度信号成像方程。在此需要特别说明,卷积运算是角度信号函数中的折射角、散射角和角度信号响应函数卷积运算,卷积运算和空间坐标无关;角度信号函数中的吸收是零角度信号,不参与卷积运算。For sample magnification imaging, the imaging of each pixel of the detector to each point of the object plane is parallel and independent of each other. Therefore, the X-ray differential phase contrast can be established only by discussing the imaging of a pixel of the detector to a point in the sample. Microscopy imaging equation. Because the mechanism of refraction and scattering conversion into light intensity at each point of the sample is the same, and it has translation invariance on the object plane, the X-ray differential phase contrast microscope can be deduced according to the convolution of the above angle signal function and angle signal response function. Angle signal imaging equation. It needs to be specially explained here that the convolution operation is the convolution operation of the refraction angle, scattering angle and angle signal response function in the angle signal function, and the convolution operation has nothing to do with the spatial coordinates; the absorption in the angle signal function is a zero-angle signal and does not participate in convolution operation.
放入样品后,物面上一物点(xo,yo)对通过该点的焦深细光束的吸收会在物镜后焦面附近引起环形光栅像发生亮度下降,折射会在物镜后焦面附近引起环形光栅像发生位置偏移,散射会在物镜后焦面附近引起环形光栅像发生条纹模糊。经过环形分析光栅滤波后,像面上像点(xi,yi)的亮度是角度信号函数和角度信号响应函数的卷积。由于分束光栅和环形分析光栅仅能对垂直光栅栅条方向角度信号产生光强响应,而对平行于光栅栅条方向的角度信号没有响应,所以只须考虑垂直于光栅条方向的角度信号函数和角度信号响应函数的卷积。因而有After the sample is put in, the absorption of the focal depth thin beam passing through this point by an object point (x o , y o ) on the object surface will cause the brightness of the annular grating image to decrease near the rear focal plane of the objective lens, and the refraction will cause The position of the ring grating image will be shifted near the surface, and the scattering will cause the ring grating image to appear fringe blur near the back focal plane of the objective lens. After the ring analysis grating filter, the brightness of the image point ( xi , y i ) on the image plane is the convolution of the angle signal function and the angle signal response function. Since the beam splitting grating and the annular analysis grating can only produce light intensity responses to the angle signal perpendicular to the direction of the grating bars, but have no response to the angle signal parallel to the direction of the grating bars, so only the angle signal function perpendicular to the direction of the grating bars needs to be considered and the convolution of the angle signal response function. Therefore there is
其中(xo,yo)为物面坐标,(xi,yi)为像面坐标,B(xi,yi,ψx)为通过物点(xo,yo)的焦深细光束在像点(xi,yi)的亮度。式(17)即为X射线微分相位衬度显微镜建立的角度信号成像方程。Where (x o , y o ) is the coordinates of the object plane, ( xi , y i ) is the coordinates of the image plane, B( xi , y i ,ψ x ) is the depth of focus passing through the object point (x o , y o ) The brightness of the thin beam at the image point ( xi , yi ). Equation (17) is the angle signal imaging equation established by the X-ray differential phase contrast microscope.
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线的谷位,在样品台上放置样品,成像探测器拍摄到样品的谷位放大像,具体在实验中拍摄二维放大像的方法如下所述:In the step S8 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the annular grating image of the beam splitting grating is fixed at the valley position of the angle signal response curve, the sample is placed on the sample stage, and the imaging The detector captures the magnified image of the valley of the sample. The specific method of taking a two-dimensional magnified image in the experiment is as follows:
放入样品前,把环形光栅像固定在角度信号响应曲线的谷位,即令ψx=0,n=0,然后放置好样品,探测器可以拍摄到样品的谷位放大像,其表达式为Before placing the sample, the ring grating image is fixed at the valley of the angle signal response curve, that is, ψ x = 0, n = 0, and then the sample is placed, and the detector can capture the enlarged image of the valley of the sample, and its expression is
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线的上坡位,拍摄样品上坡位放大像的方法:放入样品前,把环形光栅像固定在角度信号响应曲线的上坡位,即令ψx=0,n=1,然后放置好样品,探测器可以拍摄到样品的上坡放大像,其表达式为In the step S8 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the annular grating image of the beam splitting grating is fixed at the uphill position of the angular signal response curve, and the enlarged image of the sample uphill position is taken The method: before placing the sample, fix the ring grating image at the uphill position of the angle signal response curve, that is, set ψ x = 0, n = 1, then place the sample, the detector can take the uphill magnified image of the sample, Its expression is
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线的峰位,拍摄样品峰位放大像的方法:放入样品前,把环形光栅像固定在角度信号响应曲线的峰位,即令ψx=0,n=2,然后放置好样品,探测器可以拍摄到样品的峰位放大像,其表达式为In the step S8 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the method of fixing the ring grating image of the beam splitting grating at the peak position of the angle signal response curve, and taking the enlarged image of the peak position of the sample : Before placing the sample, fix the ring grating image at the peak position of the angle signal response curve, that is, set ψ x = 0, n = 2, and then place the sample, the detector can capture the peak magnified image of the sample, its expression for
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S8中,所述的把分束光栅的环形光栅像固定在角度信号响应曲线的下坡位,拍摄样品下坡位放大像的方法:放入样品前,把环形光栅像固定在角度信号响应曲线的下坡位,即令ψx=0,n=3,然后放置好样品,探测器可以拍摄到样品的下坡放大像,其表达式为In the step S8 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the annular grating image of the beam splitting grating is fixed at the downhill position of the angular signal response curve, and the downhill magnified image of the sample is taken The method: before placing the sample, fix the ring grating image at the downhill position of the angle signal response curve, that is, set ψ x = 0, n = 3, and then place the sample, the detector can take a downhill magnified image of the sample, Its expression is
在所述X射线微分相位衬度显微镜系统的二维成像方法的步骤S9中,所述的从样品的谷位放大像、上坡放大像、峰位放大像、下坡放大像中提取样品在物空间的吸收像、折射像和散射像方法如下所述:提取样品在物空间吸收像的公式为:In the step S9 of the two-dimensional imaging method of the X-ray differential phase contrast microscope system, the sample is extracted from the valley magnified image, the uphill magnified image, the peak magnified image, and the downhill magnified image of the sample. The absorption image, refraction image and scattering image method of the object space are as follows: the formula for extracting the absorption image of the sample in the object space is:
提取样品在物空间折射像的公式为:The formula for extracting the refraction image of the sample in object space is:
在散射可以忽略的条件下,提取样品在物空间折射像的公式简化为:Under the condition that the scattering can be ignored, the formula for extracting the refraction image of the sample in the object space is simplified as:
提取样品在物空间散射像的公式为:The formula for extracting the scattering image of the sample in the object space is:
在折射可以忽略的条件下,提取样品在物空间散射像的公式简化为:Under the condition that the refraction can be ignored, the formula for extracting the scattering image of the sample in the object space is simplified as:
另外,作为一种可选的实施方式,X射线光源10、聚光镜1、分束光栅3’可以集成为一个元件,即周期为微米量级的X射线环形栅源,如图7所示,此时在X射线环形栅源出口处沿光轴方向依次放置两个环形通光孔20和30,使来自X射线环形栅源的光束经过两个环形通光孔20和30后,形成照明样品的并具有空间周期结构的空心锥光束,并在物镜5后焦面附近形成环形栅源像及环形栅源像光束,环形分析光栅6’的形状和尺寸和环形栅源像相同。这种基于环形栅源的X射线微分相位衬度显微镜系统100’使得整个X射线微分相位衬度显微镜系统的长度可以进一步缩短,不仅可以降低X射线显微镜系统的制造成本,而且光的利用效率也能进一步提高,具有广阔的应用前景。In addition, as an optional implementation, the X-ray light source 10, the condenser lens 1, and the beam-splitting grating 3' can be integrated into one component, that is, an X-ray ring grating source with a period of micron order, as shown in FIG. 7 . At the same time, two annular light-through holes 20 and 30 are placed in sequence along the optical axis at the exit of the X-ray annular grid source, so that the light beam from the X-ray annular grid source passes through the two annular light-through holes 20 and 30 to form an illumination sample. It also has a hollow cone beam with a spatially periodic structure, and forms a ring source image and a ring source image beam near the rear focal plane of the objective lens 5. The shape and size of the ring analysis grating 6' are the same as that of the ring source image. The X-ray differential phase contrast microscope system 100' based on the ring grid source can further shorten the length of the entire X-ray differential phase contrast microscope system, which can not only reduce the manufacturing cost of the X-ray microscope system, but also improve the light utilization efficiency. It can be further improved and has broad application prospects.
由于上述基于环形栅源、环形栅源像和环形分析光栅的X射线微分相位衬度显微镜系统的二维成像方法,与前述基于分束光栅、分束光栅的环形光栅像和环形分析光栅的X射线微分相位衬度显微镜系统的二维成像方法类似,此处不再赘述。Due to the above-mentioned two-dimensional imaging method based on the ring grating source, the ring grating source image and the ring analysis grating X-ray differential phase contrast microscope system, compared with the aforementioned X-ray based on the beam splitting grating, the ring grating image of the beam splitting grating and the ring analysis grating The two-dimensional imaging method of the ray differential phase contrast microscope system is similar and will not be repeated here.
以上描述仅为本申请的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本申请中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离所述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。The above description is only a preferred embodiment of the present application and an illustration of the applied technical principle. Those skilled in the art should understand that the scope of the invention involved in this application is not limited to the technical solution formed by the specific combination of the above-mentioned technical features, but should also cover the technical solution formed by the above-mentioned technical features without departing from the inventive concept. Other technical solutions formed by any combination of or equivalent features thereof. For example, a technical solution formed by replacing the above-mentioned features with technical features with similar functions disclosed in (but not limited to) this application.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610617865.1A CN107664648B (en) | 2016-07-29 | 2016-07-29 | A kind of X-ray differential phase contrast microscopic system and its two-dimensional imaging method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610617865.1A CN107664648B (en) | 2016-07-29 | 2016-07-29 | A kind of X-ray differential phase contrast microscopic system and its two-dimensional imaging method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107664648A CN107664648A (en) | 2018-02-06 |
| CN107664648B true CN107664648B (en) | 2019-08-13 |
Family
ID=61122057
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610617865.1A Active CN107664648B (en) | 2016-07-29 | 2016-07-29 | A kind of X-ray differential phase contrast microscopic system and its two-dimensional imaging method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107664648B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108469443A (en) * | 2018-04-18 | 2018-08-31 | 北京航空航天大学 | X-ray grating differential phase contrast imaging method and device based on two dimension dislocation absorption grating |
| CN109801733B (en) * | 2018-12-29 | 2020-10-27 | 深圳大学 | Manufacturing method of X-ray absorption grating and X-ray absorption grating |
| CN111595877B (en) * | 2020-05-27 | 2022-03-29 | 合肥工业大学 | Multi-contrast image extraction method for X-ray diffraction enhanced imaging |
| CN111721784A (en) * | 2020-06-04 | 2020-09-29 | 中国科学院高能物理研究所 | A method and device for measuring imaging resolution of nanoscale zone plate |
| CN115308202B (en) * | 2021-05-07 | 2025-05-27 | 中国科学院高能物理研究所 | Phase contrast microscope, phase shift quantitative imaging method, system and storage medium |
| CN116297578A (en) * | 2021-12-20 | 2023-06-23 | 中国科学院深圳先进技术研究院 | X-ray Phase Quantitative Imaging Technology and Measurement Method |
| CN115656232B (en) * | 2022-12-05 | 2023-04-14 | 济南汉江光电科技有限公司 | X-ray differential phase contrast nano CT system and imaging method thereof |
| CN118328914B (en) * | 2023-07-24 | 2025-05-30 | 张江国家实验室 | Small-angle X-ray scattering measurement device and method based on multi-beam incidence |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103364416A (en) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | X ray differential phase contrast microscopic imaging system and imaging method |
| CN104622492A (en) * | 2013-11-11 | 2015-05-20 | 中国科学技术大学 | X-ray grating phase-contrast imaging device and method |
| CN105675631A (en) * | 2016-01-05 | 2016-06-15 | 合肥泰禾光电科技股份有限公司 | Fast fan-beam geometric phase contrast CT imaging device and method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006015358B4 (en) * | 2006-02-01 | 2019-08-22 | Paul Scherer Institut | Focus / detector system of an X-ray apparatus for producing phase-contrast images, associated X-ray system and storage medium and method for producing tomographic images |
| DE102008048683A1 (en) * | 2008-09-24 | 2010-04-08 | Siemens Aktiengesellschaft | Method for determining phase and / or amplitude between interfering adjacent X-rays in a detector pixel in a Talbot interferometer |
-
2016
- 2016-07-29 CN CN201610617865.1A patent/CN107664648B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103364416A (en) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | X ray differential phase contrast microscopic imaging system and imaging method |
| CN104622492A (en) * | 2013-11-11 | 2015-05-20 | 中国科学技术大学 | X-ray grating phase-contrast imaging device and method |
| CN105675631A (en) * | 2016-01-05 | 2016-06-15 | 合肥泰禾光电科技股份有限公司 | Fast fan-beam geometric phase contrast CT imaging device and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107664648A (en) | 2018-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107664648B (en) | A kind of X-ray differential phase contrast microscopic system and its two-dimensional imaging method | |
| CN103377746B (en) | Method for realizing super-resolution imaging of microscope system | |
| US10352880B2 (en) | Method and apparatus for x-ray microscopy | |
| CN103364416B (en) | X ray differential phase micro imaging system and formation method | |
| Giavazzi et al. | Scattering information obtained by optical microscopy: differential dynamic microscopy and beyond | |
| US20150055745A1 (en) | Phase Contrast Imaging Using Patterned Illumination/Detector and Phase Mask | |
| CN102830102B (en) | Method and device for hollow focused light spot excitation-based confocal microscopy | |
| WO2017213996A1 (en) | Method and apparatus for x-ray microscopy | |
| CN108169257B (en) | There are three types of the nanometer CT three-D imaging methods of quantitative imaging mechanism for tool | |
| CN101632134A (en) | Holographic microscopy system and method for optical trapping and material inspection | |
| CN101303302A (en) | Method and device for realizing quasi-confocal fluorescence microscopy with dynamic speckle illumination | |
| CN114113019A (en) | A kind of array scanning super-resolution microscopy imaging device, method, equipment and storage medium based on multiple signal classification algorithm | |
| JP5981443B2 (en) | Observation device | |
| US9297990B2 (en) | Confocal microscope | |
| EP3722861A1 (en) | Observation method using microscope and transmission-type microscopic device | |
| CN101893755A (en) | Fluorescence microscopy method and apparatus using a quadrangular pyramid mirror to produce structured illumination | |
| US20200014837A1 (en) | Method for observing a sample by lensless imaging, with a spatial dispersion in the sample taken into account | |
| US20250314870A1 (en) | Microscope device and data generation method using microscope | |
| JP2009169420A (en) | Polarization phase microscope | |
| JP2020137429A (en) | Cell observation system and cell observation method | |
| CN105319195A (en) | Super-resolution structure detection array confocal fluorescence imaging device and imaging method thereof | |
| US20210348998A1 (en) | Method and apparatus for detecting nanoparticles and biological molecules | |
| Sung | Snapshot three-dimensional absorption imaging of microscopic specimens | |
| Geilhufe et al. | Extracting depth information of 3-dimensional structures from a single-view X-ray Fourier-transform hologram | |
| CN105319196A (en) | Super-resolution structure detection confocal fluorescence imaging device and imaging method thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230805 Address after: 25F03B, Phase 2, Zhongyin Square, Luoyuan Street, Lixia District, Jinan City, Shandong Province, 250012 Patentee after: Jinan Hanjiang Photoelectric Technology Co.,Ltd. Address before: 100049 Shijingshan District, Yuquanlu Road, Beijing No. 19 (b) Patentee before: INSTITUTE OF HIGH ENERGY PHYSICS, CHINESE ACADEMY OF SCIENCES |