Time-dependent density functional theory (TD-DFT) and correlated ab initio methods have been applied to explore the electronically excited states of vitamin B(12) (cyanocobalamin or CNCbl). Different experimental techniques have been used to probe the excited states of CNCbl, revealing many issues that remain poorly understood from an electronic structure point of view. Due to its efficient scaling with size, TD-DFT emerges as one of the most practical tools that can be used to study the electronic properties of these fairly complex molecules. However, the description of excited states is strongly dependent on the type of functional used in the calculations. In the present contribution, the choice of a proper functional for vitamin B(12) was evaluated in terms of its agreement with both experimental results and correlated ab initio calculations. Three different functionals, i.e., B3LYP, BP86, and LC-BLYP, were tested. In addition, the effect of the relative contributions of DFT and HF to the exchange-correlation functional was investigated as a function of the range-separation parameter, μ. The issues related to the underestimation of charge-transfer excitation energies by TD-DFT were validated by the Λ diagnostic, which measures the spatial overlap between occupied and virtual orbitals involved in the particular excitation. The nature of the low-lying excited states was also analyzed based on a comparison of TD-DFT and ab initio results. Based on an extensive comparison with experimental results and ab initio benchmark calculations, the BP86 functional was found to be the most appropriate in describing the electronic properties of CNCbl. Finally, an analysis of electronic transitions and reassignment of some excitations were discussed.

Linear and quadratic response time-dependent density functional theory (TD-DFT) has been applied to investigate absorption (Abs), circular dichroism (CD), and magnetic CD (MCD) spectra of cyanocobalamin (CNCbl) and methylcobalamin (MeCbl). Although electronically excited states of both cobalamins have been probed by applying different experimental techniques, their exact nature remains poorly understood from an electronic structure point of view. Recent theoretical studies have revealed a lot of relevant information about their properties but also left some unresolved issues related to the nature of individual transitions. In this contribution, not only Abs but also CD and MCD spectra of both cobalamins were computed for direct comparison with experiment. The results were evaluated with respect to the choice of exchange-correlation functional, basis set, and the environment (gas phase or solvent) used in the calculation. Taking into account the complexity of the CNCbl and MeCbl systems, reliable agreement between theory and experiment was achieved based on calculations employing the BP86 functional, particularly for the low-energy α/β bands. This spectral range has been traditionally interpreted as a vibrational progression associated with a single electronic excitation, but according to the present analysis for both cobalamins, these bands are best interpreted as consisting of multiple electronic transitions.

Time-dependent density functional theory (TD-DFT) has been applied to explore electronically excited states of vitamin B(12) (cyanocobalamin or CNCbl). To explain why the Co-C bond in CNCbl does not undergo photodissociation under conditions of simple photon excitation, electronically excited states have been computed along the Co-C(CN) stretched coordinate. It was found that the repulsive (3)(σ(Co-C) → σ*(Co-C)) triplet state drops in energy as the Co-C(CN) bond lengthens, but it does not become dissociative. Low-lying excited states were also computed as function of two axial bond lengths. Two energy minima have been located on the S(1)/CNCbl, as well as T(1)/CNCbl, surfaces. The full geometry optimization was carried out for each minimum and electronic properties associated with each optimized structure were analyzed in details. One minimum was described as excitation having mixed ππ*/MLCT (metal-to-ligand charge transfer) character, while the second as ligand-to-metal charge transfer (LMCT) transition. Neither of them, however, can be viewed as pure MLCT or LMCT transitions since additional excitation to or from σ-bonds (SB) of N-Co-C unit have also noticeable contributions. Inclusion of solvent altered the character of one of the excitations from ππ*/MLCT/SBLCT to ππ*/LMCT/LSBCT-type, and therefore, both of them gained significant contribution from LMCT/LSBCT transition. Finally, the nature of S(1) electronic state has been comparatively analyzed in CNCbl and MeCbl cobalamins.

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 C=C 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 Co-C(Me) 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.