# Mathematical Models of Dissipative Systems in Quantum Engineering

Authors: Andreea Sterian; Paul Sterian
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

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 parametersain 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.


---

## Body

## 1. 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, coordinateq 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 parametersain 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 forN-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].


## 2. 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 coefficientsain in the Lindblad equation:
(2.1)ρ̇≡-iħ[H,ρ]+12ħ∑n{[Xnρ,Xn+]+[Xn,ρXn+]}depending on the system HamiltonianH and the operators of openingXn:
(2.2)Xn≡∑iainsi,wheresi are system operators, ain are complex coupling  coefficients or amplitudes and Γi operators of environment defined using the interaction Hamiltonian as(2.3)HSE=ħ∑isiΓi.These relationships have been established under the form [10](2.4)∑nainajn*=2ħ〈ΓiΓj〉,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.


## 3. 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]:(3.1)ddtρ(t)=-iħ[H,ρ(t)]+∑i,jλij{[ci+cjρ(t),cj+ci]+[ci+cj,ρ(t)cj+ci]}with dissipative coefficients:(3.2)λij=λijF+λijBincluding a componentλijF for a dissipative environment of fermions and a component λijB for a dissipative environment of bosons.Equation (3.1) is of Lindblad’s form, with dissipative operators depending on the transition/population operators ci+cj. For a system with N levels, the total number of these operators is N2-1 the number of the independent operators defined.If we denote byVF and VB the interaction dissipative potentials of the environment containing YF fermions and YB bosons, respectively, it is possible to write the expressions of the coefficients λijF şi λijB 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(εα) şi fβB(εβ). 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)λijF=πħYF∫|〈αi|VF|βj〉|2[1-fαF(εα)]fβF(εβ)gαF(εα)gβF(εβ)dεβ,εα-εβ=εj-εi,λijB=πħYB∫|〈αi|VB|βj〉|2[1+fαB(εα)]fβB(εβ)gαB(εα)gβB(εβ)dεβ,εα-εβ=εj-εi.


## 4. The Environment Dynamics Correlated with that of a Fermion Systems Coupled with Free Electromagnetic Field

We consider a system ofZ 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 two-particle 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] (4.1)HT=∑n=1Z(p̅n-eA̅B)22m+∑n=1ZU(1)(r̅n)+12∑n,m=1ZU(2)(r̅n,r̅m)+HB.In the total hamiltonian (4.1),(4.2)V=-em∑n=1Zp̅nA̅Bis the system-field interaction potential, while(4.3)HS=∑n=1Zp̅n22m+∑n=1ZU(1)(r̅n)+12∑n,m=1ZU(2)(r̅n,r̅m)is the fermion system hamiltonian, and(4.4)HB=∑νH(ν)is the field hamiltonian, where(4.5)H(ν)=ħων(aν+aν+12)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(4.6)ρ(t)=Tr⁡B{χ(t)}over the environment states.The total density operatorχ(t) satisfies the equation of motion:(4.7)dχ̃dt=-iħ[εṼR(t)+εṼ(t),χ̃(t)],where the sign aboveχ designs operators within the framework of interaction picture of the system and environment(4.8)χ̃(t)=e(i/ħ)(HB+H0S)tχ(t)e-(i/ħ)(H0S+HB)t,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:(4.9)χ̃(t)=R⊗ρ̃(t)+εχ̃(1)(t)+ε2χ̃(2)(t)+⋯,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(4.10)Tr⁡B{χ(1)(t)}=Tr⁡B{χ(2)(t)}=⋯=0,describe the environment dynamics that is correlated to the system dynamics. For an equation of motion of the form(4.11)dρ̃dt=εB̃(1)[ρ̃(t),t]+ε2B̃(2)[ρ̃(t),t]From (4.7), (4.9), and (4.11) we get a system of coupled equations:(4.12)R⊗B̃(1)[ρ̃(t),t]+dχ̃(1)dt=-iħ[ṼR(t)+Ṽ(t),R⊗ρ̃(t)],R⊗B̃(2)[ρ̃(t),t]+dχ̃(2)dt=-iħ[ṼR(t)+Ṽ(t),χ̃(1)(t)].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):(4.13)B̃(1)[ρ̃(t),t]=-iħTr⁡B[ṼR(t)+Ṽ(t),R⊗ρ̃(t)],B̃(2)[ρ̃(t),t]=-iħTr⁡B[ṼR(t)+Ṽ(t),χ̃(1)(t)],while, integrating by time, we get "excitation" terms of the total density operator (4.9):(4.14)χ̃(1)(t)=∫0t{-iħ[ṼR(t′)+Ṽ(t′),R⊗ρ̃(t′)]-R⊗B̃(1)[ρ̃(t′),t′]}dt′,χ̃(2)(t)=∫0t{-iħ[ṼR(t′)+Ṽ(t′),χ̃(1)(t′)]-R⊗B̃(2)[ρ̃(t′),t′]}dt′.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.


## 5. 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,Q12 is the field operator, P12 is the polarization operator, and N2 is population operator:(5.1)ddt〈Q12〉=-γ⊥〈Q12〉+ω21〈P12〉,ddt〈P12〉=-ω21〈Q12〉-γ⊥〈P12〉,ddt〈N2〉=-γ||[〈N2〉-N2(0)],with microscopic coefficientsγ⊥ şi γ|| expressed by dissipative coefficients λij of the master equation:(5.2)γ⊥=λ12+λ21+λ11+λ22,(5.3)γ||=2(λ12+λ21),(5.4)N2(0)=λ21λ12+λ21.The condition2γ⊥≥γ|| is a confirmation of master equation (4.7) which led to the establishment of Bloch equations-Feynman, because this condition is verified experimentally.


## 6. 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(6.1)dwdt=-γL(w+1)+μħ(Ẽ*Rab+ẼRab*),(6.2)dRabdt=-[γT+i(Δ+εw)]Rab-μ2ħẼw.In the above equations,w is the inversion of population, Rab nondiagonal elements of density matrix slow variable, indices a and b refer to lower and higher energy states, with the gap ħω0, EL 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εwRab, 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 γ||-1; this enables us to neglect dissipation processes). In these conditions, the matrix element Rab is decomposed into its real and imaginary parts Rab=1/2(ν+iu), resulting in system(6.3)dudt′=-(ετp)νw,dνdt′=(ετp)(u+Ωε)w,dwdt′=-(ετp)(Ωε)νwhose outcome is possible only numerically.In the above equationst′=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.In Figure1, we present the final population inversion function of maximum Rabi frequency for hyperbolic secant pulses Ẽ(t)=Ẽ0sech(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.Figure 1
Final state of population inversion, depending on the hyperbolic secant pulse maximum valueẼ(t)=Ẽ0sech(t/τp)(solid line), ετp=20 (continuous line) and ετp=30.In Figure2, 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).Figure 2
Temporal evolution of the system for a hyperbolic secant pulse,Ẽ(t)=Ẽ0sech(t/τp).


## 7. Conclusions

General expressions which relate the coupling parametersain 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 forN-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.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.



---

*Source: 347674-2012-05-24.xml*