Mathematical Models of Dissipative Systems in Quantum Engineering

The paper shows the results of theoretical research concerning the modeling and characterization of the dissipative structures generally, the dissipation being an essential property of the system with self-organization which include the laser-type systems also. The most important results presented are new formulae which relate the coupling parameters ain from Lindblad equation with environment operators Γi; microscopic quantitative expressions for the dissipative coefficients of the master equations; explicit expressions which describe the changes of the environment density operator during the system evolution for fermion systems coupled with free electromagnetic field; the generalized Bloch-Feynman equations for N-level systems with microscopic coefficients in agreement with generally accepted physical interpretations. Based on Maxwell-Bloch equations with consideration of the interactions between nearing atomic dipoles, for the dense optical media we have shown that in the presence of the short optical pulses, the population inversion oscillates between two extreme values, depending on the strength of the interaction and the optical pulse energy.


Introduction
An essential problem of the quantum information systems is the controllability and observability of the quantum systems.In this context, Fermi systems are essential for several important physical effects in quantum engineering as the dynamics of semiconductor nanostructures and high temperature superconductivity, nuclear resonances, fusion-fission reactions, and analysis of optical quantum systems.These effects are essentially determined by the dissipative coupling of the system.
Dissipation in quantum systems is a complex phenomenon which raises important theoretical investigations.A dissipative system is a system of interest, coupled with another system usually considered as being of much larger-environment.Fundamental and difficult problem of dissipative quantum theory is to design the total system system of interest environment on the space system of interest.In this way obtain a quantum master equation describing the evolution of the system using two terms: 1 a hamiltonian term for processes with energy conservation and 2 a nonhamiltonian term with coefficients that depend on the dissipative coupling.A master equation is based on approximations that consist in mediating rapid oscillations of reduced density matrix describing the interaction.
Such an approximation is the assumption that the evolution operators of a dissipative system forms a semigroup, not a group like for isolated systems.In this framework was derived a quantum master equation with dissipative terms which is consistent with all principles of quantum mechanics.Considering two operators, coordinate q and momentum p, master equation was used to describe the harmonic oscillator.In this theoretical framework, dissipation is described by the friction and diffusion coefficients that satisfy certain conditions called basic restrictions and Heisenberg's uncertainty relations are observed during the whole evolution of the system.
A rigorous method for deducting the master equation with microscopic expressions of the dissipative coefficients is developed in the literature.
For a weak dissipative coupling one obtains a master equation of Lindblad form 1 , but with the microscopic expressions of the dissipative coefficients.
In the development of quantum theory of dissipative systems an important step was the connection between Lindblad's generator and the previous phenomenological descriptions, realized by Sȃndulescu and Scutaru 2 .Besides, we must mention Isar et al.'s contributions 3 .This school developed by the above-mentioned researchers in the field are well recognized in the scientific world 4-8 .
Firstly, in the paper general expressions which relate the coupling parameters a in in Lindblad equation with environment operators Γ i have been established 9-11 .In this way, became possible deeper causality understanding of processes of friction and diffusion and of related quantum effects: broadening and shift of spectral lines, tunneling rates, bifurcations and instability 12, 13 .
Secondly, for a system of fermions, coupled with a dissipative environment quantitative microscopic expressions for the coefficients of the dissipative master equation depending on the potential matrix elements, the densities of states of the environment and the occupation probabilities of these states are presented 14-19 .
The study continue with the systems of fermions coupled by electric dipole interactions of free electromagnetic field for which has established general explicit expressions which describe the changes of the environment density operator during the system evolution.This description is not restricted to the Born approximation, taking into account the environment time evolution as a function of the system evolution.The results of the dissipative dynamics of the system of fermions in the presence of laser field are applicable to the dissipative structures 14, 20-29 .
Next, generalized Bloch-Feynman equations for N-level systems with microscopic coefficients in agreement with generally accepted physical interpretations are presented.
In the last part, we study the dynamics of dense media under the action of ultrafast optical pulses using Maxwell-Bloch formalism to include interaction between close atomic dipoles 30-34 .It is shown that, in a system initially without inversion, in the presence of optical pulses, the final population has two extreme values, results which contribute to understanding the specific mechanisms of switching for applications, with specific examples concerning the coherent radiation generation and amplification 35-43 .A computational specific software, to verify the experimental and numerical existing models and in the same time to discover new important situations for operative systems design and implementation, was developed 44-48 .

Relationship between Coupling Coefficients in Lindblad Master Equation and Environment Observables
Research on dissipative processes has led to evidence for the first time concerning the relationship between coupling coefficients a in in the Lindblad equation: depending on the system Hamiltonian H and the operators of opening X n : where s i are system operators, a in are complex coupling coefficients or amplitudes and Γ i operators of environment defined using the interaction Hamiltonian as These relationships have been established under the form 10 and allow an understanding of the physical causes of quantum processes of friction and diffusion, with their known effects: broadening and shift of spectral lines 11 , increased rates of tunneling, nonlinear characteristics, leading to bifurcation, instability, and chaos.

Microscopic Quantitative Expressions of the Dissipative Coefficients in Master Equations
A general quantum master equation for a many-level many-particle system, with microscopic coefficients, that preserves the quantum-mechanical properties of the density matrix was obtained 12 : with dissipative coefficients: including a component λ F ij for a dissipative environment of fermions and a component λ B ij for a dissipative environment of bosons.Equation 3.1 is of Lindblad's form, with dissipative operators depending on the transition/population operators c i c j .For a system with N levels, the total number of these operators is N 2 − 1 the number of the independent operators defined.
If we denote by V F and V B the interaction dissipative potentials of the environment containing Y F fermions and Y B bosons, respectively, it is possible to write the expressions of the coefficients λ F ij s ¸i λ B ij for the resonant transition |j → |i of the system coupled with |β → |α environmental transition, with fermionic state having densitiesaly g F α , g F β and populations f F α ε α , f F β ε β , and bosonic states with densities g B α , g B β and population f B α ε α s ¸i The probability that the final state |α of the environment to be free is 1 − f ε α while the probability the initial state of the environment to be occupied is f ε β .
General expressions of the dissipative coefficients are written for this type of interaction in the form: 3.3

The Environment Dynamics Correlated with that of a Fermion Systems Coupled with Free Electromagnetic Field
We consider a system of Z charged fermions with the coordinates r n and momenta p n n 1, 2, . . ., Z in a single-particle potential U 1 r n , while U 2 r n , r m represents the twoparticle residual potential.This system is coupled to the modes ν of the free electromagnetic field.In order to describe the dynamics of this system, for simplicity, we neglect the particle spin and its dimensions with respect to the electromagnetic field wavelength the electric dipole approximation .In this case, the total hamiltonian is of the form 14 In the total hamiltonian 4.1 , is the system-field interaction potential, while is the fermion system hamiltonian, and is the field hamiltonian, where is the field mode ν hamiltonian.
Let us take the density operator χ t of the total system with hamitonian 4.1 and the reduced density matrix ρ t Tr B χ t 4.6 over the environment states.The total density operator χ t satisfies the equation of motion: where the sign above χ designs operators within the framework of interaction picture of the system and environment χ t e i/ħ H B H S 0 t χ t e − i/ħ H S 0 H B t , 4.8 while ε is an intensity parameter used to show the orders of the series expansion of this density.Considering the radiation field of the black body in the initial state R, the total density operator of the system can be taken under the form: where χ 1 t , χ 2 t represent modifications of the field during the system evolution.The first term of this expression corresponds to the Born approximation when the environment state is a constant state R, while the higher-order terms, which satisfy the normalization relations describe the environment dynamics that is correlated to the system dynamics.For an equation of motion of the form Mathematical Problems in Engineering From 4.7 , 4.9 , and 4.11 we get a system of coupled equations:

4.12
By calculating the partial traces over the environment states and using the normalization conditions 4.10 , from these equations we get successively the terms of the equation of motion 4.11 : while, integrating by time, we get "excitation" terms of the total density operator 4.9 :

4.14
The first-order equation 4.13 represents the system evolution when the environment is considered as being in a constant state R, while for the higher-order term 28 , we take into consideration some changes of the environment matrix 4.14 .Further on, we will show that the first-order terms 4.13 describe the hamiltonian dynamics of the system, while the second-order term 28 describes system one-particle transitions related to environment.

The Generalized Bloch-Feynman Equations
An alternative description of dissipative system dynamics is given by Bloch-Feynman equations for systems of fermions obtained by defining the pseudo-spin operators 14 .
In particular, for a system with two-level known Bloch-Feynman, equations are obtained, where, Q 12 is the field operator, P 12 is the polarization operator, and N 2 is population operator: .

5.4
The condition 2γ ⊥ ≥ γ || is a confirmation of master equation 4.7 which led to the establishment of Bloch equations-Feynman, because this condition is verified experimentally.

Dynamics of Dense Media under the Action of Short Optical Pulses
Maxwell-Bloch equations of a two-level atomic medium generalized to include interactions between the dipoles approach 15, 16, 32 have been used to describe the system dynamics under the action of ultrafast optical pulses.These equations, for systems with homogeneous broadening of spectral lines in about semiclassical treating, were established using the density matrix formalism as In the above equations, w is the inversion of population, R ab nondiagonal elements of density matrix slow variable, indices a and b refer to lower and higher energy states, with the gap ħω 0 , E L is slowly varying local field Δ ω 0 − ω is the frequency deviation in relation to the center frequency of the field resonance frequency, μ is the transition matrix element of the electric dipole, and γ || , γ ⊥ are longitudinal and transverse relaxation rates.Contributions of the dipole-dipole interactions occur in 6.2 by term iεwR ab , where ε nμ 2 /3ħε d ω 0 is the strength parameter of dipole-dipole interactions having a dimension of a frequency.Equations 6.1 and 6.2 for the atomic variables and for field variables realise the description Maxwell-Bloch of optically dense environment.These equations were generalized and used to study intrinsic optical bistability, propagation effects in nonlinear media, and so forth.
For numerical simulation, we considered the case resonant Δ 0 , a characteristic distance between dipoles much smaller than the wavelength of the central field propagation effects are negligible and ultrafast pulses pulses much shorter than γ − In Figure 1, we present the final population inversion function of maximum Rabi frequency for hyperbolic secant pulses E t E 0 sech t/τ p .As long as the Rabi frequency has a value so that Ω 0 /ε < 1, the final population inversion is w −1.In the region Ω 0 /ε > 1, the final population inversion has an oscillatory behavior, almost rectangular wave.As the parameter ετ p value is greater, the oscillation period decreases, the transitions become abrupt, and the first half cycle of the rectangular wave becomes more centered to Ω 0 /ε 1.In Figure 2, temporal evolution of the system is presented for a hyperbolic secant pulse with a peak higher than one when t → ∞, the population inversion performs a number of oscillations before reaching a value 1; under certain conditions when t → ∞, after a number of oscillations, the system remains in the ground state .

Conclusions
General expressions which relate the coupling parameters a in in Lindblad equation with environment operators Γ i have been established.These expressions allow deeper understanding of causal processes of friction-and diffusion-related quantum effects: broadening and shift of spectral lines, tunneling rates, bifurcations, and instability.For a system of fermions coupled with a dissipative environment quantitative microscopic expressions for the coefficients of the dissipative master equation are presented.
These coefficients depend on the potential matrix elements, the densities of states of the environment, and the occupation probabilities of these states.
Expressions of the dependence of the particle distributions on temperature are taken into account.It can be shown that a system of fermions located in a dissipative environment of bosons tends to a Bose-Einstein distribution.
Studying the systems of fermions coupled by electric dipole interactions of free electromagnetic field, has established general explicit expressions which describe the changes of the environment density operator during the system evolution for fermion systems coupled with free electromagnetic field.This description is not restricted to the Born approximation, taking into account the environment time evolution as a function of the system evolution.The study can be continued with the calculation of the higher-order term of the reduced matrix equation in order to describe the correlated transition of the system particles.The results of the dissipative dynamics of the system of fermions in the presence of laser field are applicable to the dissipative structures.
Generalized Bloch-Feynman equations for N-level systems with microscopic coefficients in agreement with generally accepted physical interpretations are presented.On this basis, the problem of a quantum system control is explicitly formulated in terms of microscopic quantities: matrix elements of the dissipative two-body potential, densities of the environment states, and occupation probabilities of these states.

Mathematical Problems in Engineering
Studying the dynamics of dense media under the action of ultrafast optical pulses using Maxwell-Bloch formalism to include interaction between close atomic dipoles showed that, in a system initially without inversion, in the presence of optical pulses, the final population has two extreme values, the ratio of Rabi frequency and the parameter that describes the interactions between close dipoles, which contribute to understanding the specific mechanisms of switching.

Figure 2 :
Figure 2: Temporal evolution of the system for a hyperbolic secant pulse, E t E 0 sech t/τ p .
⊥ s ¸i γ || expressed by dissipative coefficients λ ij of the master equation:γ ⊥ λ 12 λ 21 λ 11 λ 22 , Final state of population inversion, depending on the hyperbolic secant pulse maximum value E t E 0 sech t/τ p solid line , ετ p 20 continuous line and ετ p 30.neglect dissipation processes .In these conditions, the matrix element R ab is decomposed into its real and imaginary parts R ab 1/2 ν iu , resulting in system In the above equations t t/τ p is the normalized time, τ p is the measured width pulse, Ω t μ E t /h is the instantaneous Rabi frequency, and E t is the intensty of electrical pulse.