It has been recently understood that mode locking of lasers has the signification of a thermodynamic phase transition in a system of many interacting light modes subject to noise. In the same framework, self starting of passive mode locking has the thermodynamic significance of a noise-activated escape process across an entropic barrier. Here we present the first dynamical study of the light mode system. While accordant with the predictions of some earlier theories, it is the first to give precise quantitative predictions for the distribution of self-start times, in closed form expressions, resolving the long standing self starting problem. Numerical simulations corroborate these results, which are also in good agreement with experiments.
©2006 Optical Society of America
Passive mode locking is a nonlinear phenomenon: It is a method where nonlinearity in the form of a saturable absorber is deliberately introduced into a laser cavity, and this nonlinearity triggers a pulsed regime in lasers. The idea is simple: an intensity fluctuation upon a uniform intensity profile will experience less loss than its surroundings, and will further grow until a pulse is formed.
However, not always does the above described process self start, and often, in order for the laser to self start, it has to operate at a sufficiently high intracavity power. Although it is plausible that a nonlinear phenomenon such as mode locking requires sufficient power, naively the equations of motion of a waveform repeatedly passing through a saturable absorber predict that pulses should always appear, only that the buildup time would scale like the inverse of the laser power, without any threshold behavior.
This puzzle has intrigued many authors in the electro-optic community, and manyworks were dedicated to explaining the power threshold behavior of passive mode locking and deriving a condition for self starting of pulsation. Ippen et al.  and later Chen et al.  suggested dynamic gain saturation as the source of the threshold behavior. Haus and Ippen  suggested reflections in the cavity to be the reason, since they introduce random dispersion to the cavity modes, and Krausz et. al. proposed  that a decoherence process in the cavity tends to disorder the phases of cavity modes, thus opposing mode locking. These theories have been further studied in a series of works [5, 6, 7, 8, 9], and other models have also been suggested [10, 11, 12].
Many of the above mentioned theories associate the threshold behavior with some sort of noise, randomness and decoherence. However the first stochastic theory of the onset of passive mode locking was presented only recently [15, 16, 17, 18, 19]. Statistical mechanics then turned out to be a very convenient tool for studying mode locking with noise, and the onset of pulsing in passively mode locked lasers has been described as a first order phase transition which occurs under variation of the intracavity power or of the noise power. The latter plays a role analogous to temperature in equilibrium statistical mechanics.
The resulting statistical-mechanics theory, which has been termed statistical lightmode dynamics (SLD), predicts that for high noise power or low intracavity power, cw is the stable state (in the thermodynamic sense) of a passively mode locked laser. When noise is decreased or laser power increased sufficiently, cw becomes metastable and the mode-locked pulsed state becomes stable, i. e., the ultimate stationary state of the laser is mode locked. However, such a metastable state can be very long lived, and the cw operation of the laser may therefore persist, even though the stable state is the mode-locked one, much as a supercooled liquid may stay unfrozen before an appropriate fluctuation drives it to the solid phase. The resulting hysteresis behavior is typical in first order phase transitions.
The self-starting problem stems therefore from the fact that the laser may remain trapped in a metastable cw state. It then needs a perturbation to drive it to the stable pulsed state, as is widely agreed. The perturbation may be provided by an external “morning wake-up kick”, but self-starting will occur only if a perturbation is provided by an internal source, i. e., by the intrinsic noise, which is also responsible for the presence of the entropic barrier which traps the laser in cw in the first place. In Fig. 1 we present a self-starting event recorded from a direct numerical simulation of the model studied in this work, see Eq. (1).
The self-starting process is therefore the noise-activated crossing of an entropic barrier. Since activation rates decrease exponentially with the strength of the barrier, the idea that there is a self-starting “threshold” has been put forward. In fact, when the entropic barrier is strong compared with the noise power, the self-starting process is a Poisson process, and as a consequence, the distribution of self-start times is exponential, a universal property of noise-activated escape processes .
In this paper we present a detailed study of the self-starting dynamics of passive mode locking, and obtain a quantitative theory of the process where a system initially in a metastable quasi-cw configuration reaches a steady state. We focus on the the parameter region where self starting occurs by a rare event of barrier crossing and the mean self-start time is considerably longer than the dynamical pulse buildup timescale. In the opposite case, the self-start process and rate are determined by the nonlinear dynamics of the saturable absorber and noise plays a negligible role. Since pulse buildup times are typically shorter than 1 msec, there is a broad class of passively mode locked lasers where self-starting is slowed down by a large factor because of the entropic barrier but can still be observed practically. The present theory applies to all such self starting laser systems.
The main result of this work is an explicit expression, Eq. (13), for the self-start time in a model of passive mode locking with a finite number N of modes. The number N stands for the number of active modes in the initial cw state, and not in the pulsed state. It is worth pointing out that the entropic barrier to mode locking exists only in such lasers, where the non mode-locked state is sometimes called “quasi-cw”. Section 2 introduces this model, reviews its steady-state properties, and identifies the self-start dimensionless parameter ε, which measures the relative strength of the noise and the entropic barrier. In Sec. 3 we derive the self-start rate using the overdamped Kramers escape rate formula [21, 22], and present corroborating numerical results. The last section presents our conclusions.
2. The self-starting problem
In this work we study the onset of mode locking in a model where the spectral filtering due to the finite bandwidth of the gain is implemented by assuming that the slow electric field envelope ψ is constant over a time interval of the order of the coherence length. The physical reasoning leading to this model and its properties are discussed in Refs. [15, 20]. The model offers the simplest setting where the mode locking phase transition occurs, and as such it is natural to use it to study self-starting dynamics. Although it does not account for important dispersive terms in the dynamics, like group velocity dispersion and Kerr nonlinearity, steady state studies which include these terms  indicate that they have a quantitative rather than a qualitative role in the dynamics of mode locking.
The number N of constant ψ intervals is the number of active laser modes when it is operating in cw. When N is large, cavity noise creates an entropic barrier which obstructs mode locking, and may make the cw configuration stable or metastable . Note that the cw operation we consider here, which is typical for multimode lasers, is characterized by the presence of many active modes with comparable amplitudes but inchoerent phases, and could be more precisely described as quasi-cw.
The master equation in our model system is
where ψm denotes the amplitude of the electric field envelope in the m-th interval, t is time, m stands for the time derivative of ψm , the constant γs is the nonlinear absorption coefficient, and g is the net gain coefficient. The cavity noise functions Γm are modeled by uncorrelated complex Gaussian white noise processes, with covariance
We also define the mean power
Slow saturable gain is responsible for stabilizing 𝒫 around some constant value P 0. Here we consider the simplest form of gain saturation, in which the intracavity power is kept at a strictly fixed value P 0. In this case one can obtain an explicit expression for the net gain
where the products of ψ*jΓj are in the Stratonovich interpretation .
2.1. The statistical steady state
Standard methods can be used to show that the steady-state distribution ρ of the ψm variables under the dynamics defined by (1) with (4) is 
It has been shown  that defining W=NT and taking the limit N→∞ keeping T, γs , and P 0 fixed makes ρ an equilibrium distribution of a nontrivial thermodynamic system where T is the effective temperature, and the dimensionless parameter plays the role of inverse temperature. As γ increases the system exhibits a phase transition between the cw configuration and a pulsed configuration. In the cw state the intracavity power is roughly evenly divided and |ψm |2=O(P 0) for all m, while in a pulsed there is a single site with power of O(NP 0) and the remaining power is divided between the other sites in a statistically homogeneous manner.
It has also been shown  that thermodynamic quantities are exactly given by mean field theory with a pulse-power dependent free energy (Landau function) F of the form
where ξ is the pulse power divided by p̄=W/(γsP) (the motivation for this scaling is explained below in Sec. 2.2). In particular the steady state pulse power is equal to p̄ times the abscissa of the global minimum of F.
The phase diagram is therefore determined by the qualitative properties of F and their dependence on γ. When γ<γc =4, F(ξ) has a single minimum at ξc =0, corresponding to a stable cw state. When γ>γc the minimum at zero persists, but another minimum appears at , corresponding to the pulsed state, along with a local maximum at ; when γc <γ<γe ≈4.91, F(ξp )>F(ξc ) and the pulsed state is metastable. The two states exchange stability at γe in the standard scenario of mean-field first-order phase transitions and the mode-locked state is the thermodynamically stable phase for larger values of γ. These properties of the free energy are presented graphically in Fig. (2).
It follows from these results that if a system is prepared in a pulsed state, and system parameters are changed to lower γ below γe , the pulsed state will persist unless a perturbation drives the system to the cw state, until the pulsed state becomes unstable when γ=γc ; this hysteresis scenario is born out in experiments . On the other hand, the cw state is always (meta)stable; it follows that the transition to the pulsed state can only occur as a crossing of the entropic barrier by a fluctuation. Self starting means that the barrier crossing is activated by an internal noise, which is also the origin of the entropic barrier. This fact is the root of the self starting problem: The lifetime of metastable states grows exponentially with the height of the barrier, and spontaneous activation will only be observed when γ is such that the strength of the barrier is not too large compared with the noise power.
2.2. The self-start parameter and the Arrhenius formula
The prototypical noise activated escape problem describes a particle or a system with a small number of degrees of freedom trapped in a local minimum of a confining potential and subject to thermal noise, a problem first studied quantitatively by Kramers . When the energy barrier ΔE required to escape the local minimum is much larger than the temperature, escape events are very rare and the metastable lifetime tesc is much longer than the typical dynamical timescale tdyn . In this case, the Arrhenius formula gives the lifetime as
up to an order one factor.
The self-starting problem is of a different nature: The number N of degrees of freedom is large, while the dynamics of any particular degree of freedom is devoid of a potential barrier. The entropic barrier is exposed when one writes down the equation of motion for pm =|ψm |2, the power in the m-th interval, using the usual rules of stochastic calculus
ηm (t) is a real Gaussian white noise process derived from Γm, which has pm -dependent covariance
where now the “multiplicative noise” implied by Eq. (9) is in the Ito interpretation . Note that the power of ηm vanishes at the endpoints of the allowed interval 0≤pm ≤NP 0, while the deterministic part of Eq. (8) is positive at the lower endpoint and negative at the upper endpoint; together, these two properties guarantee that pm indeed stays in this interval.
Unlike Eq. (1), the equation for pm contains a restoring term proportional to the intensity W of the noise. When pm is small, the linear restoring term in Eq. (8) dominates, and pm is trapped near zero, i. e. in the cw configuration. Assuming, as is argued below, that the term depending on pj , j≠m in Eq. (8) is subdominant, escape from cw is possible only with the help of the random term ηm , derived from the cavity noise. Hence in the self-starting problem, noise is responsible both for the existence of the barrier, and the activation which allows the system to overcome the barrier. It will be shown now, however, that the barrier height depends more strongly on the noise power than the activation, so that activation is harder and cw lifetime is longer for stronger noise intensity; this state of affairs is opposite to the standard Kramers-like scenario with a given potential, where stronger noise intensity implies a shorter lifetime.
We proceed to make an estimate of the relative strength of the noise and the barrier. In order to cast this problem in a form where the usual theory of noise activated escape can be applied we rescale the variables so as to make the confining “potential” largely independent of the noise power. For this purpose let ξ=pm /p̄, where p̄=W/(γsP) is the power needed to accumulate in a single degree of freedom for the pulse buildup process to commence. This rescaling transforms Eq. (8) into
where γ=γs /T as before, , τ=(W/P 0)t is the rescaled time variable, and the prime designates derivative with respect to τ. η is a normalized white noise, with 〈η(τ)η(τ′)〉=δ(τ-τ′).
Let us compare the deterministic drift and random diffusion terms in Eq. (10). Ignoring for a moment the terms proportional to ε, it is evident that there is an O(1) negative drift for small ξ, which becomes positive for larger ξ if γ is large enough, while the activating noise is of O(√ε). It follows that the parameter ε, rather than γ, determines the character of the self starting process. If ε=O(1) or larger, then the entropic barrier is too weak to inhibit the ordering interaction of the saturable absorber, and mode locking is achieved on a dynamical timescale 1/(γsP 0). On the other hand if ε≪1, then ξ remains trapped near zero by the drift force, and barrier crossing can occur only as a result of a rare event of a large fluctuation of η. The Arrhenius formula indicates that self-start times are then slowed down by a factor of . Since N is large in passively mode locked lasers, tends to be small even when γ is significantly larger than the threshold value needed to maintain a mode locked pulse. This observation is the basic reason that self-staring of passive mode locking difficult to achieve in many practical systems, and external perturbations are needed instead to drive these systems to mode locked operation. The case ε≪1 is the one considered in this paper. In this case the drift terms proportional to ε in Eq. (10) are small, and do not change the physical picture just described; nevertheless, because of the strong dependence of the self-start time on the drift, these small terms have an appreciable effect on the self-start process, and cannot be neglected. In the next section we carefully analyze the term proportional to ∑j≠m , show that it is indeed small, and find its contribution to the self-start time.
3. The self-start time distribution
3.1. The mean field approximation
In the rest of this paper it is assumed that ε≪1. In this case self-starting occurs via noise activated barrier crossing, and the cw lifetime is exponentially larger than the other timescales in the system. Since the pulse buildup process occurs on a dynamical timescale, this fact implies that the probability of self starting occuring simultaneously in two different sites is exponentially small, and it can be safely assumed that the power at no more than a single site may simultaneously reach values which are much larger than the mean power P 0. Therefore the self-starting time is determined by the fastest of a system of N-independent escape processes, each described by Eq. (10).
The calculation of the mean self-start time is therefore reduced to that of a single degree of freedom ξ, governed by Eq. (10). However, this is still not a single-variable escape problem, becuase it contains the term
which couples the variable ξ to the N-1 other degrees of freedom in the system. Fortunately, in the case of interest where ε≪1, we claim that the dynamics are such that the distribution of Q is determined to leading order by ξ, enabling us to write a statistically equivalent equation, which involves ξ alone.
For this purpose, consider the dynamics of ψj with j≠m. Since self-starting occurs predominantly through the growth of a single degree of freedom, it can be assumed that ψj remains of O(√P 0) throughout the process, and in the leading order (in ε) the nonlinear terms in its dynamics, except those involving ψm can be neglected. It follows that the probability distribution of ψj , for a given value of ψm , is gaussian. Q is therefore the sum of the fourth powers of independent centered gaussian variables, the sum of whose variances is NP 0-p̄ ξ; from this it follows that the expectation of Q is  . As a sum of many independent variables, the fluctuations in Q are smaller by a factor of O(√N) than its expectation, and will be neglected.
We now make the mean-field approximation by replacing Q in (10) with its expectation value conditional on ξ. The resulting equation in rescaled coordinates for the self-starting process at site m is then
3.2. The mean lifetime of cw
The self-starting problem has been reduced in Eq. (12) to a standard problem of noise-activated escape in one dimension, whose solution is well-known [25, 28]. When γ>4, the O(1) part of the drift force has three fixed points ξc , ξb and ξp , whose values are precisely those of the three extrema of F of Sec. 2.1. Since ε≪1, there are exact fixed points of the drift near ξc , ξb , and ξp , the first and last of which are (dynamically) stable, while the one near ξb is unstable. It follows that the drift force is derivable from a potential whose shape is well-represented by the steady-state free energy function F shown in Fig. 2, with a potential barrier near ξb that inhibits self-starting. The escape time may be defined as the first time the variable ξ with zero initial value reaches an arbitrary value strictly between ξb and ξp . For typical escape times, which are much longer than the dynamical time scale, the distribtution of escape times is Poissonic . This fundamental property follows from the fact that the dynamics are Markovian, and from the principle of separation of scales : the system reaches a quasi-steady state very fast compared to τesc , and therefore the probability that it escapes is independent of history.
The escape time τesc depends very little on the precise choice of the target pulse size , and is well-approximated by the mean first passage time to reach . The self-start time tcw is obtained from τesc by converting back to the physical variable t and dividing the result by N to take into account the fact that self start can occur via N independent escape processes. As the calculation of the mean first passage time for stochastic equations of the type of Eq. (12) is standard [25, 28], we merely cite the result, leaving details to the appendix
where F(ξ) is again the free energy function. Eq. (13) is the main quantitative result of this work. It contains an exponential dependence of the self-start time on ε, which is equal to the one which can be deduced heuristically by applying the Arrhenius formula to the steady state free energy, and an algebraic prefactor, which has a ε 5/2 dependence; it is verified by numerical simulationz in Sec. 3.4.
The present analysis and its result Eq. (13) are valid when ε≪1 for any value of γ>4; it is accompanied by a reverse process of activation from the pulsed state back to cw, which is not studied here. The steady state analysis  tells us that when γ<γe the rate of activation of mode locking is slower than the rate of cw activation from the metastable mode-locked state, while the converse is true when γ>γe , in which case the process indeed describes self-starting of mode locking.
Because tcw grows very fast as ε decreases, systems where self-starting may be practically observed involve ε values which are not too small; for example if F(ξb )/ε is larger than 100, say, the probability of observing a self-starting event is less than 10-35 per second. Since , and N is large, this means that γ may be assumed large in many cases. If γ≫ε -1, or equivalently, γ≫√N then Eq. (13) simplifies to
Rather large values of N are needed for this asymptotic expression for the mean lifetime to become precise, as can be seen in Fig. 3, where the exact and asymptotic expressions for the mean lifetime Eqs. (13) and (14) are displayed as a function of system parameters. For reference purposes we also include a version of Eq. (14) in the original unscaled variables,
3.3. Discussion and comparison with experiments
The formulation of the self-starting problem as a noise activated escape process implies that the mean self start time is an exponentially increasing function of the barrier height. Since the latter is a decreasing function of the intracavity power, a threshold-like behavior of the self-starting process has been often observed in experimental tests of self-starting with varying power.
Our theory which explains the self-starting process as an entropic barrier crossing is compatible with the idea put forward by [4, 8, 7] that the action of the saturable absorber to align the phases of the laser modes is counterbalanced by a noise-induced decoherence process, that tends to randomize them. The decoherence process and the pulse-buildup process are each characterized by an associated time scale, and self-starting occurs, according to the decoherence time theory, when the pulse-buildup time is of the order of or shorter than the decoherence time scale. The dimensionless ratio of the two time scales in our model is equal to γs/2W=1/2ε, and the condition for self starting based on the decoherence argument is therefore ε≳1. As shown above, this domain is where the noise is too weak for an entropic barrier to form, and self-starting occurs on a dynamical time scale. When ε≪1 a significant entropic barrier forms, and the simple picture of decoherence versus saturable absorber-induced ordering is insufficient; rather, a full stochastic analysis is necessary leading ultimately to Eqs. (13) and (14).
The theory and results obtained here motivated a new experimental study of self-starting behavior in an additive-pulse mode locked fiber laser  focusing on the activation regime ε≪1, and the experimental results agree well with our theoretical predictions. Firstly, as shown in the right panel of Fig. 4, the experiments have confirmed that the distribution of the self-start times has an exponential tail, as expected in this regime. Furthermore the approximately linear relation between logtcw is 1/ shown in the right panel of Fig. 4 agrees with our prediction in Eq. (14), which is remarkable, since the latter result has been obtained in the framework of a simplified model, and leads us to believe that this property of the self start time has a wide range of validity.
3.4. Numerical analysis
The theoretical analysis leading to the asymptotic expression for the mean lifetime Eq. (13) includes some reasoning which, while strongly supported by physical arguments, is hard to justify rigorously.We therefore performed full numerical simulations of Eq. (1) as an independent test of the theory. The numerical results, to be discussed presently below, confirm in a satisfactory manner the theoretical predictions, including the prefactor in Eq. (13).
We used the stochastic Euler method  and imposed the constraint (4) by normalizing 𝒫 to P 0 in every time step. For a given value of N the problem depends on γ as a single dimensionless parameter. The initial values of the ψj ’s were chosen as independent samples a random variable with a complex Gaussian distribution and normalized to the appropriate value of P 0. The time when the order parameter ∑j|ψj |4 has reached half its maximal possible value was recorded as the self-starting time.
As explained above, the theory developed in this paper has two main predictions, the first of which is that the distribution of self-start times is Poissonic beyond dynamical times, or in other words, that the probability that the system remains in the cw state decays exponentially. This prediction is confirmed by the numerical simulations, as shown in Fig. 5 where the measured probability distribution of self-start times is displayed. It is interesting to note that since this property is a consequence of the statistics of rare events, it is necessary to simulate the noise by a high-quality random number generator to obtain it correctly. In the simulations reported here we used the generator ran2 of .
The second theoretical prediction concerns the mean self-start time which is given by the asymptotic approximations (13). The theoretical prediction is shown as a solid line in Fig. 6, and compared with the numerical measurements shown as red crosses. As can be seen in the figure, there is an excellent agreement between Eq. (13) and simulation, indicating that the exponential dependence on system paramaters as well as the preexponetial factor are correctly reproduced. The simulations were carried out in a parameter regime where γ~ε -1 and Eq. (14) is not valid. Because of the ε 5/2 preexponential factor, a considerably larger number of active modes N is needed to reach the region of validity of the latter expression; on the other hand, such values of N and larger ones are quite common in actual laser systems.
The main result of this work is the theoretical demonstration that the process of self starting of passive mode locking has the dynamical significance of a noise-activated crossing of an entropic barrier, and the quantitative expression for the mean self start time.We found the dimensionless parameter that governs the process, which is equal to the the ratio of two time scales, dynamical and noise-related, in accordance with some earlier studies of self starting. Furthermore, numerical simulations corroborate the theoretical predictions with excellent accuracy, and the theory is also in good agreement with recent experiments in which the self-start times distribution was measured.
The self-starting parameter, which dictates the nature of the self-starting dynamics, is inversely proportional to the large number of active laser modes (in its non-pulse state). This fact is the origin of the self-starting problem, and implies that self starting is harder to achieve when the active bandwidth of the laser becomes larger.
The results presented in this paper have been obtained in the framework of a simplified model of a passively mode locked laser, where fine details of the nonlinear pulse buildup process and of the gain dyamics and its spectral filtering are not taken into account. Nevertheless, the fact that the effective activation barrier is determined by the maximal value of the static free energy divided by the self-starting parameter can serve as a good guiding principle for making an estimate of the self-starting behavior of more realistic laser systems.
We acknowledge support by the Israel Science Foundation and by AFOSR under contract FA9550-04-01-0011.
A. Appendix: Calculation of the mean self-start time
The calculation of the mean escape time for one-dimensinal Markovian activation problems such as Eq. (12) is standard [25, 28], and is outlined in this appendix for the sake of completeness. The mean first passage time τ (ξ) to reach the target point starting from the point ξ obeys the equation
with the boundary conditions τ()=0, and τ′(0)=0 (reflecting boundary).
It is straightforward to solve this equation in quadratures for τesc =τ(0),
where the function F is the steady state free energy defined in Sec. 2.1.
In the case of interest, ε≪1, the main contribution to the integral arises from the vicinity of the maximum of F(ξ)-F(ξ′) in the region 0≤ξ′≤ξ≤, which is reached when ξ=ξb , the abcissa of the potential barrier of F, and ξ′=0. The leading term in the asymptotic expansion of τesc in powers of ε is obtained by approximating the integrand as a function quadratic in ξ and linear in ξ′ near the maximum (Laplace’s method), giving
Eq. (13) is obtained from this result upon scaling back to the physical time variable t and division by N to take into account the N possible self-starting paths.
References and links
7. Ch. Spielman, F. Krausz, T. Brabec, E. Wintner, and A. J. Schmidt, “Experimental study of additive-pulse mode locking in an Nd:Glass laser”, IEEE J. Quantum Electron. 27, 1207 (1991) [CrossRef]
10. J. Hermann, “Starting dynamic, self-starting condition and mode-locking threshold in passive, coupled-cavity or Kerr-lens mode locked solid-state lasers”, Opt. Comm. 98, 111 (1993). [CrossRef]
11. A. K. Komarov, K. P. Komarov, and F. M. Mitschke, “Phase-modulation bistability and threshold self-start of laser passive mode locking”, Phys. Rev. A. 65, 053803
12. J. M. Soto-Crespo, N. Akhmediev, and G. Town, “Continuous-wave versus pulse regime in a passively mode-locked laser with a fast saturable absorber”, J. Opt. Soc. Am. B 19, 234 (2002) [CrossRef]
13. H. A. Haus, “Theory of mode locking with a fast saturable absorber”, J. Appl. Phys. , 46, 3049 (1975). [CrossRef]
14. H. A. Haus, “Mode-Locking of Lasers”, IEEE J. Sel. Top. Quantum Electron. 6, 1173 (2000) [CrossRef]
16. A. Gordon and B. Fischer, Phase transition theory of pulse formation in passively mode-locked lasers with dispersion and Kerr nonlinearity” Opt. Commun. 223, 151 (2003). [CrossRef]
17. A. Gordon and B. Fischer, “Inhibition of modulation instability in lasers by noise”, Opt. Lett 18, 1326 (2003). [CrossRef]
18. O. Gat, A. Gordon, and B. Fischer, “Light-mode locking - A new class of solvable statistical physics systems”, New J. Phys. 7, 151 (2005) [CrossRef]
20. O. Gat, A. Gordon, and B. Fischer, “Solution of a statistical mechanics model for pulse formation in lasers”, Phys. Rev. E. 70, 046108 (2004) [CrossRef]
21. H. A. Kramers, “Brownian motion in a field of fource and the diffusion model of chemical reactions”, Physica (Utrecht) 7, 284 (1940). [CrossRef]
22. P. Hänggi, P. Talkner, and M. Borkovec, “Reaction-rate theory: fifty years after Kramers,” Rev. Mod. Phys. 62, 251 (1990). [CrossRef]
23. B. Vodonos, A. Bekker, V. Smulakovsky, A. Gordon, O. Gat, N. K. Berger, and B. Fischer “Experimental study of the stochastic nature of the pulsation self-starting process in passive mode-locking,” Opt. Lett. 30, 2787 (2005). [CrossRef] [PubMed]
25. R. L. Stratonovich, “Some Markov methods in the theory of stochastic processes in nonlinear dynamical systems”, in “Noise in nonlinear dynamical systems, Vol. 1, edited by F. Moss and P. V. E. McClintock, Cambridge University Press (1989)
26. B. Vodonos, R. Weill, A. Gordon, A. Bekker, V. Smulakovsky, O. Gat, and Baruch Fischer “Formation and Annihilation of Laser Pulse Quanta in a Thermodynamic-like Pathway”, Phys. Rev. Lett. 93,153901 (2004). [CrossRef] [PubMed]
27. H. H. Risken, “The Fokker-Planck Equation”, Second edition, Springler-Verlag (1989, 1996).
28. C. W. GardinerHandbook of Stochastic Methods, 3rd ed., Springer, New York (2004).
29. P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, Springer New York (2000).
30. W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical recipes in C: The Art of Scientific Computing, p. 281., 2nd edition, Cambridge University Press, New York (1992).