1
|
Lee S, Shostak S, Filatov M, Choi CH. Conical Intersections in Organic Molecules: Benchmarking Mixed-Reference Spin-Flip Time-Dependent DFT (MRSF-TD-DFT) vs Spin-Flip TD-DFT. J Phys Chem A 2019; 123:6455-6462. [PMID: 31283235 DOI: 10.1021/acs.jpca.9b06142] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The mixed-reference spin-flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S1/S0 conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular single-reference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.
Collapse
Affiliation(s)
- Seunghoon Lee
- Department of Chemistry , Seoul National University , Seoul 151-747 , South Korea
| | - Svetlana Shostak
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| | - Michael Filatov
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| | - Cheol Ho Choi
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| |
Collapse
|
2
|
Nikiforov A, Gamez JA, Thiel W, Huix-Rotllant M, Filatov M. Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules. J Chem Phys 2015; 141:124122. [PMID: 25273427 DOI: 10.1063/1.4896372] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.
Collapse
Affiliation(s)
- Alexander Nikiforov
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Jose A Gamez
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt am Main, Germany
| | - Michael Filatov
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
3
|
Huix-Rotllant M, Nikiforov A, Thiel W, Filatov M. Description of Conical Intersections with Density Functional Methods. Top Curr Chem (Cham) 2015; 368:445-76. [PMID: 25896441 DOI: 10.1007/128_2015_631] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.
Collapse
Affiliation(s)
- Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438, Frankfurt am Main, Germany
| | | | | | | |
Collapse
|
4
|
Chiang SY, Fang YS, Bahou M, Sankaran K. Experiments and Calculations on Photoionization and Dissociative Photoionization of CH2CO. J CHIN CHEM SOC-TAIP 2013. [DOI: 10.1002/jccs.200400103] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
5
|
Liu Y, Yu JK, Huang XR, Sun CC. Theoretical study of photodissociation dynamics on the lowest-lying Rydberg state of ketene. J Chem Phys 2006; 125:44311. [PMID: 16942146 DOI: 10.1063/1.2220568] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
In the present study, an attempt is made to reveal the main mechanism of photodissociation on the lowest-lying Rydberg state (1)B(1) of ketene, referred to as the second singlet excited state S(2), by means of the complete active space self-consistent field and the second-order multiconfigurational perturbation theory methods. The located S(2)S(1)T(1) three-surface intersection plays an important role in the dissociation process. It is shown that the intersection permits an efficient internal conversion from S(2) to S(1) state, but prohibits the intersystem crossing from S(2) to T(1) state because of the small spin-orbital coupling value of 0.136 cm(-1). The main photodissociation process could be described as follows: after one photon absorption to the S(2) state, ketene preferentially relaxes to the minimum S(2)C(2v), and undergoes a transition state S(2)TS with small potential barrier along the C(s)-I (out-of-plane bent) symmetry, and passes through the S(2)S(1)T(1) intersection to reach S(1) surface, then arrives at the transition state S(1)TS along the minimum energy path. As is well known, S(1)-->S(0) internal conversion around the Franck-Condon region is expected to be very efficient, and eventually the hot S(0) molecule has accumulated enough energy to yield the CH(2) (a (1)A(1)) and CO (X (1)Sigma(+)) products.
Collapse
Affiliation(s)
- Yang Liu
- State Key Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, People's Republic of China
| | | | | | | |
Collapse
|
6
|
Liu J, Wang FY, Wang H, Jiang B, Yang XM. A Velocity Map Ion-imaging Study on Ketene Photodissociation at 218 nm. CHINESE J CHEM PHYS 2006. [DOI: 10.1360/cjcp2006.19(1).1.5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
7
|
Lu IC, Lee SH, Lee YT, Yang X. Photodissociation dynamics of ketene at 157.6nm. J Chem Phys 2006; 124:024324. [PMID: 16422600 DOI: 10.1063/1.2147221] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Photodissociation dynamics of ketene at 157.6 nm has been investigated using the photofragment translational spectroscopic technique based on photoionization detection using vacuum-ultraviolet synchrotron radiation. Three dissociation channels have been observed: CH2+CO, CH+HCO, and HCCO+H. The product translational energy distributions and angular anisotropy parameters were measured for all three observed dissociation channels, and the relative branching ratios for different channels were also estimated. The experimental results show that the direct C-C bond cleavage (CH2+CO) is the dominant channel, while H migration and elimination channels are very minor. The results in this work show that direct dissociation on excited electronic state is much more significant than the indirect dissociation via the ground state in the ketene photodissociation at 157.6 nm.
Collapse
Affiliation(s)
- I-Chung Lu
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| | | | | | | |
Collapse
|
8
|
Liu J, Wang F, Wang H, Jiang B, Yang X. A velocity map ion-imaging study on ketene photodissociation at 208 and 213nm: Rotational dependence of product angular anisotropy. J Chem Phys 2005; 122:104309. [PMID: 15836317 DOI: 10.1063/1.1858435] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Photodissociation dynamics of ketene following excitation at 208.59 and 213.24 nm have been investigated using the velocity map ion-imaging method. Both the angular distribution and translational energy distribution of the CO products at different rotational and vibrational states have been obtained. No significant difference in the translational energy distributions for different CO rotational state products has been observed at both excitation wavelengths. The anisotropy parameter beta is, however, noticeably different for different CO rotational state products at both excitation wavelengths. For lower rotational states of the CO product, beta is smaller than zero, while beta is larger than zero for CO at higher rotational states. The observed rotational dependence of angular anisotropy is interpreted as the dynamical influence of a peculiar conical intersection between the (1)B(1) excited state and (1)A(2) state along the C(S)-I coordinate.
Collapse
Affiliation(s)
- Jie Liu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning 116023, People's Republic of China
| | | | | | | | | |
Collapse
|
9
|
Wang S, Shi Y, Dénommée S, Simard B, Lee YP. Investigation of some Rydberg states of ketene by two-photon resonance-enhanced multiphoton ionization spectroscopy. J Chem Phys 2003. [DOI: 10.1063/1.1609392] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
|
10
|
Feltham EJ, Qadiri RH, Cottrill EEH, Cook PA, Cole JP, Balint-Kurti GG, Ashfold MNR. Ketene photodissociation in the wavelength range 193–215 nm: The H atom production channel. J Chem Phys 2003. [DOI: 10.1063/1.1601604] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|