# Numerical Simulations on Nonlinear Dynamics in Lasers as Related High Energy Physics Phenomena

Authors: Andreea Rodica Sterian
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

This paper aims to present some results on nonlinear dynamics in active  nanostructures as lasers with quantum wells and erbium doped laser systems using mathematical models, methods, and numerical simulations for some related high energy physics phenomena. We discuss nonlinear dynamics of laser with quantum wells and of fiber optics  laser and soliton interactions. The results presented have important implications in particle detection and postdetection processing of information as well as in soliton generation and amplification or in the case that these simulations are thought to be useful in the experiments concerning the high energy particles. The soliton behaviour as particle offers the possibility to use solitons for better understanding of real particles in this field. The developed numerical models concerning nonlinear dynamics in nanostructured lasers, erbium doped laser systems, the soliton interactions, and the obtained results are consistent with the existing data in the literature.


---

## Body

## 1. Introduction

Numerical simulations on the nonlinear dynamics in lasers are largely presented in the literature, but the research work for their development has determined unprecedented progress in particle physics [1–8].The research put into evidence the existence of relevant conditions and limitations which are important in the particle detectors and data processing after detection and in the experimental studies in high energy physics also. These unprecedented advances offer opportunities in theoretical and numerical predictions in high energy physics and also in the related fields. The particles having soliton behavior, as being localized and preserved under collisions, lead to a large number of applications but offer the possibility to use solitons for better understanding of the real particles, in high energy physics. On the other hand, the research developing in Romania concerning the Project “Extreme Light Infrastructure-Nuclear Physics” (ELI-NP) demonstrate the complexity of problems in the physics and engineering of high energy particles [9]. The ELI-NP facilities, imply as main tools for the research directions, and the realization of the laser systems of high power of petawatt and ultrashort pulses together with a gamma beam system, highly collimated and very brilliant and tunable [10, 11], so that numerical simulations on nonlinear dynamics are essential and inevitable.The mathematical models, the methods, and the numerical simulations from nonlinear dynamics are important in processing of data concerning particle physics. Some numerical results concerning the dynamics of a two-level medium under the action of short optical pulses are presented in the papers [12–18]. We analyse in the paper the nonlinear dynamics of the lasers with quantum wells and optic fibers as well as the soliton collisions and the behavior of particles in high energy physics as solitons. The used methods are typical for the study of quantum fields and particles, particularly cosmic rays [19].


## 2. Modulation Capabilities and Bifurcation Diagrams

The numerical simulations of the lasers with quantum wells (QWL) are based on solving the nonlinear differential equations which describe the behaviour of this kind of devices [20–31]. Study of functional characteristics of QWL is possibly based on the models present in [20], by implementing programs to integrate these equations in MathCad using appropriate algorithms. The MathCad software has the advantage of solving complicated systems starting from the usual mathematical form by comparison with other methods of integration of the nonlinear equations. The modulation studies are important, and for solitons generation in optical fibers also, the phase self-modulation phenomena are the essential mechanism of this process [26]. In general, all the analysed models had given satisfactory results for the study of transient and dynamics operational regimes at the small signal of injection. At high level of injection, the obtained numerical results had showed the limits of specific model used but are consistent with the theoretical analysis.We have considered the amplitude modulation of the current of injection [27]. The output power of the laser versus the modulation index was investigated (Figure 1).Figure 1
Bifurcation diagram: photon density dependence of modulation index.Graphical representation of the modulated signal simultaneously with the evolution of the system in the phase space highlights also the route to chaos with the increase ofm by doubling the period.We have examined the diagram of bifurcation of the photon density versus the modulation index directly (with the increase ofm) and vice versa (with decreasing of m) for a cycle of modulation as an indicator of the system progress. The chaotic region has periodic windows and a dynamic type of hysteresis.We have considered also the additional noise sources in the rate equations [28]. In the Markov approximation, the noise sources are random Gaussian variables with zero mean and delta-correlated quantities [24]. For numerical representation of the dynamics including the noise sources, we had used the random Gaussian numbers applied by MATLAB functions. Numerical calculations showed that a small noise disturbs and can also occur transitions between coexisting attractors (Figure 2).Gaussian noise influence on the system evolution in the phase plane.
(a)
(b)
(c)For a suitable modulation index, there are two coexisting stable attractors as it is showed in Figures2(a) and 2(b), but in the presence of the noise, the induced transitions from one attractor to another can occur (Figure 2(c)). The contribution of the included noise had proved to be an important factor that limits the design of the format of modulation.


## 3. Case Study on Nonlinear Dynamics

A case study in nonlinear dynamics of lasers with quantum wells is analysed systematically based on the model with two equations proposed by Anghel et al. [29–31]. These equations were solved numerically using typical parameters of a laser with multiple quantum wells. The parameter values used in numerical simulations are presented in the papers [20, 21]. Simulations were performed in MATLAB programming language using a scheme Runge-Kutta of order four. Detailed numerical investigations showed that for the lower frequency compared to the frequency of relaxation, dynamics of the system is periodic with the period of the radio-frequency modulating signal.For a sufficiently low modulation index is obtained a sinusoidal periodic regime, so that laser output power is sinusoidal. For a growth of the modulation index and the frequency of modulation near the frequency of the relaxation oscillations, a pulsing mode of operation is obtained with a duration of tens of picoseconds (ultrashort pulses), so that the laser system can be used as a source of pulses for optical communications. Higher frequency modulation causes multiple bifurcation points including a chaotic dynamics for certain parameter values (Figure3). Also, in Figure 3 critical dependence of the modulation index (which leads to the bifurcation) of frequency is observed.The laser bifurcation diagrams for different modulation frequencies.
(a)
(b)The behavior is similar for other values of the injection current above the threshold. However, an increasing dependence of the modulation index for the doubling of period with the current by the polarization can be observed.


## 4. Nonlinear Dynamics in Optical Fibers Systems

In recent years much attention has been paid to the study of nonlinear dynamics in optical fiber lasers [32–39].Numerical simulations discussed here refer to nonlinear effects in such kind of systems. The chaotic dynamics was studied numerically using the rate equations, based on phenomenon of formation of the ion pairs [40–42]. It was demonstrated the stable, nonchaotic operation of the analyzed laser system and discussed the performance of their modulation using the methods of communications theory.Numerical modeling was realized using the MATLAB programming medium. The base element of the program was the function ode 45, which realize the integration of the right side expressions of the nonlinear coupled equation using the method Runge-Kutta. The self-pulsing and chaotic functioning of the EDFLs have been reported in various experimental conditions including in the case the pumping is near the laser threshold [32–35]. A model for the single-mode laser was analysed firstly taking into account the presence of the erbium ion pairs that act as a saturable absorber. For this type of laser system don’t yet completely are known the interaction mechanisms, in spite of many published works. The nonlinear dynamics of an erbium doped fiber laser is explained based on a simple model of the ion pairs present in heavily doped fibers. The dynamics is reducible to four nonlinear differential equations. Next, the nonlinear dynamics of a two-mode laser with optical fiber is explained in terms of a classical two-modes laser model and two additional equations of the quantum states of ion pairs [40], so that the laser dynamics is described through six coupled differential equations. At a suitable level of the concentration of the ion pairs, the laser system is in a steady state or self-pulsing operation regime, depending on the value of the pumping parameter. For larger values of the ion pairs, there is a range of self-pulsation dynamics. This is showed by means of a bifurcation diagram in which the mode intensity of the laser is plotted against the pumping parameter. This dynamics has also been encountered in the case of the one-mode model, the change from a stationary stable state to a self-pulsing one being a Hopf bifurcation.Figure4(a) shows the maxima of intensity at a larger concentration of the ion pairs. There exists a region of two maximum intensities, that is, a period doubling. A small increase of the ion pairs concentration (Figure 4(b)) leads to the appearance of a window of a tripled period. Further increase of the ion-pairs concentrations leads to a more complicated bifurcation diagram, including of the regions with multiple stable states (generalized bistability) and routes to chaotic dynamics such as the period-doubling route or a quasiperiodic one.Bifurcation diagram of the mode intensity for two values of ion pairs versus the pumping strength.
(a)
(b)The used method of the fourth-order Runge-Kutta for the numerical simulation demonstrate the importance of “computer experiments” in the designing, improving, and optimization of these coherent optical systems.Some new features of the computer analysed systems and the existence of some new situations have been put into evidence for their use by the designers in different applications [42–48]. The obtained results are important not only for the optimization of the functioning conditions of this kind of devices but also for designing of the systems for the processing of data in the particle physics and in the systems engineering of high energy.


## 5. Numerical Simulation of Soliton Propagation

In the paper [49] the specific features of the soliton pulses were studied, so: the type of the nonlinearity, dispersive properties of the propagation medium and the propagation equations, for the optical and acoustic solitons. Were analysed the specific equations as quadratic nonlinear Schrodinger (QNLS), Kerr cubic (CNLS), Saturable (SNLS), discretized (DNLS), Korteveg de Vries (KdV), Boussinesq, Burgers, and sine-Gordon. Also a study of the methods of discretization used by the different numerical simulations was realised.A numerical simulation of soliton propagation, based on the Korteweg-de Vries equation, using a powerful PC program (Maple10) which permits to perform numerical calculations, plot or animate functions, and manage analytical expressions is presented in paper [50]. We discuss a model of the interaction between two solitons, taking into account both the fiber dispersion and nonlinearity. The numerical simulations show how the solitons are propagating in optical fibers and how two solitons interact, passing one through another, with only a phase change. These simulations are thought to be useful for the designers working in the digital data transmission and in particles physics as well. The program Maple10 helps us plot solutions (we can use animation or individual frames) for different values of the wave number k, for different moments t, and for an arbitrary interval for x. For the study of the interaction between two solitons, we need a solution of the KdV equation present in [3, 50]:
(1)u2∶=2(-k22k12cosh(k1(x-4k12t))2+k22k12cosh(k2(x-4k22t))2-k22k12-k14cosh(k2(x-4k22t))2+k14+cosh(k1(x-4k12t))2k24)×(∫(cosh(k1(x-4k12t))cosh(k2(x-4k22t))k2-sinh(k2(x-4k22t))×sinh(k1(x-4k12t))k1)2∫)-1.
With that solution, Maple10 helps us visualize the interaction between the two solitons as in Figure 5.Figure 5
The interaction of two solitons fork1=0.4;  k2=0.9;  x=-50⋯50;  t = −10;  −2; 0; 2; 10.We can see at both solitons the regaining of their shapes after collision (after they pass one through another). Obviously, it is better to use the animation feature of the Maple10 program by modification of the parameters and observing the changes.For designers working in digital data transmission, the simulation of the solitons propagation along optical fibers is useful and may be performed successfully also with Maple10. The animated images are some useful features of this computer program. Especially, the collision of two solitons and ofk dependence of the soliton propagation can be studied by this means, respectively.As it was stressed above, after the collision, the solitons regain their original shape and velocity [8, 50, 51]. The only remaining effect of the scattering is a phase shift (i.e., a change in the position they would have reached without interaction). For the trisolitonic model, using Maple11 similarly, the solution is plotted in Figure 6, showing the collision of three solitons [49, 51], by choosing various values for the parameters.Figure 6
The interaction of three solitons.


## 6. Solitons Like Behavior of Particles in High Energy Physics

The behavior as particle of the solitons (they are localised and preserved in the collision) leads to a large number of applications and is important in the particle physics also [50]. It is well known that we hope the solitons to be used for a better understanding of real particles. There are already soliton models of nuclei. On the other hand, there are theories which permit us to see particles such as electrons and positrons as being solitons in appropriate cases. In the paper [52] they had proposed a general semiclassical method for computing of the probability of the production of soliton-antisoliton pairs in the particle collisions. They calculated also the probability of this process at high energies. In the paper [53] they had studied the creation of the solitons from the particles, using an suitable model used as a prototype. It was considered the scattering of the small identical wave pulses, which are equivalent to a sequence of particles, and they found that kink-antikink pairs are created in a large region of the parameters space. The dissertation [54] is devoted to the study of monopoles and solitons in particle physics. The monopoles and solitons are of relevance as intermediate and final states of the particles in the collision experiments. The study of monopoles and solitons in the theory of strong interactions had shown the possible interpretation of the solitons as baryons. However, so far there is no general theory in which the particles are described as solitons. There are macroscopic phenomena such as the spontaneous transparency and the behavior of light in fiber optics that are understood only in terms of solitons. The methods specified above are typical for the quantum fields and the particles study, particularly of the cosmic rays [55–61].


## 7. Conclusions

(i)
The nonlinear dynamics in lasers was treated taking into account the complexity of the phenomena which govern the field-substance interaction with applications in particle detectors and data processing after detection, but these methods are typical for the study and for quantum fields and particles also.
(ii)
A computational software was elaborated to study the nonlinear phenomena in quantum well lasers (QWL). Mathcad proved to be a powerful tool in analyzing and designing of these models. For the QWL, we put into evidence their complex nonlinear dynamics, characterized by bifurcation points and chaos, the critical values of the parameters being determined.
(iii)
The classical Runge-Kutta method was used successfully for numerical simulation of the complex nonlinear dynamics in lasers with optical fibers. The chaotic dynamics was studied numerically based on the rate equations and on the phenomenon of formation of the ion pairs.
(iv)
Our results put into evidence the existence of the new situations which are important for the optimization of the functioning conditions for this kind of devices. The nonlinear dynamics can be controlled, and undesired situations can be avoided if it is necessary, by suitable parameters choosing.
(v)
The studies on laser modulation are important for the generation of solitons in optical fibers, the phase self-modulation phenomenon being the essential mechanism of this process. The study of the interaction between two solitons using a powerful PC program (Maple10) is of great interest for a better understanding of real particles in spite of the fact that now there is no general theory in which the particles are described as solitons.
(vi)
The developed numerical methods and the obtained results are consistent with the existing data in the literature. The methods of the computer experiments make a complex study possible taking into account the intercorrelations of the parameters by simulating experimental conditions, as it has been shown, so that the obtained results can be successfully applied for processing of data in the particle physics.



---

*Source: 516396-2013-10-27.xml*