The aim of this work is to investigate projection operator method of deriva- tion of equations of motion for reduced density matrix and apply it to a model open quantum system. We gradually pass from quantum mechanical model of a molecule with one vibrational degree of freedom to an example of open quantum system relevant in the theory of nonlinear spectroscopy. In the thesis we present results of numerical simulations of the time evolution of the open quantum system performed with a program written for this purpose. We are specially concerned with simulations of the solution of the time-convolutionless generalized master equation up to the a second order of the perturbation expan- sion, and we show that under certain conditions it provides an exact solution of the problem. The text also contains derivation of the recurrence relations for the Franck-Condon factors for the most general case of two quantum harmonic oscillators in one space dimension, i. e. transformation matrix between two bases of the L2 (R) space determined by the solutions of the time-independent Schrödinger equation appropriate for these oscillators. 1
The aim of this work is to investigate projection operator method of deriva- tion of equations of motion for reduced density matrix and apply it to a model open quantum system. We gradually pass from quantum mechanical model of a molecule with one vibrational degree of freedom to an example of open quantum system relevant in the theory of nonlinear spectroscopy. In the thesis we present results of numerical simulations of the time evolution of the open quantum system performed with a program written for this purpose. We are specially concerned with simulations of the solution of the time-convolutionless generalized master equation up to the a second order of the perturbation expan- sion, and we show that under certain conditions it provides an exact solution of the problem. The text also contains derivation of the recurrence relations for the Franck-Condon factors for the most general case of two quantum harmonic oscillators in one space dimension, i. e. transformation matrix between two bases of the L2 (R) space determined by the solutions of the time-independent Schrödinger equation appropriate for these oscillators. 1
The dynamics of finite dimension open quantum systems is studied with the help of the simplest possible form of projection operators, namely the ones which project only onto one-dimensional subspaces. The simplicity of the action of the projection operators always leads to an analytical solution of the dynamical master equation, even in the non-Markovian case, in any perturbative order. The analytical solution correctly reproduces the short-time dynamics, and can be used to recursively recover the dynamics for an arbitrary time interval with arbitrary precision. The necessary number of relevant degrees of freedom to completely characterise an open quantum system is ( n − 1 ) ( n + 2 ) / 2 , where n is the dimension of the Hilbert space of the open system. The method is illustrated by two examples, the relaxation of a qubit in a thermal bath and the dynamics of two interacting qubits in a common environment.
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.
Our main goal in this paper is to extend to any system of coupled quadratic Hamiltonians some properties known for systems of quantum harmonic oscillators related with the Brownian Quantum Motion model. In a first part we get a rather general formula for the purity (or the linear entropy) in a short time approximation. In a second part we establish a master equation (or a Fokker-Planck type equation) for the time evolution of the reduced matrix density for bilinearly coupled quadratic Hamiltonians. The Hamiltonians and the bilinear coupling can be time dependent. Moreover we give an explicit formula for the solution of this master equation so that the time evolution of the reduced density at time t is connected with the reduced density at initial time t0 for $t_0 \leq t
Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.
The quantum molecular trajectory is the deterministic trajectory, arising from the Bohm theory, that describes the instantaneous positions of the nuclei of molecules by assuring the agreement with the predictions of quantum mechanics. Therefore, it provides the suitable framework for representing the geometry and the motions of molecules without neglecting their quantum nature. However, the quantum molecular trajectory is extremely demanding from the computational point of view, and this strongly limits its applications. To overcome such a drawback, we derive a stochastic representation of the quantum molecular trajectory, through projection operator techniques, for the degrees of freedom of an open quantum system. The resulting Fokker-Planck operator is parametrically dependent upon the reduced density matrix of the open system. Because of the pilot role played by the reduced density matrix, this stochastic approach is able to represent accurately the main features of the open system motions both at equilibrium and out of equilibrium with the environment. To verify this procedure, the predictions of the stochastic and deterministic representation are compared for a model system of six interacting harmonic oscillators, where one oscillator is taken as the open quantum system of interest. The undeniable advantage of the stochastic approach is that of providing a simplified and self-contained representation of the dynamics of the open system coordinates. Furthermore, it can be employed to study the out of equilibrium dynamics and the relaxation of quantum molecular motions during photoinduced processes, like photoinduced conformational changes and proton transfers.
Abstract This chapter introduces the density operator and the Liouville equation and outlines how reduced dynamical description of open quantum systems naturally leads to the appearance of relaxation processes. (1) The density operator and the quantum Liouville equation (the density matrix for a pure system, statistical mixtures, representations, coherences, thermodynamic equilibrium). (2) An example: The time evolution of a two-level system in the density matrix formalism. (3) Reduced descriptions. (4) Time evolution equations for reduced density operators: The quantum master equation (projection operators, The Nakajima–Zwanzig equation, Derivation of the quantum master equation using the thermal projector, The quantum master equation in the interaction representation, The Markovian limit—the Redfield equation, The secular approximation and the Lindblad equation). (5) The Lindblad equation in the ‘Heisenberg representation’. (6) Relaxation of a two-level system and a harmonic oscillator revisited. (7) The optically driven two-level system in a thermal environment—the Bloch equations. (8) The oscillating dipole model of a two-level molecule. (9) Analogy of a coupled two-level system to a spin ½ system in a magnetic field.