Advanced
Alternative Method of AWG Phase Measurement Based on Fitting Interference Intensity
Alternative Method of AWG Phase Measurement Based on Fitting Interference Intensity
Journal of the Optical Society of Korea. 2012. Jun, 16(2): 91-94
Copyright ©2012, Optical Society of Korea
  • Received : April 04, 2012
  • Accepted : June 06, 2012
  • Published : June 25, 2012
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Yong Ho Oh
Sungwoo Lim
Chun Soo Go
kosoo@wku.ac.kr
Abstract
Arrayed waveguide grating (AWG) phase errors are normally assessed from the Fourier transform of the interference intensity data in the frequency domain method. However it is possible to identify the phases directly from the intensity data if one adopts a trial-and-error method. Since the functional form of the intensity profile is known, the intensities can be calculated theoretically by assuming arbitrary phase errors. Then we decide the phases that give the best fit to the experimental data. We verified this method by a simulation. We calculated the intensities for an artificial AWG which is given arbitrary phases and amplitudes. Then we extracted the phases and amplitudes from the intensity data by using our trial-and-error method. The extracted values are in good agreement with the originally given values. This approach yields better results than the analysis using Fourier transforms.
Keywords
I. INTRODUCTION
Low-crosstalk arrayed waveguide gratings (AWGs) are key components in dense wavelength-division-multiplexing (DWDM) photonic networks [1] . The crosstalk depends on the phase errors of the waveguide array and can be reduced by various compensation techniques requiring precise phase measurements [2 - 4] .
The phases of an AWG can be measured by an optical low coherence method [5 - 7] or by a frequency domain method [8 - 10] . The latter method is easier to implement because it does not require fine control over the path length in an interferometer. In this method, the phases are obtained from the Fourier transform of the interference intensities, which appear as a series of peaks that have the phase information of each waveguide. However, it is impossible to completely eliminate overlaps of neighboring peaks. Hence, another process is required to extract the phase values. The measurement accuracy may also be lowered by the error in calculating the Fourier transform.
In this paper, we propose a method to obtain phase values directly from the interference intensities, without the need for Fourier transformation. We verify this approach by a simulation in which an artificial AWG is analyzed. In section 2, we describe the method of phase measurement. The results of the simulation are reported in section 3.
II. PHASE MEASUREMENT
Figure 1 shows a schematic diagram of the experimental setup for the frequency domain method. The AWG is included in an arm of a fiber interferometer, and the interference intensity is measured as a function of the laser frequency.
When the AWG is composed of N-waveguides, the interferogram, which is defined as the normalized AC component of the interference intensity, can be expressed as follows [11] :
Lager Image
Here, Ak is the amplitude of the transmittance of the k -th waveguide, Lk is the length of the path containing the k -th waveguide including the effect of the refractive index, Lr is the path length of the reference arm, and σ is the
Lager Image
Schematic diagram of the experimental setup for frequency domain method.
reciprocal of the laser wavelength. Lk can be expressed as Lk = Lo + kΔL + ϵk , where Lo is the length of the path containing the shortest waveguide ( k = 0), ΔL is the step increment, and ϵk is the error in the step increment.
We derive the values of the amplitudes and lengths by fitting the theoretical interferogram to the experimental data through the Monte-Carlo method. The process is as follows. First, we assign initial values to the set of parameters { Ak , Lk }. Either uniform, Gaussian or sinc distributions can be assigned initially to { Ak }. The initial values of { Lk } may be given randomly, however the computational time can be reduced by starting with the ideal lengths of the AWG. Using the initial values, we calculate the cost function, which is defined as the sum of the squared differences between the theoretical and the experimental interferograms over σ . Next, one of the parameters is chosen randomly and changed slightly. If the cost function increases, the value of the changed parameter is restored. If the cost function decreases, the change is accepted. This process is repeated until the cost function converges. We wrote a program for the process of automatic optimization. The last set { Ak , Lk } includes the desired amplitudes and phases.
III. SIMULATION
We simulated the phase measurement process to verify our method. First, we created an artificial AWG with arbitrary phases and amplitudes. Then, we calculated the virtual interferogram by using Eq. (1). Considering this as an experimental interferogram, we derived the phases and compared them with the initially given values.
The number of waveguides in the array ( N ) was assumed to be 100. A random length error ( ϵk ) within a wavelength was given initially to each waveguide. Figure 2 shows the given errors, which were converted to the phase values for the central wavelength. The amplitudes were given as either Gaussian, sinc or uniform distribution. The laser wavelength of the interferometer was assumed to be scanned from 1500 nm to 1600 nm. Then, the value of Δ σ was equal to 0.042 μm -1 ). The path difference Lo - Lr was fixed at 10 mm
Figure 3 shows an example of a virtual interferogram for an AWG whose step increment ( Δ L ) is equal to 50 μm.
Lager Image
Given phase error distribution for the waveguide array in the simulation.
Lager Image
Example of virtual interferogram. The step increment of the AWG is 50 μm. The symbols in the inset indicate the results of the theoretical fit.
The symbols in the inset represent the theoretical data fitted by the Monte-Carlo method. Here, 10000 data points for σ are used to calculate the cost function. The computer running time for the Monte-Carlo-optimization is approximately 5 minute. We can see in the figure that the intensity oscillates with σ . The period depends on the value of Lo - Lr . The interval between the two peaks corresponds to the free spectral range (FSR), i.e., the reciprocal of ΔL .
Previous studies have shown that the possibility of phase measurement depends on the laser tuning range ( ΔL ) of the interferometer. Hence, we measured the phases for various ratios of the tuning range to FSR (Δσ/FSR). The results are shown in Fig. 4 . The vertical axis represents the phase measurement error, which is the difference between measured and the given lengths of each waveguide. The numbers in the box indicate the values of Δσ /FSR. As we can see in the figure, the measurement errors are nearly zero for all the waveguides if the ratio is equal to or
Lager Image
Phase measurement errors of the waveguide array. The numbers in the box are the ratios of the frequency tuning range to the FSR.
Lager Image
Standard deviations of phase errors as a function of the ratio of the frequency tuning range to the FSR.
greater than 1.0. However, as the ratio becomes less than 1.0, the errors increase rapidly. Figure 5 shows the standard deviation of the phase errors. The deviation is 0.07° when the ratio is 1.0. However, it increases to 8° when the ratio is 0.95. Hence, we can conclude that the phases can be measured correctly if Δσ is greater than the FSR.
The relative amplitudes were also measured for various ratios Δσ /FSR. Figure 6 shows the measured amplitudes when Gaussian distribution was given to an AWG. The measured amplitudes are in good agreement with the given amplitudes even when Δσ is less than the FSR. We got the similar results even if we start the Monte-Carlo process with other initial amplitude distribution. When a sinc or uniform distribution was given to the amplitudes, the measurement errors were also negligible regardless of the initial conditions. Hence, we can conclude that amplitude measurement is insensitive to the experimental condition.
When we used a trial-and-error method in Fourier space, we could measure the phases precisely only if Δσ was larger than 1.7FSR and the phase errors were given within
Lager Image
Measured(symbols) and given(line) amplitudes of the waveguide array when Gaussian distribution was given to the AWG.
50° [11] . And even in that case, the measurement errors were very large for waveguides with amplitudes lower than 20 % of the peak value. The standard deviation of the phase measurement errors was greater than 1°, which is more than ten times of that in intensity space. Above all things, computational time of the Monte-Carlo process is much longer in Fourier space than in intensity space. Hence, it would be more advantageous to measure them in the intensity space. The poor results in Fourier space are believed to be related to the errors in calculating the Fourier transform.
IV. CONCLUSION
We show that it is possible to measure the phases of an AWG directly from the interferogram by a trial-and-error method. The phase measurement error is as low as 0.1° if the laser tuning range is greater than FSR. Moreover, this method shows high performance even for waveguides at the edges of the array, which tend to have very small amplitudes. All the characteristics are superior to the analysis using the Fourier transform.
Acknowledgements
This paper was supported by Wonkwang University in 2010.
References
Doerr C. R. , Okamoto K. (2006) “Advances in silica planar lightwave circuit” J. Lightwave Technol. 24 4763 - 4789
Sim E. C. , Abbou F. M. , Faidz A. R. (2007) “System degradation due to phase error induced crosstalk in WDM optical networks employing arrayed waveguide grating multi/demultiplexer” Opt. Quant. Electron. 39 553 - 560
Takada K. , Takadav Y. , Yokota A. , Satoh S. “Metal mask fabrication with an inkjet printer for AWG phase trimming using a photosensitive refractive index change” IEEE Photon. Technol. Lett. 17 813 - 815
Takada K. , Abe M. , Shibata M. , Ishii M. , Okamoto K. (2001) “Low-corsstalk 10-GHz-spaced 512-channel arrayed-waveguide grating multi/demultiplexer fabricated on a 4-in wafer” IEEE Photon. Technol. Lett. 13 1182 - 1184
Takada K. , Yamada H. , Inoue Y. (1996) “Optical low coherence method for characterizing silica-based arrayed-waveguide gratings multiplexers” J. Lightwave Technol. 14 1677 - 1689
Lazaro J. A. , Wessel R. , Koppenborg J. , Dudziak G. , Blewett I. J. (2003) “Inverse Fourier tramsform method for characterizing arrayed-waveguide grating” IEEE Photon. Technol. Lett. 15 93 - 95
Chen W. , Chen Y. J. , Yan M. , McGinnis B. , Wu Z. (2003) “Improved techniques for the measurement of phase error in waveguide based optical devices.” J. Lightwave Technol. 21 198 - 205
Takada K. , Okamoto K. (2000) “Frequency-domain measurement of phase error distribution in narrow-channel arrayed waveguide grating” Electron. Lett. 36 160 - 161
Takada K. , Satoh S. (2006) “Method for measuring the phase error distribution of a wideband arrayed waveguide grating in the frequency domain” Opt. Lett. 31 323 - 325
Takada K. , Hirose T. (2009) “Phase modulation method for AWG phase-error measurement in the frequency domain” Opt. Lett. 34 3914 - 3916
Go C. S. , Oh Y. H. , Lim S. (2011) “Measurement of the phase errors of AWG by using the Monte-Carlo Analysis.” Korean J. Opt. Photon. 22 207 - 213