In this thesis, the influence of an environment on molecules and, in particular, on the quantum control of such systems is investigated. Different approaches to describe system-bath dynamics are implemented and applied. The inclusion of a dissipation term in the system Hamiltonian leads to energy loss and relaxation to the ground state. As a first application, the isomerisation reaction in an aromatic complex is treated. It is shown that this simple model is able to reproduce results of time-resolved spectroscopic measurements. Next, the influence of noise is investigated. The incorporation of fluctuations reveals that energy is not conserved and coherences are destroyed. As an example, the quantum control of a population transfer in Na2 is examined. The efficiency of control processes is studied in dependence on the strength of the noise and different system-bath couplings. Starting with the unperturbed system, Local Control Theory is applied to construct a field which selectively transfers population into a single excited electronic state. The coupling to the bath is then switched on to monitor the dependence of the coupling strength on the transfer efficiency. The perturbation of the bath effects the Na2 molecule in such a way that potential energy curves and transition dipole moments are distorted. An important result is that already elastic collisions lead to a substantial loss of control efficiency. The most promising approach used in this thesis is the stochastic Schrödinger equation. It is equivalent to the commonly employed descriptions of system-bath dynamics within the reduced density matrix formalism. It includes decoherences and dissipation caused by elastic and inelastic collisions. Our contribution is the incorporation of laser excitation into the kinetic Monte-Carlo scheme. Thus we are able to apply this stochastic approach to the quantum control of population transfer in the sodium dimer. Because within our description it is possible to separate pure dephasing, inelastic transitions, and coherent time-evolution, we can identify the relative influence of these processes on the control efficiency. This leads to a far more physical picture of the basic processes underlying the perturbations of an environment then what a reduced density matrix description can provide. In utilising the stochastic wave function approach instead of the density matrix formalism, the computations are quite efficient. The stochastic Schrödinger equation is realised by N independent runs, where, in our case, an ensemble size of N = 1000 gives converged results. The efficiency of the laser control process is studied as a function of temperature and collision rates. A rise in temperature (or collision rate) reeffects a stronger fluctuation and thus results in a less efficient transfer by the control field. Though the Gaussian fluctuations used here do not strictly represent 'white'- noise, since a deterministic machine is not able to produce uncorrelated random numbers, an acceptable distribution is achieved by simple procedures. An improvement of the here applied algorithms would, for instance, include a more sophisticated sampling of the dephasing rates. Only one example of a control process is studied here and an application of the developed approach to other problems of quantum control is to be performed. This thesis established a systematic approach to understand quantum control in the presence of an environment.
In this thesis, the influence of an environment on molecules and, in particular, on the quantum control of such systems is investigated. Different approaches to describe system-bath dynamics are implemented and applied. The inclusion of a dissipation term in the system Hamiltonian leads to energy loss and relaxation to the ground state. As a first application, the isomerisation reaction in an aromatic complex is treated. It is shown that this simple model is able to reproduce results of time-resolved spectroscopic measurements. Next, the influence of noise is investigated. The incorporation of fluctuations reveals that energy is not conserved and coherences are destroyed. As an example, the quantum control of a population transfer in Na2 is examined. The efficiency of control processes is studied in dependence on the strength of the noise and different system-bath couplings. Starting with the unperturbed system, Local Control Theory is applied to construct a field which selectively transfers population into a single excited electronic state. The coupling to the bath is then switched on to monitor the dependence of the coupling strength on the transfer efficiency. The perturbation of the bath effects the Na2 molecule in such a way that potential energy curves and transition dipole moments are distorted. An important result is that already elastic collisions lead to a substantial loss of control efficiency. The most promising approach used in this thesis is the stochastic Schrödinger equation. It is equivalent to the commonly employed descriptions of system-bath dynamics within the reduced density matrix formalism. It includes decoherences and dissipation caused by elastic and inelastic collisions. Our contribution is the incorporation of laser excitation into the kinetic Monte-Carlo scheme. Thus we are able to apply this stochastic approach to the quantum control of population transfer in the sodium dimer. Because within our description it is possible to separate pure dephasing, inelastic transitions, and coherent time-evolution, we can identify the relative influence of these processes on the control efficiency. This leads to a far more physical picture of the basic processes underlying the perturbations of an environment then what a reduced density matrix description can provide. In utilising the stochastic wave function approach instead of the density matrix formalism, the computations are quite efficient. The stochastic Schrödinger equation is realised by N independent runs, where, in our case, an ensemble size of N = 1000 gives converged results. The efficiency of the laser control process is studied as a function of temperature and collision rates. A rise in temperature (or collision rate) reeffects a stronger fluctuation and thus results in a less efficient transfer by the control field. Though the Gaussian fluctuations used here do not strictly represent 'white'- noise, since a deterministic machine is not able to produce uncorrelated random numbers, an acceptable distribution is achieved by simple procedures. An improvement of the here applied algorithms would, for instance, include a more sophisticated sampling of the dephasing rates. Only one example of a control process is studied here and an application of the developed approach to other problems of quantum control is to be performed. This thesis established a systematic approach to understand quantum control in the presence of an environment.
We investigate the effect of inter-mode coupling on the vibrational relaxation dynamics of molecules in weak dissipative environments. The simulations are performed within the reduced density matrix formalism in the Markovian regime, assuming a Lindblad form for the system-bath interaction. The prototypical two-dimensional model system representing two CO molecules approaching a Cu(100) surface is adapted from an ab initio potential, while the diatom-diatom vibrational coupling strength is systematically varied. In the weak system-bath coupling limit and at low temperatures, only first order non-adiabatic uni-modal coupling terms contribute to surface-mediated vibrational relaxation. Since dissipative dynamics is non-unitary, the choice of representation will affect the evolution of the reduced density matrix. Two alternative representations for computing the relaxation rates and the associated operators are thus compared: the fully coupled spectral basis, and a factorizable ansatz. The former is well-established and serves as a benchmark for the solution of Liouville-von Neumann equation. In the latter, a contracted grid basis of potential-optimized discrete variable representation is tailored to incorporate most of the inter-mode coupling, while the Lindblad operators are represented as tensor products of one-dimensional operators, for consistency. This procedure results in a marked reduction of the grid size and in a much more advantageous scaling of the computational cost with respect to the increase of the dimensionality of the system. The factorizable method is found to provide an accurate description of the dissipative quantum dynamics of the model system, specifically of the time evolution of the state populations and of the probability density distribution of the molecular wave packet. The influence of intra-molecular vibrational energy redistribution appears to be properly taken into account by the new model on the whole range of coupling strengths. It demontrates that most of the mode mixing during relaxation is due to the potential part of the Hamiltonian and not to the coupling among relaxation operators.
The thesis comprises two major themes of quantum statistical dynamics. One is the development of quantum dissipation theory (QDT). It covers the establishment of some basic relations of quantum statistical dynamics, the construction of several nonequivalent complete second-order formulations, and the development of exact QDT. Another is related to the applications of quantum statistical dynamics to a variety of research fields. In particular, unconventional but novel theories of the electron transfer in Debye solvents, quantum transport, and quantum measurement are developed on the basis of QDT formulations. The thesis is organized as follows. In Chapter 1, we present some background knowledge in relation to the afore-mentioned two themes of this thesis. The key quantity in QDT is the reduced density operator ρ(t) ≡ trBρT(t); i.e., the partial trace of the total system and bath composite ρT(t) over the bath degrees of freedom. QDT governs the evolu-tion of reduced density operator, where the effects of bath are treated in a quantum statistical manner. In principle, the reduced density operator contains all dynam-ics information of interest. However, the conventional quantum transport theory is formulated in terms of nonequilibrium Green's function. The newly emerging field of quantum measurement in relation to quantum information and quantum computing does exploit a sort of QDT formalism. Besides the background of the relevant theoretical development, some representative experiments on molecular nanojunctions are also briefly discussed. In Chapter 2, we outline some basic (including new) relations that highlight several important issues on QDT. The content includes the background of nonequi-librium quantum statistical mechanics, the general description of the total compos-ite Hamiltonian with stochastic system-bath interaction, a novel parameterization scheme for bath correlation functions, a newly developed exact theory of driven Brownian oscillator (DBO) systems, and its closely related solvation mode trans-formation of system-bath coupling Hamiltonian in general. The exact QDT of DBO systems is also used to clarify the validity of conventional QDT formulations that involve Markovian approximation. In Chapter 3, we develop three nonequivalent but all complete second-order QDT (CS-QDT) formulations. Two of them are of the conventional prescriptions in terms of time-local dissipation and memory kernel, respectively. The third one is called the correlated driving-dissipation equations of motion (CODDE). This novel CS-QDT combines the merits of the former two for its advantages in both the application and numerical implementation aspects. Also highlighted is the importance of correlated driving-dissipation effects on the dynamics of the reduced system. In Chapter 4, we construct an exact QDT formalism via the calculus on path integrals. The new theory aims at the efficient evaluation of non-Markovian dis-sipation beyond the weak system-bath interaction regime in the presence of time-dependent external field. By adopting exponential-like expansions for bath corre-lation function, hierarchical equations of motion formalism and continued fraction Liouville-space Green's function formalism are established. The latter will soon be used together with the Dyson equation technique for an efficient evaluation of non-perturbative reduced density matrix dynamics. The interplay between system-bath interaction strength, non-Markovian property, and the required level of hierarchy is also studied with the aid of simple spin-boson systems, together with the three proposed schemes to truncate the infinite hierarchy. In Chapter 5, we develop a nonperturbative theory of electron transfer (ET) in Debye solvents. The resulting exact and analytical rate expression is constructed on the basis of the aforementioned continued fraction Liouville-space Green's func-tion formalism, together with the Dyson equation technique. Not only does it re-cover the celebrated Marcus' inversion and Kramers' turnover behaviors, the new theory also shows some distinct quantum solvation effects that can alter the ET mechanism. Moreover, the present theory predicts further for the ET reaction ther-modynamics, such as equilibrium Gibbs free-energy and entropy, some interesting solvent-dependent features that are calling for experimental verification. In Chapter 6, we discuss the constructed QDTs, in terms of their unified math-ematical structure that supports a linear dynamics space, and thus facilitates their applications to various physical problems. The involving details are exemplified with the CODDE form of QDT. As the linear space is concerned, we identify the Schrödinger versus Heisenberg picture and the forward versus backward propaga-tion of the reduced, dissipative Liouville dynamics. For applications we discuss the reduced linear response theory and the optimal control problems, in which the correlated effects of non-Markovian dissipation and field driving are shown to be important. In Chapter 7, we turn to quantum transport, i.e. electric current through molecular or mesoscopic systems under finite applied voltage. By viewing the nonequilibrium transport setup as a quantum open system, we develop a reduced-density-matrix approach to quantum transport. The resulting current is explicitly expressed in terms of the molecular reduced density matrix by tracing out the degrees of freedom of the electrodes at finite bias and temperature. We propose a conditional quantum master equation theory, which is an extension of the conven-tional (or unconditional) QDT by tracing out the well-defined bath subsets indi-vidually, instead of the entire bath degrees of freedom. Both the current and the noise spectrum can be conveniently analyzed in terms of the conditional reduced density matrix dynamics. By far, the QDT (including the conditional one) has only been exploited in a second-order form. A self-consistent Born approximation for the system-electrode coupling is further proposed to recover all existing nonlinear current-voltage behaviors including the nonequilibrium Kondo effect. Transport theory based on the exact QDT formalism will be developed in future. In Chapter 8, we study the quantum measurement of a qubit with a quantum-point-contact detector. On the basis of a unified quantum master equation (a form of QDT), we study the measurement-induced relaxation and dephasing of the qubit . Our treatment pays particular attention on the detailed-balance relation, which is a consequence of properly accounting for the energy exchange between the qubit and detector during the measurement process. We also derive a conditional quantum master equation for quantum measurement in general, and study the readout characteristics of the qubit measurement. Our theory is applicable to the quantum measurement at arbitrary voltage and temperature. A number of remarkable new features are found and highlighted in concern with their possible relevance to future experiments. In Chapter 9, we discuss the further development of QDT, aiming at an effi-cient evaluation of many-electron systems. This will be carried out by reducing the many-particle (Fermion or Boson) QDT to a single-particle one by exploring, e.g. the Wick's contraction theorem. It also results in a time-dependent density func-tional theory (TDDFT) for transport through complex large-scale (e.g. molecules) systems. Primary results of the TDDFT-QDT are reported. In Chapter 10, we summary the thesis, and comment and remark on the future work on both the theoretical and application aspects of QDT.
Classical Brownian motion has well been investigated since the pioneering work of Einstein, which inspired mathematicians to lay the theoretical foundation of stochastic processes. A stochastic formulation for quantum dynamics of dissipative systems described by the system-plus-bath model has been developed and found many applications in chemical dynamics, spectroscopy, quantum transport, and other fields. This article provides a tutorial review of the stochastic formulation for quantum dissipative dynamics. The key idea is to decouple the interaction between the system and the bath by virtue of the Hubbard-Stratonovich transformation or Itô calculus so that the system and the bath are not directly entangled during evolution, rather they are correlated due to the complex white noises introduced. The influence of the bath on the system is thereby defined by an induced stochastic field, which leads to the stochastic Liouville equation for the system. The exact reduced density matrix can be calculated as the stochastic average in the presence of bath-induced fields. In general, the plain implementation of the stochastic formulation is only useful for short-time dynamics, but not efficient for long-time dynamics as the statistical errors go very fast. For linear and other specific systems, the stochastic Liouville equation is a good starting point to derive the master equation. For general systems with decomposable bath-induced processes, the hierarchical approach in the form of a set of deterministic equations of motion is derived based on the stochastic formulation and provides an effective means for simulating the dissipative dynamics. A combination of the stochastic simulation and the hierarchical approach is suggested to solve the zero-temperature dynamics of the spin-boson model. This scheme correctly describes the coherent-incoherent transition (Toulouse limit) at moderate dissipation and predicts a rate dynamics in the overdamped regime. Challenging problems such as the dynamical description of quantum phase transition (local- ization) and the numerical stability of the trace-conserving, nonlinear stochastic Liouville equation are outlined.
Next, quantum-classical correspondence in intrinsic decoherence dynamics is examined via purities and reduced density matrices in both the position and energy representations. Dependence of the correspondence on the nature of the coupling potentials and on the initial conditions are discussed. Decoherence effects on quantum and classical dynamics in reactive scattering problems are investigated. Through the study of dynamics of the isotopic analogues of the collinear H + H2 reaction and transmission over simple one-dimensional barrier potentials, the role of decoherence in improving quantum-classical correspondence is identified with suppressing resonance and tunnelling. Finally, the efficiency of 1+3 photon phase control in the presence of decoherence effects phase is carefully examined in several model cases which include from a population inversion in a two-level bound system to IBr photodissociation. Further, the ways to overcome the control loss caused by dephasing are discussed.
The work is devoted to the investigation of the influence of a heat bath on the physical processes in a quantum system. We use the density matrix theory as one of the most powerfool tool for investigation of quantum relaxation. In the beginning of the work (chapter 2) we mention and recall the most important steps of derivation of the equation of motion for the reduced density matrix (master equation) for an arbitrary quantum system in diabatic representation interacting with the environment modeled by a set of independent harmonic oscillators. Chapter 3 deals with the question of the border between classical and quantum effects and reports on a study of the environmental influence on the time evolution of a coherent state or the superposition of two coherent states of a harmonic oscillator as a simple system displaying the peculiarities of the transition from quantum to classical regime. Chapters 4 and 5 concern the electron transfer (ET) problem, namely the mathematical description of the ET in molecular zinc-porphyrin-quinone complexes modeling artificial photosynthesis (chapter 4) and photoinduced processes in the porphyrin triad (chapter 5). Each chapter starts with an introduction and ends with a brief summary. The main achievements of the present work are summarized in the Conclusions.
All fields of physics - be it nuclear, atomic and molecular, solid state, or optical - offer examples of systems which are strongly influenced by the environment of the actual system under investigation. The scope of what is called "the environment" may vary, i.e., how far from the system of interest an interaction between the two does persist. Typically, however, it is much larger than the open system itself. Hence, a fully quantum mechanical treatment of the combined system without approximations and without limitations of the type of system is currently out of reach. With the single assumption of the environment to consist of an internally thermalized set of infinitely many harmonic oscillators, the seminal work of Stockburger and Grabert [Chem. Phys., 268:249-256, 2001] introduced an open system description that captures the environmental influence by means of a stochastic driving of the reduced system. The resulting stochastic Liouville-von Neumann equation describes the full non-Markovian dynamics without explicit memory but instead accounts for it implicitly through the correlations of the complex-valued noise forces. The present thesis provides a first application of the Stockburger-Grabert stochastic Liouville-von Neumann equation to the computation of the dynamics of anharmonic, continuous open systems. In particular, it is demonstrated that trajectory based propagators allow for the construction of a numerically stable propagation scheme. With this approach it becomes possible to achieve the tremendous increase of the noise sample count necessary to stochastically converge the results when investigating such systems with continuous variables. After a test against available analytic results for the dissipative harmonic oscillator, the approach is subsequently applied to the analysis of two different realistic, physical systems. As a first example, the dynamics of a dissipative molecular oscillator is investigated. Long time propagation - until thermalization is reached - is shown to be possible with the presented approach. The properties of the thermalized density are determined and they are ascertained to be independent of the system's initial state. Furthermore, the dependence on the bath's temperature and coupling strength is analyzed and it is demonstrated how a change of the bath parameters can be used to tune the system from the dissociative to the bound regime. A second investigation is conducted for a dissipative tunneling scenario in which a wave packet impinges on a barrier. The dependence of the transmission probability on the initial state's kinetic energy as well as the bath's temperature and coupling strength is computed. For both systems, a comparison with the high-temperature Markovian quantum Brownian limit is performed. The importance of a full non-Markovian treatment is demonstrated as deviations are shown to exist between the two descriptions both in the low temperature cases where they are expected and in some of the high temperature cases where their appearance might not be anticipated as easily.:1 Introduction. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 2 Theory of Open Quantum Systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 2.1 Influence Functional Formalism . . . . . . . . . . . . . . . . . . . . . . . . . 14 2.2 Quantum Brownian Limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 2.3 Stochastic Unraveling of the Influence Functional . . . . . . . . . . . . . . . 20 2.4 Improved Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 2.4.1 Modified Dynamic Response . . . . . . . . . . . . . . . . . . . . . . . 23 2.4.2 Guide Trajectory Transformation . . . . . . . . . . . . . . . . . . . . 24 2.5 Obtaining Properly Correlated Stochastic Samples from Filtered White Noise 24 3 Unified Stochastic Trajectory Propagation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.1 Semiclassical Brownian Motion . . . . . . . . . . . . . . . . . . . . . . . . . 28 3.1.1 Guide Trajectory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 3.1.2 Real Coherent State Center Coordinates . . . . . . . . . . . . . . . . 31 3.1.3 Propagation Scheme Including Stochastic Forces . . . . . . . . . . . 32 3.2 Stochastic Bohmian Mechanics with Complex Action . . . . . . . . . . . . . 33 3.2.1 Hydrodynamic Formulation of Bohmian Mechanics . . . . . . . . . . 33 3.2.2 Bohmian Mechanics with Complex Action . . . . . . . . . . . . . . . 34 3.2.3 Stochastic BOMCA Trajectories . . . . . . . . . . . . . . . . . . . . 38 3.3 Noise Distribution Preserving Removal of Adverse Samples . . . . . . . . . . 39 4 Dissipative Harmonic Oscillator. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 4.1 Reservoir Specification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 4.2 Harmonic Oscillator Analytic Expectation Values . . . . . . . . . . . . . . . 42 4.2.1 Ohmic Bath . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 4.2.2 Drude Regularized Bath . . . . . . . . . . . . . . . . . . . . . . . . . 43 4.3 Sampling Strategies and Analytic Comparison . . . . . . . . . . . . . . . . . 44 4.4 Limits of the Markovian Approximation . . . . . . . . . . . . . . . . . . . . 45 5 Dissipative Vibrational Dynamics of Diatomics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 5.1 Molecular Morse Potential . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 5.2 Anharmonic Phenomena . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 5.3 Transient Non-Markovian Effects . . . . . . . . . . . . . . . . . . . . . . . . 53 5.4 Trapping by Dissipation and Thermalization . . . . . . . . . . . . . . . . . . 53 6 Tunneling with Dissipation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57 6.1 Eckart Barrier Scattering . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58 6.2 Dissipative Tunneling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 6.3 Investigation of Markovianity . . . . . . . . . . . . . . . . . . . . . . . . . . 61 7 Summary and Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 Appendix A Conventions for Constants, Reservoir Kernels, and Influence Phases 69 Appendix B Stochastic Calculus. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 B.1 Stochastic Differential Equations . . . . . . . . . . . . . . . . . . . . . . . . 71 B.2 Position Verlet Scheme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72 B.3 Runge-Kutta Fourth Order Scheme . . . . . . . . . . . . . . . . . . . . . . . 73 Appendix CMorse Oscillator Expectation Values . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 Appendix DPrerequisites of a Successful Stochastic Propagation . . . . . . . . . . . . . . 79 D.1 Hubbard-Stratonovich Transformation . . . . . . . . . . . . . . . . . . . . . 79 D.2 Kernels of the Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 D.2.1 Quadratic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 D.2.2 Quartic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86 D.2.3 Strictly Ohmic Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . 89 D.3 Guide Trajectory Integration . . . . . . . . . . . . . . . . . . . . . . . . . . 90 D.3.1 Quadratic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90 D.3.2 Quartic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 D.3.3 Strictly Ohmic Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . 92 D.4 Computation of Matrix Elements and Expectation Values . . . . . . . . . . 92 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97