Bound-soliton fiber laser

Similar documents
Mechanism of multisoliton formation and soliton energy quantization in passively mode-locked fiber lasers

Dissipative soliton resonance in an all-normaldispersion erbium-doped fiber laser

Observation of spectral enhancement in a soliton fiber laser with fiber Bragg grating

Self-started unidirectional operation of a fiber ring soliton. laser without an isolator

Dark Soliton Fiber Laser

Vector dark domain wall solitons in a fiber ring laser

Group interactions of dissipative solitons in a laser cavity: the case of 2+1

Bound states of gain-guided solitons in a passively modelocked

Vector dark domain wall solitons in a fiber ring laser

Induced solitons formed by cross polarization coupling. in a birefringent cavity fiber laser

arxiv: v1 [physics.optics] 26 Mar 2010

Dark Pulse Emission of a Fiber Laser

Compact graphene mode-locked wavelength-tunable erbium-doped fiber lasers: from all anomalous dispersion to all normal dispersion

SOLITON DYNAMICS IN PASSIVELY MODE-LOCKED FIBER LASERS ZHAO LUMING 2006

Solitons. Nonlinear pulses and beams

Dissipative soliton molecules with independently evolving or flipping phases in mode-locked fiber lasers

Large energy mode locking of an erbium-doped fiber. laser with atomic layer graphene

Group-velocity-locked vector soliton molecules in a birefringence-enhanced fiber laser

Fiber Gratings p. 1 Basic Concepts p. 1 Bragg Diffraction p. 2 Photosensitivity p. 3 Fabrication Techniques p. 4 Single-Beam Internal Technique p.

Supplementary Figure 1: The simulated feedback-defined evolution of the intra-cavity

DISTRIBUTION A: Distribution approved for public release.

HARMONICALLY mode-locked lasers are attractive as

Stable soliton pairs in optical transmission lines and fiber lasers

Regimes of Passive Mode-Locking of Fiber Lasers

Nonlinear dynamics of mode-locking optical fiber ring lasers

Optical solitons and its applications

Group velocity locked vector dissipative solitons in a high repetition rate fiber laser

Nonlinear Transmission of a NOLM with a Variable Wave Retarder inside the Loop

Linear pulse propagation

Multipulse Operation and Limits of the Kerr-Lens Mode-Locking Stability

Geometrical description of the onset of multipulsing in mode-locked laser cavities

Moving fronts for complex Ginzburg-Landau equation with Raman term

OPTI510R: Photonics. Khanh Kieu College of Optical Sciences, University of Arizona Meinel building R.626

Developing spatiotemporal solitons in step-index. multimode fibers

Dissipative soliton interactions inside a fiber laser cavity

Deterministic chaos in an ytterbium-doped mode-locked fiber laser

Statistical light-mode dynamics of multipulse passive mode locking

Pulse-spacing manipulation in a passively mode-locked multipulse fiber laser

Ultrafast Laser Physics

Effect of multiple time delays on intensity fluctuation dynamics in fiber ring lasers

Control of Multiple Solitons in Mode Locked Figure-of-Eight Micro Structured Fiber Lasers to Improve Soliton Stability

Vector solitons with locked and precessing states of polarization

Nonlinearity management: a route to high-energy soliton fiber lasers

Supplementary Information. Temporal tweezing of light through the trapping and manipulation of temporal cavity solitons.

Self-Phase Modulation in Optical Fiber Communications: Good or Bad?

Soliton Molecules. Fedor Mitschke Universität Rostock, Institut für Physik. Benasque, October

Computational Study of Amplitude-to-Phase Conversion in a Modified Unitraveling Carrier Photodetector

Roadmap to ultra-short record high-energy pulses out of laser oscillators

Performance Limits of Delay Lines Based on "Slow" Light. Robert W. Boyd

Fiber lasers and amplifiers. by: Khanh Kieu

Multipulse operation of a Ti:sapphire laser mode locked by an ion-implanted semiconductor saturable-absorber mirror

Strongly enhanced negative dispersion from thermal lensing or other focusing effects in femtosecond laser cavities

Pulsating solitons, chaotic solitons, period doubling, and pulse coexistence in mode-locked lasers: Complex Ginzburg-Landau equation approach

A short tutorial on optical rogue waves

Mode-locking pulse dynamics in a fiber laser with a saturable Bragg reflector

Anticipating Synchronization Based on Optical Injection-Locking in Chaotic Semiconductor Lasers

150 Zhang Sheng-Hai et al Vol. 12 doped fibre, and the two rings are coupled with each other by a coupler C 0. I pa and I pb are the pump intensities

Breakup of Ring Beams Carrying Orbital Angular Momentum in Sodium Vapor

System optimization of a long-range Brillouin-loss-based distributed fiber sensor

Nonlinear Optics (WiSe 2016/17) Lecture 9: December 16, 2016 Continue 9 Optical Parametric Amplifiers and Oscillators

Pulse interaction via gain and loss dynamics in passive mode locking

MEFT / Quantum Optics and Lasers. Suggested problems Set 4 Gonçalo Figueira, spring 2015

Temporal modulation instabilities of counterpropagating waves in a finite dispersive Kerr medium. II. Application to Fabry Perot cavities

Highly Nonlinear Fibers and Their Applications

Orthogonally polarized bright dark pulse pair. generation in mode-locked fiber laser with a large-angle

A tunable corner-pumped Nd:YAG/YAG composite slab CW laser

New Concept of DPSSL

ROLE OF HYSTERESIS ON THE MODE-SHIFT CHARACTERISTICS OF INJECTION LOCKING A LASER DIODE

Final Report for AOARD grant FA Measurement of the third-order nonlinear susceptibility of graphene and its derivatives

Nonlinear pulse shaping and polarization dynamics in mode-locked fiber lasers

1 Mathematical description of ultrashort laser pulses

An Efficient Method to Simulate the Pulse Propagation and Switching Effects of a Fiber Bragg Grating

Feedback-controlled Raman dissipative solitons in a fiber laser

A Guide to Experiments in Quantum Optics

Krauskopf, B., Erzgraber, H., & Lenstra, D. (2006). Dynamics of semiconductor lasers with filtered optical feedback.

Contents Three Sources and Three Component Parts of the Concept of Dissipative Solitons Solitons in Viscous Flows

Ho:YLF pumped HBr laser

Ultrafast Laser Physics

Mode-locked laser pulse fluctuations

Differential Brillouin gain for improving the temperature accuracy and spatial resolution in a long-distance distributed fiber sensor

Optics, Optoelectronics and Photonics

VECTOR SOLITON FIBER LASERS

2. THE RATE EQUATION MODEL 2.1 Laser Rate Equations The laser rate equations can be stated as follows. [23] dn dt

Numerical Analysis of Soft-Aperture Kerr-Lens Mode Locking in Ti:Sapphire Laser Cavities by Using Nonlinear ABCD Matrices

Computer Modelling and Numerical Simulation of the Solid State Diode Pumped Nd 3+ :YAG Laser with Intracavity Saturable Absorber

Drop-port study of microresonator frequency combs: power transfer, spectra and time-domain characterization

Laser Physics OXFORD UNIVERSITY PRESS SIMON HOOKER COLIN WEBB. and. Department of Physics, University of Oxford

Effect of cross-phase modulation on supercontinuum generated in microstructured fibers with sub-30 fs pulses

Sintec Optronics Pte Ltd

Light bullets and dynamic pattern formation in nonlinear dissipative systems

Spatiotemporal coupling in dispersive nonlinear planar waveguides

Effects of self-steepening and self-frequency shifting on short-pulse splitting in dispersive nonlinear media

LIST OF TOPICS BASIC LASER PHYSICS. Preface xiii Units and Notation xv List of Symbols xvii

Stable One-Dimensional Dissipative Solitons in Complex Cubic-Quintic Ginzburg Landau Equation

Applications of Nonlinear Fiber Optics. Second Edition

Enhanced Detection of Gas Absorption Using an Erbium-Doped Fiber Ring Laser

Optoelectronic Applications. Injection Locked Oscillators. Injection Locked Oscillators. Q 2, ω 2. Q 1, ω 1

Graphene-based passively mode-locked bidirectional fiber ring laser

Coherent Combining and Phase Locking of Fiber Lasers

Time resolved optical spectroscopy methods for organic photovoltaics. Enrico Da Como. Department of Physics, University of Bath

Transcription:

PHYSICAL REVIEW A, 66, 033806 2002 Bound-soliton fiber laser D. Y. Tang, B. Zhao, D. Y. Shen, and C. Lu School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore W. S. Man and H. Y. Tam Department of Electrical Engineering, Hong Kong Polytechnic University, Hung Hom, Hong Kong Received 6 May 2002; published 16 September 2002 Experimental study of the soliton dynamics of a passively mode-locked fiber ring laser revealed a state of bound-soliton operation in the laser, where two solitons bind together tightly with fixed pulse separation. We further report on the properties of the bound-soliton emission of the laser. In particular, we demonstrate both experimentally and numerically that, like the single-pulse soliton operation of the laser, the bound soliton emission is another intrinsic feature of the laser. DOI: 10.1103/PhysRevA.66.033806 PACS number s : 42.55.Wd, 42.81.Dp, 42.60.Fc, 42.65.Sf Self-started, passively mode-locked fiber lasers as a potential source of ultrashort optical pulses have been intensively investigated 1 7. A generic feature of the lasers found is that under suitable operation conditions, they can emit the so-called soliton pulses optical pulses that are sech-form shaped and have nearly transform-limited bandwidth-duration product. Under a soliton operation, not only the output pulses of the lasers become ultrashort, but also the pulse-to-pulse energy and peak power become ultrastable, which was found to be quantum noise limited 8. Soliton emission of the lasers is a natural consequence of the nonlinear pulse propagation in the fiber cavity, where due to the balanced action between the fiber-optical Kerr effect and the cavity dispersion on a pulse, the shape and duration remain unchanged with the propagation. Although in a passively mode-locked fiber laser, apart from the optical fiber other optical components such as the gain medium and output coupler coexist in the cavity, which affect the detailed dynamics of the formed pulses. It is nevertheless demonstrated that under weak influence of them, the average dynamics of the solitons could be well described by the nonlinear Schrödinger equation 9. In this paper we report on the states of bound-soliton emission and their properties in a passively mode-locked fiber ring laser. We show both experimentally and numerically that apart from the single-pulse soliton emission, the laser can also emit stable, closely spaced soliton pairs. In particular, when operating in the regime, the bound-soliton pair is the only stable structure of the solitary wave in the laser, as a unit, it has exactly the same features as those of a single-pulse soliton. For this reason, we refer the fiber laser as a bound-soliton fiber laser. Bound states of solitons have recently been predicted in the coupled nonlinear Schrödinger equations 10 and the quintic complex Ginzburg-Landau equation 11,12. Formation of bound solitons was explained as a result of direct soliton interaction. Solitons formed in these systems have an oscillating tail; when they interact, their effective interaction potential has spatial local minima, which give rise to stable, bound solitons. Bound solitons thus formed are characterized as having fixed, discrete pulse separations independent of soliton propagation. For a passively mode-locked fiber laser, when the influences of laser gain and cavity losses and/or fiber birefringence become strong, its dynamics is actually described by the Ginzburg-Landan equation or its coupled forms 13. The passively mode-locked fiber ring laser used in our experiments is schematically shown in Fig. 1. It has a ring cavity of about 5.5 m in length. The cavity is comprised of a 3.5-m-long, 2000 ppm erbium-doped fiber, with a group velocity dispersion of about 10 ps/nm km, a piece of a 1-mlong single-mode dispersion shifted fiber, whose groupvelocity dispersion is 2 ps/nm km, and a piece of a 1-mlong standard single-mode fiber SM-28. The nonlinear polarization rotation technique 5 is used to achieve the selfstarted mode locking in the laser. To this end, a polarization dependent isolator, together with two polarization controllers, one consisting of two quarter-wave plates and the other of two quarter-wave plates and one half-wave plate, is used to adjust the polarization of light in the cavity. The polarization dependent isolator and the polarization controllers are mounted on a 7-cm-long fiber bench to achieve accurate polarization adjustments. The laser is pumped by a pigtailed In x Ga 1 x AsP semiconductor diode of wavelength 1480 nm. The output of the laser is taken via a 10% fiber coupler and analyzed with an optical spectrum analyzer Ando AQ- 6315B and a commercial optical autocorrelator Inrad 5-14- LDA. A 50-GHz-wide bandwidth sampling oscilloscope Agilent 86100A and a 25-GHz high-speed photodetector FIG. 1. A schematic of the passively mode-locked fiber laser. /4 represents quarter-wave plate; /2 represents half-wave plate; PI represents polarization dependent isolator. WDM represents wavelength-dividend multiplexer. 1050-2947/2002/66 3 /033806 6 /$20.00 66 033806-1 2002 The American Physical Society

TANG, ZHAO, SHEN, LU, MAN, AND TAM FIG. 2. Optical spectra of the bound states of the observed solitons. a With spectral modulation period of about 9.1 nm, b 4.7 nm, c 2.6 nm, and d 1.3 nm. The soliton pulse width is about 325 fs and the separations between the bound solitons can be calculated from the corresponding spectral modulation periods. New Focus 1414FC were used to study the evolution of the bound solitons in the cavity. In two previous papers 14,15 we have reported on how to experimentally achieve the bound-soliton operation and some of the basic properties of the observed bound solitons, respectively. The key technique used to get stable boundsoliton operation of the laser is the intracavity mode locking of the dynamical laser modes. With this technique the influences caused by the random environmental noise and the dispersive waves on the soliton interaction could be suppressed. In our previous experiment we have also confirmed that the observed bound solitons have the characteristics of discrete, fixed soliton separations that are independent of soliton propagation, as predicted theoretically. The states of bound solitons are usually abruptly formed in our laser from a single-pulse soliton state by adjusting the orientations of the wave plates. However, if their orientations are already appropriately set, bound-soliton operation can also be directly obtained by simply increasing the pump power beyond the mode-locking threshold. For the completeness of the paper we have shown again in Fig. 2 the optical spectra of the bound solitons observed in our laser. In order to check if the formation of bound states of solitons is a generic feature of the passively mode-locked fiber lasers, in the current experiment we have deliberately changed the cavity property of our laser. Our experimental results confirmed that with exactly the same technique and procedure, the bound states of solitons could be reobtained in the laser. With the new fiber laser setup, we have obtained four discrete states of bound solitons as shown. Comparing them with those obtained with our previous fiber laser setup, except the different pulse separation and soliton pulse width, they have exactly the same properties. Namely, they are ultrastable and have discrete, fixed pulse separations that are neither dependent on the soliton propagation nor on the exact laser operation conditions, such as the pump intensity and the cavity detuning. The change in the soliton pulse width is a natural result of the total cavity dispersion change. In turn it causes the pulse separation change of the bound solitons. Nevertheless, once a laser setup is fixed, the pulse separations of the obtainable bound solitons are fixed. Analyzing the pulse separations among the four bound-soliton states, surprisingly, they again have the relationship that one is roughly twice of the other, as can be easily checked from the spectral modu- 033806-2

BOUND-SOLITON FIBER LASER FIG. 4. Harmonic mode locking of the bound solitons. Eight bound solitons are equally spaced in the cavity forming the special mode-locked state. FIG. 3. Oscilloscope traces of multiple bound solitons in the cavity. The cavity round-trip time is about 26 ns for the laser. a With twelve bound solitons in the cavity, b with eleven bound solitons in the cavity. From a to b the pump laser power is slightly reduced. lation of the bound solitons shown in Fig. 2. This indicates that this relationship could be a universal property of the bound solitons in the fiber laser system, and it is independent of the concrete laser setup. Like the single-pulse soliton operation of the laser, depending on the pumping strength, multiple bound solitons can also coexist in the laser cavity. As an example Fig. 3 shows two oscilloscope traces of such states. Figure 3 a shows that twelve bound-soliton pulses coexist in the laser cavity. Note that the pulse separation between the bound solitons is about 920 fs, and the soliton pulse width of the laser is about 325 fs when a sech-form pulse shape is assumed. As the detector cannot separate the pulse separation in the oscilloscope traces, the bound solitons look like a single pulse. When multiple bound solitons coexist in the laser cavity and they are separated far apart, except that the spectral or signal strength is increased, there is no difference observable on the optical spectrum and the measured autocorrelation traces, compared to those of single bound-soliton pair in the cavity. Based on this experimental evidence together with the result shown in Fig. 3, we conclude that all the bound solitons in the cavity have exactly the same property, namely, the same soliton separation and same pulse energy. There are small peak intensity fluctuations on the pulses shown in the measured oscilloscope traces. We have checked it and determined that they were due to the sampling problem of the sampling oscilloscope. By reducing the recorded time range of the trace, the fluctuation disappeared and all the pulses had exactly the same height. In all of our experiments in a stable state no bound solitons of different pulse separations or different pulse energies have been observed to coexist, and no unbound solitons with the bound solitons, or bound solitons of higher order with more than two soliton pulses binding together have been observed to coexist, which suggests strongly that in the parameter regime, the only stable pulse structure is the bound-soliton pairs of the fixed pulse separation. Carefully increasing or decreasing the pump power, the number of bound-soliton pairs can be increased or reduced one by one. Figures 3 a and 3 b show a case of two consecutive states when the pump power is reduced. Initially there were twelve bound solitons in the cavity. Reducing the pump power to a certain value, one bound-soliton pair was suddenly destroyed. After losing one pair of the bound solitons, the rest of the bound solitons rearranged their relative positions in the cavity and stabilized finally at a new pattern as shown in Fig. 3 b. In our experiments we found that the dispersive waves play an essential role in the long-range interaction among the bound solitons. Mediated through the dispersive waves, the bound solitons can interact among themselves and form various states of multiple bound-soliton operation. In fact, Fig. 3 shows two cases of the so-called quasiharmonic mode locking by some authors in the singlepulse soliton operation of fiber lasers. States of harmonic mode locking of the bound solitons have also been frequently observed. Figure 4 shows as an example one of these states. In the present case, eight bound solitons are equally spaced in the cavity. We have also observed situations where several tens of bound solitons coexist in the cavity and form harmonic mode locking. Depending on the interaction between the bound solitons, bunches of bound solitons have also been revealed in the laser, as shown in Fig. 5. Several bound solitons can tightly couple together locally to form a bound soliton bunch, and the bunch moves with the fundamental repetition rate in the cavity, exactly as also observed in the single-pulse soliton case 1,2,4. When the bound-soliton pairs interact directly with each other, they behave as a unit; either colliding elastically or 033806-3

TANG, ZHAO, SHEN, LU, MAN, AND TAM u z i u u t i 2 2 u t 2 i u 2 2 3 v 2 u i 3 v2 u* g 2 u g 2 g 2 2 u t 2, v v i v z t i 2 2 v t i u 2 2 v 2 2 v 3 i 3 u2 v* g 2 v g 2 g 2 2 v t 2. 1 FIG. 5. Bunching of the bound solitons. Seven bound solitons are coupled tightly together, with no relative movement among the solitons. simply crossing each other without damaging their bindings or altering their pulse separations. With the high-speed sampling oscilloscope and detector, we could monitor in the real time the interaction between two bound-soliton pairs, as shown in Fig. 6. The interaction between the bound-soliton pairs clearly exhibits the particlelike nature of the soliton interaction. Bound states of the bound solitons mediated through the dispersive waves were occasionally observed. However, these bound states are unstable, a small perturbation can destroy it. All these experimental findings demonstrate that the bound-soliton pairs as a unit has exactly all the same properties as those of a single-pulse soliton in the laser. Based on their properties we conclude that the observed bound states of solitons could be another type of solitary wave in the laser and the bound-soliton operation could be an intrinsic feature of the system. To confirm our conclusions, we have also numerically simulated the operation of our laser. For this purpose we have built up a theoretical model based on the coupled extended nonlinear Schrödinger equations. Our model takes into account the birefringence of the optical fiber and the perturbations caused by the laser gain and cavity loss, and has a form FIG. 6. Collision between two bound-soliton pairs. After collision they still remain as bound solitons of the same property. Here u and v are the two normalized slowly varying pulse envelopes along the slow and the fast axes, respectively. 2 2 n/ is the wave-number difference, 2 2 /2 c is the inverse group-velocity difference. is the dispersion parameter and is the nonlinearity of the fiber. g is the laser gain coefficient and g is the bandwidth of the laser gain. Gain saturation of the laser is considered by writing g g 0 1 u 2 v 2 dt E s, 2 where g 0 is the small signal gain and E s is the saturation energy. To also account the effects of other cavity components and the cavity feedback, we simulate the laser by simply allowing the light to circulate in the cavity; whenever it meets a cavity component, we then discretely multiply the component s transformation matrix to the light. In all of our simulations where it is possible, we have used the actual value of the laser parameters. Using the same theoretical model we have previously simulated the dynamics of single-pulse soliton lasing of passively mode-locked fiber lasers, and correctly reproduced the features such as the intrinsic soliton wavelength tuning 16, soliton sideband asymmetry, and subsideband generation 13. In the current simulation we found that by appropriately choosing the linear cavity transmission of the laser and the polarization orientation of the polarization dependent isolator, bound states of solitons could actually be reproduced in our model. In particular, we found that with the parameter settings, bound solitons are the only stable state. Figure 7 shows results of our simulations. Three stable bound-soliton states with different pulse separations have been revealed in our simulations under the laser parameter settings. In our simulations, in order to find out the bound-soliton states we simply fix the laser parameters and input either random initial noise wave form or pulses of different shapes to start our simulations, we then wait for until a stable state is achieved. We found that independent of the initial conditions the laser output always settled down to one of the bound soliton states. Similar to the experimental results, the obtained stable bound solitons exhibit discrete, fixed pulse separations. Once a state of bound solitons is obtained, it will remain there for several thousand rounds of calculation, even under slightly 033806-4

BOUND-SOLITON FIBER LASER FIG. 7. States of bound solitons calculated from the theoretical model. a c show bound solitons with different pulse separations, c d show the corresponding optical spectra. Parameters used are 3 W 1 km 1, 2 ps/nm km for dispersion shifted fiber, 10 ps/nm km for erbium-doped fiber, g 0 254, g 2 10 THz, E s 110, cavity length L 11 m, and beat length L b L/2. changed parameters such as the small signal gain and saturation energy, indicating that it is a stable state of the system. Our numerical simulations have fully confirmed our experimental observations. In summary, we have experimentally studied properties of bound-soliton emission in a passively mode-locked fiber ring laser, and observed features such as the bound-soliton energy quantization, quasibound and harmonically bound soliton mode locking, bound-soliton bunching, and particlelike interactions between bound solitons. Our experimental results have also revealed that no single-pulse soliton can coexist with the bound solitons, and when multiple bound solitons coexist, all bound solitons have exactly the same properties. Exactly the same features have also been observed for the single-pulse soliton in the fiber soliton lasers, which suggests strongly that, like the single-pulse soliton, the observed bound solitons are another fundamental structure of stable light pulse in the laser. The bound-soliton operation of the laser is, in fact, an intrinsic feature of the fiber laser. Finally, we note that numerical simulations on the laser operation have also confirmed the existence of bound solitons in the laser and its property of discrete, fixed soliton separations, which further support our conclusions. W. S. Man and H. Y. Tam acknowledge financial support by The Hong Kong Polytechnic University. 033806-5

TANG, ZHAO, SHEN, LU, MAN, AND TAM 1 D. J. Richardson, R. I. Laming, D. N. Payne, V. J. Matsas, and M. W. Phillips, Electron. Lett. 27, 1451 1991. 2 A. B. Grudinin, D. J. Richardson, and D. N. Payne, Electron. Lett. 27, 1860 1993. 3 C. J. Chen, P. K. A. Wai, and C. R. Menyuk, Opt. Lett. 17, 417 1992. 4 A. B. Grudinin, D. J. Richardson, and D. N. Payne, Electron. Lett. 28, 67 1992. 5 V. J. Matsas, D. J. Richardson, T. P. Newson, and D. N. Payne, Opt. Lett. 18, 358 1993. 6 M. J. Guy, D. U. Noske, and J. R. Taylor, Opt. Lett. 18, 1447 1993. 7 K. Tamura, E. P. Ippen, H. A. Haus, and L. E. Nelson, Opt. Lett. 18, 1080 1993. 8 H. A. Haus and A. Mecozzi, IEEE J. Quantum Electron. 29, 983 1993. 9 S. M. J. Kelly, K. Smith, K. J. Blow, and N. J. Doran, Opt. Lett. 16, 1337 1991. 10 B. A. Malomed, Phys. Rev. A 45, R8321 1992. 11 B. A. Malomed, Phys. Rev. A 44, 6954 1991. 12 N. N. Akhmediev, A. Ankiewicz, and J. M. Soto-Crespo, Phys. Rev. Lett. 79, 4047 1997. 13 D. Y. Tang, S. Fleming, W. S. Man, H. Y. Tam, and M. S. Demokan, J. Opt. Soc. Am. B 18, 1443 2001. 14 D. Y. Tang, W. S. Man, H. Y. Tam, and P. D. Drummond, Phys. Rev. A 64, 033814 2001. 15 D. Y. Tang, B. Zhao, D. Y. Shen, W. S. Man, and H. Y. Tam, unpublished. 16 W. S. Man, H. Y. Tam, M. S. Demokan, P. K. A. Wai, and D. Y. Tang, J. Opt. Soc. Am. B 17, 28 2000. 033806-6