0.54 μm resolution two-photon interference with dispersion cancellation for quantum optical coherence tomography.

Quantum information technologies harness the intrinsic nature of quantum theory to beat the limitations of the classical methods for information processing and communication. Recently, the application of quantum features to metrology has attracted much attention. Quantum optical coherence tomography (QOCT), which utilizes two-photon interference between entangled photon pairs, is a promising approach to overcome the problem with optical coherence tomography (OCT): As the resolution of OCT becomes higher, degradation of the resolution due to dispersion within the medium becomes more critical. Here we report on the realization of 0.54 μm resolution two-photon interference, which surpasses the current record resolution 0.75 μm of low-coherence interference for OCT. In addition, the resolution for QOCT showed almost no change against the dispersion of a 1 mm thickness of water inserted in the optical path, whereas the resolution for OCT dramatically degrades. For this experiment, a highly-efficient chirped quasi-phase-matched lithium tantalate device was developed using a novel 'nano-electrode-poling' technique. The results presented here represent a breakthrough for the realization of quantum protocols, including QOCT, quantum clock synchronization, and more. Our work will open up possibilities for medical and biological applications.

One of the most distinct feature of quantum physics is quantum entanglement. Entanglement attracted attention first in the test of nonlocality of quantum mechanics123, and started to be considered as an essential resource for quantum information protocols4, including quantum key distribution56, quantum teleportation78, and quantum computation91011. Recently, the application of quantum entanglement for metrology and sensing is attracting attention1213. One recent example is an entanglement-enhanced microscope, where a photon-number entangled state is used as probe light to enhance the sensitivity14. Here we focus on another example: an application of the frequency entangled state of photons for optical coherence tomography (OCT)1516.

OCT based on low-coherence interference (LCI)17 has been widely used in various fields, including medical applications such as imaging of the retina and cardiovascular system1516. Figure 1a shows a schematic diagram of OCT. Broadband light from a source is divided at a beam splitter (BS). One beam is incident on a sample after passing through a dispersive medium, while the other beam is reflected from a mirror with a temporal delay τ. The OCT interference fringe I(τ) is obtained by the detection of the interfered light intensity with varying delay τ (Fig. 1a inset). When the bandwidth of the source is made broader to achieve higher resolution, the resolution, far from being improved, degrades due to dispersion in the medium. This constitutes a severe problem in OCT18. Although dispersion effect can be compensated by inserting a ‘phantom’, a medium with the same dispersion, in the reference path19, it requires a priori knowledge of both the structure and the frequency-dependent refractive index of the target object. Furthermore, as the resolution becomes higher, a slight difference between the target object and the ‘phantom’ becomes a crucial problem.

Figure 1

OCT and QOCT schemes.

(a) Schematic diagram of OCT. is the interfered light intensity measured at a detector with varying delay τ (inset). BS is a beam splitter. (b) Schematic diagram of QOCT. is the coincidence count rate counted at two single photon detectors with varying delay τ (inset).

As an alternative method, quantum optical coherence tomography (QOCT)2021222324 based on the two-photon interference (TPI) of frequency-entangled photon pairs25 has been proposed20 (Fig. 1b). Broadband entangled photon pairs can be generated via a spontaneous parametric down-conversion process from a nonlinear crystal. Signal photons reflected at a sample through a dispersive medium and idler photons with a temporal delay τ, interfere quantum mechanically at a BS, and coincidence detection events are counted at two single photon detectors. The QOCT interference dip , which is so called Hong-Ou-Mandel (HOM) dip25, is obtained by the coincidence count rate with varying delay τ (Fig. 1b inset). Due to the frequency correlation of entangled photon pairs, the resolution of QOCT (the width of the HOM dip) does not change even with group velocity dispersion (GVD) in the medium2627. This ‘dispersion cancellation’ of TPI was first demonstrated with 19 μm resolution21 and very recently with 3 μm resolution, where the GVD effect becomes significant28.

In this work, we report 0.54 μm resolution TPI with dispersion cancellation for ultra-high resolution QOCT. For highly efficient generation of ultra-broadband (166 THz, λ = 660–1040 nm) frequency-entangled photon pairs, we developed a 1st-order chirped quasi-phase-matched (QPM)29 lithium tantalate30 device using the nanofabrication technology for fine electrode patterns. The device was pumped using a narrowband (~100 kHz) pump laser to ensure the dispersion cancellation of TPI28. We also constructed stable interferometer setup with hybrid ultra-broadband detection systems (HUBDeS) operated at room temperature. In addition, dispersion cancellation in TPI was demonstrated against a 1 mm thickness of water inserted in the optical path. Almost no degradation in resolution was observed for TPI (from 0.54 μm to 0.56 μm), which is in contrast to the significant degradation in LCI resolution (from 1.5 μm to 7.8 μm). The 0.54 μm resolution in air, which corresponds to the resolution of 0.40 μm in water, is the best among the previously achieved LCI resolution 0.75 μm31 for OCT and also TPI resolutions for QOCT, including that (0.85 μm23) where dispersion cancellation was not verified.

Results

The 1st-order QPM device was developed using a ‘nano-electrode-poling’ technique. Note that the conversion efficiency of the 1st-order QPM32 is supposed to be almost one order magnitude (9 times) larger than that of the 3rd-order QPM2223; however, it is impossible to realize using conventional photolithography techniques. For high-power ultraviolet pumping, Mg-doped stoichiometric lithium tantalate (Mg:SLT)3033 was selected because it has a short absorption edge around 270 nm34 and high thermal conductivity3536. Furthermore, Mg:SLT is free from visible-light-induced nonlinear infrared absorption, which is observed in Mg:LiNbO337. The new device fabrication process that was employed is shown in Fig. 2a. A 500 μm thick Mg:SLT wafer doped with 1.0 mol% Mg was used to suppress photorefractive damage and a 100 nm thick Al film was evaporated onto the wafer surface. A nanoscale resist pattern was defined using electron beam lithography (EBL) and then transferred to an Al pattern by dry etching on the surface of the Mg:SLT wafer. EBL free from the diffraction limit of light provides over 10 times higher accuracy than photolithography, but special care is required to avoid charge-up in the ferroelectric Mg:SLT substrate. The width of the Al electrode was 400 nm and the duty ratio of the period was approximately 10% to the period assuming reasonable sidewise expansion (scanning electron microscopy (SEM) images are shown in Fig. 2b). An electric field of 2 kV/mm was then applied in the vacuum chamber to achieve a high electric field contrast by suppression of the surface screening charge for precise control of the domain duty ratio. Domain sidewise motion was successfully controlled over 6,000 domains along the 20 mm device length with a 6.7% chirped poling period varying from 3.12 to 3.34 μm (optical microscopy images are shown in Fig. 2c). Although each EBL scanning area is limited to a width of 200 μm  mm device length), the extremely accurate translation stage enables a 20 nm connection error, which results in a low phase slip between scanned areas. The fluctuation in the domain duty ratio (typically 0.65) was less than 10% to the period (Fig. 2c), which suggests uniform conversion efficiency along the entire spectral range.

Figure 2

Fabrication of chirped QPM lithium tantalate device using nano-electrode-poling technique.

(a) Device fabrication process. Black arrows indicate the direction of ferroelectric spontaneous polarization in the periodically poled device. Mg:SLT is Mg-doped stoichiometric lithium tantalate. (b) Scanning electron microscopy images (upper and lower) of the 400 nm wide Al electrodes fabricated for the poling period of 3.2 μm. (c) Optical Microscopy images of periodically poled structures. The QPM period varies from 3.12 μm (left) to 3.34 μm (right) along the 20 mm device length.

The experimental setup for the TPI is shown in Fig. 3a (Details are given in Methods). The chirped QPM Mg:SLT device is set in a temperature controlled metal holder (Fig. 3a upper inset). The temperature T is controlled by a thermoelectric cooler with an accuracy of 0.1 K. We have found that in ultra-high resolution regime the linewidth of a pump laser degrades the dispersion tolerance28. To ensure the dispersion tolerance, we used a narrowband pump laser (~100 kHz) with a wavelength of 401 nm. The pump beam is focused to a diameter of 80 μm with a lens; this diameter is sufficiently smaller than the cross-section of the device (500 μm (vertical) × 800 μm (horizontal)), and then cut by filters after the device. Collinearly emitted photons had ultra-broad bandwidth with a device temperature T = 351 K as shown in Fig. 3b. The spectrum spans from 660 nm with a sharp rise. We think that the offset in the shorter wavelength region (≤650 nm) is stray fluorescence from the device. Considering the frequency correlation with the center wavelength of 802 nm, the spectrum should span up to 1040 nm with a bandwidth of 166 THz. The low intensity in the longer wavelength region (≥950 nm) may be due to the low coupling efficiency to the optical fiber for the spectrometer we used. Figure 3c shows non-collinearly emitted photon pairs (signal photons and idler photons) with an emission angle of 1.5 degrees relative to the pump beam with a device temperature T = 353 K. With the detection event of a signal photon as a trigger, the coincidence count rate of the idler photons was typically 5% (5 × 105 Hz) that of the single photon count rate (~1 × 107 Hz). Signal and idler photons coupled to fiber couplers (FCs) are transferred to the TPI interferometer through polarization-maintaining fibers. The delay τ is determined by the physical position of the FC. The coincidence count rate is obtained by two HUBDeSs, which consist of a Si avalanche photodiode (APD) with a detection bandwidth spanning from 400 to 1060 nm and an InGaAs APD with a bandwidth spanning from 950 to 1150 nm. For the LCI experiments, we used signal photons as a low-coherence light source with the exactly same bandwidth of the source for TPI and the HUBDeS 1 for detection (Fig. 3a lower inset). For the check of the dispersion effect, a 1 mm thickness of water enclosed by two thin glass cover plates is inserted in the optical paths of both TPI and LCI interferometers. Note that according to the conventional definition for OCT16, the interferograms below are expressed in units of the ‘delay’ cτ/2, where c is the speed of light, considering the physical displacement of the delay mirrors in OCT and QOCT systems (Fig. 1).

Figure 3

Experimental setup.

(a) TPI and LCI interferometers with hybrid ultra-broadband detection systems (HUBDeS). The upper inset shows the chirped QPM device set in the temperature controlled metal holder. HUBDeS consists of a Si avalanche photodiode (APD) and an InGaAs APD. The lower inset shows the LCI interferometer. A 1 mm thickness of water can be inserted in the optical path in TPI and LCI interferometers. BS, beam splitter; FC, fiber coupler; PMF, polarization-maintaining fiber; Mg:SLT, Mg-doped stoichiometric lithium tantalate. (b) Frequency spectrum of collinearly emitted photons from the device. Experimental data (red dots) and the theoretical curve (black line) are plotted. (c) Frequency spectra of photon pairs generated from the device in non-collinear emission. The observed data are plotted for signal photons (red dots) and idler photons (blue dots). The transmission efficiency and detection efficiency of the spectrometer were calibrated (b,c).

The experimental results are shown in Fig. 4. First, the LCI signal was obtained as plotted by the red dots in Fig. 4a. The full width at half maximum (FWHM) of the interference fringe is 1.5 μm, which is slightly larger than the FWHM of 1.1 μm for theoretical calculation (blue line in Fig. 4a) assuming a rectangular spectral shape with a 166 THz bandwidth. We think that this degradation is induced by the GVD: When the thicknesses of the glass of BS in the optical paths of the probe beam and the reference beam (Fig. 3a lower inset) are different, the degradation occurs by the GVD of the glass. We have found that the difference of the FWHMs can be explained when we assume that the difference of the thicknesses of the glass is just 150 μm.

Figure 4

Obtained LCI and TPI signals.

(a,b) LCI fringes obtained in the LCI interferometer using signal photons as the source. The fringes without (a) and with (b) a 1 mm thickness of water inserted in the optical path. The experimental data (red dots) and the theoretical curves (blue line) are plotted in units of the delay. The integration time was 1 second per point. (c,d) TPI dips obtained in the TPI interferometer. The dips without (c) and with (d) the 1 mm thickness of water. The experimental data (red dots), the theoretical curves (blue lines) and the Gaussian fitting curves (black dashed lines) are plotted in units of the delay. The integration time was 10 seconds per point. The red lines connecting the data points are a guide to the eye (a–d).

Figure 4b shows the LCI signal when a 1 mm thickness of water is inserted in the optical path. Due to the GVD of the 1 mm thickness of water38, the interference fringe becomes much broader. The FWHM of the fringe is 7.8 μm, which is more than 5 times larger than that without the 1 mm thickness of water sample. These results illustrate how GVD effect becomes crucial in ultra-high resolution regime.

The main results of this paper are shown in Fig. 4c,d. Figure 4c shows the TPI signal without the water sample. The FWHM of the observed HOM dip is 0.54 ± 0.05 μm, which surpasses 0.75 μm of LCI31 for OCT using a broadband laser light with a center wavelength of 725 nm, which is the current best record resolution as far as we know. The resolution of 0.54 μm in air corresponds to the resolution of 0.40 μm in water or biological tissue. The theoretical curve (blue line in Fig. 4c) assuming the observed spectrum of the signal photons (Fig. 3c) fits well with the experimental data. Note that for the same bandwidth of the photon source with a rectangular shaped spectrum, the FWHM of the TPI is a half of the LCI28, which is also an advantage of QOCT especially when the bandwidth of the optical window is limited.

Finally, experimental demonstration of the dispersion cancellation of an ultra-high resolution TPI for QOCT is shown in Fig. 4d. The HOM dip is almost unchanged from Fig. 4c even when a 1 mm thickness of water is inserted in the optical path, which is in striking contrast to the case of LCI (Fig. 4a,b). The FWHM of the HOM dip is 0.56 ± 0.04 μm and the difference between the FWHMs with and without the 1 mm thickness of water is 0.02 μm, which is smaller than the margins of errors. The theoretical curve (blue line in Fig. 4d) calculated taking the higher-order dispersion of the 1 mm thickness of water38 into account fits also well with the experimental data. Note that the asymmetricity of the HOM dip observed in the experimental data is due to the higher odd-order dispersion2839, which can be also observed in the theoretical curve. Our theoretical calculation suggests that the resolution will be still kept 0.60 μm even when a thickness of water is increased to 3 mm. For thicker medium, ‘phantom’ can also be used to compensate the dispersion effect for QOCT, similar to OCT. In case of QOCT, the thickness of the phantom does not have to be exactly the same with the medium due to the inherit dispersion cancellation by frequency entanglement demonstrated here. In case of OCT, on the other hand, the thickness or the GVD of the phantom has to be exactly the same with the medium as we discussed on Fig. 4a, which is technically very difficult. The detail of the theoretical estimation of FWHMs and the calculation are given in Methods.

Conclusion

In conclusion, we have achieved 0.54 μm resolution TPI with dispersion cancellation for ultra-high resolution QOCT, which surpasses the current record resolution 0.75 μm of LCI for OCT. We developed a 1st-order chirped QPM Mg:SLT device using a ‘nano-electrode-poling’ technique with EBL for ultra-broadband entangled photon pair generation (166 THz, λ = 660–1040 nm). A 6.7% chirped poling period varying from 3.12 to 3.34 μm was fabricated in over 6,000 domains along a 20 mm long device. We constructed stable interferometers with HUBDeSs that consist of commercially available single photon detectors operated at room temperature. In addition, dispersion cancellation was demonstrated in TPI against a 1 mm thickness of water inserted in the optical path. Almost no degradation in resolution was observed in TPI (from 0.54 μm to 0.56 μm), whereas the LCI resolution was significantly degraded (from 1.5 μm to 7.8 μm). These results will open the door to ultra-high resolution QOCT imaging with depth resolution less than half a micrometer. Such ultra-high resolution QOCT will be beneficial for many different areas, for example, following up the change of the retinal thickness with the ultra-high resolution will greatly help the early detection of glaucoma40. For this end, the increase in the flux of entangled photons is very important. The flux of 0.3 μW has already been realized using bulk QPM devices414243. Further increase in flux up to tens of μW can be expected using slab or ridge type waveguide structures4445. A promising direction may be an OCT/QOCT hybrid system where OCT is used for a quick wide-range scan and QOCT is used when ultra-high resolution/high precision observation is required. We also note that the ultra-high resolution TPI with dispersion cancellation demonstrated here is useful for quantum protocols, including not only QOCT but also quantum clock synchronization46, time-frequency entanglement measurement47 and multimode frequency entanglement48.

Methods

Details of experimental setup

The pump laser system consists of a single-frequency continuous wave Ti:sapphire laser (MBR-110, Coherent) excited by a diode-pumped solid state laser (Verdi G-10, Coherent) and a resonant frequency doubling unit (MBD-200, Coherent). The output of the Ti:sapphire laser (wavelength: 802 nm; linewidth: approximately 100 kHz) is frequency doubled by second-harmonic generation. It is used as the pump beam with a power of 100 mW. This narrow linewidth (~100 kHz) ensures the dispersion tolerance of TPI in ultra-high resolution regime28. We treated the narrowband pump beam as a monochromatic pump for numerical calculations based on the theory2839. The focused pump beam at the chirped quasi-phase-matched device has a confocal parameter (i.e., twice the Rayleigh length) of 18 mm.

For the measurements of frequency spectra, generated entangled photons from the device are coupled to the the polarization-maintaining fiber and sent to a 300-mm spectrometer with a 150-grooves/mm grating blazed at 80 nm (SP-2358, Princeton Instruments) and a charge coupled device (CCD) camera (Pixis:100BRX, Princeton Instruments). The transmission efficiency of the spectrometer and the quantum efficiency of the CCD camera were calibrated in Fig. 3b,c.

Each of HUBDeS shown in Fig. 3a consists of a Si APD (SPCM-AQRH-14, Perkin Elmer) and an InGaAs APD (id401, idQuantique). The coincidence count rate C between two HUBDeSs is obtained by the sum of coincidence counts between two APDs as , where for example is the coincidence count rate between Si APD 1 in HUBDeS 1 and InGaAs APD 2 in HUBDeS 2. The coincidence time window of the coincidence counter (id800, idQuantique) was 2 ns.

Analysis of interference signals

For LCI experimental data (plotted by red dots in Fig. 4a,b), a FWHM of interference fringe is determined by the full width at the middle between a base line and a peak height. The base line is the average of the whole interference fringe along the delay. The peak height is calculated using the smallest or largest data point. In theoretical calculations for LCI (blue lines in Fig. 4a,b), we calculated based on the theory28 assuming a rectangular spectral shape with a center wavelength of 802 nm and a bandwidth of 166 THz for the light source. In the case with a 1 mm thickness of water inserted in the optical path, we assumed the 2nd-order dispersion and the 3rd-order dispersion of the water38. For TPI experimental data (red dots in Fig. 4c,d), a FWHM of interference dip is determined by the FWHM of a Gaussian fitting curve (black dashed lines in Fig. 4c,d). The Gaussian fit uses the whole range of the dip along the delay and assumes the experimental visibility. The experimental visibilities of dips with and without the water dispersion are 0.67 ± 0.05 and 0.73 ± 0.04, respectively. In theoretical calculations for TPI (blue lines in Fig. 4c,d), we assumed the observed spectrum of the signal photons (Fig. 3c), experimental coincidence count rate and experimental visibilities. We also assumed the GVD and the 3rd-order dispersion of the water38 for the case with the 1 mm thickness of water inserted.

Additional Information

How to cite this article: Okano, M. et al. 0.54 μm resolution two-photon interference with dispersion cancellation for quantum optical coherence tomography. Sci. Rep. 5, 18042; doi: 10.1038/srep18042 (2015).