Transition rate using density matrices

So far, one can be much more successful in calculating a rate constant when one knows in advance that it exists, than in answering the question of whether it exists. A considerable breakthrough in this area was the solution of the spin-boson problem, which, however, has only limited relevance to any problem in chemistry because it neglects the effects of intrawell dynamics (vibrational relaxation) and does not describe thermally activated transitions. A number of attempts have been made to go beyond the two-level system approximation, but the basic question of how vibrational relaxation affects the transition from coherent oscillations to exponential decay awaits a quantitative solution. Such a solution might be obtained by numerical computation of real-time path integrals for the density matrix using the influence functional technique. 

Warshel and Chu [42] and Hwang et al. [60] were the first to calculate the contribution of tunneling and other nuclear quantum effects to PT in solution and enzyme catalysis, respectively. Since then, and in particular in the past few years, there has been a significant increase in simulations of quantum mechanical-nuclear effects in enzyme and in solution reactions [16]. The approaches used range from the quantized classical path (QCP) (for example. Refs. [4, 58, 95]), the centroid path integral approach [54, 55], and variational transition state theory [96], to the molecular dynamics with quantum transition (MDQT) surface hopping method [31] and density matrix evolution [97-99]. Most studies of enzymatic reactions did not yet examine the reference water reaction, and thus could only evaluate the quantum mechanical contribution to the enzyme rate constant, rather than the corresponding catalytic effect. However, studies that explored the actual catalytic contributions (for example. Refs. [4, 58, 95]) concluded that the quantum mechanical contributions are similar for the reaction in the enzyme and in solution, and thus, do not contribute to catalysis. 

This system in many cases can be simplified further. For example, if we have a broad spectral line excitation with a not very intense laser radiation, we have a situation for an open transition when 7 Ti, H. In practical cases this condition is often fulfilled at excitation with cw lasers operating in a multimode regime. If the homogeneous width of spectral transition usually is in the range of 10 MHz, then the laser radiation spectral width broader than 100 MHz usually can be considered as a broad line excitation. In this case we can use a procedure known as adiabatic elimination. It means that we are assuming that optical coherence pi2 decays much faster than the populations of the levels puJ = 1,2. Then we can find stationary solution for off-diagonal elements for the density matrix and afterwards find a rate equations for populations in this limit. For the two level system we will have... 

In order to evaluate the effects of the approach to the localization transition upon the reaction rate we first define the rate quantum mechanically using the dynamics of the density matrix (37). Following Zwanzig s method (38) one can partition the space into reactant states, indicated by a and product states, indicated by b (39,40). Then a projector... 

Chachaty and co-workers [8.20, 8.37, 8.38] were first to describe correlated internal motions in alkyl chains of surfactant molecules that form lyotropic liquid crystals. The last section described an extension of the master equation method of Wittebort and Szabo [8.4] to treat spin relaxation of deuterons on a chain undergoing trans-gauche jump rotations in liquid crystals. This method was also followed by Chachaty et al. to deal with spin relaxation of nuclei in surfactants. However, they assumed that the conformational changes occur by trans-gauche isomerization about one bond at a time. In their spectral density calculations (see Section 8.3.1), they used a transition rate matrix that was constructed from the jump rate Wi, W2, and Ws about each bond. Since W3 is much smaller than Wi and W2, the time scale of internal motions was practically governed by Wi and W2 of each C-C bond. Since... 

The high external magnetic field can mix the zero-field spin sublevels and therefore will change the rate constants of the T-S transitions [65]. This concept can be applied to triplet chromophores in biopolymers with slow rotational diffusion. Moreover, the Zeemann perturbation can induce additional mechanisms for the T-S mixing [65, 66]. The standard Liouville equation for spin density matrix [65] can be used for prediction of magnetic field effects on spin-selective photoprocesses and enzymatic reactions if the ZFS and hyperfine coupling parameters are properly interpreted in connection with the T-S transitions rate constants. 

The transition energies are not sharp but rather have finite widths due to thermal vibrations in the solid. Thus, the wave functions and the density of states can be treated as functions of energy. The wave functions can be normalized with respect to energy and a Dirac delta function used for the density of final states to insure conservation of energy. Then Einstein s A and B coefficients can be used to relate the transition matrix elements to experimentally measurable quantities such as oscillator strengths and luminescence lifetimes. For electric dipole-dipole interaction the energy transfer rate becomes... 

The spontaneous emission rate into the optical mode l, Wspontw, contains the dipole matrix element of the two electronic states involved in the transition [21]. Thus WspontM will not be changed by placing the optically active medium inside an optical cavity. However, the optical mode density, g(ve), is strongly modified by the cavity. Next, the changes in optical mode density will be used to calculate the changes in spontaneous emission rate. 

The IVR rates in the aniline(Ar)i case are rather sensitive to the densities of states. In this case, the values of pi (=N(E used in the construction of the IVR transition probability matrix are calculated with a direct count method, Vibrational frequencies of 45 cm l (stretch), and IS cm l (both bends) and anharmonicities of 3% (i.e., Avi, i+i /Avi-i,i =. 97) are assumed. 

It is possible to evaluate the rate of energy transfer using the Fermi golden rule [3] of quantum mechanics. At the lowest level of calculation, the matrix element for the process is just the electrostatic interaction between the transition densities of the donor and of the acceptor. At large distances, this may be approximated by the interaction between the corresponding transition dipoles, which is proportional to where R is the distance between the donor and the acceptor [3]. The rate of energy transfer is proportional to the square of it, and it is usual to write the rate as... 

---