Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Two-photon X-ray ghost microscope

Open Access Open Access

Abstract

This article presents a non-classical imaging mechanism that produces a diffraction-limited and magnified ghost image of the internal structure of an object through the measurement of intensity fluctuation correlation formed by two-photon interference. In principle, the lensless X-ray ghost imaging mechanism may achieve a spatial resolution determined by the wavelength and the angular diameter of the X-ray source, ∼λθs, with possible reduction caused by additional optics. In addition, it has the ability to image select “slices” deep within an object, which can be used for constructing 3D view of its internal structure.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Due to the properties of X rays, such as higher penetrating power and shorter wavelengths compared to visible light, the field of X-ray microscopy has attracted a lot of research and has a wide range of applications in various fields of study [13]. As a result, many different types of techniques have been developed that achieve micron and nano-level spatial resolution [46]. A well known type of X-ray imaging is standard projectional radiography with hard X-rays ($> 10$ keV) due to its simplicity and common use in the medical field [7,8]. This technique is also commonly used in X-ray microscopy, including the slightly adjusted form of projection-based phase contrast imaging [9]. The resolution of this projection is closely dependent on how small of a source is used, meaning a “perfect,” magnified image can be obtained with a “perfect” point source. With the help of the state-of-the-art techniques, focusing the X-ray beam down to 30 nm has been achieved; producing a projection with an equivalent resolution [10]. Unlike lens-based imaging mechanisms, this technique cannot take advantage of the wave nature of radiation to achieve diffraction-limited point-to-point imaging resolution.

Focusing optics such as imaging lenses play a critical role in producing a diffraction-limited point-to-point correlation between the object plane and the image plane, forming a magnified or demagnified image of the object [1113]. If desired, an additional lens system is then able to map the primary image onto a secondary image plane for further magnification or demagnification; notably making an optical microscope possible. While lens-based optical microscopes are commonly used to obtain high-resolution images, the relatively longer wavelengths of visible light limits this to $\sim$200 nm resolution [11,12]. Unfortunately, unlike visible light, X rays experience a refractive index close to 1 for many materials such as those used in traditional imaging lenses, making such materials useless for focusing the X-ray beam [7,8]. To overcome this, some alternative focusing X-ray optics have been developed such as compound refractive lenses, focusing mirrors, and zone plates [1416]. While zone plates and focusing mirrors are typically limited to soft X rays ($< 10$ keV), compound refractive lenses can be designed for hard X rays ($> 10$ keV). Unlike projectional radiography, which is a projection, or “shadow,” of the X rays that pass through the object, focusing X-ray optics produce high-resolution images by taking full advantage of the short wavelength of the X rays used. Most notably, this has been demonstrated by producing 10 nm resolution images with soft X rays [16]. We would expect even higher resolution with higher energy levels but limitations in the fabrication of the focusing X-ray optics are the limiting factor. Even as advances are made in focusing X-ray optics, it is difficult for classic imaging devices such as these to obtain magnified, high-resolution images of the interior structure of a thicker object. Recording the resulting image on a photodetector array relies on magnification produced by these focusing X-ray optics or an optical lens following the scintillator. Due to the fact that object planes further away from the focal plane do not achieve as much magnification, the magnification may not be enough for a sensor array to capture the desired resolution.

For the sake of comparison, it is worth noting that other X-ray microscopy techniques exist. This mainly includes coherent X-ray diffraction imaging (CXDI) and a similar technique of X-ray ptychography, both of which can theoretically provide diffraction-limited resolution without the need of a lens [6,17,18]. These techniques utilize the measurement of a diffraction pattern following the object that then requires post processing to recover the image. Similar to these techniques is lensless Fourier-transform imaging which also requires post processing [19]. Projection-based X-ray phase contrast imaging was already mentioned, but the larger, encompassing field of phase contrast microscopy is well studied [20,21].

In this article, we discuss the physics of two-photon ghost imaging to produce a diffraction-limited, lensless, image-forming correlation between the image plane and object plane, for which X rays allow imaging of the internal structure of the object [2225]. Due to the diffraction-limited nature of this imaging process, sub-nanometer resolution can theoretically be achieved but will be limited by physical constraints of the imaging system. One key example of this is the fact that the ghost imaging process itself does not provide any magnification to allow for photodetectors (pixels of a CCD or CMOS) with relatively significant physical size to capture the full resolving power of the imaging process. To account for this, we propose the help of a secondary imaging device, either focusing X-ray optics or a scintillator paired with a visible-light lens assembly. With secondary imaging, the primary lensless ghost image can be mapped onto a secondary image plane with significant magnification to be resolvable by a standard CCD or CMOS. Although the addition of secondary imaging will introduce its own restrictions on resolution (among other fabrication restrictions) compared to X-ray ghost imaging on its own, we may label this setup as an X-ray ghost microscope as it produces a magnified, high-resolution image of the internal structure of an object. One interesting result of the X-ray ghost microscope is that the secondary imaging system is focused on the primary ghost image rather than the physical object. This removes the common constraint encountered in classic X-ray microscopy with focusing X-ray optics: in order to get high magnification from the setup it is desired to have the object as close to the focal point as possible, but the short focal lengths of focusing X-ray optics makes it impossible to do this for the internal structure of larger (thicker) objects [16].

To avoid confusion with other discussions of ghost imaging, it is worth discussing the different types of ghost imaging that have been developed. Visible-light ghost imaging was first demonstrated in 1995 using two-photon interference from entangled photon pairs through the measurement of the second-order coherence function [26]. In addition, this has been demonstrated recently with entangled X-ray photons [27]. 10 years after the initial, 1995 demonstration, it was demonstrated that two-photon interference of randomly created and randomly paired photons in a thermal state can produce a similar point-to-point image-forming correlation without the need of an imaging lens, commonly called lensless ghost imaging [2225]. Interestingly, the first lensless ghost imaging of thermal light was observed in the secondary ghost imaging plane with the help of a lens, providing the inspiration for this paper. An unfolded schematic diagram of the 2005 ghost imaging experiment of thermal light is illustrated in Fig. 1. X-ray ghost imaging has now been demonstrated multiple times [2835]; however, most demonstrations did not take advantage of two-photon interference to produce an image. The most commonly used technique for current X-ray ghost imaging is to introduce a spatially varied material in the X-ray beam to produce an artificial “speckle” pattern instead of quantum speckles inherently present in the beam through two-photon interference. This results in a classical speckle-to-speckle ghost imaging with resolution dependent on the artificial “speckle” size as opposed to being diffraction-limited.

 figure: Fig. 1.

Fig. 1. Experimentally achieved with a beam splitter, this “unfolded” schematic of the 2005 thermal light ghost imaging demonstration of Valencia et al. helps depict the symmetry of the two-photon amplitudes. The one-to-one lensless ghost image and the magnified secondary ghost image are results of two-photon interference. Of the many two-photon amplitudes, such as the red and blue two-photon amplitudes of the random $m$th and $n$th pair of photons shown here, each superpose constructively at a corresponding $\vec {\rho }_{i}$ to from an image.

Download Full Size | PDF

Classified according to their experimental setup and working mechanism, we can group the two classes of ghost imaging: (I) The observed ghost image is produced from a natural point-to-point image forming correlation that is the result of two-photon interference. This class of ghost imaging follows the mechanism of the entangled state ghost imaging of Pittman et al. [26] and the thermal light ghost imaging of Valencia et al. [22]. With thermal light ghost imaging, the point-to-point image-forming correlation is able to achieve a diffraction-limited spatial resolution of $\lambda / \Delta \theta _s$, where $\Delta \theta _s$ is the angular diameter of the light source, and also has the potential to be turbulence-free [13,3638]. Pelliccia et al. have used type (I) ghost imaging in their first demonstration of X-ray ghost imaging [29]. Alternatively, (II) the observed ghost image is produced from a classical speckle-to-speckle correlation shown with visible light [3941] and with a variety of X-ray sources [30,3235]. Here, a ghost image of the object is obtained from the coincidences between two sets of identical, artificial “speckles” formed by either spatially correlated laser beams or an aperture mask following the light source to produce shadows and bright spots distributed on the object plane and on the ghost image (detector) plane. Included as a classical speckle-to-speckle correlation is computational ghost imaging which removes the need of a beam splitter by measuring the object plane and ghost image plane at separate times with an identical speckle pattern [28,31,42]. These artificial speckles are different from the two-photon interference that produces a point-to-point image-forming correlation and, correspondingly, this type of ghost imaging has a resolution dependent on the size of the classically formed speckles. Through the correlation measurement, this type of ghost imaging observes a projection or a shadow of the object, comparable to how classic X-ray imaging technology is a projection of the object. Both classes of ghost imaging use the measurement of intensity fluctuation correlation to obtain an image, utilizing either changes in the artificial speckle distribution over time or fluctuations resulting from two-photon interference.

The X-ray ghost microscope studied in this article belongs to type (I) ghost imaging via two-photon interference. As with lensless ghost imaging, the spatial resolution of X-ray ghost imaging before magnification is determined by $\lambda / \Delta \theta _s$, where $\lambda$ is the wavelength and $\Delta \theta _s$ is the angular diameter of the X-ray source. This unique characteristic of lensless ghost imaging makes the resolving potential of an X-ray ghost imaging a point of interest. As an example, ignoring potential physical constraints, using a high-energy ($>$20 keV) X-ray source with a relatively large angular diameter may produce a one-to-one (no magnification or demagnification) lensless ghost image with sub-nanometer resolution. It should be noted that achieving higher resolution with a source with a larger angular diameter is the opposite of what is desired for projectional X-ray imaging, phase contrast imaging, and CXDI which all see an increase in resolution when using a smaller, point-like source, such as modern X-ray tubes and synchrotron X-ray sources [7,43]. Type (I) ghost imaging is still possible with these small sources, but the resolution will be reduced. When expecting a ghost image with even sub-micrometer resolution, the image becomes unresolvable by any state-of-the-art CCD or CMOS detector arrays. However, this primary ghost image can be magnified significantly onto a secondary ghost image plane through either focusing X-ray optics or a scintillator-lens system, making sub-micron (potentially even 10 nm [16]) sized features resolvable by a CCD or CMOS with micrometer pixels. Importantly, the secondary ghost image is produced directly by the point-to-point image-forming correlation between the object plane, $\vec {\rho }_o$, and the secondary ghost image plane, $\vec {\rho }_i$, where $\vec {\rho }_o$ and $\vec {\rho }_i$ represent the transverse coordinates of the object plane and image plane, respectively. An additional benefit of two-photon interference is that this X-ray ghost microscope can be set up such that it is insensitive to any rapid phase variations along the optical path due to random changes in composition, density, length, index of refraction, or medium vibration, namely “turbulence-free” [3638]. Overcoming vibrations is especially important for the extremely high resolution imaging obtainable with the X-ray ghost microscope.

This article is organized in four sections: Sec. 2 we introduce the concept of two-photon interference induced intensity fluctuation correlation, Sec. 3 we show a point-to-point image-forming correlation for producing a lensless ghost image, Sec. 4 we extend this ghost image to a secondary image using either focusing X-ray optics or a scintillator-lens system, and Sec. 5 we discuss how to obtain observable intensity fluctuation correlation when using an X-ray source with a broad spectrum.

2. Two-photon interference induced X-ray intensity fluctuation correlation

Recall, in 1905, Einstein introduced a granularity to radiation, abandoning the continuum interpretation of Maxwell [44]. This led to a microscopic picture of radiation and a statistical view of light, including X rays. In Einstein’s picture, a light source consists of many point-like sub-sources, each of which emit their own subfields. Originally labeled in German by “strahlenbündel,” which translates to “bundle of ray” in English, we now label these bundles of light as subfields or the quantum term of photons. For a thermal light source, these subfields (photons) are emitted in a random manner such that the $m$th subfield emitted from the $m$th point-like sub-source may propagate in all possible directions with a random initial phase. It has been shown that an effective wavefunction can be defined from the quantum theory of optical coherence [45,46] to specify the space-time behavior of a photon [13,47]. The effective wavefunction of a photon in a thermal state is mathematically the same function as Einstein’s subfield. In Einstein’s picture, the radiation measured at coordinate $(\boldsymbol{r}, t)$ is the result of a superposition of a large number of subfields,

$$E(\boldsymbol{r}, t) = \sum_{m}E_m(\boldsymbol{r}, t) = \sum_{m}E_m(\boldsymbol{r}_m, t_m) \, g_m(\boldsymbol{r}_m, t_m; \boldsymbol{r}, t),$$
where $E_m(\boldsymbol{r}_m, t_m)$ labels the subfield emitted from the $m$th sub-source at coordinate ($\boldsymbol{r}_m, t_m$), and $g_m(\boldsymbol{r}_m, t_m; \boldsymbol{r}, t)$ represents the field propagator or Green’s function that propagates the $m$th subfield from coordinate $(\boldsymbol{r}_m, t_m)$ to coordinate $(\boldsymbol{r}, t)$. Shortening the notation, we will replace $E_m(\boldsymbol{r}_m, t_m)$ with $E_m$ and $g_m(\boldsymbol{r}_m, t_m; \boldsymbol{r}, t)$ with $g_m(\boldsymbol{r}, t)$. In Einstein’s picture of light, the expectation value of the intensity corresponds to a statistical ensemble average which takes into account all possible realizations of the field or, more specifically, takes into account all possible relative phases of the subfields:
$$\begin{aligned} \langle I(\boldsymbol{r}, t) \rangle & = \langle E^*(\boldsymbol{r}, t) E(\boldsymbol{r}, t) \rangle \\ & = \big{\langle} \sum_{m} E^*_m(\boldsymbol{r}, t) \sum_n E_n(\boldsymbol{r}, t) \big{\rangle} \\ & = \big{\langle} \sum_{m} \big{|} E_m(\boldsymbol{r}, t) \big{|}^2 \big{\rangle} + \big{\langle} \sum_{m\neq n} E^*_m(\boldsymbol{r}, t) E_n(\boldsymbol{r}, t) \big{\rangle} \\ & = \sum_{m} \big{|} E_m(\boldsymbol{r}, t) \big{|}^2. \end{aligned}$$
The expectation value of $\sum _{m\neq n} E^*_m(\boldsymbol{r}, t) E_n(\boldsymbol{r}, t)$ goes to zero when taking into account all possible relative phases of the subfields. We may conclude the theoretical expectation value of intensity $\langle I(\boldsymbol{r}, t) \rangle = \sum _{m} \big {|} E_m(\boldsymbol{r}, t) \big {|}^2$ is the result of the $m$th subfield interfering with the $m$th subfield itself, while the $m\neq n$ term is a result of the $m$th subfield interfering with a different $n$th subfield. In a realistic measurement in which only a limited number of subfields contribute to the measurement, all possible phases may not be present, meaning the $m\neq n$ term may not vanish and contribute noise to the measurement. We may name this term a two-photon interference induced intensity fluctuation,
$$\Delta I(\boldsymbol{r}, t) = \sum_{m\neq n} E^*_m(\boldsymbol{r}, t) E_n(\boldsymbol{r}, t).$$
The two-photon induced intensity fluctuation is different from other traditional “intensity fluctuations.” For instance, the number of subfields contributing to a measurement may vary from measurement to measurement and thus the value of $\sum _{m} \big {|} E_m(\boldsymbol{r}, t) \big {|}^2$ may change from measurement to measurement. These variations may be from a classical mask in the path of the light or variations in the intensity of the source. These would be the fluctuations used in type (II) ghost imaging. Note that in the following discussion, no classical intensity fluctuations are involved, neither from a spatial intensity distribution nor temporal intensity distribution, but only taking the two-photon induced intensity fluctuations of thermal light.

The X-ray ghost microscope studied in this article utilizes the two-photon induced intensity fluctuation correlation $\langle \Delta I(\boldsymbol{r}_1, t_1) \Delta I(\boldsymbol{r}_2, t_2) \rangle$ for type (I) ghost imaging. Although the expectation or ensemble average of the intensity fluctuation measured by a single detector $D_j$, $j = 1, 2$, is zero, $\langle \Delta I(\boldsymbol{r}_j, t_j) \rangle = \langle \sum _{m\neq n} E^*_m(\boldsymbol{r}_j, t_j) E_n(\boldsymbol{r}_j, t_j) \rangle = 0$, the expectation or ensemble average of the correlation of the intensity fluctuations measured by $D_1$ and $D_2$, jointly, may not equal zero,

$$\begin{aligned} \langle \Delta I(\boldsymbol{r}_1, t_1) \Delta I(\boldsymbol{r}_2, t_2) \rangle & = {\Big{\langle}} \sum_{m\neq n} E_m^*(\boldsymbol{r}_1, t_1) E_{n}(\boldsymbol{r}_1, t_1) \sum_{p\neq q} E_p^*(\boldsymbol{r}_2, t_2) E_{q}(\boldsymbol{r}_2, t_2) {\Big{\rangle}} \\ & = \sum_{m\neq n} E_m^*(\boldsymbol{r}_1, t_1) E_{n}(\boldsymbol{r}_1, t_1) E_n^*(\boldsymbol{r}_2, t_2) E_{m}(\boldsymbol{r}_2, t_2). \end{aligned}$$
Due to the random relative phases between the subfields canceling in a specific case, when $m = q$ and $n = p$, there is a surviving term in the above summation. Mathematically, the result of Eq. (4) can be represented as the cross term of the following superposition,
$$\begin{aligned} G^{(2)}(\boldsymbol{r}_1, t_1; \boldsymbol{r}_2, t_2) & = \sum_{m \neq n} \big{|} E_m(\boldsymbol{r}_1, t_1) E_n(\boldsymbol{r}_2, t_2) + E_n(\boldsymbol{r}_1, t_1) E_m(\boldsymbol{r}_2, t_2) \big{|}^2 \\ & = \langle I(\boldsymbol{r}_1, t_1) \rangle \langle I(\boldsymbol{r}_2, t_2) \rangle + \langle \Delta I(\boldsymbol{r}_1, t_1) \Delta I(\boldsymbol{r}_2, t_2) \rangle, \end{aligned}$$
corresponding to the superposition of two different yet indistinguishable alternatives of joint photodetection: (1) the $m$th subfield (photon) is measured at $(\boldsymbol{r}_1, t_1)$ while the $n$th subfield (photon) is measured at $(\boldsymbol{r}_2, t_2)$; (2) the $n$th subfield (photon) is measured at $(\boldsymbol{r}_1, t_1)$ while the $m$th subfield (photon) is measured at $(\boldsymbol{r}_2, t_2)$. The cross term of this superposition models the concept of randomly paired photons interfering with the pair itself, namely two-photon interference.

3. Lensless ghost imaging with X-rays

To better visualize the X-ray ghost microscope it is best to start with the working mechanism of X-ray lensless ghost imaging. This simple experimental setup, schematically illustrated in Fig. 2, consists of an X-ray beam splitter [48, 49] that divides the X-ray beam from a disk-like source into two beams. A 2D array of X-ray detectors, $D_1$, is placed in beam-one at a selected plane of $z_1 = z_{gi}$ in the Fresnel near-field. Following the output of beam-two, we place an object followed by a bucket X-ray photodetector, $D_2$, which collects all X rays transmitted from the object. The Green’s function that propagates a subfield from point $(\rho_0,z_0,t_0)$ to point $(\rho, z, t)$ is thus

$$g_{0}(\rho,z,t)= \frac{1}{r}\, e^{-i (\omega t - k r)},$$
where $r = \sqrt {(z^2 + |\rho - \rho_0|^2)}$. In the following, we use the Fresnel near-field approximation to provide an analytical solution for a proof-of-principle analysis. In certain cases, the Fresnel approximation may not be accurate enough representing the spatial correlation and a numerical solution may be necessary. In the Fresnel approximation, when $|\vec {\rho } - \vec {\rho }_0|^2 \ll z^2$ we take the first-order expansion of $r$ in terms of $z$ and $\vec {\rho }$
$$r = \sqrt{z^2 + |\vec{\rho} - \vec{\rho}_0|^2} \simeq z(1 + \frac{|\vec{\rho} - \vec{\rho}_0|^2}{2 z^2}).$$
We will use the above Fresnel approximation to calculate the X-ray ghost image forming correlation by means of propagating the field from one space-time location such as the $m$th sub-source, $(\rho_m, z_m, t_m)$, to another space-time location at $(\rho_j, z_j, t_j)$ [13],
$$g_m(\rho_j, z_j, t_j) = \frac{c_o}{|z_j - z_m|} e^{-i \omega \tau_j} e^{i \frac{\omega}{2 c |z_j - z_m|} |\vec{\rho}_j -\vec{\rho}_m |^2},$$
where $c_0$ is a normalization constant and $\tau _j \equiv t_j - |z_j - z_m|/c$. It may be necessary to emphasize that the Fresnel propagator used here is just an approximation for a point-to-point propagation which is very different from the “beam propagator” commonly used in projection radiography. For the following discussion, we will assume a single frequency source, and will address concerns of temporal coherence in broad spectrum sources in Section 5. Assuming a disk-like source and randomly distributed and randomly radiated point-like sub-sources, we can approximate the sum of $m$ in Eq. (4) into an integral of $\vec {\rho }_s$ on the source plane of $z_s = 0$. Before considering the detectors, we can look at the object plane on beam-two, $z_o$, and the corresponding plane on beam-one, which we will label $z_{gi}$. The result of intensity fluctuation correlation is
$$\begin{aligned} \langle \Delta I(\vec{\rho}_{gi}, z_{gi}) \Delta I(\vec{\rho}_{o}, z_o) \rangle \big{|}_{z_{gi} = z_o} & \propto \Big{\langle} \Big{|} \int d \vec{\rho}_s \,\big{[} g^*_s(\vec{\rho}_{gi}, z_{gi}) \big{]} \big{[} g_s(\vec{\rho}_o, z_o) \big{]} \Big{|}^2 \Big{\rangle} \\ & \propto \textrm{somb}^2 \frac{\pi \Delta \theta_s}{\lambda} |\vec{\rho}_{gi} - \vec{\rho}_{o}|, \end{aligned}$$
where the somb-function is defined as $2J_1(x)/x$, where $J_1(x)$ is the Bessel function of the first kind, and $\Delta \theta _s \approx 2R/z_{gi}$ is the angular diameter of the radiation source. This point-to-spot correlation is similar to that of the original Hanbury Brown-Twiss experiment that started the practice of optical correlation measurements, except it is measured in the Fresnel near-field instead of the Fraunhofer far-field [5052]. This point-to-spot correlation has also been demonstrated at X-ray synchrotron sources [5355]. The width of this correlation is $\sim \lambda / \Delta \theta _s$ so, due to the high energy (short wavelength) nature of X-rays, this point-to-spot correlation is much narrower than that of visible light. For instance, ignoring experimental limitations such as those introduced by the beam splitter used and the pixel size on the photodetector array, a realistic setup with a $\sim 5 \times 10^{-11}$ meter wavelength ($\sim 25$ keV) X-ray source with an angular diameter of 10 mrad would achieve $\sim 5$ nm correlation and the same photon energies with an angular diameter of 100 mrad, perhaps with a specially designed source, could achieve $\sim 0.5$ nm correlation. And this of course can be enhanced further by using higher photon energies (shorter wavelengths) depending on what is appropriate for the object of interest.

 figure: Fig. 2.

Fig. 2. X-ray ghost imaging. A beam splitter (most likely a crystal aligned to utilize Laue diffraction which would not provide 90$^{\circ }$ separation as depicted) creates two paths for the beam, one directed passing through the object and followed by a bucket detector, and the second directed at a 2D photodetection array (CCD or CMOS). The array is placed on the desired ghost image plane, $z_1 = z_{gi}$, corresponding to different “slices”, or cross sections inside the object: $z_{gi} = z_o, z_{gi}' = z_o', z_{gi}'' = z_o'',$ etc.

Download Full Size | PDF

Unless extreme magnification or point-like pixels are used, we can approximate this point-to-spot correlation (somb-function) as a point-to-point correlation (delta-function) between $\vec {\rho }_{gi}$ and $\vec {\rho }_o$ to simplify the math. This point-to-point correlation, or ghost image forming function, identifies a unique internal plane of the object, $z_{gi} = z_o$, with the X-ray ghost image observed from the X-ray intensity fluctuation, or photon number fluctuation, correlation. With the point-to-point approximation and $z_{gi} = z_o$, we may approximate the X-ray intensity fluctuation correlation as

$$\langle \Delta I(\vec{\rho}_{gi}, z_{gi}) \Delta I(\vec{\rho}_{o}, z_o) \rangle \propto \delta (z_{gi} - z_o) \, \delta (|\vec{\rho}_{gi} - \vec{\rho}_{o}|).$$
We now calculate the X-ray lensless ghost image by including a 2D photodetection array on the $z_1 = z_{gi}$ plane and a bucket detector, $D_2$, which integrates all possible X rays transmitted through each transverse coordinate point $\vec {\rho }_o$ of the object. The internal “aperture function” of an object cannot be simply represented as a 2D function $A(\vec {\rho }_{o}, z_o)$ as typically done for visible-light imaging. In visible light imaging, $A(\vec {\rho }_{o})$ is usually used to represent the surface plane of an object; however, for X-ray ghost imaging it is reasonable to model a 3D aperture function representing the internal structure of an object,
$$A(\vec{\rho}_{o}) \approx \int d z_o A(\vec{\rho}_{o}, z_o).$$
Assuming perfect temporal correlation is satisfied experimentally, the intensity fluctuation correlation results in
$$\begin{aligned}\langle \Delta I(\vec{\rho}_{1}, z_1) \Delta I_2 \rangle & \propto \Big{|} \int d \vec{\rho}_{o} \int d z_o \, A(\vec{\rho}_{o}, z_o) \delta (z_1 - z_o)\, \textrm{somb} \frac{\pi \Delta \theta_s}{\lambda} |\vec{\rho}_{1} - \vec{\rho}_{o}| \Big{|}^2 \\ & \approx \Big{|} \int d \vec{\rho}_{o} \int d z_o \, A(\vec{\rho}_{o}, z_o) \delta (z_1 - z_o) \, \delta (|\vec{\rho}_{1} - \vec{\rho}_{o}|) \Big{|}^2 \\ & \approx \big{|} A(\vec{\rho}_{1} = \vec{\rho}_{o}, z_1 = z_o) \big{|}^2, \end{aligned}$$
which represents image formation. In the first line, prior to the point-to-point approximation, we see that the somb-function is the image-forming correlation function that represents the resolving capabilities of the system. As discussed before with the width of the correlation, this can theoretically provide sub-nanometer resolution, but will be limited by physical constraints in the system such as those introduced by the beam splitter used and the pixel size on the photodetector array. It is also worth noting that, even though the selected $z_o = z_{gi}$ plane is “in focus,” there may still be effects present from the out-of-focus planes. These effects may average to zero with enough statistics or may result in a certain level of background “noise” from the blurred images of the out-of-focus planes.

This result also indicates the reproduction of a 2D X-ray ghost image of the internal transverse cross section of $z_o = z_{gi} = z_1$. In other words, the longitudinal position of the 2D X-ray detector array, $z_1$, selected an object plane, $z_o = z_1$, of the internal structure of the object to image. Interestingly, when the 2D X-ray detector array is linearly scanned from $z_1$ to $z'_1$ to $z''_1$ along the optical axis, the selected object plane will be changed from $z_o = z_1$ to $z'_o = z'_1$ to $z''_o = z''_1$, respectively, as illustrated in Fig. 2. The imaging resolution along the $z$-axis is the coherence time of the light, similar to optical coherence tomography [56]. By scanning the position of the 2D X-ray detector array, a set of slices of the internal structure of the object can be grouped together to form a 3D ghost image of the object with sub-nanometer resolution. This differs from traditional X-ray computerized tomography (CT) imaging [57] and ghost tomography (GT) demonstrated by Kingston et al. [34] which rely on rotating the object or revolving the detectors around the object 360$^{\circ }$ [58, 59].

4. X-ray ghost microscope

Unfortunately, the sub-nanometer resolution of lensless ghost imaging is unresolvable by any state-of-the-art 2D photodetector array, such as CCD or CMOS sensors which may have micrometer sized pixels, unless photodetector array technology ever advances to sub-nanometer-sized pixels. In order to resolve the two-photon X-ray ghost image with a standard 2D photodetector array, significant magnification is necessary.

4.1 Configuration I

To magnify the X-ray ghost image to a secondary ghost image plane we may follow the optical setup demonstrated by Valencia et al. in 2005 for visible light ghost imaging. The schematic design of this standard X-ray ghost microscope matches what an optical setup would be and is shown in Fig. 3. For this configuration of the X-ray ghost microscope, focusing X-ray optics such as compound refractive lenses, X-ray zone plates, etc. can be used to produce a magnified secondary ghost image [1416]. Note that focusing X-ray optics may limit the usable photon energies of the X-ray source unless the technology advances to accommodate higher photon energies. As with the primary ghost image, this secondary ghost image is observable from the measurement of the intensity fluctuation correlation. To confirm the nonlocal point-to-point correlation between the object plane, $z_o$, and the secondary image plane of the lens, $z_i$ (where $D_1$ will be placed, $z_1 = z_i$), we adjust Eq. (9) to include the lens system,

$$\begin{aligned}\langle \Delta I(\vec{\rho}_{i}, z_i) \Delta I(\vec{\rho}_{o}, z_o) \rangle & \propto \Big{|} \int d \vec{\rho}_s \big{[} g^*_s(\vec{\rho}_o, z_o) \big{]} \big{[} g_s(\vec{\rho}_i, z_i) \big{]} \Big{|}^2 \\ & \propto \Big{|} \int d \vec{\rho}_s \big{[} g^*_s(\vec{\rho}_o, z_o) \big{]} \big{[} g_s(\vec{\rho}_{gi}, z_{gi}) \int d \vec{\rho}_{gi} \, g_{gi}(\vec{\rho}_L, z_L) \\ & \qquad \times g_{Lens} \int d \vec{\rho}_L \, g_L(\vec{\rho}_i, z_i) \big{]} \Big{|}^2. \end{aligned}$$
Here $g_s(\vec {\rho }_{gi}, z_{gi})$, $g_{gi}(\vec {\rho }_L, z_L)$, $g_{Lens}$, and $g_L(\vec {\rho }_i, z_i)$ are the Green’s functions propagating the field from the source plane to the one-to-one primary ghost image plane, from the primary ghost image plane to the lens plane, from the input plane of the lens to the output plane of the lens, and from the lens plane to the secondary ghost image plane, respectively. This results in,
$$\begin{aligned} \langle \Delta I(\vec{\rho}_{i}, z_i) \Delta I(\vec{\rho}_{o}, z_o) \rangle & \propto \Big{\langle} \Big{|} \int d \vec{\rho}_{gi} \, \delta(|\vec{\rho}_{gi} - \vec{\rho}_o|) \, \textrm{somb} \frac{\pi D }{s_o \lambda} |\vec{\rho}_{gi} - \vec{\rho}_{i} / \mu | \Big{|}^2 \Big{\rangle} \\ & \propto \textrm{somb}^2 \frac{\pi D }{s_o \lambda} |\vec{\rho}_{o} - \vec{\rho}_{i} / \mu |, \end{aligned}$$
indicating a point-to-spot secondary ghost image forming function [22,13]. In Eq. (14), $D$ is the diameter of the lens, $s_o$ is the distance from the primary ghost image to the lens, $s_i$ is the distance from the lens to the secondary ghost image, satisfying the Gaussian thin lens equation $1/s_o + 1/s_i = 1/f$, where $f$ is the focal length of the lens and $\mu = s_i / s_o$ is the magnification factor of the secondary ghost image. The 2D detector array, $D_1$, is now placed on the image plane of the lens (secondary ghost image plane), $z_1 = z_i$. The position of $D_1$ defines the position of $z_{gi}$ and thus defines the “slice” at $z_o = z_{gi}$ of the internal cross section of the object by means of the Gaussian thin-lens equation. We then take into account the complex aperture function $\tilde {A}(\vec {\rho }_o, z_{o} = z_{gi})$ and the bucket detector $D_2$ for the calculation. A magnified secondary ghost image of the aperture function is then observed from the joint detection between the CCD (CMOS), $D_1$, and the bucket detector $D_2$,
$$\begin{aligned}\langle \Delta I(\vec{\rho}_{1}, z_1) \Delta I_2 \rangle & \propto \Big{\langle} \Big{|} \int d \vec{\rho}_{o} \, \tilde{A}(\vec{\rho}_o, z_o = z_{gi}) \, \textrm{somb} \frac{\pi D }{s_o \lambda} |\vec{\rho}_{o} - \vec{\rho}_{1} / \mu| \Big{|}^2 \Big{\rangle} \\ & \propto \int d \vec{\rho}_{o} \big{|} \tilde{A}(\vec{\rho}_o, z_o = z_{gi}) \big{|}^2 \textrm{somb}^2 \frac{\pi D }{s_o \lambda} |\vec{\rho}_{o} - \vec{\rho}_{1} / \mu| \\ & \approx \big{|} \tilde{A}(\vec{\rho}_1 / \mu, z_o = z_{gi}) \big{|}^2. \end{aligned}$$
Scanning $D_1$ from one position to another along the optical axis, or refocusing the microscope from one ghost image plane of $z_{gi}$ to $z'_{gi}$ etc., we obtain a magnified 3D ghost image of the internal structure of the object.

 figure: Fig. 3.

Fig. 3. The X-ray ghost microscope. Focusing X-ray optics are used to produce a magnified secondary ghost image prior to the 2D photodetection array which allows for a high resolution image of the object to be resolvable by a standard CCD or CMOS. Different primary ghost image planes correspond to different “slices” of the object. By scanning the lens and $D_1$ along the optical axis, we select new $z_{gi} = z_o = d$ planes to be in focus, obtaining a magnified 3D image of the internal structure of the object.

Download Full Size | PDF

Due to the high resolution of the primary ghost image, assuming the use of a source with a large angular diameter, the result of the magnified secondary ghost image has identical resolution to an image from a classical X-ray microscope that was obtained with the focusing X-ray optics. Currently, further advancing the fabrication of focusing X-ray optics is a difficult dilemma to overcome (with 10 nm resolution when using a zone plate being state-of-the-art [16]), but as the the fabrication capabilities advance, the potential of this ghost microscope setup will advance as well. In addition to matching the resolution of simply using the focusing X-ray optics to image the object itself, one major benefit to the ghost microscope setup is that the X-ray optics are imaging the ghost image plane, which isn’t physically present, as opposed to imaging the physical object (hence the name “ghost” microscope). Imaging small or thin objects with a classical X-ray microscope may not be an issue, but often it may be desirable to image the deeper internal structure (e.g. a bone) of a thicker object (e.g. a body) for which the focusing X-ray optics can not physically be close enough to the internal structure to image it with proper magnification and resolution. This would not be an issue for the X-ray ghost microscope as the physical object is placed “nonlocally” on the other path following the beam splitter and is simply followed by a bucket detector.

4.2 Configuration II

If the use of X-ray optics is not possible or not preferred, one can place a scintillator on the primary ghost image plane, $z_{gi}$ to convert it into the visible spectrum, allowing it to then be magnified by a visible-light lens system onto a secondary image plane. This is inspired by many current X-ray imaging setups that use scintillator and visible light lenses. Unlike configuration I which has limitations from the focusing X-ray optics, this design is more suitable for higher energy X-ray imaging, such as $\geq 20$ keV, depending on the scintillator used. The schematic design of the adjusted X-ray microscope is shown in Fig. 4. Unlike the standard X-ray ghost microscope which has a clear path for the two-photon amplitudes from the light source to the detectors, here it is not clear that the scintillator-lens system preserves the result of the two-photon interference, and thus the secondary ghost image. To understand how the secondary ghost image is preserved, we can say that the scintillator essentially “detects” the X rays, thus establishing the presence of the X-ray ghost image on the scintillator plane; although, prior to correlation, this plane is simply a distribution of quantum speckles from interfering photon pairs. The scintillator converts these X-ray fluctuations to the visible spectrum, $\Delta I_{\textrm {X}}(\vec {\rho }_{gi}, z_{gi}) \approx \Delta I_{\textrm {v}}(\vec {\rho }_{gi}, z_{gi})$. The lens system then images this distribution of fluctuations and produces a diffraction-limited magnified image of them on the image plane,

$$\begin{aligned}\Delta I(\vec{\rho}_i, z_i) & = \sum_{m \neq n} E^*_m(\vec{\rho}_i, z_i) E_n(\vec{\rho}_i, z_i) \\ & = \sum_{m \neq n} \big{[} \int_{lens} d \vec{\rho}_l \, E^*_m(\vec{\rho}_{gi}) g^*_{gi}(\vec{\rho}_l, z_l) g_{Lens} g^*_l(\vec{\rho}_i, z_i) \big{]} \\ & \qquad \times \big{[} \int_{lens} d \vec{\rho'}_l \, E_n(\vec{\rho}_{gi}) g_{gi}(\vec{\rho'}_l, z'_l) g_{Lens} g_{l'}(\vec{\rho}_i, z_i) \big{]} \\ & = \sum_{m \neq n} E^*_m(\vec{\rho}_{gi}) E_n(\vec{\rho}_{gi}) \, \textrm{somb}^2\frac{\pi D}{s_o \lambda_{\textrm{v}}} |\vec{\rho}_{o} + \vec{\rho}_{i}/\mu| \\ & = \Delta I(\vec{\rho}_{gi}, z_{gi}) \, \textrm{somb}^2\frac{\pi D }{s_o \lambda_{\textrm{v}}} |\vec{\rho}_{o} + \vec{\rho}_{i}/\mu|. \end{aligned}$$
where $\lambda _{\textrm {v}}$ is the center wavelength emitted from the scintillator. It is interesting that the X-ray intensity fluctuations at $(\vec {\rho }_{gi}, z_{gi})$ of the primary ghost image plane are “propagated” to an unique point in the secondary image plane $(\vec {\rho }_i, z_i)$, where $(\vec {\rho }_i, z_i)$ is defined by the somb-function and the Gaussian thin-lens equation of the visible-light microscope. Unlike configuration I, this process allows the use of type (II) ghost imaging as well, which would involve “detection” and imaging of the classical speckle-pattern; although, type (II) ghost imaging does not have the capability of tomography via a linear scan (linear ghost tomography) [58]. Returning to the application of type (I) ghost imaging, a magnified secondary ghost image is observable from the correlation measurement between the intensity fluctuations of the X-ray and the intensity fluctuations of the visible light,
$$\begin{aligned}\langle \Delta I(\vec{\rho}_{1}, z_1) \Delta I_2 \rangle & \propto \Big{\langle} \Big{|} \int d \vec{\rho}_{o} \, \tilde{A}(\vec{\rho}_o, z_o = z_{gi}) \, \textrm{somb} \frac{\pi D }{s_o \lambda_{\textrm{v}}} |\vec{\rho}_{o} - \vec{\rho}_{1} / \mu| \Big{|}^2 \Big{\rangle} \\ & \propto \int d \vec{\rho}_{o} \, |\tilde{A}(\vec{\rho}_o, z_o = z_{gi})|^2 \, \textrm{somb}^2 \frac{\pi D }{s_o \lambda_{\textrm{v}}} |\vec{\rho}_{o} - \vec{\rho}_{1} / \mu|. \end{aligned}$$
Interestingly, because each ghost image plane along this beam represents a different “slice” of the object (as discussed in Section 3.), a thick scintillator may contain multiple ghost image planes depending on the penetrating depth of the X rays. The resolution along the $z$-axis of the image would be more dependent on the depth of field of the lens system as it is focused onto specific planes of the scintillator. It should be noted, however, that this scintillator model assumes an exact point-to-point conversion of the X-rays to visible light. While this may not be actually the case and loss of resolution will be introduced, this has not been taken into account for this proof-of-principle discussion. With a new dependence on visible light, pairing the scintillator with a standard compound microscope will be the best accessible option with a theoretical resolution limit of around 200 nm but, in general, a more practical limit may be micron resolution. However, even with these limitations in resolution, this setup would still have benefits of standard optical microscopy as it allows for imaging the “deeper” internal structure of the object. It allows the scintillator to be placed directly on the desired ghost image plane as opposed to behind the physical object in projectional radiography. If preferred, it would also be possible to combine configurations I and II.

 figure: Fig. 4.

Fig. 4. The adjusted X-ray ghost microscope. Nearly identical to the standard X-ray ghost microscope, but now a scintillator is placed on the ghost image plane to convert the X-ray ghost image into the visible spectrum. Now a lens (or lens system) operating in the visible spectrum produces a magnified secondary ghost image.

Download Full Size | PDF

One factor that may aid in the high resolving capabilities of the X-ray ghost microscope is the optical stability of two-photon interference. It has been shown that ghost imaging and other two-photon interference phenomena can be set up to achieve turbulence-free measurements [3638,51]. This is achieved when the superposed two-photon amplitudes experience the same turbulence and medium vibrations along their optical paths, meaning any composition, density, length, refractive index, or medium vibration induced random phase variations along the optical paths do not have any effect on each individual two-photon interference. While optical turbulence may not be prominent in the X-ray regime, this technique also includes small vibrations of the object or detectors that would become more likely to blur a classical image at the resolutions discussed.

5. Achieving observable two-photon interference with broad-spectrum X-ray sources

The intensity fluctuation correlation utilized for the X-ray ghost microscope is directly related to the second-order correlation function, $G^{(2)}(\boldsymbol{r}_1, t_1; \boldsymbol{r}_2, t_2)$, of X-rays in a thermal state [13,4547]. This second-order correlation function follows a quantum description of light, so we use the quantized notation of photon numbers, $n(\boldsymbol{r}_j, t_j) \propto I(\boldsymbol{r}_j, t_j)$ and $\Delta n(\boldsymbol{r}_j, t_j) \propto \Delta I(\boldsymbol{r}_j, t_j)$ for $j = 1, 2$. Up to this point, we have assumed perfect temporal correlation to allow us to emphasize the spatial portion of this function. Now we will do the opposite and focus on the temporal aspect of the measurement. Similar to ghost imaging with sunlight, we have to face the problem caused by the extremely broad spectrum of X-ray sources and relatively slow detectors. Due to the limited ability of the photodetectors in determining the registration time of a photoelectron, the response time of the photodetectors and the associated electronics will affect the measurement of $G^{(2)}(t_{1}-t_{2})$, where $t_j$ is the registration time of photodetectors $D_1$ and $D_2$. For instance, due to the slow response time of the photodetectors, relatively speaking, the measured second-order correlation may have a much wider temporal width and much smaller amplitude. Note that this also holds true with the relatively slow decay time of the scintillator in configuration II. We may characterize this uncertainty as a response function of the photodetector, $D(\tilde {t}_j - t_j)$, where $t_j$ is the photon annihilation time and $\tilde {t}_j$ is the electronic registration time. Thus, the joint photodetection measurement of $D_1$ and $D_2$ can be treated as a convolution between the response functions and the second-order correlation function $G^{(2)}(t_{1} - t_{2})$,

$$G^{(2)}(\tilde{t}_{1} -\tilde{t}_{2}) = \frac{1}{t^{2}_{c}} \, \int dt_{1} \int dt_{2} \, G^{(2)}(t_{1} - t_{2})\, D(\tilde{t}_{1} - t_{1}) \, D(\tilde{t}_{2} - t_{2}),$$
where we have normalized the function with the response time, or characteristic time, of the photodetector, $t_{c}$. When using fast detectors, the width of the response functions are much narrower then the temporal width of $G^{(2)}(t_{1} -t_{2})$ so the response functions can be treated as delta functions, $D(\tilde {t}-t) \sim t_{c}\delta (\tilde {t}-t)$. In this extreme case, the measured second-order correlation function will reveal the theoretical expectation of $G^{(2)}(t_{1} - t_{2})$. However, when the response times are larger, the situation is different; especially when the temporal widths of the response functions are much wider then that of $G^{(2)}(t_{1} - t_{2})$. In this case, the second-order correlation function itself can then be treated as a delta function and Eq. (18) turns into the following convolution between the response functions of the two photon counting detectors,
$$G^{(2)}(\tilde{t}_{1} -\tilde{t}_{2}) = G^{(2)}(0) \, \Big{(}\frac{\tau_{0}}{t^{2}_{c}}\Big{)} \int d\tau \, D(\tilde{t}_{1} - \tilde{t}_{2}-\tau) \, D(\tau),$$
where $\tau = t_{1} - t_{2}$ and we have normalized the delta function with $G^{(2)}(0) \tau _{0}$. The factor $\tau _{0}$ is the second-order correlation time and is a constant defined as the inverse of the bandwidth of the source, $\tau _{0} = 1/\Delta \nu$. Equation (19) indicates two things: (1) the width of the observed $G^{(2)}(\tilde {t}_{1} -\tilde {t}_{2})$ is now determined by the response function of the photodetectors, which could be significantly broadened compared to the original correlation and (2) the relatively slow response time of the photodetectors may reduce the magnitude of the measured second-order correlation. For Gaussian response functions the convolution yields a Gaussian function of $G^{(2)}(\tilde {t}_{1} -\tilde {t}_{2})$ with a reduced central value of $G^{(2)}(0)$. The reduction factor is roughly $\tau _0 / t_c$.

The above result for slow detectors may not be a problem in the measurement of entangled states; however, it may affect the measurement of a $G^{(2)}$ function for thermal or pseudo-thermal light significantly. The $G^{(2)}$ function of thermal or pseudo-thermal field has two terms, the product of two measured mean photon numbers (trivial part) and the photon number fluctuation correlation (nontrivial part) [13],

$$\begin{aligned} G^{(2)}(\tilde{t}_{1} -\tilde{t}_{2}) & = \frac{1}{t^{2}_{c}} \int dt_{1} \int dt_{2} \, \langle n(t_1) n(t_2) \rangle D(\tilde{t}_{1} - t_{1}) \, D(\tilde{t}_{2} - t_{2})\\ & = \frac{1}{t^{2}_{c}} \int dt_{1} \int dt_{2} \, \big{[} \langle n(t_1) \rangle \langle n(t_2) \rangle + \langle \Delta n(t_1) \Delta n(t_2) \rangle \big{]} D(\tilde{t}_{1} - t_{1}) \, D(\tilde{t}_{2} - t_{2}). \end{aligned}$$
The time average over $t_c$ has no effect on the first term (product of mean photon numbers) for a CW thermal or pseudo-thermal field with broad spectrum, because $\langle n(t_1) \rangle \langle n(t_2) \rangle = \bar {n}_1 \bar {n}_2$ is a constant. However, it reduces the magnitude of the second term (photon number fluctuation correlation) significantly when $\tau _{0} \ll t_{c}$,
$$\begin{aligned} & \frac{1}{t^{2}_{c}} \int dt_{1} \int dt_{2} \, \langle \Delta n(t_1) \Delta n(t_2) \rangle D(\tilde{t}_{1} - t_{1}) D(\tilde{t}_{2} - t_{2}) \\ & \qquad \qquad \approx \frac{\bar{n}_1 \bar{n}_2}{t^{2}_{c}} \int dt_{1} \int dt_{2} \, \tau_0 \delta (t_1 - t_2) D(\tilde{t}_{1} - t_{1}) D(\tilde{t}_{2} - t_{2}) \\ & \qquad \qquad \approx \frac{\bar{n}_1 \bar{n}_2}{t^{2}_{c}} \tau_0 \int d\tau \, D(\tilde{t}_{1} - \tilde{t}_{2}-\tau) \, D(\tau). \end{aligned}$$
For Gaussian response functions, the magnitude of the photon number fluctuation term is thus roughly $\tau _0 / t_c$ times that of the product of mean photon numbers, which may reach $10^{-6}$ (one part of a million) for an X-ray source with $10^{15}$Hz bandwidth ($\tau _0 \sim 10^{-15}$) and a nanosecond photodetector ($t_c \sim 10^{-9}$).

How does one distinguish the relatively small, photon number fluctuations correlation part of the second-order correlation of a thermal field with short coherence time $\tau _0$? The obvious approach is to use fast photodetectors with $t_c \sim \tau _0$. However, even if faster, individual photodetectors are used instead of a full detector array, the state-of-the-art technology has not been able to produce photodetectors fast enough to achieve this for broadband light sources yet. One realistic approach that could utilize current photodetector technology is using a pulsed thermal or pseudo-thermal radiation source with pulse width of $\tau _p$ and recording a single pulse per frame. In this case, the time integral in Eq. (20) shall have similar effects on the first term of $G^{(2)}$ (product of mean photon numbers) as it had in Eq. (21),

$$\begin{aligned} & \frac{1}{t^{2}_{c}} \int dt_{1} \int dt_{2} \langle n(t_1) \rangle \langle n(t_2) \rangle D(\tilde{t}_{1} - t_{1}) D(\tilde{t}_{2} - t_{2}) \\ & \qquad \qquad \approx \frac{\bar{n }_{1} \bar{n}_2}{t^{2}_{c}} \int dt_{1} \int dt_{2} \, \tau_p \delta (t_1 - t_2) D(\tilde{t}_{1} - t_{1}) D(\tilde{t}_{2} - t_{2}) \\ & \qquad \qquad \approx \frac{\bar{n }_{1} \bar{n}_2}{t^{2}_{c}} \tau_p \int d\tau \, D(\tilde{t}_{1} - \tilde{t}_{2}-\tau) \, D(\tau). \end{aligned}$$
where we have approximated the short pulsed $\langle n(t_1) \rangle \langle n(t_2) \rangle$ as a delta-function such that $\langle n(t_1) \rangle \langle n(t_2) \rangle \sim \bar {n}_1 \bar {n}_2 \tau _p \delta (t_1 - t_2)$. For Gaussian response functions, this results in the pulse being broadened and reduced by a factor of $\tau _p / t_c$. In the case where the pulse width is approximately the same as the coherence time, $\tau _p \approx \tau _0$, all of the photons in the pulse are considered second-order coherent. This results in,
$$\frac{1}{t^{2}_{c}} \int dt_{1} \int dt_{2} \langle n(t_1) \rangle \langle n(t_2) \rangle D(\tilde{t}_{1} - t_{1}) D(\tilde{t}_{2} - t_{2}) \approx \frac{\bar{n}_1 \bar{n}_2}{t^{2}_{c}} \tau_0 \int d\tau \, D(\tilde{t}_{1} - \tilde{t}_{2}-\tau) \, D(\tau),$$
which matches Eq. (21), meaning the two terms in the second-order correlation are reduced by the same factor. This allows the photon number fluctuation correlation (intensity fluctuation correlation) to become more distinguishable and is suitable for measuring the second-order correlation of X rays.

Based on the above discussions, we may conclude that a short pulsed X-ray source with large angular diameter is a better solution for X-ray ghost imaging and the X-ray ghost microscope. Although synchrotrons have been chosen for many X-ray imaging techniques (projectional radiography, phase-contrast radiography, and CXDI) along with various other X-ray experiments in which point-like sources are desired, synchrotron X-ray sources are not necessarily the ideal candidate for the X-ray ghost microscope setups because higher-resolution is expected from a larger angular diameter X-ray source. An X-ray source with a large diameter and short pulsed emission is currently in development at Los Alamos National Laboratory and other institutions.

6. Summary

In summary, we have analyzed the working mechanism of the X-ray ghost microscope. By applying Einstein’s granularity picture of radiation, which includes X rays, we found that two-photon interference is able to produce an image-forming correlation that forms a lensless X-ray ghost image with high spatial resolution, but no magnification. Utilizing a high-energy ($>$ 20 keV) source with a relatively large angular diameter, a ghost image with sub-nanometer resolution could be achieved but is not obtainable with current photodetector array technology. Either with the addition of focusing X-ray optics or the help of a scintillator and an optical lens system, the primary ghost image can be mapped onto a secondary image plane with significant magnification and reduced (but still high) resolution. Because the primary ghost image is along a different beam than the physical object, the focusing X-ray optics or scintillator can be placed anywhere along the beam, allowing for high-resolution, magnified imaging of the deep internal structure of an object without direct physical constraints. Linearly scanning the secondary imaging setup along the propagation axis, a 3D image can be obtained with Linear Ghost Tomography. We have also found that a short-pulsed X-ray source and measuring a single pulse per frame is preferred for observing the two-photon interference induced intensity fluctuation correlation of X rays due to the wide bandwidth of the source. Even with measuring a single pulse per frame, in order to achieve a higher degree of second-order correlation it is preferred to have longer coherence time, $\tau _0$, that is more comparable to the pulse width, $\tau _p$. In addition to this, the desire for a large angular diameter makes synchrotron X-ray sources less favorable as the resolution of the primary ghost image would be greatly reduced. As sensor arrays and X-ray sources see technological advancements, this X-ray ghost microscope will open up new capabilities that would be of interest to the fields of physics, materials science, and biological imaging.

Funding

Los Alamos National Laboratory.

Acknowledgments

The authors thank J. N. Sprigg and B. Joshi for helpful discussions.

Disclosures

The authors declare no conflicts of interest related to this article.

References

1. A. Sakdinawat and D. Attwood, “Nanoscale x-ray imaging,” Nat. Photonics 4(12), 840–848 (2010). [CrossRef]  

2. M. Holt, R. Harder, R. Winarski, and V. Rose, “Nanoscale hard x-Ray microscopy methods for materials studies,” Annu. Rev. Mater. Res. 43(1), 183–211 (2013). [CrossRef]  

3. V. Cnudde and M. N. Boone, “High-resolution x-ray computed tomography in geosciences: A review of the current technology and applications,” Earth-Sci. Rev. 123, 1–17 (2013). [CrossRef]  

4. P. J. Withers, “X-ray nanotomagraphy,” Mater. Today 10(12), 26–34 (2007). [CrossRef]  

5. V. V. Lider, “X-ray microscopy,” Phys.-Usp. 60(2), 187–203 (2017). [CrossRef]  

6. H. N. Chapman and K. A. Nugent, “Coherent lensless x-ray imaging,” Nat. Photonics 4(12), 833–839 (2010). [CrossRef]  

7. E. A. Spiller, Soft X-ray Optics (SPIE Optical Engineering, 1994).

8. M. A. Flower, Webb’s Physics of Medical Imaging (CRC, 2012).

9. A. Snigirev, I. Snigireva, V. Kohn, S. Kuznetsov, and I. Schelokov, “On the possibilities of x-ray phase contrast microimaging by coherent high-energy synchrotron radiation,” Rev. Sci. Instrum. 66(12), 5486–5492 (1995). [CrossRef]  

10. B. Yu, L. Weber, A. Pacureanu, M. Langer, C. Olivier, P. Cloetens, and F. Peyrin, “Evaluation of phase retrieval approaches in magnified x-ray phase nano computerized tomography applied to bone tissue,” Opt. Express 26(9), 11110–11124 (2018). [CrossRef]  

11. M. Born and E. Wolf, Principles of Optics (Pergamon, 1980).

12. E. Hecht, Optics (Addison-Wesley, 2002).

13. Y. H. Shih, An Introduction to Quantum Optics: Photon and Biphoton Physics (CRC, 2011).

14. C. G. Schroer, M. Kuhlmann, U. T. Hunger, T. F. Günzler, O. Kurapova, S. Feste, F. Frehse, B. Lengeler, M. Drakopoulos, A. Somogyi, A. S. Simionovici, A. Snigirev, I. Snigireva, C. Schug, and W. H. Schröder, “Nanofocusing parabolic refractive x-ray lenses,” Appl. Phys. Lett. 82(9), 1485–1487 (2003). [CrossRef]  

15. H. Mimura, S. Handa, T. Kimura, H. Yumoto, D. Yamakawa, H. Yokoyama, S. Matsuyama, K. Inagaki, K. Yamamura, Y. Sano, K. Tamasaku, Y. Nishino, M. Yabashi, T. Ishikawa, and K. Yamauchi, “Breaking the 10 nm barrier in hard-x-ray focusing,” Nat. Phys. 6(2), 122–125 (2010). [CrossRef]  

16. W. Chao, P. Fischer, T. Tyliszczak, S. Rekawa, E. Anderson, and P. Naulleau, “Real space soft x-ray imaging at 10 nm spatial resolution,” Opt. Express 20(9), 9777–9783 (2012). [CrossRef]  

17. F. Pfeiffer, “X-ray ptychography,” Nat. Photonics 12(1), 9–17 (2018). [CrossRef]  

18. I. Robinson and R. Harder, “Coherent x-ray diffraction imaging of strain at the nanoscale,” Nat. Mater. 8(4), 291–298 (2009). [CrossRef]  

19. J. Cheng and S. Han, “Incoherent coincidence imaging and its applicability in x-ray diffraction,” Phys. Rev. Lett. 92(9), 093903 (2004). [CrossRef]  

20. M. Stampanoni, R. Mokso, F. Marone, J. Vila-Comamala, S. Gorelick, P. Trtik, K. Jefimovs, and C. David, “Phase-contrast tomography at the nanoscale using hard x rays,” Phys. Rev. B 81(14), 140105 (2010). [CrossRef]  

21. A. Momose, “Demonstration of phase-contrast x-ray computed tomography using an x-ray interferometer,” Nucl. Instrum. Methods Phys. Res., Sect. A 352(3), 622–628 (1995). [CrossRef]  

22. A. Valencia, G. Scarcelli, M. D’Angelo, and Y. H. Shih, “Two-photon imaging with thermal light,” Phys. Rev. Lett. 94(6), 063601 (2005). [CrossRef]  

23. G. Scarcelli, V. Berardi, and Y. H. Shih, “Can two-photon correlation of chaotic light be considered as correlation of intensity fluctuations?” Phys. Rev. Lett. 96(6), 063602 (2006). [CrossRef]  

24. R. Meyers, K. S. Deacon, and Y. H. Shih, “Ghost-imaging experiment by measuring reflected photons,” Phys. Rev. A 77(4), 041801 (2008). [CrossRef]  

25. T. A. Smith, Y. H. Shih, Z. Wang, X. Li, B. Adams, M. Demarteau, R. Wagner, J. Xie, L. Xie, R. Y. Zhu, L. Zhang, and C. Hu, “From optical to x-ray ghost imaging,” Nucl. Instrum. Methods Phys. Res., Sect. A 935, 173–177 (2019). [CrossRef]  

26. T. B. Pittman, Y. H. Shih, D. V. Strekalov, and A. V. Sergienko, “Optical imaging by means of two-photon quantum entanglement,” Phys. Rev. A 52(5), R3429–R3432 (1995). [CrossRef]  

27. A. Schori, D. Borodin, K. Tamasaku, and S. Shwartz, “Ghost imaging with paired x-ray photons,” Phys. Rev. A 97(6), 063804 (2018). [CrossRef]  

28. H. Yu, R. Lu, S. Han, H. Xie, G. Du, T. Xiao, and D. Zhu, “Fourier-transform ghost imaging with hard x rays,” Phys. Rev. Lett. 117(11), 113901 (2016). [CrossRef]  

29. D. Pelliccia, A. Rack, M. Scheel, V. Cantelli, and D. M. Paganin, “Experimental x-ray ghost imaging,” Phys. Rev. Lett. 117(11), 113902 (2016). [CrossRef]  

30. A. Schori and S. Shwartz, “X-ray ghost imaging with a laboratory source,” Opt. Express 25(13), 14822–14828 (2017). [CrossRef]  

31. A. X. Zhang, Y. H. He, L. A. Wu, L. M. Chen, and B. B. Wang, “Tabletop x-ray ghost imaging with ultra-low radiation,” Optica 5(4), 374–377 (2018). [CrossRef]  

32. D. Pelliccia, M. P. Olbinado, A. Rack, A. M. Kingston, G. R. Myers, and D. M. Paganin, “Towards a practical implementation of x-ray ghost imaging with synchrotron light,” IUCrJ 5(4), 428–438 (2018). [CrossRef]  

33. D. Ceddia and D. M. Paganin, “Random-matrix bases, ghost imaging, and x-ray phase contrast computational ghost imaging,” Phys. Rev. A 97(6), 062119 (2018). [CrossRef]  

34. A. M. Kingston, D. Pelliccia, A. Rack, M. P. Olbinado, Y. Cheng, G. R. Myers, and D. M. Paganin, “Ghost tomography,” Optica 5(12), 1516–1520 (2018). [CrossRef]  

35. Y. Y. Kim, L. Gelisio, G. Mercurio, S. Dziarzhytski, M. Beye, L. Bocklage, A. Classen, C. David, O. Y. Gorobtsov, R. Khubbutdinov, S. Lazarev, N. Mukharamova, Y. N. Obukhov, B. Röner, K. Schlage, I. A. Zaluzhnyy, G. Brenner, R. Röhlsberger, J. von Zanthier, W. Wurth, and I. A. Vartanyants, “Ghost imaging at an XUV free-electron laser,” Phys. Rev. A 101(1), 013820 (2020). [CrossRef]  

36. R. E. Meyers, K. S. Deacon, and Y. H. Shih, “Turbulence-free ghost imaging,” Appl. Phys. Lett. 98(11), 111115 (2011). [CrossRef]  

37. T. A. Smith and Y. H. Shih, “Turbulence-free double-slit interferometer,” Phys. Rev. Lett. 120(6), 063606 (2018). [CrossRef]  

38. T. A. Smith and Y. H. Shih, “Turbulence-free two-photon double-slit interference with coherent and incoherent light,” Opt. Express 27(23), 33282–33297 (2019). [CrossRef]  

39. R. S. Bennink, S. J. Bentley, R. W. Boyd, and J. C. Howell, “Quantum and classical coincidence imaging,” Phys. Rev. Lett. 92(3), 033601 (2004). [CrossRef]  

40. A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, “Correlation imaging, quantum and classical,” Phys. Rev. A 70(1), 013802 (2004). [CrossRef]  

41. A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, “Ghost imaging with thermal light: Comparing entanglement and classical correlation,” Phys. Rev. Lett. 93(9), 093602 (2004). [CrossRef]  

42. J. H. Shapiro, “Computational ghost imaging,” Phys. Rev. A 78(6), 061802 (2008). [CrossRef]  

43. Additionally, X-ray free-electron lasers (XFELs) do not produce X rays in a thermal state, so type (I) ghost imaging is not possible unless the coherent X-ray beam is thermalized into a pseudo-thermal state; although, type (II) ghost imaging is possible when a beam diffuser is used to introduce speckles.

44. A. Einstein, “On a heuristic viewpoint concerning the production and transformation of light,” Ann. Phys. 322(6), 132–148 (1905). [CrossRef]  

45. R. J. Glauber, “Photon correlations,” Phys. Rev. 10(3), 84–86 (1963). [CrossRef]  

46. R. J. Glauber, “The quantum theory of optical coherence,” Phys. Rev. 130(6), 2529–2539 (1963). [CrossRef]  

47. M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University, 1997).

48. This is trivial when using visible light but becomes more difficult with high-energy X-rays. So far, Laue diffraction of crystals has been used with success but advancement in kinoform X-ray beam splitters may prove useful in the future.

49. M. Lebugle, G. Seniutinas, F. Marschall, V. A. Guzenko, D. Grolimund, and C. David, “Tunable kinoform x-ray beam splitter,” Opt. Lett. 42(21), 4327–4330 (2017). [CrossRef]  

50. R. Hanbury Brown and R. Q. Twiss, “Correlation between photons in two coherent beams of light,” Nature 177(4497), 27–29 (1956). [CrossRef]  

51. R. Hanbury Brown and R. Q. Twiss, “A test of a new type of stellar interferometer on sirius,” Nature 178(4541), 1046–1048 (1956). [CrossRef]  

52. R. Hanbury Brown and R. Q. Twiss, “The question of correlation between photons in coherent light rays,” Nature 178(4548), 1447–1448 (1956). [CrossRef]  

53. Y. Kunimune, Y. Yoda, K. Izumi, M. Yabashi, X. W. Zhang, T. Harami, M. Ando, and S. Kikuta, “Two-photon correlations in x-rays from a synchrotron radiation source,” J. Synchrotron Radiat. 4(4), 199–203 (1997). [CrossRef]  

54. E. Gluskin, E. E. Alp, I. McNulty, W. Sturhahn, and J. Sutter, “A classical Hanbury Brown-Twiss experiment with hard x-rays,” J. Synchrotron Radiat. 6(5), 1065–1066 (1999). [CrossRef]  

55. M. Yabashi, K. Tamasaku, and T. Ishikawa, “Characterization of the transverse coherence of hard synchrotron radiation by intensity interferometry,” Phys. Rev. Lett. 87(14), 140801 (2001). [CrossRef]  

56. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254(5035), 1178–1181 (1991). [CrossRef]  

57. G. N. Hounsfield, “Computerized transverse axial scanning (tomography): Part 1. Description of system,” Br. J. Radiol. 46(552), 1016–1022 (1973). [CrossRef]  

58. A benefit of “linear” ghost tomography as opposed to “rotational” ghost tomography is the theoretical capability of having a 3D array of photodetectors such that each 2D plane captures a slice of the object simultaneously, removing the need to move the object or photodetectors. Of course the technology of 3D photodetector arrays does not exist yet, but initial assessments have recently been completed for multilayer, ultrafast silicon sensor development.

59. X. Li, P. Chu, Z. Wang, C. M. O’Shaughnessy, C. Morris, M. Demarteau, R. Wagner, J. Xie, L. Xia, R. Y. Zhu, L. Zhang, C. Hu, B. Adams, J. Katsoudas, Y. Ding, C. Segre, T. A. Smith, and Y. H. Shih, “Nucl. Instrum. Methods Phys. Res., Sect. A,” Nucl. Instrum. Methods Phys. Res., Sect. A 942, 162414 (2019). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1.
Fig. 1. Experimentally achieved with a beam splitter, this “unfolded” schematic of the 2005 thermal light ghost imaging demonstration of Valencia et al. helps depict the symmetry of the two-photon amplitudes. The one-to-one lensless ghost image and the magnified secondary ghost image are results of two-photon interference. Of the many two-photon amplitudes, such as the red and blue two-photon amplitudes of the random $m$th and $n$th pair of photons shown here, each superpose constructively at a corresponding $\vec {\rho }_{i}$ to from an image.
Fig. 2.
Fig. 2. X-ray ghost imaging. A beam splitter (most likely a crystal aligned to utilize Laue diffraction which would not provide 90$^{\circ }$ separation as depicted) creates two paths for the beam, one directed passing through the object and followed by a bucket detector, and the second directed at a 2D photodetection array (CCD or CMOS). The array is placed on the desired ghost image plane, $z_1 = z_{gi}$, corresponding to different “slices”, or cross sections inside the object: $z_{gi} = z_o, z_{gi}' = z_o', z_{gi}'' = z_o'',$ etc.
Fig. 3.
Fig. 3. The X-ray ghost microscope. Focusing X-ray optics are used to produce a magnified secondary ghost image prior to the 2D photodetection array which allows for a high resolution image of the object to be resolvable by a standard CCD or CMOS. Different primary ghost image planes correspond to different “slices” of the object. By scanning the lens and $D_1$ along the optical axis, we select new $z_{gi} = z_o = d$ planes to be in focus, obtaining a magnified 3D image of the internal structure of the object.
Fig. 4.
Fig. 4. The adjusted X-ray ghost microscope. Nearly identical to the standard X-ray ghost microscope, but now a scintillator is placed on the ghost image plane to convert the X-ray ghost image into the visible spectrum. Now a lens (or lens system) operating in the visible spectrum produces a magnified secondary ghost image.

Equations (23)

Equations on this page are rendered with MathJax. Learn more.

E ( r , t ) = m E m ( r , t ) = m E m ( r m , t m ) g m ( r m , t m ; r , t ) ,
I ( r , t ) = E ( r , t ) E ( r , t ) = m E m ( r , t ) n E n ( r , t ) = m | E m ( r , t ) | 2 + m n E m ( r , t ) E n ( r , t ) = m | E m ( r , t ) | 2 .
Δ I ( r , t ) = m n E m ( r , t ) E n ( r , t ) .
Δ I ( r 1 , t 1 ) Δ I ( r 2 , t 2 ) = m n E m ( r 1 , t 1 ) E n ( r 1 , t 1 ) p q E p ( r 2 , t 2 ) E q ( r 2 , t 2 ) = m n E m ( r 1 , t 1 ) E n ( r 1 , t 1 ) E n ( r 2 , t 2 ) E m ( r 2 , t 2 ) .
G ( 2 ) ( r 1 , t 1 ; r 2 , t 2 ) = m n | E m ( r 1 , t 1 ) E n ( r 2 , t 2 ) + E n ( r 1 , t 1 ) E m ( r 2 , t 2 ) | 2 = I ( r 1 , t 1 ) I ( r 2 , t 2 ) + Δ I ( r 1 , t 1 ) Δ I ( r 2 , t 2 ) ,
g 0 ( ρ , z , t ) = 1 r e i ( ω t k r ) ,
r = z 2 + | ρ ρ 0 | 2 z ( 1 + | ρ ρ 0 | 2 2 z 2 ) .
g m ( ρ j , z j , t j ) = c o | z j z m | e i ω τ j e i ω 2 c | z j z m | | ρ j ρ m | 2 ,
Δ I ( ρ g i , z g i ) Δ I ( ρ o , z o ) | z g i = z o | d ρ s [ g s ( ρ g i , z g i ) ] [ g s ( ρ o , z o ) ] | 2 somb 2 π Δ θ s λ | ρ g i ρ o | ,
Δ I ( ρ g i , z g i ) Δ I ( ρ o , z o ) δ ( z g i z o ) δ ( | ρ g i ρ o | ) .
A ( ρ o ) d z o A ( ρ o , z o ) .
Δ I ( ρ 1 , z 1 ) Δ I 2 | d ρ o d z o A ( ρ o , z o ) δ ( z 1 z o ) somb π Δ θ s λ | ρ 1 ρ o | | 2 | d ρ o d z o A ( ρ o , z o ) δ ( z 1 z o ) δ ( | ρ 1 ρ o | ) | 2 | A ( ρ 1 = ρ o , z 1 = z o ) | 2 ,
Δ I ( ρ i , z i ) Δ I ( ρ o , z o ) | d ρ s [ g s ( ρ o , z o ) ] [ g s ( ρ i , z i ) ] | 2 | d ρ s [ g s ( ρ o , z o ) ] [ g s ( ρ g i , z g i ) d ρ g i g g i ( ρ L , z L ) × g L e n s d ρ L g L ( ρ i , z i ) ] | 2 .
Δ I ( ρ i , z i ) Δ I ( ρ o , z o ) | d ρ g i δ ( | ρ g i ρ o | ) somb π D s o λ | ρ g i ρ i / μ | | 2 somb 2 π D s o λ | ρ o ρ i / μ | ,
Δ I ( ρ 1 , z 1 ) Δ I 2 | d ρ o A ~ ( ρ o , z o = z g i ) somb π D s o λ | ρ o ρ 1 / μ | | 2 d ρ o | A ~ ( ρ o , z o = z g i ) | 2 somb 2 π D s o λ | ρ o ρ 1 / μ | | A ~ ( ρ 1 / μ , z o = z g i ) | 2 .
Δ I ( ρ i , z i ) = m n E m ( ρ i , z i ) E n ( ρ i , z i ) = m n [ l e n s d ρ l E m ( ρ g i ) g g i ( ρ l , z l ) g L e n s g l ( ρ i , z i ) ] × [ l e n s d ρ l E n ( ρ g i ) g g i ( ρ l , z l ) g L e n s g l ( ρ i , z i ) ] = m n E m ( ρ g i ) E n ( ρ g i ) somb 2 π D s o λ v | ρ o + ρ i / μ | = Δ I ( ρ g i , z g i ) somb 2 π D s o λ v | ρ o + ρ i / μ | .
Δ I ( ρ 1 , z 1 ) Δ I 2 | d ρ o A ~ ( ρ o , z o = z g i ) somb π D s o λ v | ρ o ρ 1 / μ | | 2 d ρ o | A ~ ( ρ o , z o = z g i ) | 2 somb 2 π D s o λ v | ρ o ρ 1 / μ | .
G ( 2 ) ( t ~ 1 t ~ 2 ) = 1 t c 2 d t 1 d t 2 G ( 2 ) ( t 1 t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) ,
G ( 2 ) ( t ~ 1 t ~ 2 ) = G ( 2 ) ( 0 ) ( τ 0 t c 2 ) d τ D ( t ~ 1 t ~ 2 τ ) D ( τ ) ,
G ( 2 ) ( t ~ 1 t ~ 2 ) = 1 t c 2 d t 1 d t 2 n ( t 1 ) n ( t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) = 1 t c 2 d t 1 d t 2 [ n ( t 1 ) n ( t 2 ) + Δ n ( t 1 ) Δ n ( t 2 ) ] D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) .
1 t c 2 d t 1 d t 2 Δ n ( t 1 ) Δ n ( t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) n ¯ 1 n ¯ 2 t c 2 d t 1 d t 2 τ 0 δ ( t 1 t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) n ¯ 1 n ¯ 2 t c 2 τ 0 d τ D ( t ~ 1 t ~ 2 τ ) D ( τ ) .
1 t c 2 d t 1 d t 2 n ( t 1 ) n ( t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) n ¯ 1 n ¯ 2 t c 2 d t 1 d t 2 τ p δ ( t 1 t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) n ¯ 1 n ¯ 2 t c 2 τ p d τ D ( t ~ 1 t ~ 2 τ ) D ( τ ) .
1 t c 2 d t 1 d t 2 n ( t 1 ) n ( t 2 ) D ( t ~ 1 t 1 ) D ( t ~ 2 t 2 ) n ¯ 1 n ¯ 2 t c 2 τ 0 d τ D ( t ~ 1 t ~ 2 τ ) D ( τ ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.