Mathematical ultrashort-pulse laser physics
V.L. Kalashnikov,
Institut fuer Photonik, TU Wien, Gusshausstrasse 27/387, A-1040 Vienna, Austria
vladimir.kalashnikov@tuwien.ac.at
http://www.geocities.com/optomaplev
Abstract: The analytical and numerical approaches to the analysis of the ultrashort-pulse solid-state lasers are presented. The unique self-consistent method of the laser dynamics analysis is based on the symbolical, numerical, programming, and graphical capacities of Maple 6. The algorithmization of basic conceptions as well as sophisticated research methods is of interest to both students and experts in the laser physics and the nonlinear dynamics.
Application Areas/Subjects : Optics, Laser Physics, Nonlinear Physics, Differential Equations, and Programming
Keywords: soliton, ultrashort pulse, mode locking, Q-switching, solid-state laser, nonlinear Schroedinger equation, nonlinear Landau-Ginzburg equation, nonlinear Klein-Gordon equation, harmonic oscillator, nonlinear oscillator, self-phase modulation, group-delay dispersion, self-induced transparency, stimulated Raman scattering
Introduction
The ultrashort laser pulses, i.e. the pulses with the durations ~ - sec, have a lot of the applications, which range ultrafast spectroscopy, tracing chemical reactions, precision processing of materials, optical networks and computing, nuclear fusion and X - ray lasing, ophthalmology and surgery (for review see T. Brabec, F. Krausz , "Intense few-cycle laser fields: Frontiers of nonlinear optics", Rev. Mod. Phys. 72 , 545 (2000)). The mechanisms of the ultrashort pulse generation are active or passive loss switching (so-called Q-switching, part I) and locking of the longitudinal laser modes (part II) due to the active (part IV) or passive ultrafast modulation resulting in the laser quasi-soliton formation. Such quasi-soliton is very similar to the well-known Schroedinger soliton, which runs in the optical networks (part V). As a matter of fact, the model describing active mode locking is based on the usual equation of the harmonic oscillator or its nonlinear modifications, while for the passive mode locking description we have primordially nonlinear Landau-Ginzburg equation (part VI). This equation is the dissipative analog of the nonlinear Schroedinger equation and, as a result of the nonlinear dissipation, there exist a lot of nonstationary regimes of the ultrashort pulse generation (part VII). This requires the generalization of the model, which leads to the numerical simulations on the basis of FORTRAN (or C) codes generated by Maple (part VIII). Simultaneously, the obtained numerical results are supported by the analytical modelling in the framework of the computer algebra approach. The last takes into account the main features of the nonlinear dissipation in the mode-locked laser, viz. the power- (part VI) or energy-dependent (part IX) response of the loss to the generation field. In the latter case, there is the possibility of the so-called self-induced transparency formation, which is described by the nonlinear Klein-Gordon equation (part X).
Our considerations are based on the analytical or semi-analytical search of the steady-state soliton-like solutions of the laser dynamical equations and on the investigation of their stability in the framework of linear stability analysis. Also, the breezer-like solutions are considered using the aberrationless approximation. The computer algebra analysis is supported by the numerical simulations on the basis of the Maple generated FORTRAN-code. We present the analysis of these topics by means of the powerful capacities of Maple 6 in the analytical and numerical computations. This worksheet contains some numerical blocks, which can take about of 12 Mb and 18 min of computation (1 GHz Athlon).
Author is Lise Meitner Fellow at TU Vienna and would like to acknowledge the support of Austrian Science Fund's grant #M611. I am grateful to Dr. I.T. Sorokina and Dr. E. Sorokin for their hospitality at the Photonics Institute (TU Vienna), for stimulating discussions and experimental support of Part VIII. Also, I thank D.O. Krimer for his help in programming of part V.
Contents:
1. Nonstationary lasing: passive Q-switching
2. Conception of mode locking
3. Basic model
4. Active mode locking: harmonic oscillator
5. Nonlinear Schroedinger equation: construction of the soliton solution using the direct Hirota's method
6. Nonlinear Landau-Ginzburg equation: quasi-soliton solution
7. Autooscillations of quasi-solitons in the laser
8. Numerical approaches: ultrashort pulse spectrum, stability and multipulsing
9. Mode locking due to a "slow" saturable absorber
10. Coherent pulses: self-induced transparency in the laser
Part I. Nonstationary lasing: passive Q-switching
Continuous-wave oscillation
The basic principle of Q-switching is rather simple, but in the beginning let's consider the steady-state oscillation of laser. The near-steady-state laser containing an active medium and pumped by an external source of the energy (lamp, other laser or diode, for example) obeys the following coupled equations:
Warning, the protected names norm and trace have been redefined and unprotected
Here is the time-dependent field intensity, is the dimensionless gain coefficient, P is the time-independent (for simplicity sake) pump intensity, and are the frequencies of the pump and generation fields, respectively, and are the absorption and generation cross-sections, respectively, is the gain relaxation time, is the linear loss inclusive the output loss of the laser cavity, and, at last, is the gain coefficient for the full population inversion in the active medium. The pump increases the gain coefficient (first term in eq2 ), that results in the laser field growth (first term in eq1 ). But the latter causes the gain saturation (second term), which can result in the steady-state operation (so-called continuous-wave, or simply cw, oscillation):
The second solution defines the cw intensity, which is the linear function of pump intensity:
The pump corresponding to =0 defines so-called generation threshold. Now let's consider the character of the steady-state points of our system { eq1 , eq2 } presented by sol . The Jacobian of the system { eq1 , eq2 } is:
For the cw-solution the eigenvalues of the perturbations are:
So, cw oscillations is stable in our simple case because the eigenvalues are negative. For the zero field solution the perturbation eigenvalues are:
The existence of the positive eigenvalue suggests the instability of the zero-field steady-state solution. Hence there is the spontaneous generation of the cw oscillation above threshold in the model under consideration.
Q-switching
The situation changes radically due to insertion of the saturable absorber into laser cavity. In this case, in addition to the gain saturation, the loss saturation appears. This breaks the steady-state operation and produces the short laser pulses. As a result of the additional absorption, Q-factor of laser is comparatively low (high threshold). This suppresses the generation. When is small, the gain increases in the absence of the gain saturation (see eq2 from the previous subsection). This causes the field growth. The last saturates the absorption and abruptly increases Q-factor. The absorption "switching off" leads to the explosive generation, when the most part of the energy, which is stored in the active medium during pumping process, converts into laser field. The increased field saturates the gain and this finishes the generation.
As the reference for the model in question see, for example, J.J. Degnan, "Theory of the optimally coupled Q-switched laser", IEEE J. Quant. Elect. 25 , 214 (1989). To formulate the quantitative model of the laser pulse formation let's use the next approximations: 1) the pulse width is much larger than the cavity period, and 2) is less than the relaxation time, 3) the pump action during the stage of the pulse generation is negligible. We shall use the quasi-two level schemes for the gain and loss media (the relaxation from the intermediate levels is fast). Also, the excited-state absorption in absorber will be taken into account.
The system of equation describing the evolution of the photon density is
Warning, the name changecoords has been redefined
Now we shall search the ground state population in the absorber as a function of the initial population inversion in amplifier:
The similar manipulation allows to find the photon density as a function of inversion in amplifier:
Hence the photon density is:
So, we have:
( x )= ( Eq. 1 )
Now let's define the key Q-switching parameters:
So, the initial inversion defining the gain at Q-switching start is
( Eq. 2 )
Additionally, we define the inversion at the pulse maximum, when tends to infinity:
So, we have the expressions for (initial inversion, sol_3 ), (the inversion at pulse maximum, sol_4 and e5 ), (the inversion at pulse maximum when --> , sol_6 ), (the final inversion, sol_5 and e6 ) and the photon density as function of inversion x ( sol_2 ).
As an example, we consider the real situation of Yb/Er-glass laser with the crystalline Co:MALO saturable absorber. The obtained expressions allow to plot the typical dependencies for the pulse parameters:
This numerical procedure plots the dependence of the output pulse energy on the reflectivity of the output mirror:
Warning, `delta` is implicitly declared local to procedure `En`
Warning, `sigma` is implicitly declared local to procedure `En`
Warning, `alpha` is implicitly declared local to procedure `En`
Warning, `h` is implicitly declared local to procedure `En`
Warning, `nu` is implicitly declared local to procedure `En`
Warning, `S` is implicitly declared local to procedure `En`
Warning, `R` is implicitly declared local to procedure `En`
Warning, `x_i` is implicitly declared local to procedure `En`
Warning, `eq` is implicitly declared local to procedure `En`
Warning, `sol_f` is implicitly declared local to procedure `En`
Warning, `sol_En` is implicitly declared local to procedure `En`
For the comparison we use the experimental data (crosses in Figure):
Similarly, for the output power we have:
Warning, `delta` is implicitly declared local to procedure `Pow`
Warning, `sigma` is implicitly declared local to procedure `Pow`
Warning, `alpha` is implicitly declared local to procedure `Pow`
Warning, `h` is implicitly declared local to procedure `Pow`
Warning, `nu` is implicitly declared local to procedure `Pow`
Warning, `S` is implicitly declared local to procedure `Pow`
Warning, `R` is implicitly declared local to procedure `Pow`
Warning, `c` is implicitly declared local to procedure `Pow`
Warning, `refractivity` is implicitly declared local to procedure `Pow`
Warning, `t_cav` is implicitly declared local to procedure `Pow`
Warning, `x_i` is implicitly declared local to procedure `Pow`
Warning, `x_t0` is implicitly declared local to procedure `Pow`
Warning, `eq` is implicitly declared local to procedure `Pow`
Warning, `sol_t` is implicitly declared local to procedure `Pow`
Warning, `sol_Pow` is implicitly declared local to procedure `Pow`
And, at last, the pulse durations are:
Warning, `delta` is implicitly declared local to procedure `Width`
Warning, `sigma` is implicitly declared local to procedure `Width`
Warning, `alpha` is implicitly declared local to procedure `Width`
Warning, `h` is implicitly declared local to procedure `Width`
Warning, `nu` is implicitly declared local to procedure `Width`
Warning, `S` is implicitly declared local to procedure `Width`
Warning, `R` is implicitly declared local to procedure `Width`
Warning, `c` is implicitly declared local to procedure `Width`
Warning, `refractivity` is implicitly declared local to procedure `Width`
Warning, `t_cav` is implicitly declared local to procedure `Width`
Warning, `x_i` is implicitly declared local to procedure `Width`
Warning, `x_t0` is implicitly declared local to procedure `Width`
Warning, `eq1` is implicitly declared local to procedure `Width`
Warning, `eq2` is implicitly declared local to procedure `Width`
Warning, `sol_f` is implicitly declared local to procedure `Width`
Warning, `sol_t` is implicitly declared local to procedure `Width`
Warning, `sol_Width` is implicitly declared local to procedure `Width`
The worse agreement with the experimental data for the pulse durations is caused by the deviation of the pulse shape from the Gaussian profile, which was used for the analytical estimations. More precise consideration is presented in J. Liu, D. Shen, S.-Ch. Tam, and Y.-L. Lam, "Modelling pulse shape of Q-switched lasers", IEEE J. Quantum Electr. 37 , 888 (2001).
The main advantage of the analytical model in question is the potential of the Q-switched laser optimization without any cumbersome numerical simulations. Let's slightly transform the above obtained expressions:
Now, let's plot some typical dependence for Q-switching.
Warning, `sol` is implicitly declared local to procedure `xt`
Warning, `x` is implicitly declared local to procedure `xt`
From the definition of this dependence for growing tends to the linear one, which correspond to the maximal efficiency of the population inversion utilization. For the fixed absorption and emission cross-sections, the parameter increases as a result of the laser beam focusing in the saturable absorber. Then ( --> ), we have:
and the output energy optimization can be realized by this simple way:
From the first figure the optimization can be performed by the simple graphical way. We have to define the appropriate for our scheme laser's net-loss and to determine (measure or calculate) the intracavity linear loss. This gives the value of b , which is the relation of the linear loss to the net-loss. The upper group of curves gives the value of c , the lower curves give a (for the different contribution of the excited-state absorption, i. e. ).
Two-color pulsing
Now let consider a more complicated situation, which corresponds to the two-color Q-switching due to presence of the stimulated Raman scattering in the active medium (for example, Yb(3+):KGd(W0(4))(2), see A.A. Lagatsky, A. Abdolvand, N.V. Kuleshov, Opt. Lett. 25 , 616 (2000)). In this case the analytical modelling is not possible, but we can use the numerical capacities of Maple.
Let's the gain medium length is and the Raman gain coefficient is g . We assume the exact Raman resonance and neglect the phase and group-velocity effects. Then the evolutional equations for the laser and scattered intensities and , respectively, are:
The integration produces:
There are the amplification of the scattered field and the depletion of the laser field, which plays a role of the pump for the Stokes component.
Now let's transit from the intensities to the photon densities:
Hence the photon densities evolution due to the stimulated Raman scattering obeys:
The inverse and ground-state populations for the gain medium and the absorber (Cr(4+):YAG) obey (see previous subsection):
Contribution of gain, saturable, output and linear intracavity loss to the field's evolution results in:
where and are the linear loss for the laser and Stokes components, respectively, and are the output mirror reflectivity at the laser and Stokes wavelengths, respectively, is the reduction factor taking into account the decrease of the loss cross-section at the Stokes wavelength relatively to the lasing one.
Hence, the field densities evaluate by virtue of:
In the agreement with the results of the previous subsection we can transform the equation for the populations in the gain and absorption media:
Let' introduce new variables:
Then from Eqs. ( eq5, eq10, eq11, eq12 ):
The next procedures solve the obtained system numerically to obtain the dependences of the normalized photon densities at laser and Stokes wavelengths and the relative inversion vs. cavity roundtrip, n is the duration of the simulation in the cavity roundtrips):
Warning, `gam` is implicitly declared local to procedure `ODE_plot1`
Warning, `rho` is implicitly declared local to procedure `ODE_plot1`
Warning, `kappa` is implicitly declared local to procedure `ODE_plot1`
Warning, `alpha` is implicitly declared local to procedure `ODE_plot1`
Warning, `x` is implicitly declared local to procedure `ODE_plot1`
Warning, `N` is implicitly declared local to procedure `ODE_plot1`
Warning, `G` is implicitly declared local to procedure `ODE_plot1`
Warning, `Xi` is implicitly declared local to procedure `ODE_plot1`
Warning, `sys` is implicitly declared local to procedure `ODE_plot1`
Warning, `gam` is implicitly declared local to procedure `ODE_plot2`
Warning, `rho` is implicitly declared local to procedure `ODE_plot2`
Warning, `kappa` is implicitly declared local to procedure `ODE_plot2`
Warning, `alpha` is implicitly declared local to procedure `ODE_plot2`
Warning, `x` is implicitly declared local to procedure `ODE_plot2`
Warning, `N` is implicitly declared local to procedure `ODE_plot2`
Warning, `G` is implicitly declared local to procedure `ODE_plot2`
Warning, `Xi` is implicitly declared local to procedure `ODE_plot2`
Warning, `sys` is implicitly declared local to procedure `ODE_plot2`
Warning, `gam` is implicitly declared local to procedure `ODE_plot3`
Warning, `rho` is implicitly declared local to procedure `ODE_plot3`
Warning, `kappa` is implicitly declared local to procedure `ODE_plot3`
Warning, `alpha` is implicitly declared local to procedure `ODE_plot3`
Warning, `x` is implicitly declared local to procedure `ODE_plot3`
Warning, `N` is implicitly declared local to procedure `ODE_plot3`
Warning, `G` is implicitly declared local to procedure `ODE_plot3`
Warning, `Xi` is implicitly declared local to procedure `ODE_plot3`
Warning, `sys` is implicitly declared local to procedure `ODE_plot3`
So, we can obtain a quite efficient two-color pulsing in the nanosecond time domain without any additional wavelength conversion.
The nonstationary lasing in question produces the pulses with the durations, which are larger than the cavity period. The cavity length shortening, i. e. the use of the microchip lasers, decreases the pulse durations down to ten-hundred picoseconds. But there is a method allowing the fundamental pulse width reduction, viz. mode locking.
Part II. Conception of mode locking
The laser cavity is, in fact, interferometer, which supports the propagation of only defined light waves. Let consider a plane wave, which is reflected from a laser mirror. The initial wave is ( cc is the complex conjugated term):
Then the reflected wave (normal incidence and full reflectivity are supposed) is:
where is the phase shift due to reflection. An interference between incident and reflected waves results in
This is the so-called standing wave:
Note, that a wave node lies on a surface (point x =0). The similar situation takes place in the laser resonator. But the laser resonator consists of two (or more) mirrors and the standing wave is formed due to reflection from the each mirror. So, the wave in the resonator is the standing wave with the nodes placed on the mirrors. Such waves are called as the longitudinal laser modes . The laser resonator can contain a lot of modes with the different frequencies (but its nodes have to lie on the mirrors!) and these modes can interfere.
Let suppose, that the longitudinal modes are numbered by the index m . In fact, we have M harmonic oscillators with the phase and frequency differences dphi and domega , correspondingly. Let the amplitude of modes is A0 .
Here 0 and 0 are the phase and frequency of the central mode, respectively.
The interference between these modes produces the wave packet:
Now, we can extract the term describing the fast oscillation on the central ("carrier") frequency 0 from the previous expression. The obtained result is the packet's envelope (its slowly varying amplitude):
It is obviously, that this expression can be converted into following form:
The squared envelope's amplitude (i. e. a field intensity) is depicted in the next figure for the different M .
One can see, that the interference of modes results in the generation of short pulses. The interval between pulses is equal to 2 / domega . The growth of M decreases the pulse duration 2 /( M*domega ) and to increases the pulse intensity * . The last is the consequence of the following relation:
In this example, the phase difference between neighboring modes is constant. Such mode locking causes the generation of short and intense pulses. But in the reality, the laser modes are not locked, i. e. the modes are the oscillations with the independent and accidental phases. In this case:
Thus, the interference of the unlocked modes produces the irregular field beatings, i. e. the noise spikes with a duration ~ 1/(M*domega) .
What are the methods for the mode locking? Firstly, let consider the simplest model of the harmonic oscillation in the presence of the periodical force.
We can see, that the external force causes the oscillations with the additional frequencies: + and - . If is equal to the intermode interval, the additional oscillation of mode plays a role of the resonance external force for the neighboring modes. Let consider such resonant oscillation in the presence of the resonant external force:
The term, which is proportional to t ("secular" term), equalizes the phase of the oscillations to the phase of the external force. It is the simplest model of a so-called active mode locking . Here the role of the external force can be played by the external amplitude or phase modulator. Main condition for this modulator is the equality of the modulation frequency to the intermode interval that causes the resonant interaction between modes and, as consequence, the mode locking (part 3 ).
The different mechanism, a passive mode locking , is produced by the nonlinear interaction of modes with an optical medium. Such nonlinearity can be caused by saturable absorption, self-focusing etc. (see further parts of worksheet). Now we shall consider the simplest model of the passive mode locking. Let suppose, that there are two modes, which oscillate with different phases and frequencies in the cubic nonlinear medium:
We can see, that the initial oscillations with the different phases are locked due to nonlinear interaction that produces the synchronous oscillations.
As conclusion, we note that the mode locking is resulted from the interference between standing waves with constant and locked phases. Such interference forms a train of the ultrashort pulses. The mechanisms of the mode locking are the external modulation with frequency, which is equal to intermode frequency interval, or the nonlinear response of the optical medium. Later on we shall consider both methods. But firstly we have to obtain the more realistic equations describing the ultrashort pulse generation.
Part III. Basic model
The models describing the laser field evolution are based usually on the so-called semi-classical approximation. In the framework of this approximation the field obeys the classical Maxwell equation and the medium evolution has the quantum-mechanical character. Here we shall consider the wave equation without concretization of the medium evolution.
The Maxwell equation for the light wave can be written as:
where E(z,t) is the field strength, P(t) is the medium polarization, z is the longitudinal coordinate, t is the time, c is the light velocity. The change of the variables z --> z* , t - z/c --> t* produces
We does not take into account the effects connected with a wave propagation in thin medium layer, that allows to eliminate the second-order derivation on z* . Then
So, we reduced the order of wave equation. The inverse transformation of the coordinates leads to the so-called shortened wave equation:
Next step is the transition to the slowly-varying amplitude approximation. We shall consider field envelope (z,t) and polarization P(t), which are filled by the fast oscillation with frequency ( k is the wave number, N is the concentration of the atoms, d is the medium length):
Then the wave equation can be transformed as:
The obtained result is the system of the shortened wave equation and the dispersion condition. The right-hand side of the field equation (material part) will be different for the different applications (see below).
Part IV. Active mode locking: harmonic oscillator
Amplitude modulation
We start our consideration with a relatively simple technique named as active mode locking due to amplitude modulation. The modulator, which is governed by external signal and changes the intracavity loss periodically, plays the role of the external "force" (see part 1). Let consider the situation, when the modulation period is equal and the pulse duration is much less than the cavity period. If the modulation curve is close to cosine, then the master equation for the ultrashort pulse evolution can be written as [ D. Kuizenga, A. Siegman,"FM and AM mode locking of homogeneous laser", IEEE J. Quant. Electr., 6 , 694 (1970) ]:
= g + ,
where g is the net-gain in the laser accounting for the saturated gain and linear loss l (including output loss), is the inverse bandwidth of the spectral filter, M is the characteristic of the modulation strength taking into account curvature of the modulation curve at the point of maximal loss.
Let try to solve this equation in the case of the steady-state propagation of ultrashort pulse (when =0) and in the absence of detuning of modulation period from cavity period. If the time is normalized to , then the obtained equation is a well-known equation of harmonic oscillator:
The next step is suggested by the asymptotic behavior of the prospective solution: = 0. But previously it is convenient to transit from Whittaker functions to hypergeometric and Kummer functions:
WhittakerM( , ,z)= hypergeom( , 1+2 , x ),
WhittakerW( , ,z)= KummerU( , 1+2 , x ).
As result, we have ( = , = , x = ):
(t)= t hypergeom( , , ) +
t KummerU( , , )=
t hypergeom(( )+ , , ) +
hypergeom( , , ).
Now, the asymptotic condition = 0, which is similar to condition for quantum states in harmonic oscillator, results in (see, for example, S. Flugge, Practical Quantum Mechanics I, Springer-Verlag (1971) )
(t) = C ( t ) for n= - ( ) && n is integer,
where and are Hermite polynomials. Value of constant C can be obtained from the energy balance condition, which results from the equation of gain saturation:
= .
Here is the inverse flux of the gain saturation energy, is the gain for small signal defined by gain medium properties and pump intensity (note, that g = ).
Now let investigate the parameters of the steady-state solution of the master equation.
With that end in view, we have to search the generation field energy:
The use of normalization of intensity to
and energy balance condition (see above) results in
Hence the ultrashort pulse obtained as result of active mode locking can be represented as
Now we calculate the pulse duration measured on half-level of maximal pulse intensity:
One can see, that the increase of modulation parameter decreases the pulse width, but this decrease is slow (~1/ ).
Next step is the taking into account the detuning of cavity and modulation periods. The corresponding normalized parameter can be introduced in the following form (see corresponding Doppler transformation: t-->t - z and, as consequence, -->- for steady-state pulse, here is, in fact, inverse relative velocity with dimension [time/cavity transit], i. e. time delay on the cavity round-trip):
The comparison with above obtained result leads to
and we can repeat our previous analysis
Note, that there is the maximal permitting the "ground state" ultrashort pulse generation:
We see, that the upper permitted level of detuning parameter is increased due to pump growth (rise of ) and is decreased by M growth.
It is of interest, that the growth of does not leads to generation of the "excite state" pulses because of the corresponding limitation for these solutions is more strict:
The dependence of pulse width on has following form:
We can see almost linear and symmetric rise of pulse width due to detuning increase. This is an important characteristic of active mode-locked lasers. But in practice (see below), the dependence of the pulse parameters on detuning has more complicated character.
The dependence of the pulse maximum location on detuning parameter can be obtained from the solution for pulse envelope's maximum: = 0. Hence, the pulse maximum location is:
The increase (decrease) of the pulse round-trip frequency ( <0 and >0, respectively) increases positive (negative) time shift of the pulse maximum relatively modulation curve extremum.
Phase modulation
The external modulation can change not only field amplitude, but its phase. In fact, this regime (phase active modulation) causes the Doppler frequency shift of all field components with the exception of those, which are located in the vicinity of extremum of modulation curve. Hence, there exists the steady-state generation only for field located in vicinity of points, where the phase is stationary. Steady-state regime is described by equation
where is the phase delay on the cavity round-trip.
This equation looks like previous one, but has complex character. This suggests to search its partial solution in the form ( t ) = :
So, we have ( t ) = and = , g= . The time-dependent parabolic phase of ultrashort pulse is called chirp and is caused by phase modulation. Pulse width for this pulse is:
that is 2 . As one can see, that result is equal to one for amplitude modulation. The main difference is the appearance of the chirp.
Let take into account the modulation detuning. In this case we may to suppose the modification of the steady-state solution:
We see, that there is not frequency shift of the pulse ( d= 0), but, as it was for amplitude modulation, the time delay appears ( ) that changes the pulse duration as result of modulation detuning. The rise of the detuning prevents from the pulse generation due to saturated net-gain coefficient increase ( ). The pulse parameters behavior coincides with one for amplitude modulation.
Ultrashort pulse stability
Now we shall investigate the ultrashort pulse stability against low perturbation ( t ). The substitution of the perturbed steady-state solution in dynamical equation with subsequent linearization on ( t ) results in
= + - l + - ,
where is the perturbed saturated gain, which is obtained from the assumption about small contribution of perturbation to gain saturation process:
Here A= , B= and we neglected the high-order terms relatively perturbation amplitude. Note, that this is negative quantity.
Let the dependence of perturbation on z has exponential form with increment . Then
Now we introduce * where B*= , that allows to eliminate the exponent from right-hand side of (we shall eliminate the asterix below).
This is equation for eigenvalues and eigenfunctions of the perturbed laser operator. Stable generation of the ultrashort pulse corresponds to decaying of perturbations, i.e. <0. For Gaussian pulse:
As result (see above), we have
(t) = + C ( t ) for n= - ( ) && n is integer,
The glance on the solution is evidence of absence of solution corresponding to . For others modes we take into account, that = l + . As result, all perturbation modes are unstable because of
< 0 ( n= 1, 2, ...),
< 0 ( n= 0, 1, ...).
It should be noted, that for the amplified pulse
,
Here we see the decrease of increment as result of n rise. Therefore, only "ground state" with will be amplified predominantly. This fact provides for Gaussian pulse generation (see, F.X. Kartner, D.M. Zumbuhl, and N. Matuschek, "Turbulence in Mode-Locked Lasers", Phys. Rev. Lett. 82 , 4428 (1999)).
In the presence of detuning we have:
It is surprisingly, but, as it was above, our linear analysis predicts the pulse stability regardless of detuning (because of = l + ). But for the pulse the detuning growth decreases the increment that does not favor the single pulse generation.
Nonlinear processes: self-phase modulation and dynamical gain saturation
Among above considered effects only gain saturation by full pulse energy can be considered as nonlinear process, which, however, does not affect on the pulse envelope, but governs its energy. The time-dependent nonlinear effects, which can transform pulse profile, are self-phase modulation (SPM) and dynamical gain saturation. First one is the dependence of the field phase on its intensity and can play essential role in solid-state lasers [ V.L. Kalashnikov, I.G. Poloyko, V.P. Mikhailov, "Phase modulation of radiation of solid-state lasers in the presence of Kerr optical nonlinearity", Optics and Spectriscopy, 84 , 104 (1998) ]. Second effect is caused by the change of the gain along pulse profile and is essential in lasers with large gain cross-sections and comparatively narrow gain band [ V.L. Kalashnikov, I.G. Poloyko, V.P. Mikhailov, "Generation of ultrashort pulses in lasers with external frequency modulation", Quantum Electronics, 28 , 264 (1998) ].
At first, let analyze the presence of SPM, which can be considered as perturbation of our master equation describing active phase modulation:
Hence we have solution with perturbed parameters:
( t ) = and g= .
As one can see, this solution has enlarged chirp and reduced pulse duration due to SPM.
Now we shall investigate the influence of dynamical gain saturation on the pulse characteristics in the case of active mode locking due to amplitude modulation. Let the contribution of the dynamical gain saturation can be considered as perturbation for Gaussian pulse (see above). The instant energy flux of such pulse is:
The approximation of the small contribution of the gain saturation allows the expansion of energy into series on t up to second order:
Hence we have the modified gain coefficient ( t ), where is the gain coefficient at pulse peak, is the pulse intensity for unperturbed solution. Note, that the additional term in brackets is resulted from the shift of pulse maximum.
Then the master equation for perturbed solution (in our approximation!) is
The comparison with above obtained result gives
(t) = C ( t ) for n= - ( ) && n is integer.
In future we shall omit -terms. If we take the unperturbed intensity for calculation of perturbation action, the perturbed pulse energy is
Now we plot the dependencies of the pulse intensity, width and maximum location versus detuning parameter .
We see, that the main peculiarity here is the asymmetric dependence of the pulse parameters on . The pulse width minimum and intensity maximum don't coincide with =0 and the detuning characteristics have sharper behavior in negative domain of detuning.
Now, as it was made in previous subsection, we estimate the condition of the ultrashort pulse stability.
Here we take into consideration -term, but this does not fail our analysis because of this term contribute only to pulse energy without shift pulse inside modulation window.
For the sake of the simplification, we shall consider the contribution of the destabilizing field to dynamical gain saturation, but only in the form of the unperturbed peak intensity variation and perturbation of the saturated gain coefficient (see above). Then the stability condition:
We can see, that the perturbation rise can destabilize the pulse as result of increase (compare with subsection "Amplitude modulation" ). Also, there is the possibility of ultrashort pulse destabilization near =0. So, the presence of dynamical gain saturation gives the behavior of the ultrashort pulse parameters and stability condition, which is close to the experimentally observed and numerically obtained (see, for example, J.M. Dudley, C.M. Loh and J.D. Harvey, Stable and unstable operation of a mode-locked argon laser, Quantum and Semiclass. Opt. , 8 , 1029 ( 1996 )).
Now we try to investigate the influence of dynamical gain saturation in detail by using so-called aberrationless approximation . Let the pulse profile keeps its form with accuracy up to n -order of time-series expansion, but n pulse parameters are modified as result of pulse propagation. Then the substitution of the expression for pulse profile in master equation with subsequent expansion in t -series produces the system of n ODE describing the evolution of pulse parameters.
Now try to find the steady-state points of ODE-system, which correspond to stationary pulse parameters. For this aim let introduce substitution :
Hence , ( x is the pulse intensity). So, the shift is that differs from the usual result (see above) as result of dynamical gain saturation (last term), which shifts the pulse maximum in negative side. This additional shift has obvious explanation. The gain at the pulse front is greater than one at pulse tail due to gain saturation. This shifts the pulse forward as hole .
Pulse width is:
We can see, that there is the minimum of the pulse duration in negative domain of that corresponds to result, which was obtained on the basis of perturbation theory. The pulse intensity :
So, we have the following dependencies for the pulse duration
and pulse intensity
The last step is the stability analysis. The stability of our solutions can be estimated from the eigenvalues of Jacobian of sol .
Now we find the eigenvalues of Jacobian directly by calculation of determinant.
One can see that the pulse with Gaussian-like form is stable in the region of its existence (see two previous figures). We have to note that the considered here perturbations belong to limited class therefore this criterion is necessary but not sufficient condition of pulse stability (compare with previous consideration on the basis of perturbation theory, where we analyzed not only pulse peak variation but also gain coefficient change).
Part V. Nonlinear Schroedinger equation: construction of the soliton solution using the direct Hirota's method
The Schroedinger equation is the well-known nonlinear equation describing the weak nonlinear waves, in the particular, the laser pulse propagation in fibers. In the last case, a pulse can propagate without decaying over large distance due to balance between two factors: SPM and group delay dispersion (GDD). These pulses are named as optical solitons [ M. J. Ablowitz, H. Segur, "Solitons and the Inverse Scattering Transform", SIAM Philadelphia, 1980 ]. The ultrashort pulse evolution obeys to the next master equation:
which is the consequence of eq_field from part 2 in the case of transition to local time t-->t-z/c . The right-hand side terms describe GDD (with coefficient ) and SPM (with coefficient ), respectively .
It is very important to obtain the exact soliton solutions of nonlinear equations. There are the inverse and direct methods to obtain such solutions. One of the direct methods is the so-called Hirota's method. The main steps of this method are: 1) the selection of the suitable substitution instead of the function (see the master equation), that allows to obtain the bilinear form of the evolution equation; 2) the consideration of the formal series of perturbation theory for this bilinear equation. In the case of soliton solutions these series are terminated.
The useful substitution for the nonlinear Schroedinger equation is (z,t)=G(z,t)/F(z,t). Let suppose that F is the real function. It should be noted that we can make any assumption about to satisfy the assumptions 1) and 2). Hirota proposed to introduce a new D - operator in following way:
After substitution of in the terms of functions G and F we obtain two bilinear differential equations with regard to the new operator D :
(1)
The functions G and F can be expanded into the series of the formal parameter :
Let substitute G and F into Eq. (1) and treat the terms with powers of as independent, to get the infinite set of the differential equations relatively G1, G3, ...; F2, F4, ... . These formal series are terminated only in the case when the master equation has exact N -soliton solution. For instance, the set of first six differential equations in our case is:
For sake of the simplification of the very cumbersome manipulations we introduce the procedure for operator , which acts on the functions a and b . The lasts are the exponents (or linear combination of the exponents) in the form , where is linear function.
Procedure
The next procedure will be used for calculation of the derivative of (or combination of exponents) on t or z with further simplification of the obtained expression.
Procedure der
The next procedure is used to calculate an integral of (or combination of exponents) on t or z with further simplification of the expression.
Procedure Integr
Now, let try to obtain a first-order soliton for nonlinear Schroedinger equation.
To obtain this result we use the trivial relationships:
As was shown above a= - i , hence the last term in third equation of set is equal to 0 . So, we are to choose G3 = 0 to satisfy third equation. Furthermore F2 from fourth equation of the set is ( ) exp( + _s). But in concordance with above obtained relationships this expression is equal to zero. So we can choose F4 = 0. Thus to satisfy other equations we can keep in the expansion of functions G and F only G1 , G3 and F2 .
So, the formal series are terminated. Since is independent parameter we can take =1.
All right! This is the exact solution of the Schroedinger equation. Physically b1 has a sense of the inverse pulse duration. So it is real parameter. But what is the free parameter 10 ? Let 10 is real and = b1. Then