Open access peer-reviewed chapter

Nonequilibrium Thermodynamic and Quantum Model of a Damped Oscillator

By Gyula Vincze and Andras Szasz

Submitted: November 21st 2014Reviewed: June 9th 2015Published: December 21st 2015

DOI: 10.5772/61010

Downloaded: 867

Abstract

We describe the linearly damped harmonic quantum oscillator in Heisenberg’s interpretation by Onsager’s thermodynamic equations. Ehrenfest’s theorem is also discussed in this framework. We have also shown that the quantum mechanics of the dissipative processes exponentially decay to classical statistical theory.

Keywords

  • quantum oscillator
  • damping
  • Onsager’s theory
  • Ehrenfest’s theory
  • classic limit

1. Introduction

Dissipation is essential for the evolution of a quantum-damped oscillator. It is responsible for the decay of quantum states, the broadening of the spectral line, and the shifting the resonance frequency. This has been a persistent challenge for a long time since dissipation causes difficulties in the quantization of the damped oscillator [1, 2, 3]. This problem has remained under intensive investigation [4, 5]. There are some widely accepted Hamilton-like variation theories about the treatment of a linearly damped classic or quantum-damped oscillator. One of these theories is Bateman’s mirror-image model [1], which consists of two different damped oscillators, where one of them represents the main linearly damped oscillator. The energy dissipated by the main oscillator will be absorbed by the other amplified oscillator, and thus the energy of the total system will be conserved. The fundamental commutation relations of this model are time independent; however, the time-dependent uncertainty products, obtained in this way, vanish as time tends to infinity [6]. The Caldirola–Kanai theory with an explicit time-dependent Hamiltonian is another kind of variation theory [7, 8, 9]. In the quantum version of this theory, both the canonical commutation rules and the uncertainty products tend to zero as time tends to infinity. The system-plus-reservoir model [10, 11] is another damped oscillator model. It is coupled linearly to a fluctuating bath. If the bath is weakly perturbed by the system, then it can be modeled with a continuous bath of the harmonic oscillator. A quantum Langevin equation in the form of a Heisenberg operator differential equation can be deduced in this model. However, this equation in general does not obey Onsager’s regression hypothesis [12], i.e., only in case when 0[13]. A direct consequence of this fact is that the expected value of the fundamental observable does not satisfy the equation of the classic linearly damped oscillator. Another consequence is that no spontaneous dissipative process exists in this theory. The above models are based on the Heisenberg’s mechanical reinterpretation model [14].

A possible reinterpretation model based on irreversible thermodynamics was recently published [15]. This model started from the Rosen–Chambers restricted variation principle of the nonequilibrium thermodynamics [16, 17, 18] and used a Hamilton-like variation approach to the linearly damped oscillator. The usual formalisms of classical mechanics, such as the Lagrangian, Hamiltonian, and Poisson brackets, were also covered by this variational principle. By means of canonical quantization, the quantum mechanical equations of the linearly damped oscillator are given. The resulting Heisenberg operator differential equations of the damped oscillator are consistent with the classical equations of motion and can be solved by using ladder operators, which are time dependent. By this theory, the exponential decay of quantum states, the natural width of the spectral line, and the shifts in the resonance frequency can be explained. This work describes the quantum theory of a linearly damped oscillator, which could be reinterpreted in terms of a classical model based on Onsager’s nonequilibrium thermodynamic theory, corresponding to the Heisenberg reinterpretation principle. The first chapters are devoted to Onsager’s thermodynamic theory and the quantum theory of a damped oscillator. The dissipative quantum theory given in the Heisenberg picture is deduced from the general evolution equation of a Hermitian observable by means of two system-specific constitutive equations. The first of the constitutive equations belongs to unitary dynamics, while the second belongs to the dissipative dynamics of the observable. The fundamental commutators, which are a consequence of the constitutive equations, are time dependent. The quantum mechanical equations of motion of the oscillator in the Heisenberg picture, the Ehrenfest theorem, and the uncertainty principle of that oscillator are given. A significant part of the work deals with applications such as the expected value of the main operators of the damped oscillator, the probability description of the wave packet motion belonging to the damped oscillator, the calculation of the wave function by matrix calculus, the spectral density of the energy dissipation, and the natural width of the spectral line. Another significant part of this work deals the quantum statistics of the damped oscillator. By a generalization of the Liouville–von Neumann equation, the statistical thermodynamic theory of the ensemble of the damped oscillators in contact with a thermal bath is given. By introducing the quantum entropy of the ensemble, it is shown that the entropy of the ensemble grows in a dissipative process and in thermal equilibrium for the probability distribution of the quantum states, such that Gibbs’ canonical distribution is valid. Finally, a wave equation of the linearly damped oscillator is given.

2. Nonequilibrium thermodynamic theory of the linearly damped oscillator

Meixner was the first to propose a nonequilibrium thermodynamic theory for linear dissipative networks [19, 20]; for a general overview on network thermodynamics, see [21]. In this theory, it can be shown that, for example, electrical networks are thermodynamic systems, and it is possible to derive the network equations (Kirchhoff equations) by application of the principles of nonequilibrium thermodynamics. In what follows, we give an Onsagerian thermodynamic theory of the linearly damped harmonic oscillator. A damped oscillator, as a primitive network, is considered under the isotherm condition, which is maintained by removing the irreversible heat as it developed in damping resistance, or in other words, by placing the damping resistance of the oscillator in a temperature bath. In this case, it is possible to speak not only of the entropy of the damped oscillator but also of the free energy, and not of entropy production but of energy dissipation instead. To show these, we give the actual form of the first law of thermodynamics in the case of a damped oscillator. To do this, let us introduce from the energy conservation law of the oscillator + thermal bath system (Figure 1), such that we obtain the following:

ddt(U+Ubath)=0dUdt=dUbathdtE1

where U and Ubathare the internal energies of the oscillator and the thermal bath. Since the oscillator is isolated by a rigid diathermal wall, the energy exchange between the oscillator and the bath must be heat transfer only. Thus, the first law of thermodynamics of the oscillator and the thermal bath has the following forms

dUdt=dQdt,dUbathdt=dQdtE2

where dQdtis the power of exchanged heat on the rigid diathermal wall.

Figure 1.

Mechanical equivalent circuit of a linearly damped oscillator.

Assume that the entropy S of the oscillator exists, and it is a state function of the variables U, q. By the second law of nonequilibrium thermodynamics, the entropy S(U,q)satisfies the balanced equation:

dSdt=dSrdt+dSidt,dSidt0E3

where dSrdtis the reversible rate of change of the entropy or entropy flux on the rigid diathermal wall and dSidtis the so-called entropy production. From the second law of the thermodynamics follows the well-known expression for entropy flux:

dSrdt=1TdQdtE4

On the other hand, the rate of change in the entropy of the oscillator can be written as

dSdt=SUdUdt+SqdqdtE5

From this equation and Equations (2) and (4), the following relations can result:

SU=1T,dSidt=Sqdqdt0E6

Thus, the entropy balance equation has the form

dSdt=1TdQdt+SqdqdtE7

Because the temperature of the oscillator is constant, we could introduce

F=UTSE8

i.e., the free energy of the oscillator, and from Equations (2) and (7), we can obtain the simple balanced equation for free energy

dFdt=TSqdqdt0E9

Also, entropy production and the so-called rate of energy dissipation

R=TdSidt=TSqdqdt0E10

decrease the free energy of the oscillator. In what follows, we shall give the actual form of the rate of energy dissipation. Next, we see from Equation (10) that the rate of energy dissipation must be some explicit function of dqdtand may depend implicitly on U, q and the following equation

R=R(dqdt;U,q)0E11

such that

R(0;U,q)=0E12

is therefore quite general. We now expand the energy dissipation Equation (11) in a Taylor series, i.e.,

R=A0+A1dqdt+12A2(dqdt)2+.....E13

The sufficient condition of nonnegativity of entropy production, which will always be satisfied, is the complete exclusion of all odd terms in dqdtin Equation (13), and A0must be zero to exclude entropy production in an equilibrium state. Thus, one need only neglect the fourth-order term in the Taylor series (Equation (13)) to obtain the Rayleigh dissipation function

R=c(dqdt)2,c=12A2E14

In addition, the so-called damping constant of oscillator c must be positive to satisfy the nonnegativity condition in Equation (13). If we assume that the nondissipative elements of the damped oscillator are linear, then the free energy in Equation (8) can be identified with the energy stored in the mass m and the spring of the oscillator (Figure 1.). Also, the free energy is equal to the Hamiltonian of the oscillator, i.e.,

F=H=12p2m+12kq2E15

where q is the displacement of the mass from its equilibrium position, p is the momentum of the mass of the oscillator, and k is the constant of the spring. A direct consequence of the above results is that R can be constructed as a bilinear form, namely,

R=dHdt=dqdt(Hq)+dpdt(Hp)0E16

Here we now interpret, in the usual nonequilibrium thermodynamic fashion, the quantities (dqdt,dpdt)in terms of thermodynamic fluxes and the quantities (Hq,Hp)in terms of thermodynamic forces. Since we now have Equation (14) of the rate of energy dissipation of the damped oscillator, Equation (16) can be written in the form

dqdt(Hq)+dpdt(Hp)=c(dqdt)20E17

In Onsagerian thermodynamics, the constitutive (kinetic) equations between fluxes and forces are linear

[dqdtdpdt]=[LqqLqpLpqLpp][HqHp]E18

where the kinetic matrix

L=[LqqLqpLpqLpp]E19

can be split into nondissipative (so-called reactive) and dissipative parts. To do this, take these kinetic equations into Equation (17), and using Equation (15), with a simple calculation, we obtain for these kinetic matrices

[dqdtdpdt]=([0aa0]+[000c])[HqHp]E20

Here, a is an arbitrary constant. Now, we see that the dissipative part of the kinetic matrix satisfies the Onsager symmetry relation and the positivity of the damping constant c trivially. The constant a can be evaluated as follows. In the case of zero for the dissipative part of the kinetic matrix, these equations must be transformed into Hamilton equations of a simple harmonic oscillator. From this fact, it follows that a is a universal constant and a=1. Also, the final form of the Onsagerian constitutive (kinetic) equation the of damped oscillator is

[dqdtdpdt]=([0110]+[000c])[HqHp]E21

It is easy to show that these kinetic equations are equivalent to the Newtonian equations of motion of the linearly damped oscillator, namely,

dqdt=Hp=pm,dpdt=HqcHp=kqcpmE22

From this equivalence, it follows that the velocity of oscillator as a generalized thermodynamic flux has only a reactive part, while the rate of momentum, as another thermodynamic flux, has both reactive and dissipative constituents. The presented thermodynamic deduction of equations of a linearly damped oscillator enables us to build a stochastic force Fsinto equations of motion (Equation (22)). This stochastic force and the thermodynamic forces introduced above are statistically independent and take into account the effect of the temperature bath. In this case, the thermodynamic fluxes are fluctuations, which obey another type of equation (Equation (22)). The correct form of these equations follows from the regression hypothesis of Onsager [22], which states that “the average regression of fluctuations will obey the same laws as the corresponding macroscopic irreversible process.” From this, we give the Langevin-type equations for a linearly damped oscillator,

dqdt=Hp,dpdt=HqcHp+FsE23

A consequence of these equations is that the dissipative kinetic coefficient c can be related to the correlation coefficient p(t),p(t+t')via the fluctuation dissipation theorem [23, 24]. According to Equation (22), we can conclude that the fluctuations of thermodynamic fluxes are similar to an impressed macroscopic deviation, except they appear spontaneously.

2.1. Bohlin’s first integral

By means of Equation (21) or (22), we can deduce the time rate of change of any observable defined in the phase space of the damped oscillator. Let O(q,p,t)be an observable, such that its time rate of change can be expressed as

dOdt=Ot+Odqdqt+Opdpdt=Ot+OdqHpOpHqcOpHpE24

where we take into account the Onsagerian equations (Equation (22)). An observable O(q,p,t)is the first integral of the (Onsagerian equation (Equation (21)) of the damped oscillator if

dOdt=0E25

and if it is a constant of the motion. Now, let us give the constant of the motion of the linearly damped harmonic oscillator. It was Bohlin [25] who first dealt with the problem of the constants of motion for a damped linear oscillator. It is easy to prove that Bolin’s observable, defined as

B(q,p,t)=m2e2βt(dqdtγq)(dqdtγq)==e2βt(p22m+β2qp+β2pq+mω022q2)E26

where

β=c2m,ω0=km,γ=β+iω,γ=βiω,ω=ω02β2E27

is the first integral of the damped oscillator. Now, we may see that the Bolin’s observable in the case of the undamped oscillator is equal to the Hamiltonian of the oscillator.

3. Quantum theory of linearly damped oscillator

In the standard theory of quantum mechanics, two kinds of evolution processes are introduced, which are qualitatively different from each other. One is the spontaneous process, which is a reactive (unitary) dynamical process and is described by the Heisenberg or Schrödinger equation in an equivalent manner. The other is the measurement process, which is irreversible and described by the von Neumann projection postulate [26], which is the rigorous mathematical form of the reduction of the wave packet principle. The former process is deterministic and is uniquely described, while the latter process is essentially probabilistic and implies the statistical nature of quantum mechanics.

3.1. The general evolution equation of the Hermitian operator

Unlike classical quantum mechanics, the spontaneous processes of the damped oscillator are irreversible, so its quantum mechanical description needs changes to some instruments of classical quantum mechanics. To do this, we use the Heisenberg picture of quantum processes. In this picture, the observables are time-dependent linear Hermitian operators, and the state vector of the system is time independent. Using the terminology introduced in the first part, the infinitesimal time transformation of the Hermitian operator could happen in two ways:

  • By reactive transformation, when the orthonormal eigenvectors of the Hermitian observable turn in time, keeping the orthonormal system with unchanged eigenvalues. The eigenvectors belonging to the different moments are connected with unitary transformation, as in classic quantum mechanics. Dynamics belonging to this transformation are so-called unitary dynamics.

  • By dissipative transformation, when the real eigenvalues of the Hermitian operator change irreversibly in time.

Let us study the general evolution equation of the Hermitian operator, considering both the above time-dependent processes. For simplicity in demonstrating the derivation, we suppose a discrete eigenvalue spectrum of the Hermitian operator, although the spectra of the displacement and momentum operators could be continuous. In this case, the orthonormal eigenvectors |ΨOi(t)and eigenvalues λOi(t)are solutions of the eigenvalue equation

O(t)|ΨOi(t)=λOi(t)|ΨOi(t),(i=1,2,3,...)E28

and the eigenvectors form a complete orthonormal basis in a Hilbert space, when the eigenvalue spectrum is nondegenerate. Thus, the spectral representation of the operator would be

O(t)=O(t)δ=i0λOi(t)|ΨOi(t)Ψ0i(t)|,δ=i0|ΨOi(t)Ψ0i(t)|E29

where δis the unity operator and ΨOi(t)|is the dual of |ΨOi(t), so ΨOi(t)|ΨOi(t)=1.

Since the observable is Hermitian, the transformation of the eigenvector |ΨOi(0)at t=0to the eigenvector |ΨOi(t)at time t will be represented by an unitary operator U(t). Hence,

|ΨOi(t)=U(t)|ΨOi(0),ΨOi(t)|=U+(t)ΨOi(0)|,U(t)U+(t)=δE30

Let us substitute Equation (30) into Equation (29), then we obtain

O(t)=i0λOi(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)E31

The time derivative of the operator is

O(t)=i0λOi(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)++i0λOi(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)++i0λOi(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)E32

or by using the instantaneous eigenvectors, |ΨOi(t)can be written as

O(t)=i0λOi(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)++i0λOi(t)U(t)U+(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)++i0λOi(t)U(t)|ΨOi(0)Ψ0i(0)|U+(t)U(t)U+(t)==i0λOi(t)|ΨOi(t)Ψ0i(t)|+U(t)U+(t)(i0λOi(t)|ΨOi(t)Ψ0i(t)|)++(i0λOi(t)|ΨOi(t)Ψ0i(t)|)U(t)U+(t)=i0λOi(t)|ΨOi(t)Ψ0i(t)|+U(t)U+(t)O(t)+O(t)U(t)U+(t)E33

Using the identity U(t)U+(t)=U(t)U+(t), which is a consequence of unitarity, the general evolution equation of the Hermitian operator results

O(t)=O(t)t+[U(t)U+(t),O(t)],O(t)t:=i0λOi(t)|ΨOi(t)Ψ0i(t)|E34

where [,]is Dirac’s symbol of a quantum mechanical commutator and O(t)tis the local time rate of change of the operator in the coordinate system of the instantaneous eigenvectors. This equation is universal in the meaning of its independence of the constitutive behavior of the quantum system.

3.2. The Heisenberg equation of motion of the linearly damped oscillator

The actual form of Heisenberg’s dynamic equation can be constructed when the expression U(t)U+(t)is formed from the unitary operator and the local rate of change of the operator O(t)tcan be connected to the constitutive properties of the studied physical system. The first term will be ordered to the unitary/reactive and the second to the dissipative dynamics. To do this, we accept Heisenberg’s reinterpretation principle [14] (for the philosophical details of this principle, see [27]), which states the possibility of constructing a quantum mechanical description of a physical system whose classical description is known. In our case, the physical system is a linearly damped oscillator, for which we know its Onsagerian thermodynamic description. This description uses the Hamiltonian and the rate of energy dissipation of the system represented by the Rayleigh potential. Since the Hamiltonian is not the first integral of the motion, Bohlin’s first integral could be used for the unitary dynamics. In classical quantum theory, the Hamiltonian belongs to the unitary dynamics as a constitutive property, i.e.,

iH=U(t)U+(t)E35

where the Hamilton operator His a first integral of the system. In a closed physical system, the local time derivative of the observables is zero since the system is reactive. On this basis, the constitutive equations of classical quantum mechanics are

iH=U(t)U+(t),O(t)t=0E36

With these constitutive equations from Equation (34), the general evolution equation we give to Heisenberg’s equation of the observable is

O=i[H,O]E37

Also, the all constitutive properties of the quantum system are contained in the Hamilton operator only, which could have originated from the Hamilton function of the classical model by means of Heisenberg’s reinterpretation principle. Figure 2 shows the above-presented scheme of the deduction of Heisenberg’s equation of motion.

Figure 2.

The schema of the deduction of Heisenberg’s equation of motion of a Hermitian operator and the role of Heisenberg’s reinterpretation principle.

The Hamiltonian His a trivial nullifier of Dirac’s commutator in this approach, so His a conserved observable of motion, as was requested. To summarize, we get Heisenberg’s classical Equation (37) from the general evolution Equation (34), if the O(t)tlocal rate of change of the operator in the coordinate system of the instantaneous eigenvectors is zero. In the case of the Heisenberg’s equation (Equation (37)), the entropy of the system is constant in time, as proven by von Neumann [26]. However, the entropy cannot remain constant in dissipative processes. Consequently, in a correct description of the dissipative system, it is possible to take into account the local rate of change of the operator. In the case of a damped oscillator, this means that for the time rate of change q(t)t, p(t)tof two fundamental observables of a linearly damped oscillator must be given constitutive equations. In this way, we assume that the constitutive equations of the linearly damped oscillator are

q(t)t=βq(t),p(t)t=βp(t),iB=U(t)U+(t)E38

where the Hermitian operator Bbelonging to unitary/reactive dynamics is the quantum mechanical equivalent of Bohlin’s constant in Equation (26). Following Equation (26), this could be written as

B=e2βt(p22m+β2pq+β2qp+mω022q2)E39

Note, the third constitutive equation of (38) is the direct consequence of Stone’s theorem [28]. If we take into account these constitutive equations in the general evolution Equation (34) of the Hermitian operators, then we obtain Heisenberg’s equations of motion of a quantum-damped oscillator

q(t)=βq(t)+i[B,q(t)],p(t)=βp(t)+i[B,p(t)]E40

According to Heisenberg’s reinterpretation principle, these equations could be interpreted by means of the Onsagerian equations of the oscillator. To do this, split the Bohlin operator (Equation (39)) into two parts. The first part contains the Hamilton operator Hof the oscillator and the second part Dbelongs to the dissipation. We then obtain

q(t)=βq(t)+i[B,q(t)]={q(t)t+ie2βt[D,q(t)]}+ie2βt[H,q(t)]p(t)=βp(t)+i[B,p(t)]={p(t)t+ie2βt[D,p(t)]}+ie2βt[H,p(t)]+H=(p22m+mω022q2),D=β2(qp+pq)E41

The expression in {} is connected to dissipative thermodynamic current by analogy, while the currents outside the bracket are analogous to reactive currents. This interpretation, analogous to Equation (22), is supported by

q(t)t+ie2βt[D,q(t)]=0,p(t)t+ie2βt[D,p(t)]=2p(t)t=2βp(t)E42

where the two first equations of Equation (38) were used. In detail, we could write

qt+ie2βt[D,q]=βq+β2(ie2βt[p,q])u+β2q(ie2βt[p,q])=0pt+ie2βt[D,p]=βpβ2p(ie2βt[p,q])+β2(ie2βt[p,q])p=2βpE43

Now, we could see the desired interpretation analogy could be applied when the commutator relation e2βt[p,q]=iδis valid. Since every operator commutes with itself, we also have the fundamental brackets of our quantum theory

{p,q}=iδ,{q,q}=0,{p,p}=0,{,}:=e2βt[]E44

Consequently, the first fundamental bracket in Equation (44) ensures that the dissipative part β2qp+β2pqof the Bohlinian adds βq(t)term to the first equation and the βq(t)term to the second equation in Equation (41). This allows to us change our attention from Bohlinian to Hamiltonian in Equation (41), obtaining

q(t)=ie2βt[H(t),q(t)]=i{H(t),q(t)},p(t)=2βp(t)+ie2βt[H(t),p(t)]=2βp(t)+i{H(t),p(t)}E45

which are equivalent with the equations

q(t)=i{H(t),q(t)},p(t)=i{H(t),p(t)}ci{H(t),q(t)}E46

The quantum mechanical equations of a damped oscillator with the fundamental brackets in Equation (44), applying the rules of Lie algebra, are as follows

q(t)=i{H(t),q(t)}=p(t)m,p(t)=i{H(t),p(t)}ci{H(t),q(t)}=2βp(t)mω02q(t)E47

which are the operator differential equations version of Onsager’s equations in Equation (22).

To use the Lie algebraic method in an evaluation of the above-introduced time-dependent commutators, it is assumed that the scalar time functions must necessarily be considered as ordinary numbers (for details, see [15]). In summary, according to the Onsagerian equations of the damped oscillator by application of Heisenberg’s reinterpretation principle, the quantum mechanical equation of a damped oscillator in the Heisenberg picture can be obtained.

3.3. Ehrenfest theorem of a linearly damped oscillator

It is easy to show, similar to classical quantum mechanics, that the following operator analytics relations are valid [29]

i{H,q}=Hp,i{H,p}=HqE48

As a consequence of Equation (48), Equation (47), the quantum mechanical equations of the oscillator, could be written in the form

q(t)=Hp=p(t)m,p(t)=HqcHp=2βp(t)mω02q(t)E49

where the formal equivalence with Onsager’s equations Equation (22) is obvious.

It is well-known in the Heisenberg picture that the expectation value of an operator is defined as

O(t)=Ψ|O(t)|ΨE50

where |Ψis the time-independent state vector of the oscillator. Thus, from Equation (49), the expectation values of the time rate of change of displacement and momentum can be evaluated as

dq(t)dt:=Ψ|q|Ψ=iΨ|{H,q}|Ψ=Ψ|Hp|Ψ=1mΨ|p|Ψ=p(t)m,dp(t)dt:=Ψ|p|Ψ=iΨ|{H,p}|ΨciΨ|{H,q}|Ψ==Ψ|Hq|ΨcΨ|Hp|Ψ=mω02Ψ|q|ΨcmΨ|p|Ψ==cmp(t)mω02q(t)E51

where we take into account Equation (49). Also, the expectation values of displacement and momentum of the linearly damped oscillator obey time evolution equations, which are exactly equivalent to those of Onsager’s equations (Equations (21) and (22)). This result is Ehrenfest’s theorem.

4. Evaluation of the equations of the quantum linearly damped oscillator

The solutions of the operator differential equations (Equation (49)) are

q=2mω0(Aeiωt+A+eiωt),p=m2ω0(γAeiωt+γ*A+eiωt),A=aeβtE52

By substituting the above two expressions into the first fundamental commutation relation of Equation (44), the time-dependent and time-independent amplitude operators are used to obtain the following commutation relations

{A,A+}=e2βtδ[a,a+]=δE53

To solve the damped oscillator problem, we have to determine the operator Abecause this should be known for the specification of displacement, momentum, and the energy of the oscillator. In the case of a nondamped oscillator, the amplitude operator can be determined from the Hamilton operator of the oscillator, which is a constant of the motion. This is, however, not true for our case; thus, we will use the Bohlin operator introduced earlier. By substituting Equation (52) into the Bohlinian Equation (39), we get

B=e2βt(p22m+β2pq+β2qp+mω022q2)=ωm2e2βt(AA++A+A)==ωme2βt(A+A+12e2βt)=ωm(a+a+12),ωm=ω2ω01E54

where the commutation relations (Equation (44)) were used. Pursuant to the above two relations, it is easy to show that the time-independent amplitude operators fulfill the equations

[B,a+]=wma+,[a,B]=wmaE55

Now, we see that if we replace the operator Bby the Hamiltonian Hof a simple oscillator, these equations are identical to the corresponding equations of the simple quantum oscillator [30, 31]. According to this strong analogy, we are able to determine the amplitude matrix and the matrices of the Bohlin operator B, the displacement operator and the momentum operator. The results are as follows:

  • The operators a,a+,Band the occupation number operator N:=a+ahave the same eigenvectors and different eigenvalues. Since N:=a+ais positive definite, Bcan possess no negative eigenvalue. The lowest eigenvalue of Bbelongs to the eigenket |0of the operator afor which the relation a|0=0holds. From Equation (54), this so-called vacuum state belongs to the 12ωmzero-point Bohlinian eigenvalue and zero occupation number eigenvalue. The Bohlinian eigenvalue belonging to the eigenket |n(where the occupation number is n) can be calculated in the form of |n=a+nn!|0is (12+n)ωm, while the occupation number eigenvalue is n.

  • For the actions of the eigenket |nof the ladder operators, aand a+can be written as a|n=(n1)|n1,a+|n=(n+1)|n+1, from which it follows for the occupation number operator that N|n=n|n.

  • The matrices of the above-introduced operators are

a=[010.002.000.....],E56
N=aa+=[100.020.003.....].E57

on the basis of which is formed the orthonormal eigenkets

|0=[100.],|1=[010.],|2=[001.],etc.,E58

It is easy to see that the matrices that belong to aaand a+a+are not diagonal. From the above equations, it follows that the rules for the time-dependent ladder and occupation number operators are

A+(t)|n=eβtn+1|n,A(t)|n=eβtn1|n1,N(t)|n=e2βtn|nE59

5. Applications

5.1. Expected values of the main operators of a linearly damped oscillator

Moreover, the expected value of the occupation number in the nth energy eigenstate at time t is

N(t):=n|N(t)|n=N0e2βt,E60

where N0is the occupation number at t=0. This result agrees well with the corresponding result derived from the system-plus-reservoir model [32]. In the energy representation, since the matrices of the operators a,a+,a2,a+2have zero diagonal elements, the expected values of the operators A,A+,A2,A+2are zero in every energy eigenstate, i.e.,

A=0,A+=0,A2=0,A+2=0E61

According to these equations, the expected values q=:n|q|n,p=:n|p|nof the displacement and the momentum operator

q=2mω0(Aeiωt+A+eiωt),p=m2ω0(γAeiωt+γ*A+eiωt)E62

in the nth energy eigenstate are zero. The variance q2=n|q2|n,p2=n|p2|nof the displacement and momentum operator in the nth energy eigenstate can be evaluated as

q2=2mω0AA++A+A=2mω0δe2βt+2AA+==2mω0e2βt(1+2n)E63

and

p2=m2ω0γγ*(Aeiωt)(A+eiωi)+γγ*(A+eiωt)(Aeiωt)==mω02δe2βt+2(Aeiωt)(A+eiωt)==mω0e2βt(12+aa+)=mω0e2βt(12+n)E64

where we considered the commutation relation (Equation (53)). According to these results, we obtain the expected value of the energy of the damped oscillator

H(t)=p22m+mω02q22=ω0δe2βt+2AA+=ω0e2βt(n+12)E65

and the uncertainty relation

(Δq)(Δp)=2e2βt(2n+1),(Δq)(Δp)2e2βt(Δq)2:=(qδq)2=q2,(Δp)2:=(pδp)2¯=p2E66

where we considered that q=0,p=0. Now, we can see Heisenberg’s uncertainty relation is not fulfilled and, in the case of the simple oscillator and also when c=0, this relation transforms into Heisenberg’s relation.

5.2. Probability description of the wave packet motion of the damped oscillator

To learn something about the time dependence of our system in a certain state |, we will calculate q|(t)and |q||2(t), which represent the probability amplitude and probability of finding the damped oscillator at q at time t in that state. In particular, it is useful to study the so-called coherent state |a, which is an eigenstate of the non-Hermitian time-dependent operator A, i.e.,

Aeiωt|a=mω02A0eβteiωt|aE67

We shall also calculate the probability amplitude Ψa(q)=q|aof the wave packet |aat q. To do this, we shall start the following fact of bracket calculus

q|Aeiωt|q=q|Aeiωt|q'q'|qdq'=q|Aeiωt|q'Ψd(q')dq'==mω02A0eβteiωtΨa(q)E68

Now, we are going to express the operator Aby using the displacement and the momentum operator (52), so we get

Aeiωt=2mω0pmγqiωE69

Taking this expression into (68), we obtain

q|pmγ*qi|qΨa(q)dq=mωA0eβteiωtΨa(q)E70

According to the following

q|p|q'=e2βtidδ(qq')dq+f(q),q|q|q'=qδ(qq')E71

coordinate representation of the operators originating from the commutation relation Equation (44) (for details, see [15]), we get

(ωq+mddpe2βt)Ψd(q)=ωA0eβteiwtΨd(q)E72

an ordinary differential equation, where cqis chosen for the arbitrary function f(q). The solution of this equation in a normalized form is

Ψd(q)=g(t)(ωπ)14eωm2(qA0eβteiωt)2e2βtE73

where the function g(t)is evaluated in follows. From which the probability will be

|Ψa(q)|2=(mω2π)12eβtemωe2βt(qA0eβtcos(ωt))2E74

Now, we might see that this is the Gaussian distribution with the

|Ψa(q)|2=12π(mωeβt)e(qA0eβtcos(ωt))22(mωeβt)2E75

probability density function. Therefore, the motion of the center of the wave packet |ais a damped oscillation, and its uncertainty width decreases exponentially from the initial value of mωto zero (see Figure 3).

Figure 3.

The evolution of the wave packet |a⟩. The motion of the center of packet and its uncertainty width Δq are represented.

Now, we see that the initial uncertainty of the packet |akeeps getting smaller with the progression of time and becomes negligible as t. Also, the evolution of the wave packet continually proceeds toward the motion of a classic damped oscillator with the progression of time.

5.3. Calculation of wave function by matrix calculus

Resulting from Equation (52) using Equation (56), the matrix of the displacement operator in energy representation has the form

n'|q|n":=n'|A+A+|n"=2mω0[0100010200020300030400040etc.]eβtE76

Here, the displacement operator was used in a narrow sense. Next, we are going to solve the

q|q=q|qE77

eigenvalue problem in terms of the dn=n|qprobability amplitudes. To do this, we rewrite the above eigenvalue equation in matrix form

n0n|q|nn|q=qn|q,(n=0,1,2,)E78

where we take into account the n'0|n'n'|=δclosure relation. In explicit form, this looks like

[0100010200020300030400040etc.]eβt[d0d1d3...]=2mω0q[d0d1d3...]E79

From this, we get the difference schema

d1=qo1d0,d2=qo2d112d0,,dn=qondn1n1ndn2qo=2mω0eβtqE80

After some algebra, we obtain another form

d'n=2q2d'n12(n1)d'n2,d'n:=2nn!dnE81

By introducing a new coordinate variable, we have the difference equation

d'n=2qd'n12(n1)d'n2,q:=q2E82

This difference equation is satisfied by the Hermitian polynomials. Thus, we obtain

n|q=dn=c0(eβtmω0q)12nn!Hn(eβtmω0q),E83

where c0(eβtmω0q)is a function to be determined. Starting from the fact that q'|q''=δ(q'q''), we get

q|q=nq|nn|q=nc0(eβtmω0q)12nn!Hn(eβtmω0q)××c0*(eβtmω0q)12nn!Hn(eβtmω0q)=nc0(eβtmω0q)c0*(eβtmω0q)12nn!Hn(eβtmω0q)×Hn(eβtmω0q)=δ(qq)E84

By using the

n12nn!Hn(eβtmω0q')×Hn(eβtmω0q'')==πe(eβtmω0q')2δ(eβtmω0q'eβtmω0q'')==πmω0e(eβtmω0q')2eβtδ(q''q')E85

relationship, the final form of Equation (83) is given by

Ψn(q):=n|q=dn=e12βt(mω0π)14e12e2βtmω0q212nn!Hn(eβtmω0q)E86

The physical meaning of the strange variable in the Hermit polynomials is that the distance of the nodes of these functions keeps getting smaller with the progression of time by the exponential law eβt. According to this result, the probability density of the nth energy state with displacement q is

|Ψn(q)|2=|n|q|2=12nn!eβt(mω0π)12emω0q2e2βtHn2(eβtmω0q)=12π(2mω0eβt)eq22(2mω0eβt)2[12nn!Hn2(eβtmω0q)]E87

This result is exactly identical to the equation given by Kim and Page [33] on the basis of another theory. Now, we might see that this is the density function of a modulated Gaussian distribution, where the modulating term has finite amplitude which runs over in time, while the Gaussian distribution sharpens toward to a Dirac delta distribution. This means that the particle will get closer and closer to the equilibrium point as t. From last result, we can conclude that in the case of β0we get back to the well-known wave function of the simple oscillator.

5.4. Spectrum of the energy dissipation of the linearly damped oscillator

We are going to give the frequency spectrum of radiation and explain the natural width of the spectrum line. As an atom emits photons, its energy drops and the amplitude of transition decreases over time. Therefore, the emission is not harmonic, and a spectrum occurs. We shall see that the natural width of the spectral line can be connected to the attenuation coefficient of the damped oscillator. Inversely, from the width of the spectral line, we might determine the attenuation coefficient of the oscillator.

5.4.1. Spectral density of the energy dissipation

In the first section, the time rate of energy dissipation for a damped oscillator is introduced by the Rayleigh dissipation potential. The quantum version of this quantity, i.e., the time rate of the energy dissipation operator, can be originated from Equation (14) as

R:=cm1mp2=2βmp2E88

From this, it follows that the expected value of the operator of energy dissipation is

wdiss=0Rdt=0n|R|ndt=2βm0n|p|2ndtE89

Substituting the expression of the momentum operator (Equation (64)) into this equation, then the above equation has the form

wdiss=βω00(γγ*(Aeiωt)(A+eiωt)+γγ*(A+eiωt)(Aeiωt))dt==βω00(δe2βt+2(Aeiωt)(A+eiωt))dt=ω02+2βω0aa+0eγteγ*tdt=ω02+2βnω00eγteγ*tdt,γ=β+iωE90

where we assume that the occupation number is n, i.e., n=aa+. Evidently, the first term of this expression belongs to the vacuum fluctuation, and the second term belongs to the essential dissipative process of the damped oscillator in which the occupation number could change. We will evaluate the second term of this energy dissipation formula. According to the Parseval theorem of the Fourier transformation theory, this term of energy dissipation may be written as

Δwdiss(n)=2βnω00eγteγ*tdt=n2βω00F[eγt]F[eγt]*dω==n2βω00dω(ωω)2+β2E91

where F[eγt]is the Fourier transform of eγt. Also, the spectral density of the dissipated energy is

2βω0(ωω')2+β2E92

i.e., a Lorentz distribution about the shifted circular frequency ωof the damped oscillator. This result agrees well with corresponding result derived from the two-state atom model of Wigner and Weisskopf [34, 35] and the system-plus-reservoir model [10, 32, 36]. It is well known that the half value width Δωof this distribution is

Δω=βE93

Now, we can see the transition from the nth occupation number state to the vacuum state, in which the oscillator will emit nω0energy. Indeed, from Equation (91) it follows that

Δwdiss(n)=n2βω00dω(ωω)2+β2=nω0E94

5.4.2. Natural width of the spectral line

The natural line width of the spectral line is a significant result of the dissipative quantum process which accompanies the spontaneous emission of an atom. We will treat this emission process in a dissipative two-state model. We consider the two states of the atom as the zeroth and the first occupation number state of a linearly damped oscillator. In this case, the spontaneous emission of a photon is the consequence of the transition from the first occupations number state to the equilibrium state of the damped oscillator. In this model, the spectrum density of the emitted photon follows from Equation (92)

Δwdiss(ω')=2βω0(ωω')2+β2E95

The width of this frequency spectrum of a spontaneous emission of the atom is a direct consequence of the dissipative self-force on the atom due to the back-reaction of the emitted photon. This back-reaction of the emitted photon can be characterized by two physical quantities, namely, the frequency shift ω0ωand the half value width Δω=βof the spectrum. If we consider Δω2as the energy uncertainty ΔEof the emitted wave packet and the time constant of the emission process Δt=β1as the time uncertainty, we obtain an uncertainty relation

ΔEΔt2ββ=2E96

The quantum mechanical interpretation of the width of the natural spectral line should be based on this relation, in which the physical quantities ΔEand Δt=β1have a precise meaning. In our model, the natural line width occurs at wavelength λand can be calculated as

Δλ=|Δ2πcω|=2πcω2Δω=2πcω2βE97

where Equation (93) was used and c is the vacuum velocity of light. It is well known that in the classical dipole model of light emission, the natural line width can be calculated as

Δλ=4πε03mc2E98

where ε0is the vacuum permittivity and re:=ε03mc2=2,8181015mis the so-called classical electron radius. From the above two equations, it follows that

β=23ω2cre=4π3reλωE99

in the dipole radiation model.

6. Uncertainty relation of the linearly damped oscillator

The standard derivation of Heisenberg’s uncertainty relation neglects the possibility that two operators Aand B, say qand p,which fulfill the commutator relation

{A,B}=iE100

could have a compatible component which is the first part of the trivial identity

AB=AB+BA2+ABBA2E101

This observation has importance when we take into account the irreversibility. Due to irreversibility, the damped oscillator proceeds to thermal equilibrium with the thermal bath. This thermal equilibrium can be characterized in terms of classical statistic theory. However, in classical statistics, random variables have a joint distribution function, which could exist in the case of quantum theory if the operators are compatible. The commutator relation (Equation (100)) is compatible this physical picture, but from Equations (100) and (101), we obtain

AB=AB+BA2+δi2e2βtE102

From this relation, in the case of t, the compatibility of the operators follows, i.e.,

AB=AB+BA2AB=BAE103

In what follows, we will show that the above-mentioned arguments appear in the uncertainty relation. The variance of the Hermitian operators Aand Bcan be calculated by the norm of the following vectors

f=(AδA)|Ψ(0),g=(BδB)|Ψ(0)E104

Indeed, we can write

(ΔA)2=f2=Ψ(0)|(AδA)2|Ψ(0)(ΔB)2=g2=Ψ(0)|(BδB)2|Ψ(0)E105

where Ψ(0)is the state vector of the system, Aand Bare the expected value of the operators defined as

A=Ψ(0)|A|Ψ(0),B=Ψ(0)|B|Ψ(0)E106

Thus, the f2g2|f,g|2Schwarz inequality implies

(ΔA)2(ΔB)2|Ψ(0)|(AδA)(BδB)|Ψ(0)|2=|ABAB|2E107

By substituting into this expression the identity (Equation (101)), then we get

(ΔA)2(ΔB)2|AB+BA2+ABBA2AB|2==|AB+BA2+i2e2βtAB|2=24e4βt+(AB+BA2AB)2E108

where we take into account that the quadrate of the absolute value of a complex number is equally the sum of the quadrate of its real and imaginary parts. From the above expression, in the case of t, it follows that

ΔAΔBAB+BA2AB,ΔA:=+(ΔA)2,ΔB:=+(ΔB)2E109

On the another hand, in this case, the commutating relation (Equation (103)) is valid; thus, we can conclude that

(ΔA)(ΔB)ABAB0ABAB(ΔA)(ΔB)1E110

which is the most primitive “uncertainty relation” of classical statistic theory in which the random variables have a joint distribution function. It states the simple fact that the regression coefficient is smaller than one if the random variables are not statistically independent.

In summary, we can provide a speculative interpretation of irreversibility in quantum mechanics, namely, in an irreversible quantum process. The incompatible operators proceed to compatible ones, which are submitted to the laws of classical statistic theory.

7. Quantum statistics of the linearly damped oscillator

It was von Neumann [26] who first dealt with the problem of the quantum statistical ensemble. The density operator is the statistical operator of a quantum statistical ensemble. In our case, the statistical ensemble is a set of linearly damped oscillators of several quantum states in contact with a heat bath with temperature T. The density operator is an operator whose eigenvalues are the classical statistical probability of the chosen microstates denoted by pi. If the chosen microstates are denoted by |i, which are eigenstates of a Hermitian operator but not necessarily the eigenstate of a Bohlinian or Hamiltonian, the general density operator is written as

ρ:=ρδ=i0ρ|ii|=i0pi|ii|,δ=i0|ii|E111

From this definition, it follows that ρis Hermitian and normalized

ρ=ρ+,Trρ=1E112

In the Heisenberg picture, the density operator is time independent and is written as

ρH=i0pi(0)|i(0)i(0)|E113

The ensemble average of an operator in the Heisenberg picture AH(t)is defined as

A(t)=Tr(ρHAH(t))=i0pi(0)i(0)|AH(t)|i(0)E114

Ensemble averages of time rate of change of the displacement and the momentum of the linearly damped oscillator can be evaluated from Equation (51) as follows

dq(t)dt:=Tr(ρHq)=iTr(ρH{H,q})=Tr(ρHHp)=Hp=pm=p(t)mdp(t)dt:=Tr(ρHp)=iTr(ρH{H,p})ciTr(ρH{H,q})==Tr(ρHHq)cTr(ρHHp)=HqcHp=mω02qcpm==cmp(t)mω02q(t)E115

Here, Equation (50) was used. Now, we see that these equations are equivalent to those of the macroscopic Onsagerian equations (Equation (21) or (22)). In the Schrödinger picture, the density operator is time dependent, but the observables of the oscillator are time independent. We define this density operator as

ρS(t)=i0pi(t)|i(t)i(t)|E116

where we allowed a time-dependent probability pi(t)to microstate |i(t). Also in this picture, the occupation of microstates is not conserved. The ensemble average of an operator ASin the Schrödinger picture is defined as

A(t)=Tr(ρH(t)AS)=i0pi(t)i(t)|AS|i(t)E117

Two ensemble averages of an observable Amust be equal, i.e.,

Tr(ρS(t)AS)=Tr(ρHAH(t))E118

from which, in the case of pure unitary dynamics, follows the well-known transformation

ρS(t)=U+(t)ρH(t)U(t)E119

where the unitary operator U(t)belongs to time evolution [37]. Since, as we have seen in the case of dissipative processes, this cannot be written by unitary dynamics only, we will use this requirement in a weaker form. We require that the basic ensemble averaged equations of the damped oscillator have the same forms in each picture, i.e.,

dq(t)dt=Tr(ρHq)=Tr(ρSq)=p(t)mdp(t)dt=Tr(ρHp)=Tr(ρSp)=cmp(t)mω02q(t)E120

According to this requirement, we could give the actual form of the equation of motion for the density operator in the Schrödinger picture. We will see that this equation corresponds to the Liouville–von Neumann equation in the case of dissipative processes. From Equations (113) and (119), it follows that the density operator in the Schrödinger picture could be written by a Hermitian operator in the form

ρS(t)=i0pi(t)U+(t)|i(0)i(0)|U(t)E121

From the general evolution equation (Equation (34)) of the Hermitian operator, the equation of motion of the density operator in Schrödinger picture could be derived as follows:

ρS=dρS(t)dt=ρS(t)t+[U+(t)U(t),ρS]ρS(t)t:=i0dpi(t)dtU+(t)|i(0)i(0)|U(t)E122

where in the case of a damped oscillator, the unitary transformation belongs to the Bohlin operator of Equation (38), i.e.,

U+(t)U(t)=iBE123

Thus, the equation of motion of Schrödinger’s density operator is

ρS=dρS(t)dt=ρS(t)t+i[B,ρS]=ρS(t)t+i{H,ρS}+i{D,ρS}E124

Here, similar to the Heisenberg equations (Equation (41)), we introduced the Hamiltonian and the dissipation operator, by means of the commutator {,}=e2βt[,]. To construct a constitutive equation for the local change in the density operator, one must take into account the consequences of the commutation relation (Equation (44)) and the facts (Equations (102) and (103)) by which the incompatibility of the operators show strict fading over time. A consequence of this fading property could be increasing uncertainty in the distinction of the quantum states of the oscillators. On the other hand, in a canonical ensemble, the oscillators weakly interact with each other. By this means, the occupation of the states could change in the ensemble. Thus, we can assume that the occupation of states is not conserved in the time evolution of the ensemble. As a consequence, the statistical ensemble of the oscillators could proceed to a final state in which the classical probabilities of the microstates correspond to a classical thermal equilibrium distribution. Denoted by the density operator ρSin the instantaneous and ρSequin the equilibrium final state, then the suggested linear constitutive equation is

ρSt=β(ρSρSequ),Tr(ρS)=Tr(ρSequ)E125

Thus, the final form of the Liouville–von Neumann Equation (124) is

ρS=dρS(t)dt=β(ρSρSequ)+i[B,ρS]==β(ρSρSequ)+i{H,ρS}+i{D,ρS}E126

We will show that this evolution equation guarantees that the equivalence relation (Equation (120)) is fulfilled, the density matrix proceeds to an equilibrium state and that the entropy of the ensemble of the damped oscillator proceeds the maximum value over time, which corresponds to thermal equilibrium. Indeed, the proof of the relations in Equation (120) proceeds as follows

dq(t)dt=Tr(ρSq)=βTr((ρSρSequ)q)+iTr({H,ρS}q)+iTr({D,ρS}q)==βq(t)βqequ+iTr(ρS{H,q})+iTr(ρS{D,q})==βq(t)βqequ+Tr(ρSHp)+βTr(ρSq)=βqequ+p(t)mdp(t)dt=Tr(ρSp)=βTr((ρSρSequ)p)+iTr({H,ρS}p)+iTr({D,ρS}p)==βp(t)βpequ+iTr(ρS{H,p})+iTr(ρS{D,p})==βp(t)βpequTr(ρSHq)βTr(ρSp)==βpequcmp(t)mω02q(t)E127

where we take into account the cyclic invariance of the trace and the facts in Equation (42).

Now we see that if we choose an ensemble of a damped oscillator in which qequ=0,qequ=0, the required equivalence of the ensemble averages is fulfilled. To prove the increase in entropy, first we introduce quantum entropy. The kBρSlnρSoperator is an operator whose eigenvalues are the terms of Shannon entropy kBpi(t)lnpi(t). Thus, the Shannon entropy is the minus trace of that operator [26]

S(t):=kBi0pi(t)lnpi(t)=kBTr(ρSlnρS)E128

Here, kBis the Boltzmann constant. According to this definition, the excess entropy of an ensemble, as suggested by Bedeaux and Mazur [38], is given by

S(t)Sequ=kBTr(ln(δρSρSequ1)ρS),δρS=ρSρSequE129

The time rate of change of the entropy in that approximation is

dS(t)dt=12kBTr((δρSρSequ1+ρSequ1δρS)dρSdt)E130

Entropy production results by substituting the Liouville–von Neumann equation into this equation

dS(t)dt=12kBTr((δρSρSequ1+ρSequ1δρS)dρSdt)==12kBTr((δρSρSequ1+ρSequ1δρS)(βδρS+i{B,ρS}))==β2kBTr((δρSρSequ1+ρSequ1δρS)δρS)12kBTr((δρSρSequ1+ρSequ1δρS)(i{B,ρS}))==β2kBTr((δρSρSequ1+ρSequ1δρS)δρS)=βkBTr(ρSequ1(δρS)2)0E131

where the cyclic invariance of the trace and the fact that the Bohlin operator and the ρSequ1commutator were used. The positivity of entropy production follows from the fact that the matrices of both operators ρSequ1and (δρS)2have nonnegative elements only. Thus, the increase in entropy is demonstrated. It can be seen from the above deduction of entropy production that pure unitary dynamics (in the case of an undamped oscillator) is isentropic and that the entropy production is a direct consequence of the unconserved property of the occupation of states. In the thermal equilibrium, from Equation (126), it follows that Band ρSequcommute, i.e.,

dρS(t)dt=0,ρS=ρSequ,[B,ρSequ]=0E132

So the equilibrium density operator ρSequis a function of the Bohlinian B. At equilibrium, the entropy is at maximum. Now, we maximize S=kBn0pnlnpnunder the conditions

Tr(ρS)=n0pn=1,B=Tr(ρSB)=n0pnBn=const,Bn=n|B|n=ωm(12+n)E133

where we use Equation (54). The necessary condition of that maximum is

kBn0δpn(lnpn+1)=0E134

where the variations δpnare restricted by the conditions

n0δpn=0,n0δpnBn=0E135

Applying the method of Lagrange multipliers, we get

n0δpn[(lnpn+1)+βBn+γ]=0E136

So

pn=eβBnγ1E137

From the first equation of the conditions (133), the normalized version of the probability distribution is obtained

pn=eβBni0eβBiE138

Choosing β=1kBTas usual, then we get Gibbs’ canonical distribution for the occupation probabilities

pn=eβBni0eβBiE139

Introducing the partition function by definition

Z:=Tr(eβB)=i0eβBiE140

then we get the equilibrium density operator

ρSequ=eβBZE141

Thus, the equilibrium ensemble average of an operator AScan be written as

A=Tr(ρSequA)=Tr(eβBA)Z=n0n|A|neβEni0eβBiE142

In particular, for the ensemble average of the Bohlinian, this is

B=n0BneβEnn0eβBn=lnZβ=ωm(12+eωmkBT1eωmkBT)=ωm2cothωmkBTE143

Introducing the free energy by definition F=kBTZ, then in terms of free energy, the Gibbs distribution is written as usual

pn=eFBnkBTE144

Substituting this into the definition equation of entropy (128), then we get

S=BFTF=BTSE145

Thus, the ensemble average of the Bohlinian is the equilibrium internal energy. It is evident that the actual choice of the angular frequency ωmin the Bohlinian is a convention. It depends on the normalization of Bohlin’s constant (Equation (26)). What is the correct angular frequency? It seems from the physical aspect that the correct choice is that the angular frequency is ω0. In this case, Bohlin’s constant of motion corresponds to the maximum free energy of the linearly damped oscillator measured at time t=0, so it is the exergy of the linearly damped oscillator.

8. Wave equation of the linearly damped oscillator

From the above-presented theory, we can conclude that an ensemble from a pure state always proceeds to a mixed state a consequence of irreversibility. Thus, it is impossible to describe the evolution of the pure state of a damped oscillator in the Schrödinger picture. Consequently, it is impossible to construct a linear Schrödinger equation in which the position and the momentum operator are time independent.

However, when the operators are time dependent, the model could show similarities to Schrodinger’s interpretation, which we show below.

In the case of s linearly damped oscillator, the transformation of the Heisenberg picture into the Schrödinger picture by the method applied in classical quantum theory is impossible because the operator has a time-dependent part due to the dissipative process. Thus, a new way must be found to construct the wave equation of the oscillator. Kostin introduced a supplementary dissipation potential into his wave equation and constructed this dissipation potential by an assumption that the energy eigenvalues of the oscillator decay exponentially over time [39]. In Kostin’s version of the wave equation, the operators are time independent, but the dissipation potential is nonlinear with respect to the wave function. In our theory, it is assumed that the abstract wave equation of the linearly damped oscillator has the form

id|Ψ(t)dt=(HD)|ΨE146

where the Hamiltonian Hand the dissipative term D has the same mathematical form as in the Bohlinian, i.e.,

H(p,q)=p22m+mω022q2,D(p,q)=β2(pq+qq),E147

The operators p,qin this picture are time dependent and satisfy the classical fundamental commutators

[p,q]=iδ,[p,p]=[q,q]=0E148

The time derivative in Equation (146) is “material” (in the sense of continuum mechanics) because of the time dependence of the observable q^. To construct a wave equation, first we rewrite this abstract wave equation in the eigenbase |qof the position operator. To do this, consider the following one-dimensional eigenvalue problem

q|q'=q|q'E149

which could constitute a continuous spectrum. Thus, we must write the orthonormality condition and completeness relation for the eigenvectors as follows

q|q'=δ(qq'),|q'q'|dq'=δE150

where δ(qq')is Dirac’s distribution and δis the unity operator. The wave function is the probability amplitude that a position measurement on the damped oscillator in state |Ψ(t)will yield an eigenvalue q, mathematically

Ψ(q,t)=q|Ψ(t)E151

Inserting Equation (150) of the unity operator in the abstract wave Equation (146) and projecting from the left with q|, then we obtain

iq|a|Ψ(t)dt=dq|Ψ(t)dt=dΨ(q,t)dt=q|H(p,q)|q'q'|Ψ(t)dq'+q|D(p,q)|q'q'|Ψ(t)dq'E152

We chose the differential operator representation for the time-dependent operators in the eigenbase |qof the position operator in the form

p=ieβtq,q=qeβtE153

which resulted in the following wave equation

idΨ(eβtq,t)dt=(HD)Ψ(eβtq,t),H=22md2d(eβtq)2+mω022(eβtq)2,D=i(β2+β(eβtq)dd(eβtq))E154

which is a linear partial differential equation. To construct the eigenvalue problem that belongs to this wave equation, we chose the wave function as

Ψn=e(Eni+β2)tΦn(eβtq)E155

With this wave function, the eigenvalue equation

22md2Φnd(eβtq)2+mω022(eβtq)2Φn=EnΦnE156

is obtained from the wave equation because the equation

D(eβ2tΦn)=iβ(12+(qeβt)dd(qeβt))(eβ2tΦn)=ideβ2tΦndt==iβ2eβ2tΦn+ieβ2tΦn(qeβt)qdeβtdtE157

is satisfied identically for every eigenfunction Φn. By introducing a new variable into the eigenvalue Equation (156) defined as ξ:=mω0eβtq, then we get

d2Φndξ2+(2Enω0ξ2)Φn=0E158

We chose the eigenfunction Φnin the form

Φn=eξ22ψnE159

then we obtained the differential equation

d2ψndξ22ξψn+(2Enω01)ψn=0E160

which has a solution in terms of Hermitian polynoms if the

2Enω01=2nEn=(n+12)ω0,n=1,2,3...E161

relation is fulfilled. With these, the solution of the wave equation is obtained as follows:

Ψn(ξ,t)=e(Eni+β2)tΦn(eβtq)ei(n+12)ω0eβ2teξ22Hn(ξ),ξ:=mω0eβtqE162

from which the probability density function of the oscillator has the form

|Ψn(q)|2=12nn!12π(2mω0eβt)eq22(2mω0eβt)2[Hn2(eβtmω0q)]E163

which is exactly identical to Equation (87) resulting from the Heisenberg picture of the damped oscillator and the equation given by Kim and Page [33] using another theory. Due to this correspondence, the quantum decoherence of linearly damped oscillators could be described in the same way as done in the publication by Kim et al. [40].

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Gyula Vincze and Andras Szasz (December 21st 2015). Nonequilibrium Thermodynamic and Quantum Model of a Damped Oscillator, Recent Advances in Thermo and Fluid Dynamics, Mofid Gorji-Bandpy, IntechOpen, DOI: 10.5772/61010. Available from:

chapter statistics

867total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Linear Approximation of Efficiency for Similar Non- Endoreversible Cycles to the Carnot Cycle

By Delfino Ladino-Luna, Ricardo T. Páez-Hernández and Pedro Portillo-Díaz

Related Book

First chapter

Recent Studies on Fundamentals and Application of Microwave Processing of Materials

By Noboru Yoshikawa

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us