Computational studies on ethylene addition to nickel bis(dithiolene)

The density functionals B3LYP, B3PW91, BMK, HSE06, LC-ωPBE, M05, M06, O3LYP, TPSS, ω-B97X, and ω-B97XD are used to optimize key transition states and intermediates for ethylene addition to Ni(edt)(2) (edt = S(2)C(2)H(2)). The efficacy of the basis sets 6-31G**, 6-31++G**, cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ is also examined. The geometric parameters optimized with different basis sets and density functionals are similar and agree well with experimental values. The ω-B97XD functional gives relative energies closest to those from CCSD, while M06 and HSE06 yield results close to those from CCSD(T). CASSCF and CASSCF-PT2 calculation results are also given. Variation of the relative energies from different density functionals appears to arise, in part, from the multireference character of this system, as confirmed by the T1 diagnostic and CASSCF calculations.     

Predictive power of long-range corrected functionals on the spectroscopic properties of tetrapyrrole derivatives for photodynamic therapy

Porphyrin and chlorin based compounds possess promising properties to be utilized as photosensitizers in photodynamic therapy (PDT). However, the photosensitizers available on the market today are not ideal for use in PDT, which has emphasized the need for new photosensitizers with improved photodynamic properties to be developed. Computational drug-design can be utilized in the search for improved pharmaceutical compounds, provided that the methods used are able to reproduce experimental data. In the present study we investigated, by the use of time-dependent density functional theory (TD-DFT), the performance of the long-range corrected functionals ωB97, ωB97X and ωB97XD on their ability to predict low-lying singlet excitations (>600 nm) of a set of well-known photosensitizing compounds. It was found that ωB97X reproduced the experimental red-most absorption band most satisfactorily. The use of either B3LYP, ωB97XD or M06 in geometry optimizations has a minor effect on the spectra in most cases. Calculated energy differences between the optimized singlet ground states and optimized first excited triplet states show consistent and overall higher triplet state energies for B3LYP, M06, and PBE0 compared with ωB97, ωB97X, and ωB97XD. The calculated triplet state energies are, however, sufficient to generate singlet oxygen in most cases.     

Conformational equilibrium study of calix[4]tetrolarenes using Density Functional Theory (DFT) and Molecular Dynamics simulations

The present work deals with the theoretical study of conformational equilibrium of calix[4]tetrolarenes (1,2,3 trimethoxy substituted calix[4]arene) in gaseous and solvation phases. A total of 64 calculations (four conformations at eight level of theories) were performed using Density Functional Theory (DFT) functionals viz. wB97XD, B97D, B3LYP, CAM-B3LYP with diffused (6–31 + G(d)) and non-diffused basis sets (6-31G(d)). It has been found that the conformational profile of calix[4]tetrolarene changes from cone to 1,3-Alternate upon incorporating – OMe group. The B97D and wB97XD functionals gave the most accurate result having rmsd value ~0.73Å followed by B3LYP and CAM-B3LYP > 0.78Å. Furthermore, NBO calculations demonstrated that reduction in charges at lower rim oxygens reduces the chances of hydrogen bonding. Moreover, global reactivity parameters and molecular dynamics also complement the observed trend. Owing to the methoxy substitutions, anion binding study of these new molecules indicates towards the promising capability to bind Cl^? and F^? ions.     

High first-hyperpolarizabilities of thiobarbituric acid derivative-based donor-π-acceptor nonlinear optical-phores: Multiple theoretical investigations of substituents and conjugated bridges effect

A number of compounds including thiobarbituric acid derivative-based acceptor unit and diverse donor moieties along with two types of π-conjugated bridges have been constructed to aid in the design of the superior nonlinear optical (NLO) materials. The effect of varying the donor and bridge parts of this donor-bridge-acceptor system was analyzed in terms of structural and opto-electronic parameters such as bond length alternation and ultraviolet-visible absorption spectra. Various functionals with aug-cc-pVDZ basis set including B3LYP, PBE0, PBE38, BMK, CAM-B3LYP, and ωB97XD were employed to calculate the static and dynamic first-hyperpolarizabilities, and also the linear polarizabilities. Furthermore, the variation of the static first-hyperpolarizabilities has been explained satisfactorily in the light of the sum-over-states method and two-level model. The comprehensive study indicates that the coplanar compound D-5 with the low electronic absorption energy, strongest oscillator strength paired with the largest ground, and excited state dipole moment difference generates a dramatical increase in its static and dynamic first-hyperpolarizabilities, which would be the most worthwhile target for development as the NLO-phores.     

The relative energies of polypeptide conformers predicted by linear scaling second-order Møller-Plesset perturbation theory

We describe an implementation of the cluster-in-molecule(CIM) resolution of the identity(RI) approximation second-order M?ller–Plesset perturbation theory(CIM-RI-MP2), with the purpose of extending RI-MP2 calculations to very large systems. For typical conformers of several large polypeptides, we calculated their conformational energy differences with the CIM-RI-MP2 and the generalized energy-based fragmentation MP2(GEBF-MP2) methods, and compared these results with the density functional theory(DFT) results obtained with several popular functionals. Our calculations show that the conformational energy differences obtained with CIM-RI-MP2 and GEBF-MP2 are very close to each other. In comparison with the GEBF-MP2 and CIM-RI-MP2 relative energies, we found that the DFT functionals(CAM-B3LYP-D3, LC-?PBE-D3, M05-2X, M06-2X and ?B97XD) can give quite accurate conformational energy differences for structurally similar conformers, but provide less-accurate results for structurally very different conformers.     

A Quantum Chemical Study of the Ground and Excited State Electronic Structures of Carbazole Oligomers with and without Triarylborane Substitutes

Recent experimental investigation (Reitzenstein and Lambert,Macromolecules, 2009, 42, 773) indicated that the quite different optical properties of 2,7- and 3,6-linkage triarylboryl carbazole oligomers may arise from the different nature of their low-lying excited states: a low-lying delocalized within-backbone excitation in longer 2,7-linked oligomers vs a backbone-to-sidechain charge-transfer (CT) excitation independent of the polymerization length in 3,6-linked oligomers. Here in this paper, two long-range corrected functionals, CAM-B3LYP and ωB97X, are applied together with the traditional B3LYP functional in time-dependent density functional theory (TDDFT) calculations to systematically investigate the low-lying electronic excitations in both oligomers. Our calculations indicate that an extensive conjugation exists between monomer molecular orbitals in 2,7-linked oligomers, which is absent in those of 3,6-linked structures, resulting in a considerable narrowing of the HOMO-LUMO gap of their backbone moiety, while having little effect on the side-chains. CAM-B3LYP and ωB97x calculations confirm that the lowest-energy absorption is a within-backbone excitation in longer 2,7-linked oligomers as opposed to a backbone to side-chain charge transfer excitation in 2,7-linked oligmers of shorter length and 3,6-linked oligomers of any length. All these findings are consistent with the experimental findings and the qualitative energy diagram proposed by Reitzenstein and Lambert.     

Diarylethene-dihydroazulene multimode photochrome: a theoretical spectroscopic investigation

Using Time-Dependent Density Functional Theory, we have assessed the structural, energetic and spectroscopic properties of a hybrid diarylethene-dihydroazulene multi-addressable switch synthesized by Mrozek et al. [Chem. Commun., 1999, 1487-1488]. All possible closed/open structures have been considered and the different activation barriers along each reaction path have been evaluated. The determination of the energetic profiles allowed us to pinpoint the thermally possible reactions. To simulate the electronic absorption spectra of this compact two-way four-state hybrid compound, we relied on a PCM-TD-DFT approach combined with a molecular orbital analysis. Key insights are reached and the experimental photochromic properties of this multi-component switch are rationalised. From a more methodological point of view, this work also shows that range-separated hybrid functionals (CAM-B3LYP and ωB97XD) allow to reproduce the measured spectroscopic features with a remarkable accuracy.     

TD-CI simulation of the electronic optical response of molecules in intense fields II: comparison of DFT functionals and EOM-CCSD

Time-dependent configuration interaction (TD-CI) simulations can be used to simulate molecules in intense laser fields. TD-CI calculations use the excitation energies and transition dipoles calculated in the absence of a field. The EOM-CCSD method provides a good estimate of the field-free excited states but is rather expensive. Linear-response time-dependent density functional theory (TD-DFT) is an inexpensive alternative for computing the field-free excitation energies and transition dipoles needed for TD-CI simulations. Linear-response TD-DFT calculations were carried out with standard functionals (B3LYP, BH&HLYP, HSE2PBE (HSE03), BLYP, PBE, PW91, and TPSS) and long-range corrected functionals (LC-ωPBE, ωB97XD, CAM-B3LYP, LC-BLYP, LC-PBE, LC-PW91, and LC-TPSS). These calculations used the 6-31G(d,p) basis set augmented with three sets of diffuse sp functions on each heavy atom. Butadiene was employed as a test case, and 500 excited states were calculated with each functional. Standard functionals yield average excitation energies that are significantly lower than the EOM-CC, while long-range corrected functionals tend to produce average excitation energies slightly higher. Long-range corrected functionals also yield transition dipoles that are somewhat larger than EOM-CC on average. The TD-CI simulations were carried out with a three-cycle Gaussian pulse (ω = 0.06 au, 760 nm) with intensities up to 1.26 × 10(14) W cm(-2) directed along the vector connecting the end carbons. The nonlinear response as indicated by the residual populations of the excited states after the pulse is far too large with standard functionals, primarily because the excitation energies are too low. The LC-ωPBE, LC-PBE, LC-PW91, and LC-TPSS long-range corrected functionals produce responses comparable to EOM-CC.     

A DFT study of substituent effects in corannulene dimers

Corannulene dimers made up of corannulene monomers with different curvature and substituents were studied using M06-2X, B97D and ωB97XD functionals and 6-31+G* basis set. Corannulene molecules were substituted with five alternating Br, Cl, CH(3), C(2)H or CN units. Geometric results showed that substituents gave rise to small changes in the curvature of corannulene bowls. So, there was not a clear relationship between the curvature of bowls and the changes on interaction energy generated by addition of substituents in the bowl. Electron withdrawing substituents gave rise to a more positive molecular electrostatic potential (MEP) of the bowl, which was able to get a strong interaction with the negative MEP at the surface of a fullerene. Substitution with CN caused the largest effect, giving rise to the most positive MEP and to a large interaction energy of -24.64 kcal mol(-1), at the ωB97XD/6-31+G* level. Dispersive effects must be taken into account to explain the catching ability of the different substituted corannulenes. For unsubstituted dimers, calculations with DFT-D methods employing ωB97XD and B97D functionals led to similar results to those previously reported at the SCS-MP2/cc-pVTZ level for corannulene dimers (A. Sygula and S. Saeb?, Int. J. Quant. Chem., 2009, 109, 65). In particular, the ωB97XD functional led to a difference of only 0.35 kcal mol(-1), regarding MP2 interaction energy for corannulene dimers. On the other hand, the M06-2X functional showed a general considerable underestimation of interaction energies. This functional worked quite well to study trends, but not to obtain absolute interaction energies.     

A mixed experimental and theoretical study on chelidamate copper(II) complex with 4-methylpyrimidine

A new chelidamate complex, [Cu(chel)(H₂O)₂(mpd)] (chel = chelidamate; mpd = 4-methylpyrimidine), has been synthesized and characterized through a combination of single crystal X-ray analysis, electron paramagnetic resonance (EPR), ultraviolet-visible (UV-vis), and fourier transform infrared spectroscopy (FT-IR). The complex has six-coordinate distorted octahedral geometry around Cu(II). The theoretical vibrational frequencies and optimized geometric parameters (bond lengths and angles) have been calculated using Density Functional Theory (DFT)/B3LYP and Hartree Fock quantum chemical methods with 6-31G(d, p) basis set by Gaussian 09W software. The EPR spectrum of the compound showed that the paramagnetic center has rhombic symmetry. The EPR studies were carried out using the following unrestricted hybrid density functionals: B3LYP, CAM-B3LYP, HSEH1PBE, WB97XD, MPW1PW91, and BPV86. The UV–vis absorption spectra have been examined in different media and compared with the calculated one using TD-DFT method by applying the polarizable continuum model. Natural bond orbital property of complex has been performed by DFT/B3LYP with 6-31G (d, p) basis set.     

Electronic structure and absorption spectra of supramolecular complexes of a fullerene crown ether with a π-extended TTF derivative

The present work is a theoretical investigation on supramolecular complexes of a fullerene crown ether (A and B isomers) with a derivative of π-extended tetrathiafulvalene (T). The geometry and the electronic structure of seven different conformers of the complex of dibenzo-18-crown-6 ether of fullero-N-methylpyrrolidine with a N-benzyl-N-(4-{[9,10-bis(1,3-dithiol-2-ylidene)-9,10-dihydroanthracen-2-yl]ethynyl}benzyl)ammonium cation were determined. We calculated the complexation energies and the absorption spectra, i.e., the lowest 50 excited electronic states of the complexes have been determined at the ground state optimum geometry. All calculations were carried out employing the density functional theory (DFT) and the time-dependent DFT, using the B3LYP, CAM-B3LYP, ωB97X-D, and M06-2X functionals in conjunction with the 6-31G(d,p) basis set. Various types of van der Waals interactions are observed in the complexes. Conformer complexation energies (CE) range from 2.54 to 2.14 eV in the gas phase and from 1.75 to 1.34 eV in CHCl(3) solvent at the ωB97X-D/6-31G(d,p)//M06-2X/6-31G(d,p) level of theory. There are three major features at about 390, 330, and 290 nm in the calculated absorption spectra of all the conformers. The major peaks correspond to T→T, T→T/F (electron density in both T and the fullerene F of B) and to T→F transitions, depending on the particular conformer. Other charge transfer T→F transitions are observed close to the T→T transition, indicating the possibility of photoinduced electron transfer in all these complexes.     

How accurate are TD-DFT excited-state geometries compared to DFT ground-state geometries?

In this work, we take a different angle to the benchmarking of time-dependent density functional theory (TD-DFT) for the calculation of excited-state geometries by extensively assessing how accurate such geometries are compared to ground-state geometries calculated with ordinary DFT. To this end, we consider 20 medium-sized aromatic organic compounds whose lowest singlet excited states are ideally suited for TD-DFT modeling and are very well described by the approximate coupled-cluster singles and doubles (CC2) method, and then use this method and six different density functionals (BP86, B3LYP, PBE0, M06-2X, CAM-B3LYP, and ωB97XD) to optimize the corresponding ground- and excited-state geometries. The results show that although each hybrid functional reproduces the CC2 excited-state bond lengths very satisfactorily, achieving an overall root mean square error of 0.011 Å for all 336 bonds in the 20 molecules, these errors are distinctly larger than those of only 0.004–0.006 Å with which the hybrid functionals reproduce the CC2 ground-state bond lengths. Furthermore, for each functional employed, the variation in the error relative to CC2 between different molecules is found to be much larger (by at least a factor of 3) for the excited-state geometries than for the ground-state geometries, despite the fact that the molecules/states under investigation have rather uniform chemical and spectroscopic character. Overall, the study finds that even in favorable circumstances, TD-DFT excited-state geometries appear intrinsically and comparatively less accurate than DFT ground-state ones.     

用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.     

Structural and electronic properties of polyacetylene and polyyne from hybrid and Coulomb-attenuated density functionals

The bond length alternation (BLA), the highest-occupied-lowest-unoccupied (HO-LU) orbital energy gap, and the corresponding excitation energy are determined for trans-polyacetylene (PA) and polyyne (PY) using density functional theory. Results from the Coulomb-attenuated CAM-B3LYP functional are compared with those from the conventional BHHLYP and B3LYP hybrid functionals. BLA values and HO-LU gaps are determined using both finite oligomer and infinite chain calculations, subject to periodic boundary conditions. TDDFT excitation energies are determined for the oligomers. The oligomer excitation energies and HO-LU gaps are then used, in conjunction with the infinite chain HO-LU gap, to estimate the infinite chain excitation energy. Overall, BHHLYP and CAM-B3LYP give BLA values and excitation energies that are larger and more accurate than those obtained using B3LYP. The results highlight the degree to which excitation energies can be approximated using the HO-LU gaps-at the infinite limit, this approximation works well for B3LYP, but not for the other functionals, where the HO-LU gap is significantly larger. The study provides further evidence for the high-quality theoretical predictions that can be obtained from the CAM-B3LYP functional.     

What is the best DFT functional for vibronic calculations? A comparison of the calculated vibronic structure of the S1-S0 transition of phenylacetylene with cavity ringdown band intensities

The sensitivity of vibronic calculations to electronic structure methods and basis sets is explored and compared to accurate relative intensities of the vibrational bands of phenylacetylene in the S(1)(A(1)B(2)) ← S(0)(X(1)A(1)) transition. To provide a better measure of vibrational band intensities, the spectrum was recorded by cavity ringdown absorption spectroscopy up to energies of 2000 cm(-1) above the band origin in a slit jet sample. The sample rotational temperature was estimated to be about 30 K, but the vibrational temperature was higher, permitting the assignment of many vibrational hot bands. The vibronic structure of the electronic transition was simulated using a combination of time-dependent density functional theory (TD-DFT) electronic structure codes, Franck-Condon integral calculations, and a second-order vibronic model developed previously [Johnson, P. M.; Xu, H. F.; Sears, T. J. J. Chem. Phys. 2006, 125, 164331]. The density functional theory (DFT) functionals B3LYP, CAM-B3LYP, and LC-BLYP were explored. The long-range-corrected functionals, CAM-B3LYP and LC-BLYP, produced better values for the equilibrium geometry transition moment, but overemphasized the vibronic coupling for some normal modes, while B3LYP provided better-balanced vibronic coupling but a poor equilibrium transition moment. Enlarging the basis set made very little difference. The cavity ringdown measurements show that earlier intensities derived from resonance-enhanced multiphoton ionization (REMPI) spectra have relative intensity errors.     

Conformations and fluorescence of encapsulated stilbene

Absorption and emission spectra of free and encapsulated stilbene in two different capsules were calculated using the DFT and the TDDFT methodology at the B3LYP, CAM-B3LYP, M06-2X, PBE0, and ωB97X-D/6-31G(d,p) levels of theory. The present work is directed toward the theoretical interpretation of recent experimental results on control of stilbene conformation and fluorescence in capsules [Ams, M. R.; et al. Beilstein J. Org. Chem. 2009, 5, 79]. The results of the calculations are in agreement with experiment and show that fluorescence of trans-stilbene persists in the large cage while it is quenched in the small one. It is found that the geometry of trans-stilbene in the ground as well as in the first excited singlet state is unaffected by encapsulation in the large cage, and consequently the absorption and emission spectra are similarly unaffected. In the small cage, the ground state of encapsulated trans-stilbene is distorted, with the two phenyl groups twisted, while the geometry of the excited state, after relaxation, lies at the conical intersection with the ground state. Consequently, there is no emission similar to that of free trans-stilbene, and the state decays nonradiatively to the ground state.     

Halogen bonds with benzene: An assessment of DFT functionals

The performance of an extensive set of density functional theory functionals has been tested against CCSD(T) and MP2 results, extrapolated to the complete basis set (CBS) limit, for the interaction of either DCl or DBr (D = H, HCC, F, and NC) with the aromatic system of benzene. It was found that double hybrid functionals explicitly including dispersion, that is, B2PLYPD and mPW2PLYPD, provide the better agreement with the CCSD(T)/CBS results on both energies and equilibrium geometry, indicating the importance of dispersive contributions in determining this interaction. Among the less expensive functionals, the better performance is provided by the ωB97X and M062X functionals, while the ωB97XD and B97D functionals are shown to work very well for bromine complexes but not so well for chlorine complexes. © 2013 Wiley Periodicals, Inc.     

Density functional and spin-orbit ab initio study of CF3Br: molecular properties and electronic curve crossing

Quantum chemical calculations of CF(3)Br and the CF(3) radical are performed using density functional theory (DFT) and time-dependent DFT (TDDFT). Molecular structures, vibrational frequencies, dipole moment, bond dissociation energy, and vertical excitation energies of CF(3)Br are calculated and compared with available experimental results. The performance of six hybrid and five hybrid meta functionals in DFT and TDDFT calculations are evaluated. The ωB97X, B3PW91, and M05-2X functionals give very good results for molecular structures, vibrational frequencies, and vertical excitation energies, respectively. The ωB97X functional calculates well the dipole moment of CF(3)Br. B3LYP, one of the most widely used functionals, does not perform well for calculations of the C-Br bond length, bond dissociation energy, and vertical excitation energies. Potential energy curves of the low-lying excited states of CF(3)Br are obtained using the multiconfigurational spin-orbit ab initio method. The crossing point between 2A(1) and 3E states is located near the C-Br bond length of 2.45 ?. Comparison with CH(3)Br shows that fluorination does not alter the location of the crossing point. The relation between the calculated potential energy curves and recent experimental result is briefly discussed.