1
|
Neville SP, Schuurman MS. Calculation of quasi-diabatic states within the DFT/MRCI(2) framework: The QD-DFT/MRCI(2) method. J Chem Phys 2024; 160:234109. [PMID: 38899687 DOI: 10.1063/5.0214637] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Accepted: 05/23/2024] [Indexed: 06/21/2024] Open
Abstract
We describe a procedure for the calculation of quasi-diabatic states within the recently introduced DFT/MRCI(2) framework [S. P. Neville and M. S. Schuurman, J. Chem. Phys. 157, 164103 (2022)]. Based on an effective Hamiltonian formalism, the proposed procedure, which we term QD-DFT/MRCI(2), has the advantageous characteristics of being simultaneously highly efficient and effectively black box in nature while directly yielding both quasi-diabatic potentials and wave functions of high quality. The accuracy and efficiency of the QD-DFT/MRCI(2) formalism are demonstrated via the simulation of the vibronic absorption spectra of furan and chlorophyll a.
Collapse
Affiliation(s)
- Simon P Neville
- National Research Council Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
| | - Michael S Schuurman
- National Research Council Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| |
Collapse
|
2
|
Weike N, Viel A, Eisfeld W. Hydrogen-iodine scattering. I. Development of an accurate spin-orbit coupled diabatic potential energy model. J Chem Phys 2023; 159:244119. [PMID: 38156638 DOI: 10.1063/5.0186787] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Accepted: 12/05/2023] [Indexed: 01/03/2024] Open
Abstract
The scattering of H by I is a prototypical model system for light-heavy scattering in which relativistic coupling effects must be taken into account. Scattering calculations depend strongly on the accuracy of the potential energy surface (PES) model. The methodology to obtain such an accurate PES model suitable for scattering calculations is presented, which includes spin-orbit (SO) coupling within the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) approach. In this approach, the SO coupling is determined only for the atomic states of the heavy atom, and the geometry dependence of the SO effect is accounted for by a diabatization with respect to asymptotic states. The accuracy of the full model, composed of a Coulomb part and the SO model, is achieved in the following ways. For the SO model, the extended ERCAR approach is applied, which accounts for both intra-state and inter-state SO coupling, and an extended number of diabatic states are included. The corresponding coupling constants for the SO operator are obtained from experiments, which are more accurate than computed values. In the Coulomb Hamiltonian model, special attention is paid to the long range behavior and accurate c6 dispersion coefficients. The flexibility and accuracy of this Coulomb model are achieved by combining partial models for three different regions. These are merged via artificial neural networks, which also refine the model further. In this way, an extremely accurate PES model for hydrogen iodide is obtained, suitable for accurate scattering calculations.
Collapse
Affiliation(s)
- Nicole Weike
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Alexandra Viel
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes)-UMR 6251, F-35000 Rennes, France
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
3
|
Littlejohn R, Rawlinson J, Subotnik J. Representation and conservation of angular momentum in the Born-Oppenheimer theory of polyatomic molecules. J Chem Phys 2023; 158:104302. [PMID: 36922131 DOI: 10.1063/5.0143809] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023] Open
Abstract
This paper concerns the representation of angular momentum operators in the Born-Oppenheimer theory of polyatomic molecules and the various forms of the associated conservation laws. Topics addressed include the question of whether these conservation laws are exactly equivalent or only to some order of the Born-Oppenheimer parameter κ = (m/M)1/4 and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed in both problems involving a single potential energy surface and those with multiple, strongly coupled surfaces and in both the electrostatic model and those for which fine structure and electron spin are important. The analysis leads to an examination of the transformation laws under rotations of the electronic Hamiltonian; of the basis states, both adiabatic and diabatic, along with their phase conventions; of the potential energy matrix; and of the derivative couplings. These transformation laws are placed in the geometrical context of the structures in the nuclear configuration space that are induced by rotations, which include the rotational orbits or fibers, the surfaces upon which the orientation of the molecule changes but not its shape, and the section, an initial value surface that cuts transversally through the fibers. Finally, it is suggested that the usual Born-Oppenheimer approximation can be replaced by a dressing transformation, that is, a sequence of unitary transformations that block-diagonalize the Hamiltonian. When the dressing transformation is carried out, we find that the angular momentum operator does not change. This is a part of a system of exact equivalences among various representations of angular momentum operators in Born-Oppenheimer theory. Our analysis accommodates large-amplitude motions and is not dependent on small-amplitude expansions about an equilibrium position. Our analysis applies to noncollinear configurations of a polyatomic molecule; this covers all but a subset of measure zero (the collinear configurations) in the nuclear configuration space.
Collapse
Affiliation(s)
- Robert Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Jonathan Rawlinson
- School of Mathematics, University of Manchester, Manchester, United Kingdom
| | - Joseph Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
4
|
Weike N, Chanut E, Hoppe H, Eisfeld W. Development of a fully coupled diabatic spin-orbit model for the photodissociation of phenyl iodide. J Chem Phys 2022; 156:224109. [PMID: 35705416 DOI: 10.1063/5.0088205] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The theoretical treatment of the quantum dynamics of the phenyl iodide photodissociation requires an accurate analytical potential energy surface (PES) model. This model must also account for spin-orbit (SO) coupling. This study is the first step to construct accurate SO coupled PESs, namely, for the C-I dissociation coordinate. The model is based on the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) method developed over the past ten years. The SO-free Hamiltonian is represented in an asymptotic diabatic basis and then combined with an atomic effective relativistic coupling operator determined analytically. In contrast to the previously studied cases (HI, CH3I), the diabatic basis states are due to excitations in the phenyl fragment rather than the iodine atom. An accurate analytical model of the ab initio reference data is determined in two steps. The first step is a simple reference model describing the data qualitatively. This reference model is corrected through a trained artificial neural-network to achieve high accuracy. The SO-free and the fine structure states resulting from this ERCAR model are discussed extensively in the context of the photodissociation.
Collapse
Affiliation(s)
- Nicole Weike
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Emma Chanut
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Hannes Hoppe
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
5
|
Rey M. Novel methodology for systematically constructing global effective models from ab initio-based surfaces: A new insight into high-resolution molecular spectra analysis. J Chem Phys 2022; 156:224103. [PMID: 35705402 DOI: 10.1063/5.0089097] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, a novel methodology is presented for the construction of ab initio effective rotation-vibration spectroscopic models from potential energy and dipole moment surfaces. Non-empirical effective Hamiltonians are obtained via the block-diagonalization of selected variationally computed eigenvector matrices. For the first time, the derivation of an effective dipole moment is carried out in a systematic way. This general approach can be implemented quite easily in most of the variational computer codes and turns out to be a clear alternative to the rather involved Van Vleck perturbation method. Symmetry is exploited at all stages to translate first-principles calculations into a set of spectroscopic parameters to be further refined on experiment. We demonstrate on H2CO, PH3, CH4, C2H4, and SF6 that the proposed effective model can provide crucial information to spectroscopists within a very short time compared to empirical spectroscopic models. This approach brings a new insight into high-resolution spectrum analysis of polyatomic molecules and will be also of great help in the modeling of hot atmospheres where completeness is important.
Collapse
Affiliation(s)
- Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| |
Collapse
|
6
|
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: 31] [Impact Index Per Article: 15.5] [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
|
7
|
Richings GW, Habershon S. Analyzing Grid-Based Direct Quantum Molecular Dynamics Using Non-Linear Dimensionality Reduction. Molecules 2021; 26:molecules26247418. [PMID: 34946499 PMCID: PMC8708769 DOI: 10.3390/molecules26247418] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Revised: 11/30/2021] [Accepted: 12/02/2021] [Indexed: 11/18/2022] Open
Abstract
Grid-based schemes for simulating quantum dynamics, such as the multi-configuration time-dependent Hartree (MCTDH) method, provide highly accurate predictions of the coupled nuclear and electronic dynamics in molecular systems. Such approaches provide a multi-dimensional, time-dependent view of the system wavefunction represented on a coordinate grid; in the case of non-adiabatic simulations, additional information about the state populations adds a further layer of complexity. As such, wavepacket motion on potential energy surfaces which couple many nuclear and electronic degrees-of-freedom can be extremely challenging to analyse in order to extract physical insight beyond the usual expectation-value picture. Here, we show that non-linear dimensionality reduction (NLDR) methods, notably diffusion maps, can be adapted to extract information from grid-based wavefunction dynamics simulations, providing insight into key nuclear motions which explain the observed dynamics. This approach is demonstrated for 2-D and 9-D models of proton transfer in salicylaldimine, as well as 8-D and full 12-D simulations of cis-trans isomerization in ethene; these simulations demonstrate how NLDR can provide alternative views of wavefunction dynamics, and also highlight future developments.
Collapse
|
8
|
Bai X, Guo X, Wang L. Machine Learning Approach to Calculate Electronic Couplings between Quasi-diabatic Molecular Orbitals: The Case of DNA. J Phys Chem Lett 2021; 12:10457-10464. [PMID: 34672582 DOI: 10.1021/acs.jpclett.1c03053] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Diabatization of one-electron states in flexible molecular aggregates is a great challenge due to the presence of surface crossings between molecular orbital (MO) levels and the complex interaction between MOs of neighboring molecules. In this work, we present an efficient machine learning approach to calculate electronic couplings between quasi-diabatic MOs without the need of nonadiabatic coupling calculations. Using MOs of rigid molecules as references, the MOs that can be directly regarded to be quasi-diabatic in molecular dynamics are selected out, state tracked, and phase corrected. On the basis of this information, artificial neural networks are trained to characterize the structure-dependent onsite energies of quasi-diabatic MOs and the intermolecular electronic couplings. A representative sequence of DNA is systematically studied as an illustration. Smooth time evolution of electronic couplings in all base pairs is obtained with quasi-diabatic MOs. In particular, our method can calculate electronic couplings between different quasi-diabatic MOs independently, and thus, this possesses unique advantages in many applications.
Collapse
Affiliation(s)
- Xin Bai
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
9
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 162] [Impact Index Per Article: 54.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
10
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
11
|
Milovanović B, Novak J, Etinski M, Domcke W, Došlić N. Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters. Phys Chem Chem Phys 2021; 23:2594-2604. [PMID: 33475644 DOI: 10.1039/d0cp05618a] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.
Collapse
Affiliation(s)
| | - Jurica Novak
- Department of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia. and Higher Medical and Biological School, South Ural State University, Chelyabinsk, Russia
| | - Mihajlo Etinski
- University of Belgrade, Faculty of Physical Chemistry, Belgrade, Serbia
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching, Germany
| | - Nađa Došlić
- Department of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia.
| |
Collapse
|
12
|
Viel A, Williams DMG, Eisfeld W. Accurate quantum dynamics simulation of the photodetachment spectrum of the nitrate anion (NO 3 -) based on an artificial neural network diabatic potential model. J Chem Phys 2021; 154:084302. [PMID: 33639724 DOI: 10.1063/5.0039503] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The photodetachment spectrum of the nitrate anion (NO3 -) is simulated from first principles using wavepacket quantum dynamics propagation and a newly developed accurate full-dimensional fully coupled five state diabatic potential model. This model utilizes the recently proposed complete nuclear permutation inversion invariant artificial neural network diabatization technique [D. M. G. Williams and W. Eisfeld, J. Phys. Chem. A 124, 7608 (2020)]. The quantum dynamics simulations are designed such that temperature effects and the impact of near threshold detachment are taken into account. Thus, the two available experiments at high temperature and at cryogenic temperature using the slow electron velocity-map imaging technique can be reproduced in very good agreement. These results clearly show the relevance of hot bands and vibronic coupling between the X̃ 2A2 ' ground state and the B̃ 2E' excited state of the neutral radical. This together with the recent experiment at low temperature gives further support for the proper assignment of the ν3 fundamental, which has been debated for many years. An assignment of a not yet discussed hot band line is also proposed.
Collapse
Affiliation(s)
- Alexandra Viel
- University Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, FranceTheoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - David M G Williams
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
13
|
Westermayr J, Marquetand P. Machine learning and excited-state molecular dynamics. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab9c3e] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
14
|
Williams DMG, Eisfeld W. Complete Nuclear Permutation Inversion Invariant Artificial Neural Network (CNPI-ANN) Diabatization for the Accurate Treatment of Vibronic Coupling Problems. J Phys Chem A 2020; 124:7608-7621. [DOI: 10.1021/acs.jpca.0c05991] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David M. G. Williams
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
15
|
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
|
16
|
Neville SP, Seidu I, Schuurman MS. Propagative block diagonalization diabatization of DFT/MRCI electronic states. J Chem Phys 2020; 152:114110. [DOI: 10.1063/1.5143126] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Simon P. Neville
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Issaka Seidu
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Michael S. Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
- National Research Council of Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
17
|
Sapunar M, Domcke W, Došlić N. UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects. Phys Chem Chem Phys 2019; 21:22782-22793. [PMID: 31595896 DOI: 10.1039/c9cp04662c] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.
Collapse
Affiliation(s)
- Marin Sapunar
- Department of Physical Chemistry, Ruder Bošković Institute, 10000 Zagreb, Croatia.
| | | | | |
Collapse
|
18
|
Sapunar M, Piteša T, Davidović D, Došlić N. Highly Efficient Algorithms for CIS Type Excited State Wave Function Overlaps. J Chem Theory Comput 2019; 15:3461-3469. [DOI: 10.1021/acs.jctc.9b00235] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marin Sapunar
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| | - Tomislav Piteša
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| | - Davor Davidović
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| | - Nadja Došlić
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| |
Collapse
|
19
|
Robertson C, González-Vázquez J, Corral I, Díaz-Tendero S, Díaz C. Nonadiabatic scattering of NO off Au3
clusters: A simple and robust diabatic state manifold generation method for multiconfigurational wavefunctions. J Comput Chem 2018; 40:794-810. [DOI: 10.1002/jcc.25764] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2018] [Revised: 10/28/2018] [Accepted: 11/01/2018] [Indexed: 01/19/2023]
Affiliation(s)
- Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing; University Of Warwick, CV4 7AL; Coventry United Kingdom
| | - Jesús González-Vázquez
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
| | - Ines Corral
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
| | - Sergio Díaz-Tendero
- Condensed Matter Physics Center (IFIMAC); Universidad Autónoma de Madrid, 28049; Madrid Spain
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
| | - Cristina Díaz
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
- Condensed Matter Physics Center (IFIMAC); Universidad Autónoma de Madrid, 28049; Madrid Spain
| |
Collapse
|
20
|
Jiang S, Xie Y, Lan Z. The role of the charge-transfer states in the ultrafast excitonic dynamics of the DTDCTB dimers embedded in a crystal environment. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.07.044] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
21
|
Liu M, Chen X, Grofe A, Gao J. Diabatic States at Construction (DAC) through Generalized Singular Value Decomposition. J Phys Chem Lett 2018; 9:6038-6046. [PMID: 30277783 DOI: 10.1021/acs.jpclett.8b02472] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
A procedure, called generalized diabatic-at-construction (GDAC), is presented to transform adiabatic potential energy surfaces into a diabatic representation by generalized singular value decomposition. First, we use a set of localized, valence bond-like configuration state functions, called DAC, as the basis states. Then, the adiabatic ground and relevant excited states are determined using multistate density functional theory (MSDFT). GDAC differs in the opposite direction from traditional approaches based on adiabatic-to-diabatic transformation with certain property restraints. The method is illustrated with applications to a model first-order bond dissociation reaction of CH3OCH2Cl polarized by a solvent molecule, the ground- and first-excited-state potential energy surfaces near the minimum conical intersection for the ammonia dimer photodissociation, and the multiple avoided curve crossings in the dissociation of lithium hydride. The GDAC diabatization method may be useful for defining charge-localized states in studies of electron transfer and proton-coupled electron transfer reactions in proteins.
Collapse
Affiliation(s)
- Meiyi Liu
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
| | - Xin Chen
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
| | - Adam Grofe
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
| | - Jiali Gao
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
- Department of Chemistry and Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| |
Collapse
|
22
|
|
23
|
Wittenbrink N, Eisfeld W. Extension of the effective relativistic coupling by asymptotic representation (ERCAR) approach to multi-dimensional potential energy surfaces: 3D model for CH3I. J Chem Phys 2018. [DOI: 10.1063/1.5011757] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nils Wittenbrink
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
24
|
Xie Y, Jiang S, Zheng J, Lan Z. Construction of Vibronic Diabatic Hamiltonian for Excited-State Electron and Energy Transfer Processes. J Phys Chem A 2017; 121:9567-9578. [DOI: 10.1021/acs.jpca.7b07737] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yu Xie
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
| | - Shengshi Jiang
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jie Zheng
- Industrial Research Institute of Nonwovens & Technical Textiles, College of Textiles & Clothing, Qingdao University, Qingdao 266071, China
| | - Zhenggang Lan
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
25
|
Wittenbrink N, Eisfeld W. An improved spin-orbit coupling model for use within the effective relativistic coupling by asymptotic representation (ERCAR) method. J Chem Phys 2017; 146:144110. [PMID: 28411618 DOI: 10.1063/1.4979949] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An improved atomic spin-orbit model is presented, which is designed to be used within the framework of the effective relativistic coupling by asymptotic representation method. This method is used for the generation of highly accurate coupled potential energy surfaces (PESs) to represent the fine structure energies of appropriate systems. The approach is demonstrated using CH3I and its photodissociation as a typical example. The method is based on a specific diabatization of electronic spin-space ("spin-free") states with respect to the asymptote at which a single relativistic atom is separated from a molecular non-relativistic fragment. Thus, the relativistic coupling effects can be treated entirely within the atomic framework. So far, an effective spin-orbit coupling operator which only accounts for intra-state coupling within each atomic spin-space state was used. In the present work, this approach is extended to account for inter-state couplings among different atomic spin-space states as well. It is shown that this extended approach improves the accuracy of the PESs significantly for higher excited states and also enhances the accuracy of low energy states. In particular, it improves the representation of the spin-orbit induced conical intersection among the 3Q0 and 1Q1 states of CH3I, which is of high relevance for the nonadiabatic quantum dynamics of the photodissociation.
Collapse
Affiliation(s)
- Nils Wittenbrink
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
26
|
Affiliation(s)
- Tao Zeng
- Department of Chemistry, Carleton University, Ottawa, Ontario K1S5B6, Canada
| |
Collapse
|
27
|
Grofe A, Qu Z, Truhlar DG, Li H, Gao J. Diabatic-At-Construction Method for Diabatic and Adiabatic Ground and Excited States Based on Multistate Density Functional Theory. J Chem Theory Comput 2017; 13:1176-1187. [PMID: 28135420 PMCID: PMC5793876 DOI: 10.1021/acs.jctc.6b01176] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We describe a diabatic-at-construction (DAC) strategy for defining diabatic states to determine the adiabatic ground and excited electronic states and their potential energy surfaces using the multistate density functional theory (MSDFT). The DAC approach differs in two fundamental ways from the adiabatic-to-diabatic (ATD) procedures that transform a set of preselected adiabatic electronic states to a new representation. (1) The DAC states are defined in the first computation step to form an active space, whose configuration interaction produces the adiabatic ground and excited states in the second step of MSDFT. Thus, they do not result from a similarity transformation of the adiabatic states as in the ATD procedure; they are the basis for producing the adiabatic states. The appropriateness and completeness of the DAC active space can be validated by comparison with experimental observables of the ground and excited states. (2) The DAC diabatic states are defined using the valence bond characters of the asymptotic dissociation limits of the adiabatic states of interest, and they are strictly maintained at all molecular geometries. Consequently, DAC diabatic states have specific and well-defined physical and chemical meanings that can be used for understanding the nature of the adiabatic states and their energetic components. Here we present results for the four lowest singlet states of LiH and compare them to a well-tested ATD diabatization method, namely the 3-fold way; the comparison reveals both similarities and differences between the ATD diabatic states and the orthogonalized DAC diabatic states. Furthermore, MSDFT can provide a quantitative description of the ground and excited states for LiH with multiple strongly and weakly avoided curve crossings spanning over 10 Å of interatomic separation.
Collapse
Affiliation(s)
- Adam Grofe
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Zexing Qu
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Hui Li
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
| | - Jiali Gao
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
28
|
Medders GR, Alguire EC, Jain A, Subotnik JE. Ultrafast Electronic Relaxation through a Conical Intersection: Nonadiabatic Dynamics Disentangled through an Oscillator Strength-Based Diabatization Framework. J Phys Chem A 2017; 121:1425-1434. [DOI: 10.1021/acs.jpca.6b12120] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Gregory R. Medders
- Department
of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Ethan C. Alguire
- Department
of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amber Jain
- Department
of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E. Subotnik
- Department
of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Stanford
PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
29
|
Eisfeld W, Viel A. Vibronic eigenstates and the geometric phase effect in the2E″ state of NO3. J Chem Phys 2017; 146:034303. [DOI: 10.1063/1.4973983] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Wittenbrink N, Venghaus F, Williams D, Eisfeld W. A new approach for the development of diabatic potential energy surfaces: Hybrid block-diagonalization and diabatization by ansatz. J Chem Phys 2016; 145:184108. [DOI: 10.1063/1.4967258] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nils Wittenbrink
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Florian Venghaus
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - David Williams
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
31
|
Dawes R, Ndengué SA. Single- and multireference electronic structure calculations for constructing potential energy surfaces. INT REV PHYS CHEM 2016. [DOI: 10.1080/0144235x.2016.1195102] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|