1
|
Wang Y, Verma P, Zhang L, Li Y, Liu Z, Truhlar DG, He X. M06-SX screened-exchange density functional for chemistry and solid-state physics. Proc Natl Acad Sci U S A 2020; 117:2294-2301. [PMID: 31953258 PMCID: PMC7007546 DOI: 10.1073/pnas.1913699117] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Screened-exchange hybrid density functionals are especially recommended for solid-state systems because they combine the advantages of hybrid functionals with the correct physics and lower computational cost associated with the attenuation of Hartree-Fock exchange at long range. We present a screened-exchange hybrid functional, M06-SX, that combines the functional form of the local revM06-L functional with a percentage of short-range nonlocal Hartree-Fock exchange. The M06-SX functional gives good results not only for a large set of training data but also for several databases quite different from the training data. The mean unsigned error (MUE) of the M06-SX functional is 2.85 kcal/mol for 418 atomic and molecular energies (AME418) in Minnesota Database 2019, which is better than all five other screened-exchange hybrid functionals tested in this work. The M06-SX functional also gives especially good results for semiconductor band gaps, molecular dissociation energies, noncovalent interactions, barrier heights, and electronic excitation energies excluding long-range charge transfer excitations. For the LC18 lattice constants database, the M06-SX functional gives an MUE of only 0.034 Å. Therefore, the M06-SX functional is well suited for studying molecular chemistry as well as solid-state physics.
Collapse
Affiliation(s)
- Ying Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Pragya Verma
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455-0431
- Chemical Theory Center, University of Minnesota, Minneapolis, MN 55455-0431
- Nanoporous Materials Genome Center, University of Minnesota, Minneapolis, MN 55455-0431
- Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Lujia Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| | - Yaqi Li
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
| | - Zhonghua Liu
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
| | - Donald G Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455-0431;
- Chemical Theory Center, University of Minnesota, Minneapolis, MN 55455-0431
- Nanoporous Materials Genome Center, University of Minnesota, Minneapolis, MN 55455-0431
- Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
2
|
Verma P, Wang Y, Ghosh S, He X, Truhlar DG. Revised M11 Exchange-Correlation Functional for Electronic Excitation Energies and Ground-State Properties. J Phys Chem A 2019; 123:2966-2990. [DOI: 10.1021/acs.jpca.8b11499] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Pragya Verma
- Department of Chemistry, Chemical Theory Center, Nanoporous Materials Genome Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Ying Wang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
| | - Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, Nanoporous Materials Genome Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry, New York University, Shanghai, Shanghai 200062, China
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, Nanoporous Materials Genome Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
3
|
Grabarek D, Andruniów T. Initial excited-state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations. J Comput Chem 2018; 39:1720-1727. [PMID: 29727036 DOI: 10.1002/jcc.25346] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2018] [Revised: 04/06/2018] [Accepted: 04/12/2018] [Indexed: 11/07/2022]
Abstract
The initial S1 excited-state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight-membered ring (locked-11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS-CASPT2, MS-CASPT2) as well as single-reference coupled-cluster CC2 method. The analysis of XMCQDPT2-based geometries reveals rather weak coupling between in-plane and out-of-plane structural evolution and minor energetical relaxation of three locked-11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out-of-plane deformation resulting in a very steep S1 energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground-state structures, the excited-state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4-0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked-11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| |
Collapse
|
4
|
Wolański Ł, Grabarek D, Andruniów T. Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal? J Comput Chem 2018; 39:1470-1480. [PMID: 29635695 DOI: 10.1002/jcc.25217] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2017] [Revised: 03/09/2018] [Accepted: 03/10/2018] [Indexed: 11/07/2022]
Abstract
To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S0 →S1 and S0 →S2 energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Łukasz Wolański
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| |
Collapse
|
5
|
Lischka H, Nachtigallová D, Aquino AJA, Szalay PG, Plasser F, Machado FBC, Barbatti M. Multireference Approaches for Excited States of Molecules. Chem Rev 2018; 118:7293-7361. [DOI: 10.1021/acs.chemrev.8b00244] [Citation(s) in RCA: 197] [Impact Index Per Article: 32.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hans Lischka
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Palacký University, 78371 Olomouc, Czech Republic
| | - Adélia J. A. Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute for Soil Research, University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Laboratory of Theoretical Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Leicestershire LE11 3TU, United Kingdom
| | - Francisco B. C. Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | | |
Collapse
|
6
|
Dong SS, Gagliardi L, Truhlar DG. Excitation spectra of retinal by multiconfiguration pair-density functional theory. Phys Chem Chem Phys 2018; 20:7265-7276. [PMID: 29484326 DOI: 10.1039/c7cp07275a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Retinal is the chromophore in proteins responsible for vision. The absorption maximum of retinal is sensitive to mutations of the protein. However, it is not easy to predict the absorption spectrum of retinal accurately, and questions remain even after intensive investigation. Retinal poses a challenge for Kohn-Sham density functional theory (KS-DFT) because of the charge transfer character in its excitations, and it poses a challenge for wave function theory because the large size of the molecule makes multiconfigurational perturbation theory methods expensive. In this study, we demonstrate that multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to predict the vertical excitation energies of 11-Z retinal, and it reproduces the experimentally determined absorption band widths and peak positions better than complete active space second-order perturbation theory (CASPT2). The consistency between complete active space self-consistent field (CASSCF) and KS-DFT dipole moments is demonstrated to be a useful criterion in selecting the active space. We also found that the nature of the terminal groups and the conformations of retinal play a significant role in the absorption spectrum. By considering a thermal distribution of conformations, we predict an absorption spectrum of retinal that is consistent with the experimental gas-phase spectrum. The location of the absorption peak and the spectral broadening based on MC-PDFT calculations agree better with experiments than those of CASPT2.
Collapse
Affiliation(s)
- Sijia S Dong
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455, USA.
| | | | | |
Collapse
|
7
|
Gozem S, Luk HL, Schapiro I, Olivucci M. Theory and Simulation of the Ultrafast Double-Bond Isomerization of Biological Chromophores. Chem Rev 2017; 117:13502-13565. [DOI: 10.1021/acs.chemrev.7b00177] [Citation(s) in RCA: 175] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Samer Gozem
- Department
of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Hoi Ling Luk
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Fritz
Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro
2, 53100 Siena, Italy
| |
Collapse
|
8
|
Varsano D, Caprasecca S, Coccia E. Theoretical description of protein field effects on electronic excitations of biological chromophores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:013002. [PMID: 27830666 DOI: 10.1088/0953-8984/29/1/013002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.
Collapse
Affiliation(s)
- Daniele Varsano
- S3 Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | | | | |
Collapse
|
9
|
Guido CA, Cortona P, Adamo C. Effective electron displacements: A tool for time-dependent density functional theory computational spectroscopy. J Chem Phys 2014; 140:104101. [DOI: 10.1063/1.4867007] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
|
10
|
Valsson O, Campomanes P, Tavernelli I, Rothlisberger U, Filippi C. Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls. J Chem Theory Comput 2013; 9:2441-54. [PMID: 26583734 DOI: 10.1021/ct3010408] [Citation(s) in RCA: 71] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.
Collapse
Affiliation(s)
- Omar Valsson
- MESA+ Institute for Nanotechnology, University of Twente , P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Pablo Campomanes
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Ivano Tavernelli
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente , P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
11
|
Coccia E, Varsano D, Guidoni L. Protein Field Effect on the Dark State of 11- cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics. J Chem Theory Comput 2013; 9:8-12. [PMID: 24611033 PMCID: PMC3943175 DOI: 10.1021/ct3007502] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate determination of the geometrical details of the dark state of 11-cis retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the β-ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [J. Am. Chem. Soc. 2004, 126, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S1 and S2 states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.
Collapse
Affiliation(s)
- Emanuele Coccia
- Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, via Vetoio, 67100, L’Aquila, Italy
| | - Daniele Varsano
- Dipartimento di Fisica, “Sapienza” - Università di Roma, piazzale Aldo Moro 5, 00185, Rome, Italy
| | - Leonardo Guidoni
- Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, via Vetoio, 67100, L’Aquila, Italy
| |
Collapse
|
12
|
Coccia E, Guidoni L. Quantum monte carlo study of the retinal minimal model C5H6NH2+. J Comput Chem 2012; 33:2332-9. [DOI: 10.1002/jcc.23071] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2012] [Revised: 06/22/2012] [Accepted: 06/26/2012] [Indexed: 12/22/2022]
|
13
|
Valsson O, Angeli C, Filippi C. Excitation energies of retinal chromophores: critical role of the structural model. Phys Chem Chem Phys 2012; 14:11015-20. [DOI: 10.1039/c2cp41387f] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
14
|
Nishizaka M, Mori T, Inoue Y. Axial chirality of donor-donor, donor-acceptor, and tethered 1,1'-binaphthyls: a theoretical revisit with dynamics trajectories. J Phys Chem A 2011; 115:5488-95. [PMID: 21557592 DOI: 10.1021/jp202776g] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The circular dichroism (CD) spectra of (R)-2,2'-dimethoxy-1,1'-binaphthyl (DD) and its untethered and tethered donor-acceptor analogues (DA and DA7-DA9) were investigated experimentally and theoretically. The experimental CD spectra of DD and DA resembled each other in several aspects, displaying a positive-positive-negative Cotton effect pattern in the (1)L(b)-(1)L(a) region and a strong negative couplet at the (1)B(b) band, but significantly differed in transition energy and rotatory strength. The couplet amplitude (A) of the main band was 1.6 times larger in DA than in DD, despite the comparable extinction coefficients and seemingly analogous conformations. An additional positive Cotton effect was observed at the CT (CT) band for donor-acceptor binaphthyl DA. Our theoretical prediction of the CD spectra of binaphthyls involves three sequential first principle quantum mechanics (QM) calculations. Thus, the geometry optimizations of a series of conformers with varying dihedral angles were performed by the dispersion-corrected DFT-D method using the B97-D functional and the TZV2P basis set. The potential curve as a function of the dihedral angle (θ) was obtained by using the SCS-MP2/TZVPP single-point energy calculations with and without application of the solvent correction. The CD spectrum of each conformer was independently calculated by the second-order approximate coupled cluster calculation (CC2 method) using the TZVPP basis sets and the resolution of the identity (RI-J) approximation. The (net) theoretical CD spectrum was obtained by averaging over all possible conformers, where the dynamics trajectories based on the relative SCS-MP2 energies were taken into account. By using 17 possible conformers at θ varying from 50 to 130° by 5° intervals, the experimental CD spectra were successfully reproduced in a quantitative manner, enabling us to characterize properly almost all of the important spectral features and chiroptical properties. The two-state model, reported previously, turned out to have led to the right answer with wrong reasons. The couplet sign and amplitude A are critical functions of θ and can be used not only for (qualitatively) determining the absolute configuration but also for quantitatively analyzing the binaphthyl conformations. The angle dependence of A was already argued in the classical coupled oscillator and exciton chirality theories to provide reasonable structure elucidations but only in a qualitative or semiquantitative manner. Our method is able to predict the A value quantitatively as a function of θ. For tethered binaphthyls DA7-DA9, particular care should be exercised in the conformational assessment based on the classical treatment because the amplitude A was shown to be significantly affected by the existence of the tether itself. In the present method, the couplet amplitude A was nicely related to the dihedral angle θ of DA and DD by the state-of-the-art ab initio calculations, enabling us to gain the quantitative information about the conformation of axially chiral binaphthyls. The Cotton effect at the CT band also serves as a complementary clue for elucidating the conformation of donor-acceptor binaphthyls.
Collapse
Affiliation(s)
- Masaki Nishizaka
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-oka, Suita 565-0871, Japan
| | | | | |
Collapse
|
15
|
Li R, Zheng J, Truhlar DG. Density functional approximations for charge transfer excitations with intermediate spatial overlap. Phys Chem Chem Phys 2010; 12:12697-701. [PMID: 20733991 DOI: 10.1039/c0cp00549e] [Citation(s) in RCA: 89] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Density functional theory is now the method of choice for calculating the electronic structure of complex systems, and time-dependent density functional theory (TDDFT) is now the preferred method for calculating spectroscopic properties of large molecules. The validity of the theory depends mainly on the quality of the approximation to the unknown exchange-correlation energy. In the present paper we consider TDDFT calculations of electronic excitation energies and oscillator strengths. We show that the M06-2X and M08-HX density functionals perform as well as and better than the range-separated CAM-B3LYP functional for charge transfer excitations with intermediate spatial overlap but have better performance for bond energies, noncovalent interactions, and chemical reaction barrier heights for representative systems; we conclude that M06-2X and M08-HX should be preferred for studies requiring the exploration of potential energy surfaces as well as electronic excitation energies, provided that those excitations with the longest-range charge transfer are excluded.
Collapse
Affiliation(s)
- Ruifang Li
- Department of Chemistry, Nankai University, Tianjin, 300071, PR China.
| | | | | |
Collapse
|
16
|
Valsson O, Filippi C. Photoisomerization of Model Retinal Chromophores: Insight from Quantum Monte Carlo and Multiconfigurational Perturbation Theory. J Chem Theory Comput 2010. [DOI: 10.1021/ct900692y] [Citation(s) in RCA: 120] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Omar Valsson
- Faculty of Science and Technology and MESA+ Research Institute, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Claudia Filippi
- Faculty of Science and Technology and MESA+ Research Institute, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
17
|
Dwyer AD, Tozer DJ. Effect of chemical change on TDDFT accuracy: orbital overlap perspective of the hydrogenation of retinal. Phys Chem Chem Phys 2010; 12:2816-8. [DOI: 10.1039/c002428g] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
18
|
Send R, Sundholm D, Johansson MP, Pawłowski F. Excited State Potential Energy Surfaces of Polyenes and Protonated Schiff Bases. J Chem Theory Comput 2009; 5:2401-14. [DOI: 10.1021/ct900240s] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Robert Send
- Institut für Physikalische Chemie, Universität Karlsruhe, Kaiserstrasse 12, D-76128 Karlsruhe, Germany
| | - Dage Sundholm
- Department of Chemistry, P.O. Box 55 (A.I. Virtanens plats 1), University of Helsinki, FI-00014 Helsinki, Finland
| | - Mikael P. Johansson
- Lundbeck Foundation Centre for Theoretical Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Århus C, Denmark
| | - Filip Pawłowski
- Physics Institute, Kazimierz Wielki University, Plac Weyssenhoffa 11, PL-85-072 Bydgoszcz, Poland
| |
Collapse
|