Decoherence Allows Model Reduction in Nonadiabatic Dynamics

Publication Date (Web): July 29, 2015 ... not only provide the physical mechanism for NAMD trajectory branding and improve the accuracy of NAMD simula...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCA

Decoherence Allows Model Reduction in Nonadiabatic Dynamics Simulations Dhara J. Trivedi† and Oleg V. Prezhdo*,‡ †

Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, United States Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States



ABSTRACT: A nonadiabatic (NA) molecular dynamics (MD) simulation requires calculation of NA coupling matrix elements, the number of which scales as a square of the number of basis states. The basis size can be huge in studies of nanoscale materials, and calculation of the NA couplings can present a significant bottleneck. A quantum-classical approximation, NAMD overestimates coherence in the quantum, electronic subsystem, requiring decoherence correction. Generally, decoherence times decrease with increasing energy separation between pairs of states forming coherent superpositions. Since rapid decoherence stops quantum dynamics, one expects that decoherence-corrected NAMD can eliminate the need for calculation of NA couplings between energetically distant states, notably reducing the computational cost. Considering several types of dynamics in a semiconductor quantum dot, we demonstrate that indeed, decoherence allows one to reduce the number of needed NA coupling matrix elements. If the energy levels are spaced closer than 0.1 eV, one obtains good results while including only three nearest-neighbor couplings, and in some cases even with just the first nearest-neighbor coupling scheme. If the energy levels are spaced by about 0.4 eV, the nearest-neighbor model fails, while three or more nearest-neighbor schemes also provide good results. In comparison, the results of NAMD simulation without decoherence vary continuously with changes in the number of NA couplings. Thus, decoherence effects induced by coupling to a quantum-mechanical environment not only provide the physical mechanism for NAMD trajectory branding and improve the accuracy of NAMD simulations, but also afford significant computational savings. attention.26,27 The simulations of photoinduced dynamics in large molecular systems involving coupled electronic excited states require an “on-the-fly”28,29 calculation of energies, forces and nonadiabatic couplings (NACs). While the number of energies and forces is a linear function of the number of the electronic states involved in the dynamics, the number of the NACs scales quadratically. As the number of states increases, the NAC calculation start to pose a computational bottleneck. For example, multiple exciton generation even in small quantum dots (QDs) involves hundreds of thousands of double exciton states.30,31 Therefore, it is highly desirable to decrease the number of NAC matrix elements needed for accurate description of NAMD. Alternative schemes have been proposed in order to reduce the computational cost associated with evaluating the NAC terms.32 Although such schemes introduce error in the propagation of the electronic wave function, they help to reduce the computational cost. A recent study has shown that the standard SH technique requires evaluation of all or most NAC terms.27 It has long been recognized and appreciated33 that MQCD overestimates the magnitude of coherence effects in the

1. INTRODUCTION Mixed quantum-classical dynamics (MQCD) methodologies provide an intuitive approach to time-domain simulations of physical and chemical processes in condensed phases.1−9 Only a small subset of particles, such as electrons, holes and protons, are dealt with quantum mechanically in MQCD. The remaining parts of the system, e.g., heavy atoms, are treated classically. MQCD simulations of the nonadiabatic molecular dynamics (NAMD) provide an atomistic time-domain description of complex ultrafast photophysical and photochemical processes, including charge and energy transfer, and nonradiative relaxation.10−18 They generate a comprehensive perspective on the elastic and inelastic electron−phonon scattering dynamics of photoexcited charge carriers in nanoscale materials, directly mimicking numerous time-resolved experiments, and offering insights into the fundamental mechanisms underlying key processes in optics,19 photovoltaics,20,21 electronics,22,23 and related fields. A number of programs are capable of performing NAMD, employing the Ehrenfest24 and surface hopping (SH)25 methods. The stability of these methods with respect to changes in parameters and the number of trajectories has been extensively analyzed for simple models, while the robustness of the method applied to larger systems has received less © XXXX American Chemical Society

Received: June 19, 2015

A

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A quantum subsystem. An equivalent subsystem coupled to a quantum-mechanical rather than classical environment loses coherence quite rapidly, typically on a femtosecond time scale.34−36 Decoherence effects within a subsystem can be captured properly by considering quantum dynamics of open systems, which extends the unitary evolution of the Schrödinger equation for the wave function into nonunitary stochastic Schrödinger equation,37 or, alternatively, which generalizes the quantum Liouville−von Neumann equation for the density matrix to a quantum Fokker−Planck equation.38 A classical environment in an MQCD simulation does not modify the Schrödinger or von Neumann equation, but rather provides a classical external field, and therefore, cannot properly describe loss of coherence. It has been shown that decoherence effects have a strong influence on the calculated time scales, as exemplified by the quantum Zeno effect in the limit of infinitely fast decoherence.39−41 Further, decoherence provides a physical mechanism for trajectory branching in MQCD, leading naturally to SH algorithms.37,42−47 Inclusion of decoherence strengthens quantum-classical correlations, in the context of trajectory surface hopping simulations, and can increase in the internal stability.48,49 Generally, decoherence effects become more important as the energy gaps between pairs of states forming coherent superpositions grow. Since decoherence slows down quantum dynamics in general, rapid decoherence between energetically distant states shuts off populations transfer, suggesting that NAC between such states can be ignored. Combined with the fact that the magnitude of NAC typically decreases with increasing state separation,50 decoherence effects can be used to reduce the number of NAC evaluations and to increase the computational efficiency of NAMD. In this paper, we investigate the possibility of reducing the number of NAC terms in a NAMD simulation that incorporated decoherence effects. We compare performance of several NAC reduction schemes within fewest switches surface hopping (FSSH), which is the most popular NAMD approach, and within decoherence induced surface hopping (DISH), which uses decoherence as a SH mechanism. Focusing on a CdSe QD, we consider three different models, including a dense manifold of quantum states, a sparse state manifold, and a combination of dense and sparse states. Confirming earlier results, we find that SH without decoherence requires full consideration of NAC matrix elements. At the same time, we show that indeed, decoherence effects allow one to construct reduced models of NAMD. In the case of a dense state manifold, it is already sufficient to include three nearestneighbor NACs. Even the nearest neighbor scheme can be quite accurate for some properties. As the energy gaps between quantum states grow, one requires a more complete NAC model. Overall, methods such as DISH, which include decoherence effects, can scale better with number of states than the standard approaches, such as FSSH. The demonstration that decoherence effects can be used to reduce the number of NAC matrix elements provides an important tool for improving the efficiency of NAMD simulations of large systems.

provides a rigorous description of quantum dynamics. Here, Ψ(r, R, t) is the full wave function, and Η̂(r, R, t) is the total Hamiltonian of the system. They depend on the electronic (r) and nuclear (R) coordinates. A quantum-classical approximation greatly reduces computational requirements on the solution of the TDSE by treating nuclei as classical particles. This leads to elimination of the terms in the Hamiltonian operator associated with the nuclear degrees of freedom and introduction of parametric dependence of the electronic Hamiltonian on classical nuclear coordinates. In the Born− Oppenheimer approximation, one can write the total wave function as a product of the electronic, ψ(r; R(t)), and nuclear, χ(R(t)), components, Ψ(r , R , t ) = ψ (r ; R(t ))χ (R(t ))

Nuclei evolve according to classical equations of motion in NAMD, and the classical paths are determined by a quantum force, which can be computed as an adiabatic Hellmann− Feynman force or mean-field Ehrenfest force. The electronic wave functions are obtained as solutions of the electronic TDSE, iℏ

∂Ψ(r , R , t ) = Η̂(r , R , t )Ψ(r , R , t ) ∂t

∂ψ (r ; R(t )) = Η̂el(r ; R(t ))ψ (r ; R(t )) ∂t

(3)

where Η̂el is the electronic Hamiltonian of the system. It can be further decomposed into 2

ℏ Η̂el = − 2

1 2 ∇i + Vel(r ; R(t )) mi

∑ i

(4)

Here, Vel(r; R(t)) includes the electron−nuclear, electron− electron, and nuclear−nuclear interactions. It is evaluated for given nuclear positions, neglecting nuclear quantum-effects. The first term in eq 4 is the electron kinetic energy. The nuclear kinetic energy is ignored, as justified by the small electron-to-nucleus mass ratio. The electronic wave function can be expanded on the basis of adiabatic wave functions, ψ (r ; R(t )) =

∑ ci(t )φi(r ; R(t ))

(5)

i

obtained as eigenstates of the electronic Hamiltonian. A wide variety of methods have been developed to calculate the electronic energies, densities, and wave functions, while balancing the desired accuracy with computational cost. The less expensive methods range from semiempirical to various formulations of density functional theory (DFT).51−58 We employ ab initio DFT and its time-dependent generalization.59 The adiabatic Kohn-sham (KS) orbitals, obtained from timeindependent DFT are used as the basis in the time-dependent DFT (TDDFT) simulation.57 Substituting a basis expansion into the TDSE results in a set of coupled equations of the expansion coefficients: iℏ

2. THEORETICAL BACKGROUND The time-dependent Schrödinger equation (TDSE), ıℏ

(2)

∂ci(t ) = ∂t

∑ cj(t )(ϵjδij − iℏd ij·Ṙ ) j

(6)

Here, ϵj is the energy of the adiabatic KS orbital, and dij Ṙ is the NAC. The NAC is calculated numerically10 as the overlap of the adiabatic KS orbitals φi and φj at sequential MD time steps Δt:

(1) B

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A d ij·Ṙ = ⟨φi|∇R |φj⟩· ≈

dR = dt

φi

DISH, as compared to decoherence corrected FSSH, is that branching appears directly due to decoherence rather than with the help of an ad hoc algorithm, which has to be corrected for decoherence effects. The original idea of decoherence-driven trajectory branching was implemented within the stochastic mean-field approach.37 DISH extends the TDSE through a stochastic adaption of the time-evolution operator of quantum system coupled to a quantum environment, and provides a mechanism for surface hopping that models diverging evolutions.42 The calculations were performed with the Vienna Ab initio Simulation Package (VASP).75 The PBE functional76 with projector-augmented-wave (PAW) pseudopotentials77 were employed in a converged plane wave basis. The simulations were performed in a periodically replicated cubic cell with at least 8 A of vacuum between the QD replicas. Following geometry optimization at 0K, we heated the QD to 300 K with repeated velocity rescaling. Then, we generated a 3 ps microcanonical MD trajectory using the Verlet algorithm with the 1 fs time step and Hellman-Feynman forces. At each time step, we computed NACs. Along with KS orbital energies, this time-dependent information was stored and used to perform NAMD within the recently developed PYthon eXtension for Ab Initio Dynamics (PYXAID) program.72,78

∂ φ ∂t j

1 (⟨φi(t )|φj(t + Δt )⟩ − ⟨φi(t + Δt )|φj(t )⟩) 2Δt

(7)

The NAC describes electron−phonon interactions that arise from dependence of the adiabatic KS orbitals on the nuclei coordinates, R(t). The above equations establish evolution of the electrons subject to the classical field of nuclei. Defining the evolution of the nuclei constitutes the so-called quantum backreaction problem, which prescribes how the quantum subsystem influences the classical subsystem.60−67 NAMD provides a generalization of ordinary molecular dynamics to include transitions between electronic states. The branching of the overall nuclear wave function into wavepackets, correlated with different electronic states, is a quantum mechanical phenomenon. Correlations between the nuclear motion and electronic states are built into MQCD using SH techniques.25,42,68,69 SH can be interpreted as a master equation in which transition rates are nonperturbative and evolve in time. FSSH minimizes the number of surface hops by prescribing a probability for hopping between the electronic states that is based on the change in the electronic state populations. The probability of hopping from state i to j during the time step Δt is given in FSSH by,25 ⎧ bji ·Δt ⎫ , 0⎬ gij(t , Δt ) = max⎨ ⎩ aii(t ) ⎭

3. RESULTS AND DISCUSSION The study is performed with the Cd33Se33 “magic” size cluster. Known experimentally,79 the cluster is 1.3 nm in diameter and is one of the smallest stable CdSe QDs that support a crystalline-like core.80,81 These properties make Cd33Se33 an excellent model for quantum-mechanical and quantumdynamical studies of semiconductor QDs. While smaller stable clusters with bulk topology are known, for instance Pb16Se16,82−84 the properties of CdSe are more representative of semiconductor QDs. For instance, CdSe has a significantly higher density of valence band (VB) than conduction band (CB) states, while the band structure of PbSe is unusually symmetrical. The density of states (DOS) of Cd33Se33 is shown in Figure 1. The optimized structure of the system is given as insert. Surface reconstruction results in saturation of unsaturated valencies and removal of associated localized surface states. Consequently, this opens the energy gap between the highest occupied (HOMO) and lowest unoccupied (LUMO) Kohn−

(8) −1

where aji(t) = cj*(t)ci(t) and bji = 2ℏ Im(aji⟨φj|H|φi⟩) − 2Re(ajidji). A hop from state i to state j can occur only if the electronic occupation of state i decreases and the electronic occupation of state j increases, minimizing the number of hops. By contrast, the original SH technique70 uses the standard quantum mechanical probabilities to hop, and allows transitions even if state populations do not change. The total electron−nuclear energy is conserved by velocity rescaling,10,25 normally performed in the direction of the NAC vector, dji. If the velocity component is too small in this direction to reflect the increase in the electronic energy, the hop is rejected. This phenomenon leads to the detailed balance between upward and downward transitions.71 In order to further simplify the implementation of FSSH, it is assumed that the energy exchanged between the electronic and vibrational degrees of freedom during a hop is rapidly redistributed among all of the vibrational modes. The velocity rescaling and hop rejection procedure is replaced by multiplying the probability for transitions upward in energy by the Boltzmann factor.13,72 This allows for the use of a ground state nuclear trajectory to determine the time-dependent potential that drives the electronic dynamics. This simplification of the original FSSH technique provides great computational savings, eliminating the need to run multiple excited state trajectories. By treating nuclei classically, the original FSSH scheme excludes coherence loss that occurs if the electronic subsystem is coupled to quantum nuclei.73,74 Decoherence can be viewed as an environment induced destruction of quantum state superpositions. The effect of decoherence on the evolution of quantum subsystem is implemented here using the DISH approach.42 Originally, Rossky and co-workers introduced the decoherence correction to FSSH63 by working with the Liouville-von Neumann equation. A conceptual advantage of

Figure 1. Density of states in Cd33Se33. The optimized geometry is shown in the inset. C

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

partial coupling (PC) schemes. PC1 includes NAC only between nearest neighbor states. PC3 and PC5 include three and five nearest-neighbor NACs, respectively. The CC, PC1, PC3 and PC5 NAC schemes are applied to the full intraband relaxation, dense manifold, and sparse manifold models (Figure 2) to obtain the rates of decay of the initial electronic state population and electronic energy (Figures 3−5). Figure 6 reports values of NACs, while Figure 7

Sham orbitals of the CdSe cluster. The calculated zero temperature energy gap is about 1.6 eV. This value is underestimateda typical feature of most semilocal DFT functionals, such as PBE used in our calculations. The problem is attributed to incomplete elimination of the self-interaction energy by these functionals. Nonlocal corrections are routinely used as remedy to the problem, for instance, in hybrid functionals that include a portion of the exact Hartree−Fock exchange.81 However, hybrid functionals require significant additional computational efforts. The VB DOS is higher than the CB DOS in agreement with the effective-mass theory,85 which uses a higher effective mass for holes than for electrons. The LUMO is well separated energetically from the rest of the CB manifold, giving rise to the so-called phonon bottleneck to the electron relaxation, observed experimentally,86 and rationalized by ab initio time-domain simulation.41 We focus on intraband relaxation of the excited electron in the CdSe QD and consider three scenarios, Figure 2. First, we

Figure 3. Decay of (a) population and (b) energy of the initially excited state for complete intraband relaxation at 300 K, obtained with FSSH and DISH using complete coupling (CC), coupling with the first five nearest neighbors (PC5), coupling with the first three nearest neighbors (PC3), and coupling with nearest neighbors (PC1). Figure 2. Schematics of (a) complete intranband relaxation, (b) dense manifold relaxation, and (c) sparse manifold relaxation via nonadiabatic vibronic dynamics. Schemes (a) and (b) employ 15 excited states, from LUMO to LUMO+14, spanning about 1.3 eV of energy. Scheme (c) has 7 excited states distributed evenly in energy with gaps around 0.4 eV, spanning 2.3 eV of energy. The LUMO to LUMO+1 gap is close to 0.4 eV (Figure 1).

study the relaxation involving all CB states. Second, we restrict the dynamics to the dense state manifold, excluding the LUMO, which is separated energetically from the rest of the CB, Figure 1. Finally, we construct a sparse manifold model by retaining few CB orbitals such that they are separated by approximately 0.4 eV, which is the gap between LUMO and LUMO+1. The dense manifold model provides a good representation of excitation dynamics in nanoscale materials, such as QDs. The sparse manifold model is typical of molecular systems, which exhibit large energy gaps between electronic states. We report two sets of data. On one hand, we present timedependent population of the initially excited state. Such data correspond to pump−probe experiments that trace the signal at a fixed photoexcitation energy, for instance, by monitoring signal bleach and recovery. On the other hand, we show the evolution of the excited electron energy, which decreases due to coupling to vibrations. Energy losses are detrimental to many applications, in which nonradiative relaxation competes with productive processes, such as charge separation. The ratio of the rates of energy losses and productive processes determines the efficiency of a device. In order to investigate the possibility of reducing the number of NAC matrix elements needed to obtain reliable results, we compare the complete coupling (CC) description with three

Figure 4. Same as Figure 3, for the dense state manifold defined in Figure 2.

characterizes pure-dephasing times. Tables 1 and 2 present the time scales of decay of the initial state population and energy, respectively. The time scales are obtained by exponential fits of the data shown in Figures 3−5. The results are obtained with both the decoherence free FSSH approach,25 and the DISH method42 that uses decoherence as the physical mechanism of SH. Note that including decoherence in DISH slows down the dynamics relative to decoherence free FSSH, as should be expected.35,36,78 As demonstrated by the data, DISH allows significant reduction in the number of needed NACs, while the D

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Figure 5. Same as Figure 3, for the sparse state manifold defined in Figure 2.

Figure 7. Same as Figure 6, for the pure-dephasing (decoherence) time.

Table 1. Time Scales (ps) of Relaxation of Population of the Initially Excited Statea model

intraband relaxation

dense manifold

sparse manifold

FSSH_CC FSSH_PC5 FSSH_PC3 FSSH_PC1 DISH_CC DISH_PC5 DISH_PC3 DISH_PC1

0.33 0.40 0.46 0.79 0.97 1.03 1.08 1.20

0.33 0.40 0.45 0.80 0.98 1.04 1.09 1.19

2.58 2.65 2.91 9.80 26.0 27.2 28.8 55.6

a

The data are presented for NAMD with (DISH) and without (FSSH) decoherence. CC stands for complete coupling. PCN refers to partial coupling schemes, with N denoting the number of nearest neighbor couplings. Intraband, dense and sparse manifold relaxations are illustrated in Figure 2.

Table 2. Time Scales (ps) of Relaxation of the Excess Electronic Energya

Figure 6. Canonically averaged root-mean-squared value of the NAC for each state involved in complete intraband relaxation for different coupling schemes: complete coupling (CC), coupling with first five nearest neighbors (PC5), coupling with first three nearest neighbors (PC3), and coupling with nearest neighbor (PC1). The NAC increases with orbital energy and is largest between nearest neighbor states.

FSSH results are much more sensitive to reduction of the number of NAC matrix elements. Focusing on the numerical values presented in Tables 1 and 2, we observe that FSSH shows much more variation with the number of NAC matrix elements than DISH. The variation is larger for the energy relaxation (Table 2) than for the decay of the initial state population (Table 1). For instance, the FSSH decay of the initial state population in the dense manifold model changes from 0.33 ps for the CC scheme to 0.80 ps for PC1, almost by a factor of 3. The corresponding variation is insignificant in DISH, 0.98 and 1.19 ps, respectively. With both FSSH and DISH, the initial state population in the sparse manifold model show little variance for the CC, PC5, and PC3 schemes, while PC1 gives notably different results. The relaxation of energy in the dense manifold model differs between the CC and PC1 schemes by an order of magnitude in FSSH, 1.28 ps vs 12.7 ps. The corresponding difference in

a

model

intraband relaxation

dense manifold

sparse manifold

FSSH_CC FSSH_PC5 FSSH_PC3 FSSH_PC1 DISH_CC DISH_PC5 DISH_PC3 DISH_PC1

0.76 2.32 4.61 12.2 10.5 18.7 22.3 30.5

1.28 2.56 4.71 12.7 11.9 20.0 24.3 32.4

4.94 5.77 7.00 54.7 88.5 83.2 95.6 234

The notation is the same as in Table 1.

DISH is only a factor of 3: 11.9 ps versus 32.4 ps. Similarly to the initial state population decay, the energy relaxation times agree between the CC, PC5, and PC3 schemes applied to the sparse manifold model, while the PC1 scheme shows significant deviations, with both FSSH and DISH. Just as with the dense manifold, here the difference between CC and PC1 is a factor of 10 in FSSH, and a factor of 3 in DISH. The reported results lead to the overall conclusion that reduction in the number of NAC matrix elements is much more justified in DISH than FSSH. It is dangerous to apply the nearest neighbor scheme PC1 with sparse states manifold, since E

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

decoherence time can be understood qualitatively. The NAC depends on overlap of wave functions at different timesteps. Such overlap generally decreases for wave functions with a larger difference in the number of nodes. The pure-dephasing (decoherence) time depends on fluctuation in the energy gap between a given pair of states. More distant states have larger energy gaps, and typically, larger gap fluctuations. Because quantum transition rates are more sensitive to the coherence time than the NAC, methods including decoherence effects, such as DISH, are more suitable for reduction of NAC to a few nearest neighbors.

one may miss an important NAC matrix element. Including several nearest neighbors gives good results, especially with DISH. The energy relaxation time scale is much more sensitive to elimination of distant NAC matrix elements than decay of the initial state population. For practical purposes, one can employ the PC3 scheme with DISH and expect to obtain results that differ from the CC scheme by a factor of 2 or less. The possibility of model reduction in NAMD arises due to the general properties of NAC matrix elements and decoherence times. Figure 6 shows root-mean-squared values of the NAC between a given state and all other states in a particular scheme, averaged over the MD trajectory. Similarly, Figure 7 shows pure-dephasing (decoherence) times between a given state and all other states involved in the four coupling schemes. NAC grows with energy. More importantly, NACs between nearest neighbors are significantly larger than NACs between distant states. This fact forms the basis for model reduction within FSSH. The behavior of pure-dephasing times gives additional motivation of model reduction within DISH. It is important to note that the frequency of quantum transitions depends on the decoherence time much more strongly than on the NAC. As follows from Fermi’s golden rule,87,88 the transition rate is a quadratic function of NAC and a Gaussian function of the decoherence time. Just as NAC, the puredephasing time is larger for nearest neighbor states than for distant states. Longer coherence favors quantum transitions.50 As a result, the nearest neighbor pathways dominate the dynamics much more strongly in DISH than FSSH. Qualitatively, once can argue that NACs are larger for nearest neighbor states because these states have more “similar” wave functions, resulting in larger overlaps in the NAC expression (eq 7). For instance, distant wave functions have very different nodal structure. A shift in the nodes induced by a nuclear displacement will give both positive and negative contributions, which will cancel each other. Decoherence time decreases with increasing fluctuation of the energy gap between a given pair of states.34 It is natural to expect that larger energy gaps between more distant states also have larger fluctuations.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Financial support of the U.S. National Science Foundation, Grant CHE-1300118, is gratefully acknowledged.



REFERENCES

(1) Herman, M. F. Dynamics by Semiclassical Methods. Annu. Rev. Phys. Chem. 1994, 45, 83−111. (2) Schatz, G. C. Scattering Theory and Dynamics: Time-Dependent and Time-Independent Methods. J. Phys. Chem. 1996, 100, 12839− 12847. (3) Borgis, D.; Hynes, J. T. Curve Crossing Formulation for Proton Transfer Reactions in Solution. J. Phys. Chem. 1996, 100, 1118−1128. (4) Ben-Nun, M.; Quenneville, J.; Martinez, T. J. Ab Initio Multiple Spawning: Photochemistry from First Principles Quantum Molecular Dynamics. J. Phys. Chem. A 2000, 104, 5161−5175. (5) Yarkony, D. R. Conical Intersections: The New Conventional Wisdom. J. Phys. Chem. A 2001, 105, 6277−6293. (6) Miller, W. H. The Semiclassical Initial Value Representation: A Potentially Practical Way for Adding Quantum Effects to Classical Molecular Dynamics Simulations. J. Phys. Chem. A 2001, 105, 2942− 2955. (7) Kapral, R. Progress in the Theory of Mixed Quantum-Classical Dynamics. Annu. Rev. Phys. Chem. 2006, 57, 129−157. (8) Hammes-Schiffer, S.; Soudackov, A. V. Proton-Coupled Electron Transfer in Solution, Proteins, and Electrochemistry. J. Phys. Chem. B 2008, 112, 14108−14123. (9) Gorshkov, V. N.; Tretiak, S.; Mozyrsky, D. Semiclassical MonteCarlo Approach for Modelling Non-Adiabatic Dynamics in Extended Molecules. Nat. Commun. 2013, 4, 2144. (10) Hammes-Schiffer, S.; Tully, J. C. Proton Transfer in Solution: Molecular Dynamics with Quantum Transitions. J. Chem. Phys. 1994, 101, 4657−4667. (11) Ben-Nun, M.; Quenneville, J.; Martínez, T. J. Ab Initio Multiple Spawning: Photochemistry from First Principles Quantum Molecular Dynamics. J. Phys. Chem. A 2000, 104, 5161−5175. (12) Craig, C. F.; Duncan, W. R.; Prezhdo, O. V. Trajectory Surface Hopping in the Time-Dependent Kohn-Sham Approach for ElectronNuclear Dynamics. Phys. Rev. Lett. 2005, 95, 163001. (13) Prezhdo, O. V.; Duncan, W. R.; Prezhdo, V. V. Photoinduced Electron Dynamics at the Chromophore−Semiconductor Interface: A Time-Domain Ab Initio Perspective. Prog. Surf. Sci. 2009, 84, 30−68. (14) Bedard-Hearn, M. J.; Sterpone, F.; Rossky, P. J. Nonadiabatic Simulations of Exciton Dissociation in Poly-P-Phenylenevinylene Oligomers. J. Phys. Chem. A 2010, 114, 7661−7670. (15) Akimov, A. V.; Neukirch, A. J.; Prezhdo, O. V. Theoretical Insights into Photoinduced Charge Transfer and Catalysis at Oxide Interfaces. Chem. Rev. 2013, 113, 4496−4565.

4. CONCLUSIONS To recapitulate, we have demonstrated that decoherence effects, which improve accuracy of MQCD and provide the physical basis for trajectory branching, also help to reduce the computational effort. In particular, a major computational bottleneck in NAMD of large systems, computation of NAC matrix elements, which scales as the square of the number of states, can be reduced to a few nearest neighbors, provided that a NAMD includes decoherence. Denser state manifolds allow one to consider fewer nearest neighbors. This is a particularly favorable feature, since dynamics involving denser state manifolds require more basis states. It is not advisable to restrict the calculation to just the first nearest neighbor, since one may easily miss an important dynamical pathway. Investigation of electron−phonon energy relaxation over a broad energy range requires a more complete NAC scheme than study of decay of the initial state population. Both NAC and decoherence time exhibit the properties needed for reduction of the computational effort. They both decrease for pairs of more distant states. Larger NAC and longer coherence times accelerate quantum dynamics in general. Therefore, decrease in the NAC and coherence time with state separation reduces the importance of quantum transitions between distant states. This feature of the NAC and F

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A (16) Neukirch, A. J.; Hyeon-Deuk, K.; Prezhdo, O. V. Time-Domain Ab Initio Modeling of Excitation Dynamics in Quantum Dots. Coord. Chem. Rev. 2014, 263−264, 161−181. (17) Trivedi, D. J.; Wang, L.; Prezhdo, O. V. Auger-Mediated Electron Relaxation Is Robust to Deep Hole Traps: Time-Domain Ab Initio Study of Cdse Quantum Dots. Nano Lett. 2015, 15, 2086−2091. (18) Tafen, D. N.; Prezhdo, O. V. Size and Temperature Dependence of Electron Transfer between Cdse Quantum Dots and a Tio2 Nanobelt. J. Phys. Chem. C 2015, 119, 5639−5647. (19) Lal, S.; Hafner, J. H.; Halas, N. J.; Link, S.; Nordlander, P. Noble Metal Nanowires: From Plasmon Waveguides to Passive and Active Devices. Acc. Chem. Res. 2012, 45, 1887−1895. (20) Galoppini, E. Linkers for Anchoring Sensitizers to Semiconductor Nanoparticles. Coord. Chem. Rev. 2004, 248, 1283−1297. (21) Kamat, P. V. Boosting the Efficiency of Quantum Dot Sensitized Solar Cells through Modulation of Interfacial Charge Transfer. Acc. Chem. Res. 2012, 45, 1906−1915. (22) Kambhampati, P. Multiexcitons in Semiconductor Nanocrystals: A Platform for Optoelectronics at High Carrier Concentration. J. Phys. Chem. Lett. 2012, 3, 1182−1190. (23) Kim, D. K.; Lai, Y.; Diroll, B. T.; Murray, C. B.; Kagan, C. R. Flexible and Low-Voltage Integrated Circuits Constructed from HighPerformance Nanocrystal Transistors. Nat. Commun. 2012, 3, 1216. (24) Li, X.; Tully, J. C.; Schlegel, H. B.; Frisch, M. J. Ab Initio Ehrenfest Dynamics. J. Chem. Phys. 2005, 123, 084106. (25) Tully, J. C. Molecular Dynamics with Electronic Transitions. J. Chem. Phys. 1990, 93, 1061−1071. (26) Fabiano, E.; Keal, T. W.; Thiel, W. Implementation of Surface Hopping Molecular Dynamics Using Semiempirical Methods. Chem. Phys. 2008, 349, 334−347. (27) Nelson, T.; Fernandez-Alberti, S.; Chernyak, V.; Roitberg, A. E.; Tretiak, S. Nonadiabatic Excited-State Molecular Dynamics: Numerical Tests of Convergence and Parameters. J. Chem. Phys. 2012, 136, 054108. (28) Lasorne, B.; Robb, M. A.; Worth, G. A. Direct Quantum Dynamics Using Variational Multi-Configuration Gaussian Wavepackets. Implementation Details and Test Case. Phys. Chem. Chem. Phys. 2007, 9, 3210−3227. (29) Worth, G. A.; Robb, M. A.; Lasorne, B. Solving the TimeDependent Schrödinger Equation for Nuclear Motion in One Step: Direct Dynamics of Non-Adiabatic Systems. Mol. Phys. 2008, 106, 2077−2091. (30) Hyeon-Deuk, K.; Prezhdo, O. V. Multiple Exciton Generation and Recombination Dynamics in Small Si and Cdse Quantum Dots: An Ab Initio Time-Domain Study. ACS Nano 2012, 6, 1239−1250. (31) Hyeon-Deuk, K.; Prezhdo, O. V. Photoexcited Electron and Hole Dynamics in Semiconductor Quantum Dots: Phonon-Induced Relaxation, Dephasing, Multiple Exciton Generation and Recombination. J. Phys.: Condens. Matter 2012, 24, 363201. (32) Pittner, J.; Lischka, H.; Barbatti, M. Optimization of Mixed Quantum-Classical Dynamics: Time-Derivative Coupling Terms and Selected Couplings. Chem. Phys. 2009, 356, 147−152. (33) Hänggi, P.; Talkner, P.; Borkovec, M. Reaction-Rate Theory: Fifty Years after Kramers. Rev. Mod. Phys. 1990, 62, 251−341. (34) Akimov, A. V.; Prezhdo, O. V. Persistent Electronic Coherence Despite Rapid Loss of Electron−Nuclear Correlation. J. Phys. Chem. Lett. 2013, 4, 3857−3864. (35) Jaeger, H. M.; Green, J. R.; Prezhdo, O. V. Decoherence Reduces Thermal Energy Loss in Graphene Quantum Dots. Appl. Phys. Lett. 2013, 103, 073111. (36) Liu, J.; Neukirch, A. J.; Prezhdo, O. V. Phonon-Induced PureDephasing of Luminescence, Multiple Exciton Generation, and Fission in Silicon Clusters. J. Chem. Phys. 2013, 139, 164303. (37) Prezhdo, O. V. Mean Field Approximation for the Stochastic Schrödinger Equation. J. Chem. Phys. 1999, 111, 8366−8377. (38) Bittner, E. R.; Rossky, P. J. Quantum Decoherence in Mixed Quantum-Classical Systems: Nonadiabatic Processes. J. Chem. Phys. 1995, 103, 8130−8143.

(39) Onuchic, J. N.; Wolynes, P. G. Classical and Quantum Pictures of Reaction Dynamics in Condensed Matter: Resonances, Dephasing, and All That. J. Phys. Chem. 1988, 92, 6495−6503. (40) Skourtis, S. S.; Balabin, I. A.; Kawatsu, T.; Beratan, D. N. Protein Dynamics and Electron Transfer: Electronic Decoherence and NonCondon Effects. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 3552−3557. (41) Kilina, S. V.; Neukirch, A. J.; Habenicht, B. F.; Kilin, D. S.; Prezhdo, O. V. Quantum Zeno Effect Rationalizes the Phonon Bottleneck in Semiconductor Quantum Dots. Phys. Rev. Lett. 2013, 110, 180404. (42) Jaeger, H. M.; Fischer, S.; Prezhdo, O. V. Decoherence-Induced Surface Hopping. J. Chem. Phys. 2012, 137, 22A545. (43) Schwartz, B. J.; Bittner, E. R.; Prezhdo, O. V.; Rossky, P. J. Quantum Decoherence and the Isotope Effect in Condensed Phase Nonadiabatic Molecular Dynamics Simulations. J. Chem. Phys. 1996, 104, 5942−5955. (44) Hack, M. D.; Truhlar, D. G. A Natural Decay of Mixing Algorithm for Non-Born−Oppenheimer Trajectories. J. Chem. Phys. 2001, 114, 9305−9314. (45) Zhu, C.; Jasper, A. W.; Truhlar, D. G. Non-Born−Oppenheimer Trajectories with Self-Consistent Decay of Mixing. J. Chem. Phys. 2004, 120, 5543−5557. (46) Wong, K. F.; Rossky, P. J. Solvent-Induced Electronic Decoherence: Configuration Dependent Dissipative Evolution for Solvated Electron Systems. J. Chem. Phys. 2002, 116, 8429−8438. (47) Bedard-Hearn, M. J.; Larsen, R. E.; Schwartz, B. J. Mean-Field Dynamics with Stochastic Decoherence (Mf-Sd): A New Algorithm for Nonadiabatic Mixed Quantum/Classical Molecular-Dynamics Simulations with Nuclear-Induced Decoherence. J. Chem. Phys. 2005, 123, 234106. (48) Fang, J.-Y.; Hammes-Schiffer, S. Improvement of the Internal Consistency in Trajectory Surface Hopping. J. Phys. Chem. A 1999, 103, 9399−9407. (49) Granucci, G.; Persico, M. Critical Appraisal of the Fewest Switches Algorithm for Surface Hopping. J. Chem. Phys. 2007, 126, 134114. (50) Kilina, S. V.; Kilin, D. S.; Prezhdo, V. V.; Prezhdo, O. V. Theoretical Study of Electron−Phonon Relaxation in Pbse and Cdse Quantum Dots: Evidence for Phonon Memory. J. Phys. Chem. C 2011, 115, 21641−21651. (51) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A: At., Mol., Opt. Phys. 1988, 38, 3098−3100. (52) Becke, A. D. Density-Functional Thermochemistry. Iii. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (53) Burke, K. Perspective on Density Functional Theory. J. Chem. Phys. 2012, 136, 150901. (54) Burke, K.; Perdew, J. P.; Ernzerhof, M. Why the Generalized Gradient Approximation Works and How to Go Beyond It. Int. J. Quantum Chem. 1997, 61, 287−293. (55) Elliott, P.; Lee, D.; Cangi, A.; Burke, K. Semiclassical Origins of Density Functionals. Phys. Rev. Lett. 2008, 100, 256406. (56) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864−B871. (57) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133− A1138. (58) Scuseria, G. E. Comparison of Coupled-Cluster Results with a Hybrid of Hartree−Fock and Density Functional Theory. J. Chem. Phys. 1992, 97, 7528−7530. (59) Marques, M. A. L.; Gross, E. K. U. Time-Dependent Density Functional Theory. Annu. Rev. Phys. Chem. 2004, 55, 427−455. (60) Aleksandrov, I. The Statistical Dynamics of a System Consisting of a Classical and a Quantum Subsystem. Z. Naturforsch., A: Phys. Sci. 1981, 36, 902. (61) Anderson, A. Quantum Backreaction on ″Classical″ Variables. Phys. Rev. Lett. 1995, 74, 621−625. (62) Prezhdo, O. V.; Kisil, V. V. Mixing Quantum and Classical Mechanics. Phys. Rev. A: At., Mol., Opt. Phys. 1997, 56, 162−175. G

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A (63) Prezhdo, O. V.; Rossky, P. J. Mean-Field Molecular Dynamics with Surface Hopping. J. Chem. Phys. 1997, 107, 825−834. (64) Prezhdo, O. V.; Rossky, P. J. Evaluation of Quantum Transition Rates from Quantum-Classical Molecular Dynamics Simulations. J. Chem. Phys. 1997, 107, 5863−5878. (65) Kapral, R.; Ciccotti, G. Mixed Quantum-Classical Dynamics. J. Chem. Phys. 1999, 110, 8919−8929. (66) Prezhdo, O. V. A Quantum-Classical Bracket That Satisfies the Jacobi Identity. J. Chem. Phys. 2006, 124, 201104. (67) Prezhdo, O. V.; Brooksby, C. Quantum Backreaction through the Bohmian Particle. Phys. Rev. Lett. 2001, 86, 3215−3219. (68) Wang, L.; Prezhdo, O. V. A Simple Solution to the Trivial Crossing Problem in Surface Hopping. J. Phys. Chem. Lett. 2014, 5, 713−719. (69) Wang, L.; Trivedi, D.; Prezhdo, O. V. Global Flux Surface Hopping Approach for Mixed Quantum-Classical Dynamics. J. Chem. Theory Comput. 2014, 10, 3598−3605. (70) Tully, J. C.; Preston, R. K. Trajectory Surface Hopping Approach to Nonadiabatic Molecular Collisions: The Reaction of H+ with D2. J. Chem. Phys. 1971, 55, 562−572. (71) Parandekar, P. V.; Tully, J. C. Mixed Quantum-Classical Equilibrium. J. Chem. Phys. 2005, 122, 094102. (72) Akimov, A. V.; Prezhdo, O. V. The Pyxaid Program for NonAdiabatic Molecular Dynamics in Condensed Matter Systems. J. Chem. Theory Comput. 2013, 9, 4959−4972. (73) Zurek, W. H. Decoherence, Einselection, and the Quantum Origins of the Classical. Rev. Mod. Phys. 2003, 75, 715−775. (74) Leggett, A. J.; Chakravarty, S.; Dorsey, A. T.; Fisher, M. P. A.; Garg, A.; Zwerger, W. Dynamics of the Dissipative Two-State System. Rev. Mod. Phys. 1987, 59, 1−85. (75) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (76) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (77) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (78) Akimov, A. V.; Prezhdo, O. V. Advanced Capabilities of the Pyxaid Program: Integration Schemes, Decoherence Effects, Multiexcitonic States, and Field-Matter Interaction. J. Chem. Theory Comput. 2014, 10, 789−804. (79) Kasuya, A.; et al. Ultra-Stable Nanoparticles of Cdse Revealed from Mass Spectrometry. Nat. Mater. 2004, 3, 99−102. (80) Puzder, A.; Williamson, A. J.; Gygi, F.; Galli, G. Self-Healing of Cdse Nanocrystals: First-Principles Calculations. Phys. Rev. Lett. 2004, 92, 217401. (81) Kilina, S.; Ivanov, S.; Tretiak, S. Effect of Surface Ligands on Optical and Electronic Spectra of Semiconductor Nanoclusters. J. Am. Chem. Soc. 2009, 131, 7717−7726. (82) Long, R.; Prezhdo, O. V. Ab Initio Nonadiabatic Molecular Dynamics of the Ultrafast Electron Injection from a Pbse Quantum Dot into the Tio2 Surface. J. Am. Chem. Soc. 2011, 133, 19240−19249. (83) Bao, H.; Habenicht, B. F.; Prezhdo, O. V.; Ruan, X. Temperature Dependence of Hot-Carrier Relaxation in Pbse Nanocrystals: An Ab Initio Study. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 79, 235306. (84) Kilina, S. V.; Craig, C. F.; Kilin, D. S.; Prezhdo, O. V. Ab Initio Time-Domain Study of Phonon-Assisted Relaxation of Charge Carriers in a Pbse Quantum Dot. J. Phys. Chem. C 2007, 111, 4871−4878. (85) Nozik, A. J. Spectroscopy and Hot Electron Relaxation Dynamics in Semiconductor Quantum Wells and Quantum Dots. Annu. Rev. Phys. Chem. 2001, 52, 193−231. (86) Pandey, A.; Guyot-Sionnest, P. Slow Electron Cooling in Colloidal Quantum Dots. Science 2008, 322, 929−932. (87) Hyeon-Deuk, K.; Madrid, A. B.; Prezhdo, O. V. Symmetric Band Structures and Asymmetric Ultrafast Electron and Hole Relaxations in

Silicon and Germanium Quantum Dots: Time-Domain Ab Initio Simulation. Dalton Transactions 2009, 10069−10077. (88) Chen, L.; Bao, H.; Tan, T.; Prezhdo, O. V.; Ruan, X. Shape and Temperature Dependence of Hot Carrier Relaxation Dynamics in Spherical and Elongated Cdse Quantum Dots. J. Phys. Chem. C 2011, 115, 11400−11406.

H

DOI: 10.1021/acs.jpca.5b05869 J. Phys. Chem. A XXXX, XXX, XXX−XXX