1
|
Shu Y, Varga Z, Kanchanakungwankul S, Zhang L, Truhlar DG. Diabatic States of Molecules. J Phys Chem A 2022; 126:992-1018. [PMID: 35138102 DOI: 10.1021/acs.jpca.1c10583] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
2
|
Dresselhaus T, Bungey CBA, Knowles PJ, Manby FR. Coupling electrons and vibrations in molecular quantum chemistry. J Chem Phys 2020; 153:214114. [PMID: 33291918 DOI: 10.1063/5.0032900] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
We derive an electron-vibration model Hamiltonian in a quantum chemical framework and explore the extent to which such a Hamiltonian can capture key effects of nonadiabatic dynamics. The model Hamiltonian is a simple two-body operator, and we make preliminary steps at applying standard quantum chemical methods to evaluate its properties, including mean-field theory, linear response, and a primitive correlated model. The Hamiltonian can be compared to standard vibronic Hamiltonians, but it is constructed without reference to potential energy surfaces through direct differentiation of the one- and two-electron integrals at a single reference geometry. The nature of the model Hamiltonian in the harmonic and linear-coupling regime is investigated for pyrazine, where a simple time-dependent calculation including electron-vibration correlation is demonstrated to exhibit the well-studied population transfer between the S2 and S1 excited states.
Collapse
Affiliation(s)
- Thomas Dresselhaus
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Callum B A Bungey
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
3
|
Sibaev M, Polyak I, Manby FR, Knowles PJ. Molecular second-quantized Hamiltonian: Electron correlation and non-adiabatic coupling treated on an equal footing. J Chem Phys 2020; 153:124102. [DOI: 10.1063/5.0018930] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R. Manby
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Peter J. Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
4
|
Richings GW, Habershon S. A new diabatization scheme for direct quantum dynamics: Procrustes diabatization. J Chem Phys 2020; 152:154108. [DOI: 10.1063/5.0003254] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
5
|
Accomasso D, Persico M, Granucci G. Diabatization by Localization in the Framework of Configuration Interaction Based on Floating Occupation Molecular Orbitals (FOMO−CI). CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900056] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Davide Accomasso
- Dipartimento di Chimica e Chimica IndustrialeUniversitá di Pisa v. G. Moruzzi 13 I-56124 Pisa Italy
| | - Maurizio Persico
- Dipartimento di Chimica e Chimica IndustrialeUniversitá di Pisa v. G. Moruzzi 13 I-56124 Pisa Italy
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica IndustrialeUniversitá di Pisa v. G. Moruzzi 13 I-56124 Pisa Italy
| |
Collapse
|
6
|
Ryabinkin IG, Hsieh CY, Kapral R, Izmaylov AF. Analysis of geometric phase effects in the quantum-classical Liouville formalism. J Chem Phys 2014; 140:084104. [PMID: 24588145 DOI: 10.1063/1.4866366] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. As a consequence, GP effects in WA-QCL can be associated with a dynamical term in the corresponding equation of motion. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.
Collapse
Affiliation(s)
- Ilya G Ryabinkin
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| | - Chang-Yu Hsieh
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Raymond Kapral
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| |
Collapse
|
7
|
An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1526-1] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
8
|
|
9
|
Sajeev Y, Moiseyev N. Enhancement of charge transition from donor to acceptor in organic and biochemical reactions by the intermediate-mode-assisted tunneling mechanism. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.01.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
10
|
Paul AK, Adhikari S, Baer M. Space-time contours to treat intense field-dressed molecular states. J Chem Phys 2010; 132:034303. [PMID: 20095735 DOI: 10.1063/1.3282333] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this article we consider a molecular system exposed to an intense short-pulsed external field. It is a continuation of a previous publication [A. K. Paul, S. Adhikari, D. Mukhopadhyay et al., J. Phys. Chem. A 113, 7331 (2009)] in which a theory is presented that treats quantum effects due to nonclassical photon states (known also as Fock states). Since these states became recently a subject of intense experimental efforts we thought that they can be treated properly within the existing quantum formulation of dynamical processes. This was achieved by incorporating them in the Born-Oppenheimer (BO) treatment with time-dependent coefficients. The extension of the BO treatment to include the Fock states results in a formidable enhancement in numerical efforts expressed, in particular, in a significant increase in CPU time. In the present article we discuss an approach that yields an efficient and reliable approximation with only negligible losses in accuracy. The approximation is tested in detail for the dissociation process of H(2) (+) as caused by a laser field.
Collapse
Affiliation(s)
- Amit K Paul
- Department of Physical Chemistry, Indian Association for Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | | | | |
Collapse
|
11
|
Nagae H, Kakitani Y, Koyama Y. Theoretical description of diabatic mixing and coherent excitation in singlet-excited states of carotenoids. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.04.039] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
12
|
Sarkar B, Adhikari S, Baer M. Space-time contours to treat intense field-dressed molecular states. I. Theory. J Chem Phys 2007; 127:014301. [PMID: 17627340 DOI: 10.1063/1.2743437] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A molecular system exposed to an intense external field is considered. The strength of the field is measured by the number L of electronic states that become populated during this process. In the present article the authors discuss a rigorous way, based on the recently introduced space-time contours [R. Baer, et al., J. Chem. Phys. 119, 6998 (2003)], to form N coupled Schrodinger equations where N<L, which maintains the effects due to the remaining (L-N) populated states. It is shown that whereas the size of L is unlimited, the main requirement concerning N is that the original group of N field-free states forms a Hilbert subspace in the spatial region of interest. From previous studies it is known that a group of states forms a Hilbert subspace if and only if the corresponding topological D matrix is diagonal [M. Baer, et al., Farad, Discuss 127, 337 (2004)].
Collapse
Affiliation(s)
- Biplab Sarkar
- Department of Chemistry, Indian Institute of Technology, Guwahati, North Guwahati, Guwahati 781039, India
| | | | | |
Collapse
|
13
|
Valero R, Truhlar DG. Nonadiabatic effects in C–Br bond scission in the photodissociation of bromoacetyl chloride. J Chem Phys 2006; 125:194305. [PMID: 17129101 DOI: 10.1063/1.2363991] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Bromoacetyl chloride photodissociation has been interpreted as a paradigmatic example of a process in which nonadiabatic effects play a major role. In molecular beam experiments by Butler and co-workers [J. Chem. Phys. 95, 3848 (1991); J. Chem. Phys. 97, 355 (1992)], BrCH2C(O)Cl was prepared in its ground electronic state (S0) and excited with a laser at 248 nm to its first excited singlet state (S1). The two main ensuing photoreactions are the ruptures of the C-Cl bond and of the C-Br bond. A nonadiabatic model was proposed in which the C-Br scission is strongly suppressed due to nonadiabatic recrossing at the barrier formed by the avoided crossing between the S1 and S2 states. Recent reduced-dimensional dynamical studies lend support to this model. However, another interpretation that has been given for the experimental results is that the reduced probability of C-Br scission is a consequence of incomplete intramolecular energy redistribution. To provide further insight into this problem, we have studied the energetically lowest six singlet electronic states of bromoacetyl chloride by using an ab initio multiconfigurational perturbative electronic structure method. Stationary points (minima and saddle points) and minimum energy paths have been characterized on the S0 and S1 potential energy surfaces. The fourfold way diabatization method has been applied to transform five adiabatic excited electronic states to a diabatic representation. The diabatic potential energy matrix of the first five excited singlet states has been constructed along several cuts of the potential energy hypersurfaces. The thermochemistry of the photodissociation reactions and a comparison with experimental translational energy distributions strongly suggest that nonadiabatic effects dominate the C-Br scission, but that the reaction proceeds along the energetically allowed diabatic pathway to excited-state products instead of being nonadiabatically suppressed. This conclusion is also supported by the low values of the diabatic couplings on the C-Br scission reaction path. The methodology established in the present study will be used for the construction of global potential energy surfaces suitable for multidimensional dynamics simulations to test these preliminary interpretations.
Collapse
Affiliation(s)
- Rosendo Valero
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | |
Collapse
|
14
|
Neugebauer J, Baerends EJ, Nooijen M. Vibronic coupling and double excitations in linear response time-dependent density functional calculations: Dipole-allowed states of N2. J Chem Phys 2004; 121:6155-66. [PMID: 15446909 DOI: 10.1063/1.1785775] [Citation(s) in RCA: 73] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The present study serves two purposes. First, we evaluate the ability of present time-dependent density functional response theory (TDDFRT) methods to deal with avoided crossings, i.e., vibronic coupling effects. In the second place, taking the vibronic coupling effects into account enables us, by comparison to the configuration analysis in a recent ab initio study [J. Chem. Phys. 115, 6438 (2001)], to identify the neglect of double excitations as the prime cause of limited accuracy of these linear response based TDDFRT calculations for specific states. The "statistical averaging of (model) orbital potentials (SAOP)" Kohn-Sham potential is used together with the standard adiabatic local-density approximation (ALDA) for the exchange-correlation kernel. We use the N2 molecule as prototype, since the TDDFRT/SAOP calculations have already been shown to be accurate for the vertical excitations, while this molecule has a well-studied, intricate vibronic structure as well as significant double excitation nature in the lowest 1Pi(u) state at elongated bond lengths. A simple diabatizing scheme is employed to obtain a diabatic potential energy matrix, from which we obtain the absorption spectrum of N2 including vibronic coupling effects. Considering the six lowest dipole allowed transitions of 1Sigma(u)+ and 1Pi(u) symmetry, we observe a good general agreement and conclude that avoided crossings and vibronic coupling can indeed be treated satisfactorily on the basis of TDDFRT excitation energies. However, there is one state for which the accuracy of TDDFRT/ALDA clearly breaks down. This is the state for which the ab initio calculations find significant double excitation character. To deal with double excitation character is an important challenge for time-dependent density functional theory.
Collapse
Affiliation(s)
- Johannes Neugebauer
- Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands.
| | | | | |
Collapse
|
15
|
Baer M, Vértesi T, Halász GJ, Vibók Á. Electronic Diabatic Framework: Restrictions Due to Quantization of the Nonadiabatic Coupling Matrix. J Phys Chem A 2004. [DOI: 10.1021/jp0487051] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- M. Baer
- Soreq Nuclear Research Center, Yavne 81800, Israel, Department of Theoretical Physics and Institute of Informatics, University of Debrecen, Debrecen, Hungary
| | - T. Vértesi
- Soreq Nuclear Research Center, Yavne 81800, Israel, Department of Theoretical Physics and Institute of Informatics, University of Debrecen, Debrecen, Hungary
| | - G. J. Halász
- Soreq Nuclear Research Center, Yavne 81800, Israel, Department of Theoretical Physics and Institute of Informatics, University of Debrecen, Debrecen, Hungary
| | - Á. Vibók
- Soreq Nuclear Research Center, Yavne 81800, Israel, Department of Theoretical Physics and Institute of Informatics, University of Debrecen, Debrecen, Hungary
| |
Collapse
|
16
|
Vértesi T, Vibók A, Halász GJ, Baer M. A field theoretical approach to calculate electronic Born-Oppenheimer coupling terms. J Chem Phys 2004; 121:4000-13. [PMID: 15332946 DOI: 10.1063/1.1778691] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper we suggest to consider the spatial distribution of the Born-Oppenheimer nonadiabatic coupling terms as fields which are created by sources, located at degeneracy points, and which can be derived using the ordinary mathematical tools of field theory. It is shown that the curl-divergence equations as formed within a given Hilbert space [M. Baer, Chem. Phys. Lett. 35, 112 (1975)] can be converted into a set of inhomogeneous coupled Poisson equations which are solved for a given set of boundary conditions. The method is applied to the three-state Hilbert subspace of the H(3) system. The numerical results are compared with ab initio calculations for which a very encouraging fit is found.
Collapse
Affiliation(s)
- T Vértesi
- Department of Theoretical Physics, University of Debrecen, Debrecen, Hungary
| | | | | | | |
Collapse
|