Time-convolutionless master equation for mesoscopic electron-phonon systems

The time-convolutionless master equation for the electronic populations is derived for a generic electron-phonon Hamiltonian. The equation can be used in the regimes where the golden rule approach is not applicable. The equation is applied to study the electronic relaxation in several models with the finite number of normal modes. For such mesoscopic systems the relaxation behavior differs substantially from the simple exponential relaxation. In particular, the equation shows the appearance of the recurrence phenomena on a time scale determined by the slowest mode of the system. The formal results are quite general and can be used for a wide range of physical systems. Numerical results are presented for a two level system coupled to Ohmic and super-Ohmic baths, as well as for a model of charge-transfer dynamics between semiconducting organic polymers.     

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.     

Stochastic surrogate Hamiltonian

The surrogate Hamiltonian is a general scheme to simulate the many body quantum dynamics composed of a primary system coupled to a bath. The method has been based on a representative bath Hamiltonian composed of two-level systems that is able to mimic the true system-bath dynamics up to a prespecified time. The original surrogate Hamiltonian method is limited to short time dynamics since the size of the Hilbert space required to obtain convergence grows exponentially with time. By randomly swapping bath modes with a secondary thermal reservoir, the method can simulate quantum dynamics of the primary system from short times to thermal equilibrium. By averaging a small number of realizations converged values of the system observables are obtained avoiding the exponential increase in resources. The method is demonstrated for the equilibration of a molecular oscillator with a thermal bath.     

Correction of apparent finite size effects in the area per lipid of lipid membranes simulations

Molecular dynamics simulations of lipids bilayers have reported that the average area per lipid increases with the size of the simulated unit cell under constant temperature, pressure, and number of molecules. Here we show that the cause of this finite size effect are artifacts associated with the heat bath coupling. This can be corrected by coupling individually each degree of freedom to the heat bath, instead of coupling globally the system. We present the results of the investigation on three aspects of molecular dynamics simulations and their effect on the computed average area per lipid: (I) the accuracy in the computation of electrostatic interactions, the energy, and the virial, (II) long range Lennard-Jones interactions for systems with symmetry in one plane, and (III) thermodynamic baths. We show that the average area per lipid remains constant for simulations of systems containing 32, 64, and 256 lipids.     

Condensed-phase relaxation of multilevel quantum systems. I. An exactly solvable model

An analytically solvable model of multilevel condensed-phase quantum dynamics relevant to vibrational relaxation and electron transfer is presented. Exact solutions are derived for the reduced system density matrix dynamics of a degenerate N-level quantum system characterized by nearest-neighbor hopping and off-diagonal coupling (which is linear in the bath coordinates) to a harmonic oscillator bath. We demonstrate that for N> 2 the long-time steady-state system site occupation probabilities are not the same for all sites; that is, they are distributed in a non-Boltzmann manner, which depends on the initial conditions and the number of levels in the system. Although the system-bath Hamiltonian considered here is restricted in form, the availability of an exact solution enables us to study the model in all regions of an extensive parameter space.     

Excitation spectrum of heisenberg spin chains

In the framework of the Heisenberg model of interacting spins we prove the existence of a continuum of zero-activation excitations in the spectrum of alternant spin chains with a non-singular ground state. It is shown that for certain weak restrictions on the spin moment M, the excitations of the lower part of the spectrum with given M are continuous. The results are extended to the case of spin anisotropy, which effectively appears in the spin Hamiltonian method with the use of the Hartree-Fock approximation. For one-dimensional molecular systems described by the Parizer-Papp-Pople (PPP) Hamiltonian, the existence of a continuum of zero-activation excitations is established, provided that the ground state of the system is nondegenerate and nonsingular.Translated from Teoreticheskaya i Éksperimental' naya Khimiya, Vol. 24, No. 2, pp. 204–208, March–April, 1988.The author thanks A. A. Ovchinnikov and V. Ya. Krivnov for discussions of the results of the study.     

Benchmark calculations for dissipative dynamics of a system coupled to an anharmonic bath with the multiconfiguration time-dependent Hartree method

In this paper, we present benchmark results for dissipative dynamics of a harmonic oscillator coupled to an anharmonic bath of Morse oscillators. The microscopic Hamiltonian has been chosen so that the anharmonicity can be adjusted as a free parameter, and its effect can be isolated. This leads to a temperature dependent spectral density of the bath, which is studied for ohmic and lorentzian cases. Also, we compare numerically exact multiconfiguration time-dependent Hartree results with approximate solutions using continuous configuration time-dependent self-consistent field and local coherent state approximation.     

Continuum limit semiclassical initial value representation for dissipative systems

In this paper, we consider a dissipative system in which the system is coupled linearly to a harmonic bath. In the continuum limit, the bath is defined via a spectral density and the classical system dynamics is given in terms of a generalized Langevin equation. Using the path integral formulation and factorized initial conditions, it is well known that one can integrate out the harmonic bath, leaving only a path integral over the system degrees of freedom. However, the semiclassical initial value representation treatment of dissipative systems has usually been limited to a discretized treatment of the bath in terms of a finite number of bath oscillators. In this paper, the continuum limit of the semiclassical initial value representation is derived for dissipative systems. As in the path integral, the action is modified with an added nonlocal term, which expresses the influence of the bath on the dynamics. The first order correction term to the semiclassical initial value approximation is also derived in the continuum limit.     

A local coherent-state approximation to system-bath quantum dynamics

A novel quantum method to deal with typical system-bath dynamical problems is introduced. Subsystem discrete variable representation and bath coherent-state sets are used to write down a multiconfigurational expansion of the wave function of the whole system. With the help of the Dirac-Frenkel variational principle, simple equations of motion--a kind of Schrodinger-Langevin equation for the subsystem coupled to (pseudo) classical equations for the bath--are derived. True dissipative dynamics at all times is obtained by coupling the bath to a secondary, classical Ohmic bath, which is modeled by adding a friction coefficient in the derived pseudoclassical bath equations. The resulting equations are then solved for a number of model problems, ranging from tunneling to vibrational relaxation dynamics. Comparison of the results with those of exact, multiconfiguration time-dependent Hartree calculations in systems with up to 80 bath oscillators shows that the proposed method can be very accurate and might be of help in studying realistic problems with very large baths. To this end, its linear scaling behavior with respect to the number of bath degrees of freedom is shown in practice with model calculations using tens of thousands of bath oscillators.     

Multidimensional reactive scattering with quantum trajectories: dynamics with 50-200 vibrational modes

The dynamics of ensembles containing thousands of quantum trajectories are studied for multidimensional systems undergoing reactive scattering. The Hamiltonian and equations of motion are formulated in curvilinear reaction path coordinates, for the case of a planar (zero-torsion) reaction path. In order to enhance the computational efficiency, an improved least squares fitting procedure is introduced. This scheme involves contracted basis sets and the use of inner and outer stencils around points where fitting is performed. This method is applied to reactive systems with 50-200 harmonic vibrational modes which are coupled to motion along the reaction coordinate. Dynamical results, including trajectory evolution and time-dependent reaction probabilities, are presented and power law scaling of computation time with the number of vibrational modes is described.     

Theory of quantum energy transfer in spin chains: superexchange and ballistic motion

Quantum energy transfer in a chain of two-level (spin) units, connected at its ends to two thermal reservoirs, is analyzed in two limits: (i) in the off-resonance regime, when the characteristic subsystem excitation energy gaps are larger than the reservoirs frequencies, or the baths temperatures are low and (ii) in the resonance regime, when the chain excitation gaps match populated bath modes. In the latter case, the model is studied using a master equation approach, showing that the dynamics is ballistic for the particular chain model explored. In the former case, we analytically study the system dynamics utilizing the recently developed Energy-Transfer Born-Oppenheimer formalism [L.-A. Wu and D. Segal, Phys. Rev. E 83, 051114 (2011)], demonstrating that energy transfers across the chain in a superexchange (bridge assisted tunneling) mechanism, with the energy current decreasing exponentially with distance. This behavior is insensitive to the chain details. Since at low temperatures the excitation spectrum of molecular systems can be truncated to resemble a spin chain model, we argue that the superexchange behavior obtained here should be observed in widespread systems satisfying the off-resonance condition.     

Optimal control simulation of the Deutsch-Jozsa algorithm in a two-dimensional double well coupled to an environment

The quantum Deutsch-Jozsa algorithm is implemented by using vibrational modes of a two-dimensional double well. The laser fields realizing the different gates (NOT, CNOT, and HADAMARD) on the two-qubit space are computed by the multitarget optimal control theory. The stability of the performance index is checked by coupling the system to an environment. Firstly, the two-dimensional subspace is coupled to a small number Nb of oscillators in order to simulate intramolecular vibrational energy redistribution. The complete (2+Nb)D problem is solved by the coupled harmonic adiabatic channel method which allows including coupled modes up to Nb=5. Secondly, the computational subspace is coupled to a continuous bath of oscillators in order to simulate a confined environment expected to be favorable to achieve molecular computing, for instance, molecules confined in matrices or in a fullerene. The spectral density of the bath is approximated by an Ohmic law with a cutoff for some hundreds of cm(-1). The time scale of the bath dynamics (of the order of 10 fs) is then smaller than the relaxation time and the controlled dynamics (2 ps) so that Markovian dissipative dynamics is used.     

Symmetric double proton tunneling in formic acid dimer: a diabatic basis approach

A model of double proton tunneling in formic acid dimer is developed using a reaction surface Hamiltonian. The surface includes the symmetric OH stretch plus the in-plane stretch and bend interdimer vibrations. The surface Hamiltonian is coupled to a bath of five A1g and B3g normal modes obtained at the D2h transition state structure. Eigenstates are calculated using Davidson and block-Davidson iterative methods. Strong mode specific effects are found in the tunneling splittings for the reaction surface, where splittings are enhanced upon excitation of the interdimer bend motion. The results are interpreted within the framework of a diabatic representation of reaction surface modes. The splitting patterns observed for the reaction surface eigenstates are only slightly modified upon coupling to the bath states. Splitting patterns for the bath states are also determined. It is found that predicting these splittings is greatly complicated by subtle mixings with the inter-dimer bend states.     

Theoretical study of electronically excited radical cations of naphthalene and anthracene as archetypal models for astrophysical observations. Part I. Static aspects

Motivated by the recent discovery of new diffuse interstellar bands and results from laboratory experiments, ab initio quantum chemistry calculations are carried out for the lowest six electronic states of naphthalene and anthracene radical cations. The calculated adiabatic electronic energies are utilized to construct suitable diabatic electronic Hamiltonians in order to perform nuclear dynamics studies in Part II. Complex entanglement of the electronic states is established for both the radical cations and the coupling surfaces among them are also derived in accordance with the symmetry selection rules. Critical examination of the coupling parameters of the Hamiltonian suggests that 29 (out of 48) and 31 (out of 66) vibrational modes are relevant in the nuclear dynamics in the six lowest electronic states of naphthalene and anthracene radical cations, respectively.     

Spin dynamics in sodium clusters

We investigate the dynamics of the electron cloud in a metal cluster from the regime of small oscillations up to high excitations in the multi-plasmon regime. Particular attention is paid to the effect of spin modes. Test cases are a spin saturated cluster, Na₉ ⁺, and a strongly spin polarized cluster, the spherical isomer of Na₆ ⁺. The spectral distributions are dominated by the Mie plasmon resonance and a strong collective spin mode in all cases considered. Nonlinear effects as cross talk between multi-plasmon states or chaotic pattern cannot be observed in the investigated energy range.     

Simple model for transport phenomena: microscopic construction of Maxwell demonlike engine

We present a microscopic Hamiltonian framework to develop Maxwell demonlike engine. Our model consists of an equilibrium thermal bath and a nonequilibrium bath, latter generated by driving with an external stationary, Gaussian noise. The engine we develop can be considered as a device to extract work by modifying internal fluctuations. Our theoretical analysis focuses on finding the essential ingredients necessary for generating fluctuation induced transport under nonequilibrium condition. An important outcome of our model is that the net motion occurs when the nonlinear bath is modulated by the external noise, creating the nonzero effective temperature even when the temperature of both the baths are the same.