1
|
Sisodiya DS, Chattopadhyay A. The photochemical trans → cis and thermal cis → trans isomerization pathways of azobenzo-13-crown ether: A computational study on a strained cyclic azobenzene system. J Chem Phys 2024; 161:034307. [PMID: 39017425 DOI: 10.1063/5.0206946] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 06/26/2024] [Indexed: 07/18/2024] Open
Abstract
The isomerization of azobenzo-13-crown ether can be expected to be hindered due to the polyoxyethylene linkage connecting the 2,2'-positions of azobenzene. The mixed reference spin-flip time-dependent density functional theory results reveal that the planar and rotational minima of the first photo-excited singlet state (S1) of the trans-isomer pass through a barrier (2.5-5.0 kcal/mol) as it goes toward the torsional conical intersection (S0/S1) geometry (
Collapse
Affiliation(s)
- Dilawar Singh Sisodiya
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani, K. K. Birla Goa Campus, Zuarinagar, India
| | - Anjan Chattopadhyay
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani, K. K. Birla Goa Campus, Zuarinagar, India
| |
Collapse
|
2
|
Park W, Lashkaripour A, Komarov K, Lee S, Huix-Rotllant M, Choi CH. Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT. J Chem Theory Comput 2024; 20:5679-5694. [PMID: 38902891 DOI: 10.1021/acs.jctc.4c00640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2024]
Abstract
Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Alireza Lashkaripour
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
- Department of Chemistry, University of Zürich, Zürich 8057, Switzerland
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
3
|
Alías-Rodríguez M, Bonfrate S, Park W, Ferré N, Choi CH, Huix-Rotllant M. Solvent Effects and pH Dependence of the X-ray Absorption Spectra of Proline from Electrostatic Embedding Quantum Mechanics/Molecular Mechanics and Mixed-Reference Spin-Flip Time-dependent Density-Functional Theory. J Phys Chem A 2023. [PMID: 38019644 DOI: 10.1021/acs.jpca.3c05070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
The accurate description of solvent effects on X-ray absorption spectra (XAS) is fundamental for comparing the simulated spectra with experiments in solution. Currently, few protocols exist that can efficiently reproduce the effects of the solute/solvent interactions on XAS. Here, we develop an efficient and accurate theoretical protocol for simulating the solvent effects on XAS. The protocol combines electrostatic embedding QM/MM based on electrostatic potential fitted operators for describing the solute/solvent interactions and mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) for simulating accurate XAS spectra. To demonstrate the capabilities of our protocol, we compute the X-ray absorption of neutral proline in the gas phase and ionic proline in water in all relevant K-edges, showing excellent agreement with experiments. We show that states represented by core to π* transitions are almost unaffected by the interaction with water, whereas the core to σ* transitions are more impacted by the fluctuation of proline structure and the electrostatic interaction with the solvent. Finally, we reconstruct the pH-dependent XAS of proline in solution, determining that the N K-edge can be used to distinguish its three protonation states.
Collapse
Affiliation(s)
| | | | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Nicolas Ferré
- Aix-Marseille Univ, CNRS, ICR, Marseille 13013, France
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | | |
Collapse
|
4
|
Komarov K, Park W, Lee S, Huix-Rotllant M, Choi CH. Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory. J Chem Theory Comput 2023; 19:7671-7684. [PMID: 37844129 DOI: 10.1021/acs.jctc.3c00884] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.
Collapse
Affiliation(s)
- Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
| | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul, 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
5
|
Park W, Komarov K, Lee S, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory. J Phys Chem Lett 2023; 14:8896-8908. [PMID: 37767969 PMCID: PMC10561896 DOI: 10.1021/acs.jpclett.3c02296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 09/21/2023] [Indexed: 09/29/2023]
Abstract
The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.
Collapse
Affiliation(s)
- Woojin Park
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center
for Quantum Dynamics, Pohang University
of Science and Technology, Pohang 37673, South Korea
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Cheol Ho Choi
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
6
|
Chen X, Wang W, Xiao D, Xia SH, Zhang Y. Non-adiabatic dynamics simulations of the S 1 excited-state relaxation of diacetyl phenylenediamine. Phys Chem Chem Phys 2023. [PMID: 37427748 DOI: 10.1039/d3cp01826a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/11/2023]
Abstract
The small molecule built around the benzene ring, diacetyl phenylenediamine (DAPA), has attracted much attention due to its synthesis accessibility, large Stokes shift, etc. However, its meta structure m-DAPA does not fluoresce. In a previous investigation, it was found that such a property is due to the fact that it undergoes an energy-reasonable double proton transfer conical intersection during the deactivation of the S1 excited-state, then returns to the ground state by a nonradiative relaxation process eventually. However, our static electronic structure calculations and non-adiabatic dynamics analysis results indicate that only one reasonable non-adiabatic deactivation channel exists: after being excited to the S1 state, m-DAPA undergoes an ultrafast and barrierless ESIPT process and reaches the single-proton-transfer conical intersection. Subsequently, the system either returns to the keto-form S0 state minimum with proton reversion or returns to the single-proton-transfer S0 minimum after undergoing a slight twist of the acetyl group. The dynamics results show that the S1 excited-state lifetime of m-DAPA is 139 fs. In other words, we propose an efficient single-proton-transfer non-adiabatic deactivation channel of m-DAPA that is different from previous work, which can provide important mechanistic information of similar fluorescent materials.
Collapse
Affiliation(s)
- Xiaohang Chen
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China.
| | - Wei Wang
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China.
| | - Dongyi Xiao
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China.
| | - Shu-Hua Xia
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China.
| | - Yan Zhang
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China.
| |
Collapse
|
7
|
Komarov K, Mironov V, Lee S, Pham BQ, Gordon MS, Choi CH. High-performance strategies for the recent MRSF-TDDFT in GAMESS. J Chem Phys 2023; 158:2890476. [PMID: 37184015 DOI: 10.1063/5.0148005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 05/02/2023] [Indexed: 05/16/2023] Open
Abstract
Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.
Collapse
Affiliation(s)
- Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
| | - Vladimir Mironov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Buu Q Pham
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
8
|
Nakata H, Fedorov DG. Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method. J Chem Theory Comput 2023; 19:1276-1285. [PMID: 36753486 DOI: 10.1021/acs.jctc.2c01177] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The analytic energy gradient of energy with respect to nuclear coordinates is derived for the fragment molecular orbital (FMO) method combined with time-dependent density functional theory (TDDFT). The response terms arising from the use of a polarizable embedding are derived. The obtained analytic FMO-TDDFT gradient is shown to be accurate in comparison to both numerical FMO-TDDFT and unfragmented TDDFT gradients, at the level of two- and three-body expansions. The gradients are used for geometry optimizations, molecular dynamics, vibrational calculations, and simulations of IR and Raman spectra of excited states. The developed method is used to optimize the geometry of the ground and excited electronic states of the photoactive yellow protein (PDB: 2PHY).
Collapse
Affiliation(s)
- Hiroya Nakata
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
9
|
Komarov K, Park W, Lee S, Zeng T, Choi CH. Accurate Spin-Orbit Coupling by Relativistic Mixed-Reference Spin-Flip-TDDFT. J Chem Theory Comput 2023; 19:953-964. [PMID: 36655271 DOI: 10.1021/acs.jctc.2c01036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Relativistic mixed-reference spin-flip (MRSF)-TDDFT is developed considering the spin-orbit coupling (SOC) within the mean-field approximation. The resulting SOC-MRSF faithfully reproduces the experiments with very high accuracy, which is also consistent with the values by four-component (4c) relativistic CASSCF and 4c-CASPT2 in the spin-orbit-energy splitting calculations of the C, Si, and Ge atoms. Even for the fifth-row element Sn, the SOC-MRSF yielded accurate splittings (∼ 3 % error). In the SOC calculations of the molecular 4-thiothymine with a third-row element, SOC-MRSF values are in excellent agreement with those of the SO-GMC-QDPT2 level, regardless of geometries and exchange-correlation functionals. The same SOC-MRSF predicted the anticipated chance of S1 (nπ*) → T1 (ππ*) intersystem crossing, even in thymine with only second-row elements. With its accuracy and practicality, thus, SOC-MRSF is a promising electronic structure protocol in challenging situations such as nonadiabatic molecular dynamics (NAMD) incorporating both internal conversions and intersystem crossings in large systems.
Collapse
Affiliation(s)
- Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang37673, South Korea
| | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California91125, USA
| | - Tao Zeng
- Department of Chemistry, York University, Toronto, ONM3J 1P3, Canada
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu41566, South Korea
| |
Collapse
|
10
|
Shostak S, Park W, Oh J, Kim J, Lee S, Nam H, Filatov M, Kim D, Choi CH. Ultrafast Excited State Aromatization in Dihydroazulene. J Am Chem Soc 2023; 145:1638-1648. [PMID: 36633597 DOI: 10.1021/jacs.2c09800] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Excited-state aromatization dynamics in the photochemical ring opening of dihydroazulene (DHA) is investigated by nonadiabatic molecular dynamics simulations in connection with the mixed-reference spin-flip (MRSF)-TDDFT method. It is found that, in the main reaction channel, the ring opening occurs in the excited state in a sequence of steps with increasing aromaticity. The first stage lasting ca. 200 fs produces an 8π semiaromatic S1 minimum (S1, min) through an ultrafast damped bond length alternation (BLA) movement synchronized with a partial planarization of the cycloheptatriene ring. An additional ca. 200 fs are required to gain the vibrational energy needed to overcome a ring-opening transition state characterized by an enhanced Baird aromaticity. Unlike other BLA motions of ππ* state, it was shown that their damping is a characteristic feature of aromatic bond-equalization process. In addition, some minor channels of the reaction have also been discovered, where noticeably higher barriers of the S1 non/antiaromatic transition structures must be surmounted. These anti-Baird channels led to reformation of DHA or other closed-ring products. The observed competition between the Baird and anti-Baird channels suggests that the quantum yield of photochemical products can be controllable by tipping their balance. Hence, here we suggest including the concept of anti-Baird, which would expand the applicability of Baird rule to much broader situations.
Collapse
Affiliation(s)
- Svetlana Shostak
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Juwon Oh
- Department of ICT Environmental Health System (Graduate school) and Department of Chemistry, Soonchunhyang University, Asan, Chungnam 31538, Republic of Korea
| | - Jinseok Kim
- Spectroscopy Laboratory for Functional π-Electronic Systems and Department of Chemistry, Yonsei University, Seoul 03722, Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Hyeongwoo Nam
- Department of Chemistry, Yonsei University, Seoul 03722, Korea
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Dongho Kim
- Spectroscopy Laboratory for Functional π-Electronic Systems and Department of Chemistry, Yonsei University, Seoul 03722, Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
11
|
Huix-Rotllant M, Schwinn K, Pomogaev V, Farmani M, Ferré N, Lee S, Choi CH. Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT. J Chem Theory Comput 2023; 19:147-156. [PMID: 36574493 DOI: 10.1021/acs.jctc.2c01010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fittingfitted (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) nonadiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed accelerated deactivation decay pathways, which is consistent with the experimental decay time of ∼400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.
Collapse
Affiliation(s)
| | - Karno Schwinn
- Aix-Marseille Univ, CNRS, ICR, Marseille13013, France
| | - Vladimir Pomogaev
- Department of Chemistry, Kyungpook National University, Daegu41566, South Korea
| | - Maryam Farmani
- Department of Chemistry, Kyungpook National University, Daegu41566, South Korea
| | - Nicolas Ferré
- Aix-Marseille Univ, CNRS, ICR, Marseille13013, France
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California91125, United States
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu41566, South Korea
| |
Collapse
|
12
|
Park W, Alías-Rodríguez M, Cho D, Lee S, Huix-Rotllant M, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory for Accurate X-ray Absorption Spectroscopy. J Chem Theory Comput 2022; 18:6240-6250. [PMID: 36166346 DOI: 10.1021/acs.jctc.2c00746] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
It is demonstrated that the challenging core-hole particle (CHP) orbital relaxation for core electron spectra can be readily achieved by the mixed-reference spin-flip (MRSF)-time-dependent density functional theory (TDDFT). With the additional scalar relativistic effects on K-edge excitation energies of 24 second- and 17 third-row molecules, the particular ΔCHP-MRSF(R) exhibited near perfect predictions with RMSE ∼0.5 eV, featuring a median value of 0.3 and an interquartile range of 0.4. Overall, the CHP effect is 2-4 times stronger than relativistic ones, contributing more than 20 eV in the cases of sulfur and chlorine third-row atoms. Such high precision allows to explain the splitting and spectral shapes of O, N, and C atom K-edges in the ground state of thymine with atom as well as orbital specific accuracy. The same protocol with a double hole particle relaxation also produced remarkably accurate K-edge spectra of core to valence hole excitation energies from the first (nO8π*) and second (ππ*) excited states of thymine, confirming the assignment of 1s → n excitation for the experimentally observed 526.4 eV peak. Regarding both accuracy and practicality, therefore, MRSF-TDDFT provides a promising protocol for core electron spectra of both ground and excited electronic states alike.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Marc Alías-Rodríguez
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, Marseille 13284, France
| | - Daeheum Cho
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Miquel Huix-Rotllant
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, Marseille 13284, France
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
13
|
Park W, Filatov (Gulak) M, Sadiq S, Gerasimov I, Lee S, Joo T, Choi CH. A Plausible Mechanism of Uracil Photohydration Involves an Unusual Intermediate. J Phys Chem Lett 2022; 13:7072-7080. [PMID: 35900137 PMCID: PMC9358713 DOI: 10.1021/acs.jpclett.2c01694] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2022] [Accepted: 07/11/2022] [Indexed: 05/28/2023]
Abstract
It is well-known that photolysis of pyrimidine nucleobases, such as uracil, in an aqueous environment results in the formation of hydrate as one of the main products. Although several hypotheses regarding photohydration have been proposed in the past, e.g., the zwitterionic and "hot" ground-state mechanisms, its detailed mechanism remains elusive. Here, theoretical nonadiabatic simulations of the uracil photodynamics reveal the formation of a highly energetic but kinetically stable intermediate that features a half-chair puckered pyrimidine ring and a strongly twisted intracyclic double bond. The existence and the kinetic stability of the intermediate are confirmed by a variety of computational chemistry methods. According to the simulations, the unusual intermediate is mainly formed almost immediately (∼50-200 fs) upon photoabsorption and survives long enough to engage in a hydration reaction with a neighboring water. A plausible mechanism of uracil photohydration is proposed on the basis of the modeling of nucleophilic insertion of water into the twisted double bond of the intermediate.
Collapse
Affiliation(s)
- Woojin Park
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | | | - Saima Sadiq
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Igor Gerasimov
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Taiha Joo
- Department
of Chemistry, Pohang University of Science
and Technology (POSTECH), Pohang 37673, South Korea
| | - Cheol Ho Choi
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
14
|
David G, Ben Amor N, Zeng T, Suaud N, Trinquier G, Malrieu JP. Difficulty of the evaluation of the barrier height of an open-shell transition state between closed shell minima: The case of small C 4n rings. J Chem Phys 2022; 156:224104. [PMID: 35705394 DOI: 10.1063/5.0090129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
C4n cyclacenes exhibit strong bond-alternation in their equilibrium geometry. In the two equivalent geometries, the system keeps an essentially closed-shell character. The two energy minima are separated by a transition state suppressing the bond-alternation, where the wave function is strongly diradical. This paper discusses the physical factors involved in this energy difference and possible evaluations of the barrier height. The barrier given as the energy difference between the restricted density functional theory (DFT)/B3LYP for the equilibrium and the broken symmetry DFT/B3LYP of the transition state is either negative or small, in contradiction with the most reliable Wave Function Theory calculations. The minimal (two electrons in two molecular orbitals) Complete Active Space self-consistent field (CASSCF) overestimates the barrier, and the subsequent second-order perturbation cancels it. Due to the collective character of the spin-polarization effect, it is necessary to perform a full π CASSCF + second-order perturbation to reach a reasonable value of the barrier, but this type of treatment cannot be applied to large molecules. DFT procedures treating on an equal foot the closed-shell and open-shell geometries have been explored, such as Mixed-Reference Spin-Flip Time-dependent-DFT and a new spin-decontamination proposal, namely, DFT-dressed configuration interaction, but the results still depend on the density functional. M06-2X without or with spin-decontamination gives the best agreement with the accurate wave function results.
Collapse
Affiliation(s)
- Grégoire David
- University Rennes, CNRS, ISCR (Institut des Sciences Chimiques de Rennes)-UMR 6226, F-35000 Rennes, France
| | - Nadia Ben Amor
- Laboratoire de Chimie et Physique Quantiques, IRSAMC-CNRS-UMR5626, Université Paul-Sabatier (Toulouse III), 31062 Toulouse Cedex 4, France
| | - Tao Zeng
- Department of Chemistry, York University, Toronto, Ontario M3J1P3, Canada
| | - Nicolas Suaud
- Laboratoire de Chimie et Physique Quantiques, IRSAMC-CNRS-UMR5626, Université Paul-Sabatier (Toulouse III), 31062 Toulouse Cedex 4, France
| | - Georges Trinquier
- Laboratoire de Chimie et Physique Quantiques, IRSAMC-CNRS-UMR5626, Université Paul-Sabatier (Toulouse III), 31062 Toulouse Cedex 4, France
| | - Jean-Paul Malrieu
- Laboratoire de Chimie et Physique Quantiques, IRSAMC-CNRS-UMR5626, Université Paul-Sabatier (Toulouse III), 31062 Toulouse Cedex 4, France
| |
Collapse
|
15
|
Liang W, Pei Z, Mao Y, Shao Y. Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges. J Chem Phys 2022; 156:210901. [DOI: 10.1063/5.0088271] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin–orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm–Dancoff approximation to TDDFT, spin–flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine–Ag20 complex), we emphasize that caution is needed to properly account for system–environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.
Collapse
Affiliation(s)
- WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Zheng Pei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| |
Collapse
|
16
|
Nam JS, Hong Y, Lee CG, Kim TI, Lee C, Roh DH, Lee IS, Kweon S, Ahn G, Min SK, Kim BS, Kwon TH. Singlet Oxygen Generation from Polyaminoglycerol by Spin-Flip-Based Electron Transfer. JACS AU 2022; 2:933-942. [PMID: 35557761 PMCID: PMC9088781 DOI: 10.1021/jacsau.2c00050] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 03/15/2022] [Accepted: 03/15/2022] [Indexed: 06/15/2023]
Abstract
Reactive oxygen species have drawn attention owing to their strong oxidation ability. In particular, the singlet oxygen (1O2) produced by energy transfer is the predominant species for controlling oxidation reactions efficiently. However, conventional 1O2 generators, which rely on enhanced energy transfer, frequently suffer from poor solubility, low stability, and low biocompatibility. Herein, we introduce a hyperbranched aliphatic polyaminoglycerol (hPAG) as a 1O2 generator, which relies on spin-flip-based electron transfer. The coexistence of a lone pair electron on the nitrogen atom and a hydrogen-bonding donor (the protonated form of nitrogen and hydroxyl group) affords proximity between hPAG and O2. Subsequent direct electron transfer after photo-irradiation induces hPAG•+-O2 •- formation, and the following spin-flip process generates 1O2. The spin-flip-based electron transfer pathway is analyzed by a series of photophysical, electrochemical, and computational studies. The 1O2 generator, hPAG, is successfully employed in photodynamic therapy and as an antimicrobial reagent.
Collapse
Affiliation(s)
- Jung Seung Nam
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
- Center
for Wave Energy Materials, Ulsan National
Institute of Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Youngjoo Hong
- Department
of Chemistry, Yonsei University, Seoul 03722, Republic of Korea
| | - Chae Gyu Lee
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
- Center
for Wave Energy Materials, Ulsan National
Institute of Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Tae In Kim
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Chaiheon Lee
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
- Center
for Wave Energy Materials, Ulsan National
Institute of Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Deok-Ho Roh
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
- Center
for Wave Energy Materials, Ulsan National
Institute of Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - In Seong Lee
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Songa Kweon
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Gyunhyeok Ahn
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Seung Kyu Min
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| | - Byeong-Su Kim
- Department
of Chemistry, Yonsei University, Seoul 03722, Republic of Korea
| | - Tae-Hyuk Kwon
- Department
of Chemistry, Ulsan National Institute of
Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
- Center
for Wave Energy Materials, Ulsan National
Institute of Science and Technology (UNIST), Ulsan 44919, Republic
of Korea
| |
Collapse
|
17
|
James D, Pradhan E, Lee S, Choi CH, Zeng T. Dicarbonyl anthracenes and phenanthrenes as singlet fission chromophores. CAN J CHEM 2022. [DOI: 10.1139/cjc-2021-0241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Singlet fission is a highly desired process in photovoltaic devices as it can significantly enhance photoelectric conversion efficiency. Exploitation of this process in photovoltaics is hindered by the lack of appropriate chromophores. We used mixed-reference spin-flipping time-dependent density functional theory to investigate five di-carbonyl anthracenes and phenanthrenes, with the purpose to design singlet fission chromophores. Two molecules were found to be promising candidates. For all the dicarbonyl molecules, the oxygen lone pair orbitals were found to be involved in the excited states that are relevant to singlet fission.
Collapse
Affiliation(s)
- Dylan James
- York University, 7991, Department of Chemistry, Toronto, Ontario, Canada
| | - Ekadashi Pradhan
- York University, 7991, Department of Chemistry, Toronto, Ontario, Canada
| | - Seunghoon Lee
- California Institute of Technology, 6469, Division of Chemistry and Chemical Engineering, Pasadena, California, United States
| | - Cheol Ho Choi
- Kyungpook National University, 34986, Department of Chemistry, Daegu, Daegu, Korea (the Republic of)
| | - Tao Zeng
- York University, 7991, Department of Chemistry, Toronto, Canada, M3J 1P3
| |
Collapse
|
18
|
Lee S, Park W, Nakata H, Filatov M, Choi CH. Recent advances in ensemble density functional theory and linear response theory for strong correlation. B KOREAN CHEM SOC 2021. [DOI: 10.1002/bkcs.12429] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Seunghoon Lee
- Division of Chemistry and Chemical Engineering California Institute of Technology Pasadena California USA
| | - Woojin Park
- Department of Chemistry Kyungpook National University Daegu South Korea
| | - Hiroya Nakata
- Department of Chemistry Kyungpook National University Daegu South Korea
| | - Michael Filatov
- Department of Chemistry Kyungpook National University Daegu South Korea
| | - Cheol Ho Choi
- Department of Chemistry Kyungpook National University Daegu South Korea
| |
Collapse
|
19
|
Pomogaev V, Lee S, Shaik S, Filatov M, Choi CH. Exploring Dyson's Orbitals and Their Electron Binding Energies for Conceptualizing Excited States from Response Methodology. J Phys Chem Lett 2021; 12:9963-9972. [PMID: 34617764 DOI: 10.1021/acs.jpclett.1c02494] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The molecular orbital (MO) concept is a useful tool, which relates the molecular ground-state energy with the energies (and occupations) of the individual orbitals. However, analysis of the excited states from linear response computations is performed in terms of the initial state MOs or some other forms of orbitals, e.g., natural or natural transition orbitals. Because these orbitals lack the respective energies, they do not allow developing a consistent orbital picture of the excited states. Herein, we argue that Dyson's orbitals enable description of the response states compatible with the concepts of molecular orbital theory. The Dyson orbitals and their energies obtained by mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) for the response ground state are remarkably similar to the canonical MOs obtained by the usual DFT calculation. For excited states, the Dyson orbitals provide a chemically sensible picture of the electronic transitions, thus bridging the chasm between orbital theory and response computations.
Collapse
Affiliation(s)
- Vladimir Pomogaev
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Sason Shaik
- The Lise Meitner-Minerva Center for Computational Quantum Chemistry, Institute of Chemistry, The Hebrew University of Jerusalem, Givat Ram Campus, Jerusalem 91904, Israel
| | - 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
|
20
|
Park W, Shen J, Lee S, Piecuch P, Filatov M, Choi CH. Internal Conversion between Bright (1 1Bu+) and Dark (2 1Ag-) States in s- trans-Butadiene and s- trans-Hexatriene. J Phys Chem Lett 2021; 12:9720-9729. [PMID: 34590847 DOI: 10.1021/acs.jpclett.1c02707] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Internal conversion (IC) between the two lowest singlet excited states, 11Bu+ and 21Ag-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){No}, and DEA-EOMCC(4p2h){Nu} or DEA-EOMCC(3p1h,4p2h){Nu}, consistently predict that the 11Bu+/21Ag- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 11Bu+ and 21Ag- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
21
|
Zhang X, Herbert JM. Nonadiabatic dynamics with spin-flip vs linear-response time-dependent density functional theory: A case study for the protonated Schiff base C 5H 6NH 2. J Chem Phys 2021; 155:124111. [PMID: 34598550 DOI: 10.1063/5.0062757] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
Nonadiabatic trajectory surface hopping simulations are reported for trans-C5H6NH2 +, a model of the rhodopsin chromophore, using the augmented fewest-switches algorithm. Electronic structure calculations were performed using time-dependent density functional theory (TDDFT) in both its conventional linear-response (LR) and its spin-flip (SF) formulations. In the SF-TDDFT case, spin contamination in the low-lying singlet states is removed by projecting out the lowest triplet component during iterative solution of the TDDFT eigenvalue problem. The results show that SF-TDDFT qualitatively describes the photoisomerization of trans-C5H6NH2 +, with favorable comparison to previous studies using multireference electronic structure methods. In contrast, conventional LR-TDDFT affords qualitatively different photodynamics due to an incorrect excited-state potential surface near the Franck-Condon region. In addition, the photochemistry (involving pre-twisting of the central double bond) appears to be different for SF- and LR-TDDFT, which may be a consequence of different conical intersection topographies afforded by these two methods. The present results contrast with previous surface-hopping studies suggesting that the LR-TDDFT method's incorrect topology around S1/S0 conical intersections is immaterial to the photodynamics.
Collapse
Affiliation(s)
- Xing Zhang
- Department of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
22
|
Relief of excited-state antiaromaticity enables the smallest red emitter. Nat Commun 2021; 12:5409. [PMID: 34518551 PMCID: PMC8438045 DOI: 10.1038/s41467-021-25677-2] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Accepted: 08/18/2021] [Indexed: 11/08/2022] Open
Abstract
It is commonly accepted that a large π-conjugated system is necessary to realize low-energy electronic transitions. Contrary to this prevailing notion, we present a new class of light-emitters utilizing a simple benzene core. Among different isomeric forms of diacetylphenylenediamine (DAPA), o- and p-DAPA are fluorescent, whereas m-DAPA is not. Remarkably, p-DAPA is the lightest (FW = 192) molecule displaying red emission. A systematic modification of the DAPA system allows the construction of a library of emitters covering the entire visible color spectrum. Theoretical analysis shows that their large Stokes shifts originate from the relief of excited-state antiaromaticity, rather than the typically assumed intramolecular charge transfer or proton transfer. A delicate interplay of the excited-state antiaromaticity and hydrogen bonding defines the photophysics of this new class of single benzene fluorophores. The formulated molecular design rules suggest that an extended π-conjugation is no longer a prerequisite for a long-wavelength light emission.
Collapse
|
23
|
Lee S, Horbatenko Y, Filatov M, Choi CH. Fast and Accurate Computation of Nonadiabatic Coupling Matrix Elements Using the Truncated Leibniz Formula and Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory. J Phys Chem Lett 2021; 12:4722-4728. [PMID: 33983029 DOI: 10.1021/acs.jpclett.1c00932] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10-10 a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H2 for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.
Collapse
Affiliation(s)
- Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Yevhen Horbatenko
- 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
|
24
|
Park W, Lee S, Huix-Rotllant M, Filatov M, Choi CH. Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine. J Phys Chem Lett 2021; 12:4339-4346. [PMID: 33929858 DOI: 10.1021/acs.jpclett.1c00712] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S1-trapping, S2-trapping, and S1&S2-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S2-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S1-trapping mechanism with two lifetimes, τ1 = 30 ± 1 fs and τ2 = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S2 (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S2→S1 internal conversion and resides around the minimum on the S1 (nOπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S2→S1 internal conversion is mediated by fast bond length alternation distortion, the subsequent S1→S0 occurs through several conical intersections, involving a slow puckering motion.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
25
|
Baek YS, Lee S, Filatov M, Choi CH. Optimization of Three State Conical Intersections by Adaptive Penalty Function Algorithm in Connection with the Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory Method (MRSF-TDDFT). J Phys Chem A 2021; 125:1994-2006. [PMID: 33651623 DOI: 10.1021/acs.jpca.0c11294] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A new adaptive algorithm for penalty function optimization for minimum-energy three-states conical intersections (ME3CI) is suggested. The new algorithm differs from the original penalty function algorithm by (a) removing the redundancy in the target function, (b) using an adaptive increment for the penalty function weighting factor, and (c) using tighter convergence criteria for the energy gap. The latter was introduced to guarantee convergence to a true conical intersection rather than to a narrowly avoided crossing geometry. The new algorithm was tested in the optimization of the ME3CI geometries in butadiene and malonaldehyde, where all of the previously found true ME3CI geometries were recovered. The previously found butadiene's CI3/2/1 turned out to be a narrowly avoided crossing. For butadiene, seven new ME3CI geometries have been located. Because of the removal of the redundancy and the use of the adaptive weighting factor, the convergence rate of the new algorithm is noticeably improved as compared to that of the previously proposed penalty function algorithm. The application to malonaldehyde and butadiene demonstrates that the three-state conical intersections may be more abundant and hence more involved in the photochemistry than previously thought. The recently developed mixed-reference spin flip (MRSF)-TDDFT method yields ME3CI geometries and relative energies quantitatively consistent with the previously reported calculations at a much reduced computational cost.
Collapse
Affiliation(s)
- Yong Su Baek
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - 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
|
26
|
Horbatenko Y, Sadiq S, Lee S, Filatov M, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory (MRSF-TDDFT) as a Simple yet Accurate Method for Diradicals and Diradicaloids. J Chem Theory Comput 2021; 17:848-859. [PMID: 33401894 DOI: 10.1021/acs.jctc.0c01074] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Due to their multiconfigurational nature featuring strong electron correlation, accurate description of diradicals and diradicaloids is a challenge for quantum chemical methods. The recently developed mixed-reference spin-flip (MRSF)-TDDFT method is capable of describing the multiconfigurational electronic states of these systems while avoiding the spin-contamination pitfalls of SF-TDDFT. Here, we apply MRSF-TDDFT to study the adiabatic singlet-triplet (ST) gaps in a series of well-known diradicals and diradicaloids. On average, MRSF displays a very high prediction accuracy of the adiabatic ST gaps with the mean absolute error (MAE) amounting to 0.14 eV. In addition, MRSF is capable of accurately describing the effect of the Jahn-Teller distortion occurring in the trimethylenemethane diradical, the violation of the Hund rule in a series of the didehydrotoluene diradicals, and the potential energy surfaces of the didehydrobenzene (benzyne) diradicals. A convenient criterion for distinguishing diradicals and diradicaloids is suggested on the basis of the easily obtainable quantities. In all of these cases, which are difficult for the conventional methods of density functional theory (DFT), MRSF shows results consistent with the experiment and the high-level ab initio computations. Hence, the present study documents the reliability and accuracy of MRSF and lays out the guidelines for its application to strongly correlated molecular systems.
Collapse
Affiliation(s)
- Yevhen Horbatenko
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Saima Sadiq
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - 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
|
27
|
Horbatenko Y, Lee S, Filatov M, Choi CH. How Beneficial Is the Explicit Account of Doubly-Excited Configurations in Linear Response Theory? J Chem Theory Comput 2021; 17:975-984. [DOI: 10.1021/acs.jctc.0c01214] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Yevhen Horbatenko
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - 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
|
28
|
Wang Z, Li Z, Zhang Y, Liu W. Analytic energy gradients of spin-adapted open-shell time-dependent density functional theory. J Chem Phys 2020; 153:164109. [DOI: 10.1063/5.0025428] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Zikuan Wang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People’s Republic of China
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Yong Zhang
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| |
Collapse
|
29
|
Pradhan E, Lee S, Choi CH, Zeng T. Diboron- and Diaza-Doped Anthracenes and Phenanthrenes: Their Electronic Structures for Being Singlet Fission Chromophores. J Phys Chem A 2020; 124:8159-8172. [PMID: 32902270 DOI: 10.1021/acs.jpca.0c06915] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We used quantum chemistry methods at the levels of mixed-reference spin-flip time-dependent density functional theory and multireference perturbation theory to study diboron- and diaza-doped anthracenes and phenanthrenes. This class of structures recently surged as potential singlet fission chromophores. We studied electronic structures of their excited states and clarified the reasons why they satisfy or fail to satisfy the energy criteria for singlet fission chromophores. Many studied structures have their S1 states not dominated by HOMO → LUMO excitation, so they cannot be described using the conventional two site model. This is attributed to frontier orbital energy shifts induced by the doping and different charge-transfer energies in different one-electron singlet excitations or, in other words, different polarizations of hole and/or particle orbitals in their S1 and T1 states. There is a mirror relation between the orbital energy shifts induced by diboron- and diaza-dopings, which together with alternant hydrocarbon pairings of occupied and unoccupied orbitals, leads to more mirror relations between the excited states of the two types of doped structures.
Collapse
Affiliation(s)
- Ekadashi Pradhan
- Department of Chemistry, York University, Toronto, Ontario M3J1P3, Canada
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Tao Zeng
- Department of Chemistry, York University, Toronto, Ontario M3J1P3, Canada
| |
Collapse
|
30
|
Pradhan E, Bentley JN, Caputo CB, Zeng T. Designs of Singlet Fission Chromophores with a Diazadiborinine Framework**. CHEMPHOTOCHEM 2020. [DOI: 10.1002/cptc.202000081] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Ekadashi Pradhan
- Department of Chemistry York University Toronto Ontario M3 J 1P3 Canada
| | - Jordan N. Bentley
- Department of Chemistry York University Toronto Ontario M3 J 1P3 Canada
| | | | - Tao Zeng
- Department of Chemistry York University Toronto Ontario M3 J 1P3 Canada
| |
Collapse
|
31
|
Photoinduced Double Proton Transfer in the Glyoxal-Methanol Complex Revisited: The Role of the Excited States. J Chem Theory Comput 2020; 16:3273-3286. [PMID: 32275423 DOI: 10.1021/acs.jctc.0c00007] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Under irradiation in the visible range, the glyoxal-methanol complex in a cryogenic argon matrix undergoes a double proton transfer (DPT) reaction through which the glyoxal molecule isomerizes into hydroxyketene. In this work, we employ electronic structure simulations in order to shed more light on the underlying mechanism. Rewardingly, we find that the lowest singlet excited state (S1) of the complex acts as a gateway to two previously unknown isomerization pathways, of which one takes place entirely in the singlet manifold and the other also involves the lowest triplet state (T1). Both of these pathways are fully compatible with the available experimental data, implying that either or both are operative under experimental conditions. In either pathway, the methanol molecule acts as a proton shuttle between the proton-donating and proton-accepting sites of glyoxal, resulting in a dramatic lowering of the potential energy barrier to isomerization with respect to the case of isolated glyoxal. The occurrence of DPT in the singlet manifold is demonstrated directly with the use of nonadiabatic molecular dynamics simulations at the spin-flip time-dependent density functional theory level.
Collapse
|
32
|
Minezawa N, Nakajima T. Quantum mechanical/molecular mechanical trajectory surface hopping molecular dynamics simulation by spin-flip time-dependent density functional theory. J Chem Phys 2020; 152:024119. [DOI: 10.1063/1.5132879] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Noriyuki Minezawa
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Takahito Nakajima
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
33
|
Horbatenko Y, Lee S, Filatov M, Choi CH. Performance Analysis and Optimization of Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory (MRSF-TDDFT) for Vertical Excitation Energies and Singlet–Triplet Energy Gaps. J Phys Chem A 2019; 123:7991-8000. [DOI: 10.1021/acs.jpca.9b07556] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yevhen Horbatenko
- Department of Chemistry, Kyungpook National University, Daegu 41566, Republic of Korea
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul 151747, Republic of Korea
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, Republic of Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, Republic of Korea
| |
Collapse
|
34
|
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
|
35
|
Japahuge A, Lee S, Choi CH, Zeng T. Design of singlet fission chromophores with cyclic (alkyl)(amino) carbene building blocks. J Chem Phys 2019; 150:234306. [PMID: 31228896 DOI: 10.1063/1.5099062] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
We use MRSF-TDDFT and NEVPT2 methods to design singlet fission chromophores with the building blocks of cyclic (alkyl)(amino)carbenes (CAACs). CAAC dimers with C2, C4, and p-phenylene spacers are considered. The substitutions with trifluoromethyls and fluorine atoms at the α C position are investigated. The electronegative substituents enhance the π accepting capability of the α C while maintaining it as a quaternary C atom. The phenylene-connected dimers with the two substitutions are identified as promising candidates for singlet fission chromophores. The cylindrically symmetric C2 and C4 spacers allow for substantial structural reorganizations in the S0-to-S1 and S0-to-T1 excitations. Although the two substituted dimers with the C4 spacer satisfy (or very close to satisfy) the primary thermodynamics criterion for singlet fission, the significant structural reorganizations result in high barriers so that the fission is kinetically unfavorable.
Collapse
Affiliation(s)
- Achini Japahuge
- Department of Chemistry, Carleton University, Ottawa, Ontario K1S5B6, Canada
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul 151-747, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Tao Zeng
- Department of Chemistry, Carleton University, Ottawa, Ontario K1S5B6, Canada
| |
Collapse
|