Clarification on Generalized Lau condition for X-ray interferometers based on dual phase gratings

TL;DR

This paper corrects misconceptions about the Lau-condition in dual phase grating x-ray interferometry, deriving a new generalized condition through geometrical analysis to improve design and implementation.

Contribution

It provides a new generalized Lau-condition specific to dual phase grating interferometry, correcting previous assumptions and aiding practical design.

Findings

The traditional Lau-condition does not directly apply to dual phase grating interferometry.

A new generalized Lau-condition was derived based on geometrical analysis.

The new condition improves the design accuracy for x-ray interferometers.

Abstract

To implement dual phase grating x-ray interferometry with x-ray tubes, one needs to incorporate an absorbing source grating. In order to attain good fringe visibility, the period of a source grating should be subject to a stringent condition. In literature some authors claim that the Lau-condition in Talbot-Lau interferometry can be literally transferred to dual phase grating interferometry. In this work we show that this statement in literature is incorrect. Instead, through an intuitive geometrical analysis of fringe formation, we derived a new generalized Lau-condition that provides a useful design tool for implementation of dual phase grating interferometry.

Figures

Equations

$$p_{{}_0}=\frac{R_1}{R_2}{p}_{{}_2}.$$

$$p_{{}_0}^{\ast }=\frac{R_s}{R_g}{p}_{\mathrm{proj}}=\frac{R_s+R_g}{R_g}{p}_{{}_1}.$$

$$M_{g_1}=\frac{R_s+R_g+R_d}{R_s};M_{g_2}=\frac{R_s+R_g+R_d}{R_s+R_g},$$

$$\exp \left[i2\pi l\frac{x}{M_{g_1}{p}_{{}_1}}\right]\cdot \exp \left[i2\pi r\frac{x}{M_{g_2}{p}_{{}_2}}\right]=\exp \left[i2\pi x\cdot \left(\frac{l}{M_{g_1}{p}_{{}_1}}+\frac{r}{M_{g_2}{p}_{{}_2}}\right)\right],$$

$$p_{{}_{\mathrm{fr}}}={⟮\frac{-1}{M_{g_1}{p}_{{}_1}}+\frac{1}{M_{g_2}{p}_{{}_2}}⟯}^{-1}=\frac{R_s+R_g+R_d}{R_g/p_{{}_2}+R_s(1/p_{{}_2}-1/p_{{}_1})}.$$

$$⟮\frac{R_g+R_d}{M_{g_1}{R}_s{p}_{{}_1}}-\frac{R_d}{M_{g_2}(R_s+R_g)p_{{}_2}}⟯\cdot p_{{}_0}=1.$$

$$p_{{}_0}=\frac{R_s+R_g+R_d}{R_g/p_{{}_1}+R_d(1/p_{{}_1}-1/p_{{}_2})}$$

Full text

OCIS codes: (110.6760) Talbot and self-imaging effects; (110.7440) X-ray imaging; (340.7440) X-ray imaging; (340.7450) X-ray interferometry.

References and links

• [1]

A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, H. Takai, and Y. Suzuki, “Demonstration of x-ray talbot interferometry,” Jpn. J. Appl. Phys. 42, L866–L868 (2003).

• [2]

T. Weitkamp, A. Diaz, C. David, F. Pfeiffer, M. Stampanoni, P. Cloetens, and E. Ziegler, “X–ray phase imaging with a grating interferometer,” Opt. Express 13, 6296–6304 (2005).

• [3]

A. Momose, “Recent advances in x-ray phase imaging,” Jpn. J. Appl. Phys. 44, 6355–6367 (2005).

• [4]

F. Pfeiffer, T. Weitkamp, O. Bunk, and C. David, “Phase retrieval and differential phase-contrast imaging with low-brilliance x-ray sources,” Nat. Phys. 2, 258–261 (2006).

• [5]

W. Yashiro, Y. Terui, K. Kawabata, and A. Momose, “On the origin of visibility contrast in x-ray talbot interferometry,” Opt. Express 18, 16890–16901 (2010).

• [6]

P. Zhu, K. Zhang, Z. Wang, Y. Liu, X. Liu, Z. Wu, S. McDonald, F. Marone, and M. Stampanoni, “Low-dose, simple, and fast grating-based x-ray phase-contrast imaging,” Proc Natl. Acad. Sci. USA 107, 13576–13581 (2010).

• [7]

N. Bevins, J. Zambelli, K. Li, Z. Qi, and G.-H. Chen, “Multicontrast x-ray computed tomography imaging using talbot-lau interferometry without phase stepping,” Med. Phys. 39, 424–428 (2012).

• [8]

X. Tang, Y. Yang, and S. Tang, “Characterization of imaging performance in differential phase contrast ct compared with the conventional ct: Spectrum of noise equivalent quanta neq(k),” Med. Phys. 39, 4367–4382 (2012).

• [9]

T. Teshima, Y. Setomoto, and T. Den, “Two-dimensional gratings-based phase-contrast imaging using a conventional x-ray tube,” Opt. Lett. 36, 3551–3553 (2011).

• [10]

E. Bennett, R. Kopace, A. Stein, and H. Wen, “A grating-based single shot x-ray phase contrast and diffraction method for in vivo imaging,” Med. Phys. 37, 6047–6054 (2010).

• [11]

M. Jiang, C. Wyatt, and G. Wang, “X-ray phase-contrast imaging with three 2d gratings,” Int. J. Biomed. Imaging pp. 827152–927159 (2008).

• [12]

H. Itoh, K. Nagai, G. Sato, K. Yamaguchi, T. Nakamura, T. Kondoh, C. Ouchi, T. Teshima, Y. Setomoto, and T. Den, “Two-dimensional grating-based x-ray phase-contrast imaging using fourier transform phase retrieval,” Opt. Express 19, 3339–3346 (2011).

• [13]

A. Bravin, P. Coan, and P. Suortti, “X-ray phase-contrast imaging: from pre-clinical applications towards clinics,” Physics in Medicine and Biology 58, R1–R35 (2013).

• [14]

A. Momose, H. Kuwabara, and W. Yashiro, “X-ray phase imaging using lau effects,” Appl. Phys. Express 4, 066603 (2011).

• [15]

N. Morimoto, S. Fujino, K. Ohshima, J. Harada, T. Hosoi, H. Watanabe, and T. Shimura, “X-ray phase contrast imaging by compact talbot-lau interferometer with a single transmission grating,” Opt. Lett. 39, 4297–4300 (2014).

• [16]

N. Morimoto, S. Fujino, A. Yamazaki, Y. Ito, T. Hosoi, H. Watanabe, and T. Shimura, “Two dimensional x-ray phase imaging using single grating interferometer with embedded x-ray targets,” Opt. Express 23, 16582–16588 (2015).

• [17]

A. Yan, X. Wu, and H. Liu, “A general theory of interference fringes in x-ray phase grating imaging,” Med. Phys. 42, 3036–3047 (2015).

• [18]

A. Yan, X. Wu, and H. Liu, “Predicting visibility of interference fringes in x-ray grating interferometry,” Opt. Express 24, 15927–15939 (2016).

• [19]

A. Yan, X. Wu, and H. Liu, “Polychromatic x-ray effects on fringe phase shifts in grating interferometry,” Opt. Express 25, 6053–6068 (2017).

• [20]

A. Yan, X. Wu, and H. Liu, “Beam hardening correction in polychromatic x-ray grating interferometry,” Opt. Express 25, 24690–24704 (2017).

• [21]

T. Donath, M. Chabior, F. Pfeiffer, O. Bunk, E. Reznikova, J. Mohr, E. Hempel, S. Popescu, M. Hoheisel, M. Schuster, J. Baumann, and C. David, “Inverse geometry for grating-based x-ray phase-contrast imaging,” J. Appl. Phys. 106, 054703 (2009).

• [22]

H. Miao, A. Panna, A. Gomella, E. Bennett, S. Znati, L. Chen, and H. Wen, “A universal moiré effect and application in x-ray phase-contrast imaging,” Nature Physics 12, 830–836 (2016).

• [23]

M. Kagias, Z. Wang, K. Jefimovs, and M. Stampanoni, “Dual phase grating interferometer for tunable dark-field sensitivity,” Appl. Phys. Lett. 110, 014105 (2017).

• [24]

A. Yan, X. Wu, and H. Liu, “Quantitative theory of x-ray interferometers based on dual phase grating: fringe period and visibility,” Opt. Express 26, 23142–23155 (2018).

• [25]

J. Bopp, V. Ludwig, M. Seifert, G. Pelzer, A. Maier, G. Anton, and C. Riess, “Simulation study on x-ray phase contrast imaging with dual-phase gratings,” International Journal of Computer Assisted Radiology and Surgery 14, 3–10 (2019).

1 Introduction

Currently, the Talbot-Lau x-ray interferometry is widely used for x-ray differential phase contrast imaging[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. In Talbot-Lau grating-based interferometry, a phase grating is employed as a beam splitter to split x-ray into diffraction orders. The interference between the diffracted orders forms intensity fringes. The sample imprints x-ray beam with phase shifts and distorts the intensity fringes. Analyzing the fringes one reconstructs sample attenuation, phase gradient and dark field images[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. To increase the grating interferometer’s sensitivity, one needs to use fine-pitch phase gratings of periods as small as few micrometers. But in medical imaging and material science applications, it is only feasible to utilize common imaging detectors, whose pixels are of a few tens of micrometers. To enable fringe detection with common image detectors, one way is to use a fine absorbing grating placed as shown in Fig. 1(a) to indirectly detect fringe patterns through grating scanning, which is also called phase stepping procedure[1, 2, 3, 4]. However, the absorbing grating blocks more than half of transmitting x-ray, and will increase radiation dose in imaging exams.

Recently, dual phase grating x-ray interferometry demonstrated its attractive advantages[22, 23]. A typical dual phase grating interferometer employs two phase gratings $G_{1}$ and $G_{2}$ as the beam splitters, as is shown in Fig. 1(b). The split waves transmitting through the phase gratings interfere with each other, creating different diffraction orders, including a beat pattern in the intensity fringe pattern[24]. The imaging detector $D$, as shown in Fig. 1(b), has a pixel size much larger than periods of both the phase gratings. Due to pixel averaging effects, the detector just resolves the beat patterns of large periodicities. In this way, a common imaging detector with resolution of tens micrometers can directly resolve the fringes generated by fine phase gratings. Hence, different from Talbot-Lau setups[1, 2, 3, 4, 5, 6, 7, 8, 9, 10], dual phase grating enables direct fringe detection without need of absorbing analyzer grating. This advantage brings significant radiation dose reduction as compared to Talbot-Lau interferometry. Moreover, different from the inverse geometry setups of Talbot-Lau interferometry[21], a dual phase grating interferometer keeps the system’s length compact, since the fringe period can be conveniently tuned by adjusting the grating-spacing, without increase of system size[22, 23, 24].

In order to attain good fringe visibility in grating-based x-ray interferometry, a necessary condition is to provide spatially coherent illumination of phase gratings, such as that achieved by use of synchrotron radiation or micro-focus x-ray tubes. But for potential medical imaging applications of dual grating interferometry, it is challenging to find robust x-ray sources to provide sufficient spatial coherence and adequate x-ray flux. As a solution, one may employ an x-ray tube equipped with an absorbing source grating $G_{0}$, which serves as an aperture mask to divide the focus spot into an array of mutually incoherent virtual line sources[4, 14]. Similarly, one may utilize a periodic array $G_{0}$ of micro-anode-sources embedded in the anode for the same purpose[15, 16]. For Talbot-Lau interferometry with a source grating, the well-known Lau condition on the period of a source grating should be satisfied for attaining good fringe visibility[4, 14]. The Lau condition for incorporating a source grating has an intuitive geometric interpretation. When Lau condition is satisfied, each displayed line-source of a source grating should shift the fringe pattern by an integral multiples of the fringe period. With this condition satisfied, the fringe pattern generated by all the line-sources in the source grating are constructively superimposed[4, 14, 15, 16]. Figure 2(a) shows a schematic of geometric configuration of a typical Talbot-Lau interferometer. It is easy to see from Fig. 2(a) that the fringe generated by this line source is displaced by $\Delta x=(R_{2}/R_{1})p_{{}_{0}}$, where $p_{{}_{0}}$ is the pitch of the source grating, $R_{1}$ is source-grating $G_{0}$-phase grating $G_{1}$ distance, and $R_{2}$ is the distance between $G_{1}$-grating and the analyzer grating $G_{2}$, which is placed in front of the detector. For constructive superposition one requires $\Delta x$ should be equal to fringe period $p_{{}_{\mathrm{fr}}}$. In Talbot-Lau interferometry the fringe period $p_{{}_{\mathrm{fr}}}$ is generally equal to the projection $p_{\mathrm{proj}}$ of the period of phase grating $G_{1}$, that is, $p_{{}_{\mathrm{fr}}}=p_{\mathrm{proj}}=(R_{1}+R_{2})/R_{1}p_{{}_{1}}$. Note that for $\pi$-grating, the fringe period is reduced to $p_{{}_{\mathrm{fr}}}=p_{\mathrm{proj}}/2$. But with polychromatic x-ray, phase shift of a grating varies with photon energy. Hence, the fringe period of a Talbot-Lau interferometry with polychromatic x-ray is equal to the projection $p_{\mathrm{proj}}$. Note that in Talbot-lau interferometry the period $p_{{}_{2}}$ of the absorbing grating is always set to $p_{{}_{2}}=p_{{}_{\mathrm{fr}}}$. Since $\Delta x=p_{{}_{\mathrm{fr}}}$ is required for coherent superposition of fringes, so the source grating pitch should be set to

[TABLE]

This is the so-called Lau condition for constructive fringe superposition in Talbot-Lau interferometry.

But for dual phase grating interferometry, the fringes are generated by two phase gratings, so the the original Lau condition of Eq. (1) is not applicable anymore. In order to implement the dual phase grating interferometry for wide applications in medical imaging, there is a pressing need to generalize the original Lau condition to dual phase grating interferometry. In literature a recently published paper presented an analysis on the validity of the Lau condition for dual phase grating interferometry[25]. In that paper the authors postulated that the intensity fringe pattern results from low frequency Moire pattern generated by $G_{1}$ and $G_{2}$ at the $G_{2}$-plane[25]. Under this hypothesis, the paper claims that the original Lau condition on the source grating is literally transferred to dual phase grating interferometry. Rewriting Eq. (1) with the notations employed in Fig. 1(b) for a dual phase grating interferometer, the authors of that paper claimed that the source grating pitch should satisfy following condition for achieving coherent fringe superposition in dual phase grating interferometry:

[TABLE]

Here $R_{s}$ is the distance from source grating to the first grating $G_{1}$, and $R_{g}$ is the spacing between the two phase gratings, as is shown in Fig. 1(b). For sake of convenience in discussion we denote the source grating period determined by above equation (Eq. (2)) as $p_{{}_{0}}^{\ast}$. Note that this $p_{{}_{0}}^{\ast}$ given by Eq. (2) is independent of the distance $R_{d}$ between $G_{2}$-grating and detector. Unfortunately, the Lau condition of Eq. (2) for dual phase grating interferometry is incorrect. It fails to take into account the wave propagation from $G_{2}$ grating to the detector. Setting out to clarify on this issue, in this work we present an intuitive derivation of the condition of constructive fringe superposition in dual phase grating interferometry with a source grating. In section 2, we present an intuitive derivation of the generalized Lau condition for dual phase grating interferometry incorporating a source grating. In section 3, we present simulation results that validate our generalized Lau condition for dual phase grating interferometry. We conclude the work in section 4. We hope that our clarification on the generalized Lau condition may help researchers in design and implementation of dual phase grating x-ray interferometry incorporating a source grating.

2 Methods

We consider a typical setup of dual phase x-ray interferometry, as is shown in Fig. 1(b). Different from conventional Talbot-Lau X-ray interferometry, which uses only a single phase-grating, the new technique employs two phase gratings $G_{1}$ and $G_{2}$ as the beam splitters. For high-sensitivity the phase gratings $G_{1}$ and $G_{2}$ are made of small periods of close to one micrometer, and we assume that $G_{1}$’s period $p_{{}_{1}}$ is equal or close to $G_{2}$’s period $p_{{}_{2}}$. In addition, the detector has pixel size $p_{\!{}_{{}_{\mathrm{D}}}}$ that is about few tens of micrometers, so detector pixel size $p_{\!{}_{{}_{\mathrm{D}}}}\!\!\gg\!\!p_{{}_{1}},p_{{}_{2}}$. In order to employ sources of large focal spots, a strongly absorbing source grating $G_{0}$ of narrow slits is used to break source into an array of mutually incoherent virtual line sources. But in order to have the fringes generated by these line sources constructively superimposed, the period $p_{{}_{0}}$ of the source grating should satisfy a certain condition, which we called the generalized Lau condition in dual phase-grating interferometry. Our task in this work is to derive the generalized Lau condition by using a simple geometric analysis.

Let us consider a setup with a single narrow-slit source first. For sake of convenience in discussion, we define several geometric magnification factors $M_{g_{2}}$ and $M_{g_{1}}$ as:

[TABLE]

where $M_{g_{1}}$ is the geometric magnification factor from the $G_{1}$ plane to detector plane, and $M_{g_{2}}$ is the geometric magnification factor from the $G_{2}$ plane to detector plane. Our analysis starts from the x-ray irradiance $I_{R_{s}+R_{g}+R_{d}}(x,y)$ at detector entrance. Obviously, in absence of $G_{2}$ grating, an intensity diffraction order generated by $G_{1}$ phase grating alone would be represented by $\exp\left[i2\pi lx/(M_{g_{1}}p_{{}_{1}})\right]$, where $l$ is an integer and $p_{{}_{1}}$ is the period of $G_{1}$ grating. Similarly, an intensity diffraction order generated by the second phase grating $G_{2}$ alone would be denoted by $\exp\left[i2\pi rx/(M_{g_{2}}p_{{}_{2}})\right]$, where $r$ is an integer and $p_{{}_{2}}$ is the period of the second phase grating. As is shown from a quantitative theory of dual phase grating interferometry[24], x-ray irradiance at detector entrance is a result of cross-modulation between the intensity patterns generated by the two phase-gratings $G_{1}$ and $G_{2}$ respectively. Specifically, the intensity pattern is a weighted sum of different diffracted orders. Each of the diffraction orders in dual phase grating interferometry is represented by a product of

[TABLE]

and each of the diffraction orders is indexed by two integers $(l,r)$. It is easy to see from above equation that the intensity fringe of order $(l,r)$ has a periodicity of $\left(l/(M_{g_{1}}p_{{}_{1}})+r/(M_{g_{2}}p_{{}_{2}})\right)^{-1}$. Since grating periods are few micrometers only, so the interference fringes are generally so fine that most of them are too fine to be detected by common imaging detector. However, among the fringe patterns, there are beat patterns formed by those diffraction orders characterized by $l=-r$. These beat patterns are indeed generated by the interference beat patterns. The beat pattern consists of fundamental frequency and its harmonics resolved by the imaging detector. The period of the first harmonics of beat patterns is[24]

[TABLE]

On the other hand, period of the $l$-th harmonics is $p_{{}_{\mathrm{fr}}}/l$. Hence the period of the intensity fringe is $p_{{}_{\mathrm{fr}}}$ of Eq. (5). According to Eq. (5), the beat pattern periods can be much larger than grating periods $p_{{}_{1}}$ and $p_{{}_{2}}$. For example, if the system geometry is set in such a way that $R_{s}+R_{d}\!\gg\!R_{g}$ and $p_{{}_{1}}\approx p_{{}_{2}}$ then the resulting fringe period $p_{{}_{\mathrm{fr}}}\!\gg\!p_{{}_{1}},p_{{}_{2}}$. Consequently, the beat pattern may be resolved by a common imaging detector. On the other hand, as long as detector period $p_{\!{}_{{}_{\mathrm{D}}}}\!\gg\!p_{{}_{1}},p_{{}_{2}}$, the detector renders all other fine fringes with $l\neq-r$ orders to a constant background, owing to detector-pixel averaging effect.

Once the fringe formation mechanism is understood, we are now ready to derive the generalized Lau condition for a dual phase grating setup incorporating a source grating. Consider an off-center line source that is displaced by $p_{{}_{0}}$ from the center. As is demonstrated in Fig. 2(b), from the similar triangles relationship, this source-displacement $p_{{}_{0}}$ will cause a fringe shift $\Delta x_{1}=p_{{}_{0}}(R_{g}+R_{d})/R_{s}$ for the $G_{1}$-associated fringes, so its diffraction order will become $\exp\left[i2\pi l(x+\Delta x_{1})/(M_{g_{1}}p_{{}_{1}})\right]$. By the same reasoning, this off-center line source causes a fringe shift $\Delta x_{2}=p_{{}_{0}}R_{d}/(R_{s}+R_{g})$ for the $G_{2}$-associated fringes, so its diffraction order will become $\exp\left[i2\pi r(x+\Delta x_{2})/(M_{g_{2}}p_{{}_{2}})\right]$. Since the fringes resolved by the detector are the beat patterns formed by cross-modulation between the $G_{1}$ and $G_{2}$-associated fringes with diffracted orders $l=-r$, so this off-center line source generates a total fringe phase shift (in radians) $\Delta\Phi=2\pi l\left(\Delta x_{1}/(M_{g_{1}}p_{{}_{1}})-\Delta x_{2}/(M_{g_{2}}p_{{}_{2}})\right)$ in the resolved intensity fringe. To achieve constructive fringe superposition for all diffraction orders, one should make $\Delta x_{1}/(M_{g_{1}}p_{{}_{1}})-\Delta x_{2}/(M_{g_{2}}p_{{}_{2}})=1$, or equivalently, one should require:

[TABLE]

Using Eq. (3), we can rewrite the condition of fringe constructive superposition as:

[TABLE]

Equation (7) is the generalized Lau condition, for dual phase grating interferometry incorporating a source grating. This condition determines the source grating periodicity for achieving constructive superposition of fringes generated by all the line sources in a source grating.

3 Results

In order to validate the generalized Lau condition derived above we conducted numerical simulations, in which four interferometer setups were simulated. The four setups employ the same phase gratings and identical geometric configurations, such as $R_{s}$, $R_{g}$, and $R_{d}$, but differ in source configurations. In the first setup, the interferometer consists of a 20-keV point source, two $\pi$-phase gratings of periods $p_{{}_{1}}=1\mu\mathrm{m}$, and $p_{{}_{2}}=1.1\mu\mathrm{m}$ respectively, and a detector of $19.3\mu$m pixels. The geometry configuration of this interferometer was set to $R_{s}=45$cm, $R_{g}=4$cm, and $R_{d}=38$cm. This geometric configuration was selected according to the fringe visibility formulas for good fringe visibility[24]. This setup serves as a reference for fringe visibility comparison, as the point source is spatially coherent. For the other three interferometer setups, we replaced the point source with a focal spot of 0.35 mm in width, while keeping otherwise the same geometric configuration as the first setup. In the second setup no source grating was used. In the third and fourth setups different source gratings were placed in front of the focal spot respectively.

In the third setup a source grating of period $p_{0}=11.67\mu$m was incorporated. Here the period $p_{{}_{0}}$ was set according to Eq. (7), the generalized Lau condition derived in this work. As a comparison, in the $4^{\mathrm{th}}$ setup another source grating $G_{0}^{\ast}$ of period $p_{{}_{0}}^{\ast}=12.25\mu$m was used. This period $p_{{}_{0}}^{\ast}$ was set according to Eq. (2), the Lau condition used in[25]. The aperture-sizes of both source gratings were set to the same as $a=2\mu$m. Intensity fringe formation was simulated numerically as results of wave Fresnel diffraction propagating from the source grating, through the two phase gratings, and finally to the imaging detector.

The curves in Fig. 3(a), are the plots of intensity patterns associated with different source settings. In Fig. 3(a) four curves are presented, they all have the same period of $p_{{}_{\mathrm{fr}}}=191.4\mu$m but different fringe visibility. The blue curve is the plot for the setup with the point source, while the black curve is the plot for the setup with a 0.35 mm-wide source but without source grating. Compared to the blue curve, the black curve shows diminishing fringe visibility. This fringe visibility loss is due to lack of sufficient spatial coherence because of the source width. Incorporating a source grating can help but only if its period is set according to the generalized Lau condition of Eq. (7). The red and green curves in Fig. 3(a) correspond to the setups with the same $0.35$mm wide source but with different source grating periods. The red curve is for the source grating setup based on Eq. (7), that is, the generalized Lau condition. The red curve demonstrates good fringe visibility, which is close to that with the point source. The difference lies in finite aperture size ($a=2\mu$m) of the source grating. The smaller the aperture size is, the better the fringe visibility. Hence, the red curve clearly validates the generalized Lau condition for determination of source grating period. In stark contrast, the green curve, which is associated with the source grating setup based on Eq. (2), depicts as poor fringe visibility as the setup without source grating does. The black and red/green curves demonstrate that incorporating a source grating can restore fringe visibility for setups with finite focal spots, provided that source grating period is set according to the generalized Lau condition of Eq. (7). The green curve shows clearly that a source grating set with the Lau condition of Eq. (2) is unable to restore fringe visibility. Hence the Lau condition of Eq. (2) is not applicable for dual phase grating interferometry. Several other interferometer setups have also been tested, all confirmed the validity of the generalized Lau condition of Eq. (7) for incorporating source gratings.

In addition, in Fig. 3(b), (c) and (d) we present the intensity maps obtained when detector pixel size is set to $19.33\mu$m. These intensity maps are corresponding to setups associated with the blue, red and green curves in Fig. 3(a), respectively.

4 Discussion and conclusions

Dual phase grating x-ray interferometry is a promising new technique of grating based x-ray differential phase contrast imaging[22, 23, 24]. To implement this new X-ray interferometry technique for a setup with an X-ray tube of finite focal spot, an absorbing source grating is required. A source grating serves as an aperture mask to divide the focus spot into an array of mutually incoherent virtual line sources[4, 14]. In order to attain good fringe visibility, a necessary condition is that the fringe pattern generated by all the line-sources in the source grating are constructively superimposed. This requirement imposes a stringent design criterion on the period of a source grating in dual phase grating x-ray interferometry. Unfortunately, currently there is confusion in literature about this important design criterion. As is mentioned earlier, some authors argue that the well-known Lau-condition in Talbot-Lau interferometry can be literally transferred to dual phase grating interferometry[25]. In their reasoning they fail to consider effects on intensity fringes of x-ray wave propagation from the second grating to the imaging detector. Based on this flawed reasoning, they claim that the period of a source grating should be set according Eq. (2), the Lau condition for dual phase grating interferometry. Note that the Lau-condition Eq. (2) is independent of the distance between the second grating and detector. However, as is shown in section 3, a source grating setup based on Eq. (2) fails to provide good fringes visibility. In this work, we set out to clarify the design criterion on source grating periods in dual phase grating x-ray interferometry. We call this design criterion as the generalized Lau-condition. In the derivation we noted that it is important to consider the effects on intensity fringe of full wave propagation from the source to the detector. Hence it is incorrect to ignore the effects of wave propagation from second phase grating to the detector. Although one can derive the generalized Lau-condition from tedious derivation of coherence degree of x-ray illumination[24, 22], but in this work we rather give a simpler and more intuitive derivation of the generalized Lau-condition of Eq. (7), as is presented in section 2. In fact, Eq. (7) is a necessary condition under which fringe patterns generated by all the line-sources in the source grating are constructively superimposed. The simulation results presented in section 3 validated generalized Lau-condition of Eq. (7) as a design criterion on source grating periods in dual phase grating interferometry.

One remark on Eq. (7) is on order. In practice one usually has $R_{s}+R_{d}\!\gg\!R_{g}$, hence the pitch $p_{0}$ of the source grating is much larger than the phase grating pitches $p_{1}$ and $p_{2}$. This is a significant advantage of the dual phase grating setup as compared to the inverse geometry setup in Talbot interferometry[21]. Although both allow to use large detector pixels for intensity fringe detection, but in the inverse geometry the pitch of source grating required to be as small as of few micrometers[21]. Such narrow-pitch absorbing source grating of high-aspect ratios is hard to fabricate. With dual phase grating setups, Eq. (7) shows that the pitches of the source grating can be as large as few tens micrometers, thereby such source gratings are easier to fabricate.

In conclusion, in this work we clarify the issue of what is the condition on source grating period for good fringe visibility. In literature some authors claim that the Lau-condition in Talbot-Lau interferometry can be literally transferred to dual phase grating interferometry. In this work we show that this statement in literature is incorrect and the Lau-condition is not applicable to dual phase grating interferometry. Instead, through an intuitive geometrical analysis of fringe formation, we derived a new generalized Lau-condition for dual phase grating interferometry. The generalized Lau-condition provides a useful design tool for implementation of dual phase grating interferometry.

Funding

National Institutes of Health (NIH) (1R01CA193378).