Modelling of NOLM Demultiplexers Employing Optical Soliton Control Pulse

Similar documents
Soliton-Effect Optical Pulse Compression in Bulk Media with χ (3) Nonlinearity. 1 Introduction

3D FE Modeling Simulation of Cold Rotary Forging with Double Symmetry Rolls X. H. Han 1, a, L. Hua 1, b, Y. M. Zhao 1, c

Engineering Approach to Modelling Metal THz Structures

Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion

Beam Combining and Atmospheric Propagation of High Power Lasers

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322

Optical property of Few-Mode Fiber with Non-uniform Refractive Index for Cylindrical Vector Beam Generation

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function

Hz transmission peak, is : F 100

ChE 471: LECTURE 4 Fall 2003

Least Squares Optimal Filtering with Multirate Observations

Electric Current and Resistance

Dead-beat controller design

1 The limitations of Hartree Fock approximation

Figure 1a. A planar mechanism.

A study on GPS PDOP and its impact on position error

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Introduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem

rcrit (r C + t m ) 2 ] crit + t o crit The critical radius is evaluated at a given axial location z from the equation + (1 , and D = 4D = 555.

Displacement and Deflection Sensitivity of Gas-coupled Laser Acoustic. Detector

General Chemistry II, Unit II: Study Guide (part 1)

Name: Period: Date: ATOMIC STRUCTURE NOTES ADVANCED CHEMISTRY

Design and Simulation of Dc-Dc Voltage Converters Using Matlab/Simulink

Dispersion Ref Feynman Vol-I, Ch-31

Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018

Introduction to Smith Charts

Determining the Accuracy of Modal Parameter Estimation Methods

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

( ) kt. Solution. From kinetic theory (visualized in Figure 1Q9-1), 1 2 rms = 2. = 1368 m/s

EHed of Curvature on the Temperature Profiles

Linearization of the Output of a Wheatstone Bridge for Single Active Sensor. Madhu Mohan N., Geetha T., Sankaran P. and Jagadeesh Kumar V.

Aircraft Performance - Drag

Comparison of hybrid ensemble-4dvar with EnKF and 4DVar for regional-scale data assimilation

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell

, which yields. where z1. and z2

ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES

LEARNING : At the end of the lesson, students should be able to: OUTCOMES a) state trigonometric ratios of sin,cos, tan, cosec, sec and cot

Chapter 30. Inductance

Theoretical study of third virial coefficient with Kihara potential

Support-Vector Machines

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India

the results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must

TOPPER SAMPLE PAPER 2 Class XII- Physics

Heat Management Methodology for Successful UV Processing on Heat Sensitive Substrates

Q1. In figure 1, Q = 60 µc, q = 20 µc, a = 3.0 m, and b = 4.0 m. Calculate the total electric force on q due to the other 2 charges.

Spectral Gain-Carrier Density Distribution of SQW GaAs/AlGaAs Laser

Verification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance

Bicycle Generator Dump Load Control Circuit: An Op Amp Comparator with Hysteresis

EXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE

Lecture 17: Free Energy of Multi-phase Solutions at Equilibrium

Kinetics of Particles. Chapter 3

A PLETHORA OF MULTI-PULSED SOLUTIONS FOR A BOUSSINESQ SYSTEM. Department of Mathematics, Penn State University University Park, PA16802, USA.

Computational modeling techniques

THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

ON THE COMPUTATIONAL DESIGN METHODS FOR IMPROOVING THE GEAR TRANSMISSION PERFORMANCES

Department of Electrical Engineering, University of Waterloo. Introduction

Lab 11 LRC Circuits, Damped Forced Harmonic Motion

Math Foundations 20 Work Plan

ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT

Chapter 3: Cluster Analysis

Chapter 4. Unsteady State Conduction

Chapter 16. Capacitance. Capacitance, cont. Parallel-Plate Capacitor, Example 1/20/2011. Electric Energy and Capacitance

Examiner: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data

Aerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison

Progress In Electromagnetics Research M, Vol. 9, 9 20, 2009

CBSE Board Class XII Physics Set 1 Board Paper 2008 (Solution)

Differentiation Applications 1: Related Rates

MATHEMATICS SYLLABUS SECONDARY 5th YEAR

Series and Parallel Resonances

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

ENSC Discrete Time Systems. Project Outline. Semester

1/2 and e0 e s ' 1+ imm w 4 M s 3 πρ0 r 3 m. n 0 ktr. .Also,since n 0 ktr 1,wehave. 4 3 M sπρ 0 r 3. ktr. 3 M sπρ 0

[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )

Synchronous Motor V-Curves

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

AP Physics Kinematic Wrap Up

Performance Bounds for Detect and Avoid Signal Sensing

Characterisation of a Pyramid WFS: an experimental study

Lead/Lag Compensator Frequency Domain Properties and Design Methods

Phys102 Final-061 Zero Version Coordinator: Nasser Wednesday, January 24, 2007 Page: 1

Optimization of frequency quantization. VN Tibabishev. Keywords: optimization, sampling frequency, the substitution frequencies.

COMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d)

GENERAL FORMULAS FOR FLAT-TOPPED WAVEFORMS. J.e. Sprott. Plasma Studies. University of Wisconsin

CLASS XI SET A PHYSICS

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b

ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS

On Boussinesq's problem

Function notation & composite functions Factoring Dividing polynomials Remainder theorem & factor property

Enhancing Performance of MLP/RBF Neural Classifiers via an Multivariate Data Distribution Scheme

Power Flow in Electromagnetic Waves. The time-dependent power flow density of an electromagnetic wave is given by the instantaneous Poynting vector

Phung LeCong, Al Joseph, Erie Udd, and Paul Theriault MeDonnell Douglas Astronauties Company 5301 Bolsa Avenue Huntington Beaeh, California 92647

Multiband retardation control using multi-twist retarders

CHAPTER 2 Algebraic Expressions and Fundamental Operations

Solution to HW14 Fall-2002

8 th Grade Math: Pre-Algebra

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

Transcription:

Micwave and Optical Technlgy Letters, Vl. 1, N. 3, 1999. pp. 05-08 Mdelling f NOLM Demultiplexers Emplying Optical Slitn Cntrl Pulse Z. Ghassemly, C. Y. Cheung & A. K. Ray Electrnics Research Grup, Schl f Engineering, Sheffield Hallam university, Pnd St., Sheffield, S1 1WB, UK E-mail: z.f.ghassemly@shu.ac.uk ABSTRACT: An ptical slitn pulse may be used as a cntrl signal in nn-linear ptical lp mirrr demultiplexers in rder t reduce timing-jitter nise and crss-talk. A mathematical mdel fr calculating the width f the slitn switching windw and ptimizing critical parameters such as walk-ff time and pulse width has been presented, and its accuracy is verified by slving the nnlinear Schrdinger equatin. Key wrd: Semicnductr laser amplifier, ptical switching, slitn, ptical time divisin demultiplexer, ptical lp mirrrs. 1- INTRODUCTION The nn-linear ptical lp mirrr (NOLM) demultiplexer is a prmising cnfiguratin fr achieving all-ptical time divisin demultiplexing because f its high perating speed [1]. Channel demultiplexing is realised by the phase difference beeen the clckwise (CW) and cunter-clckwise (CCW) signal pulses prpagating within the fibre lp, see Figure 1. In ultra-high speed systems, there are majr prblems: 1) difficulty in achieving cmplete switching f the signal pulses; and ) timing jitter beeen the cprpagating cntrl and signal pulses. The latter is reduced t a certain extent by intrductin f walk-ff time T w beeen the cntrl and signal pulses. This results in a square switching windw shape, with a width f T w L, where L is the fibre lp length, rather than bell shape, thus enabling imprved switching and cnsequently better tlerance t timing jitter []. A 100% switching is difficult t achieve because f the cntrl pulse experiencing a slitn cmpressin effect due t interactin beeen the fibre dispersin and self phase mdulatin (SPM) [4]. This will result in an increased cntrl signal peak pwer which in turn leads t an asymmetrical switching windw with an increasing phase shift twards the end f the switching prfile. Maximum switching may be achieved by emplying ptical slitn as a signal r cntrl pulse, since the fundamental prperties f the slitn are unifrm phase ver the entire pulse and cnstant pulse shape ver the entire prpagatin length. Here, we investigate the later ptin, and shw that the prblem f pulse defrmatin during prpagatin within the lp is slved, thus resulting in a much imprved switching windw with high transmittance and reduced timing jitter effects. - THEORY T study the utput transmittance f the NOLM demultiplexer it is best t cnsider the cmbined impact f the fllwing effects: 1) crss phase mdulatin (XPM), ) prpagatin f ptical slitn pulse within the fibre lp; and 3i) walk-ff time beeen the cntrl and signal pulses. The switching windw is created by the phase difference beeen the CW and CCW signals. The phase f the signal pulse cprpagating with the cntrl pulse is substantially changed by the XPM effect. With the inclusin f walk-ff time, the phase change can be represented by L φ = γ PT ( Tw x ) dx (1) 0 where γ is the nn-linear cefficient, P is the ptical pwer prfile, and T w = the walk-ff time per unit length. The cntrl pulse has t meet the fllwing slitn cnditins in rder t prpagate undisrtedly alng a lssless fibre [6]:

1- the prpagatin equatin f a fundamental slitn wave u( ξ, τ) = sec h( τ)exp( iξ/ ) - the peak pwer P = β / γt where u is the nrmalised amplitude, τ = T/T and ξ = z/l D, z is the travelling distance, L D is the dispersin length = T / β, T FWHM 1763. is the pulse width f the ptical slitn cntrl pulse, and β is the first rder dispersin cefficient. The cntrl pulse walks thrugh the signal pulse due t difference in the grup velcity, and with the phase change f the signal depending n the value f the ttal walk-ff time, the time varying ptical pwer prfile in Eq. (1) can be replaced by the pulse average pwer ver the duratin f ttal walk-ff time. The maximum average pwer f the cntrl pulse, resulting in XPM, is given as P ave(max) = T PT ( ) dt T / / T () Fr slitn pulse, P(T) = P. u, therefre the peak phase change can be shwn t be T / T φ p = [ γlp T sec h ( T / T ) d( T / T )]/ T (3) T / T where T = T w L is the ttal walk-ff time. Substituting fr P and replacing φ p by π fr maximum transmittance and slving fr sech (x), we btain an expressin fr lss less fibre lp as ( T / T ) T T π = L β tanh (4) 4 The utput transmittance wuld be equal t 1 if the cnditin set by Eq.(4) is satisfied. In Eq. (4) it is assumed that fibre lp is lss free. Hwever, in NOLM demultiplexers where a lng length f fibre is used, it is necessary t include fibre lss α in the analysis. Assuming that the rate f decrease f pwer is s small, the increase in chrmatic dispersin can be neglected, and therefre the peak pwer can be expressed as P 1 L = L P ( 0)exp( α z ) dz (5) 0 Slving Eq. (5) and substituting it int Eq. (3), we can btain a mdified expressin fr maximum transmittance as [ 1 exp ( L )] tanh ( T / T ) T T πα = 4 β α (6) The utput transmittance wuld be equal t 1 if the cnditin set by Eq. (4) is satisfied. In NOLM emplying a Gaussian pulse, the switching windw prfile is a bell shape with a nn-flat tp. Therefre, relative timing jitter beeen the cntrl and signal pulses and cntrl pulse defrmatin within the lp will induce intensity fluctuatin f the demultiplexed utput signal, thus resulting in switching pwer penalty [5], whereas with an ptical slitn pulse, the shape f the switching windw becmes symmetrical with a flat-tp transmittance prfile, thus resulting in reduced timing jitter nise

and inter-channel crss-talk. T achieve the frmer, the variatin f transmittance at the tp f the switching windw has t be kept within the perturbatin-tlerance factr F p ; see Figure. The shaded areas represent the prtin f the cntrl pulse prfile walking thrugh the signal pulse alng the whle prpagating distance. At time t 1, the transmittance level is belw the threshld level. Frm t 1 t t 3, it remains abve the threshld level and reaches peak value at time t. After t 3, the transmittance level begin t drp. Therefre, the width f the transmissin windw is equal t t 3 - t 1. By cmparing the walk-ff regin f the pulse prfile at t 1, t, and t 3, the walk-ff regin f the pulse prfile at t 1 is just left shifted with t = (T /) -T 1 and t 3 is just right shifted with t = T -(T /) frm the symmetric walk-ff regin f the pulse prfile at t, respectively. T -(T /) = (T /) -T 1 because f the symmetrical shape f a slitn pulse. The transmissin windw s perturbatin-tlerance width is defined as the perid f time at which the nrmalised transmittance level remains abve the threshld value 1 - F p, which is given as W = ) ( T / T1 ) = ( T T / ) = ( T / Ti (7) where T 1 and T are referred t in Figure. Fr F p = 0.5, W = FWHM width. A mdel was develped fr calculating the width f a slitn switching windw assuming that α = 0. Fr a given set f system parameters, F p can be related t T 1 and T by the fllwing equatin: 1 ( 4cs ( F )/ ) tanh( T / T ) = tanh( T / T ) tanh[( T T ) / T ] p π (8) i Equatin (8) yields slutins fr T i : ne is T 1, and the ther is fr T. W can then be calculated by substituting either f the slutins fr T i in Eq. (8) int Eq. (7). The expressin fr W hlds true n matter which ne f the slutins f T i is used in the substitutin. i 3- SIMULATION RESULTS AND DISCUSSION Equatin (4) was slved numerically by the Newtn-Raphsn methd (NRM) [3], and the slutin was checked by applying the crrespnding values f T and T t Eq. (4) fr verificatin. Fr L = 3 km, the results fr the slitn pulse width versus the ttal walk-ff time fr different values f β and α fr a slitn switching windw at maximum transmittance are shwn in Figure 3. The slitn pulse width decreases with the increase f the ttal walk-ff time. This is due t the decrease in the interactive time fr XPM beeen signal pulses and the central peak f cntrl pulses. The reduced pulse width cmpresses mre ptical pwer arund the central peak f the cntrl pulses, and as a result, the phase change φ is increased. The increased phase change due t reduced values f the slitn pulse width is exactly equal t the decreased phase change due t the larger value f the ttal walk-ff time, thus maintaining maximum transmittance. Figure 3 als shws the affects f the first rder dispersin cefficient β n the slitn pulse width and the ttal walk-ff time. Fr a fixed value f ttal walk-ff time, a smaller dispersin cefficient results in a reduced slitn pulse width. This is because the peak pwer f a ptical slitn pulse is directly prprtinal t the magnitude f the first rder dispersin cefficient β and inversely prprtinal t the square f the slitn pulse width. A reductin in peak pwer, due t lwer values f β, is partially cmpensated by the decrease f the slitn pulse width. Fr lssy fibre Eq. (6) is slved numerically and the results fr a range f α are als shwn in Figure 3. With the inclusin f α, the peak pwer f a pulse is gradually reduced as the pulse prpagates alng the fibre lp and, as a result, the phase change due t XPM is als reduced. Fr α = 0. db/km, a narrwer pulse width may be used t cmpensate fr the lss f phase due t α. Equatin (8) was slved numerically by a cmbinatin f the functin iteratin and NRM and the results fr T i were used t calculate the transmissin windw width. Values f the slitn pulse width and the ttal walk-ff time were btained frm Eq. (4), and the results fr the transmissin windw width fr different values f β are illustrated in Figure 4, shwing a threshld level abve which the slpe increases significantly. The threshld level dependence n the β can be used as an imprtant parameter t predict the width f a slitn switching windw. T create a wider switching windw, if the ttal walk-ff time is greater than the threshld level, then a smaller value f β can be used, and if the ttal walk-ff time is smaller than the threshld level, then a larger value f β may be used. By

increasing the width f the switching windw, the timing jitter nise can be reduced, but at the cst f increased inter-channel crss-talk by allwing sme prtin f the adjacent pulses t be demultiplexed t the utput prt, which is the subject f further study. Finally, the mathematical mdel develped is verified by slving the nn-linear Schrdinger equatin with beam prpagatin methd (BPM) [], and the results are shwn in Table 1. An interactive fibre length f 3 km and first-rder dispersin parameters f 5 ps /km are used fr calculatin, and the pulse width fr lss-free and lssy fibre are btained frm Eqs. (4) and (6), respectively with T = 10 ps. The results btained agree very clsely, thus verifying the accuracy f the mdel develped. Hwever, the deviatin f the windw width increases as the fibre lss per kilmetre increases. 4- CONCLUSIONS A mathematical mdel fr NOLM emplying a slitn cntrl pulse has been presented, where the prblem f pulse defrmatin during prpagatin within the lp is vercme. It has been shwn that the mdel can be used t ptimise critical parameters, such as walk-ff time and pulse width, fr maximising utput transmittance with much reduced timing jitter effects. The accuracy f the mathematical mdel develped is verified by slving the nn-linear Schrdinger equatin. Acknwledgement One f the authrs Mr C. Y. Cheung is financially supprted by Overseas Research Students Award. REFERENCES 1) M. Jinn, and T. Matsumt, Nnlinear Sagnac Interfermeter switch and its applicatins, IEEE J. f Quantum electrnics, 8 (199), 875-88. ) K. Uchiyama, et-al., Effects f cntrl-signal pulse walk-ff n BER perfrmance f nnlinear ptical lp mirrr demultiplexer, Electrn. Lett 9 (1993), 1870-1871. 3) P. R. Turner, Numerical analysis, Macmillan, New Yrk, 1994. 4) L. P. Barry, et-al., Effect f cntrl pulse defrmatin n the switching characteristics f a NOLM, Prcs. f the 1st Australian Cnf. n Optical Fibre Technlgy, Dec. 1996, pp. 197-00. 5) K. Uchiyama, et-al., Signal-t-nise rati analysis f 100 Gb/s demultiplexing using nnlinear ptical lp mirrr, IEEE J., Lighave Technl 15 (1997), 194-01. 6) G. P. Agrawel, Nnlinear fiber ptics, Academic Press, New Yrk, 1989. 1999 Jhn Wiley & Sns, Inc. CCC 0895-477/99

Cntrl cupler Cntrl Data in (signal CW Lng Fibre Lp CCW I/O Cupler Data ut Figure 1 A typical cnfiguratin f NOLM demultiplexer Transmissin windw At time t 1 At time t 3 F p T W T W T 1 T 1 -T W T -T W T At time t -T W / T W T W / Figure Typical slitn transmissin windw: the shaded areas representing the prtin f the cntrl pulse verlapping with the signal pulse thrughut the cmplete prpagatin

β = -0 ps /km β =-10 ps /km β = -5 ps /km Fibre Lss α 0 db/km 0. db/km 0.5 db/km 0.8 db/km Figure 3 Slitn pulse width versus ttal walk-ff time fr different values f β and α β (ps /km) -5-10 -15-0 Figure 4 Perturbatin-tlerance width versus ttal walk-ff time fr different values f β

TABLE 1 Results Using BPM fr L = 3 km and β = -5 ps /km Fibre lss (db/km) Parameter BPM Mathematic al Mdel 0 Transmittance 0.9998 1.0000 Windw width (F p = 5%) ps 6.73 6.75 0. Transmittance 1.0000 1.0000 Windw width (F p = 5%) ps 6.66 6.9 0.5 Transmittance 1.0000 1.0000 Windw width (F p = 5%) ps 6.3 7.17