Abstract

We study the effects of laser phase noise on a phase diversity coherent optical frequency domain (C-OFD) technique that has been recently proposed to measure passive devices used in dense wavelength division multiplexing (DWDM) systems. Theoretical expressions are provided to calculate the laser phase-noise to intensity-noise conversion in this technique under simplified circumstances. Obtained simulation results for a realistic measurement set-up show the validity of the approximate expressions. It is concluded that this effect is one of the limiting source of error for this measurement technique.

], with promising possibilities. This method is based in the six-port measurement technique, a well established and accurate method to measure complex reflection coefficients at microwave and millimeter-wave frequencies which is used by several metrological institutes to determine the complex reflection coefficient of commercial working standards [7

]. Compared to other C-OFD techniques, which rely on frequency sweeping and interferogram demodulation, the six-port technique is intrinsically a single frequency technique, i.e. no sweeping of the laser source is necessary to completely determine the modulus and absolute phase of the reflection coefficient of the load at that specific frequency. This makes its precision basically dependent on the calibration process and not on the sweeping characteristics of the laser source. The core of the proposed technique is a small and rugged planar lightwave circuit (PLC) six-port junction, based on multi-mode interference (MMI) couplers, which was designed and simulated showing an excellent performance in a wide wavelength range. Although this germinal work has with no doubt open the possibility of extending the six-port measurement technique, so successful at microwave frequencies, up to the optical range, assessment of the optical six-port measurement technique as a whole, must still be carried out.

It is a well recognized fact that, aside from calibration errors, the accuracy of the six-port technique is determined by the power measurement errors [9

]. Thus assessment of this technique at optical frequencies requires a careful study of the possible factors causing measuring errors at the four power detectors of the six-port setup. In doing so, special attention should be paid to those uncertainty mechanisms which are not encountered on microwave measurement setups and that, consequently, have not yet been quantified.

In this work the optical six-port technique is first briefly reviewed and the possible sources of power measurement errors are summarized and compared with those existing at microwave frequencies. From this analysis, optical coherence issues related to the laser phase noise (i.e., linewidth), appear as one of the main differences with the six-port microwave counterpart. Then the effects of laser phase noise on six-port power measurement errors is studied and simulated, and simplified expressions are obtained to calculate the power measurement uncertainties from the constituent parts of the setup.

2. The six-port measurement technique

Figure 1 shows the schematic of the proposed six-port optical reflectometer for measuring the complex reflection coefficient of a load (also known as device under test or DUT). The figure corresponds to a single polarization case. Of course, in a practical situation, polarization related measurements would require additional hardware (polarization controller at the laser port and polarization diversity detection in each of the power meter outputs), but the study of phase noise effects can be carried out in a simpler manner based on this simplified single polarization assumption. The set-up is comprised of a narrow linewidth tunable laser source (TLS), a passive six-port junction, and four power meters (PM). The core of the measurement set-up is the passive six-port junction. This is a recently proposed [5

], planar lightwave circuit (PLC) based on multimode interference couplers (MMIs) on a Δn=0.75% contrast Silica-on-Silicon buried channel waveguides technology. Although completely described in the aforementioned references, in Fig. 2 a schematic representation of the proposed six-port PLC is presented for the sake of completeness. It is basically composed of three interconnected MMI’s: two identical 2×2 couplers (M1 and M2) and one 2×3 120° coupler (M3). M1 acts as a power divider; M2 is used as a directional coupler, i.e. to separate the incident and reflected waves in the DUT and M3 is the one in charged of setting the 120° phase difference between the three measurement signal at ports 3, 5 and 6, as required in a idealized six-port design [7

]. The circuit schematic also contains two interconnecting waveguides WD and WC which are important in the six-port balance. TLS and PMs are considered to be high performance commercially available equipment. Interconnection between the different subsystems in the measurement set-up is done by short pieces of optical fibre, and low reflectivity connectors, such as straight PC ones (return losses better than -40 dB).

where, for a certain frequency ω, Ki
(ω) and qi
(ω) are real and complex constant, respectively, which only depend on the six-port PLC and the reflection coefficient of the PMs (in real situations in which the PMs are well matched, these constants are almost completely determined by the six-port junction itself), b2
(ω) is the incident wave on the DUT, and Γ
L
(ω)=a2(ω)/b2(ω) is the DUT reflection coefficient to be measured.

The complex values qi
(ω) play a crucial role in six-port theory. Six-port design goal is to get q4−1
(ω)=0 (so that port 4 acts as a power reference port) and q3
(ω), q5
(ω), and q6
(ω) located symmetrically at 120° over a circumference of radius approximately equal to 1.5. These design values have shown to give the best accuracy for passive load measurement [7

], that the proposed optical PLC approximately fulfils these design criteria in a great wavelength span of 100nm.

Normalizing power equations Eq. (1) with respect to the power of the reference port (habitually P4) we come up with the three power ratio equations

pi(ω)=ki(ω)∣1−qi−1(ω)ΓL(ω)∣2∣1−q4−1(ω)ΓL(ω)∣2i=3,5,6

(2)

where pi
(ω)=Pi
(ω)/P4(ω) and ki
(ω)=Ki
(ω)/K4
(ω) which are the basis of six-port technique. For each frequency ω, Eq. (2) define three circles in the complex plane whose intersection is the reflection coefficient of the load Γ
L
(ω). In these equations, there are four complex constants (q3, q4, q5
, and q6
) and three real constants (k3, k5
, and k6
) which depend on the six-port junction and reflection coefficient of the PMs, but not of the load, that are determined by calibration [10

3. Power measurement uncertainties

In practical situations, due to power measurement errors, the three circles defined by Eq. (2) will not intersect in a point, but still six-port theory provides an statistically quasi-optimum value for the unknown reflection coefficient Γ
L
that can be estimated from the erroneous measurements [11

], that these power measurement errors will, aside from calibration errors, determine the accuracy of the measurement technique. Thus, the causes of power measurement errors in the optical six-port set-up should be analyzed.

Main causes of power measurement errors in microwave six-port technique have been found to be related to thermal noise, zero drift, and nonlinearities [9

], all of them directly linked to the PM performance. Possible effects due to the non-zero linewidth of the source have not been accounted for probably because, at microwave frequencies, the electrical size of the involved circuits is small, thus making this effects unimportant. At optical frequencies this assert no longer holds: the circuit size is typically several thousands of wavelengths long and finite coherence effects must be analyzed.

It must be noticed that, power measurement errors due to nonlinearity of the PMs are expected to be better at optical frequencies than at microwaves. Pin diodes, typically used as optical power sensing devices, are highly linear devices exhibiting dynamic ranges in excess of 60 dB, that are greater than their microwave counterparts, usually limited to approximately 50 dB (this is, for example, the case of Schottky diodes or thermocouples). In fact, it is usually recognized that optical PM nonlinearity is mainly due to the electronic circuitry that follows the photodetector, and not to the pin diode itself [1

]. Of course this kind of nonlinearity is common to both optical and microwave setups.

On the other hand, power measurement noise sources can be divided into photodetector shot (quantum) noise and electronic noise. Electronic noise sources are thermal and low frequency flicker noise and are linked to the amplifier stages following the photodetector, thus being also common to microwave and optical setups. This noise levels do not depend on the received optical power and thus establish a noise floor at low received powers. On the contrary, shot noise, which increases with received power will be the limiting noise source at high optical levels. The effects of all these noise contributions can be easily calculated through well known expressions.

Besides this easily recognizable noise sources, there still exists another effect that is not so evident: it is the one that could follow from laser phase noise (related to the non-zero linewidth of the laser or, in other words, to finite coherence of the source) to intensity noise conversion. This is a well known effect that can be troublesome in C-OFDM techniques [1

], and that is related to multiple reflections that cause different delayed optical signals to interfere with each other at the photoreceiver. The sixport technique being a phase diversity C-OFD technique in which each of the PMs (photodetectors) receives a combination, with different phases, of the two monochromatic waves to be compared (incident and reflected waves at the DUT) will with no doubt suffer from this detrimental effect.

4.1 Frequency domain interpretation of coherent homodyne detection

The influence of laser phase noise in coherent detection schemes, in which two optical beams interfere on a photodiode, has been analyzed in various situations including heterodyne and homodyne detection with different degrees of correlation between the interfering fields [13

]. In that work, closed expressions were given for the power spectral density of the received photocurrent under the following approximations: a Lorentzian laser lineshape was assumed and only two incident waves, with certain relative level of delay (τ0) and attenuation (α), were supposed to be present at the photodiode. Although it seems that no relation can be established between this simplified studies and the situation in the optical six-port, it will be shown that a close linkage exists among them. It must be realized that, phase noise to intensity noise conversion will happen whenever an optical noisy carrier passes through a linear system with frequency dependant amplitude response. In fact, from this point of view, the homodyne coherent detection scheme described in [13

], in which a field and its delayed (τ0) and attenuated (α) replica simultaneously impinge on a photodiode (with responsivity R), can be represented with the help of Fig. 3.

Fig. 3. Equivalent model for the homodyne coherent detection scheme

In this figure the laser output, which is modeled as a quasi-monochromatic amplitude-stabilized wave with random phase fluctuations

A(t)=A0·ej[ω0t+ϕ(t)]

(3)

is filtered by a linear time invariant (LTI) system with frequency response

H(ω)=TF[h(t)]=TF[δ(t)+αδ(t−τ0)]=1+αe−jτ0ω

(4)

which represents the effect of interfering the two optical beams, a direct one and a delayed and attenuated replica, giving up to the interference field at the input of the photoreceiver

B(t)=h(t)*A(t)=A(t)+αA(t−τ0).

(5)

The laser output spectrum is of Lorentzian type, thus its power spectral density is given by

SAA(ω)=A022πΔν(πΔν)2+(ω−ω0)2

(6)

where Δν is the full linewidth at half maximum (FWHM), ω0 is the central angular frequency and A0 is the amplitude of the laser source. Elementary LTI system theory dictates that the spectrum of the wave entering the photodiode will be given by

SB(ω)=SAA(ω)∣H(ω)∣2

(7)

The two sided spectrum of the detected photocurrent I(t) was calculated to be [13

This filtering effect of the linear time invariant system is represented in Fig. 4 for a typical situation.

It can be seen that the ripple in the frequency response causes the spectrum of the optical field entering the photodiode to be distorted and this, once power detected, is converted into intensity noise. From this analysis it is clear that if the frequency response of the six-port setup (from the input laser port to any of the PM output ports) is somehow approximated, inside the spectral width of the laser source, in a way similar to that in Eq. (4), then the results in [13

]) for the designed six-port PLC 5,6 show that internal reflections between the different discontinuities inside the circuit are very small and thus the following approximations can be done: a) S22
(ω) is always below -60dB ; b) qi
(ω) are slowly varying with ω; c) the modulus of Si1
(ω) is also a slowly varying function of ω, and its phase can be locally approximated by a line. Thus when applying (9) with a1
(ω) being the laser optical output centered at ω0 and with a linewidth Δν (typically below 10MHz) the following approximations hold

In Fig. 5 simulated results for these constants are presented showing that the aforementioned approximations are valid: Fig. 5(a) shows the S2,2 modulus, Fig. 5(b) shows the evolution of the qi
points in the complex plane in a total frequency span of 125 GHz (1 nm) and Figs. 5(c) and 5(d) show, respectively, the modulus and phase of S3,1 in the same frequency span. Insets are included in Figs. 5(c) and 5(d) showing the detailed behaviour in a reduced frequency span of ±10MHz correspondent to the maximum expected value of the laser linewidth. Although not plotted, for the sake of clarity, the same behaviour has been observed for the other Si1 parameters.

Putting these approximations into (9) the reduced expression

bi(ω)=Si1(ω0)e−jωτi(ω0)[1−qi−1(ω0)ΓL(ω)]a1(ω)i=3,4,5,6

(12)

is obtained. The frequency response Hi
(ω) from input port 1 to output port i (i=3,4,5,6) of the six-port can then be identified as

Hi(ω)=Si1(ω0)e−jωτi(ω0)[1−qi−1(ω0)ΓL(ω)]i=3,4,5,6

(13)

Fig. 5. Simulated frequency response of the optical six-port PLC. Simulation results have been calculated in 4000 frequency points in 1nm bandwidth centred at the laser central wavelength of 1550nm. The results shown in the insets of Fig. 5(c) and 5(d) have been calculated with a finer mesh of 4000 points in a ±10 MHz bandwidth around 1550nm.

It shows that, for the laser linewidth frequency scale, the effect of the idealized six-port junction is only to introduce a constant delay and attenuation, and thus the possible frequency domain ripple will only be due to the reflection coefficient of the DUT.

If the DUT is selected to be a perfect reflector at the end of a piece of idealized fiber (offset short), the reflection coefficient of the DUT will be

ΓL(ω)=−e−j2τLω

(14)

where τL is the one-way delay of the piece of fiber (related to its length by τ
L
=(neff LDUT)/c where neff is the equivalent refractive index of the fiber and LDUT the length of the fiber). Putting Eq. (14) into Eq. (13) yields

Hi(ω)=Si1(ω0)e−jωτi(ω0)[1−qi−1(ω0)e−j2τLω]i=3,4,5,6

(15)

It must be noticed that Eq. (15) is identical to Eq. (4) when α=q−i1
(ω0) and τ0=2τL (except for α being complex and for the constant scaling (Si1(ω0)) and delay (τi(ω0)) of the preceding terms in (15)). Thus, under these approximations, Eq. (8) can be used to calculate the output photocurrent spectrum at the PMs of the six-port setup S
II,,i
(ω) yielding

where τ¯L
=Δντ
L
is the DUT-delay to laser-coherence time ratio, θi is now given by θi=ω02τ
L
-∠qi
, and the rest of parameters are defined as before. The first term in Eq. (16) is the squared DC photocurrent due to the beating between incident and reflected waves at the DUT (we will refer to it as SII_DC=IDC2, with IDC being the DC photocurrent). It contains a exponential contrast loss term due to the finite coherence of the laser. The second term in Eq. (16) accounts for the phase noise to intensity noise conversion (it will we referred as SII_NOISE).

In practical situations τ¯L
≪1 is a must for reduced contrast loss. Furthermore, the bandwidth of the electronic amplifiers, following the PM photodiodes, is sufficiently small so that in its pass-band the approximation ω¯
≪1 holds. Then the intensity noise spectrum can be considered flat and Eq. (16) can be calculated as

where the second term gives the photocurrent intensity-noise power density. Thus the relative intensity noise (RINi), due to the laser phase-noise to intensity-noise conversion at output port ‘i’, can be calculated as

In the designed PLC |q-1i(ω0)| is almost cero for the reference port 4, while it can be approximated by 0.6 for ports 3, 5 and 6 (see Fig 5(b)) in the analysed frequency span of 1nm. Thus the reference port will not suffer from phase noise conversion and, for the remaining ports, Eq. (18) shows that the RIN only depends on the DUT length through τL and θ
i
.

5. Simulation results

In the previous section approximate analytical expressions where obtained to calculate the phase noise induced RIN in a PLC based six-port setup under idealized circumstances, i.e. perfectly matched PMs. In a real situation, PMs will not be perfectly matched and beating of incident and reflected waves on the interconnecting fibers will cause rapidly changing ripple in the frequency response which can increase the noise level.

To test the validity of the aforementioned approximations in a realistic situation a simulator has been developed to calculate the photocurrent power spectral density at the power detectors for the complete system of Fig. 3. The following procedure is applied: i) a realization of the laser source stochastic process is generated in the time domain as a quasi-monochromatic constant amplitude wave Eq. (3) with band-limited Gaussian random phase fluctuations corresponding to a Lorentzian spectrum (Eq. (6)) with spectral linewidth Δν; ii) this realization is then Fourier transformed to the frequency domain and filtered by the frequency responses Hi
(ω) (i=3,4,5,6) of the complete six-port measurement set-up (including the possible beating due to incident and reflected waves in the interconnecting fibers and the response of the DUT) to get a realization of the wave entering the photodetectors in the frequency domain; iii) this is again transformed to the time domain and power detected to get the final photocurrent diode output whose power spectral density is then calculated. This process is repeated several times (Nit=300) and averaged to reduce the variance of the estimator and improve its quality.

The frequency responses of the system, from the laser input port to any of the PMs output Hi
(ω), is calculated from the scattering parameters of the six-port PLC, those of the interconnecting devices (fibers, connectors…) and the frequency response of the DUT ΓL(ω) by simple network analysis. Six-port junction PLC scattering parameters are obtained, as previously stated, by means of a Bidirectional 3D-Scalar Method of Lines Beam Propagation Method [14

] which accounts for all the reflections in the different discontinuities of the PLC. The interconnection between the PLC and the power detectors have been modeled by pieces of standard fiber (neff=1.5) of the same length (Lf=1m). Constant reflectivity (Rp=-40dB) connectors are inserted between the PLC and the fibers, simulating the pigtailing process. Also the PMs are simulated as constant reflectivity (Rpm=-40dB) blocks followed by a detector with unit responsivity (R=1A/W). A ideal reflector at the end of a piece of standard fiber, of length LDUT, has been considered as the DUT for all the simulations. The laser source wavelength, linewidth, and power have been set to λ0=1550 nm, Δν=100 KHz and A0=0 dBm, respectively (corresponding to a high quality commercial tunable laser source).

Figure 6 shows, for two different DUT lengths, a comparison of the detected photocurrent noise power spectral density at port 5 (SII_NOISE,5) obtained by simulation and evaluated by the analytical expressions obtained in previous section: the solid red curves correspond to analytical approximation (Eq. (16)) obtained when no reflections at the PM/PLC interconnection occur; the green dashed lines correspond to its narrow-band approximation (Eq. (17)); the dotted blue line are the simulation results obtained, as described in previous paragraph, for the complete system including the effects of the reflections in the PM/PLC interconnection. In this figure an excellent agreement between the theoretical curves, obtained in the previous section, and the simulation results is observed. Insets in Fig. 6 correspond to a narrower frequency span of 10 MHz (much greater than expected bandwidth of amplifier stages following the photodetectors) showing the validity of the narrow-band (flat spectrum) approximation Eq. (17). Similar results where obtained for the other measurement ports but they are not plotted for the sake of brevity.

Once verified that the noise term of the photocurrent spectrum of Eq. (17) agrees with the simulation results, we should also confirm that the same its true for the DC part of this equation. Table 1 shows the comparison between simulated and theoretically predicted DC photocurrent (Eq. (17)) for measurement ports 3, 5 and 6, and for two different DUT lengths (LDUT=0.2 m and LDUT=1 m), showing also an excellent agreement.

Table 1:. Theoretical (eq. (16)) and simulated results of the DC photocurrent in each measurement port.

From Fig. 6 and Table 1, it is clear that Eq. (17) accurately predicts the complete photocurrent spectrum for different situations. Thus, Eq. (18) can be used to calculate the laser phase-noise induced RIN for the complete six-port measurement setup. This is done in Fig. 7 where phase-noise induced RIN (Eq. (18)) versus detected DC photocurrent (DC term of Eq. (17)) is plotted and compared with the shot and thermal noise RIN. These results correspond to port 5 and different DUT lengths, but similar results have been obtained for ports 3 and 6.

From Fig. 7 it can be seen that for the analyzed example, laser phase-noise to intensity-noise conversion is the main source of power measurement errors, greatly exceeding the expected detrimental effects of other sources such as thermal or shot noise.

6. Conclusions

The influence of laser phase-noise to intensity-noise conversion has been studied in a recently proposed phase diversity C-OFD six-port measurement technique. Closed expressions have been developed to calculate the laser phase noise induced RIN at the four measurement PMs, based on simplified approximations. A stochastic simulator has been developed to simulate the complete system under realistic conditions. Simulation results show the validity of the aforementioned expressions. This study clearly shows that the laser phase-noise to intensity-noise conversion is the main source of power meter noise and thus this effect cannot be neglected when studying the performance of the six-port measurement technique at optical frequencies. Further studies must be carried out to quantify the consequences of this effect in the six-port calibration procedure and in the global accuracy of the measurement.

Table 1:. Theoretical (eq. (16)) and simulated results of the DC photocurrent in each measurement port.

Acknowledgments

This work was supported by the Spanish MCYT under project TIC2003-07860. The authors wish to thanks Dr. Wangüemert-Pérez and Mr. Ortega-Moñux for their invaluable help in the development of the algorithms to simulate the PLC.

Other

Cited By

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.