Density functional calculations on electronic circular dichroism spectra of chiral transition metal complexes

Optically active bis(aminophenols) and their metal complexes

Optically active C2-symmetric bis(aminophenols) based on (R)-2,2'-diaminobinaphthyl (BiniqH4) and (R,R)-2,3-butanediyldianthranilate (BdanH4) have been prepared by condensation of the diamines with 3,5-di-tert-butylcatechol. Group 10 bis(iminosemiquinone) complexes (R)-(Biniq)M (M = Pd, Pt) and (C,R,R)-(Bdan)Pd have been prepared by oxidatively metalating the corresponding ligands. In (R)-(Biniq)M, the C2 axis passes through the approximate square plane of the bis(iminosemiquinone)metal core, while in (C,R,R)-(Bdan)Pd the C2 axis is perpendicular to this plane. In the latter compound, the (R,R)-butanediyl strap binds selectively over one enantioface of the metal complex in a conformation where the methyl groups are anti to one another. Osmium oxo complexes with the intrinsically chiral OsO(amidophenoxide)2 chromophore are obtained by metalation of OsO(OCH2CH2O)2 with (R,R)-BdanH4. Both the (A,R,R) and (C,R,R) diastereomers can be observed, with metalation in refluxing toluene selectively giving the latter isomer. The electronic structures of the complexes are illuminated by the circular dichroism spectra, in conjuction with the optical spectra and TDDFT calculations.

Infrared Ion Spectroscopic Characterization of the Gaseous [Co(15-crown-5)(H2O)]2+ Complex

We report fingerprint infrared multiple-photon dissociation spectra of the gaseous monohydrated coordination complex of cobalt(II) and the macrocycle 1,4,7,10,13-pentaoxacyclopentadecane (or 15-crown-5), [Co(15-crown-5)(H2O)]2+. The metal-ligand complexes are generated using electrospray ionization, and their IR action spectra are recorded in a quadrupole ion trap mass spectrometer using the free-electron laser FELIX. The electronic structure and chelation motif are derived from spectral comparison with computed vibrational spectra obtained at the density functional theory level. We focus here on the gas-phase structure, addressing the question of doublet versus quartet spin multiplicity and the chelation geometry. We conclude that the gas-phase complex adopts a quartet spin state, excluding contributions of doublet species, and that the chelation geometry is pseudo-octahedral with the six oxygen centers of 15-crown-5 and H2O coordinated to the metal ion. We also address the possible presence of higher-energy conformers based on the IR spectral evidence and calculated thermodynamics.

Facial vs. meridional coordination in gaseous Ni(II)-hexacyclen complexes revealed with infrared ion spectroscopy

We report fingerprint infrared multiple-photon dissociation (IRMPD) spectra of the isolated gaseous hexa-coordinated complex of the macrocycle hexa-aza-18-crown-6 (hexacyclen, 1,4,7,10,13,16-hexaazacyclooctadecane, 18-azacrown-6) with Ni2+. The metal-ligand complexes are generated using electrospray ionization (ESI) and IR action spectra are recorded in a Fourier transform ion cyclotron resonance mass spectrometer (FTICR) MS coupled to the infrared free-electron laser FELIX. We investigate geometric structure of the complexes and in particular the chelation motif, by comparison with computed vibrational spectra, obtained using density functional theory (DFT) at the B3LYP/6-31++G(d,p) level. The quasi-octahedral chelation motif of the complex has been well documented in condensed-phase studies, and we focus here on the gas-phase structure, addressing in particular the question of a facial (fac) versus a meridional (mer) octahedral chelation geometry. Based on the good agreement between calculated linear IR spectra and experimental IRMPD spectra, we conclude that the gas-phase complex adopts a mer chelation geometry and we exclude significant contribution of the fac isomer, which is computed to lie about 10 kJ mol-1 higher in energy. We also address the possible presence of both meridional diastereomers and of higher energy conformers of meridional isomers. Finally, as expected for the d8 Ni2+-ion in an octahedral ligand environment, the IR spectrum also shows that the complexes are in a high-spin electron configuration.

Optical Activity of Spin-Forbidden Electronic Transitions in Metal Complexes from Time-Dependent Density Functional Theory with Spin-Orbit Coupling

The calculation of magnetic transition dipole moments and rotatory strengths was implemented at the zeroth‐order regular approximation (ZORA) two‐component relativistic time‐dependent density functional theory (TDDFT) level. The circular dichroism of the spin‐forbidden ligand‐field transitions of tris(ethylenediamine)cobalt(III) computed in this way agrees very well with available measurements. Phosphorescence dissymmetry factors glum and the corresponding lifetimes are evaluated for three N‐heterocyclic‐carbene‐based iridium complexes, two of which contain helicene moieties, and for two platinahelicenes. The agreement with experimental data is satisfactory. The calculations reproduce the signs and order of magnitude of glum , and the large variations of phosphorescence lifetimes among the systems. The electron spin contribution to the magnetic transition dipole moment is shown to be important in all of the computations.

Chiral Resolution of Spin-Crossover Active Iron(II) [2x2] Grid Complexes

Chiral magnetic materials are proposed for applications in second‐order non‐linear optics, magneto‐chiral dichroism, among others. Recently, we have reported a set of tetra‐nuclear Fe(II) grid complex conformers with general formula C/S‐[Fe₄L₄]⁸⁺ (L: 2,6‐bis(6‐(pyrazol‐1‐yl)pyridin‐2‐yl)‐1,5‐dihydrobenzo[1,2‐d : 4,5‐d′]diimidazole). In the grid complexes, isomerism emerges from tautomerism and conformational isomerism of the ligand L, and the S‐type grid complex is chiral, which originates from different non‐centrosymmetric spatial organization of the trans type ligand around the Fe(II) center. However, the selective preparation of an enantiomerically pure grid complex in a controlled manner is difficult due to spontaneous self‐assembly. To achieve the pre‐synthesis programmable resolution of Fe(II) grid complexes, we designed and synthesized two novel intrinsically chiral ligands by appending chiral moieties to the parent ligand. The complexation of these chiral ligands with Fe(II) salt resulted in the formation of enantiomerically pure Fe(II) grid complexes, as unambiguously elucidated by CD and XRD studies. The enantiomeric complexes exhibited similar gradual and half‐complete thermal and photo‐induced SCO characteristics. The good agreement between the experimentally obtained and calculated CD spectra further supports the enantiomeric purity of the complexes and even the magnetic studies. The chiral resolution of Fe(II)‐ [2×2] grid complexes reported in this study, for the first time, might enable the fabrication of magneto‐chiral molecular devices.

Resonance Raman Optical Activity Shows Unusual Structural Sensitivity for Systems in Resonance with Multiple Excited States: Vitamin B12 Case

In this work, cobalamins with different upper axial substituents and a cobalamin derivative with a ring modification were studied using chiroptical spectroscopies, in particular resonance Raman optical activity (RROA), to shed light on the influence of structural modifications on RROA spectra in these strongly chiral systems in resonance with multiple excited states at 532 nm excitation. We have demonstrated that for these unique systems RROA possesses augmented structural specificity, surpassing resonance Raman spectroscopy and enabling at the same time measurement of cobalamins at fairy low concentrations of ∼10^-5 mol dm^-3. The enhanced structural specificity of RROA is a result of bisignate spectra due to resonance via more than one electronic state. The observation of increased structural capability of RROA for cobalamins opens a new perspective for studying chiral properties of other biological systems incorporating d-metal ions.

Ab Initio Study of Circular Dichroism and Circularly Polarized Luminescence of Spin-Allowed and Spin-Forbidden Transitions: From Organic Ketones to Lanthanide Complexes

Complete and restricted active space self-consistent field (CAS-/RAS-SCF) wave function methods are applied for the calculation of circular dichroism (CD) and circularly polarized luminescence (CPL) of a series of molecules comprising four organic ketones, the chiral cobalt(III) complex [Co(en)₃]³⁺, and the europium(III) complex [Eu(DPA)₃]^3-. The ab initio results are in good agreement with the experimental data and previous results obtained with Kohn-Sham density functional theory in the case of the spin-allowed transitions. CD and CPL properties are calculated ab initio for the spin-forbidden transitions of both a transition metal and a lanthanide complex.

Resolution-of-identity accelerated relativistic two- and four-component electron dynamics approach to chiroptical spectroscopies

We present an implementation and application of electron dynamics based on real-time time-dependent density functional theory (RT-TDDFT) and relativistic 2-component X2C and 4-component Dirac-Coulomb (4c) Hamiltonians to the calculation of electron circular dichroism and optical rotatory dispersion spectra. In addition, the resolution-of-identity approximation for the Coulomb term (RI-J) is introduced into RT-TDDFT and formulated entirely in terms of complex quaternion algebra. The proposed methodology was assessed on the dimethylchalcogenirane series, C₄H₈X (X = O, S, Se, Te, Po, Lv), and the spectra obtained by non-relativistic and relativistic methods start to disagree for Se and Te, while dramatic differences are observed for Po and Lv. The X2C approach, even in its simplest one-particle form, reproduces the reference 4c results surprisingly well across the entire series while offering an 8-fold speed-up of the simulations. An overall acceleration of RT-TDDFT by means of X2C and RI-J increases with system size and approaches a factor of almost 25 when compared to the full 4c treatment, without compromising the accuracy of the final spectra. These results suggest that one-particle X2C electron dynamics with RI-J acceleration is an attractive method for the calculation of chiroptical spectra in the valence region.

Electronic circular dichroism of chiral alkenes: B3LYP and CAM-B3LYP calculations

Time-dependent density functional theory (TDDFT) calculations of electronic circular dichroism (ECD) are widely used to determine absolute configurations (ACs) of chiral molecules. Two very popular DFT exchange-correlation functionals, one hybrid (B3LYP) and one long-range corrected (CAM-B3LYP), along with a hierarchical sequence of basis sets were investigated, and the ECD spectra predicted for eight alkenes and compared to gas-phase experimental spectra. Little variation in predicted ECD spectra was found with the basis set size enlargement, but the sensitivity to the functional is greater. Good agreement was obtained only with the CAM-B3LYP functional, leading to the conclusion that TDDFT calculations of ECD spectra can routinely provide reliable ACs if and only if an appropriate functional is used. For camphene, twistene, syn-(E)-bisfenchylidene, and phyllocladene, solvent effects were estimated.

Chiroptical spectra of a series of tetrakis((+)-3-heptafluorobutylyrylcamphorato)lanthanide(III) with an encapsulated alkali metal ion: circularly polarized luminescence and absolute chiral structures for the Eu(III) and Sm(III) complexes

The luminescence and circularly polarized luminescence (CPL) spectra of M(I)[Eu((+)-hfbc)(4)] show a similar behavior to the exciton CD in the intraligand π-π* transitions when the alkali metal ions and solvents are manipulated. There is a difference in susceptibility in solvation toward the alkali metal ions but not toward the Eu(III) ion, as in the case of axially symmetric DOTA-type compounds. The remarkable CPL in the 4f-4f transitions provide much more information on the stereospecific formation of chiral Eu(III) complexes, since CPL spectroscopy is limited to luminescent species and reflects selectively toward helicity of the local structural environment around the lanthanide(III). While in comparison, exciton CD reveals the chiral structural information from the helical arrangement of the four bladed chelates. Of special importance, the observation of the highest CPL activities measured to date for lanthanide(III)-containing compounds (i.e., Eu and Sm) in solution supports the theory that the chirality of lanthanide(III) in the excited state corresponds to that in the ground state, which was derived from the exciton CD.

Calculation of optical rotatory dispersion and electronic circular dichroism for tris-bidentate groups 8 and 9 metal complexes, with emphasis on exciton coupling

The optical rotatory dispersion (ORD, both non-resonant and resonant) and the electronic circular dichroism (CD) of tris-bidentate transition metal complexes of the form [M(L)(3)](n+) (M = Fe, Ru, Os, Co, Rh, Ir; n = 2, 3; L = 1,10-phenanthroline, 2,2'-bipyridine) are calculated using time-dependent density functional theory (TDDFT). The exciton CD band resulting from the coupling of ligand π-to-π* transitions is investigated in detail and analyzed in terms of exciton coupling of long-axis transitions using a dipole coupling model that takes TDDFT data for a single ligand as input. Results of the coupling model agree well with the full TDDFT CD spectra. The usefulness and reliability of this model is discussed. The resonant ORDs calculated directly from analytical damped linear TDDFT response compare well with Kramers-Kronig transformations of the calculated CD spectra. For comparisons of resonant ORD with experiment, one needs to consider wavelength shifts.

Alkali metals in ethylenediamine: a computational study of the optical absorption spectra and NMR parameters of [M(en)3(δ+)·M(δ-)] ion pairs

The optical absorption spectra of alkali metals in ethylenediamine have provided evidence for a third oxidation state, -1, of all of the alkali metals heavier than lithium. Experimentally determined NMR parameters have supported this interpretation, further indicating that whereas Na(-) is a genuine metal anion, the interaction of the alkali anion with the medium becomes progressively stronger for the larger metals. Herein, first-principles computations based upon density functional theory are carried out on various species which may be present in solutions composed of alkali metals and ethylenediamine. The energies of a number of hypothetical reactions computed with a continuum solvation model indicate that neither free metal anions, M(-), nor solvated electrons are the most stable species. Instead, [Li(en)(3)](2) and [M(en)(3)(δ+)·M(δ-)] (M = Na, K, Rb, Cs) are predicted to have enhanced stability. The M(en)(3) complexes can be viewed as superalkalis or expanded alkalis, ones in which the valence electron density is pulled out to a greater extent than in the alkali metals alone. The computed optical absorption spectra and NMR parameters of the [Li(en)(3)](2) superalkali dimer and the [M(en)(3)(δ+)·M(δ-)] superalkali-alkali mixed dimers are in good agreement with the aforementioned experimental results, providing further evidence that these may be the dominant species in solution. The latter can also be thought of as an ion pair formed from an alkali metal anion (M(-)) and solvated cation (M(en)(3)(+)).

Excitation and circular dichroism spectra of (+)-(S,S)-bis(2-methylbutyl)chalcogenides

Theoretical electronic spectra and natural circular dichroism (CD) spectra of (+)-(S,S)-bis(2-methylbutyl)chalcogenides, Ch[CH₂CH(CH₃)C₂H₅]₂ (Ch = S, Se, and Te), were calculated by the symmetry adapted cluster (SAC) and SAC-configuration interaction (SAC-CI) methods. Whereas the calculated CD spectrum for each stable conformation itself did not reproduce the corresponding experimental one, their Boltzmann-averaged spectra showed good agreement with the experimental results. We provided the assignment for each spectral band according to our calculation results. For the telluride compound, temperature dependence of the CD spectra was experimentally observed due to variation in the Boltzmann factor, and our calculations reproduced it qualitatively. The spectral features that we could not reproduce can be attributed to triplet transitions through the spin-orbit interaction effects as well as accuracy incompleteness on the calculation conditions.

Time-dependent density functional response theory for electronic chiroptical properties of chiral molecules

Methodology to calculate electronic chiroptical properties from time-dependent density functional theory (TDDFT) is outlined. Applications of TDDFT to computations of electronic circular dichroism, optical rotation, and optical rotatory dispersion are reviewed. Emphasis is put on publications from 2005 to 2010, but much of the older literature is also cited and discussed. The determination of the absolute configuration of chiral molecules by combined measurements and computations is an important application of TDDFT chiroptical methods and discussed in some detail. Raman optical activity (ROA) spectra are obtained from normal-mode derivatives of the optical rotation tensor and other linear response tensors. A few selected (ROA) benchmarks are reviewed.

Racemic atropisomeric N,N-chelate ligands for recognizing chiral carboxylates via Zn(II) coordination: structure, fluorescence, and circular dichroism

Two racemic atropisomeric N,N'-chelate ligands, bis{3,3'-[N-Ph-2-(2'-py)indolyl]} (1) and bis{3,3'-N-4-[N-2-(2'-py)indolyl]phenyl-2-(2'-py)indolyl} (2), have been found to be able to distinguish the enantiomers of Zn((R)-BrMeBu)2 and Zn((S)-BrMeBu)2 where BrMeBu = O2CCH(Br)CHMe2, with a distinct and intense CD spectral response at approximately the 10 microM concentration range. Computational studies established that the (R)-1-Zn((R)-BrMeBu)2 or (S)-1-Zn((S)-BrMeBu)2 diastereomer is more stable than (R)-1-Zn((S)-BrMeBu)2 or (S)-1-Zn((R)-BrMeBu)2. In addition, computational studies showed that the CD spectra of (S)-1-Zn((S)-BrMeBu)2 and (S)-1-Zn((R)-BrMeBu)2 are similar. (1)H NMR spectra confirmed that these two diastereomers exist in solution in about a 2:1 ratio for both complexes of 1 and 2. The distinct CD response of the racemic ligands 1 and 2 toward the chiral zinc(II) carboxylate is therefore attributed to the preferential formation of one diastereomer. The binding modes of the zinc(II) salt with ligands 1 and 2 were established by the crystal structures of the model compounds 1-Zn(tfa)2 and 2-Zn(tfa)2 (tfa = CF3CO2(-)), where the Zn(II) ion is chelated by the two central pyridyl groups in the ligand. Fluorescent titration experiments with various zinc(II) salts showed that the fluorescent spectrum of the atropisomeric ligand displays an anion-dependent change. The zinc(II) binding strength to the N,N'-chelate site of the atropisomeric ligand has been found to play a key role in the selective recognition of different chiral zinc(II) carboxylate derivatives by the racemic atropisomeric ligands.