Arbitrary optical waveform generation using

Arbitrary optical waveform generation using2D ring resonator arraysBing Xia and Lawrence R. ChenPhotonic Systems Group, Department of Electrical and Computer Engineering,McGill University, Montreal, QC H3A 2A7, Canadabing.xia@mail.mcgill.caAbstract: The direct temporal domain approach can be applied forarbitrary optical waveform generation using 2D ring resonator arrays(RRAs). To demonstrate the approach, we provide numerical exampleswhich show the generation of two very different waveforms from the sameinput pulse. In particular, we consider a hyperbolic secant input pulse with8 ps full width half maximum and generate (1) a 50 ps square-likewaveform with 5 ps rising and falling times and a 40 ps flat-top as well as(2) a 60 ps triangular waveform with 30 ps rising and falling times, bothwith a 5×5 RRA. Simulations show that the generated waveforms are well-matched to their targets.©2006 Optical Society of AmericaOCIS codes: (070.0070) Optical signal processing; (320.5540) Pulse shapingReferences and links1. A. M. Weiner, “Femtosecond optical pulse shaping and processing,” Prog. Quantum Electron. 19, 161-237(1995).2.J. Azana and M. A. Muriel, “Temporal self-imaging effects: Theory and application for multiplying pulserepetition rates,” IEEE J. Sel. Top. Quantum Electron. 7, 728-744 (2001).3. D. E. Leaird, S. Shen, A. M. Weiner, A. Sugita, S. Kamei, M. Ishii, and K. Okamoto, “Generation of high-repetition rate WDM pulse trains from an arrayed-waveguide grating,” IEEE Photon. Technol. Lett. 13, 221-223 (2001).4.Z. Jiang, D. E. Leaird and A. M. Weiner, “Line-by-line pulse shaping control for optical arbitrarywaveform generation,” Opt. Express 13, 10431-10439 (2005).5.J. D. McKinney, D. Seo, D. E. Leaird, A. M. Weiner, “Photonically assisted generation of arbitrarymillimeter-wave and microwave electromagnetic waveforms via direct space-to-time optical pulse shaping,” J. Lightwave Technol. 21, 3020-3028 (2003).6. F. Parmigiani, P. Petropoulos, M. Ibsen and D. J. Richardson, “All-optical pulse reshaping and retimingsystems incorporating pulse shaping fiber Bragg grating,” J. Lightwave Technol, 24, 357-364 (2006).7.P. Petropoulos, M. Ibsen, A. D. Ellis and D. J. Richardson, “Rectangular pulse generation based on pulsereshaping using a superstructured fiber Bragg grating,” J. Lightwave Technol. 19, 746-752 (2001).8. A. Rundquist, A. Efimov and D. H. Reitze, “Pulse shaping with the Gerchberg-Saxton algorithm,” J. Opt.Soc. Am. B 19, 2468-2478 (2002).9. A. M. Weiner, S. Oudin, D. E. Leaird and D. H. Reitze, “Shaping of femtosecond pulses using phase-onlyfilters designed by simulated annealing,” J. Opt. Soc. Am. A 10, 1112-1120 (1993).10. B. Xia and L. R. Chen, “A direct temporal domain approach for pulse-repetition rate multiplication witharbitrary envelope shaping,” IEEE J. Sel. Top. Quantum Electron. 1, 165-172 (2005).11. C. K. Madsen, C. Lenz, A. J. Bruce, M. A. Capuzzo, L. T. Gomez, T. N. Nielsen, and I. Brener,“Multistage dispersion compensator using ring resonators,” Opt. Lett. 24, 1555-1557 (1999).12. A. Rostami, and G. Rostami, “All-optical implementation of tunable low-pass, high-pass, and band-passoptical fitlers using ring resonators,” J. Lightwave Technol. 23, 446-460 (2005).13. C. K. Madsen, J. H. Zhao, Optical filter design and analysis-A signal processing approach ( John Wiley &Sons, 1999), Chap.5.14.Y. M. Landobasa, S. Darmawan and M. K. Chin, “Matrix analysis of 2-D microresonator lattice opticalfilters,” IEEE J. Quantum Electron. 41, 1410-1418 (2005).15. A. Agarwal, P. Toliver, R. Menendez, S. Etemad, J. Jackel, J. Yong, T. Banwell, B.E. Little, S. T. Chu, W.Chen, W. Chen, J. Hryniewicz, F. Johnson, D. Gill, O. King, R. Davidson, K. Donovan and P. J. Delfyett, “Fully programmable ring-resonator-based integrated photonic circuit for phase coherent applications,” J.Lightwave Technol. 24, 77-87 (2006).#71361 - $15.00 USD Received 26 May 2006; revised 14 July 2006; accepted 17 July 2006 (C) 2006 OSA24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 661916.T. Sakamoto, F. Futami, K. Kikuchi, S. Takeda, Y. Sugaya and S. Watanabe, “All-optical wavelengthconversion of 500-fs pulse trains by using a nonlinear-optical loop mirror composed of a highly nonlinear DSF,” IEEE Photon. Technol. Lett. 13, 502-504 (2001).1. IntroductionTechniques for the generation, control, and manipulation of optical pulses attract considerable interest for numerous applications and have become increasingly important in many scientific areas [1]. Of specific interest are techniques for pulse repetition rate multiplication (PRRM), which are used to generate ultra-fast optical pulse trains from a low repetition rate input pulse train [2, 3], as well as for the generation of arbitrary waveforms [4]. Arbitrary or user-defined waveforms at high repetition rates can be used to generate wide-band RF signals [5], as switching windows or gating signals in all-optical pulse reshaping and retiming systems in optical communications [6], and as probing signals for various studies on spectroscopy and the properties of materials or molecules [7].Traditional pulse shaping methods are based on frequency domain processing (i.e. spectral filtering) in which we specifically manipulate the different spectral components of the input pulse in amplitude and/or phase. However, the relationship between the input pulse spectrum and the target temporal waveform is not straightforward, especially for phase-only filtering processes. For example, if we need to generate an arbitrary waveform from a given input pulse, say Gaussian, it is difficult to determine the phase response and the corresponding filter that generates the waveform. Efficient Fourier synthesis algorithms are necessary to synthesize the required spectral response as well as the filter technologies for their implementation. Several optimization and synthesis algorithms have been used to generate arbitrary waveforms, including genetic, Gerchberg-Saxton [8], and simulated annealing [9, 10] algorithms. In terms of hardware, one of the most popular implementations is based on a bulk optical 4f pulse shaper, which separates spatially the frequency components of an input pulse and uses amplitude and/or phase masks to process the signal [1].In this paper, we propose an alternate method to synthesize arbitrary waveforms using the direct temporal domain approach. In this approach, traditional frequency domain processing is not used and the waveform is synthesized purely in the temporal domain. While simple and compact optical filters structures are used, the waveform is synthesized by optimizing the filter parameters to obtain a specific temporal impulse response, without concern for the corresponding frequency response. The direct temporal domain approach was initially proposed to generate an ultra-high repetition rate pulse train from a low repetition rate pulse train using spectrally-periodic (SP) filters, details on the theory can be found in [10]. In this paper, we show that this approach can also be applied for general pulse shaping purposes, for example to generate a square waveform from a hyperbolic secant pulse train or a Gaussian pulse train. This approach is very different from the traditional spectral filtering approach which uses spectral filters to shape the input pulse spectrum into a sinc function (with amplitude and phase information preserved) in order to achieve a square waveform in the temporal domain. For implementation, we choose 2D ring resonator arrays (RRAs) as the optical filter to perform the temporal waveform generation since ring resonators are a class of SP filters.2. Theory and configurationHigh repetition rate pulse trains with uniform or even arbitrary envelopes can be generated using an SP filter and the direct temporal domain approach, as long as the input pulse width is sufficiently narrow [10]. When the free spectral range (FSR) of the SP filter satisfies certain conditions with respect to the new (output) repetition rate, pulse repetition rate multiplication (PRRM) is achieved; moreover, the output multiplied pulses (within the original period) can have very different amplitudes and phases due to the simultaneous amplitude and/or phase filtering process. Moreover, these amplitudes and phases can be set independently of the input #71361 - $15.00 USD Received 26 May 2006; revised 14 July 2006; accepted 17 July 2006 (C) 2006 OSA24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6620pulse shape, see Eqs. 4-7 in [10]. We have found that this characteristic can be used to generate and synthesize arbitrary waveforms at the rate of the original input pulse train. The basic principle of the direct temporal domain approach for arbitrary waveform generation is illustrated in Fig. 1. First, we choose an SP filter to perform PRRM with a uniform output envelope. When an input pulse train with a broader pulse width is launched into the same SP filter, PRRM takes place, but there will be an interference among the output pulses. Next, we can manipulate and optimize the amplitude and phase of each individual output pulse (in the multiplied train) to generate the specified waveform (at the input repetition rate). The key in this approach is that we must find an SP filter which can simultaneously perform PRRM and manipulate the amplitude and phase of each individual output pulse with a range of freedom. The amplitude and phase of each individual output pulse is controlled by specifying the value of the temporal impulse response of the filter. This approach is very different from general spectral filtering approach which normally uses spectral filters to shape the input pulse spectrum into a special shape (with amplitude and phase information preserved) in order to obtain the desired waveform in the temporal domain.Multi-stage ring resonators are a powerful class of optical filters which have been widely investigated for WDM add/drop filters, dispersion compensation, and many other applications that are focused on their spectral characteristics [11, 12]. However, their temporal domain characteristics are equally interesting. Since ring resonators are SP filters, they can be used as key components to perform arbitrary waveform generation by implementing the direct temporal domain approach.Figure 2(a) shows a typical configuration of an M×N 2D RRA, where M represents the number of rings in the vertical direction and N represents the number in the horizontal direction. Figure 2(b) shows the details of the individual ring elements which incorporate a directional coupler with splitting ratio κand a phase shifter with an additional phase shift φ. Note that the rings are assumed to be identical in size and are coupled in the vertical direction only. While there is no coupling between the rings in the horizontal direction, the signals propagate in the horizontal direction through the two waveguide buses that are placed at the top and the bottom of the RRA. Since the waveguide bus introduces time delays above and beyond those of the rings, the length of the bus must satisfy specific conditions in order to generate a pair of pulse trains at the two output ports simultaneously. In particular, if L Fig. 1. Schematic of arbitrary waveform generation using the directtemporal domain approach; (a) PRRM output with a uniform envelopefrom an input pulse train with narrow pulse widths; (b) PRRM output froman input pulse train with wide pulse widths; (c) waveform generation withan optimized SP filter and an input pulse train with wide pulse widths.P R R M w i t h a w i d en(b)r m g e n e r a t i o n w i t ht i m i z e d S P f i l t e r(c)#71361 - $15.00 USD Received 26 May 2006; revised 14 July 2006; accepted 17 July 2006(C) 2006 OSA 24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6621denotes the distance between the centers of two adjacent rings, then L must be a positive multiple of πR where R is the radius of the ring. Larger values of L will lead to more loss and an overall increase in the size of the structure; hence, we choose L = πR . The free spectral range (FSR) of the device is given by )2/(R e n c FSR ⋅=π, where R is the radius of the ring and n e is the effective index of the waveguide.The RRA transfer function can be calculated by dividing the M×N array into N horizontally-cascaded M×1 arrays, in which each M×1 array can be calculated separately using transfer matrices and Z-transforms as described in [11, 12]. For the n th column, the transfer matrix in the vertical direction can be expressed by⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⋅⋅⋅⋅=⎥⎦⎤⎢⎣⎡+22222112112212111O I O I O I m θθθθθθθ ⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−⎟⎟⎟⎠⎞⎜⎜⎜⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−=+++=−−−−−−∏221121122221111)1(O I t t t j e e γt e e γt e e γt j m m m m p R πβj φj p R πβj φj p R πβj φj p p p p where I 1, I 2, O 1 and O 2 represent the inputs and outputs in the two waveguide buses as shown in Fig. 2(a); p t is the through-amplitude coefficient which is defined as p p t κ−=1; γ is the waveguide loss and β is the propagation constant in the device. This matrix can then be converted to another transfer matrix which functions in the horizontal direction using the following transformation:⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−=⎥⎦⎤⎢⎣⎡Φ=⎥⎦⎤⎢⎣⎡211121122211221221211)(1I I I I O O n θθθθθθθ Fig. 2. (a) General configuration of an M×N RRA and (b) detailedview of the individual rings, r is the radius of the ring resonator, κ isthe coupling coefficient and Φ is an additional phase shift.(a)Input 2 (M=even)Output 2 (M=even) Output 1φφ(b) (1) (2)#71361 - $15.00 USD Received 26 May 2006; revised 14 July 2006; accepted 17 July 2006(C) 2006 OSA 24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6622where ij θare the elements in the vertical transfer matrix. The two waveguide buses can also beexpressed with a transfer matrix in the horizontal direction: ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=Φ⋅−−γγββd L j L j bus e e 00. d is the signal propagation direction and is equal to -/+ 1 when the signal counter-/co-propagates withthe input, i.e. when M is odd/even. The overall RRA transfer function is then given by the system matrix 11Φ⋅Φ⋅⋅⋅Φ⋅Φ⋅Φ=Θ−bus N bus N . By properly designing the parameters p t and p φ , the RRA can perform PRRM and simultaneously manipulate amplitudes and phases forthe newly generated output pulses, thereby generating a special waveform at output. It is important to stress that the parameters of the RRA are optimized to generate a specific temporal impulse response of the filter in order to control the amplitude and phase of the output pulses.In order to find a proper set of values of p t and p φ to generate a specified waveform, apowerful optimization algorithm is highly desirable. We use a similar optimization process as described in [10], but with some modifications. The major difference is that all sampling points within a period at the output (i.e. the input period) are compared with the target waveform and a weighted root mean square error function is used to evaluate the optimized results. Furthermore, in our optimization process, we use quantized rather than arbitrary values for p t and p φ, since values with a high precision may be difficult to realize in practice.3. Simulation resultsAs an example, we consider a 5×5 RRA for arbitrary waveform generation. A 10 GHz train of hyperbolic secant input pulses is launched at input 1; we assume no input signal at input 2. The full width half maximum (FWHM) of each hyperbolic secant pulse is 8 ps as shown in Fig. 3(a). The objective is to transform this hyperbolic secant pulse into a square waveform with 5 ps rise and fall times, and a 40 ps flat-top. Square or rectangular shaped pulses are useful for enhancing the operation of nonlinear optical switches for reshaping ultrashort optical pulses [6].In our simulation, we set the FSR of the filter to 160 GHz (corresponding ring radius is 0.199 mm for an effective index n e = 1.5) and assume that the total loss in the whole device is 30%. We use the simulated annealing algorithm and a weighted root mean square error function to optimize the filter parameters, namely the through-amplitude p t and the additionalphase shift p φfor each ring element. We quantize p t and p φin steps of 0.01 and π/16, respectively. The weighted root mean square error function is defined as⎟⎟⎠⎞⎜⎜⎝⎛−=∑=N n et t output n I I w error 12arg where I is the intensity, N is the total sampling number in one repetition period, and w n is the weight for each sampling point.The results are shown in Fig. 3(b) and (c); the corresponding filter parameters are given in Table 1. Figure 3(b) shows the generated square-like waveform. The rising and falling edges are well matched to the target pulse and there are small oscillations in the flat-top portion. Figure 3(c) is a re-plot of (b) but in dB scale. It can be seen that the pulse extinction ratio (ER), defined as the ratio of the smallest intensity in flat-top and the largest intensity in background, is around 25 dB. The peak-to-peak variation in the flat-top is 0.9 dB.#71361 - $15.00 USDReceived 26 May 2006; revised 14 July 2006; accepted 17 July 2006(C) 2006 OSA 24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6623Fig. 3. (a) Input and target waveform in one repetition period; (b) The generated square-like waveform; (c) The generated waveform in dB scale.(a)(b) (c)Fig. 4. Spectrum of the 10 GHz pulse train at the RRA input and of the square-like waveform at the output. Amplitude (a) and (c); phase (b) and (d).(a) (c) (b)(d)#71361 - $15.00 USDReceived 26 May 2006; revised 14 July 2006; accepted 17 July 2006(C) 2006 OSA 24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6624Figure 4 shows the spectra (amplitudes and phases) of the 10 GHz hyperbolic secant pulse train at the input and the square-like waveform at the output. It can be seen that while the amplitude response of the output waveform exhibits the same 10 GHz spectral spacing as the input, the magnitudes for each individual mode have been altered which shows the amplitude-filtering nature of the RRA. In addition, the input and output phases are completely difference which shows the phase-filtering nature of the RRA. Since our approach only process the waveform in the temporal domain, the spectral shape of the output (and consequently of the filter) is not important in our approach and optimization algorithm; both magnitude and phase can be manipulated by RRA with a large scale of freedom (note in particular that the output spectrum does not exhibit any sinc-like features). This demonstrates the flexibility of the temporal domain approach over amplitude and phase filtering approaches which needs to alter the spectrum to a sinc shape (with appropriate phase jumps), or the phase-only filtering which is constrained to maintaining a constant magnitude or a uniform loss for all wavelengths.Table 1. Parameters of the RRA configuration for generation of a square waveform and a triangularwaveform from a 10 GHz hyperbolic secant pulse train.Fig. 5. (a) The input hyperbolic secant pulse train at 10 GHz and the target waveform;(b) the target waveform and the generated triangular waveform form. the RRA.(b) (a)) #71361 - $15.00 USD Received 26 May 2006; revised 14 July 2006; accepted 17 July 2006(C) 2006 OSA 24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6625As a second example to demonstrate the capabilities of the approach, we use the RRA to generate a triangular waveform which has 30 ps rising and falling times from the same input hyperbolic secant pulse train. The FSR of the RRA is still 160 GHz. The optimized RRA parameters p t and p φ are also given in Table 1. Figure 5 shows the generated triangularwaveform, which is well matched to the target. We obtain similar results when generating triangular waveforms of different widths, which demonstrate that this approach is extremely powerful for generating arbitrary waveformsWe have also investigated the impact of the input pulse width (FWHM) on the quality of the generated waveform. If the FWHM varies in a range of ± 0.8 ps from its ideal value (8 ps), a good square-like waveform can be maintained with the same set of filter parameters. Otherwise, we need to re-optimize the ring resonator parameters. In particular, when the FWHM of the input pulse is between 6 ps to 9 ps, a good square-like waveform can be synthesized using the optimization algorithm. However, if the FWHM < 6 ps, we cannot limit the peak-to-peak variations to < 1 dB in the flat-top, even after re-optimization. Moreover, if the FWHM > 9 ps, the rising and falling times (5 ps) cannot be achieved. We also investigated the use of Gaussian input pulses to generate the square-like waveform. The result is similar to that obtained using hyperbolic secant pulses, which demonstrate that 2D RRAs have a big potential to synthesize square-like waveforms and even arbitrarily defined waveforms using the direct temporal domain approach.Next, we investigate the impact of fabrication errors in p t and p φ of each ring element.For illustrative purposes, we consider the same 5×5 RRA that generates the square-like waveform. We assume that the fabrication errors are uniformly distributed in a range of 0 to2.5% for p t and 0 to 1% for p φ of all rings in the RRA and obtain both ER and the peak-to-peak intensity variations in the flat-top portion based on these imperfect parameters. We repeat the simulations 10,000 times and determine the average ER for the device. Figure 6(a) shows a contour plot for the average extinction ratio as the fabrication errors in p t and p φvary from 0 to 2.5% and 0 to 1%, respectively; Fig. 6(b) shows the contour plot for thepeak-to-peak variations in the flat-top portion. It can be seen that in order to maintain both an ER greater than 15 dB and peak-to-peak variations less than 2.5 dB, the fabrication error of p t and p φ should be within 2% and 0.8%, respectively. We note that programmable ring-resonator-based integrated photonic circuit with an accurate vertical coupling structure for p t and independent controls for phase shift p φ have been demonstrated [15], which indicate theFig. 6. Contour plots for fabrication errors in the RRA for the square waveformgeneration; (a) the average extinction ratio; (b) peak-to-peak intensity variation in theflat-top portion.(a)(b)#71361 - $15.00 USDReceived 26 May 2006; revised 14 July 2006; accepted 17 July 2006(C) 2006 OSA 24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6626practicality and feasibility of our temporal domain approach for arbitrary waveform generation.4. SummaryWe have demonstrated the use of 2D RRAs to synthesize square and triangular waveforms with the direct temporal domain approach. A square-like waveform with a 5 ps rising time, a 40 ps flat-top period, and a 5 ps falling time is generated by this filter from a hyperbolic secant input pulse with a FWHM of 8 ps. A high extinction ratio (25 dB) is observed in the generated square-like waveform and the peak-to-peak variation in the flat-top portion is 0.9 dB. In addition, we have shown the generation of a triangular waveform with a 30 ps rising time and a 30 ps falling time. These results show that the direct temporal domain approach can be applied to perform optical signal processing and pulse shaping and arbitrary waveforms using 2D RRAs.AcknowledgmentsThis research was supported in part by the Natural Sciences and Engineering Research Council of Canada and industrial partners through the Agile All-Photonic Networks (AAPN) Research Program.#71361 - $15.00 USD Received 26 May 2006; revised 14 July 2006; accepted 17 July 2006 (C) 2006 OSA24 July 2006 / Vol. 14, No. 15 / OPTICS EXPRESS 6627。