Impulse Spatial-Temporal Domains in Semiconductor Laser with Feedback

An initial value boundary problem for system of diffusion equations with delay arguments and dynamic nonlinear boundary conditions is considered. The problem describes evolution of the carrier density and the radiation density in the semiconductor laser or laser diodes with “memory” and with feedback. It is shown that the boundary problem can be reduced to a system of difference equations with continuous time. For large times, solutions of these equations tend to piecewise constant asymptotic periodic wave functions which represent chain of shock waves with finite or infinite points of discontinuities on a period. Applications to the optical systems with linear media and nonlinear surface optical properties with feedback have been done. The results are compared with the experiment.


Introduction
We consider semiconductor lasers or laser diodes [1]- [11].The laser is inverted carrier density system.There exist carrier generation and recombination that is when electrons interact with hole, they recombine.The energy released can be produced by thermal recombination or optical photon recombination which is used in semiconductor lasers.An electronic oscillator is an electronic circuit that produces a periodic signal.Oscillators convert direct current to an alternative current signal.We use the feedback oscillator which can increase amplitudes of signal.
We consider "surface" oscillator which can be described mathematically by the functional or dynamic boundary conditions with feedback of the following form: or , where input u is the corresponding amplitude of input signal, and output u is the amplitude of output signal.Here, , : I I Φ Ψ → are given nonlinear functions which model the transform of signal with help of laser diodes or bipolar junction transistor.
Next, we can use semiconductor laser with nonlinear layer which produces the process of recombination of electrons and holes, and the density of radiation.For surface inverse system, this process of recombination can be given by nonlinear function with feedback.Such nonlinearity may be produced by heating.
There are different physical mechanisms which can convert an optical signal.We consider an "deal" resonator when conversion of signal with feedback takes place at walls which confines the resonator.It will be shown that nonlinear surface optical properties of material produce in the volume of ideal resonator asymptotic periodic piecewise constant wave structures with finite or infinite number of "jumps" of amplitude.Such distributions of electron, holes and radiation together represent shock waves.If a number of "jumps" are finite on a period, we have dealt asymptotic distributions of relaxation type as time t → ∞ .If this number is infinite countable or infinite uncountable on a period, we get limit shock waves of pre-turbulent and turbulent type, correspondingly.Such periodic shock waves take place in n p − -type semiconductor lasers.
A corresponding mathematical model can be described for semiconductors with "memory" by a system of diffusion equations with delay arguments with dynamic nonlinear boundary conditions.To be more precise, the structure of shock waves depends not only on surface structure of material, which is described mathematically by the boundary conditions, but also from the initial data of the boundary problem that is from initial distribution of electrons, holes and density of radiation in the semiconductor laser.
In last years, in physics studied the nonlinear interaction of light which can mimic the physics at so called an event horizon.As shown in [6], this analogue arises when a weak probe wave is unable to pass through an intense soliton, despite propagating at a different velocity.These dynamics arise as a soliton-induced refractive index barrier.In all papers this barrier characterizes the volume optic properties of linear boundary conditions.In this paper, we consider the opposite problem when the optical medium is ideal or linear, but boundaries of the medium have the nonlinear optic properties, and describe, for example, the all-optical transistor [5].It may be also a case when a bright soliton is passing through the soliton.In this case, the intensity of light depends on the refractive index.Thus, the soliton creates a moving refractive index's perturbations which passage through the another soliton [1] [8].This interaction between such surface solitons plays the main role of distributions of the light in the linear medium with nonlinear surface interaction.Thus for the ideal medium the main role plays the surface nonlinear refractive index.In ( [5], Figure 1) has been mimicked two spectral modes of solitons when the mode-locked laser diode generate picosecond solitons.This generation will be described as a functional boundary conditions with feadback generated these solitons.As noted in [1], "intense pulse of light may be used to create an effective flowing medium which mimics certain properties of black hole physics".Of course, analogues models can be realized in very different physical systems.
In this paper, an initial value boundary problem (IVBP) of system of two linear wave equations with nonlinear boundary conditions has been considered.Solutions of this problem describes the propagation of density of optical radiation (on a given frequency) and electron carriers in one-dimensional in the ideal semiconductor rod with optical defects at the ends.For example, the corresponding mathematical model describes the wave distributions of the density of radiation of photons and electron density in an ideal semiconductor.The semiconductor is confined by two flat walls which emit or absorb light.The probability of absorption or emission of photons depend on the surface density of the radiation and the surface density of electrons in a nonlinear manner.Diffusion in semiconductors is one-dimensional and one is directed orthogonal to the flat walls.Thus, the initial boundary value problem for the two linear wave type equations with nonlinear boundary conditions will be considered.This problem models different optical phenomena as white and black solitons, propagation of light in resonator with feedback connection between beams at walls and so on.
In this paper, we study the structure of attractor of the IVBP.The IVBP admits a Figure 1.Solutions of turbulent type with the limit cantor set of points of discontinuities on a period.Phase portrait for simplest solutions of relaxation type (white and black solitons).
reduction to a differential-difference equations (DDE).We restrict ourselves to the case when the corresponding DDE are completely integrable.Indeed, among the differentialfunctional equations [12] ( ) ( where A is a differential operator, ( ), 1,   n k k n τ = are delay arguments, there are a class of equations which admit decomposition on the finite product of differential and functional operators.Such equations are called equations with splitting operators [12].
The study of such equations reduces to a serial study of differential and functional equations.To this class belong to the so-called completely integrable equations [12] such that A DJ = , where D is a differential operator and J is a functional operator.For example, the equation ( ) ( ) ( ) is a completely integrable.Here, , where F is some primitive of the function f.By integrating this equation can be reduced to one-parameter family of difference equations (DE) with continuous time: ( ) ( ) where c is some integrating constant.In general, the completely integrable equation can be reduced to the family of functional non-autonomous equations ( ) ( ) where ( ) We consider distributions in a semiconductor in the region 0 , 0 x l t < < > , where : is the density of radiation on a given frequency ω .The boundary conditions at have the form: , , 0 at 0 where , , P Q R are continuously differentiable functions.Define ( ) , n I µ one of the integrating factor of Equation (8).Then, after multiplying on µ , Equation ( 8) can be written as: Since all the integrating factor described by the formula ( ) and the initial conditions Here, n τ and I τ are response times of the carrier density radiation density corre- sponding to external perturbations in the semiconductor.α is an absorption coefficient, β is an emission coefficient; n D and I D are diffusion coefficient.Such approach previously applied to the study of binary alloys on the example of Kahn-Hilliard equation with delay argument to applications to binary mixtures [13]- [16] and to binary polymer blends [17]- [19].Such problem describes by hyperbolic equations.Evolution of distributions satisfies to the so-called Non-Fickian diffusion and and this distributions represent "tau approximation" for numerical turbulence [20].Below these idea will be applied to the modelling of carriers distributions in semiconductor lasers.At first, such approach has been considered by Maxwell (see, [21]).
Thus, we consider asymptotic behavior of solutions for this of linear uncoupled equations with nonlinear dynamic boundary conditions.We define conditions on parameters, boundary functions and the initial conditions where there are asymptotic periodic piecewise constant solutions with finite or infinite points of discontinuities on a period.The paper is organized as follows.

Formulation of Problem
Consider a uniform rod such that the axis (Oz-axis) extends monochromatic radiation, characterized by the frequency ω so that so that the density of the radiation is S x be a density of the light flow at a point x and let the rod has an active and nonactive centers that are capable of absorbing radiation at the frequency ω .Then the active centers emission process is accompanied by opposite process so that the resultant absorption is determined by the competition of these processes.
Nonactive absorption centers (on frequency ω ) are not accompanied by the emission: for such absorbtion we can formally include the scattering of radiation in the rod, which leads to energy loss through the side surface (see, [22], p. 11).
So we write the Burger law ( ) where l α is the linear absorption coefficient on active centers, n α is the corre- sponding coefficient in the inactive centers.Note that the competition transitions (absorption and stimulated emission) can lead to a negative value l α .In this case, the optical parameter of the medium describe the negative absorption, i.e., the increase in radiation.Due to the nonlinear optical effect (in the centers),  , where G is a given nonlinear function (see, [22], p. 11).
The corresponding model equation for densities of electrons and radiation is: , where e n and p n are the diffusion coefficients for electrons and holes.We assume , where l is the size of system and , n I τ τ are times of relaxations, correspondingly.Then these equations can be written as , , . Now, instead of diffusion Equation (19) and Equation (20), we consider the diffusion equations with delay arguments: , , , , where 1 r n t Here, r t and s t are times of response of the system on outer perturbations, correspondingly.These equations describes distributions of the number of electrons and radiation in semiconductor.We assume that the number of electrons and holes is constant.The rise times r t and s t corresponds to impulse response that is the reaction of dynamic system on some external change.We assume that 1 1 τ  and 2 1 τ  .Then these equations, with accuracy ( ) O τ and ( ) O τ can be written as: , , , , .
= , these equations can be written in the unification form: , , , , .
Bellow index in 1 t will be omitted.
Note that the rise time can be determined by the formula in the quadratic approximation, and by the formula in the common case.Here, ς is the damping ratio and r ω is the natural frequency of the network.The value ς is a dimensionless measure which describes how oscillations decay after perturbations.
Note that the rise time s t is derived from the assumption that the response of the medium at a time t determined by the field of the light wave at the same time t so that ( ) ( ) ( ) , where ( ) E t is an electric field.However, in fact it should be taken into account (inevitable) "inertness" of the medium.This means that the response of the polarization on an action of the field take place with delay time.Generally speaking, the time is not precisely defined Libenson (see, [22], p. 28).So the polarization of the medium at a time t to be determined by the wave field in all previous times t t t ′ = − ∆ , where the delay is positive.
Let us consider for Equation (25) the dynamic boundary conditions and the Neumann boundary conditions and the initial conditions Let the solutions are , , where . Substituting (33), (33) into (25), (26), we obtain the equations From ( 33), (33) it follows that Equations (35) can be written as Without loss of generality, we assume that From the Neumann boundary conditions it follows that Integrating this equation from 0 t , we obtain that Next, integrating the equation , we obtain that From ( 44), (46) it follows that Then from the dynamic boundary conditions it follows that .
From (42) it follows that Equation (48) can be written as where . Without loss of generality, we assume that 2 : F f f → .Then Equ- ation (49) can be written as : where I is an interval.In 2 C structural stable maps form an open dense subset.
is an initial function for Equation (50).Then from (50) we can find a solution of this equation on the interval [ ) 0, l V and so on.This solution is bounded for is closed and nowhere dense in the interval I.The Lebesgue measure of the set Γ is zero on I.For almost all h, the solutions of Equation (50) tends to periodic piecewise constat function with finite, countable or uncountable points of discontinuities on periods.

Example 1
Let 1 : where ( ) . Hence, solutions exists for any 0 t > .If 4 λ > , then all solutions of Equation (52) tends to infinity as t → ∞ .If , then ( ) +  , then D is countable set with with a finite set of points of accumulation.
In this case, we have deal with solutions of pre-turbulent type.If

λ >
, then D is the Cantor set and we have deal with solutions of turbulent type.

Reduction to Difference Equations
From the boundary conditions and the representation of solution in the form From the boundary conditions From the Neumann boundary conditions it follows that for a special choice of the initial conditions at the points Then functional Equations ( 57), ( 58) and ( 61), (62) can be written as: where and 1 , : . Then from (67), (68) it follows that We got the system of difference equations with continuous time.This system produce a map We assume that there is a bounded open subset where Φ is a continuous map such that: 1) differential DΦ is continuous on the set G; 2) a set ( ) , where φ is the empty set.
We will consider continuous solutions 69), (70).Such solutions are produced by the set of initial functions From these assumptions it follows that ( ) Φ for some natural N and ( ) we define N-periodic piecewise constant function ( ) where Then any solution of system is asymptotic periodic and such that for each fixed point t R + ∈ , where

Example 2
Consider the system where e , 0 a b a = < .This system produce the map Let us define the initial curve where the vector-function ( ) h t is determined by the initial data of the initial value boundary problem.
If 1 4 λ > the map λ ϕ has not fixed points and, hence, for each initial curve ( ) , solutions of the problem is such that m n is irreducible fraction).In this case, there are numbers k and i such that . Then, beginning from some number, we obtain a circle.For almost all (with respect of the Lebesque measure), the sequence is uniformly distributed in the interval.
There is a set Λ such that for almost u ∈ Λ trajectories { } < < =  the map λ ϕ has an attractive circle of the period 2 n , but all another circles are repelling.For the system of difference equations it means that ( ) u t tends to a 2 2 nl V -periodic piecewise constant function and ( ) w t tends to zero.

Physical Dynamic Boundary Conditions
Let us consider the following model [2]: , where ( ) is the strength of the nonlinearities, t I is the threshold intensity.We , where E is the electric field, * E is the value which is conjugate to the value E. Now we show how can be postulated nonlinear boundary conditions.For plane waves the density is defined as where c is velocity of light, l n is linear refractive index at low fields.The effective refractive index can be written as Indeed, as noted in [10]: "the nonlinear effects in optical fiber occur either due to intensity dependence of refractive index of the medium or due to inelastic-scattering phenomenon".There are various types of nonlinear effects based on self-phase modulation, cross-phase modulations and four-wave mixing.We consider the effect when the refractive index is the nonlinear function of the optical intensity.The power dependence describes the Kerr-effect [10].Origin of nonlinearity is the unharmonic motion of electrons under the action of the polarization  , which is induced by electric dipoles so that ( ) ( ) where 0 ε is the permittivity of vacuum and ( ) ,  is susceptibility.We assume that there is the surface polarization at a walls which confine the laser fiber.Then defects and color centers at the walls produce the second order generalization.We confined itself by the study only the quadratic boundary conditions which describes wave type distributions of relaxation, pre-turbulent and turbulent type.
Then the surface refractive index is ( ) . Then one of the boundary conditions can be written as 0, 0, , , , .
eff eff In similar form, we can consider more common boundary conditions.Then from (80) it follows that ( ) ( ) , χ χ are bifurcations parameters.Formally, their changing leads to appearing of period doubling distributions bifurcations and, correspondingly, to existence of distributions of relaxation, pre-turbulent and turbulent type.We can consider also the values Thus we obtain the problem with nonlinear boundary conditions which can be considered as above.Here, k l n is a bifurcation parameter.

Applications to Semiconductor Lasers
In semiconductor lasers, the intersection field matter is realized through the carrier density n which satisfies to the equation [3]: where Ĩ E , j is the injection current, q is the electron charge, v is active volume of the device, s τ is lifetime of the carriers, g is optical gain in cm −1 , and D is the diffusion coefficient.
We assume that device is placed at a point 0 x = , and the Neumann boundary condition Next, we assume that the media of semiconductor laser is ideal that is light propagates throw the media without changing of amplitude along lines d d where V is velocity of propagation of light.The we can consider phenomenological dynamic boundary condition We assume that ( ) where a n is the carrier density at transparency,  is the gain compress factor.As noted in [3]: "the actual dependence of the gain on the field density is still object of controversy".Note that in the paper [3] has been considered the case Thus we get the dynamic boundary conditions.Applying the method of reduction corresponding initial value boundary problem to a system of difference equations, we obtain that the problem has solutions of relaxation, pre-turbulent and turbulent type for a special case of the coefficients 1 1 , , a a n γ and so on.

Comparison with Experiment
In [9] has beeb discusses the dynamics of modulation surface directed instability and periodic waves in the coupled linear equations which describes light propagation in dispersed Kerr media.New spatial-temporal periodic solutions are found for these equations.As noted in [9]: there is "the fundamental link which exists between the phenomenon of polarization instability and the so-called black and white vector solitons".From this point of view, the asymptotic periodic spatial-temporal periodic solutions obtained in our paper are not the solitons, because solitons exists due to a balance between surface induced injection of radiation into optic media and volume radiation.Surface radiation with feedback can be modeled by laser diode.The balance can be achived only asymptotically as t → ∞ .We obtain white and black asymptotic impulse the form of which determined by the form of surface induced impulse.In this case, parameters of surface absorption and emission are bifurcation parameter which define the period and the frequency oscillations on the period white and black "solitons" [11].
Spatial solitons were experimentally observed in 1997 by Mitchel and Segev (see, [11]).Such solitons arise from incoherent white light and these solitons appears as a dynamic balance between the tendency for the beam to expend as a result of diffraction and the property of beam to contract as a result of self focusing.In [23] has been demonstrated solitons which can be obtained from both the spatially and temporally incoherent light.But, as noted in [11], the corresponding theory can not explain spatial-temporal coherence properties and also properties of temporal spectral density of such solitons.We believe that the initial value boundary problem, which is considered in the present paper, may be useful to describe formal properties of spatialtemporal solitons and their spectral properties.

Conclusion
In this paper, we analyze the dynamic of surface, and induce nonlinear instability in ideal optic volume producing light propagation, reflected from walls, which confines an ideal resonator.The problem is described by two linear difference equations with nonlinear dynamic boundary conditions for density of radiations and density of photons in the resonator.It is shown that deriving from the form of the boundary conditions with feedback we can obtain distributions of the light of relaxation, preturbulent and turbulent type.
)where α and β are positive or negative values.The diffusion coefficient n D , given that the charged particles are diffused and the Coulomb forces prevent diffusion can be associated with a carrier concentration by the formula: on parameters α and β , where ˆ1 α non-wandering points of the map Φ is finite and hyperbolic.Further, let ( ) FixΦ is a set of fixed points of the map Φ ; for all points of the map Φ we have that ( The set of non-wandering points of the map the interval I λ under an action of iterations of the map λ ϕ .Here, fixed point of the map λ ϕ .It means that each component of the solution tends to infinity as t → ∞ .The solutions are bounded if and only if 2

Λ
everywhere dense.The trajectories on Λ are unstable, but the set Λ are stable generally.It means that Λ attracts almost all trajectories from its neighbourhoods.For 2 λ = − the set Λ is the interval number of oscillations increase infinitely as t → ∞ .Such behaviour of trajectories exists not only for 2 λ = − , but for continuum values of λ .If 3 the fixed attractive point 1 β on this interval.It means that t → ∞ as t → ∞ .If 5 fixed point 1 β become repelling, but instead on Λ appears an attractive circle of the period 2 which consists from the two points 2

Figure 2 .,
Figure 2. Phase portrait for simplest solutions of relaxation type (white and black solitons).( ) ( ) 1 1 3 3 , , , p p p p -attractive − .Then the boundary conditions can be written as the other hand, from the of the injection locking it follows that are constants.Assuming that 2 1 I  at x l = , we obtain the right boundary condition Then a set Per α Φ