Abstract
X-ray grating interferometry (XGI) can provide access to unresolved sub-pixel information by utilizing the so-called dark-field or visibility reduction contrast. A recently developed variant of conventional XGI named dual-phase grating interferometer, based only on phase-shifting structures, has allowed for straightforward micro-structural investigations over multiple length scales with conventional X-ray sources. Nonetheless, the theoretical framework of the image formation for the dark-field signal has not been fully developed yet, thus hindering the quantification of unresolved micro-structures. In this work, we expand the current theoretical formulation of dual-phase grating interferometers taking into account polychromatic sources and beam hardening effects. We propose a model that considers the contribution of beam hardening to the visibility reduction and accounts for it. Finally, the method is applied to previously acquired and new experimental data showing that discrimination between actual micro-structures and beam hardening effects can be achieved.
© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
X-ray grating interferometry (XGI) [1,2] is an imaging technique that simultaneously provides three complementary contrasts; absorption, differential phase, and visibility reduction also known as dark-field contrast. There has been a great interest, particularly in the dark-field contrast [3], as it provides information about underlying scattering structures, which are much smaller than the pixel size of the utilized detector [4–6]. The utility of the dark-field contrast has been demonstrated by applications both in the industrial and medical fields. For example, it can be used to obtain the orientation of unresolved carbon fibers in composite materials [7–9] which determine their mechanical properties. Furthermore, the observations of the time-evolving hydration process of mortars were enabled by the high sensitivity of the dark-field contrast to unresolved pores [10,11]. Recently, dark-field was also used to detect unresolved structural defects in composite materials [12]. In the medical research community, it was demonstrated that the dark-field signal could be used to differentiate between various types of microcalcifications which are associated with different levels of breast cancer malignancy [13,14]. Similarly, distinct types of kidney stones can also be discriminated [15]. Another highly promising medical application is lung imaging, by exploiting the strong scattering taking place in the microstructure of the lung, high contrast images can be obtained, opening the potential for improved diagnosis [16,17]. All of these applications are aiming at utilizing the high contrast of the dark-field signal compared to conventional transmission, and do not exploit the fact that quantitative information about the scattering structure can be obtained. The reason for this is that a GI setup has a specific length scale sensitivity which is defined by the autocorrelation length of the interferometer [6]. This effectively means that a single dark-field measurement corresponds to a single point of the autocorrelation function of the scattering structure which is insufficient for microstructural characterization. Therefore, to obtain quantitative information about the unresolved microstructures in a sample, multiple measurements over a range of autocorrelation lengths needs to be performed [18–20].
In general, the autocorrelation length $\xi$ depends on the period of the interference pattern at the detector plane, the distance between the sample and the detector, and lastly, the wavelength. In a conventional XGI (Talbot-Lau interferometer, with three gratings) changing autocorrelation length can be achieved by either imaging the sample with phase gratings of different periods, an impractical and cost-inefficient solution, or by moving the sample along the optical axis, which in cone-beam geometries has the undesirable effect of varying magnification hindering further analysis. Changing the wavelength also has its disadvantages since, ideally, for each energy, a new set of gratings would be required. A so-called dual-phase interferometer was recently proposed to overcome these limitations of conventional XGI [21,22]. By using two phase gratings of small periods, enlarged fringes are generated at the detector plane, which can be recorded directly without the need for an analyzer grating [22]. Moreover, the relative distance between the two gratings can be tuned to change the fringe period, correspondingly changing the autocorrelation length. When XGI is performed with a polychromatic X-ray source, beam hardening effects from the sample can induce an apparent visibility reduction which does not correspond to small angle scattering [23–26], hence hindering any structural quantification. Nonetheless, in typical XGI setups, namely Talbot and Talbot-Lau interferometers, changing the autocorrelation length by moving the sample in the beam path neither changes the amount of beam hardening that takes place nor the spectral sensitivity (defined by the gratings and geometry) of the interferometer, hence leading to a constant bias across all measurements. On the contrary, in a dual-phase XGI, changing the relative distance between the gratings changes the geometry of the setup and thus the spectral acceptance of the system as well. Therefore, a varying effect of the sample beam hardening across different autocorrelation lengths is incorporated in the observed dark-field signal, which means that a usual imaging protocol where autocorrelation length is changed without accounting for beam hardening would not provide quantifiable structural information.
Although some recent works have presented the theoretical calculations for determining the fringe period seen at the detector and visibility [27] in the dual-phase XGI grating system, a complete theory incorporating the scattering and beam hardening effect of the sample is not available. This letter addresses the effect of the sample on the visibility reduction in the case of a polychromatic dual-phase XGI. We introduce a generalized formulation to calculate the interference fringe visibility at the detector when a sample is placed in the beam path using wave optics formulations. Furthermore, we demonstrate the varying effect of beam hardening on the visibility reduction across different autocorrelation lengths and eventually propose a method to account for it. Finally, we demonstrate the utility of this correction to quantify scattering structures and discriminate between visibility induced by beam hardening and actual scattering in unknown samples.
2. Theoretical description for visibility reduction in a dual-phase grating interferometer
Let us consider a standard dual-phase grating setup, as shown in Fig. 1. The first phase grating $G_1$ with a period of $p_1$ is placed at a distance of $L_1$ from the X-ray source. A second phase grating $G_2$ with a period of $p_2$ is placed at a distance $d$ from the first grating. Both gratings have a duty cycle of 0.5. The detector is placed at a distance $L$ from the source. The transmission functions of these phase gratings can be described by Fourier series as following
3. Results and discussion
In order to evaluate our model, we conducted a series of experiments demonstrating that the spectral visibility reduction can be predicted and accounted for. We started with imaging a water phantom, where visibility reduction would solely be due to beam hardening effects which our model can predict. For our second experiment, we showed that the already known real space correlation function of microspheres could be measured after the proposed correction is applied. Finally, we imaged a sample containing both areas with and without microstructure and demonstrated that with the proposed model, we could successfully distinguish those areas.
The setup we used for all the reported experiments is the one proposed in [21]. Specifically, we utilized two identical Si gratings with a period of 1.2 $\mu \textrm {m}$ and etched to a depth of 23 $\mu \textrm {m}$, introducing a $\pi$ phase shift at 17 keV. The gratings were fabricated by interference Talbot lithography and deep reactive ion etching [28]. The source we utilized was a HAMAMATSU L0101 microfocal X-ray tube, operated at 70 kVp and 100 $\mu \textrm {A}$ which resulted in a focal spot size of $9.5~\mu \textrm {m}$. The detection of the generated fringe was performed with a Princeton Instruments PI-SCX:4300 X-ray detector. The pixel size of the detector was $24~\mu \textrm {m}$ allowing to record fringes with a period down to approximately $100~\mu \textrm {m}$. The PI-SCX:4300 X-ray detector is equipped with a Gd2O2S:Tb (Gadox) scintillator, which is very efficient at energies up to 20 keV (50$\%$ efficiency). The first phase grating $G_{1}$ was placed at a distance of 0.5 m from the X-ray source. The second grating $G_{2}$ was placed at a distance of $d$ downstream $G_{1}$, and this distance was varied in order to change the autocorrelation length. Specifically, for all imaging examples, we changed $d$ from 2 mm to 18 mm in steps of 1 mm. This resulted in a range of autocorrelation lengths between 104 nm and 998 nm. The X-ray detector was placed at a distance of 0.5 m from $G_{1}$. Finally, the samples were placed right before the first phase grating $G_{1}$ resulting in an effective magnification of approximately two, which determined the effective imaging field of view (FOV) to $2.5 \times 2.5~\textrm {cm}^{2}$.
For all the corrections, we utilized the measured X-ray spectrum at the operating current and acceleration voltage as reported in the Appendix. The grating parameters were estimated from cross sectional scanning electron images (SEM) of the gratings (SEM images are shown in Appendix).
The system can be operated both in phase stepping and in single shot mode. Since in the single shot mode the image resolution is defined by the fringe period, which in turn varies across measurements at different autocorrelation lengths, we decided to work with the phase stepping mode. In this mode, one of the two gratings is laterally shifted by one or more periods in multiple steps which are smaller than the period. At each step, an image is recorded which eventually allows recording a period of the interference fringe at each pixel. For all the imaging experiments we utilized five phase steps with an exposure time of 60 sec each. The retrieval of the three complementary signals is identical to the methods used in standard grating interferometry when phase stepping is used [29].
3.1 Water phantom
To demonstrate that beam hardening due to the sample can result in an observed visibility reduction even in lack of small angle X-ray scattering events, we imaged an absorbing sample (non-scattering) at different inter-grating distances. As a sample, we chose a Poly(methyl methacrylate) (PMMA) cuvette filled with water; the width of the cuvette was 0.5 cm. In Figs. 2(a), (b), and (c) the retrieved visibility reduction images for autocorrelation lengths of 104 nm, 525 nm, and 998 nm are shown respectively. A signal both for the PMMA container and the water is present. Nevertheless, no microstructure is present in the sample, and therefore this varying signal is due to the beam hardening as explained by our theory. Given the spectrum of the X-ray source (shown in the Appendix), the detector efficiency, the grating parameters, and the attenuation coefficient of water and PMMA, we can theoretically calculate the expected visibility reduction of the sample. In Fig. 2(d) the average visibility reduction for the marked area is compared to the theoretical calculation for the range of probed autocorrelation lengths. In general, a strong agreement is observed validating our theoretical model. Nonetheless, small discrepancies are observed due to uncertainties regarding the exact parameters utilized such as the X-ray spectrum and grating parameters.
3.2 Glass microspheres
As a next step, we wanted to demonstrate that quantitative visibility reduction due to small angle X-ray scattering can be retrieved by correcting for beam hardening effects. As samples, we utilized glass spheres of known diameters and volume fractions whose small angle scattering signals can be estimated. Two different microsphere powders [30] with diameters of 0.26 $\mu \textrm {m}$ and 1.7 $\mu \textrm {m}$ were placed on top of each other in a cylindrical container as shown in Fig. 3(a). These samples were also utilized for the demonstration of the tunability of the dual-phase grating interferometer in [21]. In light of our theoretical analysis, we will reinterpret the previous experimental results in order to demonstrate quantitativeness. In Figs. 3(b) and (c) representative visibility reduction images at two different autocorrelation lengths are shown. The two powders have the same absorption coefficient which would lead to similar beam hardening; therefore, the difference in the observed visibility reduction between the two powders is due to their difference in the size of microstructures. In Fig. 3(d) the average visibility reduction for the measured autocorrelation length range in the two marked areas is shown.
Given that the material of the samples is known, we can estimate the spectral term $\left (\frac {V}{V_0}\right )_{spectral}$ which we use to correct the experimental data. This results in the curves shown in Fig. 3(e) which are compared to the expected visibility reduction curves from the numerical calculation of the small angle scattering from SASview [31]. In the Appendix, the scattering intensity profiles for both powders are shown, these scattering intensities can be translated to visibility reduction by applying the Hankel transform as explained here [18]. The calculated and experimental curves show a high accuracy which corroborates our model. As explained in the previous experiment a few deviations are present. In particular, the high frequency oscillations in the theoretical model of the $0.26~\mu \textrm {m}$ powder are blurred out in the experimental curve. This is due to the polychromaticity of the utilized source which is not modeled in the SAXS calculations.
3.3 Generalised sample
In this section, we show the applicability of the model for a sample whose sample properties (absorption coefficient, microstructured or not) are unknown. As samples we used, a glass block and polydispersed glass powder of unknown granularity (see Fig. 4(a) for a sketch). Assuming, the sample is made of single material, we can calculate the absorption coefficient. The absorption coefficient of the samples was estimated as explained in the Appendix in subsection 5.3. Both samples had a thickness of 0.5 cm. The visibility reduction images of the samples at two different autocorrelation lengths are shown in Figs. 4(c) and (e). The corresponding images obtained after correction are shown in Figs. 4(b) and (d). The average visibility reduction within the marked areas of Fig. 4(b) of the uncorrected and corrected data is shown in Figs. 5(a) and (b) for the glass powder and the glass block respectively. Interpreting the uncorrected visibility reduction signals would lead us to conclude that both samples contain subpixel structures since they both show a visibility reduction. Furthermore, one would possibly claim that the glass powder contains structures finer than 104 nm since the visibility reduction signal seems to have saturated. On the other hand, one would conclude that the glass block has structures in the range of the autocorrelation length since a signal fluctuation is observed. Alas, both statements are erroneous. After correction, the visibility reduction curves dramatically change as shown in Figs. 5(a) and (b). The corrected visibility reduction for the glass block seems to be constant and equal to one, meaning that no unresolved microstructures are present which is the ground truth. However, we see a strong visibility reduction for the glass spheres which is changing as the autocorrelation length is increased. Eventually, this curve could be fitted with an appropriate model in order to retrieve structural information.
4. Conclusion
Like conventional grating interferometers, the dual-phase grating interferometer suffers from beam hardening effects when operated with polychromatic sources. These effects cause visibility reduction to be detected even in the absence of small angle scattering. Through theoretical modeling and experimental investigations, we have demonstrated that this effect strongly depends on the distance between the two phase gratings, consequently affecting measurements across different autocorrelation lengths which are performed by changing the inter-grating distance. Therefore, having a deleterious effect on any statements made on alleged micro-structural properties of the sample. With a reasonable pre knowledge of the system parameters (such as spectrum and grating geometry), a rigorous theoretical model allows for mitigation of this effect. Specifically, we have demonstrated that the measured visibility reduction can be written as the product of two terms. The first one purely depends on the beam hardening introduced by the sample, and the second corresponds to the small angle X-ray scattering of the microstructure. Our model allows us to predict the spectral term and eventually retrieve the true scattering signal. Our description was validated with a series of experiments ranging from basic proof of principle of the beam hardening induced visibility reduction to the quantitative retrieval of the real space correlation function of microparticles. To conclude, we have demonstrated that quantitative retrieval of the underlying real space autocorrelation function can be achieved with polychromatic dual-phase X-ray grating interferometers. This comprises a significant step towards the application of the method for the quantitative investigation of the (sub)microstructure of materials and devices. Taking into account the advances in grating fabrication technologies [32,33], the method could be translated to higher energy ranges, hence allowing the investigation of devices such as Polymer Electrolyte Fuel Cells (PEFCs) under operation. Such fuel cells contain multiple layers whose pore sizes distributions affect water transport and eventually their performance. With our method, we could get a close look at the water propagation through the submicrometer porous structures which could lead to the design optimization of the fuel cell.
5. Appendix
5.1. Evaluating the wave equations
5.1.1. Wave propagation
We consider that the sample is placed right in front of grating $G_1$ as shown in Fig. 1. The complex refractive index of the sample can is written as $1 - \delta + i\beta$. The electric field for a given wavelength $\lambda$ after the sample is given by $E_0 \exp [ \alpha _s \left (x_1,y_1\right ) /2 ]\exp [i\phi \left (x_1,y_1\right )]$ according to the projection approximation, where $E_0$ is the amplitude of the incident X-rays, and $\alpha _{s}$ and $\phi$ are given by $\frac {2 \pi \delta \left (x,y\right )}{\lambda }$ and $\frac {4 \pi \beta \left (x_{1},y_{1}\right )}{\lambda }$ respectively, where $(x_{1},y_{1})$ are the coordinates at the sample and $G_{1}$ plane. Additionally, we assume that the sample has only fine phase variations that cannot be resolved by the detector. Initially, we propagate the multiplied transmission functions of the $G_{1}$ and the sample. By utilising the Fresnel propagator. This results in the field at the second grating $G_{2}$ being equal to
5.1.2. Absorption contrast
To calculate the transmission of a sample, we first acquire a flat image without a sample in the beam, the DC term of the fringe is then given by
5.1.3. Visibility reduction contrast
The dark-field signal is retrieved from the visibility of the interference pattern which corresponds to the modulation of the oscillating term which is given by
5.2. Deriving effective wavelength and autocorrelation function
To calculate the autocorrelation length, we first need to calculate the effective wavelength. By starting from Eq. (40) and assuming that the sample is weakly absorbing we get
5.3. Estimating absorption contrast
In order to estimate the absorption coefficient $\mu$ of the sample, we assume that it varies with the third power of the wavelength $\lambda$, such that $\mu = r\lambda ^{4}+ q \lambda ^{3}$ is a valid approximation [35]. Eventually, we would like to find a way to estimate the parameters: $r$ and $q$. For a given spectrum $s(\lambda )$ and a sample thickness $T$ the transmission of the sample can be calculated as
whereWhen there are one or more absorption edges the estimation of attenuation coefficient may not be very accurate in the absorption edge regions (e.g. K edge). Other attenuation coefficient models could be used for accurate estimation around the edge region. However, such models require additional parameters like atomic number. Given that additional information about the sample like atomic number and density are known, such models could be used to estimate wavelength dependent absorption coefficient [35].
5.4. Theoretical SAXS curves
We were able to obtain the SAXS curves shown in Fig. 7 for Silicon dioxide microsphere samples using SASview. These curves were transformed into visibility reduction curves, as described in [18].
5.5. SEM images of the gratings used in the experiments
We used two phase gratings of period 1.2 $\mu m$, with duty cycle of $0.5$ and etched to depth of $23 \mu m$. The SEM image of the gratings is shown in Fig. 8.
Funding
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (Sinergia Grant CRSII5_18356 Clinical GI-BCT); Eurostars (E!1106 INFORMAT); Paul Scherrer Institut (PSI-CROSS Proposal).
Acknowledgments
We thank our lab technician Gordan Mikuljan for his technical support.
Disclosures
The authors declare that there are no conflicts of interest related to this article.
References
1. C. David, B. Nöhammer, H. Solak, and E. Ziegler, “Differential x-ray phase contrast imaging using a shearing interferometer,” Appl. Phys. Lett. 81(17), 3287–3289 (2002). [CrossRef]
2. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of x-ray talbot interferometry,” Jpn. J. Appl. Phys. 42(Part 2, No. 7B), L866–L868 (2003). [CrossRef]
3. F. Pfeiffer, M. Bech, O. Bunk, P. Kraft, E. F. Eikenberry, C. Brönnimann, C. Grünzweig, and C. David, “Hard-X-ray dark-field imaging using a grating interferometer,” Nat. Mater. 7(2), 134–137 (2008). [CrossRef]
4. W. Yashiro, Y. Terui, K. Kawabata, and A. Momose, “On the origin of visibility contrast in x-ray talbot interferometry,” Opt. Express 18(16), 16890–16901 (2010). [CrossRef]
5. S. K. Lynch, V. Pai, J. Auxier, A. F. Stein, E. E. Bennett, C. K. Kemble, X. Xiao, W.-K. Lee, N. Y. Morgan, and H. H. Wen, “Interpretation of dark-field contrast and particle-size selectivity in grating interferometers,” Appl. Opt. 50(22), 4310–4319 (2011). [CrossRef]
6. M. Strobl, “General solution for quantitative dark-field contrast imaging with grating interferometers,” Sci. Rep. 4(1), 7243 (2015). [CrossRef]
7. T. H. Jensen, M. Bech, I. Zanette, T. Weitkamp, C. David, H. Deyhle, S. Rutishauser, E. Reznikova, J. Mohr, R. Feidenhans, and F. Pfeiffer, “Directional x-ray dark-field imaging of strongly ordered systems,” Phys. Rev. B 82(21), 214103 (2010). [CrossRef]
8. T. Lauridsen, E. M. Lauridsen, and R. Feidenhans, “Mapping misoriented fibers using x-ray dark field tomography,” Appl. Phys. A 115(3), 741–745 (2014). [CrossRef]
9. M. Kagias, Z. Wang, M. E. Birkbak, E. Lauridsen, M. Abis, G. Lovric, K. Jefimovs, and M. Stampanoni, “Diffractive small angle X-ray scattering imaging for anisotropic structures,” Nat. Commun. 10(1), 5130 (2019). [CrossRef]
10. F. Yang, F. Prade, M. Griffa, I. Jerjen, C. Di Bella, J. Herzen, A. Sarapata, F. Pfeiffer, and P. Lura, “Dark-field x-ray imaging of unsaturated water transport in porous materials,” Appl. Phys. Lett. 105(15), 154105 (2014). [CrossRef]
11. F. Prade, A. Yaroshenko, J. Herzen, and F. Pfeiffer, “Short-range order in mesoscale systems probed by x-ray grating interferometry,” Europhys. Lett. 112(6), 68002 (2015). [CrossRef]
12. S. Senck, M. Scheerer, V. Revol, K. Dobes, B. Plank, and J. Kastner, “Non-destructive evaluation of defects in polymer matrix composites for aerospace applications using x-ray talbot-lau interferometry and micro ct,” in 58th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, (2017), p. 0355.
13. T. Michel, J. Rieger, G. Anton, F. Bayer, M. W. Beckmann, J. Durst, P. A. Fasching, W. Haas, A. Hartmann, G. Pelzer, M. Radicke, C. Rauh, A. Ritter, P. Sievers, R. Schulz-Wendtland, M. Uder, D. L Wachter, T. Weber, E. Wenkel, and A. Zang, “On a dark-field signal generated by micrometer-sized calcifications in phase-contrast mammography,” Phys. Med. Biol. 58(8), 2713–2732 (2013). [CrossRef]
14. N. Hauser, Z. Wang, R. A. Kubik-Huch, M. Trippel, G. Singer, M. K. Hohl, E. Roessl, T. Köhler, U. van Stevendaal, N. Wieberneit, and M. Stampanoni, “A study on mastectomy samples to evaluate breast imaging quality and potential clinical relevance of differential phase contrast mammography,” Invest. Radiol. 49(3), 131–137 (2014). [CrossRef]
15. K. Scherer, E. Braig, K. Willer, M. Willner, A. A. Fingerle, M. Chabior, J. Herzen, M. Eiber, B. Haller, M. Straub, H. Schneider, E. J. Rummeny, P. B. Noël, and F. Pfeiffer, “Non-invasive differentiation of kidney stone types using x-ray dark-field radiography,” Sci. Rep. 5(1), 9527 (2015). [CrossRef]
16. A. Yaroshenko, K. Hellbach, A. Ö. Yildirim, T. M. Conlon, I. E. Fernandez, M. Bech, A. Velroyen, F. G. Meinel, S. Auweter, M. Reiser, O. Eickelberg, and F. Pfeiffer, “Improved in vivo assessment of pulmonary fibrosis in mice using x-ray dark-field radiography,” Sci. Rep. 5(1), 17492 (2015). [CrossRef]
17. K. Willer, A. A. Fingerle, L. B. Gromann, F. De Marco, J. Herzen, K. Achterhold, B. Gleich, D. Muenzel, K. Scherer, M. Renz, B. Renger, F. Kopp, F. Kriner, F. Fischer, C. Braun, S. Auweter, K. Hellbach, M. F. Reiser, T. Schroeter, J. Mohr, A. Yaroshenko, H.-I. Maack, T. Pralow, H. van der Heijden, R. Proksa, T. Koehler, N. Wieberneit, K. Rindt, E. J. Rummeny, F. Pfeiffer, and P. B. Noël, “X-ray dark-field imaging of the human lung–A feasibility study on a deceased body,” PLoS One 13(9), e0204565 (2018). [CrossRef]
18. R. P. Harti, M. Strobl, B. Betz, K. Jefimovs, M. Kagias, and C. Grünzweig, “Sub-pixel correlation length neutron imaging: Spatially resolved scattering information of microstructures on a macroscopic scale,” Sci. Rep. 7(1), 44588 (2017). [CrossRef]
19. S. Gkoumas, P. Villanueva-Perez, Z. Wang, L. Romano, M. Abis, and M. Stampanoni, “A generalized quantitative interpretation of dark-field contrast for highly concentrated microsphere suspensions,” Sci. Rep. 6(1), 35259 (2016). [CrossRef]
20. A. J. Brooks, D. S. Hussey, H. Yao, A. Haghshenas, J. Yuan, J. M. LaManna, D. L. Jacobson, C. G. Lowery, N. Kardjilov, S. Guo, M. M. Khonsari, and L. G. Butler, “Neutron interferometry detection of early crack formation caused by bending fatigue in additively manufactured SS316 dogbones,” Mater. Des. 140, 420–430 (2018). [CrossRef]
21. M. Kagias, Z. Wang, K. Jefimovs, and M. Stampanoni, “Dual phase grating interferometer for tunable dark-field sensitivity,” Appl. Phys. Lett. 110(1), 014105 (2017). [CrossRef]
22. H. Miao, A. Panna, A. A. Gomella, E. E. Bennett, S. Znati, L. Chen, and H. Wen, “A universal moiré effect and application in x-ray phase-contrast imaging,” Nat. Phys. 12(9), 830–834 (2016). [CrossRef]
23. W. Yashiro, P. Vagovič, and A. Momose, “Effect of beam hardening on a visibility-contrast image obtained by x-ray grating interferometry,” Opt. Express 23(18), 23462–23471 (2015). [CrossRef]
24. A. Yan, X. Wu, and H. Liu, “Beam hardening correction in polychromatic x-ray grating interferometry,” Opt. Express 25(20), 24690–24704 (2017). [CrossRef]
25. D. S. Hussey, H. Miao, G. Yuan, D. Pushin, D. Sarenac, M. G. Huber, D. L. Jacobson, J. M. LaManna, and H. Wen, “Demonstration of a white beam far-field neutron interferometer for spatially resolved small angle neutron scattering,” arXiv preprint arXiv:1606.03054 (2016).
26. M. Chabior, T. Donath, C. David, O. Bunk, M. Schuster, C. Schroer, and F. Pfeiffer, “Beam hardening effects in grating-based x-ray phase-contrast imaging,” Med. Phys. 38(3), 1189–1195 (2011). [CrossRef]
27. 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(18), 23142–23155 (2018). [CrossRef]
28. K. Jefimovs, L. Romano, J. Vila-Comamala, M. Kagias, Z. Wang, L. Wang, C. Dais, H. Solak, and M. Stampanoni, “High-aspect ratio silicon structures by displacement talbot lithography and bosch etching,” in Advances in Patterning Materials and Processes XXXIV, vol. 10146 (International Society for Optics and Photonics, 2017), p. 101460L.
29. 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(16), 6296–6304 (2005). [CrossRef]
30. See http://www.cospheric.com for one-stop reliable source of high quality microspheres. http://www.cospheric.com. Accessed: 2019-09-11.
31. M. Doucet, J. H. Cho, G. Alina, J. Bakker, W. Bouwman, P. Butler, K. Campbell, M. Gonzales, R. Heenan, A. Jackson, P. Juhas, S. King, P. Kienzle, J. Krzywon, A. Markvardsen, T. Nielsen, L. O’Driscoll, W. Potrzebowski, R. Ferraz Leal, T. Richter, P. Rozycko, and A. Washington, “Sasview version 4.1,” Zenodo and http://www.sasview.org (2017).
32. L. Romano, M. Kagias, J. Vila-Comamala, K. Jefimovs, L.-T. Tseng, V. A. Guzenko, and M. Stampanoni, “Metal assisted chemical etching of silicon in gas phase: a nanofabrication platform for X-ray optics,” Nanoscale Horiz. 5(5), 869–879 (2020). [CrossRef]
33. M. Kagias, Z. Wang, V. A. Guzenko, C. David, M. Stampanoni, and K. Jefimovs, “Fabrication of Au gratings by seedless electroplating for X-ray grating interferometry,” Mater. Sci. Semicond. Process. 92, 73–79 (2019). [CrossRef]
34. Y. I. Nesterets, “On the origins of decoherence and extinction contrast in phase-contrast imaging,” Opt. Commun. 281(4), 533–542 (2008). [CrossRef]
35. L. Gerward, “Analytical approximations for X-ray attenuation coefficients,” Nucl. Instrum. Methods 181(1-3), 11–14 (1981). [CrossRef]