A combined semiempirical and DFT computational protocol for studying bioorganometallic complexes: Application to molybdocene–cysteine complexes

Computational methods can help in the design of new bioorganometallic compounds. However, the presence of multihapto or σ/π metal‐ligand bonding still precludes the direct application of either pure molecular mechanics (MM) or hybrid quantum mechanics‐MM methods to study the flexibility of biomolecules in complex with organometallics. Herein, we present a computational protocol aimed to the evaluation of the relative free energies of bioorganometallic compounds, which explores the conformational space by means of Molecular Dynamics simulations using the semiempirical PM6 method coupled with the COnductor‐like Screening MOdel solvation model followed by density functional theory (DFT) calculations including the DFT‐D3 dispersion energy correction on the most stable conformers. This protocol is applied to investigate the complexes formed between cysteine and the molybdocene entity Cp₂Mo²⁺ (Cp?η⁵? C₅H₅), which has received considerable attention due to its potential antitumor activity and chemical reactivity. The calculated structures and free energies are in agreement with those of the experimentally detected molybdocene‐cysteine adducts and allow us to investigate the reaction mechanism for the formation of the most stable 1:1 and 1:2 adducts. © 2013 Wiley Periodicals, Inc.     

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     

Dynamics and structural changes of small water clusters on ionization

Despite utmost importance in understanding water ionization process, reliable theoretical results of structural changes and molecular dynamics (MD) of water clusters on ionization have hardly been reported yet. Here, we investigate the water cations [(H₂O)_{n = 2–6}⁺] with density functional theory (DFT), Möller–Plesset second‐order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The complete basis set limits of interaction energies at the CCSD(T) level are reported, and the geometrical structures, electronic properties, and infrared spectra are investigated. The characteristics of structures and spectra of the water cluster cations reflect the formation of the hydronium cation moiety (H₃O⁺) and the hydroxyl radical. Although most density functionals fail to predict reasonable energetics of the water cations, some functionals are found to be reliable, in reasonable agreement with high‐level ab initio results. To understand the ionization process of water clusters, DFT‐ and MP2‐based Born‐Oppenheimer MD (BOMD) simulations are performed on ionization. On ionization, the water clusters tend to have an Eigen‐like form with the hydronium cation instead of a Zundel‐like form, based on reliable BOMD simulations. For the vertically ionized water hexamer, the relatively stable (H₂O)₅⁺ (5sL4A) cluster tends to form with a detached water molecule (H₂O). © 2013 Wiley Periodicals, Inc.     

Theoretical insights into the absorption and emission properties of blue luminescent copper(I) complexes based on the pyrazolyl‐pyridine ligands

The ground geometrical and electronic structures, charge transfer (CT) behaviors, absorption, and emission properties of the three copper(I) complexes [Cu(pypz)(POP)]⁺ (1) , [Cu(pympz)(POP)]⁺ (2) , and [Cu(pytfmpz)(POP)]⁺ (3) (pypz=1‐(2‐pyridyl)pyrazole, pympz=3‐methyl‐1‐(2‐pyridyl)pyrazole, and pytfmpz=3‐trifluoromethyl‐1‐(2‐pyridyl)pyrazole), have been investigated using density functional theory (DFT) and time‐dependent density functional theory (TD‐DFT). The vertical absorption energies of the all copper(I) complexes are well reproduced by TD‐DFT calculations based on the CT amount calculations. The triplet emission properties of the all copper(I) complexes were correctly evaluated at BMK/LANL2DZ/6‐31G* level of theory. In addition, the thermally activated delayed fluorescence properties of 1–3 were discussed in detail based on the spatial separation of the HOMO and LUMO and vertical excited energies. These theoretical insights should be expected to provide some guides for the design and synthesis of efficient luminescent copper(I) complexes. © 2014 Wiley Periodicals, Inc.     

Relationship between orbital energy gaps and excitation energies for long‐chain systems

The difference between the excitation energies and corresponding orbital energy gaps, the exciton binding energy, is investigated based on time‐dependent (TD) density functional theory (DFT) for long‐chain systems: all‐trans polyacetylenes and linear oligoacenes. The optimized geometries of these systems indicate that bond length alternations significantly depend on long‐range exchange interactions. In TDDFT formalism, the exciton binding energy comes from the two‐electron interactions between occupied and unoccupied orbitals through the Coulomb‐exchange‐correlation integral kernels. TDDFT calculations show that the exciton binding energy is significant when long‐range exchange interactions are involved. Spin‐flip (SF) TDDFT calculations are then carried out to clarify double‐excitation effects in these excitation energies. The calculated SF‐TDDFT results indicate that double‐excitation effects significantly contribute to the excitations of long‐chain systems. The discrepancies between the vertical ionization potential minus electron affinity (IP–EA) values and the HOMO–LUMO excitation energies are also evaluated for the infinitely long polyacetylene and oligoacene using the least‐square fits to estimate the exciton binding energy of infinitely long systems. It is found that long‐range exchange interactions are required to give the exciton binding energy of the infinitely long systems. Consequently, it is concluded that long‐range exchange interactions neglected in many DFT calculations play a crucial role in the exciton binding energies of long‐chain systems, while double‐excitation correlation effects are also significant to hold the energy balance of the excitations. © 2016 Wiley Periodicals, Inc.     

TD‐DFT benchmarks: A review

Time‐Dependent Density Functional Theory (TD‐DFT) has become the most widely‐used theoretical approach to simulate the optical properties of both organic and inorganic molecules. In this contribution, we review TD‐DFT benchmarks that have been performed during the last decade. The aim is often to pinpoint the most accurate or adequate exchange‐correlation functional(s). We present both the different strategies used to assess the functionals and the main results obtained in terms of accuracy. In particular, we discuss both vertical and adiabatic benchmarks and comparisons with both experimental and theoretical reference transition energies. More specific benchmarks (oscillator strengths, excited‐state geometries, dipole moments, vibronic shapes, etc.) are summarized as well. © 2013 Wiley Periodicals, Inc.     

Bis‐μ‐oxo and μ‐η2:η2‐peroxo dicopper complexes studied within (time‐dependent) density‐functional and many‐body perturbation theory

Based on the equilibrium geometries of [Cu₂(dbdmed)₂O₂]²⁺ and [Cu₂(en)₂O₂]²⁺ obtained within density‐functional theory, we investigate their molecular electronic structure and optical response. Thereby results from occupation‐constrained as well as time‐dependent DFT (ΔSCF and TDDFT) are compared with Green's function‐based approaches within many‐body perturbation theory such as the GW approximation (GWA) to the quasiparticle energies and the Bethe‐Salpeter equation (BSE) approach to the optical absorption. Concerning the ground‐state energies and geometries, no clear trend with respect to the amount of exact exchange in the DFT calculations is found, and a strong dependence on the basis sets is to be noted. They affect the energy difference between bis‐μ‐oxo and μ‐η²:η²‐peroxo complexes by as much as 0.8 eV (18 kcal/mol). Even stronger, up to 5 eV is the influence of the exchange‐correlation functional on the gap values obtained from the Kohn‐Sham eigenvalues. Not only the value itself but also the trends observed upon the bis‐μ‐oxo to μ‐η²:η²‐peroxo transition are affected. In contrast, excitation energies obtained from ΔSCF and TDDFT are comparatively robust with respect to the details of the calculations. Noteworthy, in particular, is the near quantitative agreement between TDDFT and GWA+BSE for the optical spectra of [Cu₂(en)₂O₂]²⁺. © 2013 Wiley Periodicals, Inc.     

Bis(N‐Confused Porphyrin) as a Semirigid Receptor with a Chirality Memory: A Two‐Way Host Enantiomerization through Point‐to‐Axial Chirality Transfer

The adduct formation of protonated bis(N‐confused porphyrin) (BNCP, 3,3′‐bis(meso‐tetratolyl‐2‐aza‐21‐carbaporphyrin) with chiral anions, carboxylic acids, and alcohols was studied in solution by means of ¹H NMR and circular dichroism (CD) spectroscopic analysis and DFT methods. The addition of enantiopure guests to the acidified BNCP resulted in optical activity that vanished after neutralization. Pairs of the ¹H NMR‐distinguishable diastereomers were formed when enantiopure guests were applied, although a single form was observed upon the addition of the racemic mixtures in each case. Unidirectional configuration change that led to diastereomeric excess was observed in several instances. Such an excess was memorized by metalation of the adducts with AgBF₄, thus resulting in optically active silver(III) complexes of BNCP with some enantiomeric excess. Absolute configurations of BNCP cations and bis(zinc) and bis(silver(III)) complexes were determined on the basis of time‐dependent (TD)‐DFT calculations of their CD spectra. It was shown that some of the chiral carboxylates induced opposite directions of enantiomerization of di‐ and tetracations or di‐/tetracation and bis(zinc) complexes. The source of the optical activity of the equimolar diastereomeric mixture of adducts is discussed.     

Calculating core‐level excitations and x‐ray absorption spectra of medium‐sized closed‐shell molecules with the algebraic‐diagrammatic construction scheme for the polarization propagator

Core‐level excitations are generated by absorption of high‐energy radiation such as X‐rays. To describe these energetically high‐lying excited states theoretically, we have implemented a variant of the algebraic‐diagrammatic construction scheme of second‐order ADC(2) by applying the core‐valence separation (CVS) approximation to the ADC(2) working equations. Besides excitation energies, the CVS‐ADC(2) method also provides access to properties of core‐excited states, thereby allowing for the calculation of X‐ray absorption spectra. To demonstrate the potential of our implementation of CVS‐ADC(2), we have chosen medium‐sized molecules as examples that have either biological importance or find application in organic electronics. The calculated results of CVS‐ADC(2) are compared with standard TD‐DFT/B3LYP values and experimental data. In particular, the extended variant, CVS‐ADC(2)‐x, provides the most accurate results, and the agreement between the calculated values and experiment is remarkable. © 2014 Wiley Periodicals, Inc.     

Chemically modified fullerene derivatives as photosensitizers in photodynamic therapy: A first‐principles study

The first‐principles density functional theory (DFT) and its time‐dependent approach (TD‐DFT) are used to characterize the electronic structures and optical spectra properties of five chemically modified fullerenes. It is revealed that the metal fullerene derivatives possess not only stronger absorption bands in visible light regions than organically modified fullerene but also the large energy gaps (ΔE_S–T > 0.98 eV) between the singlet ground state and the triplet state, which imply their significant aspect of potential candidates as a photosensitizer. We have found that a new metal‐containing bisfullerene complexes (Pt(C₆₀)₂), with the extended conjugated π‐electrons, much degenerate orbitals and a uniform electrostatic potential surface, behave more pre‐eminent photosensitizing properties than other examined fullerene derivatives. © 2012 Wiley Periodicals, Inc.     

Impact of replacement of the central benzene ring in anthracene by a heterocyclic ring on electronic excitations and reorganization energies in anthratetrathiophene molecules

The impact of changing the central benzene ring on the electronic excitations and reorganization energies (λ) of the anthratetrathiophene (ATT) molecules is studied by density functional theory (DFT) and time‐dependent DFT (TD‐DFT) quantum chemical calculations. The effect of changing the position of the sulfur atom at the periphery of anthracene on the optical and charge transfer properties is also studied. The calculated results suggest that the HOMO, LUMO, HOMO–LUMO energy gap, ionization potential (IP), electron affinity (EA), hole extraction potential (HEP), electron extraction potential (EEP), and reorganization energies (λ) are affected by replacing the central ring with different heterocyclic rings and the position of the sulfur atom. In addition, all molecules show good hole‐ and electron‐transport properties. This work may be helpful for future design and preparation of high‐performance charge‐transport materials.     

Electronic Properties of Two Series of Blue-Emitting Conjugated Oligomers and Polymers Based on Cyanostilbene： A Quantum Chemistry Study

One of the drawbacks of the electroluminescence （EL） polymers is that they are usually much better at accepting and transporting holes than electrons due to their inherent richness of π-electrons. One approach improving electron injection and transport in conjugated polymers is to incorporate moieties with high electron affinities. In this theoretical work, to gain an insight into the chemical structure-property relationships was aimed by controllable modification of the main chain structures. Two cyanovinylene derivatives with 2,7-fluorenylene and p-phenylene moieties in the main chains, namely, poly { （2,5-dimethoxy-p-phenylene- 1,4-ylene）-alt-[ 1,2-bis（p-phenylene）- 1- cyanovinylene]} （PPhCN） and poly{[9,9-dimethyl-2,7-fluorenylene]-alt-[1,2-bis（p-phenylene）-1-cyanovinylene]} （PFCN）, were studied employing density functional theory （DFT） and time dependent density functional theory （TD-DFT） with B3LYP functional. The electronic properties of the neutral molecules, extrapolated ionization potentials （IP） and electron affinities （EA）, and energy gaps were investigated in comparison with pristine poly（2,7- fluorenylene）. From comparison with poly（2,7-fluorenylene） （PF）, the 1,2-bis（p-phenylene）-1-cyanovi-nylene unit was found to be a good electron-withdrawing moiety for electronic materials and the incorporation of 1,2-bis（p- phenylene）-1-cyanovinylene resulted in a narrow band gap and a red shift of both the absorption and photoluminescence emission peaks. Most importantly, the LUMO energies of PFCN are around 1 eV lower than those of PF, which results in the decrement of EA about 0.9 eV, indicating that the 1,2-bis（p-phenylene）-1-cyanovinylene unit has significantly improved the electron-accepting properties of the copolymer PFCN. Substitution of 2,5-dimethoxy-p-phenylene for 9,9-dimethyl-2,7-fluorenylene induced larger band gaps and thus a blue-shift in absorption and emission peaks, which can be attributed to the better conjugated backbone in PFCN.     

GALAMOST: GPU‐accelerated large‐scale molecular simulation toolkit

GALAMOST [graphics processing unit (GPU)‐accelerated large‐scale molecular simulation toolkit] is a molecular simulation package designed to utilize the computational power of GPUs. Besides the common features of molecular dynamics (MD) packages, it is developed specially for the studies of self‐assembly, phase transition, and other properties of polymeric systems at mesoscopic scale by using some lately developed simulation techniques. To accelerate the simulations, GALAMOST contains a hybrid particle‐field MD technique where particle–particle interactions are replaced by interactions of particles with density fields. Moreover, the numerical potential obtained by bottom‐up coarse‐graining methods can be implemented in simulations with GALAMOST. By combining these force fields and particle‐density coupling method in GALAMOST, the simulations for polymers can be performed with very large system sizes over long simulation time. In addition, GALAMOST encompasses two specific models, that is, a soft anisotropic particle model and a chain‐growth polymerization model, by which the hierarchical self‐assembly of soft anisotropic particles and the problems related to polymerization can be studied, respectively. The optimized algorithms implemented on the GPU, package characteristics, and benchmarks of GALAMOST are reported in detail. © 2013 Wiley Periodicals, Inc.     

β‐Iminoenamine‐BF2 Complexes: Aggregation‐Induced Emission and Pronounced Effects of Aliphatic Rings on Radiationless Deactivation

The synthesis, photophysical, and electrochemical attributes of a novel class of boron difluorides containing an aromatic‐fused alicyclic/hetero‐alicyclic ring built on a β‐iminoenamine chromophoric backbone are reported. The compounds displayed large Stokes shifts (86–121 nm), and were emissive in the solid state. The quantum yields obtained in solution at room temperature were unusually lower by an order of magnitude compared to those in the solid state. Some of the tested compounds displayed aggregation‐induced emission (AIE). Single crystal XRD analyses revealed a lack of interplanar π–π interactions, which are presumed to be absent owing to non‐planarity of the alicyclic component in the molecule. For most of the studied compounds, time‐dependent DFT (TD‐DFT) calculations invariably reveal intramolecular charge transfer (π–π*) characteristics with the frontier orbitals concentrated on the boron–nitrogen heterocycle. The participation of boron and fluorine atoms was found to be negligible.