Excited state properties,fluorescence energies,and lifetimes of a poly(fluorene‐phenylene), based on TD‐DFT investigation

The structural and electronic properties of fluorene‐phenylene copolymer (FP)_n, n = 1–4 were studied by means of quantum chemical calculations based on density functional theory (DFT) and time dependent density functional theory (TD‐DFT) using B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest singlet excited state. It was found that (FP)_n is nonplanar in its ground state while the electronic excitations lead to planarity in its S₁ state. Absorption and fluorescence energies were calculated using TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods. Vertical excitation energies and fluorescence energies were obtained by extrapolating these values to infinite chain length, resulting in extrapolated values for vertical excitation energy of 2.89 and 2.87 eV, respectively. The S₁ ← S₀ electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is distinguishing in terms of oscillator strength. Fluorescence energies of (FP)_n calculated from TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods are 2.27 and 2.26 eV, respectively. Radiative lifetimes are predicted to be 0.55 and 0.51 ns for TD‐B3LYP/SVP and TD‐B3LYP/SVP+ calculations, respectively. These fundamental information are valuable data in designing and making of promising materials for LED materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Structures and excitation energies of Zn–tetraarylporphyrin analogues: A theoretical study

The photophysical properties of β-substituted Zn–tetraarylporphyrin (ZnTAP) analogues used as dyes in dye-sensitized solar cells were studied using density functional theory (DFT). Singlet-excitation energy calculations of ZnTAP analogues were performed using time-dependent DFT with B3LYP, B3PW91, PBE0 exchange–correlation functionals at 6-31G(d) and 6-31+G(d) basis sets using B3LYP/6-31G(d) geometries. The PBE0 functional at 6-31+G(d) basis set provided a better correlation with the experimental data for both B- and Q-bands. The inclusion of solvation effect in the calculations provided a good agreement in terms of B:Q_ave ratio of the oscillator strengths for both analogues with the experimental values. Analogue 2 has a higher and a more balanced charge-carrier transport rates than analogue 1. In general, the addition of an electron-donating group in the meso-substituent (analogue 2) resulted in a narrower band gap, higher oscillator strength, a more red-shifted absorption spectra, and better charge-transfer characteristics than analogue 1.     

One‐bond 29Si‐1H spin‐spin coupling constants in the series of halosilanes: benchmark SOPPA and DFT calculations,relativistic effects,and vibrational corrections

A number of most representative second order polarization propagator approach (SOPPA) based wavefunction methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), and density functional theory (DFT) based methods, B3LYP, PBE0, KT2, and KT3, have been benchmarked in the calculation of the one‐bond ²⁹Si‐¹H spin‐spin coupling constants in the series of halosilanes SiH_nX_4?n (X = F, Cl, Br, I), both at the non‐relativistic and full four‐parameter Dirac's relativistic levels taking into account vibrational corrections. At the non‐relativistic level, the wavefunction methods showed much better results as compared with those of DFT. At the DFT level, out of four tested functionals, the Perdew, Burke, and Ernzerhof's PBE0 showed best performance. Taking into account, relativistic effects and vibrational corrections noticeably improves wavefunction methods results, but generally worsens DFT results. Copyright © 2013 John Wiley & Sons, Ltd.     

Combined experimental and computational modeling studies on 3,5‐dimethyl‐pyrazole‐1‐carbodithioic acid benzyl ester

The title compound, 3,5‐Dimethyl‐pyrazole‐1‐carbodithioic acid benzyl ester, has been synthesized and structurally characterized by X‐ray single crystal diffraction, elemental analysis, IR spectra, and UV‐Vis spectrum. The crystal belongs to orthorhombic, space group P212121, with a = 5.3829(15), b = 11.193(3), c = 21.824(6) Å, V = 1315.0(6) ų, and Z = 4. The molecules are connected via intermolecular C–H···N hydrogen bonds into 1D infinite chains. The crystal structure is consolidated by the intramolecular C–H···S hydrogen bonds. Furthermore, Density functional theory (DFT) calculations of the structure, stabilities, orbital energies, composition characteristics of some frontier molecular orbitals and Mulliken charge distributions of the title compound were performed by means of Gaussian 03W package and taking B3LYP/6‐31G(d) basis set. The time‐dependent DFT (TD‐DFT) calculations have been employed to calculate the electronic spectrum of the title compound, and the UV‐Vis spectra has been discussed on this basis. The results show that DFT method at B3LYP/6‐31G(d) level can well reproduce the structure of the title compound. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Assessment of dispersion‐improved exchange‐correlation functionals for the simulation of CO2 binding by alcoholamines

In this study, 12 bound complexes were selected to construct a database for testing 15 dispersion‐improved exchange‐correlation (XC) functionals, including hybrid generalized gradient approximation (GGA), modified using the Grimme's pairwise strategy, and double hybrid XC functionals, for specifically characterizing the CO₂ binding by alcoholamines. Bound complexes were selected based on the characteristics of their hydrogen bonds, dispersion, and electrostatic (particularly between the positive charge of CO₂ and the lone pair of N of alcoholamines) interactions. The extrapolated binding energy from the aug‐cc‐pVTZ (ATZ) to aug‐cc‐pVQZ (AQZ) basis set at the CCSD(T)/CBS(MP2+DZ) level was used as the reference for the XC functional comparison. M06‐2X produced the optimal agreement if the optimized geometries at MP2/ATZ level were chosen for all the test bound complexes. However, M06‐L, ωB97X, and ωB97, and were preferred if the corresponding density functional theory (DFT) optimized geometries were adapted for the benchmark. Simple bimolecular reaction between CO₂ and monoethanolamine simulated using polarizable continuum solvation model confirmed that ωB97, ωB97X, and ωB97XD qualitatively reproduced the energetics of MP2 level. The inconsistent performance of the tested XC functionals, observed when using MP2 or DFT optimized geometries, raised concerns regarding using the single‐point ab initio correction combined with DFT optimized geometry, particularly for determining the nucleophilic attack by alcoholamines to CO₂. © 2014 Wiley Periodicals, Inc.     

Theoretical study of the magnetic exchange coupling behavior substituting Cr(III) with Mo(III) in cyano‐bridged transition metal complexes

Molecular magnetism in a series of cyano‐bridged first and second transition metal complexes has been investigated using density functional theory (DFT) combined with the broken‐symmetry (BS) approach. Several exchange‐correlation (XC) functionals in the ADF package were used to investigate complexes I [?(Me₃tacn)₂(cyclam)NiMo₂(CN)₆]²⁺, II [?(Me₃tacn)₂(cyclam)Ni‐Cr₂(CN)₆]²⁺, III [(Me₃tacn)₆MnMo₆(CN)₁₈]²⁺, and IV [(Me₃tacn)₆MnCr₆(CN)₁₈]²⁺ (Me₃tacn = N,N′,N?‐trimethyl‐1,4,7‐triazacyclononane). For models A (the molded structure of complex I) and B (the modeled structure of complex II), all the XCs given qualitatively reasonable results and predict ferromagnetic coupling character between M (M = Mo^III for A or Cr^III for B) and Ni^II in coincidence with the experimental results (see Tables I and II ). The calculated using Operdew, OPBE, O3LYP, and B3LYP functionals and experimental J values show that substituting Cr^III with Mo^III will enhance the ferromagnetic exchange coupling interactions. But VWN, PW91, PBE, VSXC, and tau‐HCTH functionals have no way to differentiate the relative strength of the intramolecular magnetic exchange coupling interactions of A and B correctly. For models C (the modeled structure of complex III) and D (the modeled structure of complex IV), all the XCs in ADF and B3LYP in Gaussian 03 with several basis sets show that substituting Cr^III with Mo^III will enhance the antiferromagnetic exchange coupling interactions. From the above calculations, the substitution of Cr^III by Mo^III will enhance the magnetic coupling interactions, whether the magnetic coupling interactions are ferro‐ or antiferromagnetic. Moreover, Kahn's model was applied to investigate the above facts. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09

A benchmark study on all possible density functional theory (DFT) methods in Gaussian09 is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (E_int) for small non‐covalently interacting molecular dimers in “dispersion‐dominated” (class 1), “dipole‐induced dipole” (class 2), and “dipole‐dipole” (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and E_int was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good E_int, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional “recommended” for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the “recommended” list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and E_int. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in Gaussian09 to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.     

A DFT study on Cu(I) coordination in Cu‐ZSM‐5: Effects of the functional choice and tuning of the ONIOM approach

The coordination of Cu⁺ at the T1 and T7 positions of the M7 ring of Cu‐ZSM‐5, and the interaction of NO with coordinated Cu⁺ were investigated by means of DFT/ONIOM calculations. The B3LYP, BLYP, PBE1PBE, PBE, M06, and M062X functionals with the def2‐TZVP (def2‐QZVP for Cu) basis set were used in the high‐level part of ONIOM calculations, with the HF/3‐21G, B3LYP/LANL2DZ, M06/LANL2DZ, and M062X/LANL2DZ methods in the low‐level part. The ability of suitable combinations of the above methods to reproduce (i) the crystallographic structure of purely siliceous ZSM‐5, (ii) the tendency of Cu⁺ to be twofold or fourfold coordinated by framework oxygen atoms of Cu‐ZSM‐5, and (iii) the interaction energy and the N? O stretching frequency of adsorbed nitrogen oxide are discussed, showing that different results are obtained depending on the adopted computational approach. With reference to above properties, some considerations about the employment of the ONIOM approximations are also included. © 2015 Wiley Periodicals, Inc.     

Accurate prediction of vertical electronic transitions of Ni(II) coordination compounds via time dependent density functional theory

Time dependent density functional theory calculations are completed for five Ni(II) complexes formed by polydentate peptides to predict the electronic absorption spectrum. The ligands examined were glycyl‐glycyl‐glycine (GGG), glycyl‐glycyl‐glycyl‐glycine (GGGG), glycyl‐glycyl‐histidine (GGH), glycyl‐glycyl‐cysteine (GGC), and triethylenetetramine (trien). Fifteen functionals and two basis sets were tested. On the basis of the mean absolute percent deviation (MAPD), the ranking among the functionals is: HSE06 ∼ MPW1PW91 ∼ PBE0 > ω‐B97x‐D ∼ B3P86 ∼ B3LYP ∼ CAM‐B3LYP > PBE ∼ BLYP ∼ BP86 > TPSS > TPSSh > BHandHLYP > M06 ≫ M06‐2X. Concerning the basis sets, the triple‐ζ def2‐TZVP performs better than the double‐ζ LANL2DZ. With the functional HSE06 and basis set def2‐TZVP the MAPD with respect to the experimental λ_max is 1.65% with a standard deviation of 1.26%. The absorption electronic spectra were interpreted in terms of vertical excitations between occupied and virtual MOs based on Ni‐d atomic orbitals. The electronic structure of the Ni(II) species is also discussed.     

Theoretical Simulation of Vibrational Sum‐Frequency Generation Spectra from Density Functional Theory: Application to p‐Nitrothiophenol and 2,4‐Dinitroaniline

The molecular orientation of adsorbed molecules forming self‐assembled monolayers can be determined by combining vibrational sum‐frequency generation (SFG) measurements with quantum chemical calculations. Herein, we present a theoretical methodology used to simulate the SFG spectra for different combinations of polarizations. These simulations are based on calculations of the IR vectors and Raman tensors, which are obtained from density functional theory computations. The dependency of the SFG vibrational signature with respect to the molecular orientation is presented for the molecules p‐nitrothiophenol and 2,4‐dinitroaniline. It is found that a suitable choice of basis set as well as of exchange‐correlation (XC) functional is mandatory to correctly simulate the SFG intensities and consequently provide an accurate estimation of the adsorbed molecule orientation. Comparison with experimental data shows that calculations performed at the B3LYP/6‐311++G(d,p) level of approximation provide good agreement with experimental frequencies, and with IR and Raman intensities. In particular, it is demonstrated that polarization and diffuse functions are compulsory for reproducing the IR and Raman spectra, and consequently vibrational SFG spectra, of systems such as p‐nitrothiophenol. Moreover, the investigated XC functionals reveal their influence on the relative intensities, which show rather systematic variations with the amount of Hartree–Fock exchange. Finally, further aspects of the modeling are revealed by considering the frequency dependence of the Raman tensors.     

Thiocyanate‐free ruthenium(II) tetrabenzoporphyrin sensitizers for photoelectrochemical cell: A DFT/TD‐DFT probe for stability of axial donor ligands

Although previously studied [(HOOC)₄(TBPor)Ru(NCS)₂]^2– ( A ; TBPor = tetrabenzoporphrin) avoided the intrinsic π‐stacking aggregation of planar metallophorphryins via incorporating two axial ligands, these isothiocyanato groups are believed to be the weakest part of the sensitizer while operating in dye‐sensitized solar cells (DSSCs). In this work, a series of thiocyanate‐free ruthenium porphyrin complexes featuring with phenyl/substituted‐phenyl axial groups, [(HOOC)₄(TBPor)Ru(L′)₂]^2– (L′ = Ph ( 1 ), PhF₂ ( 2 ), PhCl₂ ( 3 ), PhBr₂ ( 4 ), and PhI₂ ( 5 )), have been examined using density functional theory (DFT) and time‐dependent DFT (TD‐DFT). Both analyses of electronic structures and calculations of interaction energies demonstrate that the Ru‐L′ interaction in 1 – 5 is significantly enhanced relative to the Ru‐NCS in A , which will raise chemical stability of the former in DSSCs. Single‐electron oxidation mechanism has been proposed. Oxidation potentials (E⁰) are increased by 0.2–0.6 V when changing axial groups from NCS to Ph/PhX₂. The spin‐orbit coupling (SOC) relativistic effects can be negligible in computing E⁰ values. TD‐DFT calculations show that 1 – 5 have more intense Q band in the visible region than A does. Taken together, high chemical stability, suitable oxidation potential and expanding absorption spectra would allow for potential applications of the thiocyanate‐free sensitizers in DSSCs. © 2016 Wiley Periodicals, Inc.     

Theoretical Modelling of Photoswitching of Hyperpolarisabilities in Ruthenium Complexes

Static excited‐state polarisabilities and hyperpolarisabilities of three Ru^II ammine complexes are computed at the density functional theory (DFT) and several correlated ab initio levels. Most accurate modelling of the low energy electronic absorption spectrum is obtained with the hybrid functionals B3LYP, B3P86 or M06 for the complex [Ru^II(NH₃)₅(MeQ⁺)]³⁺ (MeQ⁺=N‐methyl‐4,4′‐bipyridinium, 3 ) in acetonitrile. The match with experimental data is less good for [Ru^II(NH₃)₅L]³⁺ (L=N‐methylpyrazinium, 2 ; N‐methyl‐4‐{E,E‐4‐(4‐pyridyl)buta‐1,3‐dienyl}pyridinium, 4 ). These calculations confirm that the first dipole‐ allowed excited state (FDAES) has metal‐to‐ligand charge‐transfer (MLCT) character. Both the solution and gas‐phase results obtained for 3 by using B3LYP, B3P86 or M06 are very similar to those from restricted active‐space SCF second‐order perturbation theory (RASPT2) with a very large basis set and large active space. However, the time‐dependent DFT λ_max predictions from the long‐range corrected functionals CAM‐B3LYP, LC‐ωPBE and wB97XB and also the fully ab initio resolution of identity approximate coupled‐cluster method (gas‐phase only) are less accurate for all three complexes. The ground state (GS) two‐state approximation first hyperpolarisability β_2SA for 3 from RASPT2 is very close to that derived experimentally via hyper‐Rayleigh scattering, whereas the corresponding DFT‐based values are considerably larger. The β responses calculated by using B3LYP, B3P86 or M06 increase markedly as the π‐conjugation extends on moving along the series 2 → 4 , for both the GS and FDAES species. All three functionals predict substantial FDAES β enhancements for each complex, increasing with the π‐conjugation, up to about sevenfold for 4 . Also, the computed second hyperpolarisabilities γ generally increase in the FDAES, but the results vary between the different functionals.     

Toward the complete range separation of non‐hybrid exchange–correlation functional

In this study, we use a very simple scheme to achieve range separation of a total exchange–correlation functional. We have utilized this methodology to combine a short‐range pure density functional theory (DFT) functional with a corresponding long‐range pure DFT, leading to a “Range‐separated eXchange–Correlation” (RXC) scheme. By examining the performance of a range of standard exchange–correlation functionals for prototypical short‐ and long‐range properties, we have chosen B‐LYP as the short‐range functional and PBE‐B95 as the long‐range counterpart. The results of our testing using a more diverse range of data sets show that, for properties that we deem to be short‐range in nature, the performance of this prescribed RXC‐DFT protocol does resemble that of B‐LYP in most cases, and vice versa. Thus, this RXC‐DFT protocol already provides meaningful numerical results. Furthermore, we envisage that the general RXC scheme can be easily implemented in computational chemistry software packages. This study paves a way for further refinement of such a range‐separation technique for the development of better performing DFT procedures. © 2015 Wiley Periodicals, Inc.     

Theoretical study of the electronic spectra of bi‐ and tri‐heteronuclear platinum complexes

The electronic structure and the spectroscopic properties of [Pt(NH₃)₄][Au(CN)₂]₂, [Pt(NH₃)₄][Ag(CN)₂]₂, [Pt(CNCH₃)₄][Pt(CN)₄], and [Pt(CNCH₃)₄][Pd(CN)₄] were studied at the HF, MP2, B3LYP, and PBE levels. In all the complexes, it was found that the nature of the intermetal interactions is consistent with the presence of a high‐ionic contribution (90%) and a dispersion‐type interaction (10%). The absorption spectra of these complexes were calculated by the single‐excitation time‐dependent (TD) method at the HF, B3LYP, and PBE levels. The [Pt(NH₃)₄][M(CN)₂]₂ (M ? Au, Ag) complexes showed a ¹(dσ* → pσ) transition associated with a metal–metal charge transfer. On the other hand, the [Pt(CNCH₃)₄][M(CN)₄] (M ? Pt, Pd) complexes showed a ¹(dσ* → π*) transition associated with a metal‐to‐metal and ligand charge transfer. The values obtained theoretically are in agreement with the experimental range. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Interaction of 1,3,2,4‐Benzodithiadiazines with Aromatic Phosphines and Phosphites

Although an interaction between hydrocarbon and fluorocarbon 1,3,2,4‐benzodithiadiazines ( 1 ) and P(C₆H₅)₃ continuously produces chiral 1,2,3‐benzodithiadiazol‐2‐yl iminophosporanes ( 2 ; in this work, 5,7‐difluoro derivative 2a ) via 1:1 condensation, an interaction between 1 and other PR₃ reagents gives different products. With R  OC₆H₅ and both hydrocarbon and fluorocarbon 1 , only X=P(OC₆H₅)₃ (X = S, O) were identified in the complex reaction mixtures by ¹³С and ³¹Р NMR and GC‐MS. With R = C₆F₅, no interaction with the archetypal 1 was observed but catalytic addition of atmospheric water to the heterocycle afforded 2‐amino‐N‐sulfinylbenzenesulfenamide ( 4 ). With electrophilic B(C₆F₅)₃ instead of nucleophilic P(C₆F₅)₃, only adduct H₃N→B(C₆F₅)₃ and a new polymorph of C₆F₅B(OH)₂ were isolated and identified by X‐ray diffraction (XRD). A molecular structure of 2a was confirmed by XRD, and the π‐stacked orientation of one of phenyl groups and heterocyclic moiety was observed. This structure is in general agreement with that calculated at the RI‐MP2 level of theory, as well as at three different levels of DFT theory with the PBE and B3LYP functionals. Mild thermolysis of 2a in a dilute decane solution gave persistent 5,7‐difluoro‐1,2,3‐benzodithiazolyl ( 3a ) identified by EPR in combination with DFT calculations.     

Parallel implementation of γ‐point pseudopotential plane‐wave DFT with exact exchange

Semi‐local functionals commonly used in density functional theory (DFT) studies of solids usually fail to reproduce localized states such as trapped holes, polarons, excitons, and solitons. This failure is ascribed to self‐interaction which creates a Coulomb barrier to localization. Pragmatic approaches in which the exchange correlation functionals are augmented with small amount of exact exchange (hybrid‐DFT, e.g., B3LYP and PBE0) have shown to promise in rectifying this type of failure, as well as producing more accurate band gaps and reaction barriers. The evaluation of exact exchange is challenging for large, solid state systems with periodic boundary conditions, especially when plane‐wave basis sets are used. We have developed parallel algorithms for implementing exact exchange into pseudopotential plane‐wave DFT program and we have implemented them in the NWChem program package. The technique developed can readily be employed in Γ‐point plane‐wave DFT programs. Furthermore, atomic forces and stresses are straightforward to implement, making it applicable to both confined and extended systems, as well as to Car‐Parrinello ab initio molecular dynamic simulations. This method has been applied to several systems for which conventional DFT methods do not work well, including calculations for band gaps in oxides and the electronic structure of a charge trapped state in the Fe(II) containing mica, annite. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010