Abstract This chapter treats dissipation arising from coupling to the environment. Problems treated include photoionization, spontaneous emission, and cavity damping. A short-time perturbative analysis leads to decay of a single state coupled to a continuum at a rate determined by Fermi's golden rule. For longer times a non-perturbative solution of the amplitude equations leads to Weisskopf–Wigner decay. The final value theorem provides a simple method for obtaining the final state spectrum without first solving for the dynamics. A Heisenberg picture analysis leads to operator Langevin equations in which environment operators introduce quantum noise. A parallel analysis is possible working in the Schroedinger picture and leads to a master equation for the density operator of the damped system. This is solved to provide a complete description of the influence of dissipation on the system of interest. It is sometimes simpler to solve directly for the evolution of the field statistics.
Abstract This chapter treats dissipation arising from coupling to the environment. Problems treated include photoionization, spontaneous emission, and cavity damping. A short-time perturbative analysis leads to decay of a single state coupled to a continuum at a rate determined by Fermi's golden rule. For longer times a non-perturbative solution of the amplitude equations leads to Weisskopf–Wigner decay. The final value theorem provides a simple method for obtaining the final state spectrum without first solving for the dynamics. A Heisenberg picture analysis leads to operator Langevin equations in which environment operators introduce quantum noise. A parallel analysis is possible working in the Schroedinger picture and leads to a master equation for the density operator of the damped system. This is solved to provide a complete description of the influence of dissipation on the system of interest. It is sometimes simpler to solve directly for the evolution of the field statistics.
In the previous chapter, we developed the equation of motion for a system as it evolved under the influence of an unobserved (reservoir) system. We used the density matrix approach and worked in the interaction picture. In this chapter, we consider the same problem of the system-reservoir interaction using a quantum operator approach. We again eliminate the reservoir variables. The resulting equations for the system operators include, in addition to the damping terms, the noise operators which produce fluctuations. These equations have the form of classical Langevin equations, which describe, for example, the Brownian motion of a particle suspended in a liquid. The Heisenberg–Langevin approach discussed in this chapter is particularly suitable for the calculation of two-time correlation functions of the system operator as is, for example, required for the determination of the natural linewidth of a laser. We first consider the damping of the harmonic oscillator by an interaction with a reservoir consisting of many other simple harmonic oscillators. This system describes, for example, the damping of a single-mode field inside a cavity with lossy mirrors. The reservoir, in this case, consists of a large number of phonon-like modes in the mirrors. We also consider the decay of the field due to its interaction with an atomic reservoir. An interesting application of the theory of the system–reservoir interaction is the evolution of an atom inside a damped cavity. It is shown that the spontaneous transition rate of the atom can be substantially enhanced if it is placed in a resonant cavity.
In many problems in quantum optics, damping plays an important role. These include, for example, the decay of an atom in an excited state to a lower state and the decay of the radiation field inside a cavity with partially transparent mirrors. In general, damping of a system is described by its interaction with a reservoir with a large number of degrees of freedom. We are interested, however, in the evolution of the variables associated with the system only. This requires us to obtain the equations of motion for the system of interest only after tracing over the reservoir variables. There are several different approaches to deal with this problem. In this chapter, we present a theory of damping based on the density operator in which the reservoir variables are eliminated by using the reduced density operator for the system in the Schrödinger (or interaction) picture. We also present a ‘quantum jump’ approach to damping. In the next chapter, the damping of the system will be considered using the noise operator method in the Heisenberg picture. An insight into the damping mechanism is obtained by considering the decay of an atom in an excited state inside a cavity. The atom may be considered as a single system coupled to the radiation field inside the cavity. Even in the absence of photons in the cavity, there are quantum fluctuations associated with the vacuum state. As discussed in Chapter 1, the field may be visualized as a large number of harmonic oscillators, one for each mode of the cavity.
This chapter gives an overview of the field as it stands at the time of writing in mid 1995. It shows how the picture radically changes when there is more than one atom in the cavity. The simplest light-atom interaction one can think of is the coupling between a single atom and a single mode of the electromagnetic field. This system provides considerable insight into the fundamental physics of the light-atom coupling but the reader is warned that the more general many mode-many atom interactions are much richer in their dynamics. An atom in free space couples not only to one, but to a continuum of modes. The mathematics of the Weisskopf-Wigner theory of spontaneous emission tends to obscure the relatively simple physics of the emission process. An alternative, shorter way of calculating the single atom spontaneous emission rate is to use Fermi’s Golden rule, which is a prescription for calculating general transition rates to first order in perturbation theory.
We present the complete analytical solution of the cavity QED of a two-level atom and a field mode at zero temperature. It includes both dissipation of the field due to a finite Q value of the cavity and incoherent decay mechanisms for the atom. This analytical solution is provided by a powerful method for treating general master equations that appear in quantum optics. As distinct from the usual approaches we first deal with that part of the master equation which describes the dissipative coupling of the field and the atom to their reservoirs. Rather than using number-state or dressed-state bases we expand the density operator into the eigenstates of the nonunitary parts of the master equation which model the dissipative part of the dynamics. The set of these eigenstates is the damping basis. With the aid of this expansion we find the eigenvalues and eigenstates of the total Liouville operator. The evolution of an arbitrary initial state is then known. We employ these results to give an exact solution of the dynamics of the photon field in realistic experiments with one-atom masers at very low temperatures. It includes detuning, cavity leakage effects, spontaneous decay mechanisms for the atoms, a Fizeau-type velocity distribution for the atomic beam, and a statistical parameter for the probability of the excitation of incoming atoms, covering the limits of Poissonian pumping and of regular pumping. On the same grounds one can treat the one-atom laser, consisting of a single atom which stays in permanent interaction with the field mode and which is continuously pumped by external heat baths.
An investigation is made of the effect on the dynamics of spontaneous emission of the Wigner-Weisskopf atom in an infinite-system limit of various choices of the coupling function or form factor describing the atom's interaction with the spectrum of the radiation field. This is carried out both for the exact solution to the problem of spontaneous emission, obtained in the earlier papers in the series, and for some approximate solutions, also previously considered, in particular one based on the Schrödinger equation of the problem and one based on the weak-coupling Prigogine-Résibois master equation. The details of the form factor are found, by numerical computation of the solutions, to be critical in determining the nonexponential parts of the solutions, and these parts are seen to be capable in some cases of dominating the exponential parts, which are given only by the values of the form factor near the resonance energy. The approximate solutions discussed are found to vary widely in their worth, and one, which yields the exact solution for the Wigner-Weisskopf problem, is singled out as being of probable use in the statistical-mechanical description of more complicated systems.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.