Electronic structure of the S1 state in methylcobalamin: Insight from CASSCF/MC-XQDPT2, EOM-CCSD, and TD-DFT calculations
The methylcobalamin cofactor (MeCbl), which is one of the biologically active forms of vitamin B12, has been the subject of many spectroscopic and theoretical investigations. Traditionally, the lowest‐energy part of the photoabsorption spectrum of MeCbl (the so‐called α/β band) has been interpreted...
Saved in:
Published in | Journal of computational chemistry Vol. 34; no. 12; pp. 987 - 1004 |
---|---|
Main Authors | , , , , , , , |
Format | Journal Article |
Language | English |
Published |
Hoboken
Wiley Subscription Services, Inc., A Wiley Company
05.05.2013
|
Subjects | |
Online Access | Get full text |
Cover
Loading…
Summary: | The methylcobalamin cofactor (MeCbl), which is one of the biologically active forms of vitamin B12, has been the subject of many spectroscopic and theoretical investigations. Traditionally, the lowest‐energy part of the photoabsorption spectrum of MeCbl (the so‐called α/β band) has been interpreted as an S0→S1 electronic transition dominated by π→π* excitations associated with the CC stretching of the corrin ring. However, a more quantitative band‐shape analysis of the α/β spectral region, along with circular dichroism (CD), magnetic CD, and resonance Raman data, has revealed the presence of a second electronic transition that involves the CoCMe bond weakening. Conversely, the lowest‐energy excitations based on transient absorption spectroscopy measurements have been interpreted as metal‐to‐ligand charge transfer (MLCT) transitions. To resolve the existing controversy about the interpretation of the S1 state of MeCbl, calculations have been performed using two independent ab initio wavefunction‐based methods. These include the modified variant of the second‐order multiconfigurational quasi‐degenerate perturbation theory (MC‐XQDPT2), using complete active space self‐consistent field orbitals, and the equation‐of‐motion coupled‐cluster singles and doubles (EOM‐CCSD) approach using restricted Hartree–Fock orbitals. It is shown that both ab initio methods provide a consistent description of the S1 state as having an MLCT character. In addition, the performance of different types of functionals, including hybrid (B3LYP, MPW1PW91, TPSSh), generalized‐gradient‐approximation‐type (GGA‐type) (BP86, BLYP, MPWPW91), meta‐GGA (TPSS), and range‐separated (CAM‐B3LYP, LC‐BLYP) approaches, has been examined and the results of the corresponding time‐dependent density functional theory calculations have been benchmarked against the MC‐XQDPT2 and EOM‐CCSD data. The hybrid functionals support the interpretation in which the S1 state represents a π→π* transition localized on corrin, while pure GGA, meta‐GGA, and LC‐BLYP functionals produce results consistent with the MLCT assignment. © 2013 Wiley Periodicals, Inc.
The existing controversy regarding the π → π* versus metal‐to‐ligand charge transfer (MLCT) interpretation of the lowestenergy part of the photoabsorbtion spectrum (α/β band) of methylcobalamin (MeCbl) is explored theoretically. Time‐dependent density functional theory calculations are benchmarked against two ab initio wavefunction‐based methodologies: the second‐order multireference perturbation theory and the equation‐of‐motion coupled‐cluster singles and doubles approach. The advantage of generalizedgradient‐ approximation‐type over hybrid functionals is demonstrated, indicating the MLCT nature of the S0 → S1 transition. |
---|---|
Bibliography: | istex:168C0960DFF892F34D39D08CE6088FA4B01FCED0 US Department of Energy's National Nuclear Security Administration - No. DE-AC04-94AL85000 ArticleID:JCC23204 ark:/67375/WNG-XP2S5S4R-W Office of Basic Energy Sciences, Office of Science, US Department of Energy (to PP) - No. DE-FG02-01ER15228 ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 23 |
ISSN: | 0192-8651 1096-987X |
DOI: | 10.1002/jcc.23204 |