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.
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.
In this paper, multidimensional dissipative quantum dynamics is studied within a system–bath approach in the Markovian regime using a model Lindblad operator. We report on the implementation of a Monte Carlo wave packet algorithm in the Heidelberg version of the Multi-Configuration Time-Dependent Hartree (MCTDH) program package, which is henceforth extended to treat stochastic dissipative dynamics. The Lindblad operator is represented as a sum of products of one-dimensional operators. The new form of the operator is not restricted to the MCTDH formalism and could be used with other multidimensional quantum dynamical methods. As a benchmark system, a two-dimensional coupled oscillators model representing the internal stretch and the surface–molecule distance in the O2/Pt(111) system coupled to a Markovian bath of electron–hole-pairs is used. The simulations reveal the interplay between coherent intramolecular coupling due to anharmonic terms in the potential and incoherent relaxation due to coupling to an environment. It is found that thermalization of the system can be approximately achieved when the intramolecular coupling is weak.
In this article, we explore the dissipation dynamics of a strongly coupled multidimensional system in contact with a Markovian bath, following a system-bath approach. We use in this endeavor the recently developed stochastic multi-configuration time-dependent Hartree approach within the Monte Carlo wave packet formalism [S. Mandal et al., J. Chem. Phys. 156, 094109 (2022)]. The method proved to yield thermalized ensembles of wave packets when intramolecular coupling is weak. To treat strongly coupled systems, new Lindblad dissipative operators are constructed as linear combinations of the system coordinates and associated momenta. These are obtained by a unitary transformation to a normal mode representation, which reduces intermode coupling up to second order. Additionally, we use combinations of generalized raising/lowering operators to enforce the Boltzmann distribution in the dissipation operators, which yield perfect thermalization in the harmonic limit. The two ansatz are tested using a model two-dimensional Hamiltonian, parameterized to disentangle the effects of intramolecular potential coupling, of strong mode mixing observed in Fermi resonances, and of anharmonicity.
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.
Effective descriptions accounting for the evolution of quantum systems that are acted on by a bath are desirable. As the number of bath degrees of freedom increases and full quantum simulations turn out computationally prohibitive, simpler models become essential to understand and gain an insight into the main physical mechanisms involved in the system dynamics. In this regard, vibrational decoherence of an I2 diatomics is tackled here within the framework of Markovian quantum state diffusion. The I2 dynamics are analyzed in terms of an effective decoherence rate, Λ, and the specific choice of the initial state, in particular, Gaussian wave packets and two-state superpositions. It is found that, for Markovian baths, the relevant quantity regarding decoherence is the product of friction (η) and temperature (T); there is no distinction between varying one or the other. It is also observed that decoherence becomes faster as the energy levels involved in the system state correspond to higher eigenvalues. This effect is due to a population redistribution during the dynamical process and an eventual irreversible loss of the initial coherence. These results have been compared with those available in the literature from more detailed semiclassical IVR simulations, finding a good agreement.
Modeling the dynamics of a quantum system coupled to a dissipative environment becomes particularly challenging when the system's dimensionality is too high to permit the computation of its eigenstates. This problem is addressed by introducing an eigenstate-free formalism, where the open quantum system is represented as a mixture of high-dimensional, time-dependent wave packets governed by coupled Schrödinger equations, while the environment is described by a multi-component quantum master equation. An efficient computational implementation of this formalism is presented, employing a variational mixed Gaussian/multiconfigurational time-dependent Hartree (G-MCTDH) ansatz for the wave packets and propagating the environment dynamics via hierarchical equations, truncated at the first or second level of the hierarchy. The effectiveness of the proposed methodology is demonstrated on a 61-dimensional model of phonon-driven vibrational relaxation of an adsorbate. G-MCTDH calculations on 4- and 10-dimensional reduced models, combined with truncated hierarchical equations for the mean fields, nearly quantitatively replicate the full-dimensional quantum dynamical results on vibrational relaxation while significantly reducing the computational time. This approach thus offers a promising quantum dynamical method for modeling complex system-bath interactions, where a large number of degrees of freedom must be explicitly considered.
In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.