## Multi-mode Time-delay Interferometer for Free-space Quantum Communication

Clinton Cahall, 1, ∗ Nurul T. Islam, 2 Daniel J. Gauthier, 2 and Jungsang Kim 1, 3

1 Department of Electrical and Computer Engineering, Duke University, Durham, NC 27708, USA 2 Department of Physics, The Ohio State University, Columbus, OH 43210, USA 3 IonQ, Inc., College Park, MD 20740, USA

(Dated: August 5, 2019)

Quantum communication schemes such as quantum key distribution (QKD) and superdense teleportation provide unique opportunities to communicate information securely. Increasingly, optical communication is being extended to free-space channels, but atmospheric turbulence in free-space channels requires optical receivers and measurement infrastructure to support many spatial modes. Here we present a multi-mode, Michelson-type time-delay interferometer using a field-widened design for the measurement of phase-encoded states in free-space communication schemes. The interferometer is constructed using glass beam paths to provide thermal stability, a field-widened angular tolerance, and a compact footprint. The performance of the interferometer is highlighted by measured visibilities of 99 02 . ± 0 05 %, . and 98 38 . ± 0 01 % for single- and multi-mode inputs, . respectively. Additionally, high quality multi-mode interference is demonstrated for arbitrary spatial mode structures and for temperature changes of ± 1 0 . ◦ C. The interferometer has a measured optical path-length drift of 130 nm / ◦ C near room temperature. With this setup, we demonstrate the measurement of a two-peaked, multi-mode, single-photon state used in time-phase QKD with a visibility of 97 37 . ± 0 01 %. .

## I. INTRODUCTION

Optical communication protocols using time-bin or pulse-position modulation encoding are a good fit for free-space communication systems. Temporal states having bin widths on the order of a few hundred pico-seconds are robust against errors caused by dispersion from the atmosphere [1], while achieving high data rates. One significant challenge of establishing a free-space optical link is the distortion of the spatial mode of the beam at the receiver due to atmospheric turbulence. A freespace quantum communication scheme, such as quantum key distribution (QKD) using a time-phase encoding [2, 3], requires a time-delay interferometer designed to support many spatial modes to perform measurements in the phase-basis.

promising design for enabling high-quality multi-mode interference is a field-widened Michelson-type interferometer without relay optics [7, 8]. One such design has reported 91 % interference visibility and 85 % entanglement visibility [6]. In this article we discuss the design, construction, and characterization of a multi-mode timedelay interferometer operating at 1550 nm. The design presented here uses a Michelson-type layout without relay optics to achieve robust multi-mode interference. Our design features different glass materials in each beam path to achieve improved spatial mode and thermal performance in a compact layout. We demonstrate the application of the interferometer to a single-photon phasestate measurement used in a time-phase QKD scheme.

Previous work has demonstrated classical and quantum interference with multi-mode time-delay interferometers that leverage various methods to enable a high quality interference [4-6]. An optical beam with many spatial modes could potentially suffer degraded interference due to each mode having a slightly different wavevector k , and therefore a different propagation path length through the interferometer. Robust multi-mode interferometers are designed to minimize the path length difference for each mode to allow a complete interference for an arbitrary spatial mode.

One technique used to reduce the mode-dependent path length uses a 4 f - imaging system in each arm to relay the mode during propagation and therefore remove the k -vector dependence. This approach has demonstrated fringe visibility as high as 93% [5]. Another

∗ Corresponding author: clinton.cahall@duke.edu

## II. MULTI-MODE DELAY INTERFEROMETER DESIGN

Asimplified interferometer layout is shown in Fig. 1 [8]. A ray incident from the left at an angle θ is split by a non-polarizing beamsplitter and each beam is directed to one of two pathways. One path has a short propagation length relative to the other. The short arm is made of material with index n 1 and thickness d 1 , and the long arm is made of material with index n 2 and thickness d 2 . Optionally, an air gap of thickness d 3 can be introduced between the end of the material in the long arm and the reflector. The two beams are reflected at the end of each arm and recombine at the beamsplitter.

The optical path length difference ∆ between the two arms of the interferometer as a function of incident angle θ is given by [7]

FIG. 1. Simplified schematic of a Michelson-type interferometer. This setup includes a small air gap at the end of the long arm ( d 2 ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-b138f477-bc6b-43eb-b813-f033968ca306.jpg)

$$\Delta & = 2 [ n _ { 3 } d _ { 3 } ( 1 - \sin ^ { 2 } \theta / n _ { 3 } ^ { 2 } ) ^ { 1 / 2 } \\ & \quad + n _ { 2 } d _ { 2 } ( 1 - \sin ^ { 2 } \theta / n _ { 2 } ^ { 2 } ) ^ { 1 / 2 } \\ & \quad - n _ { 1 } d _ { 1 } ( 1 - \sin ^ { 2 } \theta / n _ { 1 } ^ { 2 } ) ^ { 1 / 2 } ],$$

where n 3 = 1 but is kept in the equation for completeness. Typically, a Michelson interferometer is designed for the chief ray to enter at normal incidence and small angular deviations from the chief ray are analyzed. In the design presented here, the input ray must enter at a finite angle θ 0 to spatially separate each output beam from the input beam. Access to both outputs are required for the time-phase QKD protocol, discussed later. Unlike the analysis performed in [7], where Eq. 1 is expanded about a normal incidence, we perform the expansion of Eq. 1 about a finite input angle θ 0 . The first three terms of this expansion are given by

$$\Delta & = 2 \left ( d _ { 3 } \beta _ { 3 } + d _ { 2 } \beta _ { 2 } - d _ { 1 } \beta _ { 1 } \right ) \\ - 2 \Phi ( \sin \theta - \Phi ) \left ( \frac { d _ { 3 } } { \beta _ { 3 } } + \frac { d _ { 2 } } { \beta _ { 2 } } - \frac { d _ { 1 } } { \beta _ { 1 } } \right ) \quad ( 2 ) \\ - ( \sin \theta - \Phi ) ^ { 2 } \left ( \frac { n _ { 3 } ^ { 2 } d _ { 3 } } { \beta _ { 3 } ^ { 3 } } + \frac { n _ { 2 } ^ { 2 } d _ { 2 } } { \beta _ { 2 } ^ { 3 } } - \frac { n _ { 1 } ^ { 2 } d _ { 1 } } { \beta _ { 1 } ^ { 3 } } \right ), \quad \text{ when} \quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{
\end{cases}$$

where Φ = sin θ 0 is the input angle and β 2 i = n 2 i -Φ 2 is the effective index of refraction. The goal of our analysis is to satisfy three conditions that are derived from Eq. 2. A full mathematical derivation can be found in Ref. [9].

## III. DESIGN REQUIREMENTS

The interferometer is designed to achieve the following; (1) an optical path-length difference (OPD) of 200 ps (5 GHz free spectral range), (2) satisfy the field-widened condition (discussed more below) for a large angular bandwidth, and (3) be thermally compensated. A successful design will feature high visibility near unity, longterm temperature stability over a modest temperature range ( ± 0 2 . ◦ C), and have phase-tuning ability over a complete interference fringe.

The OPD of the chief ray is determined by setting sin θ = Φ, leaving only the first term of Eq. 2. The path length difference ∆ is in units of length, and therefore must be converted to a time delay ∆ t by dividing it by the speed of light in vacuum c and given by

$$\Delta _ { t } = 2 ( d _ { 3 } \beta _ { 3 } + d _ { 2 } \beta _ { 2 } - d _ { 1 } \beta _ { 1 } ) / c. \text{ \quad \ \ } ( 3 )$$

A wide angular bandwidth is necessary for a multi-mode interferometer design because a multi-mode beam has a range of wavevectors and each will take a slightly different path through the interferometer. To achieve high visibility, the change in OPD as a function of input angle must be minimized. Field-widening is achieved by proper material choices and device geometry that make the second term in Eq. 2 vanish. The field-widening condition is

$$\frac { d _ { 3 } } { \beta _ { 3 } } + \frac { d _ { 2 } } { \beta _ { 2 } } - \frac { d _ { 1 } } { \beta _ { 1 } } = 0.$$

The result of satisfying Eq. 4 is that the OPD in Eq. 2 becomes first-order insensitive to deviations from the input angle of the chief ray [7, 8].

Lastly, thermal compensation provides long-term stability in an imperfectly controlled environment. Temperature fluctuations cause OPD changes due to expansion and contraction of the physical lengths of the glass material, as well as temperature-dependent changes in the index of refraction [7]. The thermal drift is characterized by taking the partial derivative of Eq. 2 with respect to temperature T , as

$$\frac { \partial \Delta } { \partial T } & = 2 d _ { 2 } \left [ \beta _ { 2 } \alpha _ { 2 } + \frac { n _ { 2 } } { \beta _ { 2 } } \left ( \frac { \partial n _ { 2 } } { \partial T } \right ) \right ] \\ & - 2 d _ { 1 } \left [ \beta _ { 1 } \alpha _ { 1 } + \frac { n _ { 1 } } { \beta _ { 1 } } \left ( \frac { \partial n _ { 1 } } { \partial T } \right ) \right ],$$

where α i = (1 /d i )( ∂d /∂T i ) is the coefficient of thermal expansion (CTE). This analysis assumes that the index and length change of the air gap with temperature is negligible and therefore those terms vanish when taking the derivative. The interferometer is made a-thermal by setting ∂ ∆ /∂T = 0, imposing a constraint on the path lengths d 1 and d 2 , given by

$$\begin{array} { c } \left [ \beta _ { 1 } \alpha _ { 1 } + \frac { n _ { 1 } } { \beta _ { 1 } } \left ( \frac { d n _ { 1 } } { d T } \right ) \right ] d _ { 1 } = \left [ \beta _ { 2 } \alpha _ { 2 } + \frac { n _ { 2 } } { \beta _ { 2 } } \left ( \frac { d n _ { 2 } } { d T } \right ) \right ] d _ { 2 }. \ ( 6 ) \end{array}$$

## IV. PERFORMANCE SIMULATION

We conducted a survey of optical glass with the goal of identifying two glass types to construct the interferometer (See Appendix A) . The basic guidelines we used for

making the glass choices are: (1) availability, (2) a large difference in index of refraction between the two materials, and (3) small relative changes in refractive index and length with temperature. A large difference in the index of refraction will enable smaller physical lengths to satisfy a given time-delay.

The finalized design uses N-SF66 and N-BK10 (Schott AG) for the long arm and short arm respectively. At the design wavelength of 1550 nm, N-BK10 has an index of refraction of n 1 = 1 4823, . ∂n /∂T 1 = 1 4154 . × 10 -6 K -1 , and a CTE of α 1 = 5 8 . × 10 -6 K -1 . N-SF66 has an index of refraction of n 2 = 1 8660, . ∂n /∂T 1 = -1 5575 . × 10 -6 K -1 , and a CTE of α 1 = 5 9 . × 10 -6 K -1 . The index and change in index for each material are calculated at λ = 1550nm using the Sellmeier equation and its derivative with respect to temperature. The glass rods are mounted to a substrate of Kovar nickel alloy. The CTE of Kovar is α Kovar = 5 86 . × 10 -6 K -1 and is closely matched to the glasses used in our design. A good match between the glass and substrate, rather than using an ultra-low expansion (ULE) substrate with a much smaller CTE, is advantageous because matching CTE values between the glass and substrate will minimize the stress on the glass during temperature changes, which can cause birefringence and refraction.

The lengths d 1 and d 2 are optimized for an input angle θ 0 = 3 ◦ and satisfy the OPD, field-widening, and athermal conditions. However, satisfying these three conditions simultaneously leads to large physical dimensions of the interferometer ( &gt; 50 mm). The negative consequence of large glass pathways are: (1) large glass path lengths increase the thermal mass to be stabilized and can increase the settling time for phase tuning: (2) as the path length increases so does the clear aperture required for a beam input at a fixed input angle; and (3) longer path-lengths will be more sensitive to overlap misalignment at the beamsplitter. Misalignment is especially sensitive for multi-mode interference because the spatial frequency of the mode can easily be on the order of the alignment tolerance, whereas in single-mode interference, the mode size is large compared to the alignment tolerance.

The goal of the design presented here is to minimize the physical dimensions of the glass beam paths. Minimizing the glass length will allow for a more robust and manufacturable setup. To meet this goal, we relax the thermal requirements of the design so that the a-thermal condition is satisfied only approximately, while the OPD and field-widening conditions are satisfied exactly. We find a unique solution that satisfies Eq. 3 and Eq. 4 by eliminating the air gap d 3 . The OPD of the chief ray and field-widening condition for an interferometer without an air gap d 3 are given by

$$\Delta _ { t } = 2 ( d _ { 2 } \beta _ { 2 } - d _ { 1 } \beta _ { 1 } ) / c, \text{ \quad \ \ } ( 7 ) \text{ \ \ sys }$$

and

$$\frac { d _ { 2 } } { \beta _ { 2 } } - \frac { d _ { 1 } } { \beta _ { 1 } } = 0,$$

respectively. The a-thermal condition remains unchanged when removing the air gap.

FIG. 2. Simulated change in the optical path difference as a function of deviations of the input about the chief ray of the multi-mode interferometer.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-2036cc7f-f5fa-40a8-b73f-77bdc9c8a4c0.jpg)

The absence of an air gap leaves two equations and two free parameters, d 1 and d 2 . The path lengths that satisfy Eq. 7 and Eq. 8 are d 1 = 34 5 mm and . d 2 = 43 5 mm. . At an input angle of θ 0 = 3 ◦ and a beam diameter ∼ 2 mm (discussed in more detail below), the required clear aperture is ∼ 10 mm. A clear aperture of this size allows the square cross-section of the glass rods to be 12 7 mm . × 12 7 mm, and leads to manageable total ther-. mal mass. Additionally, a 10 mm clear aperture is &lt; 80 % of a 12.7 mm surface, which allows each surface to be polished easily with standard glass processing capabilities.

The angular bandwidth of the interferometer is characterized by plotting the OPD in Eq. 2 for a range of input angles θ around the chief ray Φ = sin 3 ◦ , shown in Fig. 2. Recall that satisfying the field-widening condition given in Eq. 8 implies that the first term in the angular dependence of the OPD (Eq. 2) is proportional to sin 2 θ .

The thermal stability of the setup is estimated by calculating Eq. 5 for small temperature changes ∆ T , shown in Fig. 3(a). The simulation shows that thermal drift in the OPD of ≤ λ/ 50 is achievable for modest temperature fluctuations of ± 0 1 . ◦ C. The phase of the interference is adjusted by controlling the temperature in one of the glass arms, the short arm for example, to ± 1 ◦ C. Temperature control over this range will result in a change in the OPD from -λ/ 2 to λ/ 2, allowing for a tuning range of 0 to 2 π for the phase φ , as shown in Fig. 3(b). Thermal tuning is preferred over mechanical tuning because mechanical tuning techniques such as piezo-driven mirrors can be expensive, require an extra control and feedback system, and have large thermal drift.

FIG. 3. Simulated thermal performance of the multi-mode interferometer. (a) Thermal stability for small temperature drifts of the complete system. (b) Interference phase tuning via control over the temperature of the short arm d 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-1ae7e512-5dbc-43bb-b2aa-d826ac2054dd.jpg)

## V. INTERFERENCE CHARACTERIZATION

The performance of the interferometer is verified using optical beams with single-mode (SM) and multimode (MM) intensity profiles. Atmospheric turbulence is emulated by mode mixing in a perturbed multi-mode fiber [6, 10]. The experimental setup, shown schematically in Fig. 4, consists of a continuous-wave laser diode (Fitel F0L15DCWC) operating at 1550 nm with a laser controller (Arroyo 6310 Combo Source) that controls the diode current and the temperature of the package. The wavelength of the output for diodes of this type is a function of the diode current and temperature. The controller has an analog voltage input for manual adjustment of the current supplied to the laser diode. A controllable voltage source is used as the scan voltage input for diode current sweeps. Adhesive heater tape (Omega KHLVA101/5-P) is used on the short arm to provide phase tuning of the interferometer. The light from the laser is directed to the interferometer via a single or multi-mode optical fiber. The output of the interferometer is coupled in to an optical fiber where the output power is measured with an optical power meter (Newport Optical Power Meter 1803C with Model 818-IG sensor head), or the intensity profile of the beam can be viewed with a beam profiling camera (Ophir-Spricon SP907-1550).

In a practical setting, the temperature of the interferometer is meant to be passively controlled by modest isolation from the surrounding lab environment. For the

FIG. 4. Schematic of the experimental setup used to characterize the interference quality and thermal performance of the interferometer (gray box). SM: single-mode, MM: multimode, TEC: thermoelectric cooler, ADC: analog-to-digital converter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-785d55e0-f220-4c4d-8b08-7b27fdb2e890.jpg)

purposes of thermal characterization, the experimental setup includes precise temperature control with thermoelectric cooling elements (TEC, Laird 64975-502) and a TEC controller (Thorlabs TED200C) with proportionalintegral-derivative (PID) feedback control. Three TEC elements are placed between the bottom of the interferometer base plate and a large thermal bath, which is the optical table in this case. The temperature of the base plate is sensed via high-accuracy thermistors (Vishay NTCALUG02A103F161) placed at various locations on the plate. A thermistor near the center of the plate is used for the TEC feedback control. The analog output from the TEC controller is digitized with an analog-to-digital converter (ADC, National Instruments NI-9239) and recorded. A representative example of the settling time and temperature stability of the interferometer is shown in Fig. 5.

FIG. 5. Temperature control of the interferometer bottom plate. (a) Settling time of the plate temperature after a step change of 2 ◦ C. The target temperature is stable after ∼ 3 minutes. (b) Long-term temperature stability at a fixed set point. The temperature is controlled with an accuracy of ± 0 0015 . ◦ C.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-b8322d4f-bac3-47dd-b2af-7cceabad241e.jpg)

FIG. 6. Interference phase tuning with short arm temperature control. The output from the interferometer is plotted as a function of the square of the voltage supplied to the heater element on the short arm showing a tuning range over a full interference fringe.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-479df380-f1e8-4336-ab5b-2e8529eb960b.jpg)

## A. Fringe Visibility

The phase of the interference fringe is scanned by heating one of glass arms and the result is shown in Fig. 6. To maintain repeatability and execute measurements on faster timescales than thermal settling allows, interference fringes are recorded by scanning the input wavelength rather than scanning the phase of the interferometer. The wavelength of the laser diode is adjusted with an analog voltage input to the laser controller, which controls the current supplied to the diode. As the laser current is scanned, the wavelength changes but so does the total output power. The output power as a function of the laser current is first calibrated by measuring the power directly from the laser. The change in optical power as a function of diode current is compensated for by dividing the interference fringe measurements by the calibration measurement. The power output of the interferometer is recorded as a function of the input scan voltage value. The quality of the interference is quantified by the visibility V , given by

$$\mathcal { V } = \frac { \mathcal { P } _ { + } - \mathcal { P } _ { - } } { \mathcal { P } _ { + } + \mathcal { P } _ { - } } \, \quad \, \quad \, ( 9 )$$

where P + is the maximum and P -is the minimum value of the interference fringe.

The power P out , ± measured at the + / -port of the interferometer as a function of small changes of the wavelength δλ glyph[lessmuch] ∆, when the phase of the interference is set to φ , is given by [11]

$$\frac { P _ { \text{out}, \pm } } { P _ { 0 } } = \frac { a } { 2 } \left [ 1 \mp \sin \left ( k \delta \lambda + \phi \right ) \right ] + b, \quad \left ( 1 0 \right ) \quad \mathcal { V } _ { \text{M} }$$

FIG. 7. Interference fringes for SM and MM (50 µ mcore fiber) input beams as a function of the input voltage to scan the laser diode current. Each data point (black) is the average of 10 measurements and error bars indicate the standard deviation of those 10 measurements. Red lines show a fit to the data using the function given by Eq. 10. (a) SM interference fringe. The visibility of the interference is V SM = 99 02 . ± 0 05 %. The . high order fringe structure in the constructive interference could be due to reflections from multiple surfaces within the interferometer. (b) Zoomed view of the SM fringe minimum. (c) MM interference fringe. The visibility of the interference is V MM = 98 38 . ± 0 01 %. . (d) Zoomed view of the MM fringe minimum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-8f85b253-a9e7-47f3-be52-af641d97264d.jpg)

where P 0 is the input optical power, a ∈ { 0 1 , } is the insertion loss of the interferometer, b accounts for imperfect destructive interference in the minimum of the fringe, and k = 2 π/λ is the wavenumber of the light. The output power is recorded as a function of the laser diode scan voltage, which serves as the wavelength adjustment parameter. The fringe measurement is fit with a curve given by Eq. 10 with a , b , and φ as fitting parameters. The visibility of the normalized fringe is extracted from this fit where

$$\mathcal { V } = \frac { a } { a + 2 b } \, \quad \quad \quad ( 1 1 )$$

The measurement results from a single-mode input fiber (Corning SMF-28-Ultra) and multi-mode input fiber (50 µ m core, 0.1 NA, Coastal Connections MFUFU-50sL) are shown in Fig. 7. The interferometer demonstrates single-mode interference visibility V SM = 99 02 . ± 0 05 % . and multi-mode interference visibility V MM = 98 38 . ± 0 01 % extracted from the fitted . functions. The single-mode results are close to those reported

by high-performing commercial interferometers such as those used in a recent QKD demonstration [3, 11]. Additionally, the multi-mode results presented here are significantly better than those recently reported for multi-mode interferometers [5, 6].

## B. Spatial Mode Performance

A robust multi-mode interferometer is required to maintain high visibility independent of the input spatial mode structure. The multi-mode results presented above are optimized and performed on a single beam profile. To investigate the robustness of the performance, we characterize the visibility of the interference fringe as the spatial profile is changed. The visibility is maximized for the first spatial profile by tuning the beam alignment such that there is good overlap at the beamsplitter. Then, interference fringes are recorded for subsequent spatial profiles after the mode has been changed. The difference between spatial mode profiles is quantified by the two dimensional Pearson correlation coefficient R defined as

$$R = \frac { \sum _ { i = 1 } ^ { N _ { x } } \sum _ { j = 1 } ^ { N _ { y } } \left [ I _ { i j } \left ( t _ { 1 } \right ) - \mu \left ( t _ { 1 } \right ) \right ] \left [ I _ { i j } \left ( t _ { 2 } \right ) - \mu \left ( t _ { 2 } \right ) \right ] } { \sigma ( t _ { 1 } ) \sigma ( t _ { 2 } ) }, \quad \begin{array} { c } \underset { \left ( 1 2 \right ) } { \sigma ( t _ { 1 } ) } \\ \end{array} \text{(using)} }, \quad \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \text{(using)} } \end{cases}$$

where I ij ( t n ) is the optical intensity collected at time t n (a frame), indexed by a discrete array of positions ( i, j ) that correspond to the pixels of the camera, µ t ( n ) is the mean intensity, and σ t ( n ) is the standard deviation of the intensity distribution. The correlation coefficient R has a value between +1 and -1, where +1( -1) indicates perfect positive (negative) correlation, and values close to 0 indicate no correlation.

Optical intensity profiles are collected by the beam profiling camera and the correlation coefficient is processed in MATLAB. The spatial mode of the beam is changed by introducing a controlled bend in the optical fiber with a micrometer and positioning stage. A schematic depicting the setup to induce changes in the mode is shown in Fig. 8(a). The fiber is anchored at three different points with the middle point on the positioning stage. As the position of the stage is adjusted, the fiber is bent at two different points that leads to the change in the mode profile. Special care is taken to ensure that all other points in the fiber are kept stationary throughout the experiment. When the fiber is kept secure apart from the micrometer adjustment, the spatial modes can be reproduced as a function of the micrometer position.

The profile of the beam is measured at the output of the fiber (without collimation) using the beam profiler. Examples of the intensity patterns from 50 µ m and 105 µ mcore (Thorlabs FG105LVA) fibers for different micrometer positions are shown in Fig. 8(b) and (c). Note that the experiment is carried out by first recording the mode profiles with the beam profiler at each micrometer

FIG. 8. Spatial mode characterization from a 50 µ m and 105 µ m core fiber. (a) Schematic showing the setup used to adjust the mode profile of the beam. The spatial mode is changed by inducing bends in the fiber controlled by a positioning stage and a micrometer. (b) Spatial mode examples for different micrometer settings from a 50 µ m core fiber. (c) Spatial mode examples for different micrometer settings from a 105 µ m core fiber.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-13376d9d-8893-4df5-b84b-17772750e80e.jpg)

position, then performing the interference measurements at each micrometer position. Between the beam profiler measurements and the interference measurements, the output of the fiber must be moved from the beam profiler to the interferometer input. The movement of the fiber will unavoidably change the actual spatial mode that travels through the interferometer from those in the recorded images in Fig. 8(b) and (c). However, the intensity profile and the correlation between profiles at different micrometer positions recorded by the beam profiler are representative of the actual mode passing through the interferometer. Therefore, the mode analysis using the beam profiler adequately represents the spatial mode at the input of the interferometer.

The correlation coefficient and interference visibility as a function of micrometer position are shown in Fig. 9. The correlation between the first measurement and subsequent measurements decays to a value &lt; 0 5 . as the spatial mode becomes more distinguishable. As the spatial mode is changed, the interferometer maintains excellent fringe visibility for each measurement. The quality of the interference using a 50 µ m core fiber has an average value V avg = 98 05 . ± 0 01% and a maximum . value V max = 98 38 . ± 0 01%. . Additionally, the quality of the interference using a 105 µ m core fiber has an average value V avg = 97 55 . ± 0 04% and a maximum value .

10 -5

4

3

2

1

0

1

0.99998

0.99996

0.99994

FIG. 9. Multi-mode interference as a function of the input spatial mode. Correlation coefficient between the first frame and subsequent frames as a function of the micrometer setting for 50 µ m (circles) and 105 µ m (triangles) core fibers is shown against the left vertical axis. Also plotted is the visibility of the interference fringe over the same micrometer range for each fiber, shown in the right vertical axis. Error in the fringe visibility is ≤ 5 × 10 -4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-a4babd0c-6b70-4020-b6ab-2e149e470e93.jpg)

V max = 98 14 . ± 0 02%. The quality of these results indi-. cate an interferometer design that is robust for arbitrary input spatial modes.

## C. Thermal Performance

The thermal performance of the interferometer is characterized by measuring the visibility of the interference fringe and the path length shift as a function of temperature. The temperature of the bottom plate is set and maintained by the TEC controller that supplies current to the TEC elements in contact with the substrate of the interferometer. When the setpoint is changed, the system is allowed to come to thermal equilibrium over a period of ∼ 30 minutes. As thermal equilibrium is reached, the path length will drift as a function of time, causing a change in the measured output power. Therefore, when a new setpoint is investigated, the output of the interferometer is monitored as equilibrium is reached. We infer that the interferometer has come to thermal equilibrium once the output power is stabilized, shown in Fig. 10. During setpoint changes, the phase of the interference is set such that it lies near the point of the steepest slope to have the maximum sensitivity in determining the equilibrium.

Once the system has come to thermal equilibrium, a full interference fringe is scanned and fit with the function given by Eq. 10 in the same manner as discussed above. A fringe is recorded at each temperature setting and plotted as a function of the scanning voltage applied to the diode controller, corresponding to the change in

FIG. 10. Example measurement of the output power from the interferometer while the system reaches thermal equilibrium. This particular set of data shows a temperature change from 21 0 . ◦ C to 21 5 . ◦ C

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-8f34c438-722a-40f9-ab6f-1e76501c1c7f.jpg)

the laser wavelength. The results of the curve fitting are shown in Fig. 11. Path length drift brought on by the temperature change causes the fringe to move horizontally. While operating at a fixed value of φ (shown by the 'X' on the yellow curve in the figure) a horizontal shift of the fringe manifests in a change in the power measured at the output port of the interferometer. The change in output power describes change in path length.

The path length change and visibility are measured for a temperature range of ± 1 0 . ◦ C centered at 22 0 . ◦ C. The results of the measurements are shown in Fig. 11(b) and (c). Ahigh fringe visibility of ∼ 98 % is maintained across the entire temperature range. The experimentally determined path-length shift as a function of temperature is 130 nm/ ◦ C. This value matches well with the theoretically predicted performance. Therefore, stability on the order of λ/ 50 can be achieved with moderate temperature control in a lab environment on the order of 0 25 . ◦ C.

## VI. APPLICATION TO TIME-PHASE QKD

Experimental realizations of QKD systems are often hindered by non-idealalities in the equipment used to generate and measure the photons. Rate limitations are dominated by detector saturation and dead-time, and so it is advantageous to increase the channel capacity by encoding bits in a high dimensional Hilbert space. We use a time-phase encoding [2] because it can support a high dimensional encoding. An encoding scheme having d -dimensions is often referred to as a qudit protocol, where d = 2 is the well-known qubit protocol. Qudits can transmit up to log 2 d bits per photon, whereas qubit schemes, such as BB84 [12], are limited to a single bit per photon. Additionally, high-dimensional schemes can

FIG. 11. Multi-mode interference as a function of temperature. (a) Fitted functions to the measured interference fringe at each temperature setting. The temperature settings for the curves moving from left to right are 21 0 . , 21 5 21 9 . , . , 22 0 22 5 . , . , and 23 0 . ◦ C. The phase of the interference at 22 0 . ◦ C is set at data point marked with an 'X.' The path length shift is calculated using points where the fringe intersects the dotted line. (b) Experimentally measured (black) and theoretically predicted (red) path length change as a function of temperature. (c) Interference visibility as a function of temperature.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-edb75f2f-5c22-4689-ad54-d91eb64392a7.jpg)

tolerate a higher bit-error-rate, which enable longer communication distances [13]. Recently, a high-dimensional QKD scheme using a d = 4 time-phase encoding has achieved a record high maximum secret-key-rate (SKR) of 26.2 Mbit/s, and an SKR of 1.07 Mbit/s at an equivalent fiber distance of 83 km [3].

A time-phase QKD protocol consist of two mutually unbiased bases that are referred to as the temporal (or time) basis and the phase (or frequency) basis. Data is encoded in the temporal basis and transmitted from Alice to Bob, and the phase basis is used to check for the presence of an eavesdropper, Eve. In general, d -1 cascaded interferometers are needed to perform a complete measurement of the phase basis [14]. A full description of such a protocol, including a demonstration utilizing a high-dimensional encoding, can be found in Ref. [3]. Here, we use a d = 2 encoding to demonstrate the performance of our multi-mode interferometer.

The states for a d = 2 time-phase protocol are shown in Fig. 12. Temporal states are constructed by intensity modulation of a continuous-wave (CW) laser source to create short duration wavepackets. The states are highly attenuated to produce weak coherent states with a mean photon number per pulse of µ ≈ 1. Phase states are constructed using intensity modulators to create the superposition of time states and an additional phase modulator to adjust the phase of the wavepackets in each time bin. The phase states are attenuated so that the probability amplitudes of the entire frame sum to unity. It is important to use a modulation scheme with a high extinction ratio to prevent errors caused by background photons leaking through the modulators.

A measurement in the temporal basis is performed with a single-photon detector, where the electrical signal generated from a detection event is recorded with a high resolution time-tagger. The single-photon detector and processing electronics, including amplifiers and time-tagging electronics, must have low timing uncertainty (detection jitter) to localize the time-of-arrival of the photon with high accuracy [15, 16]. The arrival time of the single-photon detection event is compared to the shared global clock to declare the result of the measurement.

Phase-state measurements require a time-delay interferometer, shown schematically in Fig. 12(c). A phase state entering the interferometer is split by a 50/50 beamsplitter with half of the wavepacket traveling through a long path and half through a short path. The difference in path lengths is chosen such that the time-delay is equal to τ , or one time bin width. Therefore, the wavepackets recombining at the second beamsplitter are shifted by one time bin relative to one another. An additional phase adjustment φ is available to tune the phase of the interference fringe, where φ = 0 in the d = 2 scheme.

A phase state that passes through the interferometer will result in constructive or destructive interference in the shaded time-bin shown in Fig. 12(e). There is a one-to-one mapping between the phase-state input and the detector in which constructive interference is observed. Therefore, constructive interference is detected at a unique detector, with destructive interference at all other detectors. In the d = 2 protocol, a photon detection event in the shaded time bin of detector D 0 ( D 1) indicates a measurement outcome of | f 0 〉 ( | f 1 〉 ). The lobes in the interference pattern lying outside the shaded time bin are considered losses and do not contribute to the outcome of the measurement. The phase-basis measurement succeeds at a rate of 1 /d , or 50% for the d = 2 scheme. The security of the channel is assessed at the end of all communication by calculating the visibility of the constructive and destructive interference fringes. The visibility of the interference is given by

$$\mathcal { V } = \frac { \mathcal { P } _ { 0 } - \mathcal { P } _ { 1 } } { \mathcal { P } _ { 0 } + \mathcal { P } _ { 1 } }$$

where P 0 ( P 1 ) is the summation of counts recorded by detector D 0 ( D 1) in the shaded time bin. The visibility of interference has a direct impact on the length of the secret key that is distilled from the raw key. For the d = 2 time-phase scheme discussed here, the quantum bit-error-rate (QBER) of a phase-basis measurement is

FIG. 12. Time-phase protocol with d = 2. (a) Temporal and (b) phase basis states. (c) Time-delay interferometer and it's action on a single wavepacket. (d) Interferometer setup for a d = 2 phase basis measurement. (e) Interference patterns for a measurement on the | f 0 〉 state.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-06d71e8f-575b-4c7a-b14b-0fe783f206ca.jpg)

1.81607

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-004b6f7e-b8e3-496c-929e-6847fc73ca1b.jpg)

1.81608

1.81609

1.8161

FIG. 13. Interference characterization and spatial mode for phase-state measurement. (a) Output from the plus (minus) port of the interferometer shown in blue (black) with Eq. 10 fitted to each curve shown in red. The visibility between the plus and minus port with the phase set at the maximum of the plus port is 98 7 %. . (b) Example of the spatial mode under test measured at the output of a 25 µ m core fiber.

given by [13]

$$Q B E R = \frac { 1 - \mathcal { V } } { 2 }. \quad \quad \ \ ( 1 4 ) \quad \ t h \cdot \quad \ \ \\ \text{an}$$

As the visibility of the interference increases the QBER decreases and therefore a larger secret key can be distilled from the raw key.

Experimentally, the generation of the states shown in Fig. 12(a) and (b) are carried out in a similar manner to the method discussed in Ref. [11]. A continuouswave laser operating at 1550 nm is intensity modulated with electro-optic modulators to create short duration ( ∼ 80 ps) wavepackets. The modulators are driven by a 10 GHz transceiver from a field-programmable gate array (FPGA). The mean photon number for each wavepacket is controlled using a variable optical attenuator (VOA). The modulation and attenuation are performed in singlemode fiber, then the light is coupled to a multi-mode fiber with a 25 µ mcore (Thorlabs FG025LJA) to produce wavepackets with multiple spatial modes. Mode mixing is achieved by stressing the fiber in the same manner as shown in Fig. 8(a). A 25 µ m core fiber is chosen over the 50 µ m and 105 µ m core fibers mentioned previously because of limitations in our single-photon detection system. The light is sent through the interferometer, where each output is collected with a 25 µ m core fiber. The fibers are then coupled to our single-photon detection system.

1.81607 1.81608 1.81609 1.8161 Our detection system is comprised of amorphous superconducting nanowire single-photon detectors (SNSPDs) produced by Quantum Opus. SNSPDs are high speed and high efficiency single photon counting devices (see Ref. [17] for a review) that are widely used in quantum communication [3, 18] and quantum optics experiments [19]. Though large area devices [20] and pixel arrays [21, 22] designed to support free-space or multi-mode fiber coupling exist, our devices have an active area with a diamter of ∼ 15 µ m and are designed to couple to a single-mode telecommunications fiber (eg. Corning SMF-28) [23]. However, in this experiment, we couple the 25 µ m core fibers from the output of the interferometer. There will be some amount of aperture loss at the detector due to the mismatch in the size of the detector and the fiber core. We choose to suffer the extra loss at the detector for the purpose of demonstrating multi-mode, single-photon interference. The fringes from the plus and minus port of the interferometer and an example of the spatial mode profile of the beam under test at the output of the fiber are shown in Fig. 13.

The preparation of the | f 0 〉 state is shown with singlephoton detection events in Fig. 14(a). The error rate of the preparation, calculated from the overlap of each wavepacket in to the opposite time-bin, is 0 4 % and is . largely due the limited bandwidth of the modulators used to create the wavepackets. The result of a measurement on the | f 0 〉 state using the time-delay interferometer is shown in Fig. 14(b). As expected, constructive interference is observed in the middle time-bin of detector D0, with destructive interference observed in the middle time-bin of detector D1. Additionally, the heights of the lobes to the left and right of the middle time-bin are approximately 1 / 4 the height of the constructive interference fringe, and match well with the predicted measurement outcome. The visibility between the two detectors in the middle time-bin is determined by summing all the counts falling in the time-bin window and using the formula given by Eq. 13. The resulting visibility is 97 37 . ± 0 01 %. . This value is within 1 % of the visibility reported for a similar phase-state measurement using

FIG. 14. Phase-state | f 0 〉 preparation and measurement using a 25 µ m core fiber. (a) Single-photon data (blue) with a sum of Gaussian functions fitted to the data (red) for the | f 0 〉 state preparation. Zeros above the data denote the phase of each bin, indicating a | f 0 〉 state. Time bin boundaries are shown with black dotted lines. (b) Measurement of the | f 0 〉 state using the multi-mode interferometer showing the plus (+) output port collected at detector D0 (blue) and minus ( -) output port collected at detector D1 (black). Time bin boundaries are shown with dotted lines. The resulting visibility between the bright and dark fringe in the middle time bin is 97 37 . ± 0 01 %. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-c91076a3-4ec3-4a44-85ec-4c115e6a276d.jpg)

commercial single-mode time-delay interferometers [11].

## VII. CONCLUSION

In this article, we detail the design methodology, performance simulations, construction and characterization of a time-delay interferometer with a 5 GHz free-spectralrange (200 ps time-delay). The interferometer is set-up in a Michelson-type layout with glass beam paths in each arm. A wide field-of-view and excellent thermal stability can be achieved with a proper material choice and optical design, while satisfying the target time-delay. The performance results are highlighted by a maximum singlemode visibility of V SM = 99 02 . ± 0 05 %, comparable to . commercial single-mode products, and multi-mode visibility of V MM = 98 38 . ± 0 01 %, . several percent higher than recent reports. Additionally, multi-mode visibility V MM ≈ 98 % is maintained independent of spatial mode profile for both 50 µ m and 105 µ m core fibers. Highquality interference is demonstrated over a modest temperature range of ± 1 ◦ C where the path-length shift of the interferometer is 130 nm/ ◦ C, enabling good thermal stability with modest temperature control.

Finally, we demonstrate the application of our interferometer to a time-phase QKD protocol, where the presence of an eavesdropper is quantified by the interference visibility between the two output ports of a time-delay interferometer. The high visibility of our interferometer enables phase-state measurements with a low QBER, allowing efficient extraction of a secret key. We demonstrate a phase-state measurement using weak coherent states with visibility of 97 37 . ± 0 01 % for 1550 nm pho-. tons from a 25 µ m core fiber.

Further improvements can be made to the design and construction of the interferometer to achieve better performance. In our estimation the largest contributing factor to the degradation of the multi-mode interference visibility is an error in the length of the glass rods during polishing. Each piece of glass is undersized by 400 -600 µ m. The implication of this error is that the actual dimensions of the glass no longer satisfy the field-widening condition and therefore result in degraded interference, especially for multi-mode beam profiles. We were able to compensate for this error by a coarse adjustment of the position of one of the reflecting mirrors. Using a more precise method for measuring the position of the mirror, we expect the quality of multi-mode visbility could be improved to over 99 %.

## ACKNOWLEDGMENTS

This work is supported by the Office of Naval Research (ONR) (N00014-13-1-0627) and the National Aeronautics and Space Administration (NASA) (NNX13AP35A).

## Appendix A: Material Choice

TABLE I. Selection of glass options investigated for the construction of the interferometer.

| Glass Name                                   | Index n                                                        | CTE (10 - 6 /K                              | dn/dT (10 - 6 /K                                      | Mfr.                                                                  |
|----------------------------------------------|----------------------------------------------------------------|---------------------------------------------|-------------------------------------------------------|-----------------------------------------------------------------------|
| S-FSL5 N-BK10 N-BK7 S-NSL3 S-TIH6 N-SF6 SF57 | 1.4732 1.4823 1.5009 1.5025 1.5512 1.7628 1.7631 1.7997 1.8019 | 9.0 5.8 7.1 9.0 8.0 8.9 9.0 7.4 8.3 6.7 6.1 | -1.1 1.4 1.1 0.3 1.3 -0.4 -2.3 -0.8 6.0 -0.2 4.1 -1.5 | Ohara Schott Schott Ohara Ohara Ohara Schott Ohara Schott Ohara Ohara |
| S-BAL14                                      |                                                                |                                             |                                                       |                                                                       |
| S-NPH53                                      |                                                                |                                             |                                                       |                                                                       |
| S-NPH2                                       | 1.8614                                                         |                                             |                                                       |                                                                       |
| L-LAH86                                      | 1.8615                                                         |                                             |                                                       |                                                                       |
| N-SF66                                       | 1.8660                                                         | 5.9                                         |                                                       | Schott                                                                |
| S-LAH79                                      | 1.9552                                                         | 6.0                                         | 7.0                                                   | Ohara                                                                 |
| N-LASF35                                     | 1.9747                                                         | 8.5                                         | 2.7                                                   | Schott                                                                |

A survey of glass options was conducted and narrowed

to the list of materials given in Table I. The factors that led to the final decision of materials used are a combination of performance, availability, and cost.

## Appendix B: Interferometer Design and Construction

FIG. 15. CAD schematic of the interferometer with input and output beam steering and coupling optics. Overhead CAD view (a) and perspective view (b) of the interferometer layout. (c) Photograph of the assembled interferometer. The heater tape is shown on the short arm glass and several temperature sensors are placed at various positions on the Kovar plate. (d) Detailed view of the beamsplitter and glass rod alignment features, and the substrate cutout for thermal isolation of the short arm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00852v1-figure-f7162559-0515-4482-895b-10d7991b083a.jpg)

- [1] L. Kral, I. Prochazka, and K. Hamal, Optical signal path delay fluctuations caused by atmospheric turbulence, Optics letters 30 , 1767 (2005).
- [2] T. Brougham, C. F. Wildfeuer, S. M. Barnett, and D. J. Gauthier, The information of high-dimensional time-bin encoded photons, European Physics Journal D 70 , 214 (2016).
- [3] N. T. Islam, C. C. W. Lim, C. Cahall, J. Kim, and D. J. Gauthier, Provably secure and high-rate quantum key distribution with time-bin qudits, Science Advances 3 , e1701491 (2017).
- [4] G. Vallone, D. Dequal, M. Tomasin, F. Vedovato, M. Schiavon, V. Luceri, G. Bianco, and P. Villoresi, Interference at the single photon level along satellite-ground

Construction of the interferometer is carried out with N-BK10 and N-SF66 glass rods measuring 12 7 mm . × 12 7 mm . × 34 5 mm and 12 7 mm . . × 12 7 mm . × 43 5 mm . respectively, where the long dimension is the optical direction. The end faces are optically polished and coated with an anti-reflective (AR) dielectric coating that is optimized for 1550 nm. The glass rods are mounted to the Kovar nickel alloy substrate.

A computer-aided design (CAD) model of the complete interferometer system with input and output fiber coupling and beam steering optics is shown in Fig. 15. The input fiber is collimated with a fixed collimating package (f = 8mm asphere, Thorlabs F240-1550) that is mounted in a custom-machined mount made of stainless steel. The optical plane is defined by the center of the collimating lenses and the center of the mirrors, both of which are 12.7 mm from the surface of the Kovar plate. The non-polarizing beamsplitter (Newport 05BC16NP.11) and glass rods are mounted on a platform to place the center of these pieces on the optical plane. The beamsplitter and glass rods are aligned on the substrate via 1 mm-height, 1 mm-width raised ridges.

## Appendix C: Deviation from Design

There are two important differences between the schematic in Fig. 1 and the simulated system discussed previously. The assembled system has 1 mm gaps between the beamsplitter and glass rods, and 1mm gaps between the end of the glass rods and the reflecting mirrors in each arm. Air gaps are inserted to avoid the need for optical bonds at each interface to reduce complexity and increase the production quality of the assembly. Additionally, air gaps between the beamsplitter and glass rods allow each face to be AR coated and reduce the loss at these interfaces. Symmetric air gaps at the interfaces do not contribute to the OPD. Expansion and contraction of the substrate due to temperature changes will also be symmetric across each spacer gap and do not invalidate the mathematical analysis presented.

channels, Physical review letters 116 , 253601 (2016).

- [5] C. Zeitler, T. Graham, H. Bernstein, P. Kwiat, and J. Chapman, Super-dense teleportation for space applications, in SPIE LASE (International Society for Optics and Photonics, 2016) pp. 973912-973912.
- [6] J. Jin, S. Agne, J.-P. Bourgoin, Y. Zhang, N. L¨tkenhaus, u and T. Jennewein, Demonstration of analyzers for multimode photonic time-bin qubits, Physical Review A 97 , 043847 (2018).
- [7] G. G. Shepherd, W. A. Gault, D. Miller, Z. Pasturczyk, S. F. Johnston, P. Kosteniuk, J. Haslett, D. J. Kendall, and J. Wimperis, Wamdii: wide-angle michelson doppler imaging interferometer for spacelab, Applied Optics 24 , 1571 (1985).

- [8] W. A. Gault, S. F. Johnston, and D. J. Kendall, Optimization of a field-widened michelson interferometer, Applied Optics 24 , 1604 (1985).
- [9] C. Cahall, Enabling Technologies for High-Rate, FreeSpace Quantum Communication , Ph.D. thesis, Duke University (2019).
- [10] R. Ursin, F. Tiefenbacher, T. Schmitt-Manderbach, H. Weier, T. Scheidl, M. Lindenthal, B. Blauensteiner, T. Jennewein, J. Perdigues, P. Trojek, B. Omer, M. Furst, M. Meyenburg, J. Rarity, Z. Sodnik, C. Barbieri, H. Weinfurter, and A. Zeilinger, Entanglement-based quantum communication over 144 km, Nature Physics 3 , 481 (2007).
- [11] N. T. Islam, C. Cahall, A. Aragoneses, A. Lezama, J. Kim, and D. J. Gauthier, Robust and stable delay interferometers with application to d-dimensional timefrequency quantum key distribution, Physical Review Applied 7 , 044010 (2017).
- [12] C. H. Bennettand G. Brassard, Quantum cryptography: Public key distribution and coin tossing, in Proceedings of IEEE International Conference on Computers, Systems and Signal Processing (1984) pp. 175-179.
- [13] V. Scarani, H. Bechmann-Pasquinucci, N. J. Cerf, M. Duˇ sek, N. L¨tkenhaus, and M. Peev, The security of u practical quantum key distribution, Reviews of Modern Physics 81 , 1301 (2009).
- [14] T. Brougham, S. M. Barnett, K. T. McCusker, P. G. Kwiat, and D. J. Gauthier, Security of high-dimensional quantum key distribution protocols using franson interferometers, Journal of Physics B: Atomic, Molecular and Optical Physics 46 , 104010 (2013).
- [15] F. Marsili, V. B. Verma, J. A. Stern, S. Harrington, A. E. Lita, T. Gerrits, I. Vayshenker, B. Baek, M. D. Shaw, R. P. Mirin, and S. W. Nam, Detecting single infrared photons with 93% system efficiency, Nature Photonics 7 , 210 (2013).
- [16] C. Cahall, D. J. Gauthier, and J. Kim, Scalable cryogenic readout circuit for a superconducting nanowire singlephoton detector system, Review of Scientific Instruments

89 , 063117 (2018).

- [17] E. A. Dauler, M. E. Grein, A. J. Kerman, F. Marsili, S. Miki, S. W. Nam, M. D. Shaw, H. Terai, V. B. Verma, and T. Yamashita, Review of superconducting nanowire single-photon detector system design options and demonstrated performance, Optical Engineering 53 , 081907 (2014).
- [18] D. M. Boroson, B. S. Robinson, D. V. Murphy, D. A. Burianek, F. Khatri, J. M. Kovalik, Z. Sodnik, and D. M. Cornwell, Overview and results of the lunar laser communication demonstration, in Free-Space Laser Communication and Atmospheric Propagation XXVI , Vol. 8971 (International Society for Optics and Photonics, 2014) p. 89710S.
- [19] L. K. Shalm, E. Meyer-Scott, B. G. Christensen, P. Bierhorst, M. A. Wayne, M. J. Stevens, T. Gerrits, S. Glancy, D. R. Hamel, M. S. Allman, et al. , Strong loophole-free test of local realism, Phys. Rev. Lett. 115 , 250402 (2015).
- [20] E. E. Wollman, V. B. Verma, A. D. Beyer, R. M. Briggs, B. Korzh, J. P. Allmaras, F. Marsili, A. Lita, R. Mirin, S. Nam, et al. , Uv superconducting nanowire singlephoton detectors with high efficiency, low noise, and 4 k operating temperature, Optics express 25 , 26792 (2017).
- [21] V. B. Verma, R. Horansky, F. Marsili, J. A. Stern, M. D. Shaw, A. E. Lita, R. P. Mirin, and S. W. Nam, A fourpixel single-photon pulse-position array fabricated from wsi superconducting nanowire single-photon detectors, Applied Physics Letters 104 , 051115 (2014).
- [22] M. S. Allman, V. B. Verma, M. Stevens, T. Gerrits, R. D. Horansky, A. E. Lita, F. Marsili, A. Beyer, M. D. Shaw, D. Kumor, R. P. Mirin, and S. W. Nam, A near-infrared 64-pixel superconducting nanowire single photon detector array with integrated multiplexed readout, Applied Physics Letters 106 , 192601 (2015).
- [23] A. J. Miller, A. E. Lita, B. Calkins, I. Vayshenker, S. M. Gruber, and S. W. Nam, Compact cryogenic self-aligning fiber-to-detector coupling with losses below one percent, Opt. Express 19 , 9102 (2011).
Clinton Cahall Nurul T. Islam Daniel J. Gauthier Jungsang Kim
Multi-mode Time-delay Interferometer for Free-space Quantum Communication
Quantum communication schemes such as quantum key distribution (QKD) and
superdense teleportation provide unique opportunities to communicate
information securely. Increasingly, optical communication is being extended to
free-space channels, but atmospheric turbulence in free-space channels requires
optical receivers and measurement infrastructure to support many spatial modes.
Here we present a multi-mode, Michelson-type time-delay interferometer using a
field-widened design for the measurement of phase-encoded states in free-space
communication schemes. The interferometer is constructed using glass beam paths
to provide thermal stability, a field-widened angular tolerance, and a compact
footprint. The performance of the interferometer is highlighted by measured
visibilities of $99.02\pm0.05\,\%$, and $98.38\pm0.01\,\%$ for single- and
multi-mode inputs, respectively. Additionally, high quality multi-mode
interference is demonstrated for arbitrary spatial mode structures and for
temperature changes of $\pm1.0\,^{\circ}$C. The interferometer has a measured
optical path-length drift of $130\,$nm$/\,^{\circ}$C near room temperature.
With this setup, we demonstrate the measurement of a two-peaked, multi-mode,
single-photon state used in time-phase QKD with a visibility of $97.37\pm
0.01\,\%$.