Generalized trajectory surface hopping method based on the Zhu-Nakamura theory

We present a generalized formulation of the trajectory surface hopping method applicable to a general multidimensional system. The method is based on the Zhu-Nakamura theory of a nonadiabatic transition and therefore includes the treatment of classically forbidden hops. The method uses a generalized recipe for the conservation of angular momentum after forbidden hops and an approximation for determining a nonadiabatic transition direction which is crucial when the coupling vector is unavailable. This method also eliminates the need for a rigorous location of the seam surface, thereby ensuring its applicability to a wide class of chemical systems. In a test calculation, we implement the method for the DH(2) (+) system, and it shows a remarkable agreement with the previous results of C. Zhu, H. Kamisaka, and H. Nakamura, [J. Chem. Phys. 116, 3234 (2002)]. We then apply it to a diatomic-in-molecule model system with a conical intersection, and the results compare well with exact quantum calculations. The successful application to the conical intersection system confirms the possibility of directly extending the present method to an arbitrary potential of general topology.     

Simulation of the photodynamics of azobenzene on its first excited state: comparison of full multiple spawning and surface hopping treatments

We have studied the cis-->trans and trans-->cis photoisomerization of azobenzene after n-->pi* excitation using the full multiple spawning (FMS) method for nonadiabatic wave-packet dynamics with potential-energy surfaces and couplings determined "on the fly" from a reparametrized multiconfigurational semiempirical method. We compare the FMS results with a previous direct dynamics treatment using the same potential-energy surfaces and couplings, but with the nonadiabatic dynamics modeled using a semiclassical surface hopping (SH) method. We concentrate on the dynamical effects that determine the photoisomerization quantum yields, namely, the rate of radiationless electronic relaxation and the character of motion along the reaction coordinate. The quantal and semiclassical results are in good general agreement, confirming our previous analysis of the photodynamics. The SH method slightly overestimates the rate of excited state decay, leading in this case to lower quantum yields.     

Excited state dynamics of Michler's ketone: a laser flash photolysis study

Steady state absorption and fluorescence as well as the time resolved absorption studies in the pico and subpicosecond time domain have been performed to characterize the excited singlet and triplet states of Michler's ketone (MK). The nature of the lowest excited singlet (S₁) and triplet (T₁) states depends on the polarity of the solvent - in nonpolar solvents they have either pure nπ * character or mixed character of nπ * and ππ * states but in more polar solvents the states have CT character. Concentration dependence of the shapes of the fluorescence as well the excited singlet and triplet absorption spectra provide the evidence for the association of the MK molecules in the ground state.     

Regarding the validity of the time-dependent Kohn-Sham approach for electron-nuclear dynamics via trajectory surface hopping

The implementation of fewest-switches surface-hopping (FSSH) within time-dependent Kohn-Sham (TDKS) theory [Phys. Rev. Lett. 95, 163001 (2005)] has allowed us to study successfully excited state dynamics involving many electronic states in a variety of molecular and nanoscale systems, including chromophore-semiconductor interfaces, semiconductor and metallic quantum dots, carbon nanotubes and graphene nanoribbons, etc. At the same time, a concern has been raised that the KS orbital basis used in the calculation provides only approximate potential energy surfaces [J. Chem. Phys. 125, 014110 (2006)]. While this approximation does exist in our method, we show here that FSSH-TDKS is a viable option for computationally efficient calculations in large systems with straightforward excited state dynamics. We demonstrate that the potential energy surfaces and nonadiabatic transition probabilities obtained within the TDKS and linear response (LR) time-dependent density functional theories (TDDFT) agree semiquantitatively for three different systems, including an organic chromophore ligating a transition metal, a quantum dot, and a small molecule. Further, in the latter case the FSSH-TDKS procedure generates results that are in line with FSSH implemented within LR-TDDFT. The FSSH-TDKS approach is successful for several reasons. First, single-particle KS excitations often give a good representation of LR excitations. In this regard, DFT compares favorably with the Hartree-Fock theory, for which LR excitations are typically combinations of multiple single-particle excitations. Second, the majority of the FSSH-TDKS applications have been performed with large systems involving simple excitations types. Excitation of a single electron in such systems creates a relatively small perturbation to the total electron density summed over all electrons, and it has a small effect on the nuclear dynamics compared, for instance, with thermal nuclear fluctuations. In such cases an additional, classical-path approximation can be made. Third, typical observables measured in time-resolved experiments involve averaging over many initial conditions. Such averaging tends to cancel out random errors that may be encountered in individual simulated trajectories. Finally, if the flow of energy between electronic and nuclear subsystems is insignificant, the ad hoc FSSH procedure is not required, and a straightforward mean-field, Ehrenfest approach is sufficient. Then, the KS representation provides rigorously a convenient and efficient basis for numerically solving the TDDFT equations of motion.     

Three-state trajectory surface hopping studies of the photodissociation dynamics of formaldehyde on ab initio potential energy surfaces

Full-dimensional, three-state, surface hopping calculations of the photodissociation dynamics of formaldehyde are reported on ab initio potential energy surfaces (PESs) for electronic states S(1), T(1), and S(0). This is the first such study initiated on S(1) with ab initio-calculated spin-orbit couplings among the three states. We employ previous PESs for S(0) and T(1), and a new PES for S(1), which we describe here, as well as new spin-orbit couplings. The time-dependent electronic state populations and the branching ratio of radical products produced from S(0) and T(1) states and that of total radical products and molecular products at three total energies are calculated. Details of the surface hopping dynamics are described, and a novel pathway for isomerization on T(1) via S(0) is reported. Final translational energy distributions of H + HCO products from S(0) and T(1) are also reported as well as the translational energy distribution and final rovibrational distributions of H(2) products from the molecular channel. The present results are compared to previous trajectory calculations initiated from the global minimum of S(0). The roaming pathway leading to low rotational distribution of CO and high vibrational population of H(2) is observed in the present calculations.     

Theoretical study of the isomerization mechanism of azobenzene and disubstituted azobenzene derivatives

A series of azobenzenes was studied using ab initio methods to determine the substituent effects on the isomerization pathways. Energy barriers were determined from three-dimensional potential energy surfaces of the ground and electronically excited states. In the ground state (S(0)), the inversion pathway was found to be preferred. Our results show that electron donating substituents increase the isomerization barrier along the inversion pathway, whereas electron withdrawing substituents decrease it. The inversion pathway of the first excited state (S(1)) showed trans --> cis barriers with no curve crossing between S(0) and S(1). In contrast, a conical intersection was found between the ground and first excited states along the rotation pathway for each of the azobenzenes studied. No barriers were found in this pathway, and we therefore postulate that after n --> pi (S(1) <-- S(0)) excitation, the rotation mechanism dominates. Upon pi --> pi (S(2) <-- S(0)) excitation, there may be sufficient energy to open an additional pathway (concerted-inversion) as proposed by Diau. Our potential energy surface explains the experimentally observed difference in trans-to-cis quantum yields between S(1) and S(2) excitations. The concerted inversion channel is not available to the remaining azobenzenes, and so they must employ the rotation pathway for both n --> pi and pi --> pi excitations.     

Nonradiative deexcitation dynamics of 9H-adenine: an OM2 surface hopping study

The nonradiative relaxation of 9H-adenine was studied at the semiempirical OM2/MR-CI level using the surface-hopping approach. Geometry optimizations of energy minima and conical intersections as well as single-point calculations of excitation energies at critical points were performed to characterize the relevant potential energy surfaces of 9H-adenine and to assess the accuracy of the OM2 results. Surface-hopping calculations were performed to describe the nonradiative dynamics of 9H-adenine after vertical excitation into the optically active state. They showed that the deexcitation process is mainly governed by a two-step relaxation consisting of an ultrashort component and a longer component. These findings compare well with experimental results from time-resolved photoelectron spectroscopy.     

Carbene formation in its lower singlet state from photoexcited 3H-diazirine or diazomethane. A combined CASPT2 and ab initio direct dynamics trajectory study

The potential energy surfaces of the ground and valence excited states of both 3H-diazirine and diazomethane have been studied computationally by mean of the CASSCF method in conjunction with the cc-pVTZ basis set. The energies of the critical points found on such surfaces have been recomputed at the CASPT2/cc-pVTZ level. Additionally, ab initio direct dynamic trajectory calculations have been carried out on the S(1) and S(2) surfaces, starting each trajectory run at the region dominated by the conformational molecular rearrangement of diazomethane. It is found that both isomers are interconnected along a C(s)() reaction coordinate on each potential surface. Radiationless deactivation of the corresponding S(1) state of each isomer occurs through the same point on the surface, an S(1)/S(0) conical intersection. Thereafter, the system has enough energy to surmount the barrier which leads to dissociation products (CH(2) + N(2)) on S(0) state. Therefore, photoexcitation to S(1) state of either diazirine of diazomethane produces methylene in its lower singlet state on a very short time scale (ca. 100 fs). Furthermore, both isomers can generate excited singlet carbene when they are excited onto the S(2) surface; in this case, they lose the activation energy passing through another common S(2)/S(1) conical intersection and then proceed to dissociation into carbene and N(2) on the S(1) surface. For the special case of methylene, it rapidly experiences deexcitation to S(0) state.     

Simultaneous-trajectory surface hopping: a parameter-free algorithm for implementing decoherence in nonadiabatic dynamics

In this paper, we introduce a trajectory-based nonadiabatic dynamics algorithm which aims to correct the well-known overcoherence problem in Tully's popular fewest-switches surface hopping algorithm. Our simultaneous-trajectory surface hopping algorithm propagates a separate classical trajectory on each energetically accessible adiabatic surface. The divergence of these trajectories generates decoherence, which collapses the particle wavefunction onto a single adiabatic state. Decoherence is implemented without the need for any parameters, either empirical or adjustable. We apply our algorithm to several model problems and find a significant improvement over the traditional algorithm.     

Nonadiabatic molecular dynamics study of the cis-trans photoisomerization of azobenzene excited to the S1 state

Ab initio nonadiabatic dynamics simulations of cis-to-trans isomerization of azobenzene upon S(1) (n-π*) excitation are carried out employing the fewest-switches surface hopping method. Azobenzene photoisomerization occurs purely as a rotational motion of the central CNNC moiety. Two nonequivalent rotational pathways corresponding to clockwise or counterclockwise rotation are available. The course of the rotational motion is strongly dependent on the initial conditions. The internal conversion occurs via an S(0)/S(1) crossing seam located near the midpoint of both of these rotational pathways. Based on statistical analysis, it is shown that the occurrence of one or other pathway can be completely controlled by selecting adequate initial conditions.     

Excited states of porphyrin isomers and porphycene derivatives: a SAC-CI study

Excited states of free-base porphyrin isomers, porphycene (Pc), corrphycene (Cor), and hemiporphycene (hPc), were studied by the Symmetry-Adapted Cluster (SAC)/SAC-Configuration Interaction (CI) method. The absorption peaks of the porphyrin isomers were assigned on the basis of the SAC-CI spectra. The X, Y, X', and Y' bands of the porphyrin isomers, which have weak intensities, are identified. The differences in the Q-band absorptions among the isomers were clearly explained by the four-orbital model. In Cor and hPc, the wave function of the B-band corresponds to the mixture of the four-orbital excitations and the optically forbidden excitation of free-base porphin (P), due to the molecular symmetry lowering in the isomers. The B-band character is described by the five-orbital model in Pc and the six-orbital model in Cor and hPc. Two tetrazaporphycenes and two ring-extended (dibenzo) porphycenes were designed, and the Q-band transition moment was successfully controlled. These examples show that the control of the four-orbital energy levels is the guiding principle for pigment design in porphyrin compounds.     

The three-dimensional nonadiabatic dynamics calculation of DH(2)(+) and HD(2)(+) systems by using the trajectory surface hopping method based on the Zhu-Nakamura theory

A theoretical investigation on the nonadiabatic processes of the full three-dimensional D(+)+H(2) and H(+)+D(2) reaction systems has been performed by using trajectory surface hopping (TSH) method based on the Zhu-Nakamura (ZN) theory. This ZN-TSH method refers to not only classically allowed hops but also classically forbidden hops. The potential energy surface constructed by Kamisaka et al. is employed in the calculation. A new iterative method is proposed to yield the two-dimensional seam surface from the topography of the adiabatic potential surfaces, in which the inconvenience of directly solving the first-order partial differential equation is avoided. The cross sections of these two systems are calculated for three competing channels of the reactive charge transfer, the nonreactive charge transfer, and the reactive noncharge transfer, for ground rovibrational state of H(2) or D(2). Also, this study provides reaction probabilities of these three processes for the total angular momentum J=0 and ground initial vibrational state of H(2) or D(2). The calculated results from ZN-TSH method are in good agreement with the exact quantum calculations and the experimental measurements.     

Quasiclassical trajectory calculations of acetaldehyde dissociation on a global potential energy surface indicate significant non-transition state dynamics

A recent experimental study [Houston, P. L.; Kable, S. H. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 16079] of the photodissociation of acetaldehyde (CH(3)CHO) has suggested two distinct mechanisms for the production of the molecular products CH(4) + CO. One corresponds to the traditional transition state mechanism and the other to a transition state-skirting path similar to the roaming channel previously reported in formaldehyde. To investigate this theoretically, a full-dimensional potential energy surface (PES) has been constructed. The PES was fit with permutationally invariant polynomials to 135,000 points calculated using coupled cluster theory with single and double excitations and a perturbative treatment of triple excitations [CCSD(T)] and correlation consistent basis sets of double- and triple-zeta quality. To test the accuracy of the PES additional CCSD(T) and multireference configuration interaction calculations were carried out. Quasiclassical trajectory calculations were run on the PES starting at the acetaldehyde equilibrium geometry and also at the conventional transition state (TS) for the molecular products CH(4) + CO. The former calculations agree well with the experimental results of Houston and Kable; however, those from the TS do not. The implications for a non-transition state, roaming mechanism in this molecule are discussed.     

A direct CI method with a multiconfigurational reference state

A method for direct configuration–interaction (CI ) calculations with a multiconfigurational reference function is described. The reference state can contain several closed-shell electronic configurations and the CI expansion comprises all single and double replacements out of all these configurations. The resulting secular problem is solved using a variation-perturbation method. A number of examples are given showing the efficiency of the method. The largest CI expansion used in calculations with this program so far contains 76,471 spin- and space-symmetrized configurations.     

Diffusion dynamics of the Li+ ion on a model surface of amorphous carbon: a direct molecular orbital dynamics study

Diffusion processes of the Li+ ion on a model surface of amorphous carbon (Li+C96H24 system) have been investigated by means of the direct molecular orbital (MO) dynamics method at the semiempirical AM1 level. The total energy and energy gradient on the full-dimensional AM1 potential energy surface were calculated at each time step in the dynamics calculation. The optimized structure, where Li+ is located in the center of the cluster, was used as the initial structure at time zero. The dynamics calculation was carried out in the temperature range 100-1000 K. The calculations showed that the Li+ ion vibrates around the equilibrium point below 200 K, while the Li+ ion moves on the surface above 250 K. At intermediate temperatures (300 K < T < 400 K), the ion moves on the surface and falls in the edge regions of the cluster. At higher temperatures (600 K < T), the Li+ ion transfers freely on the surface and edge regions. The diffusion pathway of the Li+ ion was discussed on the basis of theoretical results.     

Classical trajectory study of collisions of Ar with alkanethiolate self-assembled monolayers: potential-energy surface effects on dynamics

We have investigated collisions between Ar and alkanethiolate self-assembled monolayers (SAMs) using classical trajectory calculations with several potential-energy surfaces. The legitimacy of the potential-energy surfaces is established through comparison with molecular-beam data and ab initio calculations. Potential-energy surfaces used in previous work overestimate the binding of Ar to the SAM, leading to larger energy transfer than found in the experiments. New calculations, based on empirical force fields that better reproduce ab initio calculations, exhibit improved agreement with the experiments. In particular, polar-angle-dependent average energies calculated with explicit-atom potential-energy surfaces are in excellent agreement with the experiments. Polar- and azimuthal-angle-dependent product translational energies are examined to gain deeper insight into the dynamics of Ar+SAM collisions.