1
|
Myers CA, Miyazaki K, Trepl T, Isborn CM, Ananth N. GPU-accelerated on-the-fly nonadiabatic semiclassical dynamics. J Chem Phys 2024; 161:084114. [PMID: 39193942 DOI: 10.1063/5.0223628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2024] [Accepted: 08/11/2024] [Indexed: 08/29/2024] Open
Abstract
GPU-accelerated on-the-fly nonadiabatic dynamics is enabled by interfacing the linearized semiclassical dynamics approach with the TeraChem electronic structure program. We describe the computational workflow of the "PySCES" code interface, a Python code for semiclassical dynamics with on-the-fly electronic structure, including parallelization over multiple GPU nodes. We showcase the abilities of this code and present timings for two benchmark systems: fulvene solvated in acetonitrile and a charge transfer system in which a photoexcited zinc-phthalocyanine donor transfers charge to a fullerene acceptor through multiple electronic states on an ultrafast timescale. Our implementation paves the way for an efficient semiclassical approach to model the nonadiabatic excited state dynamics of complex molecules, materials, and condensed phase systems.
Collapse
Affiliation(s)
- Christopher A Myers
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Ken Miyazaki
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Thomas Trepl
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
2
|
Kossoski F, Boggio-Pasqua M, Loos PF, Jacquemin D. Reference Energies for Double Excitations: Improvement and Extension. J Chem Theory Comput 2024; 20:5655-5678. [PMID: 38885174 DOI: 10.1021/acs.jctc.4c00410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.
Collapse
Affiliation(s)
- Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
3
|
Xu J, Hao J, Bu C, Meng Y, Xiao H, Zhang M, Li C. XMECP: Reaching State-of-the-Art MECP Optimization in Multiscale Complex Systems. J Chem Theory Comput 2024; 20:3590-3600. [PMID: 38651739 DOI: 10.1021/acs.jctc.4c00033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
The Python-based program, XMECP, is developed for realizing robust, efficient, and state-of-the-art minimum energy crossing point (MECP) optimization in multiscale complex systems. This article introduces the basic capabilities of the XMECP program by theoretically investigating the MECP mechanism of several example systems including (1) the photosensitization mechanism of benzophenone, (2) photoinduced proton-coupled electron transfer in the cytosine-guanine base pair in DNA, (3) the spin-flip process in oxygen activation catalyzed by an iron-containing 2-oxoglutarate-dependent oxygenase (Fe/2OGX), and (4) the photochemical pathway of flavoprotein adjusted by the intensity of an external electric field. MECPs related to multistate reaction and multistate reactivity in large-scale complex biochemical systems can be well-treated by workflows suggested by the XMECP program. The branching plane updating the MECP optimization algorithm is strongly recommended as it provides derivative coupling vector (DCV) with explicit calculation and can equivalently evaluate contributions from non-QM residues to DCV, which can be nonadiabatic coupling or spin-orbit coupling in different cases. In the discussed QM/MM examples, we also found that the influence on the QM region by DCV can occur through noncovalent interactions and decay with distance. In the example of DNA base pairs, the nonadiabatic coupling occurs across the π-π stacking structure formed in the double-helix system. In contrast to general intuition, in the example of Fe/2OGX, the central ferrous and oxygen part contribute little to the spin-orbit coupling; however, a nearby arginine residue, which is treated by molecular mechanics in the QM/MM method, contributes significantly via two hydrogen bonds formed with α-ketoglutarate (α-KG). This indicates that the arginine residue plays a significant role in oxygen activation, driving the initial triplet state toward the productive quintet state, which is more than the previous knowledge that the arginine residue can bind α-KG at the reaction site by hydrogen bonds.
Collapse
Affiliation(s)
- Jiawei Xu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Jian Hao
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Caijie Bu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- College of Chemistry and Materials Science, Fujian Normal University, Fuzhou 350117, Fujian, P. R. China
| | - Yajie Meng
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Han Xiao
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
| | - Minyi Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
| | - Chunsen Li
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen 361005, Fujian, P. R. China
| |
Collapse
|
4
|
Garcia MB, Singh M, Miller E, Neuenswander S, Douglas J, Boskovic Z. Twisted Intramolecular Charge-Transfer State Addition to Electron-Poor Olefins. J Org Chem 2024; 89:3058-3064. [PMID: 38354334 PMCID: PMC11006016 DOI: 10.1021/acs.joc.3c02521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2024]
Abstract
When electron-rich arylpyrrolinium salts are irradiated with ultraviolet light in the presence of Michael acceptors, the pyrrolinyl and aryl fragments add to the activated and polarized double bond in a regioselective manner, forming two C-C bonds and fragmenting the substrate. In this paper, we present a model for this intriguing reaction, supported by spectroscopy and computational analyses, and provide evidence for rectifying previously misassigned structures. We postulate that the photochemical reaction is inefficient because the reaction between the twisted intramolecular charge-transfer state and the olefin competes with fluorescence from this state upon photon absorption. We also discuss the practical advantages of performing this photochemical reaction in a continuous flow setup. Additionally, we explore several subsequent reactions that allow us to further modify the products of the photochemical step, ultimately leading to the creation of new chemical structures.
Collapse
Affiliation(s)
- Mauricio Bahena Garcia
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Manvendra Singh
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Elizabeth Miller
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Sarah Neuenswander
- Nuclear Magnetic Resonance Laboratory, University of Kansas, Lawrence, Kansas 66045, United States
| | - Justin Douglas
- Nuclear Magnetic Resonance Laboratory, University of Kansas, Lawrence, Kansas 66045, United States
| | - Zarko Boskovic
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
5
|
Ortiz-Rodríguez LA, Caldero-Rodríguez NE, Seth SK, Díaz-González K, Crespo-Hernández CE. Electronic relaxation mechanism of 9-methyl-2,6-diaminopurine and 2,6-diaminopurine-2'-deoxyribose in solution. Photochem Photobiol 2024; 100:393-403. [PMID: 38018292 DOI: 10.1111/php.13887] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 10/16/2023] [Accepted: 11/05/2023] [Indexed: 11/30/2023]
Abstract
Prolonged ultraviolet exposure results in the formation of cyclobutane pyrimidine dimers (CPDs) in RNA. Consequently, prebiotic photolesion repair mechanisms should have played an important role in the maintenance of the structural integrity of primitive nucleic acids. 2,6-Diaminopurine is a prebiotic nucleobase that repairs CPDs with high efficiency when incorporated into polymers. We investigate the electronic deactivation pathways of 2,6-diaminopurine-2'-deoxyribose and 9-methyl-2,6-diaminopurine in acetonitrile and aqueous solution to shed light on the photophysical and excited state properties of the 2,6-diaminopurine chromophore. Evidence is presented that both are photostable compounds exhibiting similar deactivation mechanisms upon the population of the S1 (ππ* La ) state at 290 nm. The mechanism involves deactivation through the C2- and C6-reaction coordinates and >99% of the excited state population decays through nonradiative pathways involving two conical intersections with the ground state. The radiative and nonradiative lifetimes are longer in aqueous solution compared to acetonitrile. While τ1 is similar in both derivatives, τ2 is ca. 1.5-fold longer in 2,6-diaminopurine-2'-deoxyribose due to a more efficient trapping in the S1 (ππ* La ) minimum. Therefore, 2,6-diaminopurine could have accumulated in significant quantities during prebiotic times to be incorporated into non-canonical RNA and play a significant role in its photoprotection.
Collapse
Affiliation(s)
| | | | - Sourav Kanti Seth
- Department of Chemistry, Case Western Reserve University, Cleveland, Ohio, USA
| | | | | |
Collapse
|
6
|
Fan J, Yue L, Liu C, Rao B, Zhou G, Li A, Su B. Isolation of Fluorescent 2π-Aromatic 1,3-Disiladiboretenes. J Am Chem Soc 2024; 146:39-44. [PMID: 38117532 DOI: 10.1021/jacs.3c09747] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
Herein, we reported the isolation of 2π-aromatic disiladiboretenes (L2Si2B2Ph2) [L = ArC(NtBu)2, Ar = Ph (1), Mes (2)], which have been synthesized from the straightforward reduction of silylene-borane adducts (LSiX → BX2Ph) [X = Cl, Br] with potassium graphite (KC8). X-ray diffraction analysis of 1 and 2 revealed that the Si2B2 units are completely planar, and DFT calculations suggested delocalization of 2π-electrons over the Si2B2 rings. Moreover, their photophysical properties and reactivity toward sulfur were also investigated in detail.
Collapse
Affiliation(s)
- Jiawei Fan
- School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices, Ministry of Education, Xi'an Key Laboratory of Sustainable Energy Materials Chemistry, Xi'an Jiaotong University, Xi'an 710049, China
| | - Ling Yue
- School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices, Ministry of Education, Xi'an Key Laboratory of Sustainable Energy Materials Chemistry, Xi'an Jiaotong University, Xi'an 710049, China
| | - Ce Liu
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, China
| | - Bin Rao
- School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices, Ministry of Education, Xi'an Key Laboratory of Sustainable Energy Materials Chemistry, Xi'an Jiaotong University, Xi'an 710049, China
| | - Guijiang Zhou
- School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices, Ministry of Education, Xi'an Key Laboratory of Sustainable Energy Materials Chemistry, Xi'an Jiaotong University, Xi'an 710049, China
| | - Anyang Li
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, China
| | - Bochao Su
- School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices, Ministry of Education, Xi'an Key Laboratory of Sustainable Energy Materials Chemistry, Xi'an Jiaotong University, Xi'an 710049, China
- State Key Laboratory of Elemento-Organic Chemistry, Nankai University, Tianjin 300071, China
| |
Collapse
|
7
|
Wang PY, Hsu YC, Chen PH, Chen GY, Liao YK, Cheng PY. Solvent-polarity dependence of ultrafast excited-state dynamics of trans-4-nitrostilbene. Phys Chem Chem Phys 2024; 26:788-807. [PMID: 38088777 DOI: 10.1039/d3cp05245a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
Ultrafast excited-state dynamics of the simplest nitrostilbenes, namely trans-4-nitrostilbene (t-NSB), was studied in solvents of various polarities with ultrafast broadband time-resolved fluorescence and transient absorption spectroscopies, and by quantum-chemical computations. The results revealed that the initially excited S1(ππ*) state deactivation dynamics is strongly influenced by the solvent polarity. Specifically, the t-NSB S1-state lifetime decreases by three orders of magnitude from ∼60 ps in high-polarity solvents to ∼60 fs in nonpolar solvents. The strong solvent-polarity dependence arises from the differences in dipole moments among the S1 and relevant states, including the major intersystem crossing (ISC) receiver triplet states, and therefore, the solvent polarity can modulate their relative energies and ISC rates. In nonpolar solvents, the sub-100 fs lifetime is due to a combination of efficient ISC and internal conversion. In medium-polarity solvents, the S1-state population decays via a competing ISC relaxation mechanism in a biphasic manner, and the ISC rates are found to obey the inverse energy gap law of the strong coupling case. In high-polarity solvents, the S1 state is stabilized to a much lower energy such that ISC becomes energetically infeasible, and the S1 state decays via barrier crossing along the torsion angle of the central ethylenic bond to the nonfluorescent perpendicular configuration. Regardless of the initial S1-state deactivation pathways in various solvents, the excited-state population is ultimately trapped in the metastable T1-state perpendicular configuration, at which a slower ISC occurs to bring the system to the ground state and bifurcate into either trans or cis form of NSB.
Collapse
Affiliation(s)
- Peng-Yun Wang
- Department of Chemistry, National Tsing Hua University, Hsinchu, Taiwan, 30043, Republic of China.
| | - Yu-Cheng Hsu
- Department of Chemistry, National Tsing Hua University, Hsinchu, Taiwan, 30043, Republic of China.
| | - Pin-Hsun Chen
- Department of Chemistry, National Tsing Hua University, Hsinchu, Taiwan, 30043, Republic of China.
| | - Guan-Yu Chen
- Department of Chemistry, National Tsing Hua University, Hsinchu, Taiwan, 30043, Republic of China.
| | - Yi-Kai Liao
- Department of Chemistry, National Tsing Hua University, Hsinchu, Taiwan, 30043, Republic of China.
| | - Po-Yuan Cheng
- Department of Chemistry, National Tsing Hua University, Hsinchu, Taiwan, 30043, Republic of China.
| |
Collapse
|
8
|
Sisodiya DS, Ali SM, Chattopadhyay A. Unexplored Isomerization Pathways of Azobis(benzo-15-crown-5): Computational Studies on a Butterfly Crown Ether. J Phys Chem A 2023; 127:7080-7093. [PMID: 37526572 DOI: 10.1021/acs.jpca.3c02363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Computational studies on trans → cis and cis → trans isomerizations of photoresponsive azobis(benzo-15-crown-5) have been reported in this work. The photoexcited ππ* state (S2) of the trans isomer relaxes through the planar S2 minimum and the planar S2/S1 conical intersection (both situated around 9 kcal/mol below the vertically excited S2 state) arising along the N═N stretching coordinate. The nπ* state (S1) of this isomer has both planar and rotated (clockwise and anticlockwise) minima, which may lead to a torsional conical intersection (S0/S1) geometry having a
Collapse
Affiliation(s)
- Dilawar Singh Sisodiya
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani, K.K. Birla Goa Campus, Zuarinagar 403726, India
| | - Sk Musharaf Ali
- Chemical Engineering Division, Bhabha Atomic Research Centre (BARC), Mumbai 400085, India
| | - Anjan Chattopadhyay
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani, K.K. Birla Goa Campus, Zuarinagar 403726, India
| |
Collapse
|
9
|
Intermolecular-Type Conical Intersections in Benzene Dimer. Int J Mol Sci 2023; 24:ijms24032906. [PMID: 36769227 PMCID: PMC9917476 DOI: 10.3390/ijms24032906] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2023] [Revised: 01/27/2023] [Accepted: 01/30/2023] [Indexed: 02/05/2023] Open
Abstract
The equilibrium and conical intersection geometries of the benzene dimer were computed in the framework of the conventional, linear-response time-dependent and spin-flipped time-dependent density functional theories (known as DFT, TDDFT and SF-TDDFT) as well as using the multiconfigurational complete active space self-consistent field (CASSCF) method considering the minimally augmented def2-TZVPP and the 6-31G(d,p) basis sets. It was found that the stacking distance between the benzene monomers decreases by about 0.5 Å in the first electronic excited state, due to the stronger intermolecular interaction energy, bringing the two monomers closer together. Intermolecular-type conical intersection (CI) geometries can be formed between the two benzene molecules, when (i) both monomer rings show planar deformation and (ii) weaker (approximately 1.6-1.8 Å long) C-C bonds are formed between the two monomers, with parallel and antiparallel orientation with respect to the monomer. These intermolecular-type CIs look energetically more favorable than dimeric CIs containing only one deformed monomer. The validity of the dimer-type CI geometries obtained by SF-TDDFT was confirmed by the CASSCF method. The nudged elastic band method used for finding the optimal relaxation path has confirmed both the accessibility of these intermolecular-type CIs and the possibility of the radiationless deactivation of the electronic excited states through these CI geometries. Although not as energetically favorable as the previous two CI geometries, there are other CI geometries characterized by the relative rotation of monomers at different angles around a vertical C-C axis.
Collapse
|
10
|
Kotaru S, Kähler S, Alessio M, Krylov AI. Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin-flip DFT and Heisenberg effective Hamiltonians. J Comput Chem 2023; 44:367-380. [PMID: 35699152 PMCID: PMC10084445 DOI: 10.1002/jcc.26941] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 04/07/2022] [Accepted: 04/22/2022] [Indexed: 12/31/2022]
Abstract
Low-energy spectra of single-molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For 12 binuclear complexes with J-values ranging from -6 to -132 cm-1 , our benchmark calculations using the short-range hybrid ωPBEh functional and 6-31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm-1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm-1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.
Collapse
Affiliation(s)
- Saikiran Kotaru
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| | - Sven Kähler
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| | - Maristella Alessio
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| | - Anna I. Krylov
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| |
Collapse
|
11
|
Ziani Z, Cobo S, Loiseau F, Jouvenot D, Lognon E, Boggio-Pasqua M, Royal G. All Visible Light Photoswitch Based on the Dimethyldihydropyrene Unit Operating in Aqueous Solutions with High Quantum Yields. JACS AU 2023; 3:131-142. [PMID: 36711101 PMCID: PMC9875246 DOI: 10.1021/jacsau.2c00552] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 11/15/2022] [Accepted: 11/15/2022] [Indexed: 06/18/2023]
Abstract
Molecular systems and devices whose properties can be modulated using light as an external stimulus are the subject of numerous research studies in the fields of materials and life sciences. In this context, the use of photochromic compounds that reversibly switch upon light irradiation is particularly attractive. However, for many envisioned applications, and in particular for biological purposes, illumination with harmful UV light must be avoided and these photoactivable systems must operate in aqueous media. In this context, we have designed a benzo[e]-fused dimethyldihydropyrene compound bearing a methyl-pyridinium electroacceptor group that meets these requirements. This compound (closed state) is able to reversibly isomerize under aerobic conditions into its corresponding cyclophanediene form (open isomer) through the opening of its central carbon-carbon bond. Both the photo-opening and the reverse photoclosing processes are triggered by visible light illumination and proceed with high quantum yields (respectively 14.5% yield at λ = 680 nm and quantitative quantum yield at λ = 470 nm, in water). This system has been investigated by nuclear magnetic resonance and absorption spectroscopy, and the efficient photoswitching behavior was rationalized by spin-flip time-dependent density functional theory calculations. In addition, it is demonstrated that the isomerization from the open to the closed form can be electrocatalytically triggered.
Collapse
Affiliation(s)
- Zakaria Ziani
- Univ.
Grenoble Alpes, CNRS, DCM, Grenoble38000, France
| | - Saioa Cobo
- Univ.
Grenoble Alpes, CNRS, DCM, Grenoble38000, France
| | | | | | - Elise Lognon
- LCPQ
UMR 5626, CNRS et Université Toulouse
III − Paul Sabatier, 118 route de Narbonne, Toulouse31062, France
| | - Martial Boggio-Pasqua
- LCPQ
UMR 5626, CNRS et Université Toulouse
III − Paul Sabatier, 118 route de Narbonne, Toulouse31062, France
| | - Guy Royal
- Univ.
Grenoble Alpes, CNRS, DCM, Grenoble38000, France
| |
Collapse
|
12
|
Abstract
Chemiluminescence (CL) utilizing chemiexcitation for energy transformation is one of the most highly sensitive and useful analytical techniques. The chemiexcitation is a chemical process of a ground-state reactant producing an excited-state product, in which a nonadiabatic event is facilitated by conical intersections (CIs), the specific molecular geometries where electronic states are degenerated. Cyclic peroxides, especially 1,2-dioxetane/dioxetanone derivatives, are the iconic chemiluminescent substances. In this Perspective, we concentrated on the CIs in the CL of cyclic peroxides. We first present a computational overview on the role of CIs between the ground (S0) state and the lowest singlet excited (S1) state in the thermolysis of cyclic peroxides. Subsequently, we discuss the role of the S0/S1 CI in the CL efficiency and point out misunderstandings in some theoretical studies on the singlet chemiexcitations of cyclic peroxides. Finally, we address the challenges and future prospects in theoretically calculating S0/S1 CIs and simulating the dynamics and chemiexcitation efficiency in the CL of cyclic peroxides.
Collapse
Affiliation(s)
- Ling Yue
- Key Laboratory for Non-equilibrium Synthesis and Modulation of Condensed Matter, Ministry of Education, School of Chemistry, Xi'an Jiaotong University, Xi'an, Shaanxi710049, China
| | - Ya-Jun Liu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai519087, China
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing100875, China
| |
Collapse
|
13
|
Ni Z, Chen Z, Zhang X, Yang X, Zhou L. Photolysis of Fungicide Triadimefon: A Combined Experimental and Theoretical Investigation on Homolytic C-O and C-N Bonds Dissociation Mechanisms. J Photochem Photobiol A Chem 2022. [DOI: 10.1016/j.jphotochem.2022.114402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
14
|
Salazar E, Reinink S, Faraji S. Providing theoretical insight into the role of symmetry in the photoisomerization mechanism of a non-symmetric dithienylethene photoswitch. Phys Chem Chem Phys 2022; 24:11592-11602. [PMID: 35531648 PMCID: PMC9116444 DOI: 10.1039/d2cp00550f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Dithienylethene (DTE) molecular photoswitches have shown to be excellent candidates in the design of efficient optoelectronic devices, due to their high photoisomerization quantum yield (QY), for which symmetry is suggested to play a crucial role. Here, we present a theoretical study on the photochemistry of a non-symmetric dithienylethene photoswitch, with a special emphasis on the effect of asymmetric substitution on the photocyclization and photoreversion mechanisms. We used the Spin-Flip Time Dependent Density Functional Theory (SF-TDDFT) method to locate and characterize the main structures (conical intersections and minima) of the ground state and the first two excited states, S1 and S2, along the ring-opening/closure reaction coordinate of the photocyclization and photoreversion processes, and to identify the important coordinates governing the radiationless decay pathways. Our results suggest that while the main features that characterize the photoisomerization of symmetric DTEs are also present for the photoisomerization of the non-symmetric DTE, the lower energy barrier on S1 along the cycloreversion reaction speaks in favor of a more efficient and therefore a higher cycloreversion QY for the non-symmetric DTEs, making them a better candidate for molecular optoelectronic devices than their symmetric counterparts.
Collapse
Affiliation(s)
- Edison Salazar
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| | - Suzanne Reinink
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| | - Shirin Faraji
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| |
Collapse
|
15
|
Lee IS, Min SK. Generalized Formulation of the Density Functional Tight Binding-Based Restricted Ensemble Kohn-Sham Method with Onsite Correction to Long-Range Correction. J Chem Theory Comput 2022; 18:3391-3409. [PMID: 35549266 DOI: 10.1021/acs.jctc.2c00037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a generalized formulation for the combination of the density functional tight binding (DFTB) approach and the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method by considering onsite correction (OC) as well as the long-range corrected (LC) functional. The OC contribution provides more accurate energies and analytic gradients for individual microstates, while the multireference character of the SSR provides the correct description for conical intersections. We benchmark the LC-OC-DFTB/SSR method against various DFTB calculation methods for excitation energies and conical intersection structures with π/π* or n/π* characters. Furthermore, we perform excited-state molecular dynamics simulations with a molecular rotary motor with variations of LC-OC-DFTB/SSR approaches. We show that the OC contribution to the LC functional is crucial to obtain the correct geometry of conical intersections.
Collapse
Affiliation(s)
- In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| |
Collapse
|
16
|
On the Photostability of Cyanuric Acid and Its Candidature as a Prebiotic Nucleobase. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27041184. [PMID: 35208973 PMCID: PMC8875432 DOI: 10.3390/molecules27041184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/19/2021] [Revised: 02/02/2022] [Accepted: 02/05/2022] [Indexed: 11/28/2022]
Abstract
Cyanuric acid is a triazine derivative that has been identified from reactions performed under prebiotic conditions and has been proposed as a prospective precursor of ancestral RNA. For cyanuric acid to have played a key role during the prebiotic era, it would have needed to survive the harsh electromagnetic radiation conditions reaching the Earth’s surface during prebiotic times (≥200 nm). Therefore, the photostability of cyanuric acid would have been crucial for its accumulation during the prebiotic era. To evaluate the putative photostability of cyanuric acid in water, in this contribution, we employed density functional theory (DFT) and its time-dependent variant (TD-DFT) including implicit and explicit solvent effects. The calculations predict that cyanuric acid has an absorption maximum at ca. 160 nm (7.73 eV), with the lowest-energy absorption band extending to ca. 200 nm in an aqueous solution and exhibiting negligible absorption at longer wavelengths. Excitation of cyanuric acid at 160 nm or longer wavelengths leads to the population of S5,6 singlet states, which have ππ* character and large oscillator strengths (0.8). The population reaching the S5,6 states is expected to internally convert to the S1,2 states in an ultrafast time scale. The S1,2 states, which have nπ* character, are predicted to access a conical intersection with the ground state in a nearly barrierless fashion (ca. ≤ 0.13 eV), thus efficiently returning the population to the ground state. Furthermore, based on calculated spin–orbit coupling elements of ca. 6 to 8 cm−1, the calculations predict that intersystem crossing to the triplet manifold should play a minor role in the electronic relaxation of cyanuric acid. We have also calculated the vertical ionization energy of cyanuric acid at 8.2 eV, which predicts that direct one-photon ionization of cyanuric acid should occur at ca. 150 nm. Collectively, the quantum-chemical calculations predict that cyanuric acid would have been highly photostable under the solar radiation conditions reaching the Earth’s surface during the prebiotic era in an aqueous solution. Of relevance to the chemical origin of life and RNA-first theories, these observations lend support to the idea that cyanuric acid could have accumulated in large quantities during the prebiotic era and thus strengthens its candidature as a relevant prebiotic nucleobase.
Collapse
|
17
|
Fluorescence of carotenoids: Probing binding site interactions and conformational motion in carotenoproteins. Methods Enzymol 2022; 674:85-111. [DOI: 10.1016/bs.mie.2022.04.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
18
|
Substituents affect the mechanism of photochemical E-Z isomerization of diarylethene triazoles via adiabatic singlet excited state pathway or via triplet excited state. J Photochem Photobiol A Chem 2022. [DOI: 10.1016/j.jphotochem.2021.113567] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
19
|
Zhou Y, Maisonnenuve S, Casimiro L, Retailleau P, Xie J, Maurel F, Métivier R. Photoisomerization of a 4-dicyanomethylene-2-methyl-6-(p-dimethylaminostyryl)-4H-pyran analog dye: a combined photophysical and theoretical investigation. Phys Chem Chem Phys 2022; 24:6282-6289. [DOI: 10.1039/d1cp05170a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A combination of experimental and theoretical investigations of a photoisomerizable analog of 4-dicyanomethylene-2-methyl-6-(p-dimethylaminostyryl)-4H-pyran (DCM) dye molecule is presented. We provide evidences that the 4 main isomers and conformers of DCM...
Collapse
|
20
|
Ziani Z, Loiseau F, Lognon E, Boggio-Pasqua M, Philouze C, Cobo S, Royal G. Synthesis of a Negative Photochrome with High Switching Quantum Yields and Capable of Singlet-Oxygen Production and Storage. Chemistry 2021; 27:16642-16653. [PMID: 34677893 DOI: 10.1002/chem.202103003] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Indexed: 12/31/2022]
Abstract
A dimethyldihydropyrene (DHP) photochromic unit has been functionalized by donor (triphenylamine group) and acceptor (methylpyridinium) substituents. This compound was characterized by NMR, absorption and emission spectroscopies as well as cyclic voltammetry, and its properties were rationalized by theoretical calculations. The incorporation of both electron-donor and -withdrawing groups at the photochromic center allows i) an efficient photo-isomerization of the system when illuminated at low energy (quantum yield: Φc-o =13.3 % at λex =660 nm), ii) the reversible and quantitative formation of two endoperoxyde isomers when illuminated under aerobic conditions at room temperature, and iii) the storage and production of singlet oxygen. The photo-isomerization mechanism was also investigated by spin-flip TD-DFT (SF-TD-DFT) calculations.
Collapse
Affiliation(s)
- Zakaria Ziani
- Département de Chimie Moléculaire, Université Grenoble-Alpes, CNRS, 38400, Grenoble, France
| | - Frédérique Loiseau
- Département de Chimie Moléculaire, Université Grenoble-Alpes, CNRS, 38400, Grenoble, France
| | - Elise Lognon
- LCPQ UMR 5626, CNRS et, Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062, Toulouse, France
| | - Martial Boggio-Pasqua
- LCPQ UMR 5626, CNRS et, Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062, Toulouse, France
| | - Christian Philouze
- Département de Chimie Moléculaire, Université Grenoble-Alpes, CNRS, 38400, Grenoble, France
| | - Saioa Cobo
- Département de Chimie Moléculaire, Université Grenoble-Alpes, CNRS, 38400, Grenoble, France
| | - Guy Royal
- Département de Chimie Moléculaire, Université Grenoble-Alpes, CNRS, 38400, Grenoble, France
| |
Collapse
|
21
|
Guan Y, Xie C, Yarkony DR, Guo H. High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics. Phys Chem Chem Phys 2021; 23:24962-24983. [PMID: 34473156 DOI: 10.1039/d1cp03008f] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China.
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131, USA.
| |
Collapse
|
22
|
Fındık V, Fındık BK, Aviyente V, Monari A. Origins of the photoinitiation capacity of aromatic thiols as photoinitiatiors: a computational study. Phys Chem Chem Phys 2021; 23:24377-24385. [PMID: 34676839 DOI: 10.1039/d1cp04345e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we report the photophysical properties of three thiol derivatives, commonly used as photoinitiators in thiol-ene free radical polymerization, the ultimate goal being to rationalize the main reason behind the photoinitiation efficiency. For this aim, time dependent density functional theory is used to simulate the absorption spectra of alkyl thiol (R-SH), thiophenol (PhSH) and p-(trifluoromethyl) thiophenol (p-CF3PhSH), describe their excited state topologies, and explore their potential energy surfaces along the S-H dissociation. Excited state calculations have shown that the S-H photolysis is achieved through the triplet excited states following intersystem crossing from the originally populated singlet manifolds. More specifically, while in aromatic thiol derivatives dissociation is mainly triplet-state mediated, the first excited singlet state and first triplet state of alkyl thiol are both dissociative and hence potentially capable of generating the photoinduced radical species. We have also justified the experimental findings concerning the photoinitiator efficiency considering both their potential energy surface topologies and the absorption intensity, in the lowest energy region.
Collapse
Affiliation(s)
- Volkan Fındık
- Univesité de Lorraine, CNRS, LPCT, F54000 Nancy, France. .,Department of Chemistry, Faculty of Arts and Sciences, Marmara University, 34722 Istanbul, Turkey
| | - Basak Koca Fındık
- Department of Chemistry, Bogazici University, 34342, Bebek, Istanbul, Turkey
| | - Viktorya Aviyente
- Department of Chemistry, Bogazici University, 34342, Bebek, Istanbul, Turkey
| | - Antonio Monari
- Univesité de Lorraine, CNRS, LPCT, F54000 Nancy, France. .,Université de Paris and CNRS, ITODYS, F-75006, Paris, France.
| |
Collapse
|
23
|
Practical treatment of singlet oxygen with density-functional theory and the multiplet-sum method. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02852-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
24
|
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
|
25
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 61] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
26
|
Photoinduced phenomena in water solution of melamine explaining the photostability of the compound. J Mol Model 2021; 27:196. [PMID: 34104983 DOI: 10.1007/s00894-021-04809-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Accepted: 05/31/2021] [Indexed: 10/21/2022]
Abstract
Two tautomers of melamine (triamino and imino-diamino) were studied at the BLYP/aug-cc-pVDZ theoretical level. It was found that the two tautomers are bridged with the 1πσ* excited-state reaction path. The high photostability of melamine in water solution was explained with the mechanism of ring deformation which occurs along the 1ππ* excited-state reaction path. The two mechanisms are investigated at the TD BLYP level of theory using the linear interpolation in internal coordinates approach.
Collapse
|
27
|
Uratani H, Yoshikawa T, Nakai H. Trajectory Surface Hopping Approach to Condensed-Phase Nonradiative Relaxation Dynamics Using Divide-and-Conquer Spin-Flip Time-Dependent Density-Functional Tight Binding. J Chem Theory Comput 2021; 17:1290-1300. [PMID: 33577323 DOI: 10.1021/acs.jctc.0c01155] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nonradiative relaxation of excited molecules is central to many crucial issues in photochemistry. Condensed phases are typical contexts in which such problems are considered, and the nonradiative relaxation dynamics are expected to be significantly affected by interactions with the environment, for example, a solvent. We developed a nonadiabatic molecular dynamics simulation technique that can treat the nonradiative relaxation and explicitly include the environment in the calculations without a heavy computational burden. Specifically, we combined trajectory surface hopping with Tully's fewest-switches algorithm, a tight-binding approximated version of spin-flip time-dependent density-functional theory, and divide-and-conquer (DC) spatial fragmentation scheme. Numerical results showed that this method can treat systems with thousands of atoms within reasonable computational resources, and the error arising from DC fragmentation is negligibly small. Using this method, we obtained molecular insights into the solvent dependence of the photoexcited-state dynamics of trans-azobenzene, which demonstrate the importance of the environment for condensed-phase nonradiative relaxation.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
28
|
Romeo-Gella F, Corral I, Faraji S. Theoretical investigation of a novel xylene-based light-driven unidirectional molecular motor. J Chem Phys 2021; 154:064111. [PMID: 33588536 DOI: 10.1063/5.0038281] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
In this study, the working mechanism of the first light-driven rotary molecular motors used to control an eight-base-pair DNA hairpin has been investigated. In particular, this linker was reported to have promising photophysical properties under physiological conditions, which motivated our work at the quantum mechanical level. Cis-trans isomerization is triggered by photon absorption at wavelengths ranging 300 nm-400 nm, promoting the rotor to the first excited state, and it is mediated by an energy-accessible conical intersection from which the ground state is reached back. The interconversion between the resulting unstable isomer and its stable form occurs at physiological conditions in the ground state and is thermally activated. Here, we compare three theoretical frameworks, generally used in the quantum description of medium-size chemical systems: Linear-Response Time-Dependent Density Functional Theory (LR-TDDFT), Spin-Flip TDDFT (SF-TDDFT), and multistate complete active space second-order perturbation theory on state-averaged complete active space self consistent field wavefunctions (MS-CASPT2//SA-CASSCF). In particular, we show the importance of resorting to a multireference approach to study the rotational cycle of light-driven molecular motors due to the occurrence of geometries described by several configurations. We also assess the accuracy and computational cost of the SF-TDDFT method when compared to MS-CASPT2 and LR-TDDFT.
Collapse
Affiliation(s)
- F Romeo-Gella
- Departamento de Química (Módulo 13, Facultad de Ciencias) and Institute of Advanced Chemical Sciences (IadChem), Universidad Autónoma de Madrid, Campus de Excelencia UAM-CSIC, Cantoblanco, 28049 Madrid, Spain
| | - I Corral
- Departamento de Química (Módulo 13, Facultad de Ciencias) and Institute of Advanced Chemical Sciences (IadChem), Universidad Autónoma de Madrid, Campus de Excelencia UAM-CSIC, Cantoblanco, 28049 Madrid, Spain
| | - S Faraji
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| |
Collapse
|
29
|
Inai N, Yokogawa D, Yanai T. Investigating the Nonradiative Decay Pathway in the Excited State of Silepin Derivatives: A Study with Second-Order Multireference Perturbation Wavefunction Theory. J Phys Chem A 2021; 125:559-569. [PMID: 33416309 DOI: 10.1021/acs.jpca.0c08738] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The fluorescence quantum yield for fluorescent organic molecules is an important molecular property, and tuning it up is desired for various applications. For the computational estimation of the fluorescence quantum yield, the theoretical prediction of the nonradiative decay rate constant has become an attractive subject of study. The rate constant of thermally activated nonradiative decay is related to the activation energy in the photoreaction; thus, the accuracy and reliability of the excited-state potential energies in the quantum chemical computation are critical. In this study, we employed a second-order multireference perturbation wavefunction theory for studying the thermally activated decay via conical intersection (CI) of 1,1-dimethyldibenzo[b,f]silepin derivatives. The correlation between the computed activation energy to reach the CI geometry in the S1 state and the experimentally determined fluorescence quantum yield implied that silepins nonradiatively decay via the CI triggered by the twisting of the central C-C bond. Geometry optimization of the transition state using multireference perturbation theory drastically reduced the estimated activation energy. Our computation gave reasonable predictions of the activation free energies of photoexcited 1,1-dimethyldibenzo[b,f]silepin. The energy profiles and geometry optimizations using proper quantum chemical methods played a critical role in reliable estimation of the rate constant and fluorescence quantum yield.
Collapse
Affiliation(s)
- Naoto Inai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Daisuke Yokogawa
- Department of Basic Science, Graduate School of Arts and Science, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan.,Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| |
Collapse
|
30
|
Goerigk L, Casanova-Paéz M. The Trip to the Density Functional Theory Zoo Continues: Making a Case for Time-Dependent Double Hybrids for Excited-State Problems. Aust J Chem 2021. [DOI: 10.1071/ch20093] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This account is written for general users of time-dependent density functional theory (TD-DFT) methods as well as chemists who are unfamiliar with the field. It includes a brief overview of conventional TD-DFT approaches and recommendations for applications to organic molecules based on our own experience. The main emphasis of this work, however, lies in providing the first in-depth review of time-dependent double-hybrid density functionals. They were first established in 2007 with very promising follow-up studies in the subsequent four years before developments or applications became scarce. The topic has regained more interest since 2017, and this account reviews those latest developments led by our group. These developments have shown unprecedented robustness for a variety of different types of electronic excitations when compared to more conventional TD-DFT methods. In particular, time-dependent double hybrids do not suffer from artificial ghost states and are able to reproduce exciton-coupled absorption spectra. Our latest methods include range separation and belong to the currently best TD-DFT methods for singlet-singlet excitations in organic molecules. While there is still room for improvement and further development in this space, we hope that this account encourages users to adjust their computational protocols to such new methods to provide more real-life testing and scenarios.
Collapse
|
31
|
Uratani H, Morioka T, Yoshikawa T, Nakai H. Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings. J Chem Theory Comput 2020; 16:7299-7313. [PMID: 33197192 DOI: 10.1021/acs.jctc.0c00936] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Toshiki Morioka
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
32
|
Teh HH, Subotnik JE. Analytic gradients and derivative couplings for configuration interaction with all single excitations and one double excitation-En route to nonadiabatic dynamics. J Chem Phys 2020; 153:184106. [PMID: 33187425 DOI: 10.1063/5.0018441] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present analytic gradients and derivative couplings for the simplest possible multireference configuration interaction method, CIS-1D, an electronic structure Ansatz that includes all single excitations and one lone double excitation on top of a Hartree-Fock reference state. We show that the resulting equations are numerically stable and require the evaluation of a similar number of integrals as compared to standard CIS theory; one can easily differentiate the required frontier orbitals (h and l) with minimal cost. The resulting algorithm has been implemented within the Q-Chem electronic structure package and should be immediately useful for understanding photochemistry with S0-S1 crossings.
Collapse
Affiliation(s)
- Hung-Hsuan Teh
- University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| | - Joseph E Subotnik
- University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| |
Collapse
|
33
|
Gurchiek JK, Rose JB, Guberman-Pfeffer MJ, Tilluck RW, Ghosh S, Gascón JA, Beck WF. Fluorescence Anisotropy Detection of Barrier Crossing and Ultrafast Conformational Dynamics in the S 2 State of β-Carotene. J Phys Chem B 2020; 124:9029-9046. [PMID: 32955881 DOI: 10.1021/acs.jpcb.0c06961] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Carotenoids are usually only weakly fluorescent despite being very strong absorbers in the mid-visible region because their first two excited singlet states, S1 and S2, have very short lifetimes. To probe the structural mechanisms that promote the nonradiative decay of the S2 state to the S1 state, we have carried out a series of fluorescence lineshape and anisotropy measurements with a prototype carotenoid, β-carotene, in four aprotic solvents. The anisotropy values observed in the fluorescence emission bands originating from the S2 and S1 states reveal that the large internal rotations of the emission transition dipole moment, as much as 50° relative to that of the absorption transition dipole moment, are initiated during ultrafast evolution on the S2 state potential energy surface and persist upon nonradiative decay to the S1 state. Electronic structure calculations of the orientation of the transition dipole moment account for the anisotropy results in terms of torsional and pyramidal distortions near the center of the isoprenoid backbone. The excitation wavelength dependence of the fluorescence anisotropy indicates that these out-of-plane conformational motions are initiated by passage over a low-activation energy barrier from the Franck-Condon S2 structure. This conclusion is consistent with detection over the 80-200 K range of a broad, red-shifted fluorescence band from a dynamic intermediate evolving on a steep gradient of the S2 state potential energy surface after crossing the activation barrier. The temperature dependence of the oscillator strength and anisotropy indicate that nonadiabatic passage from S2 through a conical intersection seam to S1 is promoted by the out-of-plane motions of the isoprenoid backbone with strong hindrance by solvent friction.
Collapse
Affiliation(s)
- J K Gurchiek
- Department of Chemistry, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Justin B Rose
- Department of Chemistry, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Matthew J Guberman-Pfeffer
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06268-1712, United States
| | - Ryan W Tilluck
- Department of Chemistry, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Soumen Ghosh
- Department of Physics, Politecnico di Milano, Piazza Leonardo da Vinci, 32, Milan, Lombardy 20133, Italy
| | - José A Gascón
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06268-1712, United States
| | - Warren F Beck
- Department of Chemistry, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| |
Collapse
|
34
|
Yoshikawa T, Doi T, Nakai H. Finite-temperature-based time-dependent density-functional theory method for static electron correlation systems. J Chem Phys 2020; 152:244111. [PMID: 32610978 DOI: 10.1063/1.5144527] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this study, we developed a time-dependent density-functional theory (TDDFT) with a finite-temperature (FT) scheme, denoted as FT-TDDFT. We introduced the concept of fractional occupation numbers for random phase approximation equation and evaluated the excited-state electronic entropy terms with excited-state occupation number. The orbital occupation numbers for the excited state were evaluated from the change in the ground-state electron configuration with excitation and deexcitation coefficients. Furthermore, we extended the FT formulation to the time-dependent density-functional tight-binding (TDDFTB) method for larger systems, denoted as FT-TDDFTB. Numerical assessment for the FT-(TD)DFT method showed smooth potential curves for double-bond rotation of ethylene in both ground and excited states. Excited-state calculations based on the FT-TDDFTB method were applied to the uniform π-stacking columns composed of trioxotriangulene, possessing neutral radicals in strong correlation systems.
Collapse
Affiliation(s)
- Takeshi Yoshikawa
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Toshiki Doi
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| |
Collapse
|
35
|
Yamamoto N. Free Energy Profile Analysis for the Aggregation-Induced Emission of Diphenyldibenzofulvene. J Phys Chem A 2020; 124:4939-4945. [DOI: 10.1021/acs.jpca.0c03240] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Norifumi Yamamoto
- Department of Applied Chemistry, Faculty of Engineering, Chiba Institute of Technology, 2-17-1 Tsudanuma, Narashino, Chiba 275-0016, Japan
| |
Collapse
|
36
|
Salazar E, Faraji S. Theoretical study of cyclohexadiene/hexatriene photochemical interconversion using spin-Flip time-Dependent density functional theory. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1764120] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Edison Salazar
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Groningen, Netherlands
| | - Shirin Faraji
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Groningen, Netherlands
| |
Collapse
|
37
|
Winslow M, Cross WB, Robinson D. Comparison of Spin-Flip TDDFT-Based Conical Intersection Approaches with XMS-CASPT2. J Chem Theory Comput 2020; 16:3253-3263. [PMID: 32302484 PMCID: PMC8279405 DOI: 10.1021/acs.jctc.9b00917] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Determining conical intersection
geometries is of key importance
to understanding the photochemical reactivity of molecules. While
many small- to medium-sized molecules can be treated accurately using
multireference approaches, larger molecules require a less computationally
demanding approach. In this work, minimum energy crossing point conical
intersection geometries for a series of molecules have been studied
using spin-flip TDDFT (SF-TDDFT), within the Tamm-Dancoff Approximation,
both with and without explicit calculation of nonadiabatic coupling
terms, and compared with both XMS-CASPT2 and CASSCF calculated geometries.
The less computationally demanding algorithms, which do not require
explicit calculation of the nonadiabatic coupling terms, generally
fare well with the XMS-CASPT2 reference structures, while the relative
energetics are only reasonably replicated with the MECP structure
as
calculated with the BHHLYP functional and full nonadiabatic coupling
terms. We also demonstrate that, occasionally, CASSCF structures deviate
quantitatively from the XMS-CASPT2 structures, showing the importance
of including dynamical correlation.
Collapse
Affiliation(s)
- Max Winslow
- Department of Chemistry and Forensics, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| | - Warren B Cross
- Department of Chemistry and Forensics, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| | - David Robinson
- Department of Chemistry and Forensics, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| |
Collapse
|
38
|
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
|
39
|
Barca GMJ, Bertoni C, Carrington L, Datta D, De Silva N, Deustua JE, Fedorov DG, Gour JR, Gunina AO, Guidez E, Harville T, Irle S, Ivanic J, Kowalski K, Leang SS, Li H, Li W, Lutz JJ, Magoulas I, Mato J, Mironov V, Nakata H, Pham BQ, Piecuch P, Poole D, Pruitt SR, Rendell AP, Roskop LB, Ruedenberg K, Sattasathuchana T, Schmidt MW, Shen J, Slipchenko L, Sosonkina M, Sundriyal V, Tiwari A, Galvez Vallejo JL, Westheimer B, Włoch M, Xu P, Zahariev F, Gordon MS. Recent developments in the general atomic and molecular electronic structure system. J Chem Phys 2020; 152:154102. [PMID: 32321259 DOI: 10.1063/5.0005188] [Citation(s) in RCA: 541] [Impact Index Per Article: 135.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.
Collapse
Affiliation(s)
- Giuseppe M J Barca
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Laura Carrington
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Dipayan Datta
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Nuwan De Silva
- Department of Physical and Biological Sciences, Western New England University, Springfield, Massachusetts 01119, USA
| | - J Emiliano Deustua
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Jeffrey R Gour
- Microsoft, 15590 NE 31st St., Redmond, Washington 98052, USA
| | - Anastasia O Gunina
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Emilie Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, USA
| | - Taylor Harville
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Stephan Irle
- Computational Science and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, USA
| | - Joe Ivanic
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, USA
| | - Karol Kowalski
- Physical Sciences Division, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, USA
| | - Sarom S Leang
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Hui Li
- Department of Chemistry, University of Nebraska, Lincoln, Nebraska 68588, USA
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, People's Republic of China
| | - Jesse J Lutz
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Joani Mato
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow 119991, Russian Federation
| | - Hiroya Nakata
- Kyocera Corporation, Research Institute for Advanced Materials and Devices, 3-5-3 Hikaridai Seika-cho, Souraku-gun, Kyoto 619-0237, Japan
| | - Buu Q Pham
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - David Poole
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Spencer R Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Alistair P Rendell
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Luke B Roskop
- Cray Inc., a Hewlett Packard Enterprise Company, 2131 Lindau Ln #1000, Bloomington, Minnesota 55425, USA
| | - Klaus Ruedenberg
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | | | - Michael W Schmidt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Lyudmila Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Masha Sosonkina
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Vaibhav Sundriyal
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Ananta Tiwari
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Jorge L Galvez Vallejo
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Bryce Westheimer
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Marta Włoch
- 530 Charlesina Dr., Rochester, Michigan 48306, USA
| | - Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Federico Zahariev
- 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
| |
Collapse
|
40
|
Inamori M, Ikabata Y, Yoshikawa T, Nakai H. Unveiling controlling factors of the S0/S1 minimum energy conical intersection (2): Application to penalty function method. J Chem Phys 2020; 152:144108. [DOI: 10.1063/1.5142592] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Mayu Inamori
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Yasuhiro Ikabata
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
41
|
Locating conical intersections using the quasidegenerate partially and strongly contracted NEVPT2 methods. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137219] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
42
|
Inamori M, Yoshikawa T, Ikabata Y, Nishimura Y, Nakai H. Spin‐flip approach within time‐dependent density functional tight‐binding method: Theory and applications. J Comput Chem 2020; 41:1538-1548. [DOI: 10.1002/jcc.26197] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2020] [Revised: 03/07/2020] [Accepted: 03/10/2020] [Indexed: 12/25/2022]
Affiliation(s)
- Mayu Inamori
- Department of Chemistry and Biochemistry, School of Advanced Science and EngineeringWaseda University Tokyo Japan
| | - Takeshi Yoshikawa
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
| | - Yasuhiro Ikabata
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and EngineeringWaseda University Tokyo Japan
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB)Kyoto University Kyoto Japan
| |
Collapse
|
43
|
Marazzi M, Francés-Monerris A, Mourer M, Pasc A, Monari A. Trans-to-cis photoisomerization of cyclocurcumin in different environments rationalized by computational photochemistry. Phys Chem Chem Phys 2020; 22:4749-4757. [PMID: 32057038 DOI: 10.1039/c9cp06565b] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cyclocurcumin is a turmeric component that has attracted much less attention compared to the well-known curcumin. In spite of the less deep characterization of its properties, cyclocurcumin has shown promising anticancer effects when used in combination with curcumin. Especially, due to its peculiar molecular structure, cyclocurcumin can be regarded as an almost ideal photoswitch, whose capabilities can also be exploited for relevant biological applications. Here, by means of state-of-the-art computational methods for electronic excited-state calculations (TD-DFT, MS-CASPT2, and XMS-CASPT2), we analyze in detail the absorption and photoisomerization pathways leading from the more stable trans isomer to the cis one. The different molecular surroundings, taken into account by means of the electrostatic solvent effect and compared with available experimental data, have been found to be critical in describing the fate of irradiated cyclocurcumin: when in non-polar environments, an excited state barrier prevents photoisomerization and favours fluorescence, whereas in polar solvents, an almost barrierless path results in a striking decrease of fluorescence, opening the way toward a crossing region with the ground state and thus funneling the photoproduction of the cis isomer.
Collapse
Affiliation(s)
- Marco Marazzi
- Department of Analytical Chemistry, Physical Chemistry and Chemical Engineering, Universidad de Alcalá, Ctra. Madrid-Barcelona, Km 33, 600, E-28805, Alcalá de Henares, Madrid, Spain. and Chemical Research Institute "Andrés M. del Rio" (IQAR), Universidad de Alcalá, Ctra. Madrid-Barcelona, Km 33, 600, E-28805, Alcalá de Henares, Madrid, Spain
| | - Antonio Francés-Monerris
- Université de Lorraine and CNRS, LPCT UMR-7019, F-5400 Nancy, France. and Departement de Quimica Fisica, Universitat de Valencia, 46100 Burjassot, Spain
| | - Maxime Mourer
- Université de Lorraine and CNRS, L2CM UMR 7053, F-54000 Nancy, France
| | - Andreea Pasc
- Université de Lorraine and CNRS, L2CM UMR 7053, F-54000 Nancy, France
| | - Antonio Monari
- Université de Lorraine and CNRS, LPCT UMR-7019, F-5400 Nancy, France.
| |
Collapse
|
44
|
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
|
45
|
Abstract
This Perspective discusses salient features of the spin-flip approach to strong correlation and describes different methods that sprung from this idea. The spin-flip treatment exploits the different physics of low-spin and high-spin states and is based on the observation that correlation is small for same-spin electrons. By using a well-behaved high-spin state as a reference, one can access problematic low-spin states by deploying the same formal tools as in the excited-state treatments (i.e., linear response, propagator, or equation-of-motion theories). The Perspective reviews applications of this strategy within wave function and density functional theory frameworks as well as the extensions for molecular properties and spectroscopy. The utility of spin-flip methods is illustrated by examples. Limitations and proposed future directions are also discussed.
Collapse
Affiliation(s)
- David Casanova
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain. and IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Euskadi, Spain
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
46
|
Mato J, Gordon MS. Analytic non-adiabatic couplings for the spin-flip ORMAS method. Phys Chem Chem Phys 2020; 22:1475-1484. [DOI: 10.1039/c9cp05849d] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Analytic non-adiabatic coupling matrix elements (NACME) are derived and implemented for the spin-flip occupation restricted multiple active space configuration interaction (SF-ORMAS-CI) method.
Collapse
Affiliation(s)
- Joani Mato
- Department of Chemistry and Ames Laboratory
- Iowa State University
- Ames
- USA
| | - Mark S. Gordon
- Department of Chemistry and Ames Laboratory
- Iowa State University
- Ames
- USA
| |
Collapse
|
47
|
Marut C, Senjean B, Fromager E, Loos PF. Weight dependence of local exchange–correlation functionals in ensemble density-functional theory: double excitations in two-electron systems. Faraday Discuss 2020; 224:402-423. [DOI: 10.1039/d0fd00059k] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We discuss the construction of first-rung weight-dependent exchange–correlation density-functional approximations for He and H2 specifically designed for the computation of double excitations within Gross–Oliveira–Kohn-DFT.
Collapse
Affiliation(s)
- Clotilde Marut
- Laboratoire de Chimie et Physique Quantiques
- Université de Toulouse
- CNRS
- UPS
- France
| | - Bruno Senjean
- Instituut-Lorentz
- Universiteit Leiden
- 2300 RA Leiden
- The Netherlands
- Division of Theoretical Chemistry
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique
- Institut de Chimie
- CNRS
- Université de Strasbourg
- Strasbourg
| | | |
Collapse
|
48
|
de Silva P. Inverted Singlet-Triplet Gaps and Their Relevance to Thermally Activated Delayed Fluorescence. J Phys Chem Lett 2019; 10:5674-5679. [PMID: 31483656 DOI: 10.1021/acs.jpclett.9b02333] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The basic design principle for emitters exhibiting thermally activated delayed fluorescence (TADF) is the minimization of the singlet-triplet gap. While typically this gap is positive, a possible inversion of states has been proposed as a pathway to improve the efficiency of organic light-emitting diodes. Despite the efforts to design such emitters, there are very few reports indicating that it is at all possible. We analyze the problem of the gap inversion from the perspective of the electronic structure theory. The key result is that inversion is possible but requires a substantial contribution of double excitations and that commonly used cheap electronic structure methods would fail to predict it.
Collapse
Affiliation(s)
- Piotr de Silva
- Department of Energy Conversion and Storage , Technical University of Denmark , Fysikvej 309 , 2800 Kongens Lyngby , Denmark
| |
Collapse
|
49
|
Pederzoli M, Wasif Baig M, Kývala M, Pittner J, Cwiklik L. Photophysics of BODIPY-Based Photosensitizer for Photodynamic Therapy: Surface Hopping and Classical Molecular Dynamics. J Chem Theory Comput 2019; 15:5046-5057. [PMID: 31390517 DOI: 10.1021/acs.jctc.9b00533] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Halogenated BODIPY derivatives are emerging as important candidates for photodynamic therapy of cancer cells due to their high triplet quantum yield. We probed fundamental photophysical properties and interactions with biological environments of such photosensitizers. To this end, we employed static TD-DFT quantum chemical calculations as well as TD-DFT surface hopping molecular dynamics on potential energy surfaces resulting from the eigenstates of the total electronic Hamiltonian including the spin-orbit (SO) coupling. Matrix elements of an effective one-electron spin-orbit Hamiltonian between singlet and triplet configuration interaction singles (CIS) auxiliary wave functions are calculated using a new code capable of dealing with singlets and both restricted and unrestricted triplets built up from up to three different and independent sets of (singlet, alpha, and beta) molecular orbitals. The interaction with a biological environment was addressed by using classical molecular dynamics (MD) in a scheme that implicitly accounts for electronically excited states. For the surface hopping trajectories, an accelerated MD approach was used, in which the SO couplings are scaled up, to make the calculations computationally feasible, and the lifetimes are extrapolated back to unscaled SO couplings. The lifetime of the first excited singlet state estimated by semiclassical surface hopping simulations is 139 ± 75 ps. Classical MD demonstrates that halogenated BODIPY in the ground state, in contrast to the unsubstituted one, is stable in the headgroup region of minimalistic cell membrane models, and while in the triplet state, the molecule relocates to the membrane interior ready for further steps of photodynamic therapy.
Collapse
Affiliation(s)
- Marek Pederzoli
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences , Dolejškova 2155/3 , 18223 Prague 8 , Czech Republic.,Faculty of Science, Department of Physical and Macromolecular Chemistry , Charles University in Prague , Hlavova 8 , 12840 Prague , Czech Republic
| | - Mirza Wasif Baig
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences , Dolejškova 2155/3 , 18223 Prague 8 , Czech Republic.,Faculty of Science, Department of Physical and Macromolecular Chemistry , Charles University in Prague , Hlavova 8 , 12840 Prague , Czech Republic
| | - Mojmír Kývala
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo nám. 2 , 16610 Prague 6 , Czech Republic
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences , Dolejškova 2155/3 , 18223 Prague 8 , Czech Republic
| | - Lukasz Cwiklik
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences , Dolejškova 2155/3 , 18223 Prague 8 , Czech Republic.,Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo nám. 2 , 16610 Prague 6 , Czech Republic
| |
Collapse
|
50
|
Levine BG, Esch MP, Fales BS, Hardwick DT, Peng WT, Shu Y. Conical Intersections at the Nanoscale: Molecular Ideas for Materials. Annu Rev Phys Chem 2019; 70:21-43. [DOI: 10.1146/annurev-physchem-042018-052425] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The ability to predict and describe nonradiative processes in molecules via the identification and characterization of conical intersections is one of the greatest recent successes of theoretical chemistry. Only recently, however, has this concept been extended to materials science, where nonradiative recombination limits the efficiencies of materials for various optoelectronic applications. In this review, we present recent advances in the theoretical study of conical intersections in semiconductor nanomaterials. After briefly introducing conical intersections, we argue that specific defects in materials can induce conical intersections between the ground and first excited electronic states, thus introducing pathways for nonradiative recombination. We present recent developments in theoretical methods, computational tools, and chemical intuition for the prediction of such defect-induced conical intersections. Through examples in various nanomaterials, we illustrate the significance of conical intersections for nanoscience. We also discuss challenges facing research in this area and opportunities for progress.
Collapse
Affiliation(s)
- Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - B. Scott Fales
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Dylan T. Hardwick
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|