2009
DOI: 10.1021/jp9063565
|View full text |Cite
|
Sign up to set email alerts
|

Ab Initio Multiple Spawning Dynamics Using Multi-State Second-Order Perturbation Theory

Abstract: We have implemented multi-state second-order perturbation theory (MS-CASPT2) in the ab initio multiple spawning (AIMS) method for first-principles molecular dynamics including nonadiabatic effects. The nonadiabatic couplings between states are calculated numerically using an efficient method which requires only two extra energy calculations per time step. As a representative example, we carry out AIMS-MSPT2 calculations of the excited state dynamics of ethylene. Two distinct types of conical intersections, pre… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
1
1
1

Citation Types

11
169
0
1

Year Published

2011
2011
2020
2020

Publication Types

Select...
6
2

Relationship

1
7

Authors

Journals

citations
Cited by 149 publications
(181 citation statements)
references
References 49 publications
11
169
0
1
Order By: Relevance
“…The directions that define the branching space for intersections involving a hydrogen bridge across the double bond are presented in Figure 2b. While the gand h-directions for CHE and CHD are completely analogous, they are quite different from the [ 26 further evidence that both of these molecules would be expected to undergo similar excitedstate dynamical processes.…”
Section: Resultsmentioning
confidence: 89%
See 1 more Smart Citation
“…The directions that define the branching space for intersections involving a hydrogen bridge across the double bond are presented in Figure 2b. While the gand h-directions for CHE and CHD are completely analogous, they are quite different from the [ 26 further evidence that both of these molecules would be expected to undergo similar excitedstate dynamical processes.…”
Section: Resultsmentioning
confidence: 89%
“…The [1,2] H-shift and [1,2] H-bridge intersections in CHE and CHD correspond to the previously discussed ethylidene-like and H-bridged intersections, respectively. 21,22,25,26,53 Likewise, in closed-ring systems, twistedpyramidalized geometries may involve the breaking of the adjacent CÀC bonds, as seen in the α-and β-cleavage processes for CHE and the single-bond breaking path in CHD. Obviously, there is no clear analogue to the [1,3] H-migration pathway observed in ethylene; however, the efficacy of this intersection type is likely highly dependent both on the twist-pyramidalization displacements required to promote hydrogen bond cleavage and migration, as well as on the ring structure that provides an initial conformationally locked geometry.…”
Section: Resultsmentioning
confidence: 99%
“…As ethylene is the simplest molecule with a C=C double bond, it has been extensively studied both experimentally and computationally. 22,[50][51][52][53][54][55][56][57][58] The AIMC and AIMS simulations were carried out with a modified version of AIMS-MOLPRO, 36 which was extended to include Ehrenfest dynamics for the TBFs. Electronic structure calculations were performed with the CASSCF method using an active space of two electrons in two orbitals and equally weighted stateaveraging over the three lowest singlet states, i.e., SA-3-CAS(2/2).…”
Section: Computational Details and Resultsmentioning
confidence: 99%
“…[49][50][51] With the above prescribed equations for the evolution of the phase and Ehrenfest amplitudes, we can now obtain the equations of motion for the TBF amplitudes c n by inserting the ansatz into the time-dependent Schrödinger equation, giving…”
Section: A Working Equationsmentioning
confidence: 99%
“…In AIMS, this flexibility is attained by adaptively increasing the number of TBFs, in a process called spawning, when the coupling of two states rises above a predefined threshold. As discussed in detail previously, 15,31,34 the nonadiabatic coupling vectors are used in AIMS in three ways: (1) to construct the Hamiltonian matrix elements in eq 5, (2) to decide when spawning occurs, and (3) to determine the momentum of the newly created child TBF's following spawning. In the first two uses, one does not need the entire nonadiabatic coupling vector but rather only the dot product with a (complex) off-diagonal matrix element of the momentum operator.…”
Section: Theoretical Methodsmentioning
confidence: 99%