It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix.
It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix.
For an open, time-dependent quantum system, Lindblad derived the most general modification of the quantum Liouville equation in the Markovian approximation that models environmental effects while preserving the non-negativity of the system's density matrix. While Lindblad's modification is correct for N-electron density matrices, solution of the Liouville equation with a Lindblad operator causes the one-electron reduced density matrix (1-RDM) to violate the Pauli exclusion principle. Consequently, after a short time, the 1-RDM is not representable by an ensemble N-electron density matrix (not ensemble N-representable). In this communication, we derive the necessary and sufficient constraints on the Lindbladian matrix within the Lindblad operator to ensure that the 1-RDM remains N-representable for all time. The theory is illustrated by considering the relaxation of an excitation in several molecules F2, N2, CO, and BeH2 subject to environmental noise.
The second-quantized, many-electron, atomic, and molecular Hamiltonian is partitioned both by the identity or labeling of the spin orbitals and by the dynamics of the spin orbitals into a system coupled to a bath. The electron bath is treated by a molecular time scale generalized Langevin equation approach designed to include one-electron dynamics in the system dynamics. The bath is formulated as an equivalent chain of spin orbitals through the introduction of equivalent-chain annihilation and creation operators. Both the dynamics and the quantum grand canonical statistical properties of the electron bath are examined. Two versions for the statistical properties of the bath are pursued. Using a weak bath assumption, a bath statistical average is defined which allows one to achieve a reduced dynamics description of the electron system which is coupled to the electron bath. In a strong bath assumption effective Hamiltonians are obtained which reproduce the dynamics of the bath and which lead to the same results as found in the weak bath assumption. The effective (but exact) Hamiltonian is found to be a one-electron Hamiltonian. A reduced dynamics equation of motion for the system population matrix is derived and found to agree with a previous version. This equation of motion is useful for studying electron transfer in the system when coupled to an electron bath.
We contrast Dirac's theory of transition probabilities and the theory of nonadiabatic transition probabilities, applied to a perturbed system that is coupled to a bath. In Dirac's analysis, the presence of an excited state |k0⟩ in the time-dependent wave function constitutes a transition. In the nonadiabatic theory, a transition occurs when the wave function develops a term that is not adiabatically connected to the initial state. Landau and Lifshitz separated Dirac's excited-state coefficients into a term that follows the adiabatic theorem of Born and Fock and a nonadiabatic term that represents excitation across an energy gap. If the system remains coherent, the two approaches are equivalent. However, differences between the two approaches arise when coupling to a bath causes dephasing, a situation that was not treated by Dirac. For two-level model systems in static electric fields, we add relaxation terms to the Liouville equation for the time derivative of the density matrix. We contrast the results obtained from the two theories. In the analysis based on Dirac's transition probabilities, the steady state of the system is not an equilibrium state; also, the steady-state population ρkk,s increases with increasing strength of the perturbation and its value depends on the dephasing time T2. In the nonadiabatic theory, the system evolves to the thermal equilibrium with the bath. The difference is not simply due to the choice of basis because the difference remains when the results are transformed to a common basis.
A general procedure is given for finding the equation of motion of the density matrix of a system in contact with a thermal bath, as for example a nuclear spin system weakly coupled to a crystal lattice. The thermal bath is treated both classically and quantum mechanically, and the theory is similar to, and a generalization of, conventional theories of time proportional transition probabilities. Relaxation of the system by the thermal bath is expressed by a linear matrix operator, and it is stressed that elements of this operator can be regarded as secular or nonsecular perturbations on the equation of motion and can be treated accordingly. When the motion of the system is slow compared to that of the thermal bath, the equation of motion can be expressed in an operator form which is independent of representation. If the system has a time-dependent Hamiltonian which varies slowly compared to the motion of the thermal bath, the same equation of motion is obeyed and the system is relaxed by the bath toward a Boltzmann distribution with respect to its instantaneous Hamiltonian. If the time variation of the Hamiltonian is more rapid, higher order corrections to the equation of motion must be applied. The theory is applied to spin-lattice relaxation of a coupled nuclear spin system in a metal, for arbitrary externally applied fixed magnetic field.
We examine the exact equation of motion for the relaxation of populations of strongly correlated electrons after a nonequilibrium excitation by a pulsed field, and prove that the populations do not change when the Green's functions have no average time dependence. We show how the average time dependence enters into the equation of motion to lowest order and describe what governs the relaxation process of the electron populations in the long‐time limit. While this result may appear, on the surface, to be required by any steady‐state solution, the proof is nontrivial, and provides new critical insight into how nonequilibrium populations relax, which goes beyond the assumption that they thermalize via a simple relaxation rate determined by the imaginary part of the self‐energy, or that they can be described by a quasi‐equilibrium condition with a Fermi‐Dirac distribution and a time‐dependent temperature. We also discuss the implications of this result to approximate theories, which may not satisfy the exact relation in the equation of motion.
In a previous work [Phys. Rev. E 48, 4263 (1993)] we have derived a nonlinear one-dimensional kinetic equation for the distribution function of particles obeying an exclusion principle. In the present work, on the same grounds, we extend this kinetics to D-dimensional continuous or discrete space, in order to study the distribution function of particles obeying a generalized exclusion-inclusion Pauli principle (EIP). This exclusion or inclusion principle is introduced into the classical transition rates by means of an inhibition or an enhancement factor, which contains a parameter \ensuremath{\kappa}, whose values range between -1 and +1 and can balance the effect of the full or partial validity of EIP. After deriving the kinetic equation we obtain a general expression of the stationary distribution function, depending on the value we give to the parameter \ensuremath{\kappa}. When we limit ourselves to Brownian particles, we derive exactly for \ensuremath{\kappa}=-1 the Fermi-Dirac (FD) distribution, for \ensuremath{\kappa}=0 the Maxwell-Boltzmann distribution, and for \ensuremath{\kappa}=1 the Bose-Einstein (BE) distribution. When \ensuremath{\kappa} assumes an intermediate value, except zero, between the extreme values -1 and +1, we obtain statistical distributions different from the FD and BE ones. We attribute to the parameter \ensuremath{\kappa} the meaning of the degree of indistinguishability of identical particles, the degree of antisymmetrization, or the symmetrization of the wave function of the particle system.