手性配合物[M(bipy)_3]~(2+)(M=Fe,Ru,Os)圆二色谱的理论解析

应用密度泛函理论,在DFT/B3LYP/LanL2DZ水平上优化了铁族配合物[M(bipy)3]2+(M=Fe,Ru,Os;bipy=2,2′-联吡啶)在溶液中的基态几何结构,并用TDDFT/B3LYP方法和相同的基组计算了该类配合物的激发能、旋转强度和振子强度,绘制了相应的圆二色谱.计算的谱带位置虽有一定的红移或蓝移,但各谱带的带形和符号均与实验谱吻合.对跃迁性质的分析表明,三个配合物在长波区,除[Fe(bipy)3]2+的第一个谱带是以d-d跃迁为主外,其他谱带均是荷移跃迁为主.短波区的两个强带则是以π-π*跃迁为主的激子耦合带:对于铁配合物,TDDFT可以正确预言其正负带的强度比,但对钌和锇配合物而言,由于对激子跃迁中荷移成份的贡献估计不足,计算的比值偏小.这一结论为进一步改进有关的计算模型指明了方向.     

Time dependent density functional theory modeling of specific rotation and optical rotatory dispersion of the aromatic amino acids in solution

Time Dependent Density Functional Theory (TDDFT) along with the COnductor-like Screening MOdel (COSMO) has been applied to model the specific rotation at 589.3 nm and the optical rotatory dispersion (ORD) of the aromatic amino acids phenylalanine, tyrosine, histidine, and tryptophan. Solution structures at low, neutral, and high pH were determined. Both the anomalous dispersion absorbing (resonance) region and the lower energy (transparent) region of the ORD of the compounds were modeled. Linear response calculation of the specific rotation and ORD as well as Kramers-Kronig transformations of calculated circular dichroism spectra to model resonant ORD were compared with experimental data from the literature.     

Theory of the absorption and circular dichroism spectra of helical molecular aggregates

A theory of the electronic circular dichroism (CD) and optical rotatory dispersion (ORD) of infinite aggregates exhibiting cylindrical symmetry is presented in which, to the authors' knowledge, for the first time vibrational structure is included explicitly. It is shown that, with the coherent exciton scattering approximation in the Green function approach, the detailed vibrational structure of the aggregate absorption. CD and ORD bands can be calculated from a knowledge of the electronic coupling and the monomer absorption line shape alone. Detailed model calculations for a single helix are made and the results are used to expose the origin of different spectral features. A good reproduction of experimental J-aggregate spectra is obtained, using the same electronic interaction to fit both absorption and CD spectral line shapes. The theory allows some prediction of aggregate geometry to be made, but it is shown that an unambiguous geometrical assignment can only be made where experimental spectra for light of different propagation directions with respect to the cylinder axis are available.     

Density Functional Study of Tetraphenylporphyrin Long-Range Exciton Coupling

The performance of time-dependent density functional theory (TDDFT) for calculations of long-range exciton circular dichroism (CD) is investigated. Tetraphenylporphyrin (TPP) is used as a representative of a class of strongly absorbing chromophores for which exciton CD with chromophore separations of 50 Å and even beyond has been observed experimentally. A dimer model for TPP is set up to reproduce long-range exciton CD previously observed for a brevetoxin derivative. The calculated CD intensity is consistent with TPP separations of over 40 Å. It is found that a hybrid functional with fully long-range corrected range-separated exchange performs best for full TDDFT calculations of the dimer. The range-separation parameter is optimally tuned for TPP, resulting in a good quality TPP absorption spectrum and small DFT delocalization error (measured by the curvature of the energy calculated as a function of fractional electron numbers). Calculated TDDFT data for the absorption spectra of TPP are also used as input for a ‘matrix method’ (MM) model of the exciton CD. For long-range exciton CD, comparison of MM spectra with full TDDFT CD spectra for the dimer shows that the matrix method is capable of producing very accurate results. A MM spectrum obtained from TPP absorption data calculated with the nonhybrid Becke88–Perdew86 (BP) functional is shown to match the experimental brevetoxin spectrum ‘best’, but for the wrong reasons.     

Circular dichroism, optical rotation and absolute configuration of 2-cyclohexenone-cis-diol type phenol metabolites: redefining the role of substituents and 2-cyclohexenone conformation in electronic circular dichroism spectra

The absolute configurations of 2-cyclohexenone cis-diol metabolites resulting from the biotransformation of the corresponding phenols have been determined by comparison of their experimental and calculated circular dichroism spectra (TDDFT at the PCM/B2LYP/Aug-cc-pVTZ level), optical rotations (calculated at the PCM/B3LYP/Aug-cc-pVTZ level) and by stereochemical correlation. It is found that circular dichroism spectra and optical rotations of 2-cyclohexenone derivatives are strongly dependent on the ring conformation (M or P sofa S(5) or half-chair), enone non-planarity and the nature and positions of the hydroxy and alkyl substituents. The effect of non-planarity of the enone chromophore, including the distortion of the C=C bond, is determined for the model structures by TDDFT calculations at the PCM/B2LYP/6-311++G(2d,2p) level. Non-planarity of the C=C bond in the enone chromophore is commonly encountered in 2-cyclohexenone derivatives and it is a source of significant rotatory strength contribution to the electronic circular dichroism spectra. It is shown that the two lowest-energy transitions in acrolein and 2-cyclohexenone and its derivatives are n(C=O)-π(C=O)* and π(C=C)-π(C=O)*, as expected, while the shorter-wavelength (below 200 nm) transitions are of more complex nature. In 2-cyclohexenone and its alkyl derivatives it is predominantly a mixture of π(C=C)-π(C=C)* and π(C=C)-σ* transitions, whereas the presence of hydroxy substituent results in a dominant contribution due to the n(OH)-π(C=O)* transition. A generalized model for correlation of the CD spectra of 2-cyclohexenones with their structures is presented.     

手性络合物[M(bpy)_2(phen)]~(2+)和[M(phen)_2(bpy)]~(2+)(M=Ru,Os)圆二色谱的理论解析

应用密度泛函理论,在B3LYP/LanL2DZ水平上对C2对称性的混配络合物[M(bpy)2(phen)]2+和[M(phen)2(bpy)]2+(M=Ru、Os;bpy=2,2'-bipyridine;phen=1,10-phenanthroline)在水溶液中的几何构型进行了优化,并用TDDFT/B3LYP方法和相同的基组计算了其激发能、旋转强度和振子强度,绘制了相应的圆二色谱(CD).在分析有关跃迁性质的基础上,对实验圆二色谱的谱带进行了明确的解析和指认,同时讨论了短波区激子裂分的规律性.结果表明:四种络合物在长波区(λ>320nm)的CD吸收带主要是由d-π*跃迁产生的荷移谱带;短波区(λ<320nm)则是配体上平行于长轴的π-π*跃迁产生的激子耦合带,且对于Λ构型表现为正的手性激子裂分.其中,[M(bpy)2(phen)]2+只显示出正负两个激子带,分属于联吡啶和邻菲咯啉配体;而[M(phen)2(bpy)]2+则有三个激子带,其中左侧的两个(一负一正)属于邻菲咯啉配体,右侧的正带则属于联吡啶配体.此外,尽管激子耦合属于远程相互作用,但用TDDFT计算的激子裂分样式仍是正确的.这些结论对于深入理解有...     

A density functional theory (DFT) and time-dependent density functional theory (TDDFT) study on optical transitions in oligo(p-phenylenevinylene)-fullerene dyads and the applicability to resonant energy transfer

The optical transitions of three different size oligo(p-phenylenevinylene)-fullerene dyads (OPV(n)-MPC(60); n = 2-4) and of the corresponding separate molecules are studied using density functional theory (DFT) and time-dependent density functional theory. The DFT is used to determine the geometries and the electronic structures of the ground states. Transition energies and excited-state structures are obtained from the TDDFT calculations. Resonant energy transfer from OPV(n) to MPC(60) is also studied and the Fermi golden rule is used, along with two simple models to describe the electronic coupling to calculate the energy transfer rates. The hybrid-type PBE0 functional is used with a split-valence basis set augmented with a polarization function (SV(P)) in calculations and the calculated results are compared to the corresponding experimental results. The calculated PBE0 spectra of the OPV(n)-MPC(60) dyads correspond to the experimental spectra very well and are approximately sums of the absorption spectra of the separate OPV(n) and MPC(60) molecules. Also, the absorption energies of OPV(n) and MPC(60) and the emission energies of OPV(n) are predicted well with the PBE0 functional. The PBE0 calculated resonant energy transfer rates are in a good agreement with the experimental rates and show the existence of many possible pathways for energy transfer from the first excited singlet states of the OPV(n) molecules to the MPC(60) molecule.     

Stereoselective formation of dinuclear complexes with anomalous CD spectra

A novel chiral tetradentate Schiff base ligand forms dinuclear lambda,lambda-[M2L2Cl2]2+ complexes with high diastereoselectivity; the complexes exhibit anomalous CD spectra due to exciton coupling of chromophores located on different metal centres.     

Calculation of circular dichroism spectra from optical rotatory dispersion, and vice versa, as complementary tools for theoretical studies of optical activity using time-dependent density functional theory

A comparison of two theoretical methods based on time-dependent density functional theory for the calculation of the linear dispersive and absorptive properties of chiral molecules has been made. For this purpose, a recently proposed computational method for the calculation of circular dichroism (CD) spectra from the imaginary part of the optical rotation parameter has been applied to six rigid organic molecules. The results have been compared to the CD spectra obtained from the rotatory strengths and from the Kramers-Kronig transformation of optical rotatory dispersion (ORD) curves. We have also investigated a criterion based on the Kramers-Kronig integration formula to determine a number of excitations in truncated CD spectra which may yield a reasonable low frequency resonant ORD. It has been tested by calculating the ORD from the sum-over-states formula both in the nonresonant and resonant regions. Finally, we have applied these methods to model the resonant optical activity of proline at low pH.     

CD spectra of polynuclear complexes of diimine ligands: theoretical and experimental evidence for the importance of internuclear exciton coupling

We have recently reported on dinuclear complexes Lambda,Lambda-[Co(2)L(2)Cl(2)]CoCl(4) of two novel chiral ligands (1a and 1b) which contain pyridyl-imine chelate groups (Telfer, S. G.; Sato, T.; Kuroda, R. Chem. Commun. 2003, 1064-1065). The absolute configuration of the cobalt(II) centers was unambiguously assigned by X-ray crystallography. However, the sign of the exciton couplets in their CD spectra was opposite to that expected on the basis of the stereochemistry of the metal centers. We present a rationalization of these anomalous spectra in terms of an "internuclear" exciton coupling model which takes into account the coupling of chromophores located on different metal centers. We have performed a series of semiempirical (ZINDO) calculations which provide quantitative support to this model. These findings show that the absolute configuration of the metal centers in a polynuclear complex may be incorrectly assigned on the basis of CD data if internuclear coupling effects are not taken into consideration. We summarize the CD spectral data of number of other chiral polynuclear complexes from the literature, including dinuclear complexes bridged by the 2,2'-bipyrimidine ligand, complexes of the HAT ligand, and dinuclear triple-stranded helicates. The amplitude of the CD spectra of many of these complexes is not additive with the number of chromophores. These anomalous spectra can be accounted for by taking internuclear coupling effects into consideration.     

C60-多吡啶Ru(II)配合物电子光谱的密度泛函理论研究

应用密度泛函理论(DFT)方法对两种C60-多吡啶Ru(II)衍生物进行理论研究. 在TZP全电子基组优化构型基础上, 通过分析前线轨道组成, 探讨金属及配体对C60母体影响; 以LB及SAOP校正局域密度近似, 用含时密度泛函(TDDFT)方法, 考虑溶剂化效应, 计算化合物1和2的电子吸收光谱. 结果表明, 化合物1和2在气相与丙酮溶液中所对应的光谱值差异较为明显, 溶剂化效应使吸收光谱蓝移. 计算得到化合物1和2在丙酮溶液中电子光谱与实验值吻合较好, 低能跃迁多为金属参与的混合跃迁, 高能跃迁主要由C60与配体部分贡献.     

Spectroscopic and theoretical studies of optically active porphyrin dimers: a system uninterpretable by exciton coupling theory.

The electronic excited states of a meso-meso beta-beta doubly linked bis-porphyrin are comprehensively investigated by measuring its circular dichroism (CD) and magnetic circular dichroism (MCD) spectra. The observed spectroscopic properties are rationalized by DFT calculations. The frontier molecular orbitals (MOs) are constructed by the linear combinations of the constituent monomers' four MOs. Comparison of a theoretical CD spectrum based on time-dependent DFT (TDDFT) with the experimental spectra resulted in the assignment of the helical conformation of the dimer. This assignment is contrary to the previous assignment based on the point-dipole approximation (exciton coupling theory).     

Synthesis, structure and characterization of fac-[Re(CO)3] complexes derived from hydrazone Schiff bases: DFT-TDDFT investigation on electronic structures

The syntheses, structural and spectroscopic characterization of the complexes of general formula [ReL(CO)₃Cl] bearing bifunctional hydrazone Schiff base ligand L are presented in this paper. The structure of one of the complexes is determined by X-ray crystallography. The solid-state structure of the compound is involved in a secondary interaction in lattice forming a supramolecular array. The gas phase geometry optimization and electronic calculation have been performed using density functional theory without any symmetry constraints. On the basis of structural and theoretical studies, ligand in the complexes is considered to be in the keto, not in enol form. Experimental ground state IR and NMR data set agree with those calculated by DFT calculations. The electronic spectra of the complexes are calculated by time dependent density functional theory (TDDFT) using conductor like polarizable continuum model (CPCM). The computed vertical excitation energies in solution are in good agreement with experimental one showing that the metal-to-ligand charge transfer transitions in visible region dominate over ligand based ILCT transition. The TDDFT excited states calculation of the electronic spectra in solution provides evidence towards luminescence spectra.     

Tailoring transition metal complexes for nonlinear optics applications. 2. A theoretical investigation of the second-order nonlinear optical properties of M(CO)(5)L complexes (M = Cr, W; L = Py, PyCHO, Pyz, PyzBF(3), BPE, BPEBF(3))

In this work, we report an ab initio investigation of second-order nonlinear optical (NLO) properties and absorption electronic spectra of push-pull transition metal chromophores of the formula [M(CO)(5)L] (M = Cr, W; L = pyridine (Py), 4-formyl-pyridine (PyCHO), pyrazine (Pyz), trans-1,2-bis(4-pyridyl)ethylene (BPE)). Pyz and BPE are considered either with one nitrogen atom free or interacting with the strong acceptor BF(3). All of the molecular properties have been calculated using two different and methodologically independent approaches: the time dependent and coupled perturbed density functional theories (TDDFT and CPDFT) and the sum-over-states (SOS) approach, where the excited states are obtained via the single configuration interaction (SCI) ab initio method. DFT results are in acceptable agreement with the experimental energy values of electronic transitions (with the exception of chromophores with the large pi-delocalization, like BPE); SCI calculations overestimate excitation energies and produce an inversion in the order of d(M) --> pi(L) and d(M) --> pi(CO) transitions. The SCI-SOS approach gives first-order hyperpolarizabilities, basically in agreement as trend and values with the experiments and seems to be a tool generally suitable for the evaluation of these properties also for transition metal complexes. On the other hand, the first-order hyperpolarizabilities computed using the CPDFT approach are consistently overestimated in comparison with the experimental results, especially in the case of a ligand with large pi-delocalization. We also show that the "two-level" approximation taking into account only the lowest energy charge transfer excitation (e.g., d(M) --> pi(L)) is not applicable to chromophores with the extended pi-delocalized ligand (BPE) coordinated to a transition metal, due to significant contributions originating from intraligand pi(L) --> pi(L) transitions. This study reports a detailed analysis and comparison of electronic NLO effects of transition metal complexes computed with DFT and ab initio SCI-SOS methodology.     

Theoretical modeling of peptide α-helical circular dichroism in aqueous solution

Reliable modeling of protein and peptide circular dichroism (CD) spectra in the far UV presents a challenge for current theoretical approaches. In this study, the time-dependent density functional theory (TDDFT), configuration interaction with single excitation (CIS), and transition dipole coupling (TDC) were used to assess the most important factors contributing to the CD spectra of the α-helical secondary structure. The dependence on the peptide chain length and also the role of the flexibility and solvent environment were investigated with a model oligopeptide Ac-(Ala)(N)-NH-Me, (N = 1, ..., 18). Both the TDDFT and TDC-like methods suggest that the CD curve typical for the α-helix arises gradually, but its basic characteristic is discernible already for peptides with 4-5 amino acid residues. The calculated dependence was in a qualitative agreement with experimental spectra of short α-helices stabilized by the histidine-metal binding. The TDDFT computations of the CD were found to be unusually sensitive to the basis set and solvent model. Explicit hydration and temperature fluctuations of the peptide geometry, simulated with the aid of molecular dynamics (MD), significantly influenced the CD and absorption spectral shapes. An extensive averaging over MD configurations is thus required to obtain a converged spectral profile in cluster simulations. On the other hand, both the TDDFT and TDC models indicate only a minor influence of the alanine side chains. The CIS and TDC calculations also point toward a relatively small effect of the helix-helix interaction on the CD spectral profiles. For a model system of two helices, the CIS method predicted larger changes in the spectra than TDC. This suggests other than interactions between peptide chains, such as mutual polarization, can have a minor, but measurable, effect on the CD spectrum.     

Inherent chirality dominates the visible/near-ultraviolet CD spectrum of rhodopsin

The visible (alpha) and near-UV (beta) CD bands of rhodopsin have been studied extensively, but their source(s) have never been definitively established. Do they result from the intrinsic chirality of the polyene chromophore of the protonated Schiff base of retinal (retPSB) or from the coupling of the transitions of this chromophore with those of protein groups? We have calculated the contributions of these two mechanisms to the CD of rhodopsin. The intrinsic CD of the retPSB chromophore was calculated using time-dependent density functional theory (TDDFT) and, for comparison, the semiempirical ZINDO method. First-order perturbation theory was used to calculate the effects of coupling of the retPSB transitions with the pi pi* transitions of the aromatic chromophores and the pi pi* and n pi* transitions of the peptide groups in rhodopsin. Calculations were performed for eight structures based upon the two molecules in the asymmetric unit of four crystal structures. The most reliable results were obtained from TDDFT calculations on the structure of Okada et al. (J. Mol. Biol. 2004, 342, 571), PDB 1U19. Averaging over the two molecules in the asymmetric unit, the intrinsic rotational strengths are 0.62 +/- 0.00 DBM (Debye-Bohr magneton) and 0.90 +/- 0.03 DBM for the alpha- and beta-bands, respectively. The contributions from coupling with protein groups are, respectively, -0.32 +/- 0.05 and -0.01 +/- 0.03 DBM. Our results show that the visible/near-UV CD bands of rhodopsin are determined by the intrinsic chirality of the retPSB chromophore and that the contributions of coupling with the protein are significantly smaller for the alpha-band and negligible for the beta-band.     

Photophysical properties of a series of electron-donating and -withdrawing platinum acetylide two-photon chromophores

To explore spectroscopic structure-property relationships in platinum acetylides, we synthesized a series of complexes having the molecular formula trans-bis(tributylphosphine)-bis(4-((9,9-diethyl-7-ethynyl-9H-fluoren-2-yl)ethynyl)-R)-platinum. The substituent, R = NH(2), OCH(3), N(phenyl)(2), t-butyl, CH(3), H, F, benzothiazole, CF(3), CN, and NO(2), was chosen for a systematic variation in electron-donating and -withdrawing properties as described by the Hammett parameter σ(p). UV/vis, fluorescence, and phosphorescence spectra, transient absorption spectra on the fs-ps time scale, and longer time scale flash photolysis on the ns time scale were collected. DFT and TDDFT calculations of the T(1) and S(1) energies were performed. The E(S) and E(T) values measured from linear spectra correlate well with the calculated results, giving evidence for the delocalized MLCT character of the S(1) state and confinement of the T(1) exciton on one ligand. The calculated T(1) state dipole moment ranges from 0.5 to 14 D, showing the polar, charge-transfer character of the T(1) state. The ultrafast absorption spectra have broad absorption bands from 575 to 675 nm and long wavelength contribution, which is shown from flash photolysis measurements to be from the T(1) state. The T(1) energy obtained from phosphorescence, the T(1)-T(n) transition energy obtained from flash photolysis measurements, and the triplet-state radiative rate constant are functions of the calculated spin density distribution on the ligand. The calculations show that the triplet exciton of chromophores with electron-withdrawing substituents is localized away from the central platinum atom, red-shifting the spectra and increasing the triplet-state lifetime. Electron-donating substituents have the opposite effect on the location of the triplet exciton, the spectra, and the triplet-state lifetime. The relation between the intersystem crossing rate constant and the S(1)-T(1) energy gap shows a Marcus relationship with a reorganization energy of 0.83 eV. The calculations show that intersystem crossing occurs by conversion from a nonpolar, delocalized S(1) state to a polar, charge-transfer T(1) state confined to one ligand, accompanied by conformation changes and charge transfer, supporting the experimental evidence for Marcus behavior.     

Structure, bonding, and linear optical properties of a series of silver and gold nanorod clusters: DFT/TDDFT studies

DFT/TDDFT calculations have been carried out for a series of silver and gold nanorod clusters (Ag(n), Au(n), n = 12-120) whose structures are of cigar-type. Pentagonal Ag(n) clusters with n = 49-121 and hexagonal Au(n) clusters with n = 14-74 were also calculated for comparison. Metal-metal distances, binding energies per atom, ionization potentials, and electron affinities were determined, and their trends with cluster size were examined. The TDDFT calculated excitation energies and oscillator strengths were fit by a Lorentz line shape modification, which gives rise to the simulated absorption spectra. The significant features of the experimental spectra for actual silver and gold nanorod particles are well reproduced by the calculations on the clusters. The calculated spectral patterns are also in agreement with previous theoretical results on different-type Ag(n) clusters. Many differences in the calculated properties are found between the Ag(n) and Au(n) clusters, which can be explained by relativistic effects.     

Electronic structure and spectra of nitrosyl complexes with cobalt and manganese porphyrins

The DFT calculations for nitrosyl manganese and cobalt porphyrins were carried out with the use of several density functionals. The binding energy of nitrosyl ligand and spin state of nitrosyl-free manganese porphyrin were determined. The best values of binding energy are obtained from the OLYP functional. The NBO analysis of metal?Cnitrosyl bonding was performed. Electronic spectra of nitrosyl cobalt and manganese porphyrin were calculated with the TDDFT method. The calculated electronic transitions agree well with the experimental data except for the Soret band of (Por)Mn(NO), where they are 0.3?C0.5?eV higher in energy than the experimental ones.     

Electronically excited states and visible region photodissociation spectroscopy of Au(m)(+) . Ar(n) clusters (m=7-9): molecular dimensionality transition?

Photodissociation spectra were determined for Au(m)(+) . Ar(n) (m=7; n=0-3 and m=8,9; n=0,1) in the photon energy range of 2.14-3.02 eV. Experimental data were compared with predictions of dipole allowed transitions using time-dependent density functional theory (TDDFT) as applied to cluster structures from both DFT (B3-LYP functional) and ab initio calculations at the MP2 level. Argon adduct formation does not significantly perturb the bare metal cluster core structure, but it does change the metal cluster spectrum for highly symmetric cluster structures. The photodissociation spectra are consistent with a transition from planar to three-dimensional gold cluster core geometries between m=7 and m=8 for both n=0 and 1. TDDFT predictions for favored isomers describe experimental absorption features to within +/-0.25 eV. We also discuss size-dependent trends in TDDFT transition energies for the lowest energy two- and three-dimensional structures of Au(m)(+)(m=3-9).