[Cu6(NGuaS)6]2+ and its oxidized and reduced derivatives: Confining electrons on a torus

The hexanuclear thioguanidine mixed‐valent copper complex cation [Cu₆(NGuaS)₆]⁺² (NGuaS = o‐SC₆H₄NC(NMe₂)₂) and its oxidized/reduced states are theoretically analyzed by means of density functional theory (DFT) (TPSSh + D3BJ/def2‐TZV (p)). A detailed bonding analysis using overlap populations is performed. We find that a delocalized Cu‐based ring orbital serves as an acceptor for donated S p electrons. The formed fully delocalized orbitals give rise to a confined electron cloud within the Cu₆S₆ cage which becomes larger on reduction. The resulting strong electrostatic repulsion might prevent the fully reduced state. Experimental UV/Vis spectra are explained using time‐dependent density functional theory (TD‐DFT) and analyzed with a natural transition orbital analysis. The spectra are dominated by MLCTs within the Cu₆S₆ core over a wide range but LMCTs are also found. The experimental redshift of the reduced low energy absorption band can be explained by the clustering of the frontier orbitals. © 2017 Wiley Periodicals, Inc.     

Time‐dependent density functional calculations on the electronic spectra of the neutral nickel complex [Ni(LISQ)2] (LISQ = 3,5‐di‐tert‐butyl‐o‐diiminobenzosemiquinonate(1−)) and its monoanion and dication

The electronic spectrum of the neutral nickel complex [Ni(L^ISQ)₂] (L^ISQ = 3,5‐di‐tert‐butyl‐o‐diiminobenzosemiquinonate(1^?)) and the spectra of its anion and dication have been calculated by means of time‐dependent density functional theory. The electronic ground state of the neutral complex exhibits an open shell singlet diradical character. The mandatory multireference problem for this electronic ground state has been treated approximately by using the unrestricted and spin symmetry broken Kohn‐Sham Slater determinant as the wave function for the noninteracting reference system in the time‐dependent density functional calculations. A reasonable agreement with observed transition energies and band intensities has been achieved. This holds also for the long wavelength transitions that are shown to be of charge transfer type. The charge distributions in the electronic ground state and the corresponding low lying excited states, however, are rather similar. Thus, the known failure of standard time‐dependent density functional theory to describe improperly long range charge transfer transitions is absent in this work. The applied computational scheme might be adequate for calculating electronic spectra of transition metal complexes with noninnocent ligands. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Electronic,bonding, and optical properties of 1d [CuCN]n (n = 1–10) chains, 2d [CuCN]n (n = 2–10) nanorings,and 3d [Cun(CN)n]m (n = 4, m = 2, 3; n = 10, m = 2) tubes studied by DFT/TD‐DFT methods

The electronic, bonding, and photophysical properties of one‐dimensional [CuCN]_n (n = 1–10) chains, 2‐D [CuCN]_n (n = 2–10) nanorings, and 3‐D [Cu_n(CN)_n]_m (n = 4, m = 2, 3; n = 10, m = 2) tubes are investigated by means of a multitude of computational methodologies using density functional theory (DFT) and time‐dependent‐density‐functional theory (TD‐DFT) methods. The calculations revealed that the 2‐D [CuCN]_n (n = 2–10) nanorings are more stable than the respective 1‐D [CuCN]_n (n = 2–10) linear chains. The 2‐D [CuCN]_n (n = 2–10) nanorings are predicted to form 3‐D [Cu_n(CN)_n]_m (n = 4, m = 2, 3; n = 10, m = 2) tubes supported by weak stacking interactions, which are clearly visualized as broad regions in real space by the 3D plots of the reduced density gradient. The bonding mechanism in the 1‐D [CuCN]_n (n = 1–10) chains, 2‐D [CuCN]_n (n = 2–10) nanorings, and 3‐D [Cu_n(CN)_n]_m (n = 4, m = 2, 3; n = 10, m = 2) tubes are easily recognized by a multitude of electronic structure calculation approaches. Particular emphasis was given on the photophysical properties (absorption and emission spectra) of the [CuCN]_n chains, nanorings, and tubes which were simulated by TD‐DFT calculations. The absorption and emission bands in the simulated TD‐DFT absorption and emission spectra have thoroughly been analyzed and assignments of the contributing principal electronic transitions associated to individual excitations have been made. © 2015 Wiley Periodicals, Inc.     

Electronic structure and spectroscopic properties of 6‐aminophenanthridine and its derivatives: Insights from density functional theory

6‐Aminophenanthridine (6AP) and its derivatives show important biological activities as antiprion compounds and inhibitors of the protein folding activity of the ribosome. Both of these activities depend on the RNA binding property of these compounds, which has been recently characterized by fluorescence spectroscopy. Hence, fundamental insights into the photophysical properties of 6AP compounds are highly important to understand their biological activities. In this work, we have calculated electronic structures and optical properties of 6AP and its three derivatives 6AP8CF₃, 6AP8Cl, and 6APi by density functional theory (DFT) and time‐dependent density functional theory (TDDFT). Our calculated spectra show a good agreement with the experimental absorption and fluorescence spectra, and thus, provide deep insights into the optical properties of the compounds. Furthermore, comparing the results obtained with four different hybrid functionals, we demonstrate that the accuracy of the functionals varies in the order B3LYP > PBE0 > M062X > M06HF. © 2015 Wiley Periodicals, Inc.     

Exploring odd-even effects of simple oligomer-like DRCNnT series: a study based on density functional theory/time-dependent density functional theory calculations

Odd-even effects of short-circuit current density and power conversion efficiency (PCE) are an interesting phenomenon in some organic solar cells. Although some explanations have been given, why they behave in such a way is still an open question. In the present work, we investigate a set of acceptor-donor-acceptor simple oligomer-like small molecules, named the DRCNnT (n = 5-9) series, to give an insight into this phenomenon because the solar cells based on them have high PCE (up to 10.08%) and show strong odd-even effects in experiments. By modeling the DRCNnT series and using density functional theory, we have studied the ground-state electronic structures of the DRCNnT (n = 5-9) series in condensed phase. The calculated results reproduce the experimental trends well. Furthermore, we find that the exciton-binding energies of the DRCNnT series may be one of the key parameters to explain this phenomenon because they also show odd-even effects. In addition, by studying the effects of alkyl branch and terminal group on odd-even effects of dipole moment, we find that eliminating one or two alkyl branches does not break the odd-even effects of dipole moments, but eliminating one or two terminal groups does. Finally, we conclude that removing one alkyl branch close to the terminal group of DRCN5T can decrease highest occupied molecular orbital (HOMO) energy (thus increasing open circuit voltage) and increase dipole moment (thus enhancing charge separation and short-circuit current). This could be a new and simple method to increase the PCE of DRCN5T-based solar cells.     

The σ‐Aromatic Clusters [Zn3]+ and [Zn2Cu]: Embryonic Brass

The triangular clusters [Zn₃Cp*₃]⁺ and [Zn₂CuCp*₃] were obtained by addition of the in situ generated, electrophilic, and isolobal species [ZnCp*]⁺ and [CuCp*] to Carmona’s compound, [Cp*Zn? ZnCp*], without splitting the Zn? Zn bond. The choice of non‐coordinating fluoroaromatic solvents was crucial. The bonding situations of the all‐hydrocarbon‐ligand‐protected clusters were investigated by quantum chemical calculations revealing a high degree of σ‐aromaticity similar to the triatomic hydrogen ion [H₃]⁺. The new species serve as molecular building units of Cu_nZn_m nanobrass clusters as indicated by LIFDI mass spectrometry.     

Time‐dependent local‐density approximation in real time: Application to conjugated molecules

The time‐dependent local‐density approximation (TDLDA) is applied to the optical response of conjugated carbon molecules in the energy range of 0–30 eV, with calculations given for carbon chains, polyenes, retinal, benzene, and C₆₀. The major feature of the spectra, the collective π–π* transition, is seen at energies ranging from below 2 to 7 eV and is reproduced by the theory to a few tenths of an electron volt with a good account of systematic trends. However, there is some indication that TDLDA predicts too much fragmentation of the strength function in large molecules. Transition strengths are reproduced with a typical accuracy of 20%. The theory also predicts a broad absorption peak in the range of 15–25 eV, and this feature agrees with experiment in the one case where quantitative data is available (benzene). ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 55–66, 1999     

An open‐source framework for analyzing N‐electron dynamics. II. Hybrid density functional theory/configuration interaction methodology

In this contribution, we extend our framework for analyzing and visualizing correlated many‐electron dynamics to non‐variational, highly scalable electronic structure method. Specifically, an explicitly time‐dependent electronic wave packet is written as a linear combination of N‐electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time‐dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open‐source Python program det CI@ORBKIT, which extends the capabilities of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one‐electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser‐driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher‐level method provided a judicious choice of functional is made. Broadband excitation of a medium‐sized organic chromophore further demonstrates the scalability of the method. In addition, the time‐dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance. © 2017 Wiley Periodicals, Inc.     

Interaction of Cu(+) and Cu(2+) ions with alpha-alanine. A density functional study

The Cu(+) and Cu(2+) preferred binding sites on alpha-alanine and their affinity values for this amino acid were determined at the density functional level using three different hybrid exchange correlation potentials and the 6-311++G** basis set. The results demonstrated that the two ions both give stable complexes with alpha-alanine but the stability order of the metalated species and the coordination sites are different depending on the nature of the cation. In particular, the Cu(+)-alpha-alanine ground-state structure is characterized by an eta(2)-N,O coordination with the nitrogen and oxygen atoms belonging to the amino and carbonyl groups, respectively. In contrast, the most stable complex of the Cu(2+)-alpha-alanine system has an eta(2)-O,O coordination with the cation bonded to the -CO(2) (-) moiety of the zwitterionic form of the amino acid. Comparison with the Cu(+) and Cu(2+) affinity values for glycine, computed at the same levels of theory, demonstrated that the relative values do not change significantly as different hybrid functionals are used, although the absolute affinities are strongly influenced by the choice of the hybrid potential.     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

Diselenophosphate‐Induced Conversion of an Achiral [Cu20H11{S2P(OiPr)2}9] into a Chiral [Cu20H11{Se2P(OiPr)2}9] Polyhydrido Nanocluster

A polyhydrido copper nanocluster, [Cu₂₀H₁₁{Se₂P(OiPr)₂}₉] ( 2_H ), which exhibits an intrinsically chiral inorganic core of C₃ symmetry, was synthesized from achiral [Cu₂₀H₁₁{S₂P(OiPr)₂}₉] ( 1_H ) of C_3h symmetry by a ligand‐exchange method. The structure has a distorted cuboctahedral Cu₁₃ core, two triangular faces of which are capped along the C₃ axis, one by a Cu₆ cupola and the other by a single Cu atom. The Cu₂₀ framework is further stabilized by 9 diselenophosphate and 11 hydride ligands. The number of hydride, phosphorus, and selenium resonances and their splitting patterns in multinuclear NMR spectra of 2_H indicate that the chiral Cu₂₀H₁₁ core retains its C₃ symmetry in solution. The 11 hydride ligands were located by neutron diffraction experiments and shown to be capping μ₃‐H and interstitial μ₅‐H ligands (in square‐pyramidal and trigonal‐bipyramidal cavities), as supported by DFT calculations on [Cu₂₀H₁₁(Se₂PH₂)₉] ( 2_H′ ) as a simplified model.     

Organic-inorganic hybrids based on novel bimolecular [Si2W22Cu2O78(H2O)]12- polyoxometalates and the polynuclear complex cations [Cu(ac)(phen)(H2O)]n n+ (n=2, 3)

The reaction of a monosubstituted Keggin polyoxometalate (POM) generated in situ with copper-phenanthroline complexes in excess ammonium or rubidium acetate led to the formation of the hybrid metal organic-inorganic compounds A7[Cu2(ac)2(phen)2(H2O)2][Cu3(ac)3(phen)3(H2O)3][Si2W22Cu2O78(H2O)].approximately 18 H2O (A=NH4+ (1), Rb+ (2); ac=acetate; phen=1,10-phenanthroline). These compounds are constructed from inorganic and metalorganic interpenetrated sublattices containing the novel bimolecular Keggin POM, [Si2W22Cu2O78(H2O)]12-, and Cu-ac-phen complexes, [Cu(ac)(phen)(H2O)]n n+ (n=2, 3). The packing of compound 1 can be viewed as a stacking of open-framework layers parallel to the xy plane built of hydrogen-bonded POMs, and zigzag columns of pi-stacked Cu-ac-phen complex cations running along the [111] direction. Magnetic and EPR results are discussed with respect to the crystal structure of the compounds. DFT calculations on [Cu(ac)(phen)(H2O)]n n+ cationic complexes have been performed, to check the influence of packing in the complex geometry and determine the magnetic exchange pathways.     

Theoretical Study of Origin of Triple Fluorescences of a Squaraine Dye, Bis[4-(N,N-dimethylamino)phenyl]-squaraine

A new scheme of photo‐fluorescent emission origin, described as S₀ (relaxed state)→S_n (Frank‐Condon state)→ S_m (relaxed state)→S₀ (Frank‐Condon state), is presented to explain the multiple fluorescent emissions of squaraine dyes observed experimentally according to the configuration interaction singles calculations of relaxed excited states of a model compound, bis[4‐(N,N‐dimethylamino)phenyl]squaraine (SQ). It is exhibited that all triple fluorescent emissions of SQ have their significant origin in vertical electron transitions of different relaxed excited states. In addition, some important absorption peaks appearing in higher energy region are most likely to be responsible for the higher energy band observed in solid states of many squaraine dyes.     

Theoretical study on photooxidation mechanism of ruthenium complex [Ru(II)‐(bpy)2(TMBiimH2)]2+ with molecular oxygen

Photoinduced reactions of ruthenium complexes with molecular oxygen have attracted a lot of experimental attention; however, the reaction mechanism remains elusive. In this work, we have used the density functional theory method to scrutinize the visible‐light induced photooxidation mechanism of the ruthenium complex [Ru(II)‐(bpy)₂(TMBiimH₂)]²⁺ (bpy: 2, 2‐bipyridine and TMBiimH₂: 4, 5, 4, 5‐tetramethyl‐2, 2‐biimidazole) initiated by the attack of molecular oxygen. The present computational results not only explain very well recent experiments, also provide new mechanistic insights. We found that: (1) the triplet energy transfer process between the triplet molecular oxygen and the metal‐ligand charge transfer triplet state of the ruthenium complex, which leads to singlet molecular oxygen, is thermodynamically favorable; (2) the singlet oxygen addition process to the S₀ ruthenium complex is facile in energy; (3) the chemical transformation from endoperoxide to epidioxetane intermediates can be either two‐ or one‐step reaction (the latter is energetically favored). These findings contribute important mechanistic information to photooxidation reactions of ruthenium complexes with molecular oxygen. © 2016 Wiley Periodicals, Inc.     

Excited states of OsO4: A comprehensive time‐dependent relativistic density functional theory study

A large number of scalar as well as spinor excited states of OsO₄, in the experimentally accessible energy range of 3–11 eV, have been captured by time‐dependent relativistic density functional linear response theory based on an exact two‐component Hamiltonian resulting from the symmetrized elimination of the small component. The results are grossly in good agreement with those by the singles and doubles coupled‐cluster linear response theory in conjunction with relativistic effective core potentials. The simulated‐excitation spectrum is also in line with the available experiment. Furthermore, combined with detailed analysis of the excited states, the nature of the observed optical transitions is clearly elucidated. It is found that a few scalar states of ³T₁ and ³T₂ symmetries are split significantly by the spin‐orbit coupling. The possible source for the substantial spin‐orbit splittings of ligand molecular orbitals is carefully examined, leading to a new interpretation on the primary valence photoelectron ionization spectrum of OsO₄. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010