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

Athermal silicon photonic wavemeter for broadband and high-accuracy wavelength measurements

Open Access Open Access

Abstract

We propose and demonstrate an integrated wavemeter capable of accurate and broadband measurements without control or knowledge of the temperature. In our design, interferometers composed of silicon and silicon nitride waveguides enable accurate measurements of an input optical wavelength despite large and rapid temperature fluctuations of 20°C by leveraging the disparity in thermo-optic properties of the waveguides. We derive formulas which resolve the wavelength and temperature ambiguity of the interferometers. The fabricated wavemeter chip is found to have a mean accuracy of 11 pm over an 80 nm range near 1550 nm. To our knowledge, this is the first demonstration of an athermal silicon wavemeter and the lowest measurement error across such a broad wavelength range using silicon photonics. This result may reduce the cost and size of wavemeters used in combination with integrated lasers for optical communications, sensing, and other applications.

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

1. Introduction

Determining the precise wavelength of a laser is important for applications such as optical communications, spectroscopy, and metrology [13]. Wavelength meters (wavemeters) provide such measurements, enabling calibration and stabilization of lasers for these and other applications. Wavemeter implementations have generally relied on either interferometers, optical beating, or wavelength-dependent structures [4]. Approaches using interferometers typically send the input laser light together with light from a reference laser into an interferometer (e.g. Michelson) and measure the output fringes while the optical path length is varied by mechanically moving an element, such as a mirror [57]. Another approach uses etalons and an imaging array to measure interference fringes [8,9], and other wavemeters similarly use interferograms or speckle patterns on camera sensors [1014].

Integrated wavemeters have recently been developed to reduce instrument size and cost and to eliminate moving parts. In one demonstration, an integrated wavemeter based on silicon waveguides used an arrayed waveguide grating (AWG) and photodiode array to achieve a detection error below 20 pm across a 20 nm range [15]. Another report involved modulating the index of silicon waveguides, achieving 73 pm root mean square (RMS) accuracy over 10 nm [16]. In another work, an indium phosphide (InP) wavemeter used two integrated Mach-Zehnder interferometers (MZIs) to attain 20 pm resolution over a 40 nm range [17]. However, these demonstrations face limitations due to temperature variations influencing the index of refraction in the waveguides, as silicon and InP have high thermo-optic coefficients (TOCs) of 1.8×10−4/°C and 2×10−4/°C, respectively. Recent works attempted to address this limitation by choice of waveguide material or confinement. In one paper, an interferometer for an on-chip wavemeter included a 62 cm length silicon nitride (Si3N4, TOC ≈ 2.4×10−5/°C) waveguide in order to produce 2.4 pm period fringes and reduced temperature dependence, but the design does not resolve wavelength ambiguity over a large range [18]. In [19], an integrated silicon wavemeter covered the C-band using multiple MZIs and leverages the low-confinement TM mode to reduce temperature sensitivity, but sub-degree temperature control is still required.

Several studies have proposed methods for achieving athermal devices. Passive athermality has been shown by using a material with a negative TOC to balance the waveguide core’s positive TOC [20,21] or by choosing a low-TOC material [22]. However, such materials are not typically available in silicon photonic foundry processes. Athermal MZIs were demonstrated by precisely controlling the waveguide dimensions, but the temperature-stability is relatively narrowband [23,24]. Ideally, a wavemeter could achieve both athermal and broadband performance with potential for integration alongside lasers on a single chip [2527].

Here we demonstrate an integrated athermal silicon wavemeter requiring neither control nor knowledge of the temperature [28]. Our wavemeter includes MZIs of varying delay lengths to cover a large wavelength range with high resolution, and by using both silicon and Si3N4 waveguides, we can isolate and subtract thermal effects in order to determine a laser’s wavelength at an unknown temperature. The wavemeter includes the MZIs in photonic circuits with integrated photodiodes on a single silicon chip, enabling wafer-scale manufacturing of low-cost wavemeters. Our calibrated wavemeter has a measured average accuracy of 11 pm across over 80 nm, which remains consistent for temperatures fluctuating from 20°C-41°C. The robustness, accuracy, and range of the wavemeter make it attractive for both optical communications and sensing.

2. Theory and design

2.1 Wavemeter using interferometers with 90-degree hybrids

Our integrated wavemeter uses asymmetric MZIs to measure the input signal’s wavelength. Each MZI splits its input into two arms, with the longer one including a path length difference $\mathrm{\Delta }L$ [ Fig. 1(a)]. The two arms connect to a 90-degree hybrid coupler to measure the phase difference between them [29]. This absolute phase difference $\varphi $ depends on the wavelength $\lambda $ and temperature T and is given by

$${\varphi = \Delta L{n_{\textrm{eff}}}({\lambda ,T} )\frac{{2\mathrm{\pi }}}{\lambda },}$$
where ${n_{\textrm{eff}}}({\lambda ,T} )$ is the effective index of the MZI waveguides. To retrieve the wavelength from the phase measurement, a linear approximation of the wavelength and temperature dependence around central values ${\lambda _0}$ and ${T_0}$ can be used:
$${{n_{\textrm{eff}}}({\lambda ,T} )\approx {n_{\textrm{eff}}}({{\lambda_0},{T_0}} )+ ({T - {T_0}} ){{\left. {\frac{{\partial {n_{\textrm{eff}}}}}{{\partial T}}} \right|}_{{\lambda _0},{T_0}}} + ({\lambda - {\lambda_0}} ){{\left. {\frac{{\partial {n_{\textrm{eff}}}}}{{\partial \lambda }}} \right|}_{{\lambda _0},{T_0}}}.}$$

 figure: Fig. 1.

Fig. 1. MZI phase and wavelength measurement. (a) Diagram of an integrated MZI for phase measurement. (b) Spectrum of a first MZI with a small $\mathrm{\Delta }{L_0}$ and large FSR. The red dots indicate the wavelength approximated by the photodetector current measurements ${X_0}$ and ${Y_0}$. (c) Relative phase inferred from the first MZI’s photocurrents. (d) Spectrum of a second MZI with a large $\mathrm{\Delta }{L_1}$ and small FSR. (e) Relative phase inferred from the second MZI’s photocurrents and from the order m as determined by the phase of the first MZI.

Download Full Size | PDF

To simplify the notation, we define the following variables:

$a \equiv {n_{\textrm{eff}}}({{\lambda_0},{T_0}} )$, the effective index at ${\lambda _0},{T_0}$,

$b \equiv {\left. {\frac{{\partial {n_{\textrm{eff}}}}}{{\partial \lambda }}} \right|_{{\lambda _0},{T_0}}}$, the first order wavelength dispersion,

${n_g} \equiv a - {\lambda _0}b$, the group index,

$\theta \equiv {\left. {\frac{{\partial {\textrm{n}_{\textrm{eff}}}}}{{\partial \textrm{T}}}} \right|_{{\lambda _0},{T_0}}}$, the TOC,

$\mathrm{\Delta }T \equiv T - {T_0}$, the temperature relative to ${T_0}$, and

$\mathrm{\Delta }\lambda \equiv \lambda - {\lambda _0}$, the wavelength relative to ${\lambda _0}$.

The phase difference at each interferometer output can be expressed as

$${\varphi = \Delta L\frac{{2\pi }}{\lambda }({a + b\mathrm{\Delta }\lambda + \theta \mathrm{\Delta }T} ).}$$
This can be rewritten using the group index:
$${\varphi = \Delta L\frac{{2\pi }}{\lambda }({{n_g} + \theta \mathrm{\Delta }T} )+ 2\pi b\Delta L.}$$
At the output of the 90-degree hybrid, photodetectors are configured in balanced pairs so that the two output electric currents X and Y are proportional to the two quadratures of the optical signal (Supplement 1, section A). A relative phase difference $\delta \varphi$ belonging to the interval $({ - \pi ,\; \pi } ]$ can be retrieved from the photodetector current measurements by the expression
$${\delta \varphi = \textrm{atan}2({Y,\; X} ).}$$
This phase measurement is independent of the input power to the MZI and provides a constant accuracy across the whole wavelength range of measurement (i.e., accuracy is not lost near interferometer nulls). The absolute and relative phase differences are related by
$${\varphi = \delta \varphi + 2\pi m - \pi /4,}$$
where the integer m is the order of the MZI. Finally, the wavelength of the input signal can be obtained from the measurements by the expression
$${\lambda = \Delta L\frac{{{n_g} + \theta \mathrm{\Delta }T}}{{\mathrm{\delta }\varphi /2\pi + m - 1/8 - \; b\mathrm{\Delta }L}}.}$$
However, we can see that this expression of $\lambda $ requires knowledge of m and T. Therefore, besides requiring accurate control or measurement of the temperature, a prior knowledge of the wavelength within a certain range is needed. Otherwise, the periodic spectral response of the MZI makes identifying the wavelength ambiguous.

The order of the MZI can be determined by using a second MZI. We previously demonstrated that the unambiguous wavelength range could be extended by cascading several interferometers [17]. An MZI with a small $\mathrm{\Delta }{L_0}$ [large free-spectral range (FSR)] can be used to identify the order of another MZI with a larger $\mathrm{\Delta }{L_1}$ (smaller FSR), as shown in Fig. 1(b-e). Figure 1(b) shows that the two photocurrent measurements from the first MZI (${X_0}$ and ${Y_0}$) identify a single wavelength and phase estimation [Fig. 1(c)]. The second MZI’s photocurrent measurements (${X_1}$ and ${Y_1}$) coincide at many wavelengths [red dots in Fig. 1(d)], making the wavelength determination ambiguous. However, the first MZI’s phase can determine the order of the second MZI, identifying a single, precise wavelength [Fig. 1(e)]. This allows the wavemeter to benefit from the large range of the former MZI while leveraging the higher accuracy of the second MZI. The temperature of the device must be accurately monitored and controlled, though, in order to maintain the accuracy of the wavelength measurements.

2.2 Athermal wavemeter

We remove the temperature dependence of the wavemeter by using MZIs with waveguides having different TOCs. Figure 2 shows a set of MZIs where the first one has waveguides with TOC θ2 and the second MZI has different waveguides with TOC θ3. By comparing the responses of each MZI subject to an unknown change in temperature, it is possible not only to determine the temperature, but also to account for the induced thermo-optic spectral shift when calculating the wavelength, as explained by the following derivation.

 figure: Fig. 2.

Fig. 2. Diagram of athermal wavemeter, composed of MZIs with different TOCs.

Download Full Size | PDF

The phase measured at each interferometer can be expressed as in Eq. (4). Because the phase is proportional to the optical frequency $\nu = c/\lambda $, we will now mainly use frequency rather than wavelength in the following calculations. The phase differences ${\varphi _2}$ and ${\varphi _3}$ can then be expressed as the system of two equations

$${{\varphi _2} = \Delta {L_2}\frac{{2\pi \nu }}{c}({{n_{g2}} + {\theta_2}\mathrm{\Delta }T} )+ 2\pi {b_2}\Delta {L_2}}$$
$${{\varphi _3} = \Delta {L_3}\frac{{2\pi \nu }}{c}({{n_{g3}} + {\theta_3}\mathrm{\Delta }T} )+ 2\pi {b_3}\Delta {L_3},}$$
where $\nu $ and $\mathrm{\Delta }T$ are the two unknown quantities. The Jacobian of this equation system near the point $({\nu _0},\; {T_0})$, is
$${J \equiv {J_\varphi }{{({\nu ,T} )}_{{\nu _0},\; {T_0}}} = \frac{{2\pi }}{c}\left( {\begin{array}{{cc}} {\mathrm{\Delta }{L_2}{n_{g2}}}&{\mathrm{\Delta }{L_2}{\theta_2}{\nu_0}}\\ {\mathrm{\Delta }{L_3}{n_{g3}}}&{\mathrm{\Delta }{L_3}{\theta_3}{\nu_0}} \end{array}} \right).}$$
For the equation system to provide a well-defined solution, the Jacobian determinant must be non-zero, which is equivalent to ${n_{g2}}/{\theta _2} \ne {n_{g3}}/{\theta _3}$. The solution of this system of equations provides an expression for the optical frequency $\nu $ as a function of the phases ${\varphi _2}$ and ${\varphi _3}$:
$${\nu = \frac{c}{{2\pi \mathrm{\Delta }{L_2}\mathrm{\Delta }{L_3}}}\frac{{{\theta _2}\mathrm{\Delta }{L_2}{\varphi _3} - {\theta _3}\mathrm{\Delta }{L_3}{\varphi _2} - 2\pi \mathrm{\Delta }{L_2}\mathrm{\Delta }{L_3}({{b_3}{\theta_2} - \; {b_2}{\theta_3}} )}}{{{\theta _2}{n_{g3}} - {\theta _3}{n_{g2}}}}.}$$
Notably, this calculation of the optical frequency does not require knowledge of the temperature. The other parameters such as TOCs, group indices, dispersion and delay lengths may be determined through a one-time calibration (Section 3.1).

The ultimate accuracy of the proposed wavemeter depends mainly on the TOCs and delay lengths. Assuming the main source of uncertainty comes from the electrical measurements of the detectors’ photocurrents, the uncertainty, i.e., the accuracy, of the optical frequency measurement can be expressed as (see Supplement 1, section B):

$${{\sigma _\nu } = \frac{c}{{2\pi \mathrm{\Delta }{L_2}\mathrm{\Delta }{L_3}}}\frac{{\sqrt {{{({\mathrm{\Delta }{L_2}{\theta_2}} )}^2} + {{({\mathrm{\Delta }{L_3}{\theta_3}} )}^2}} }}{{\; {\theta _2}{n_{g3}} - \; {\theta _3}{n_{g2}}}}{\sigma _\varphi },}$$
where ${\sigma _\varphi }$ is the phase uncertainty. We can see in Eq. (12) that, counter-intuitively, a large TOC ${\theta _2}$ becomes a desirable trait in our design. The larger the difference between the TOCs of the two waveguides (${\theta _2}$ and ${\theta _3}$), the smaller the optical frequency uncertainty will be. To determine ${\sigma _\varphi }$ in cases where the detector current uncertainty ${\sigma _D}$ is limited by the analog-to-digital converter (ADC) resolution, and assuming the range of the ADC fits the signal, the phase uncertainty becomes ${\sigma _\varphi } = {2^{ - \textrm{ENOB}}}$, where ENOB is the effective number of bits. We lastly note that the accuracy may also be impacted by other experimental effects not included in Eq. (12).

2.3 Design of the athermal silicon wavemeter

We implement the proposed athermal wavemeter on an integrated platform using silicon and Si3N4 MZIs. The wavemeter design is shown in Fig. 3(a). Input light from a laser is split using multimode interferometers (MMIs) and directed to four MZIs. One MZI has a small delay ΔL0 corresponding to a large FSR0, setting the maximum wavelength range of the wavemeter. The next MZI has an intermediate delay ΔL1, and the third and fourth MZIs have a long delay ΔL2= ΔL3. The first stage resolves the ambiguity in identifying the relevant order of the intermediate stage, and then the intermediate stage is used to identify the order of the long delay MZIs. The first three MZIs are composed of silicon, while the fourth uses Si3N4 waveguides. The phase differences at the outputs of the third and fourth MZIs can then be used to make temperature-independent wavelength measurements using Eq. (11), where θ2 and θ3 correspond to the TOCs of silicon and Si3N4 respectively. The phases are measured using 90-degree hybrid MMIs followed by integrated germanium photodiodes. Note that while we vary the TOC using different materials, it is also possible to vary the width of waveguides of a single material in order to obtain two TOCs. However, this alternative may not easily achieve a large TOC difference since waveguides that are too narrow suffer greater bending losses, and waveguides that are too wide become multimode.

 figure: Fig. 3.

Fig. 3. (a) Diagram of athermal wavemeter design. (b) Photograph of assembled device. (c) Microscope image of the wavemeter chip.

Download Full Size | PDF

We design the wavemeter’s MZIs to enable wavelength measurement over a large bandwidth with high accuracy. The delay lengths and FSRs at 1550 nm of the four MZIs are listed in Table 1. To cover a wide wavelength range for communications and sensing applications, the FSR of the first MZI is designed to be 94 nm, which corresponds to a 6 µm delay length difference. The following two MZIs successively decrease in FSR by a factor of ∼8.6, which empirically ensures reliable determination of the MZI orders, accounting for possible phase errors in the first stages. The final two MZIs, composed of silicon and Si3N4, have the finest resolution and are used to measure the wavelength according to Eq. (11). We estimate that an accuracy of 10 pm is achievable using these MZI dimensions, based on Eq. (12) and assuming a phase error of 0.02 rad.

Tables Icon

Table 1. Mach-Zehnder interferometer parameters

The integrated wavemeter is fabricated on silicon-on-insulator (SOI) wafers using a standard silicon photonic foundry process. The silicon waveguides are 450 nm×220 nm in size and can couple to the 1 µm×400 nm Si3N4 waveguides using tapered evanescent couplers. These waveguide dimensions are chosen for good confinement of the fundamental transverse electric (TE) mode and low bending loss. The 1×2 mm chip is assembled for testing by attaching an optical fiber to the input and wirebonding the electrical outputs to a fan-out board [Fig. 3(b-c)]. A thermoelectric cooler (TEC) is attached to the substrate underneath the chip along with an accompanying thermistor to assist in characterizing and calibrating the wavemeter’s thermal response.

3. Results

3.1 Wavemeter calibration

To perform a one-time calibration of the wavemeter, we extract linear fits of the phase responses of each MZI across the full wavelength range using an external tunable laser. First, we inject light from the laser, sweeping its wavelength across the measurement range while measuring the photodetector currents. We estimate an insertion loss of 5.5 dB to the detectors, including 2 dB coupling loss from the fiber. Figure 4 shows the current measured at the output of each balanced detector pair as a function of wavelength from 1500 to 1580 nm in 10 pm steps, at a fixed TEC temperature of 25°C. For each MZI, we obtain two sine curves in quadrature with a periodicity (FSR) depending on the optical delay. The distortions visible in the sinusoids likely originate in imperfections of the 90-degree hybrid MMIs and the response imbalance of the photodiodes. However, when we calculate the phase from these quadrature curves, the distortions largely disappear. The four top graphs in Fig. 5(a-d) show the phases unwrapped on the calibration range $\mathrm{\Delta }\varphi $ in radians, which are related to the absolute phases by

$${\varphi = \Delta \varphi + 2\pi {m_{\textrm{ref}}} - \pi /4,}$$
where ${m_{\textrm{ref}}}$ is the MZI order at the origin of the calibration range ${\nu _0}$ (corresponding to the shortest wavelength). Even though the x-axes use wavelength for convenience, we apply a linear fit to the measurements as a function of the optical frequency, i.e., $\mathrm{\Delta }{\varphi _i}(\nu )= {A_i} + {B_i}\nu $, where i is the respective MZI. The four bottom graphs [Fig. 5(e-h)] represent the error between the measurement and the linear fit, which is typically lower than 0.2 radians. Note that the linear fit errors of the first two MZIs have no effect on the accuracy of the final wavelength measurement, as long as they do not exceed the FSR of the successive stages (which is ensured through proper ratio of ΔL).

 figure: Fig. 4.

Fig. 4. Photodetector current measurements from each pair of balanced photodiodes (X and Y) for (a) MZI-0 (b) MZI-1 (c) MZI-2 and (d) MZI-3.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. (a-d) Measured unwrapped phases $\; \Delta \varphi $ for MZI-0 to MZI-3. (e-h) Measurement deviations from the linear fits of MZI-0 to MZI-3 phases.

Download Full Size | PDF

While the wavemeter operates athermally during regular operation, we must perform a one-time calibration to extract the temperature parameters. We accomplish this by measuring linear fits at several different temperatures. Each linear fit provides two parameters ${A_i}({\mathrm{\Delta }T} )$ and ${B_i}({\mathrm{\Delta }T} )$, such that

$${\mathrm{\Delta }{\varphi _i}({\mathrm{\Delta }\nu ,\mathrm{\Delta }T} )= {A_i}({\mathrm{\Delta }T} )+ {B_i}({\mathrm{\Delta }T} )\Delta \nu ,}$$
where $\mathrm{\Delta }\nu \equiv \nu - {\nu _0}$ and $\mathrm{\Delta }T \equiv T - {T_0}$ are the deviations of the optical frequency and temperature with respect to ${\nu _0}$ and ${T_0}$. We define ${T_0}$ as the mid-point in the temperature calibration range, 30°C. The parameters ${A_i}$ and ${B_i}$ are calculated at each temperature by measuring the MZI phases while sweeping the input laser, as in Fig. 5. We then fit the variation of ${A_i}$ and ${B_i}$ as a function of temperature (Fig. 6). The four top graphs [Fig. 6(a-d)] display the ${A_i}({\mathrm{\Delta }T} )$ coefficients with their linear fits, while the four bottom graphs [Fig. 6(e-h)] display the ${B_i}({\mathrm{\Delta }T} )$ coefficients with their linear fits. These new linear fits provide a set of parameters ${C_{ij}}$ such that ${A_i}({\mathrm{\Delta }T} )= {C_{i0}} + {C_{i1}}\mathrm{\Delta }T$ and ${B_i}({\mathrm{\Delta }T} )= {C_{i2}} + {C_{i3}}\mathrm{\Delta }T$. The entire calibration of the wavemeter can then be approximated with the 16 ${C_{ij}}$ parameters, and the four MZI phases can be approximated as
$${\mathrm{\Delta }{\varphi _i}({\mathrm{\Delta }\nu ,\mathrm{\Delta }T} )= {C_{i0}} + {C_{i1}}\Delta T + {C_{i2}}\Delta \nu + {C_{i3}}\Delta T\Delta \nu .}$$

 figure: Fig. 6.

Fig. 6. (a-d) Measurement and linear fit of A parameter versus temperature for MZI-0 to MZI-3. (e-h) Measurement and linear fit of B parameter across temperature.

Download Full Size | PDF

The ${C_{ij}}$ coefficients can be related to the physical parameters as follows:

$${{C_{i0}} = \Delta \varphi ({{\nu_0},{T_0}} )= \frac{{2\pi \mathrm{\Delta }{L_i}}}{c}{\nu _0}{n_{gi}} + 2\pi {b_i}\Delta {L_i} - 2\pi {m_{ref}} + \pi /4}$$
$${{C_{i1}} = {J_{i1}} = \frac{{2\mathrm{\pi }\Delta {L_i}}}{c}{\nu _0}{\theta _\textrm{i}}}$$
$${{C_{i2}} = {J_{i2}} = \frac{{2\mathrm{\pi }\Delta {L_i}}}{c}{n_{gi}}}$$
$${{C_{i3}} = \frac{{2\mathrm{\pi }\Delta {L_i}}}{c}{\theta _\textrm{i}} = \frac{{{C_{i1}}}}{{{\nu _0}}}.}$$
As can be seen in Fig. 6(e-h), parameter ${C_{i3}}$ has miniscule values and therefore reduced accuracy, but it can be replaced by its equivalent expression ${C_{i1}}/{\nu _0}$, where ${C_{i1}}$ and ${\nu _0}$ are accurately known values. To have a completely linearly approximated calibration of the wavemeter, one only needs the remaining 12 ${C_{ij}}$ parameters and ${\nu _0}$. These 13 values are sufficient to provide a good approximation of the optical frequency measurement during wavemeter operation. These coefficients can also be used in conjunction with the values in Table 1 to calculate physical parameters such as ${\theta _\textrm{i}}$ and ${n_{gi}}$. We again note that controlling and measuring the temperature using the experimental setup is not necessary during wavemeter operation, after the one-time calibration is completed. In fact, the wavemeter itself can measure the temperature (Section 3.3).

3.2 Wavemeter operation

We obtain the optical frequency from the MZI phase measurements using successive calculations. Once the calibration parameters have been determined, the procedure first uses MZI-0, which has the largest FSR, to make a first approximation of the optical frequency and thus determine the interferometer order of MZI-1. The phase measurement of MZI-1 is then used to determine the orders of MZI-2 and MZI-3. With these last two MZI orders known, we finally use the theoretical framework from Section 2 to make a temperature-independent calculation of the optical frequency $\nu $ from the phases of MZI-2 and MZI-3.

This procedure involves calculating relative frequency from the measured phase, as well as calculating the subsequent MZI’s order from this frequency approximation. After first retrieving the relative phase $\mathrm{\delta }{\varphi _i}$ from the photodetector currents [Eq. (5)], the relative frequency $\mathrm{\delta }{\nu _i}$ is then calculated using

$${\mathrm{\delta }{\nu _i} = \frac{{\mathrm{\delta }{\varphi _i} - {C_{i0}} - {C_{i1}}\mathrm{\Delta }T}}{{{C_{i2}} + {C_{i1}}\mathrm{\Delta }T/{\nu _0}}},}$$
which follows from Eq. (15) (see Supplement 1, section C). For these first approximations, assuming $T = {T_0}$ (i.e., $\mathrm{\Delta }T = 0$) is sufficient to determine the subsequent MZI orders. For MZI-0, the order is unambiguous, so $\mathrm{\delta }{\nu _0}$ directly refers to the approximate frequency deviation $\mathrm{\Delta }{\nu _0}$ from the beginning of the range at ${\nu _0}$. For the subsequent MZI calculations, we can use
$${\mathrm{\Delta }{\nu _i} = \mathrm{\delta }{\nu _i} + \mathrm{\Delta }{m_i}\textrm{FS}{\textrm{R}_i},}$$
where $\mathrm{\Delta }{m_i}$ is the difference between ${m_i}$, the order of MZI-i at $\nu $, and ${m_{\textrm{ref},i}}$, the order at ${\nu _0}$, i.e., $\mathrm{\Delta }{m_i} = \; {m_i}\; - \; {m_{\textrm{ref},i}}$. Figure 7 illustrates the relation between the respective orders. The FSR can be calculated in terms of the calibration parameters using Eqs. (17) and (18) as
$${\textrm{FS}{\textrm{R}_i} = \frac{\textrm{c}}{{\mathrm{\Delta }{L_i}({{n_{gi}} + {\theta_i}\mathrm{\Delta }T} )}} = \frac{{2\pi }}{{{C_{\textrm{i}2}} + {C_{i1}}\mathrm{\Delta }T/{\nu _0}}}.}$$
We can then use the calculated frequency approximation to determine the order of the subsequent MZI-(i+1) using
$${\mathrm{\Delta }{m_{i + 1}} = \textrm{round}\left( {\frac{{\mathrm{\Delta }{\nu_i} - \mathrm{\delta }{\nu_{i + 1}}}}{{\textrm{FS}{\textrm{R}_{i + 1}}}}} \right),}$$
noting that in the final step, i=1 determines the orders for both i=2 and i=3.

 figure: Fig. 7.

Fig. 7. Plot of MZI phase measurement, displaying the frequency offset and interferometer order variables. The frequency offset ${\delta }\nu $ can be visualized as the sum of two offsets ${\delta }\nu {^{\prime}}$ and ${\delta }\nu ^{\prime \prime} $.

Download Full Size | PDF

We calculate the unknown laser frequency by iteratively applying Eqs. (20)–(23) to each MZI. The first approximation $\mathrm{\Delta }{\nu _0}$ from MZI-0 determines the order of MZI-1. We then compute a second, more accurate approximation of the frequency using $\mathrm{\Delta }{\nu _1}$, which we use to determine the orders of the final two MZIs with

$${\mathrm{\Delta }{m_{2,3}} = \textrm{round}\left( {\frac{{\mathrm{\Delta }{\nu_1} - \delta {\nu_{2,3}}}}{{\textrm{FS}{\textrm{R}_{2,3}}}}} \right).}$$
With the final two MZI orders known, the relative phases can be related to the phases unwrapped on the calibration range by $\mathrm{\Delta }{\varphi _{2,3}} = \mathrm{\delta }{\varphi _{2,3}} + 2\pi \mathrm{\Delta }{m_{2,3}}$. At this point, we can determine the optical frequency input to the wavemeter using these unwrapped phases measured from the two interferometers with different TOCs. Using the two expressions of $\mathrm{\Delta }{\varphi _{2,3}}$ from Eq. (15), we finally obtain the temperature-independent measurement of the optical frequency:
$${\nu \; = \; {\nu _0} + \Delta \nu = {\nu _0} + \frac{{{C_{21}}({\mathrm{\Delta }{\varphi_3} - {C_{30}}} )- {C_{31}}({\mathrm{\Delta }{\varphi_2} - {C_{20}}} )}}{{{C_{21}}{C_{32}} - {C_{31}}{C_{22}}}}.}$$
This equation is equivalent to Eq. (11) by replacing the C coefficients with their expressions as a function of physical parameters. The result can also be expressed as a wavelength $\lambda = 2\pi c/\nu $, which is the form we will use to display the experimental results.

3.3 Wavelength measurement results

We measure the wavelength of a tunable laser across an 80 nm range using the wavemeter operation procedure described above. In the experimental setup, a Photonetics Tunics external cavity tunable laser’s output is coupled to the wavemeter chip using an optical fiber. The laser is tuned in 800 steps between 1500 nm and 1580 nm, while the eight wavemeter output currents are measured using an Agilent switch unit connected to a Keithley 2000 multimeter. The TEC temperature was treated as unknown, so its value was ignored for the computation of the measured wavelength. The temperature set point (22.5°C) was chosen not to correspond to values chosen during the calibration. Applying the procedure from Section 3.2 to calculate the laser wavelength from the MZI phases, we obtain the wavemeter’s measurement of the wavelength compared to the known value from the laser [Fig. 8(a)], which closely matches a line. The difference between the measured and set wavelengths is plotted in Fig. 8(b) and also corresponds to the accuracy, with the highly accurate (<5 pm repeatability) laser serving as the reference. The deviation from the linear model has a maximum amplitude of 200 pm and mean error (accuracy) of 84.6 pm. However, we can apply a correction to increase the accuracy even further.

 figure: Fig. 8.

Fig. 8. (a) Wavelength measurement based on linear model compared to the known input laser wavelength. (b) Deviation of the measurement from the known laser wavelength.

Download Full Size | PDF

By correcting for reproducible errors in the calibration, the wavemeter accuracy improves beyond that of the initial linear model. Figure 8(b) shows that the linear model’s error has significant periodic components. This error includes detector imbalances and other fabrication defects such as spurious reflections (for example, the period of the oscillations in Fig. 8(b) corresponds to the length of the 90-degree hybrid), which can all be systematically corrected from the measurement because they are predictable. To do so, we store the measurement error of each MZI across the wavelength range at several temperatures. During wavemeter operation, this stored error can be subtracted from calculations made using the linear model, i.e., Eq. (20). Initially, an arbitrary value for the unknown temperature +/- 10°C may be used in the Section 3.2 procedure, producing a reasonably accurate estimation of the wavelength and temperature. The temperature itself can be obtained with the formula

$${T = \frac{{\mathrm{\Delta }{\varphi _3} - {C_{30}} - {C_{32}}\mathrm{\Delta }\nu }}{{{C_{31}} + {C_{33}}\mathrm{\Delta }\nu }} + {T_0},}$$
which follows from Eq. (15). Then, a second iteration of calculations is performed using the obtained temperature and subtracting the stored corrections (applying an interpolation between the values if necessary), producing an accurate wavelength measurement. Figure 9 shows the wavemeter characterization when using the correction for reproducible errors. The measurement error improves to 29 pm amplitude and 10.8 pm RMS accuracy (1.3 GHz) across the 80 nm range, reaching an accuracy near our original target of 10 pm. We can see in Fig. 9(b) that the error still has some periodic components, which suggests that we have not reached the fundamental accuracy of the wavemeter limited by random electrical noise calculated in Eq. (12).

 figure: Fig. 9.

Fig. 9. (a) Wavelength measurement using additional corrections compared to known input laser wavelength. (b) Measurement error (accuracy) with the corrections applied across the wavelength range.

Download Full Size | PDF

The wavemeter accuracy is consistent across a large temperature range, when subjected to rapidly changing temperatures, and for both high and low input power. We measure the mean error across three temperatures unknown to the wavemeter and covering over a 20°C range (Fig. 10). The RMS accuracy stays at 10–12 pm at both ends of the temperature range tested. We also vary the input optical power to the wavemeter from −11.5 dBm (71 µW) to 18 dBm (63 mW) and again find consistent accuracy of 11.0 pm on average. We next monitor the dynamic response of the wavemeter when subjected to a rapidly shifting temperature. To perform a fast acquisition, the eight photocurrents were measured simultaneously at 8.3 kS/s using a National Instruments PXIe 6238 eight-channel acquisition module. With the input laser set to a fixed wavelength, the TEC under the wavemeter chip is switched from 21°C to 39°C (Fig. 11). During this temperature shift, which takes less than five seconds, the wavemeter continues to measure the wavelength and the temperature. Figure 11(b-c) show the wavelength measurement error has negligible fluctuations even as the temperature drifts at ∼5°C/s. The mean error during the temperature ramp at this wavelength is just 4.1 pm, while after the temperature settles, the error is 1.0 pm. This robustness to environmental perturbations demonstrates the benefit of an athermal integrated wavemeter.

 figure: Fig. 10.

Fig. 10. Wavemeter measurement accuracy for various temperatures and input powers.

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. (a) Temperature calculated by the wavemeter during a large fluctuation. (b) Measured wavelength during the temperature shift with the input laser set at 1550 nm. (c) Wavelength measurement error during the temperature shift.

Download Full Size | PDF

4. Discussion and summary

We expect that the wavemeter’s capabilities may exceed the current demonstration. The measurement range was limited by the available equipment and should in fact cover the full FSR of MZI-0, which is 94 nm. Additionally, improved assembly with high-speed analog-to-digital convertors (ADCs) and transimpedance amplifiers (TIAs) should enable characterization of the ultimate measurement rate. We expect nanosecond-scale measurements are attainable with moderate input power. The germanium photodiodes are capable of GHz operation and thus should not limit the speed in this regime. Next, the wavemeter’s operational temperature range and measurement accuracy may be optimized for a particular application by properly adjusting the MZI dimensions, i.e., the FSRs. Thanks to the chip-scale integration of waveguides, much longer delay lengths may be attained in a compact area with low loss, leading to improved measurement accuracy, while adjusting the shortest delay lengths allows control of the wavelength range. In Supplement 1, section D, we derive equations explaining that there is currently a trade-off between the operational temperature range (∼20°C in our case) and the accuracy, ultimately finding that only a rough estimate of the temperature is needed to overcome the trade-off and extend the maximum temperature deviation significantly. Finally, with the corrected phase errors’ magnitude now known, a larger ratio between the FSRs would allow for either a larger wavelength range or elimination of the intermediate stage without sacrificing accuracy. Nevertheless, in the current design we have achieved a favorable combination of accuracy, wavelength range, and temperature-independence.

We have proposed and demonstrated the first athermal integrated wavemeter using silicon photonics. By including interferometers composed of two different materials, our design removes the temperature sensitivity present in previous integrated wavemeters. We show that our wavemeter chip determines the wavelength of a laser with an average accuracy of 11 pm over an 80 nm range. This represents the lowest measurement error across such a broad wavelength range using silicon photonics, to the best of our knowledge. Such accuracy is maintained for temperature changes of at least 20°C, including a rapidly fluctuating temperature. In comparison, without the athermal design, the longest silicon MZI would cause a measurement drift of ∼1,200 pm over 20°C. This design enables low-cost production of robust, integrated wavemeters for an array of applications. The broadband and accurate measurements will be useful in future integrated sensors in combination with a scanning laser. Additionally, the wavemeter could help to stabilize tunable lasers used in communication networks.

Acknowledgments

The authors thank Michael Eggleston, Po Dong, Argishti Melikyan, Andrea Blanco-Redondo, Mark Earnshaw, Jean-Guy Provost, Alexandre Shen, and Joan Ramirez for productive conversations.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

References

1. T. Shinagawa, M. Nishita, T. Sato, H. Nasu, T. Mukaihara, T. Nomura, and A. Kasukawa, “Detailed Investigation on Reliability of Wavelength-Monitor-Integrated Fixed and Tunable DFB Laser Diode Modules,” J. Lightwave Technol. 23(3), 1126–1136 (2005). [CrossRef]  

2. L. Couturier, I. Nosske, F. Hu, C. Tan, C. Qiao, Y. H. Jiang, P. Chen, and M. Weidemüller, “Laser frequency stabilization using a commercial wavelength meter,” Rev. Sci. Instrum. 89(4), 043103 (2018). [CrossRef]  

3. T. Udem, R. Holzwarth, and T. W. Hänsch, “Optical frequency metrology,” Nature 416(6877), 233–237 (2002). [CrossRef]  

4. M. Dobosz and M. Kożuchowski, “Overview of the laser-wavelength measurement methods,” Optics and Lasers in Engineering 98, 107–117 (2017). [CrossRef]  

5. F. V. Kowalski, R. T. Hawkins, and A. L. Schawlow, “Digital wavemeter for cw lasers,” J. Opt. Soc. Am. 66(9), 965–966 (1976). [CrossRef]  

6. J. Ishikawa, N. Ito, and K. Tanaka, “Accurate wavelength meter for cw lasers,” Appl. Opt. 25(5), 639–643 (1986). [CrossRef]  

7. P. J. Fox, R. E. Scholten, M. R. Walkiewicz, and R. E. Drullinger, “A reliable, compact, and low-cost Michelson wavemeter for laser wavelength measurement,” Am. J. Phys. 67(7), 624–630 (1999). [CrossRef]  

8. Bristol Instruments, “828 Series High-Speed Optical Wavelength Meter,” Bristol Instruments (n.d.).

9. D. F. Gray, K. A. Smith, and F. B. Dunning, “Simple compact Fizeau wavemeter,” Appl. Opt. 25(8), 1339–1343 (1986). [CrossRef]  

10. T. B. Jones, N. Otterstrom, J. Jackson, J. Archibald, and D. S. Durfee, “Laser wavelength metrology with color sensor chips,” Opt. Express 23(25), 32471–32480 (2015). [CrossRef]  

11. N. H. Wan, F. Meng, T. Schröder, R.-J. Shiue, E. H. Chen, and D. Englund, “High-resolution optical spectroscopy using multimode interference in a compact tapered fibre,” Nat Commun 6(1), 7762 (2015). [CrossRef]  

12. N. Han, G. N. West, A. H. Atabaki, D. Burghoff, and R. J. Ram, “Compact and high-precision wavemeters using the Talbot effect and signal processing,” Opt. Lett. 44(17), 4187–4190 (2019). [CrossRef]  

13. G. D. Bruce, L. O’Donnell, M. Chen, M. Facchin, and K. Dholakia, “Femtometer-resolved simultaneous measurement of multiple laser wavelengths in a speckle wavemeter,” Opt. Lett. 45(7), 1926–1929 (2020). [CrossRef]  

14. J. N. Porter, J. S. Jackson, J. S. Jackson, D. S. Durfee, D. S. Durfee, and R. L. Sandberg, “Laser wavelength metrology with low-finesse etalons and Bayer filters,” Opt. Express 28(25), 37788–37797 (2020). [CrossRef]  

15. M. Muneeb, A. Ruocco, A. Malik, S. Pathak, E. Ryckeboer, D. Sanchez, L. Cerutti, J. B. Rodriguez, E. Tournié, W. Bogaerts, M. K. Smit, and G. Roelkens, “Silicon-on-insulator shortwave infrared wavelength meter with integrated photodiodes for on-chip laser monitoring,” Opt. Express 22(22), 27300–27308 (2014). [CrossRef]  

16. A. Ruocco and W. Bogaerts, “Fully integrated SOI wavelength meter based on phase shift technique,” in 2015 IEEE 12th International Conference on Group IV Photonics (GFP) (2015), pp. 131–132.

17. E. Martín-López and D. Bitauld, “Increasing the Range and Precision of Integrated Wavemeters,” in Conference on Lasers and Electro-Optics (2019), Paper AF2 K.3 (Optical Society of America, 2019), p. AF2 K.3.

18. C. Xiang, M. A. Tran, T. Komljenovic, J. Hulme, M. Davenport, D. Baney, B. Szafraniec, and J. E. Bowers, “Integrated chip-scale Si3N4 wavemeter with narrow free spectral range and high stability,” Opt. Lett. 41(14), 3309–3312 (2016). [CrossRef]  

19. L. Chen, C. Doerr, S. Liu, L. Chen, and M. Xu, “Silicon-based Integrated Broadband Wavelength-meter with Low Temperature Sensitivity,” in Optical Fiber Communication Conference (OFC) 2020 (2020), Paper M1C.3 (Optical Society of America, 2020), p. M1C.3.

20. B. Guha, J. Cardenas, and M. Lipson, “Athermal silicon microring resonators with titanium oxide cladding,” Opt. Express 21(22), 26557–26563 (2013). [CrossRef]  

21. H. Podmore, A. Scott, P. Cheben, C. Sioris, P. Cameron, J. H. Schmid, A. Lohmann, Z. Corriveau, and R. Lee, “Athermal planar-waveguide Fourier-transform spectrometer for methane detection,” Opt. Express 25(26), 33018–33028 (2017). [CrossRef]  

22. D. A. Ackerman, K. M. Paget, L. F. Schneemeyer, L. J.-P. Ketelsen, F. W. Warning, O. Sjolund, J. E. Graebner, A. Kanan, V. R. Raju, L. E. Eng, E. D. Schaeffer, and P. V. Emmerik, “Low-Cost Athermal Wavelength-Locker Integrated in a Temperature-Tuned Single-Frequency Laser Package,” J. Lightwave Technol. 22(1), 166–171 (2004). [CrossRef]  

23. M. Uenuma and T. Motooka, “Temperature-independent silicon waveguide optical filter,” Opt. Lett. 34(5), 599–601 (2009). [CrossRef]  

24. A. Melikyan, G. de Valicourt, K. Kim, N. Fontaine, Y.-K. Chen, and P. Dong, “Hybrid III-V/Silicon Laser with Integrated Athermal Wavelength Locker,” in 2017 European Conference on Optical Communication (ECOC) (2017), pp. 1–3.

25. J. M. Ramirez, H. Elfaiki, T. Verolet, C. Besancon, A. Gallet, D. Néel, K. Hassan, S. Olivier, C. Jany, S. Malhouitre, K. Gradkowski, P. E. Morrissey, P. O’Brien, C. Caillaud, N. Vaissière, J. Decobert, S. Lei, R. Enright, A. Shen, and M. Achouche, “III-V-on-Silicon Integration: From Hybrid Devices to Heterogeneous Photonic Integrated Circuits,” IEEE J. Select. Topics Quantum Elec. 26, 6100213 (2020). [CrossRef]  

26. A. Hänsel and M. J. R. Heck, “Widely Tunable Laser on an InP Photonic Integrated Circuit,” IEEE Photonics Technol. Lett. 32(2), 105–108 (2020). [CrossRef]  

27. B. Stern, K. Kim, C. Calò, C. Fortin, D. Lanteri, K. Mekhazni, J.-G. Provost, A. Wilk, F. Mallécot, F. Pommereau, A. Melikyan, and P. Dong, “Broadly and finely tunable hybrid silicon laser with nanosecond-scale switching speed,” Opt. Lett. 45(22), 6198–6201 (2020). [CrossRef]  

28. D. Bitauld, “Wavelength meter,” World Intellectual Property Organization patent WO2020083626A1 (30 April 2020).

29. P. Dong, X. Liu, S. Chandrasekhar, L. L. Buhl, R. Aroca, and Y. K. Chen, “Monolithic Silicon Photonic Integrated Circuits for Compact 100+ Gb/s Coherent Optical Receivers and Transmitters,” IEEE J. Sel. Top. Quantum Electron. 20(4), 150–157 (2014). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Supplement 1

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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 (11)

Fig. 1.
Fig. 1. MZI phase and wavelength measurement. (a) Diagram of an integrated MZI for phase measurement. (b) Spectrum of a first MZI with a small $\mathrm{\Delta }{L_0}$ and large FSR. The red dots indicate the wavelength approximated by the photodetector current measurements ${X_0}$ and ${Y_0}$. (c) Relative phase inferred from the first MZI’s photocurrents. (d) Spectrum of a second MZI with a large $\mathrm{\Delta }{L_1}$ and small FSR. (e) Relative phase inferred from the second MZI’s photocurrents and from the order m as determined by the phase of the first MZI.
Fig. 2.
Fig. 2. Diagram of athermal wavemeter, composed of MZIs with different TOCs.
Fig. 3.
Fig. 3. (a) Diagram of athermal wavemeter design. (b) Photograph of assembled device. (c) Microscope image of the wavemeter chip.
Fig. 4.
Fig. 4. Photodetector current measurements from each pair of balanced photodiodes (X and Y) for (a) MZI-0 (b) MZI-1 (c) MZI-2 and (d) MZI-3.
Fig. 5.
Fig. 5. (a-d) Measured unwrapped phases $\; \Delta \varphi $ for MZI-0 to MZI-3. (e-h) Measurement deviations from the linear fits of MZI-0 to MZI-3 phases.
Fig. 6.
Fig. 6. (a-d) Measurement and linear fit of A parameter versus temperature for MZI-0 to MZI-3. (e-h) Measurement and linear fit of B parameter across temperature.
Fig. 7.
Fig. 7. Plot of MZI phase measurement, displaying the frequency offset and interferometer order variables. The frequency offset ${\delta }\nu $ can be visualized as the sum of two offsets ${\delta }\nu {^{\prime}}$ and ${\delta }\nu ^{\prime \prime} $.
Fig. 8.
Fig. 8. (a) Wavelength measurement based on linear model compared to the known input laser wavelength. (b) Deviation of the measurement from the known laser wavelength.
Fig. 9.
Fig. 9. (a) Wavelength measurement using additional corrections compared to known input laser wavelength. (b) Measurement error (accuracy) with the corrections applied across the wavelength range.
Fig. 10.
Fig. 10. Wavemeter measurement accuracy for various temperatures and input powers.
Fig. 11.
Fig. 11. (a) Temperature calculated by the wavemeter during a large fluctuation. (b) Measured wavelength during the temperature shift with the input laser set at 1550 nm. (c) Wavelength measurement error during the temperature shift.

Tables (1)

Tables Icon

Table 1. Mach-Zehnder interferometer parameters

Equations (26)

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

φ = Δ L n eff ( λ , T ) 2 π λ ,
n eff ( λ , T ) n eff ( λ 0 , T 0 ) + ( T T 0 ) n eff T | λ 0 , T 0 + ( λ λ 0 ) n eff λ | λ 0 , T 0 .
φ = Δ L 2 π λ ( a + b Δ λ + θ Δ T ) .
φ = Δ L 2 π λ ( n g + θ Δ T ) + 2 π b Δ L .
δ φ = atan 2 ( Y , X ) .
φ = δ φ + 2 π m π / 4 ,
λ = Δ L n g + θ Δ T δ φ / 2 π + m 1 / 8 b Δ L .
φ 2 = Δ L 2 2 π ν c ( n g 2 + θ 2 Δ T ) + 2 π b 2 Δ L 2
φ 3 = Δ L 3 2 π ν c ( n g 3 + θ 3 Δ T ) + 2 π b 3 Δ L 3 ,
J J φ ( ν , T ) ν 0 , T 0 = 2 π c ( Δ L 2 n g 2 Δ L 2 θ 2 ν 0 Δ L 3 n g 3 Δ L 3 θ 3 ν 0 ) .
ν = c 2 π Δ L 2 Δ L 3 θ 2 Δ L 2 φ 3 θ 3 Δ L 3 φ 2 2 π Δ L 2 Δ L 3 ( b 3 θ 2 b 2 θ 3 ) θ 2 n g 3 θ 3 n g 2 .
σ ν = c 2 π Δ L 2 Δ L 3 ( Δ L 2 θ 2 ) 2 + ( Δ L 3 θ 3 ) 2 θ 2 n g 3 θ 3 n g 2 σ φ ,
φ = Δ φ + 2 π m ref π / 4 ,
Δ φ i ( Δ ν , Δ T ) = A i ( Δ T ) + B i ( Δ T ) Δ ν ,
Δ φ i ( Δ ν , Δ T ) = C i 0 + C i 1 Δ T + C i 2 Δ ν + C i 3 Δ T Δ ν .
C i 0 = Δ φ ( ν 0 , T 0 ) = 2 π Δ L i c ν 0 n g i + 2 π b i Δ L i 2 π m r e f + π / 4
C i 1 = J i 1 = 2 π Δ L i c ν 0 θ i
C i 2 = J i 2 = 2 π Δ L i c n g i
C i 3 = 2 π Δ L i c θ i = C i 1 ν 0 .
δ ν i = δ φ i C i 0 C i 1 Δ T C i 2 + C i 1 Δ T / ν 0 ,
Δ ν i = δ ν i + Δ m i FS R i ,
FS R i = c Δ L i ( n g i + θ i Δ T ) = 2 π C i 2 + C i 1 Δ T / ν 0 .
Δ m i + 1 = round ( Δ ν i δ ν i + 1 FS R i + 1 ) ,
Δ m 2 , 3 = round ( Δ ν 1 δ ν 2 , 3 FS R 2 , 3 ) .
ν = ν 0 + Δ ν = ν 0 + C 21 ( Δ φ 3 C 30 ) C 31 ( Δ φ 2 C 20 ) C 21 C 32 C 31 C 22 .
T = Δ φ 3 C 30 C 32 Δ ν C 31 + C 33 Δ ν + T 0 ,
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.