Vibrationally resolved photoelectron spectra of lower diamondoids: A time-dependent approach.

Vibrationally resolved lowest-energy bands of the photoelectron spectra (PES) of adamantane, diamantane, and urotropine were simulated by a time-dependent correlation function approach within the harmonic approximation. Geometries and normal modes for neutral and cationic molecules were obtained from B3LYP hybrid density functional theory (DFT). It is shown that the simulated spectra reproduce the experimentally observed vibrational finestructure (or its absence) quite well. Origins of the finestructure are discussed and related to recurrences of autocorrelation functions and dominant vibrations. Remaining quantitative and qualitative errors of the DFT-derived PES spectra refer to (i) an overall redshift by ∼0.5 eV and (ii) the absence of satellites in the high-energy region of the spectra. The former error is shown to be due to the neglect of many-body corrections to ordinary Kohn-Sham methods, while the latter has been argued to be due to electron-nuclear couplings beyond the Born-Oppenheimer approximation [Gali et al., Nat. Commun. 7, 11327 (2016)].

The 1,6,7,12-Tetraazaperylene Bridging Ligand as an Electron Reservoir and Its Disulfonato Derivative as Redox Mediator in an Enzyme-Electrode Process.

The homodinuclear ruthenium(II) complex [{Ru(l-N4 Me2 )}2 (µ-tape)](PF6 )4 {[1](PF6 )4 } (l-N4 Me2 =N,N'-dimethyl-2,11-diaza[3.3](2,6)-pyridinophane, tape=1,6,7,12-tetraazaperylene) can store one or two electrons in the energetically low-lying π* orbital of the bridging ligand tape. The corresponding singly and doubly reduced complexes [{Ru(l-N4 Me2 )}2 (µ-tape.- )](PF6 )3 {[2](PF6 )3 } and [{Ru(l-N4 Me2 )}2 (µ-tape2- )](PF6 )2 {[3](PF6 )2 }, respectively, were electrochemically generated, successfully isolated and fully characterized by single-crystal X-ray crystallography, spectroscopic methods and magnetic susceptibility measurements. The singly reduced complex [2](PF6 )3 contains the π-radical tape.- and the doubly reduced [3](PF6 )2 the diamagnetic dianion tape2- as bridging ligand, respectively. Nucleophilic aromatic substitution at the bridging tape in [1]4+ by two sulfite units gave the complex [{Ru(l-N4 Me2 )}2 {µ-tape-(SO3 )2 }]2+ ([4]2+ ). Complex dication [4]2+ was exploited as a redox mediator between an anaerobic homogenous reaction solution of an enzyme system (sulfite/sulfite oxidase) and the electrode via participation of the low-energy π*-orbital of the disulfonato-substituted bridging ligand tape-(SO3 )22- (Ered1 =-0.1 V versus Ag/AgCl/1 m KCl in water).

Organic electronic materials: recent advances in the DFT description of the ground and excited states using tuned range-separated hybrid functionals.

CONSPECTUS: Density functional theory (DFT) and its time-dependent extension (TD-DFT) are powerful tools enabling the theoretical prediction of the ground- and excited-state properties of organic electronic materials with reasonable accuracy at affordable computational costs. Due to their excellent accuracy-to-numerical-costs ratio, semilocal and global hybrid functionals such as B3LYP have become the workhorse for geometry optimizations and the prediction of vibrational spectra in modern theoretical organic chemistry. Despite the overwhelming success of these out-of-the-box functionals for such applications, the computational treatment of electronic and structural properties that are of particular interest in organic electronic materials sometimes reveals severe and qualitative failures of such functionals. Important examples include the overestimation of conjugation, torsional barriers, and electronic coupling as well as the underestimation of bond-length alternations or excited-state energies in low-band-gap polymers. In this Account, we highlight how these failures can be traced back to the delocalization error inherent to semilocal and global hybrid functionals, which leads to the spurious delocalization of electron densities and an overestimation of conjugation. The delocalization error for systems and functionals of interest can be quantified by allowing for fractional occupation of the highest occupied molecular orbital. It can be minimized by using long-range corrected hybrid functionals and a nonempirical tuning procedure for the range-separation parameter. We then review the benefits and drawbacks of using tuned long-range corrected hybrid functionals for the description of the ground and excited states of π-conjugated systems. In particular, we show that this approach provides for robust and efficient means of characterizing the electronic couplings in organic mixed-valence systems, for the calculation of accurate torsional barriers at the polymer limit, and for the reliable prediction of the optical absorption spectrum of low-band-gap polymers. We also explain why the use of standard, out-of-the-box range-separation parameters is not recommended for the DFT and/or TD-DFT description of the ground and excited states of extended, pi-conjugated systems. Finally, we highlight a severe drawback of tuned range-separated hybrid functionals by discussing the example of the calculation of bond-length alternation in polyacetylene, which leads us to point out the challenges for future developments in this field.

Accurate description of torsion potentials in conjugated polymers using density functionals with reduced self-interaction error.

We investigate the torsion potentials in two prototypical π-conjugated polymers, polyacetylene and polydiacetylene, as a function of chain length using different flavors of density functional theory. Our study provides a quantitative analysis of the delocalization error in standard semilocal and hybrid density functionals and demonstrates how it can influence structural and thermodynamic properties. The delocalization error is quantified by evaluating the many-electron self-interaction error (MESIE) for fractional electron numbers, which allows us to establish a direct connection between the MESIE and the error in the torsion barriers. The use of non-empirically tuned long-range corrected hybrid functionals results in a very significant reduction of the MESIE and leads to an improved description of torsion barrier heights. In addition, we demonstrate how our analysis allows the determination of the effective conjugation length in polyacetylene and polydiacetylene chains.

A highly K(+)-selective phenylaza-[18]crown-6-lariat-ether-based fluoroionophore and its application in the sensing of K+ ions with an optical sensor film and in cells.

Herein, we report the synthesis of two phenylaza-[18]crown-6 lariat ethers with a coumarin fluorophore (1 and 2) and we reveal that compound 1 is an excellent probe for K(+) ions under simulated physiological conditions. The presence of a 2-methoxyethoxy lariat group at the ortho position of the anilino moiety is crucial to the substantially increased stability of compounds 1 and 2 over their lariat-free phenylaza-[18]crown-6 ether analogues. Probe 1 shows a high K(+)/Na(+) selectivity and a 2.5-fold fluorescence enhancement was observed in the presence of 100 mM K(+) ions. A fluorescent membrane sensor, which was prepared by incorporating probe 1 into a hydrogel, showed a fully reversible response, a response time of 150 s, and a signal change of 7.8% per 1 mM K(+) within the range 1-10 mM K(+). The membrane was easily fabricated (only a single sensing layer on a solid polyester support), yet no leaching was observed. Moreover, compound 1 rapidly permeated into cells, was cytocompatible, and was suitable for the fluorescent imaging of K(+) ions on both the extracellular and intracellular levels.

On the relationship between bond-length alternation and many-electron self-interaction error.

Predicting accurate bond-length alternations (BLAs) in long conjugated molecular chains has been a major challenge for electronic-structure theory for many decades. While Hartree-Fock (HF) overestimates BLA significantly, second-order perturbation theory and commonly used density functional theory (DFT) approaches typically underestimate it. Here, we discuss how this failure is related to the many-electron self-interaction error (MSIE), which is inherent to both HF and DFT approaches. We use tuned long-range corrected hybrids to minimize the MSIE for a series of polyenes. The key result is that the minimization of the MSIE alone does not yield accurate BLAs. On the other hand, if the range-separation parameter is tuned to yield accurate BLAs, we obtain a significant MSIE that grows with chain length. Our findings demonstrate that reducing the MSIE is one but not the only important aspect necessary to obtain accurate BLAs from density functional theory.

Collectively induced quantum-confined Stark effect in monolayers of molecules consisting of polar repeating units.

The electronic structure of terpyrimidinethiols is investigated by means of density-functional theory calculations for isolated molecules and monolayers. In the transition from molecule to self-assembled monolayer (SAM), we observe that the band gap is substantially reduced, frontier states increasingly localize on opposite sides of the SAM, and this polarization in several instances is in the direction opposite to the polarization of the overall charge density. This behavior can be analyzed by analogy to inorganic semiconductor quantum-wells, which, as the SAMs studied here, can be regarded as semiperiodic systems. There, similar observations are made under the influence of a, typically external, electric field and are known as the quantum-confined Stark effect. Without any external perturbation, in oligopyrimidine SAMs one encounters an energy gradient that is generated by the dipole moments of the pyrimidine repeat units. It is particularly strong, reaching values of about 1.6 eV/nm, which corresponds to a substantial electric field of 1.6 × 10(7) V/cm. Close-lying σ- and π-states turn out to be a particular complication for a reliable description of the present systems, as their order is influenced not only by the docking groups and bonding to the metal, but also by the chosen computational approach. In the latter context we demonstrate that deliberately picking a hybrid functional allows avoiding pitfalls due to the infamous self-interaction error. Our results show that when aiming to build a monolayer with a specific electronic structure one can not only resort to the traditional technique of modifying the molecular structure of the constituents, but also try to exploit collective electronic effects.

Long-range corrected hybrid functionals for π-conjugated systems: dependence of the range-separation parameter on conjugation length.

Long-range corrected (range-separated hybrid) functionals represent a relatively new class of functionals for generalized Kohn-Sham theory that have proven to be very successful, for instance, when it comes to predicting ionization potentials and energy gaps for a wide range of molecules and solids. The results obtained from long-range corrected density functional theory approaches can be improved dramatically, if the range-separation parameter (ω) is optimized for each system separately. In this work, we have optimized ω for a series of π-conjugated molecular systems of increasing length by forcing the resulting functionals to obey the ionization potential-theorem, i.e., that their highest occupied eigenvalue be equal to the ΔSCF ionization potential. The optimized ω values are observed to vary substantially from their default values for the functionals. For highly conjugated chains such as oligoacenes and polyenes, we find that the characteristic length scale of the range-separation, i.e., 1/ω, grows almost linearly with the number of repeat units, for saturated alkane chains, however, 1/ω quickly saturates after 5-6 repeat units. For oligothiophenes, we find that 1/ω grows linearly for the shorter oligomers but then saturates at around 10 repeat units. Our results point to a close relation between the optimal range-separation parameter and the degree of conjugation in the system.

Metal Sulfide Nanoparticle Synthesis with Ionic Liquids - State of the Art and Future Perspectives.

Metal sulfides are among the most promising materials for a wide variety of technologically relevant applications ranging from energy to environment and beyond. Incidentally, ionic liquids (ILs) have been among the top research subjects for the same applications and also for inorganic materials synthesis. As a result, the exploitation of the peculiar properties of ILs for metal sulfide synthesis could provide attractive new avenues for the generation of new, highly specific metal sulfides for numerous applications. This article therefore describes current developments in metal sulfide nanoparticle synthesis as exemplified by a number of highlight examples. Moreover, the article demonstrates how ILs have been used in metal sulfide synthesis and discusses the benefits of using ILs over more traditional approaches. Finally, the article demonstrates some technological challenges and how ILs could be used to further advance the production and specific property engineering of metal sulfide nanomaterials, again based on a number of selected examples.

Size-Dependence of Nonempirically Tuned DFT Starting Points for G0W0 Applied to π-Conjugated Molecular Chains.

G0W0 calculations for predicting vertical ionization potentials (IPs) and electron affinities of molecules and clusters are known to show a significant dependence on the density functional theory (DFT) starting point. A number of nonempirical procedures to find an optimal starting point have been proposed, typically based on tuning the amount of HF exchange in the underlying hybrid functional specifically for the system at hand. For the case of π-conjugated molecular chains, these approaches lead to a significantly different amount of HF exchange for different oligomer sizes. In this study, we analyze if and how strongly this size dependence affects the ability of nonempirical tuning approaches to predict accurate IPs for π-conjugated molecular chains of increasing chain length. To this end, we employ three different nonempirical tuning procedures for the G0W0 starting point to calculate the IP of polyene oligomers up to 22 repeat units and compare the results to highly accurate coupled-cluster calculations. We find that, despite its size dependence, using an IP-tuned hybrid functional as a starting point for G0W0 yields excellent agreement with the reference data for all chain lengths.

Accurate Valence Ionization Energies from Kohn-Sham Eigenvalues with the Help of Potential Adjustors.

An accurate yet computationally very efficient and formally well justified approach to calculate molecular ionization potentials is presented and tested. The first as well as higher ionization potentials are obtained as the negatives of the Kohn-Sham eigenvalues of the neutral molecule after adjusting the eigenvalues by a recently [ Görling Phys. Rev. B 2015 , 91 , 245120 ] introduced potential adjustor for exchange-correlation potentials. Technically the method is very simple. Besides a Kohn-Sham calculation of the neutral molecule, only a second Kohn-Sham calculation of the cation is required. The eigenvalue spectrum of the neutral molecule is shifted such that the negative of the eigenvalue of the highest occupied molecular orbital equals the energy difference of the total electronic energies of the cation minus the neutral molecule. For the first ionization potential this simply amounts to a ΔSCF calculation. Then, the higher ionization potentials are obtained as the negatives of the correspondingly shifted Kohn-Sham eigenvalues. Importantly, this shift of the Kohn-Sham eigenvalue spectrum is not just ad hoc. In fact, it is formally necessary for the physically correct energetic adjustment of the eigenvalue spectrum as it results from ensemble density-functional theory. An analogous approach for electron affinities is equally well obtained and justified. To illustrate the practical benefits of the approach, we calculate the valence ionization energies of test sets of small- and medium-sized molecules and photoelectron spectra of medium-sized electron acceptor molecules using a typical semilocal (PBE) and two typical global hybrid functionals (B3LYP and PBE0). The potential adjusted B3LYP and PBE0 eigenvalues yield valence ionization potentials that are in very good agreement with experimental values, reaching an accuracy that is as good as the best G0W0 methods, however, at much lower computational costs. The potential adjusted PBE eigenvalues result in somewhat less accurate ionization energies, which, however, are almost as accurate as those obtained from the most commonly used G0W0 variants.

How Bond Length Alternation and Thermal Disorder Affect the Optical Excitation Energies of π-Conjugated Chains: A Combined Density Functional Theory and Molecular Dynamics Study.

We dissect the sources of error leading to inaccuracies in the description of the geometry and optical excitation energies of π-conjugated polymers. While the ground-state bond length alternation is shown to be badly reproduced by standard functionals, the recently adapted functionals PBEh* and ωPBE* as well as the double hybrid functional XYGJ-OS manage to replicate results obtained at the CCSD(T) level. By analysis of the bond length alternation in the excited state, a sensitive dependence of the exciton localization on the long-range behavior of the functional and the amount of Hartree-Fock exchange present is shown. Introducing thermal disorder through molecular dynamics simulations allows the consideration of a range of thermally accessible configurations of each oligomer, including trans to cis rotations, which break the conjugation of the backbone. Thermal disorder has a considerable effect when combined with functionals that overestimate the delocalization of the excitation, such as B3LYP. For functionals with a larger amount of exact exchange such as our PBEh* and ωPBE*, however, the effect is small, as excitations are often localized enough to fit between twists in the chain.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules II: Non-Empirically Tuned Long-Range Corrected Hybrid Functionals.

The performance of non-empirically tuned long-range corrected hybrid functionals for the prediction of vertical ionization potentials (IPs) and electron affinities (EAs) is assessed for a set of 24 organic acceptor molecules. Basis set-extrapolated coupled cluster singles, doubles, and perturbative triples [CCSD(T)] calculations serve as a reference for this study. Compared to standard exchange-correlation functionals, tuned long-range corrected hybrid functionals produce highly reliable results for vertical IPs and EAs, yielding mean absolute errors on par with computationally more demanding GW calculations. In particular, it is demonstrated that long-range corrected hybrid functionals serve as ideal starting points for non-self-consistent GW calculations.

Ionization Energies, Electron Affinities, and Polarization Energies of Organic Molecular Crystals: Quantitative Estimations from a Polarizable Continuum Model (PCM)-Tuned Range-Separated Density Functional Approach.

We propose a new methodology for the first-principles description of the electronic properties relevant for charge transport in organic molecular crystals. This methodology, which is based on the combination of a nonempirical, optimally tuned range-separated hybrid functional with the polarizable continuum model, is applied to a series of eight representative molecular semiconductor crystals. We show that it provides ionization energies, electron affinities, and transport gaps in very good agreement with experimental values, as well as with the results of many-body perturbation theory within the GW approximation at a fraction of the computational costs. Hence, this approach represents an easily applicable and computationally efficient tool to estimate the gas-to-crystal phase shifts of the frontier-orbital quasiparticle energies in organic electronic materials.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods.

The performance of different GW methods is assessed for a set of 24 organic acceptors. Errors are evaluated with respect to coupled cluster singles, doubles, and perturbative triples [CCSD(T)] reference data for the vertical ionization potentials (IPs) and electron affinities (EAs), extrapolated to the complete basis set limit. Additional comparisons are made to experimental data, where available. We consider fully self-consistent GW (scGW), partial self-consistency in the Green's function (scGW0), non-self-consistent G0W0 based on several mean-field starting points, and a "beyond GW" second-order screened exchange (SOSEX) correction to G0W0. We also describe the implementation of the self-consistent Coulomb hole with screened exchange method (COHSEX), which serves as one of the mean-field starting points. The best performers overall are G0W0+SOSEX and G0W0 based on an IP-tuned long-range corrected hybrid functional with the former being more accurate for EAs and the latter for IPs. Both provide a balanced treatment of localized vs delocalized states and valence spectra in good agreement with photoemission spectroscopy (PES) experiments.

Long-Range Corrected DFT Meets GW: Vibrationally Resolved Photoelectron Spectra from First Principles.

We propose an entirely nonempirical and computationally efficient scheme to calculate highly reliable vibrationally resolved photoelectron spectra for molecules from first principles. To this end, we combine nonempirically tuned long-range corrected hybrid functionals with non-self-consistent many-body perturbation theory in the G0W0 approximation and a Franck-Condon multimode analysis based on DFT-calculated frequencies. The vibrational analysis allows for a direct comparison of the GW-calculated spectra to gas-phase ultraviolet photoelectron measurements of neutral and anionic molecules, respectively. Direct comparison of the calculated peak maxima with experiment yields mean absolute errors below 0.1 eV for ionization potentials, electron affinities, and fundamental gaps, clearly outperforming commonly used G0W0 approaches at similar numerical costs.

Calculating Optical Absorption Spectra of Thin Polycrystalline Organic Films: Structural Disorder and Site-Dependent van der Waals Interaction.

We propose a new approach for calculating the change of the absorption spectrum of a molecule when moved from the gas phase to a crystalline morphology. The so-called gas-to-crystal shift Δ[Formula: see text] m is mainly caused by dispersion effects and depends sensitively on the molecule's specific position in the nanoscopic setting. Using an extended dipole approximation, we are able to divide Δ[Formula: see text] m = -QWm in two factors, where Q depends only on the molecular species and accounts for all nonresonant electronic transitions contributing to the dispersion while Wm is a geometry factor expressing the site dependence of the shift in a given molecular structure. The ability of our approach to predict absorption spectra is demonstrated using the example of polycrystalline films of 3,4,9,10-perylenetetracarboxylic diimide (PTCDI).

Size Effects in the Interface Level Alignment of Dye-Sensitized TiO2 Clusters.

The efficiency of dye-sensitized solar cells (DSCs) depends critically on the electronic structure of the interfaces in the active region. We employ recently developed dispersion-inclusive density functional theory (DFT) and GW methods to study the electronic structure of TiO2 clusters sensitized with catechol molecules. We show that the energy level alignment at the dye-TiO2 interface is the result of an intricate interplay of quantum size effects and dynamic screening effects and that it may be manipulated by nanostructuring and functionalizing the TiO2. We demonstrate that the energy difference between the catechol LUMO and the TiO2 LUMO, which is associated with the injection loss in DSCs, may be reduced significantly by reducing the dimensions of nanostructured TiO2 and by functionalizing the TiO2 with wide-gap moieties, which contribute additional screening but do not interact strongly with the frontier orbitals of the TiO2 and the dye. Precise control of the electronic structure may be achieved via "interface engineering" in functional nanostructures.

Sequence dependence of electron-induced DNA strand breakage revealed by DNA nanoarrays.

The electronic structure of DNA is determined by its nucleotide sequence, which is for instance exploited in molecular electronics. Here we demonstrate that also the DNA strand breakage induced by low-energy electrons (18 eV) depends on the nucleotide sequence. To determine the absolute cross sections for electron induced single strand breaks in specific 13 mer oligonucleotides we used atomic force microscopy analysis of DNA origami based DNA nanoarrays. We investigated the DNA sequences 5'-TT(XYX)3TT with X = A, G, C and Y = T, BrU 5-bromouracil and found absolute strand break cross sections between 2.66 · 10(-14) cm(2) and 7.06 · 10(-14) cm(2). The highest cross section was found for 5'-TT(ATA)3TT and 5'-TT(ABrUA)3TT, respectively. BrU is a radiosensitizer, which was discussed to be used in cancer radiation therapy. The replacement of T by BrU into the investigated DNA sequences leads to a slight increase of the absolute strand break cross sections resulting in sequence-dependent enhancement factors between 1.14 and 1.66. Nevertheless, the variation of strand break cross sections due to the specific nucleotide sequence is considerably higher. Thus, the present results suggest the development of targeted radiosensitizers for cancer radiation therapy.