Lindblad Equation for Harmonic Oscillator: Uncertainty Relation Depending on Temperature

Specific nonequilibrium states of the quantum harmonic oscillator described by the Lindblad equation have been hereby suggested. This equation makes it possible to determine time-varying effects produced by statistical operator or statistical matrix. Thus, respective representation-varied equilibrium statistical matrixes have been found. Specific mean value equations have been found and their equilibrium solutions have been obtained.

Statistical operator  or density matrix is basically applied as the quantum mechanics tool, any information of the nonequilibrium process proceeding within the tested system may be gained from [1] [2] [3] [4] [5].When the process concerned proceeds within the system which fails interacting with its environment, statistical operator  will satisfy Liouville-von Neumann equa- tion as follows: With provision for the fact that the system interacts with any environment, a new equation shall be produced [4]- [18].Lindblad is the first one who offered the equation describing interaction of the system with a thermostat [11].This work is devoted to Markovian equation, which hereby describes nonequilibrium quantum harmonic oscillator performance.
We will write the kinetic equation for a quantum harmonic oscillator as follows: ( ) ( ) ˆˆˆˆˆˆˆˆˆˆˆˆˆ, , , , î i A a a H a a i B a a a a where ( ) A and B are constants.Operator â is formulated as follows: ( ) where m ω κ = .
Equation (1.2) is very precise to describe time varying state of the thermostat-interacted quantum harmonic oscillator and its equilibrium state.

Energy Representation
Now, we will define the wave functions describing specific energy state ( ) The very functions satisfy the equation as follows: ( ) ( ) where ( ) As referred to energy representation, the matrix elements of statistical operator  will be formulated by the equation as follows: Wave functions satisfy the following equations: With provision for the above formulas the following matrix-formed equation Now, we will write the equation for diagonal elements of density matrix , where n w is the probability referred to oscillator state n ϕ .The equa- tion produced has the form as follows: ( ) ( ) This kinetic equation describes particular harmonic oscillator state transitions.
In this case, there may be gained coefficients A and B as follows: ( ) Equation (2.6) has specific oscillator state equilibrium distribution, which satisfies the following equation: ( ) ( ) This equation is solved by the method as follows: ( ) under the following condition ( ) (2.10)

Mean Value of Coordinate
Mean value b assigned by operator b is defined as For gaining mean value a the respective equation may be derived from Formula (1.2).Using the equality of: we will get the equation as follows: ( ) Now, we can find the derivatives from mean values x and p .By applying Formula (1.4) we will get: ( ) ( ) Then, we will try to equate both the real and imaginary parts of this equation: If we eliminate p from this set of equations, we can obtain the mean coor- dinate equation The above Equation (3.5) provides the following solution: ( ) ( ) ( ) where 1 C and 2 C are arbitrary constants.

Mean Oscillator Energy
Now, we will find the time derivative from a a + .By applying the above equality (3.2) we will produce the following derivative from Equation (1.2): ( ) Thus, the following differential equation is derived: 2) The solution of the equation is: where C is an arbitrary constant.The Equation (4.2) has specific stationary solution: Since constants A and B are related (2.7), the stationary solution obeys the formula as follows:

Kinetic Equation Expressed in Terms of Coordinate and Momentum Operators
Let us express the Equation (1.2) in terms of operators x and p .For this purpose, we will firstly write the Equation (1.2) as follows: ) Since the energy operator is equal to: ( ) we will insert it in Equation (5.1) along with Formula (1.4) to obtain the following one:

Coordinate Representation
In coordinate representation the density matrix looks like this: ( ) The coordinate and momentum operators are: Using the above values we can write Equation (5.3) by the formula as follows: Physical interpretation of density matrix implies that the following expression is the probability density.
Let us introduce new variables ( ) In this case Referring to density matrix ( ) , , t x x  and using the above new variables we will get the equation as follows: In this case We will find the solution of Equation (6.4) as follows: Taking into account (6.6) we will obtain: Thus, in view of function (6.8) the following equation is formed: This equation has an equilibrium solution which satisfies both equations as follows: (6.12) We will write the performance equation of the above Formula (6.11): ( ) ( ) This formula implies that the general solution of Equation (6.11) takes the form as follows: ( ) where ( ) We put this function into Equation (6.12) to get the following formula: ( ) Taking into account condition (6.9) this equation has the following solution: ( ) ( ) Thus, the equilibrium solution of Equation (6.10) takes the form as We will find the equilibrium density matrix by formula (6.7): Integration brings us to the formula: where ( ) Using Formulas (6.3) we will get the equation as follows: Using Formula (6.2) we will get equilibrium probability density [19] ( ) ( )

Momentum Representation
In momentum representation the coordinate and momentum operators are: In this representation the density matrix looks like this: ( ) . This enables to write Equation (5.3) as follows: is the probability density to detect the state when an oscillator have impulse p.

Wigner Function
In order to better understand the physical meaning of various kinetic state summands we will derive the equation for Wigner function

( )
, , w w t x p = , which is a quantum analog of classical distribution function and can be defined with the use of density matrix ( ) by the relation: If the density depends on 1 x and 2 x , than ( ) ( ) ( ) ( ) , , , , exp d t x x w t x p ip x p Since there is Formula (6.9) , ,0 d From Equation (6.4) for density matrix ( ) , , t x x  we will get the equation for Wigner function The equation obtained is much different from its quantum analog of Fokker-Planck equation.Summands containing derivatives 2 x w ∂ and 2 p w ∂ can be interpreted as those describing phase space diffusion.And still, it is necessary to add that it is rather hard to find physical meaning of the formula in parentheses that follows coefficient ε .
The equilibrium solution of Equation (8.5) should at the same time be a solution for the following equations: The solution of this equation is the function: .11)Wigner function can be obtained by Formula (8.2) inserting in it equilibrium function (6.14).We have: This formula leads to the following result.If T = 0, the uncertainty is equal to

The Lindblad Equation Is a First Order Approximation
In work [17] it was proved that the Lindblad equation can be derived the quantum equation for a small system that interacts with the equilibrium system from the equation equations of Liouville-von Neumannas.The Lindblad equation can be written as ( ) here λ is the order parameter.Statistical operator ϱ write as well Let us substitute the operator (9.2) in Equation (9.1).We have ( ) If the order parameter λ = 0, we get the unperturbed statistical operator 0  : The first value of λ, gives ( ) Equilibrium values obey the equations: ( ) Examples of such equations are the Equations (6.11), (6.12) and (8.6), (8.7).

Conclusion
We considered the equation proposed by Lindblad for the statistical operator describing nonequilibrium state of quantum harmonic oscillator.From this equation, first we obtained the density matrix equation in energy representation and the equation for the diagonal elements of this matrix.We formulated the expressions defining physical meaning of Lindblad equation coefficients.Then we derived the equation for the mean value of a coordinate and found its general solution.We demonstrated that the mean coordinate value exponentially decreases in time.We obtained the equation for the mean oscillator energy and its general solution.We found the equilibrium mean energy value.This value is a monotonic decreasing function of temperature.We formulated Lindblad equation using coordinate and momentum operators.We obtained the density matrix equation in the coordinate representation.From this equation, we derived the formula for equilibrium density matrix.We wrote the density matrix equation in the momentum representation.We obtained Wigner function equation and found the respective equilibrium state function.We found the uncertainty relation for various temperatures by applying Wigner equilibrium function.