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

Multi-channel swept source optical coherence tomography concept based on photonic integrated circuits

Open Access Open Access

Abstract

In this paper, we present a novel concept for a multi-channel swept source optical coherence tomography (OCT) system based on photonic integrated circuits (PICs). At the core of this concept is a low-loss polarization dependent path routing approach allowing for lower excess loss compared to previously shown PIC-based OCT systems, facilitating a parallelization of measurement units. As a proof of concept for the low-loss path routing, a silicon nitride PIC-based single-channel swept source OCT system operating at 840 nm was implemented and used to acquire in-vivo tomograms of a human retina. The fabrication of the PIC was done via CMOS-compatible plasma-enhanced chemical vapor deposition to allow future monolithic co-integration with photodiodes and read-out electronics. A performance analysis using the results of the implemented photonic building blocks shows a potential tenfold increase of the acquisition speed for a multi-channel system compared to an ideal lossless single-channel system with the same signal-to-noise ratio.

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

1. Introduction

Optical coherence tomography (OCT) is an imaging technique able to visualize the layer composition in a sample by interfering low-amplitude light reflected from the layer boundaries with a reference light [1,2]. The non-invasive nature of OCT makes it well suited for medical applications as well as for industrial inspection. Wider utilization of OCT is, however, limited by the cost and the size of state-of-the-art systems. A reduction of both would enable point-of-care applications with light-weight portable or even handheld systems. This can be achieved by using photonic integrated circuits (PICs) to substitute bulk and/or fiber optics systems which are used for conventional OCT implementations. PICs consist of a combination of photonic building blocks substituting individual fiber-optical and bulk components of an optical setup. In addition, optoelectronic and electronic components can be monolithically co-integrated with PICs allowing for detection and electronic processing of the optical signal on a single chip. Therefore, complete optoelectronic systems can be realized on a single chip within an area of <1 cm$^{2}$. Multiple integrated circuits can be fabricated in parallel on a single wafer, which results in a significant cost reduction for each individual system.

State-of-the-art OCT systems employ the Fourier domain (FD) modality, which offers faster acquisition speeds and higher mechanical stability compared to the time domain modality, typically using a moving reference mirror [3]. In FD-OCT the layer information is encoded in the spectral response of the interference signal between sample and reference light. The layer information is extracted via Fourier transformation. The spectral response is obtained in two different ways in so-called spectral domain (SD) and swept source (SS) OCT. While in SD-OCT a broadband light source and a spectrometer is used, in SS-OCT a high-speed wavelength tunable laser source is paired with a synchronized photodiode pair. With respect to PIC-based OCT, the SS-OCT approach holds more promise of high imaging performance because it involves less demanding integrated photonic key building blocks. In particular, SD-OCT requires an integrated spectrometer, e.g. in the form of an arrayed waveguide grating (AWG), which is difficult to implement with sufficiently high quality to compete with state-of-the-art bulk optical systems. In comparison, in SS-OCT the interference signal is directly measured with photodiodes. Several groups have already reported PIC-based SS-OCT systems. Table 1 summarizes these results and highlights the most important key parameters and implemented integrated components.

Tables Icon

Table 1. List of previously reported PIC-based SS-OCT systems.

Most of the previously reported PIC-based SS-OCT systems attempted an one-to-one replication of conventional OCT concepts using photonic building blocks. This resulted in a reduced imaging performance compared to state-of-the-art OCT systems because disadvantages of PICs were not compensated by inherent strengths. Although the potential reductions in cost and size are highly beneficial, the practical relevance of PIC-based OCT systems will strongly depend on their imaging performance. The main disadvantages of PIC-based OCT systems are increased excess losses and wavelength dispersion compared to bulk/fiber components. Coupling light into and out of the PIC and on-chip waveguide propagation loss (PL) represent major loss contributions, which affect the achievable signal-to-noise ratio (SNR). As a result of the higher wavelength dispersion, several implementations of the PIC-based SS-OCT systems exhibited a reduction of axial resolution [4,5,7,8,11].

The inherent strengths of PICs besides cost and size advantages are the capability of combining multiple functions on a monolithic chip and the resulting high mechanical stability of the associated photonic building blocks. For bulk/fiber optics systems multiple functions require additional components and extra efforts for alignment and maintenance. Using photolithography for the PIC fabrication allows a high integration density of photonic building blocks without additional cost, fabrication time, or assembly steps. These advantages become increasingly prevalent with more functionality implemented on the PIC. The high integration density has been exploited to some degree in the implementation of an OCT receiver [9] allowing balanced detection, polarization sensitive and/or polarization diversity measurements. Another approach to capitalize on these strengths for OCT is the parallelization of measurements employing multiple light spots emitted from the PIC probing a sample. With multiple measurement spots the time required for a line scan (B-scan) or volumetric scan is significantly reduced. Multiple probing spots have been previously implemented with a cascade of $1\times 2$ splitters in a SS-OCT system [12]. In this system, each measurement spot has a different on-chip path length resulting in a time division multiplexed interference signal. The interference signal is detected with a single photodiode and the individual imaging positions are separated in post processing. By using only a single photodiode the read-out and processing speed is inherently limited. Moreover, no balanced detection can be used, which causes a high unbalanced relative intensity noise characteristic of SS-OCT [2]. Another downside of this implementation is the fact that the cascade of $1\times 2$ splitters is also used to recombine the signal in the optical return path, which induces significant excess loss. In conclusion, this approach does not take real advantage of parallelization, which would require a low-loss routing of the sample light towards an interferometer. A dedicated optical return path would further allow the implementation of a Mach-Zehnder configuration and therefore balanced detection, which is a prerequisite for a high-performance SS-OCT system. Complete measurement channels can be formed by co-integrating a photodiode pair for each probing spot, therefore removing the limitation in the read-out speed of a single external photodiode.

In this paper, we propose a PIC-based multi-channel SS-OCT concept that avoids high routing losses inherent to other PIC-based systems demonstrated so far. At the core of this concept is a low-loss polarization dependent path routing approach, which we experimentally apply in a single-channel configuration for OCT imaging of a healthy human retina. On the basis of these results we provide an in-depth performance analysis of a multi-channel implementation.

2. Multi-channel SS-OCT concept

 figure: Fig. 1.

Fig. 1. Conceptual drawing of a PIC-based four-channel OCT system. The input light is coupled to the PIC with a fiber via an inverted taper. The input polarization is set to transverse electric (TE)-like polarization by minimizing the light measured at the photodiode (PD) at the termination port of the first polarization beam splitter (PBS). The optical power is split into the reference and sample path via a directional coupler (DC) with an asymmetric splitting ratio and into the individual channels with a cascade of Y-branches. On both the sample and the reference path the outgoing light passes a PBS, which redirects the returning light with the orthogonal polarization towards the multimode interferometer (MMI). The polarization rotation from the forward propagating TE-like mode to the backward propagating transverse magnetic (TM)-like mode is achieved outside of the PIC with a $\lambda /4$-plate rotated by 45$^{\circ }$. On its way from the $\lambda /4$-plate to the sample and reference mirror the light is circularly polarized. The returning reflected light is converted into the orthogonal polarization relative to the outgoing light by passing the $\lambda /4$-plate a second time. The returning reference light crosses the sample path twice, which can be achieved efficiently with low-loss and low-crosstalk dual polarization crossings (DPCs). The highly divergent beams from the PIC are collimated with a microlens array (MLA). The PIC only provides short parts of the sample and reference paths avoiding high propagation losses and wavelength dispersion effects.

Download Full Size | PDF

Figure 1 shows the proposed PIC-based multi-channel OCT concept in a four-channel configuration. In this system, an orthogonal polarization orientation is introduced between light propagating from the source to the sample and light returning from the sample. Therefore, an almost lossless polarization based on-chip path routing is achieved with a polarization beam splitter (PBS). The rotation of the polarization is introduced outside the PIC via a double pass (forward and backward) through a $\lambda /4$-plate which is rotated by 45$^{\circ }$ relative to the polarization axis of the waveguide. In combination with optoelectronic and electronic components monolithically co-integrated on the PIC a multi-channel SS-OCT system can be realized that holds the promise of meeting the demands for compactness, cost-reduction, and at same time high imaging performance.The presented concept avoids high routing losses of light by efficiently separating the forward optical path (between the light source and the sample) from the backward optical path (between the sample and the photodiodes). All PIC-based OCT concepts demonstrated so far use an optical power splitter to separate the forward optical path from the backward optical path, which inherently induces significant excess loss associated with the fixed splitting ratio. In fiber optical systems, low-loss optical circulators can be used to separate the forward and backward optical path. However, integrated optical circulators require complex fabrication steps, and even then they exhibit a limited operational wavelength bandwidth in current implementations [15,16]. Broadband PIC-based polarization beam splitters can be implemented without intricate fabrication steps in various waveguide materials [17]. Therefore, using a PBS in a polarization dependent path routing allows easy fabrication while avoiding significant loss of source light and of returning light from the sample, which is crucial for a multi-channel system. Another advantage of the proposed routing concept stems from the polarization rotation outside of the PIC. Light reflected before the rotation, e.g. at the end facet, keeps its polarization state and therefore is suppressed on its way towards the interferometer. This reduces reflection artifacts in the OCT measurement, which are caused by similar path lengths of light reflected at the end facets of the sample output and the reference output, as observed in the measurement by Schneider et al. [11].

3. Retinal OCT imaging at 840 nm

The proposed PIC-based multi-channel SS-OCT concept facilitates high imaging speed, while keeping losses to a minimum, and, thus, provides advantages in many fields of application of OCT. In particular, the acquisition speed is of high importance for retinal imaging because of the rapid uncontrollable movement of the eye. Furthermore, OCT angiography is gaining increasing attention for the imaging of microvascular structures in the retina [2]. For OCT angiography, multiple tomograms have to be recorded at the same position with a high acquisition rate to reconstruct the blood flow in microvascular structures. Hence, retinal imaging is a good benchmark to evaluate the proposed PIC-based multi-channel OCT concept. With respect to the PIC part of this concept, the polarization based low-loss optical path routing is the most crucial aspect to enable the multi-channel configuration. Hence, in this work we restrict our experiments to the demonstration of the path routing concept for a single-channel SS-OCT system.

Since retinal OCT imaging is usually performed in the wavelength region around 840 nm, we use silicon nitride (SiN) waveguides with silicon dioxide (SiO$_2$) as cladding material. Wire waveguides were chosen because of their high optical confinement and versatile application range. A SiN layer thickness of 160 nm was chosen as a compromise between several aspects. One of the main aspects was the coupling efficiency from fiber/bulk optics to the PIC and vice versa. Due to the use of two polarization states at the sample port the coupling has to be optimized for both transverse electric (TE)-like and transverse magnetic (TM)-like polarized light. This requires the use of a quadratic cross section at the end facet of the PIC resulting in similar field distributions for both polarizations. According to simulation and confirmed via measurements, a cross section of 160 nm×160 nm provides a numerical aperture (NA) of about 0.13, which is compatible with a commercially available high NA fiber (Corning RC PM 850) and is still manageable with a microlens. While a smaller cross section would further expand the mode field diameter and allow the use of conventional single mode fibers, a waveguide width of 160 nm is already at the limit of the lithographic resolution. Besides, for a thinner waveguide core the coupling losses to the silicon substrate would become larger because the lower cladding thickness is limited to a given maximum value due to boundary conditions given by the fabrication process. A waveguide width of 700 nm was selected for routing to ensure single mode propagation for both polarizations. With a cross section of 700 nm×160 nm the minimum bend radius was selected as 150 µm resulting in bend losses below 0.05 dB per 180$^{\circ }$ circular bend for both polarizations [18]. This radius allows for a multi-channel configuration with a reasonable pitch of 500 µm between two adjacent sample ports. Figures 2(a)-(c) show the routing waveguide cross section and the corresponding mode fields.

 figure: Fig. 2.

Fig. 2. (a) Waveguide cross section used for routing on the PIC. Normalized electric field amplitude distribution of light with (b) TE-like and (c) TM-like polarization. (d) Corresponding measured propagation losses in the wavelength range of 790 nm to 890 nm.

Download Full Size | PDF

The SiN waveguides were fabricated at ams AG in a CMOS-compatible process using low-temperature plasma-enhanced chemical vapor deposition (PECVD), deep UV photolithography, and reactive ion etching, which allows monolithic co-integration in the back-end-of-line (BEOL) with silicon photodiodes and read-out electronics [19]. The upper and lower SiO$_2$ cladding thickness was set to 3 µm. The chip end facets were deep etched to enable efficient coupling of light from fibers or bulk optics to the waveguides. The refractive indices were determined with ellipsometry to be 1.916 for SiN and 1.455 for SiO$_2$ at 840 nm. The fabricated PICs exhibit propagation losses below 1 dB/cm and 0.5 dB/cm for TE-like and TM-like polarization, respectively, in the relevant wavelength range, as can be seen in Fig. 2(d).

There are four essential photonic building blocks required for the proof-of-concept single-channel PIC-based SS-OCT system: a mode size converter to couple light efficiently to and from the PIC, a power splitter to split the optical power into reference and sample path, a polarization beam splitter to redirect the light in the optical return path, and finally an interferometer. To avoid degradation of the axial resolution in an OCT measurement all of these photonic building blocks have to exhibit broadband wavelength behaviour. The design bandwidth was set to 100 nm. For the mode size conversion an inverted taper was designed, which is depicted in Fig. 3(a). In this photonic building block the waveguide width is reduced to lower the confinement and expand the mode into the cladding material. At the minimum width of the taper the resulting mode field diameter and the corresponding numerical aperture provide a better match with fiber and bulk optic components resulting in lower coupling losses. Due to the underlying operation principle inverted tapers have a low wavelength dependence and, thus, are well suited for broadband wavelength applications. To ensure a smooth mode transition at the beginning and the end of the taper a sinusoidal envelope function (covering a half period) is used for the width. Figure 4(a) shows the measurement results for the inverted taper using a polarization maintaining high NA fiber at the input in combination with index matching fluid. The coupling loss amounted to ~1.2 dB.The design (see Fig. 3(c)) and performance of the broadband 90/10 directional coupler for TE-like polarization at 840 nm, acting as power splitter, has been reported in a previous publication [20]. The broadband behaviour is achieved by introducing a phase shift section between two directional couplers. The broadband polarization beam splitter uses the same phase shifting concept as the broadband directional coupler with an additional point symmetry [21], as shown in Fig. 3(d). Figure 4(c) shows the measurement results in terms of polarization extinction ratio defined as the ratio of optical power in the through port to the power in the suppressed port. Photonic building blocks based on directional couplers exhibit comparatively low insertion losses since they do not have abrupt changes in the cross section, which makes an experimental evaluation of the insertion loss difficult. Hence, the insertion losses (see Fig. 4(d)) have only been determined with 3D finite difference time domain (FDTD) using a commercial tool from Lumerical. The final photonic building block required to build a single-channel SS-OCT system is the interferometer. For this purpose, a $2\times 2$ multimode interferometer (MMI) was designed (see Fig. 3(b)). MMIs are especially well suited for interference because they can be designed with a low power imbalance across a wide wavelength range, which is of utmost importance for balanced detection. Furthermore, they feature a high fabrication tolerance because of their relatively large dimension compared to the width variations caused by the lithography process. Figure 4(b) shows the measurement results indicating insertion losses (IL) below 0.7 dB and a power imbalance (IB) below 0.5 dB in the targeted wavelength range. These MMI results were experimentally determined with a symmetric Mach-Zehnder interferometer (MZI) consisting of two cascaded MMIs in a back-to-back configuration [22]. The insertion loss is defined by

$$\textrm{IL} [\textrm{dB}] = \frac{-10\,\textrm{log}_{10}\,(P_1+P_2)}{2},$$
where $P_{1,2}$ are the normalized optical powers measured at the MZI output ports. The factor 2 results from the cascaded MMI configuration. The imbalance, which is defined as the ratio between the output powers of the two MMI output ports, was derived from the measured extinction ratio (ER) using the following relationship
$$\begin{aligned}\textrm{ER} [\textrm{dB}] & = 10\,\textrm{log}_{10}\,(P_{1}/P_{2}) \\ &= 10\,\textrm{log}_{10}\,\bigg(\frac{4}{\rho+\rho^{-1}-2\textrm{cos}(2\theta)}\bigg), \end{aligned}$$
where $\rho = 10\;\hat{}\;(\textrm {IB}/10)$ is the imbalance in linear scale and $\theta$ is the phase error. Assuming $\theta =0$ the resulting maximum power imbalance is
$$\textrm{IB}_{\textrm{max}} [\textrm{dB}] = 10\,\textrm{log}_{10}\,\bigg(1+\frac{2}{\eta}\Big(1+\sqrt{1+\eta}\Big)\bigg),$$
where $\eta = 10\;\hat{}\;(\textrm {ER}/10)$ is the extinction ratio in linear scale.

 figure: Fig. 3.

Fig. 3. Schematics of the broadband photonic building blocks required for low-loss single-channel SS-OCT system: (a) inverted taper for efficient bulk/fiber optic component coupling to and from the PIC; (b) multimode interferometer (MMI) optimized for TM-like polarization to interfere incoming light from the sample and reference paths; (c) directional coupler to split the optical input path into sample and reference path with 90% of the optical power on the through port and 10% on the cross port; (d) polarization beam splitter for redirecting the returning light in the orthogonal polarization state towards the interferometer with TE-like polarization suppressed towards the cross port and vice versa for the TM-like polarization.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Measurement (crosses) and simulation (lines) results of photonic building blocks: (a) insertion loss for the inverted taper with sinusoidal width envelope function (averaged from three different samples); (b) insertion loss and imbalance of the MMI, both determined with a symmetric Mach-Zehnder interferometer comprising two MMIs; (c) polarization extinction ratio (PER), and (d) losses of the PBS for both polarizations. The simulations for the MMIs were made using the eigenmode expansion (EME) method in FIMMPROP (Photon Design), whereas the simulations for the PBS were performed with the 3D FDTD method.

Download Full Size | PDF

Figure 5 depicts the schematic of the PIC that is used in a single-channel SS-OCT system comprising the presented photonic building blocks. The lengths of the reference and sample paths on the PIC are identical for each polarization to minimize dispersion mismatch.

 figure: Fig. 5.

Fig. 5. Schematic of PIC for single-channel SS-OCT system: The input laser light is coupled to the TE-like mode on the PIC from the west side via the inverted taper. The grating coupler located at the suppression port of the PBS provides optical feedback for adjusting an external polarization controller to optimize the input polarization. The TE-like polarized light is then separated by the directional coupler with 90% of the light going into the sample path (north) and 10% going into the reference path (south). Both in the sample and the reference path a PBS redirects the returning polarization rotated light towards the MMI. The interfered light exits the output ports on the east side, which are spatially separated to avoid cross talk between the output ports. The length of the sample and the reference paths are matched to avoid wavelength dispersion mismatch on the PIC.

Download Full Size | PDF

Preliminary OCT measurements with the fabricated PIC were conducted at the Medical University of Vienna. A 840 nm swept laser source with a bandwidth of 60 nm was provided by EXALOS. Figure 6 shows the tomogram of a healthy human retina measured with the fabricated PIC. This preliminary measurement demonstrates that in-vivo tomograms of the retina layer structure can be achieved with a PIC-based SS-OCT system using the proposed routing approach. Some trade-offs were made to allow the use of available fiber/bulk components, which reduced the overall system performance independent of the PIC. Instead of the $\lambda /4$-plates in the sample and reference paths fiber-based paddle polarization controllers were used to manipulate the polarization. This results in a reduced broadband capability of the system compared to using an achromatic $\lambda /4$-plate. Furthermore, the polarization controllers use standard single mode fibers and therefore the high NA fiber had to be replaced, which resulted in a $\approx$ 1.5 dB penalty compared to the best achievable value for the coupling loss. This was the case at the input, sample, and reference ports. The interfered light exiting the PIC end facet was projected onto the photodiodes of a commercially available dual balanced detector in free-space using multiple lenses and mirrors. This induced high excess loss, which resulted in lower contrast in the tomogram. An optimized geometric configuration of the dual balanced photodiodes matching the distance of the output ports is expected to significantly improve the efficiency. In this system a sensitivity of 79 dB was measured with a sample beam power of 0.6 mW and an acquisition rate per sample position (A-scan) of 20 kHz. The axial resolution was determined as 5.5 µm in air from the FWHM of the signal peak generated with a mirror as sample. This value is close to the theoretical limit of 5.2 µm for a central wavelength of 840 nm and a bandwidth of 60 nm. This result indicates that the PIC provides good broadband performance in the investigated wavelength region.

 figure: Fig. 6.

Fig. 6. Tomogram (averaged over 400 consecutive B-scans) of a healthy human retina acquired with the PIC-based single-channel SS-OCT system. The layer structure of the retina can be clearly seen. The axial resolution appears to be lower than the determined value of 5.5 µm because of the high number of averages resulting in a blurring of the image due to insufficient motion correction of the micro movements of the eye ball.

Download Full Size | PDF

4. Performance analysis

With the knowledge of the optical characteristics of the individual building blocks the potentially achievable performance of a PIC-based multi-channel SS-OCT system can be estimated. Although parallelization brings an obvious advantage in the imaging speed, there will be no net gain in the performance of a multi-channel system if the same imaging quality as that of a single-channel system can only be obtained by averaging over multiple scans with the number of scans being larger than the number of parallel channels. As benchmark for this analysis we consider an ideal lossless single-channel system with an A-scan rate of 100 kHz. In our analysis, we assume that the sweep rate of the swept source can be reduced arbitrarily to compensate losses introduced by the photonic building blocks for achieving the same SNR as the lossless system. We define the effective A-scan rate as the adjusted source sweep rate multiplied by the number of channels and use it as key performance indicator. In our concept the returning light from the eye experiences the same losses for each OCT channel independent of the number of channels. Consequently, the SNR is not reduced by increasing the number of channels. Therefore, the multi-channel effective A-scan rate is simply the A-scan rate of the lossy single-channel system multiplied by the number of channels. The only limitation stems from the available source power, which ultimately defines the maximum number of OCT channels. By increasing the number of OCT channels further losses are introduced requiring more source power to reach the desired single spot optical power on the sample. To determine the required optical source power the losses in the forward optical path towards the eye have to be calculated.

We consider a multi-channel SS-OCT PIC with $2^{N}$ output ports, where $N$ is the number of cascading steps. As shown in Fig. 1, in a multi-channel system there are two additional photonic building blocks compared to the single-channel configuration discussed in the previous section – broadband Y-branches for equal power splitting to the individual channels and dual polarization crossings. We have already implemented both building blocks for the 840 nm wavelength region [23,24]. Table 2 provides a summary of the loss values of all building blocks in the multi-channel system for a wavelength of 840 nm. Bend losses are not listed explicitly because radii can be increased at multiple positions making the accumulated propagation loss the main contribution. Although the losses are dependent on the wavelength, this primarily affects the achievable axial resolution rather than the SNR of the OCT system.

Tables Icon

Table 2. Summary of the propagation loss (PL), insertion loss ($\Gamma$), and crosstalk (CT) values for the individual building blocks to be used in the proposed multi-channel SS-OCT system. All values are given in dB for the central wavelength of 840 nm. They have been experimentally determined with fabricated devices except the ones indicated with $^{S}$, which have been derived from 3D FDTD simulations.

In reference to Fig. 1, for the loss of source power towards the eye two types of losses can be distinguished – losses that scale with the number of cascading steps and losses that do not. The number of Y-branches and crossings along the optical path of one channel equals the number of cascading steps $N$. For the Y-branches both the insertion losses and the equal splitting (3 dB) have to be included. Concerning the propagation length scaling with the number of cascading steps, the additional length in x-direction is given by the spacing of the channels. The channel spacing is set to 500 µm, which is a typical value for microlens arrays. In y-direction, the additional length is determined by the length of the Y-branches, the crossings and the bends as shown in Fig. 1. In total, this results in an additional path length of approximately $N \times 850$ µm. The loss contributions independent of the number of channels are dominated by the single-channel photonic building blocks because the propagation lengths of most sections can be minimized by adapting the PIC dimension. The only exception is the section from the cascading tree up to the sample ports measuring about 700 µm due to the MMI length and bends from the MMI to the PDs, with the PD spacing being half of the channel spacing. The losses of photonic building blocks that have to be included are from the DC, both in terms of insertion loss and amount of light directed to the reference path (only−0.46 dB or $\approx$90% towards sample), the two PBS, in terms of insertion and crosstalk loss, and the two inverted tapers at the input and sample port. In summary, the power is reduced from the source to the sample by

$$\begin{aligned}\Gamma_{\textrm{Source}}[\textrm{dB}] & = 2\Gamma_{\textrm{Inverted taper}}+\Gamma_{\textrm{DC}}+0.46\;\textrm{dB}+2(\Gamma_{\textrm{PBS}}+\textrm{CT}_{\textrm{PBS}})+700\;\mathrm{\mu}\textrm{m} \times \textrm{PL} \\ & + N\times(\Gamma_{\textrm{Y-branch}}+3\;\textrm{dB}+\Gamma_{\textrm{Crossing}}+850\;\mathrm{\mu}\textrm{m} \times \textrm{PL}) \\ & \approx 3.313\;\textrm{dB} + N \times 3.244\;\textrm{dB}. \end{aligned}$$
The loss $\Gamma _{\textrm {Source}}$ results from light propagating in TE-like polarization. For the loss of signal light three photonic building blocks are relevant – the inverted taper, PBS and MMI. Since the components can be connected directly to each other only the connection bends from the MMI to the photodiodes are relevant, which are approximately 350 µm in length. In conclusion, the loss from the sample to the photodiode is calculated with
$$\Gamma_{\textrm{Signal}}[\textrm{dB}] = \Gamma_{\textrm{Inverted taper}}+ \Gamma_{\textrm{PBS}}+\textrm{CT}_{\textrm{PBS}}+\Gamma_{\textrm{MMI}}+350\;\mathrm{\mu}\textrm{m}\times \textrm{PL} \approx 1.502\;\textrm{dB},$$
from light propagating in TM-like polarization.

For the calculation of the required source power the intended power at the sample is divided by the power reduction $\Gamma _{\textrm {Source}}$. We assume a single beam power of 0.5 mW to guarantee laser safe operation at 840 nm according to the norms ISO 15004 and ANSI Z80.36.. In the case that the single beam is scanned over the retina a higher power of 0.75 mW would be allowed resulting in an average power of less than 0.5 mW. While the higher allowed power was used in the single-channel measurement, for a multi-channel system a more restrictive consideration in regard to safety is recommendable. For the safety of a multi-beam ophthalmic system the combined effect of the beams has to be taken into account. While spatially separated beams at the retina can be considered independently, at the cornea and lenticle the beams may overlap and the power may add up. The norms give a safe mean intensity of 100 mW/cm2 averaged over a defined area at the corneal plane. According to the ANSI norm more than sixteen single beams with 0.5 mW each can safely overlap at the cornea. Table 3 shows the values for the required source power to reach 0.5 mW at each probing position and the effective A-scan rate for an increasing number $N$ of cascading steps (channels).

As can be seen in the table, already for two channels the effective A-scan rate exceeds the 100 kHz figure of merit. For up to 32 channels the laser power stays below 50 mW, which is realistically achievable for a swept source with a booster amplifier. The effective A-scan rate shows that the imaging can be done with a sixteen-channel OCT system at approximately ten times the speed of a lossless single-channel system with the same SNR. Therefore, even in the comparison with the unrealistic case of a lossless system the multi-channel system can achieve higher acquisition rates.

Tables Icon

Table 3. Performance analysis of multi-channel configurations taking into account losses of the designed building blocks for the center wavelength of 840 nm. The first column shows the required source power to achieve an optical power of 0.5 mW for each probing beam at the retina. The second column shows the resulting effective A-scan rate. The SNR for all systems is identical to the lossless single-channel system because the sweeping rate is adapted accordingly for the lossy single-channel system and no additional loss is introduced on the way from the sample to the interferometer for higher number of channels.

5. Outlook

The next step towards a multi-channel system is to combine the designed building blocks with optoelectronics and electronics, e.g. photodiodes and analog to digital converter (ADC), on the same chip. The inclusion of the reference path in the PIC, with the potential capability to change the length, would be a further step towards integration. However, this would require significantly lower propagation losses and dispersion. The dispersion issue can be addressed by changing the dispersion properties e.g. by changing the cross section and/or the cladding material [5]. To facilitate the numerical dispersion compensation a PIC-based k-clock can be included for accurate real-time information on the dispersion properties, consisting of an asymmetric Mach-Zehnder interferometer with a well defined physical path length difference.

6. Conclusion

In conclusion, we have presented a concept for a low-loss PIC-based multi-channel SS-OCT system utilizing the intrinsic advantages of PICs such as high integration density and mechanical stability enabling parallelization. In comparison, bulk or fiber optic multi-channel implementations are significantly more intricate to achieve with growing difficulties and losses for each additional measurement channel. As a proof of concept of the polarization based optical path routing concept, which forms the basis for our low-loss multi-channel concept, we designed and fabricated a PIC for a single-channel SS-OCT system operating at 840 nm with a bandwidth of 100 nm including all major photonic building blocks. We showed that the performance of this PIC-based OCT system is sufficient to achieve in-vivo retinal OCT tomograms. This work represents the first PIC-based SS-OCT demonstration in the 840 nm wavelength domain. With the designed and tested photonic building blocks we carried out a performance analysis of multi-channel configurations. We showed that in principle the laser source power levels can be kept reasonable, while achieving a tenfold speed increase compared to an ideal single-channel system with the same SNR. This can be done even though excess losses in PICs are higher than in bulk and fiber optics. This shows that PICs have the potential to play an important role in future OCT development, not only in regard to size and cost but also to performance.

Funding

H2020 LEIT Information and Communication Technologies (688173, 871312).

Acknowledgments

This research has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 688173 (OCTCHIP) and No 871312 (HandheldOCT). Portions of this work were presented at the SPIE OPTO in 2019 ”Silicon-nitride waveguide-based integrated photonic circuits for medical diagnostic and other sensing applications”, the SPIE Optics + Optoelectronincs in 2019 ”PECVD silicon nitride optical waveguide devices for sensing applications in the visible and <1 µm near infrared wavelength region” and the ECIO in 2020 ”PECVD SiN photonic integrated circuit for swept source OCT at 840 nm”. The authors want to thank Dr. Michael Kempe and Daniel Bublitz of Carl Zeiss AG for their input on the laser safety norms.

Disclosures

The authors declare no conflicts of interest.

References

1. W. Drexler and J. G. Fujimoto, Optical Coherence Tomography (Springer International Publishing, 2015), Chap. 1.

2. J. F. de Boer, R. Leitgeb, and M. Wojtkowski, “Twenty-five years of optical coherence tomography: the paradigm shift in sensitivity and speed provided by Fourier domain OCT [Invited],” Biomed. Opt. Express 8(7), 3248–3280 (2017). [CrossRef]  

3. R. Leitgeb, C. Hitzenberger, and A. Fercher, “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Express 11(8), 889–894 (2003). [CrossRef]  

4. G. Yurtsever, P. Dumon, W. Bogaerts, and R. Baets, “Integrated photonic circuit in silicon on insulator for Fourier domain optical coherence tomography,” Proc. SPIE 7554, 75541B (2010). [CrossRef]  

5. G. Yurtsever, K. Komorowska, and R. Baets, “Low dispersion integrated Michelson interferometer on silicon on insulator for optical coherence tomography,” Proc. SPIE 8091, 80910T (2011). [CrossRef]  

6. V. D. Nguyen, N. Weiss, W. Beeker, M. Hoekman, A. Leinse, R. G. Heideman, T. G. van Leeuwen, and J. Kalkman, “Integrated-optics-based swept-source optical coherence tomography,” Opt. Lett. 37(23), 4820–4822 (2012). [CrossRef]  

7. S. Schneider, M. Lauermann, C. Weimann, W. Freude, and C. Koos, “Silicon photonic optical coherence tomography system,” in 2014 Conference on Lasers and Electro-Optics (CLEO) - Laser Science to Photonic Applications, (IEEE, 2014), pp. 1–2.

8. G. Yurtsever, N. Weiss, J. Kalkman, T. G. van Leeuwen, and R. Baets, “Ultra-compact silicon photonic integrated interferometer for swept-source optical coherence tomography,” Opt. Lett. 39(17), 5228–5231 (2014). [CrossRef]  

9. Z. Wang, H.-C. Lee, D. Vermeulen, L. Chen, T. Nielsen, S. Y. Park, A. Ghaemi, E. Swanson, C. Doerr, and J. Fujimoto, “Silicon photonic integrated circuit swept-source optical coherence tomography receiver with dual polarization, dual balanced, in-phase and quadrature detection,” Biomed. Opt. Express 6(7), 2562–2574 (2015). [CrossRef]  

10. L. Chang, N. Weiss, T. G. van Leeuwen, M. Pollnau, R. M. de Ridder, K. Wörhoff, V. Subramaniam, and J. S. Kanger, “Chip based common-path optical coherence tomography system with an on-chip microlens and multi-reference suppression algorithm,” Opt. Express 24(12), 12635–12650 (2016). [CrossRef]  

11. S. Schneider, M. Lauermann, P.-I. Dietrich, C. Weimann, W. Freude, and C. Koos, “Optical coherence tomography system mass-producible on a silicon photonic chip,” Opt. Express 24(2), 1573–1586 (2016). [CrossRef]  

12. Y. Huang, M. Badar, A. Nitkowski, A. Weinroth, N. Tansu, and C. Zhou, “Wide-field high-speed space-division multiplexing optical coherence tomography using an integrated photonic device,” Biomed. Opt. Express 8(8), 3856–3867 (2017). [CrossRef]  

13. M. S. Eggleston, F. Pardo, C. Bolle, B. Farah, N. Fontaine, H. Safar, M. Cappuzzo, C. Pollock, D. J. Bishop, and M. P. Earnshaw, “90dB sensitivity in a chip-scale swept-source optical coherence tomography system,” in Conference on Lasers and Electro-Optics, OSA Technical Digest (online) (Optical Society of America, 2018), paper JTh5C.8.

14. T. G. van Leeuwen, I. B. Akca, N. Angelou, N. Weiss, M. Hoekman, A. Leinse, and R. G. Heideman, “On-chip Mach-Zehnder interferometer for OCT systems,” Adv. Opt. Technol. 7(1-2), 103–106 (2018). [CrossRef]  

15. D. Huang, P. Pintus, C. Zhang, P. Morton, Y. Shoji, T. Mizumoto, and J. E. Bowers, “Dynamically reconfigurable integrated optical circulators,” Optica 4(1), 23–30 (2017). [CrossRef]  

16. P. Pintus, D. Huang, P. A. Morton, Y. Shoji, T. Mizumoto, and J. E. Bowers, “Broadband TE optical isolators and circulators in silicon photonics through Ce:YIG bonding,” J. Lightwave Technol. 37(5), 1463–1473 (2019). [CrossRef]  

17. D. Dai, “Advanced passive silicon photonic devices with asymmetric waveguide structures,” in Proceedings of the IEEE (IEEE, 2018), pp. 2117–2143.

18. R. Hainberger, P. Muellner, S. Nevlacsil, A. Maese-Novo, F. Vogelbacher, M. Eggeling, J. Schotter, M. Sagmeister, G. Koppitsch, and J. Kraft, “Silicon-nitride waveguide-based integrated photonic circuits for medical diagnostic and other sensing applications,” Proc. SPIE 10922, 1092204 (2019). [CrossRef]  

19. M. Sagmeister, G. Koppitsch, P. Muellner, S. Nevlacsil, A. Maese-Novo, R. Hainberger, D. Seyringer, and J. Kraft, “Monolithically integrated, CMOS-compatible SiN photonics for sensing applications,” Proceedings 2(13), 1023 (2018). [CrossRef]  

20. S. Nevlacsil, M. Eggeling, P. Muellner, G. Koppitsch, M. Sagmeister, J. Kraft, and R. Hainberger, “Broadband SiN asymmetric directional coupler for 840 nm operation,” OSA Continuum 1(4), 1324–1331 (2018). [CrossRef]  

21. Z. Lu, Y. Wang, F. Zhang, N. A. F. Jaeger, and L. Chrostowski, “Wideband silicon photonic polarization beamsplitter based on point-symmetric cascaded broadband couplers,” Opt. Express 23(23), 29413–29422 (2015). [CrossRef]  

22. R. Halir, A. Ortega-Monux, I. Molina-Fernandez, J. G. Wanguemert-Perez, P. Cheben, D.-X. Xu, B. Lamontagne, and S. Janz, “Compact high-performance multimode interference couplers in silicon-on-insulator,” in IEEE Photonics Technology Letters, (IEEE, 2009), pp. 1600–1602.

23. R. Hainberger, P. Muellner, S. Nevlacsil, F. Vogelbacher, M. Eggeling, A. Maese-Novo, M. Sagmeister, G. Koppitsch, J. Kraft, X. Zhou, J. Huang, M. Li, K.-J. Jiang, and Y. Song, “PECVD silicon nitride optical waveguide devices for sensing applications in the visible and <1µm near infrared wavelength region,” Proc. SPIE 11031, 110310A (2019). [CrossRef]  

24. S. Nevlacsil, P. Muellner, M. Sagmeister, J. Kraft, and R. Hainberger, “Broadband low loss and ultra-low crosstalk waveguide crossings based on a multimode interferometer for 840 nm operation,” OSA Continuum 3(2), 334–344 (2020). [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 (6)

Fig. 1.
Fig. 1. Conceptual drawing of a PIC-based four-channel OCT system. The input light is coupled to the PIC with a fiber via an inverted taper. The input polarization is set to transverse electric (TE)-like polarization by minimizing the light measured at the photodiode (PD) at the termination port of the first polarization beam splitter (PBS). The optical power is split into the reference and sample path via a directional coupler (DC) with an asymmetric splitting ratio and into the individual channels with a cascade of Y-branches. On both the sample and the reference path the outgoing light passes a PBS, which redirects the returning light with the orthogonal polarization towards the multimode interferometer (MMI). The polarization rotation from the forward propagating TE-like mode to the backward propagating transverse magnetic (TM)-like mode is achieved outside of the PIC with a $\lambda /4$-plate rotated by 45$^{\circ }$. On its way from the $\lambda /4$-plate to the sample and reference mirror the light is circularly polarized. The returning reflected light is converted into the orthogonal polarization relative to the outgoing light by passing the $\lambda /4$-plate a second time. The returning reference light crosses the sample path twice, which can be achieved efficiently with low-loss and low-crosstalk dual polarization crossings (DPCs). The highly divergent beams from the PIC are collimated with a microlens array (MLA). The PIC only provides short parts of the sample and reference paths avoiding high propagation losses and wavelength dispersion effects.
Fig. 2.
Fig. 2. (a) Waveguide cross section used for routing on the PIC. Normalized electric field amplitude distribution of light with (b) TE-like and (c) TM-like polarization. (d) Corresponding measured propagation losses in the wavelength range of 790 nm to 890 nm.
Fig. 3.
Fig. 3. Schematics of the broadband photonic building blocks required for low-loss single-channel SS-OCT system: (a) inverted taper for efficient bulk/fiber optic component coupling to and from the PIC; (b) multimode interferometer (MMI) optimized for TM-like polarization to interfere incoming light from the sample and reference paths; (c) directional coupler to split the optical input path into sample and reference path with 90% of the optical power on the through port and 10% on the cross port; (d) polarization beam splitter for redirecting the returning light in the orthogonal polarization state towards the interferometer with TE-like polarization suppressed towards the cross port and vice versa for the TM-like polarization.
Fig. 4.
Fig. 4. Measurement (crosses) and simulation (lines) results of photonic building blocks: (a) insertion loss for the inverted taper with sinusoidal width envelope function (averaged from three different samples); (b) insertion loss and imbalance of the MMI, both determined with a symmetric Mach-Zehnder interferometer comprising two MMIs; (c) polarization extinction ratio (PER), and (d) losses of the PBS for both polarizations. The simulations for the MMIs were made using the eigenmode expansion (EME) method in FIMMPROP (Photon Design), whereas the simulations for the PBS were performed with the 3D FDTD method.
Fig. 5.
Fig. 5. Schematic of PIC for single-channel SS-OCT system: The input laser light is coupled to the TE-like mode on the PIC from the west side via the inverted taper. The grating coupler located at the suppression port of the PBS provides optical feedback for adjusting an external polarization controller to optimize the input polarization. The TE-like polarized light is then separated by the directional coupler with 90% of the light going into the sample path (north) and 10% going into the reference path (south). Both in the sample and the reference path a PBS redirects the returning polarization rotated light towards the MMI. The interfered light exits the output ports on the east side, which are spatially separated to avoid cross talk between the output ports. The length of the sample and the reference paths are matched to avoid wavelength dispersion mismatch on the PIC.
Fig. 6.
Fig. 6. Tomogram (averaged over 400 consecutive B-scans) of a healthy human retina acquired with the PIC-based single-channel SS-OCT system. The layer structure of the retina can be clearly seen. The axial resolution appears to be lower than the determined value of 5.5 µm because of the high number of averages resulting in a blurring of the image due to insufficient motion correction of the micro movements of the eye ball.

Tables (3)

Tables Icon

Table 1. List of previously reported PIC-based SS-OCT systems.

Tables Icon

Table 2. Summary of the propagation loss (PL), insertion loss ( Γ ), and crosstalk (CT) values for the individual building blocks to be used in the proposed multi-channel SS-OCT system. All values are given in dB for the central wavelength of 840 nm. They have been experimentally determined with fabricated devices except the ones indicated with S , which have been derived from 3D FDTD simulations.

Tables Icon

Table 3. Performance analysis of multi-channel configurations taking into account losses of the designed building blocks for the center wavelength of 840 nm. The first column shows the required source power to achieve an optical power of 0.5 mW for each probing beam at the retina. The second column shows the resulting effective A-scan rate. The SNR for all systems is identical to the lossless single-channel system because the sweeping rate is adapted accordingly for the lossy single-channel system and no additional loss is introduced on the way from the sample to the interferometer for higher number of channels.

Equations (5)

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

IL [ dB ] = 10 log 10 ( P 1 + P 2 ) 2 ,
ER [ dB ] = 10 log 10 ( P 1 / P 2 ) = 10 log 10 ( 4 ρ + ρ 1 2 cos ( 2 θ ) ) ,
IB max [ dB ] = 10 log 10 ( 1 + 2 η ( 1 + 1 + η ) ) ,
Γ Source [ dB ] = 2 Γ Inverted taper + Γ DC + 0.46 dB + 2 ( Γ PBS + CT PBS ) + 700 μ m × PL + N × ( Γ Y-branch + 3 dB + Γ Crossing + 850 μ m × PL ) 3.313 dB + N × 3.244 dB .
Γ Signal [ dB ] = Γ Inverted taper + Γ PBS + CT PBS + Γ MMI + 350 μ m × PL 1.502 dB ,
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.