A time-dependent quantum dynamics investigation of the guanine-cytosine system: a six-dimensional model

The dynamics of the guanine-cytosine base pair has been studied in the time-dependent quantum approach. A six-dimensional model involving the nonlinear three hydrogen bridges has been utilized. The modifications induced in the hydrogen transfer from a base to the other by the explicit inclusion of the out-of-plane hydrogen atom position in the three bridges have been evidenced and the consequences on stacking interaction and base pair opening are considered. The relevance of these aspects in biological properties has been suggested.     

Yb valence states in YbC2: a HERFD-XANES spectroscopic investigation

The valence state of Yb in YbC(2) was analyzed using high-energy-resolution fluorescence detection (HERFD) X-ray absorption near-edge structure (XANES) spectroscopy and time-of-flight neutron powder diffraction to clarify a controversy in the literature. The unit cell volume of YbC(2) suggests a mixed Yb valence, which was formerly determined to be 2.8 by magnetization measurements and paramagnetic neutron scattering techniques. However, the nature of the intermediate valence was not clearly established. Both homogeneous and heterogeneous mixed valences were assumed in different publications. The temperature-dependent behavior of the valence state was only predicted, albeit not explicitly studied. In this work, the valence state of Yb in YbC(2) is, therefore, investigated thoroughly by HERFD-XANES spectroscopy at low and high temperatures. Our measurements result in an average Yb valence of 2.81 that is temperature-independent from 15 to 1123 K. These findings are confirmed by neutron powder diffraction experiments, which reveal a constant C-C distance of 128.7(9) pm in a temperature range from 5 to 100 K. A significant temperature dependence of the Yb valence state in YbC(2) can, therefore, be excluded by our experimental results.     

Photoinduced ethane formation from reaction of ethene with matrix-isolated Ti, V, or Nb atoms

The reactions of matrix-isolated Ti, V, or Nb atoms with ethene (C(2)H(4)) have been studied by FTIR absorption spectroscopy. Under conditions where the ethene dimer forms, metal atoms react with the ethene dimer to yield matrix-isolated ethane (C(2)H(6)) and methane. Under lower ethene concentration conditions ( approximately 1:70 ethene/Ar), hydridic intermediates of the types HMC(2)H(3) and H(2)MC(2)H(2) are also observed, and the relative yield of hydrocarbons is diminished. Reactions of these metals with perdeuterioethene, and equimolar mixtures of C(2)H(4) and C(2)D(4), yield products that are consistent with the production of ethane via a metal atom reaction involving at least two C(2)H(4) molecules. The absence of any other observed products suggests the mechanism also involves production of small, highly symmetric species such as molecular hydrogen and metal carbides. Evidence is presented suggesting that ethane production from the ethene dimer is a general photochemical process for the reaction of excited-state transition-metal atoms with ethene at high concentrations of ethene.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

Exciton dissociation at donor-acceptor polymer heterojunctions: quantum nonadiabatic dynamics and effective-mode analysis

The quantum-dynamical mechanism of photoinduced subpicosecond exciton dissociation and the concomitant formation of a charge-separated state at a semiconducting polymer heterojunction is elucidated. The analysis is based upon a two-state vibronic coupling Hamiltonian including an explicit 24-mode representation of a phonon bath comprising high-frequency (C==C stretch) and low-frequency (torsional) modes. The initial relaxation behavior is characterized by coherent oscillations, along with the decay through an extended nonadiabatic coupling region. This region is located in the vicinity of a conical intersection hypersurface. A central ingredient of the analysis is a novel effective mode representation, which highlights the role of the low-frequency modes in the nonadiabatic dynamics. Quantum dynamical simulations were carried out using the multiconfiguration time-dependent Hartree method.     

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.     

A systematic theoretical investigation of the valence excited states of the diatomic molecules B2, C2, N2 and O2

A quantitative survey on the performance of multireference (MR), configuration interaction with all singles and doubles (CISD), MRCISD with the Davidson correction and MR-average quadratic coupled cluster (AQCC) methods for a wide range of excited states of the diatomic molecules B₂, C₂, N₂ and O₂ is presented. The spectroscopic constants r _e, ω_e, T _e and D _e for a total of 60 states have been evaluated and critically compared with available experimental data. Basis set extrapolations and size-extensivity corrections are essential for highly accurate results: MR-AQCC mean-errors of 0.001 ?, 10 cm⁻¹, 300 cm⁻¹ and 300 cm⁻¹ have been obtained for r _e, ω_e, T _e and D _e, respectively. Owing to the very systematic behavior of the results depending on the basis set and the choice of method, shortcomings of the calculations, such as Rydberg state coupling or insufficient configuration spaces, can be identified independently of experimental data. On the other hand, significant discrepancies with experiment for states which indicate no shortcomings whatsoever in the theoretical treatment suggest the re-evaluation of experimental results. The broad variety of states included in our survey and the uniform quality of the results indicate that the observed systematics is a general feature of the methods and, hence, is molecule-independent. Received: 12 June 2000 / Accepted: 1 September 2000 / Published online: 21 December 2000     

Fingerprints of delocalized transition states in quantum dynamics

Reactions with delocalized transition states (plateau reactions) can be characterized statically by their energy profile along the reaction path, where they exhibit a broad, flat region instead of one or several well-defined saddle points on the potential energy surface. Employing our new, highly flexible quantum dynamics code to perform two-dimensional and effective four-dimensional quantum wave packet propagations on ab initio based model potentials, we show that plateau reactions can also be discerned from the other standard reaction types by their dynamics.     

LAND-map, a linearized approach to nonadiabatic dynamics using the mapping formalism

We present a new approach for calculating quantum time correlation functions for systems whose dynamics exhibits relevant nonadiabatic effects. The method involves partial linearization of the full quantum path-integral expression for the time correlation function written in the nonadiabatic mapping Hamiltonian formalism. Our analysis gives an algorithm which is both numerically efficient and accurate as we demonstrate in test calculations on the spin-boson model where we find results in good agreement with exact calculations. The accuracy of our new approach is comparable to that of calculations performed using other approximate methods over a relatively broad range of model parameters. However, our method converges relatively quickly when compared with most alternative schemes. These findings are very encouraging in view of the application of the new method for studying realistic nonadiabatic model problems in the condensed phase.     

Transition states for the dimerization of 1,3-cyclohexadiene: a DFT, CASPT2, and CBS-QB3 quantum mechanical investigation

Quantum mechanical calculations using restricted and unrestricted B3LYP density functional theory, CASPT2, and CBS-QB3 methods for the dimerization of 1,3-cyclohexadiene (1) reveal several highly competitive concerted and stepwise reaction pathways leading to [4 + 2] and [2 + 2] cycloadducts, as well as a novel [6 + 4] ene product. The transition state for endo-[4 + 2] cycloaddition (endo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 28.7 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 19.0 kcal/mol) is not bis-pericyclic, leading to nondegenerate primary and secondary orbital interactions. However, the C(s) symmetric second-order saddle point on the B3LYP energy surface is only 0.3 kcal/mol above endo-2TS. The activation enthalpy for the concerted exo-[4 + 2] cycloaddition (exo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 30.1 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 21.1 kcal/mol) is 1.4 kcal/mol higher than that of the endo transition state. Stepwise pathways involving diallyl radicals are formed via two different C-C forming transition states (rac-5TS and meso-5TS) and are predicted to be competitive with the concerted cycloaddition. Transition states were located for cyclization from intermediate rac-5 leading to the endo-[4 + 2] (endo-2) and exo-[2 + 2] (anti-3) cycloadducts. Only the endo-[2 + 2] (syn-3) transition state was located for cyclization of intermediate meso-5. The novel [6 + 4] "concerted" ene transition state (threo-4TS, DeltaH(double dagger)(UB3LYP(0K)) = 28.3 kcal/mol) is found to be unstable with respect to an unrestricted calculation. This diradicaloid transition state closely resembles the cyclohexadiallyl radical rather than the linked cyclohexadienyl radical. Several [3,3] sigmatropic rearrangement transition states were also located and have activation enthalpies between 27 and 31 kcal/mol.     

Identification of unavoided crossings in nonadiabatic photoexcited dynamics involving multiple electronic states in polyatomic conjugated molecules

Radiationless transitions between electronic excited states in polyatomic molecules take place through unavoided crossings of the potential energy surfaces with substantial non-adiabatic coupling between the respective adiabatic states. While the extent in time of these couplings are large enough, these transitions can be reasonably well simulated through quantum transitions using trajectory surface hopping-like methods. In addition, complex molecular systems may have multiple "trivial" unavoided crossings between noninteracting states. In these cases, the non-adiabatic couplings are described as sharp peaks strongly localized in time. Therefore, their modeling is commonly subjected to the identification of regions close to the particular instantaneous nuclear configurations for which the energy surfaces actually cross each other. Here, we present a novel procedure to identify and treat these regions of unavoided crossings between non-interacting states using the so-called Min-Cost algorithm. The method differentiates between unavoided crossings between interacting states (simulated by quantum hops), and trivial unavoided crossings between non-interacting states (detected by tracking the states in time with Min-Cost procedure). We discuss its implementation within our recently developed non-adiabatic excited state molecular dynamics framework. Fragments of two- and four-ring linear polyphenylene ethynylene chromophore units at various separations have been used as a representative molecular system to test the algorithm. Our results enable us to distinguish and analyze the main features of these different types of radiationless transitions the molecular system undertakes during internal conversion.     

The investigation of nonadiabatic effects for the N + ND → N2 + D reaction

Nonadiabatic quantum dynamical calculations have been carried out on the two coupled potential energy surfaces (1²A′ and 2²A′) (Mota et al., J Theor Comput Chem 2009, 8, 849) for the title reaction. Initial state‐resolved reaction probabilities and cross sections for ground and excited states for collision energies of 0.005–1.0 eV are determined, respectively. Nonadiabatic transition is enhanced about four times by isotopic substitution of N + NH by N + ND reaction. It turns out that the nonadiabatic effects exert no significant contribution in the N + ND → N₂ + D reaction. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Exploring the role of decoherence in condensed-phase nonadiabatic dynamics: a comparison of different mixed quantum/classical simulation algorithms for the excited hydrated electron

Mixed quantum/classical (MQC) molecular dynamics simulation has become the method of choice for simulating the dynamics of quantum mechanical objects that interact with condensed-phase systems. There are many MQC algorithms available, however, and in cases where nonadiabatic coupling is important, different algorithms may lead to different results. Thus, it has been difficult to reach definitive conclusions about relaxation dynamics using nonadiabatic MQC methods because one is never certain whether any given algorithm includes enough of the necessary physics. In this paper, we explore the physics underlying different nonadiabatic MQC algorithms by comparing and contrasting the excited-state relaxation dynamics of the prototypical condensed-phase MQC system, the hydrated electron, calculated using different algorithms, including: fewest-switches surface hopping, stationary-phase surface hopping, and mean-field dynamics with surface hopping. We also describe in detail how a new nonadiabatic algorithm, mean-field dynamics with stochastic decoherence (MF-SD), is to be implemented for condensed-phase problems, and we apply MF-SD to the excited-state relaxation of the hydrated electron. Our discussion emphasizes the different ways quantum decoherence is treated in each algorithm and the resulting implications for hydrated-electron relaxation dynamics. We find that for three MQC methods that use Tully's fewest-switches criterion to determine surface hopping probabilities, the excited-state lifetime of the electron is the same. Moreover, the nonequilibrium solvent response function of the excited hydrated electron is the same with all of the nonadiabatic MQC algorithms discussed here, so that all of the algorithms would produce similar agreement with experiment. Despite the identical solvent response predicted by each MQC algorithm, we find that MF-SD allows much more mixing of multiple basis states into the quantum wave function than do other methods. This leads to an excited-state lifetime that is longer with MF-SD than with any method that incorporates nonadiabatic effects with the fewest-switches surface hopping criterion.     

The guanine tautomer puzzle: quantum chemical investigation of ground and excited states

Combined density functional and multireference configuration interaction methods have been employed to explore the ground and low-lying electronically excited states of the most important tautomeric and rotameric forms of guanine with the purpose of resolving the conflicting assignments of IR-UV bands found in the literature. The calculations predict sharp 1(pi-->pi*) origin transitions for the RN1 rotamer of the 7H-amino-hydroxy species and the RN7 rotamer of the 9H-amino-hydroxy species. The other 9H-amino-hydroxy rotamer, RN1, undergoes ultrafast nonradiative decay and is thus missing in the UV spectra. Because of its very small Franck-Condon factor and the presence of a conical intersection close by, it appears questionable, whether the 1(pi-->pi*) origin transition of 9H-amino-oxo-guanine can be observed experimentally. Vibrational overlap is more favorable for the 1(pi-->pi*) origin transition of the 7H- amino-oxo form, but also this tautomer is predicted to undergo ultrafast nonradiative decay of the 1(pi-->pi*) population. The good agreement of calculated IR frequencies of the amino-oxo species with recent IR spectra in He droplets and their mismatch with peaks observed in IR-UV spectra indicate that none of the bands stem from 7H- or 9H-amino-oxo guanine. Instead, our results suggest that these bands originate from 7H-imino-oxo guanine tautomers. In the excited-state dynamics of the biologically relevant 9H-amino-oxo tautomer, a diffuse charge transfer state is predicted to play a significant role.     

The photoisomerization mechanism of azobenzene: a semiclassical simulation of nonadiabatic dynamics

We have simulated the photoisomerization dynamics of azobenzene, taking into account internal conversion and geometrical relaxation processes, by means of a semiclassical surface hopping approach. Both n-->pi* and pi-->pi* excitations and both cis-->trans and trans-->cis conversions have been considered. We show that in all cases the torsion around the N==N double bond is the preferred mechanism. The quantum yields measured are correctly reproduced and the observed differences are explained as a result of the competition between the inertia of the torsional motion and the premature deactivation of the excited state. Recent time-resolved spectroscopic experiments are interpreted in the light of the simulated dynamics.     

The reaction rate for dissociative adsorption of N2 on stepped Ru(0001): six-dimensional quantum calculations

Quantum-mechanical calculations of the reaction rate for dissociative adsorption of N2 on stepped Ru(0001) are presented. Converged six-dimensional quantum calculations for this heavy-atom reaction have been performed using the multiconfiguration time-dependent Hartree method. A potential-energy surface for the transition-state region is constructed from density-functional theory calculations using Shepard interpolation. The quantum results are in very good agreement with the results of the harmonic transition-state theory. In contrast to the findings of previous model calculations on similar systems, the tunneling effect is found to be small.     

Optimized Jastrow-Slater wave functions for ground and excited states: application to the lowest states of ethene

A quantum Monte Carlo method is presented for determining multideterminantal Jastrow-Slater wave functions for which the energy is stationary with respect to the simultaneous optimization of orbitals and configuration interaction coefficients. The approach is within the framework of the so-called energy fluctuation potential method which minimizes the energy in an iterative fashion based on Monte Carlo sampling and a fitting of the local energy fluctuations. The optimization of the orbitals is combined with the optimization of the configuration interaction coefficients through the use of additional single excitations to a set of external orbitals. A new set of orbitals is then obtained from the natural orbitals of this enlarged configuration interaction expansion. For excited states, the approach is extended to treat the average of several states within the same irreducible representation of the pointgroup of the molecule. The relationship of our optimization method with the stochastic reconfiguration technique by Sorella et al. is examined. Finally, the performance of our approach is illustrated with the lowest states of ethene, in particular with the difficult case of the 1(1)B(1u) state.