Search

Halide double perovskites are a chemically-diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and {\it ab initio} many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here we apply the WOT-SRSH functional to five halide double perovskites, and compare the results with those obtained from other known functionals and previous $GW$ calculations. We also use the approach as a starting point for $GW$ calculations and we compute the band structures and optical absorption spectrum for Cs\textsubscript{2}Ag{Bi}Br\textsubscript{6}, using both time-dependent DFT and the $GW$-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

Van der Waals (vdW) polytypes of broken inversion and mirror symmetries were recently shown to exhibit switchable electric polarization even at the ultimate two-layer thin limit. Their out-of-plane polarization was found to accumulate in a ladder-like fashion with each successive layer, offering 2D building blocks for the bottom-up construction of 3D ferroelectrics. Here, we demonstrate experimentally that beyond a critical stack thickness, the accumulated polarization in rhombohedral polytypes of molybdenum disulfide (r-MoS2) saturates. The underlying saturation mechanism, deciphered via DFT and self-consistent Poisson-Schr\"odinger calculations, point to a purely electronic redistribution involving: (1) polarization-induced bandgap closure that allows for cross-stack charge transfer and the emergence of free surface charge; (2) reduction of the polarization saturation value, as well as the critical thickness at which it is obtained, by the presence of free carriers. The resilience of polar layered structures to atomic surface reconstruction, which is essentially unavoidable in polar 3D crystals, potentially allows for the design of new devices with mobile surface charges. Our findings, which are of general nature, should be accounted for when designing switching and/or conductive devices based on ferroelectric layered materials.

We present a first principles strategy for developing state-specific density functional approximations for excited states. We first clarify why approaches based on conventional ground state approximations miss density-driven correlations, by considering excited state physics through the lens of ensemble density functional theory. To solve this issue we gain insights on density driven correlations by exploiting the recently understood low-density limit of electrons in excited states. The theory developments are then combined to produce a proof-of-concept excited state approximation that resolves urgent paradigmatic failures (double excitations, charge transfer excitations, piecewise linearity) of existing state-of-art density-functional approaches, directly from differences of self-consistent field calculations; i.e., $\Delta$SCF. In light of its observed impressive performance, we conclude that the approach represents a major step toward unified and accurate modelling of neutral and charged excitations.

Accurate prediction of electronic and optical excitations in van der Waals (vdW) materials is a long-standing challenge for density functional theory. The recently proposed Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functional has proven successful in non-empirical determination of electronic band gaps and optical absorption spectra for various covalent and ionic crystals. However, for vdW materials the tuning of the material- and structure-dependent functional parameters has, until now, only been attained semi-empirically. Here, we present a non-empirical WOT-SRSH approach applicable to vdW materials, with the optimal functional parameters transferable between monolayer and bulk. We apply this methodology to prototypical vdW materials: black phosphorus, molybdenum disulfide, and hexagonal boron nitride (in the latter case including zero-point renormalization). We show that the WOT-SRSH approach consistently achieves accuracy levels comparable to experiments and ab initio many-body perturbation theory (MBPT) calculations for band structures and optical absorption spectra, both on its own and as an optimal starting point for MBPT calculations., Comment: 11 pages + 4 figures; Supporting Information (15 pages + 2 figures)

The ionization potential of molecular chains is well-known to be a tunable nano-scale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nano-sized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an \textit{ansatz} that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nano-scale regime.

Using both time-dependent density functional theory (TDDFT) and the ``single-shot" $GW$ plus Bethe-Salpeter equation ($GW$-BSE) approach, we compute optical band gaps and optical absorption spectra from first principles for eight common binary and ternary closed-shell metal oxides (MgO, Al$_2$O$_3$, CaO, TiO$_2$, Cu$_2$O, ZnO, BaSnO$_3$, and BiVO$_4$), based on the non-empirical Wannier-localized optimally-tuned screened range-separated hybrid functional. Overall, we find excellent agreement between our TDDFT and $GW$-BSE results and experiment, with a mean absolute error less than 0.4 eV, including for Cu$_2$O and ZnO, traditionally considered to be challenging for both methods.

An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2 where M = Mo, W and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with reference calculations of binding energy curves and sliding potential energy surfaces. It is found to be transferable to TMD junctions outside the training set that contain the same atom types. Calculated bulk moduli agree with most previous dispersion corrected DFT predictions, which underestimate available experimental values. Calculated phonon spectra of the various junctions under consideration demonstrate the importance of appropriately treating the anisotropic nature of layered interfaces. Considering our previous parameterization for MoS2, the interlayer potential enables accurate and efficient large-scale simulations of the dynamical, tribological, and thermal transport properties of a large set of homogeneous and heterogeneous TMD interfaces.

The halide perovskite CsPbBr$_3$ belongs to the Cs-Pb-Br material system, which features two additional thermodynamically stable ternary phases, Cs$_4$PbBr$_6$ and CsPb$_2$Br$_5$. The coexistence of these phases and their reportedly similar photoluminescence have resulted in a debate on the nature of the emission in these systems. Here, we combine optical and microscopic characterization with an effective mass, correlated electron-hole model of excitons in confined systems, to investigate the emission properties of the ternary phases in the Cs-Pb-Br system. We find that all Cs-Pb-Br phases exhibit green emission and the non-perovskite phases exhibit photoluminescence quantum yields orders of magnitude larger than CsPbBr$_3$. In particular, we measure blue- and red-shifted emission for the Cs- and Pb-rich phases, respectively, stemming from embedded CsPbBr$_3$ nanocrystals. Our model reveals that the difference in emission shift is caused by the combined effects of nanocrystal size and different band mismatch. Furthermore, we demonstrate the importance of including the dielectric mismatch in the calculation of the emission energy for Cs-Pb-Br composites. Our results explain the reportedly limited blue shift in CsPbBr$_3$@Cs$_4$PbBr$_6$ composites and rationalize some of its differences with CsPb$_2$Br$_5$.

The accurate description of electronic properties and optical absorption spectra is a long-standing challenge for density functional theory. Recently, the introduction of screened range-separated hybrid (SRSH) functionals for solid-state materials has allowed for the calculation of fundamental band gaps and optical absorption spectra that are in very good agreement with many-body perturbation theory. However, since solid-state SRSH functionals are typically tuned to reproduce the properties of bulk phases, their transferability to low-dimensional structures, which experience substantially different screening than in the bulk, remains an open question. In this work, we explore the transferability of SRSH functionals to several prototypical van der Waals materials, including transition-metal sulfides and selenides, indium selenide, black phosphorus, and hexagonal boron nitride. Considering the bulk and a monolayer of these materials as limiting cases, we show that the parameters of the SRSH functional can be determined systematically, using only the band-edge quasiparticle energies of these extremal structural phases as fitting targets. The resulting SRSH functionals can describe both electronic bandstructures and optical absorption spectra with accuracy comparable to more demanding ab initio many-body perturbation theory (GW and Bethe-Salpeter equation) approaches. Selected examples also demonstrate that the SRSH parameters, obtained from the bulk and monolayer reference structures, display good accuracy for bandstructures and optical spectra of bilayers, indicating a degree of transferability that is independent of the fitting procedure., Comment: 14 pages + 8 figures; Supporting Information (21 pages + 18 figures)

A crystalline solid is a periodic sequence of identical cells, each containing one or more atoms. If the constituting unit cell is not centrosymmetric, charge may distribute unevenly between the atoms, resulting in internal electric polarization. This effect serves as the basis for numerous ferroelectric, piezoelectric, and pyroelectric phenomena. In nearly all polar materials, including multilayered van der Waals stacks that were recently found to exhibit interfacial polarization, inversion symmetry is broken by having two or more atomic species within the unit cell. Here, we show that even elemental crystals, consisting of one type of atom, and composed of non-polar centrosymmetric layers, exhibit electric polarization if arranged in an appropriate three-dimensional architecture. This concept is demonstrated here for inversion and mirror asymmetric mixed-stacking tetra-layer polytypes of non-polar graphene sheets. Furthermore, we find that the room temperature out-of-plane electric polarization increases with external electrostatic doping, rather than decreases owing to screening. Using first-principles calculations, as well as tight-binding modeling, we unveil the origin of polytype-induced polarization and its dependence on doping. Extension of this idea to graphene multilayers suggests that solely by lateral shifts of constituent monolayers one can obtain multiple meta-stable interlayer stacking sequences that may allow for even larger electrical polarization.

We present first-principles calculations of the interface between GaN and strained AlN, using a slab model in which polarization is compensated via surface fractional-charge pseudo-hydrogen atoms. We show that an interface two-dimensional carrier electron or hole gas emerges naturally in response to different compensating surface charges, but that this need not involve in-gap surface states.

Consistency between the exchange-correlation (xc) functional used during pseudopotential construction and planewave-based electronic structure calculations is important for an accurate and reliable description of the structure and properties of condensed-phase systems. In this work, we present a general scheme for constructing pseudopotentials with range-separated hybrid (RSH) xc functionals based on the solution of the all-electron radial integro-differential equation for a spherical atomic configuration. As proof-of-principle, we demonstrate pseudopotential construction with the PBE, PBE0, HSE, and sRSH (based on LC-$\omega$PBE0) xc functionals for a select set of atoms, and then investigate the importance of pseudopotential consistency when computing band gaps, equilibrium lattice parameters, bulk moduli, and atomization energies of several solid-state systems. In doing so, we find that pseudopotential consistency errors (PSCE) tend to be systematic and can be as large as $1.4\%$ when computing these properties.

Two-dimensional (2D) halide perovskites offer a unique platform to investigate the ground-state of materials possessing significant anharmonicity. In contrast to three-dimensional perovskites, their 2D counterparts offer substantially fewer degrees of freedom, resulting in multiple well-defined crystal structures. In this work, we thoroughly investigate the anharmonic ground-state of the benchmark (PEA)2PbI4 compound, using complementary information from low-temperature x-ray diffraction (XRD) and photoluminescence spectroscopy, supported by density functional theory (DFT) calculations. We extrapolate four crystallographic configurations from the low-temperature XRD. These configurations imply that the ground-state has an intrinsic disorder stemming from two coexisting chiral sub-lattices, each with a bi-oriented organic spacer molecule. We further show evidence that these chiral structures form unevenly populated ground-states, portraying uneven anharmonicity, where the state population may be tuned by surface effects. Our results uncover a disordered ground-state that may induce intrinsic grain boundaries, which cannot be ignored in practical applications., Comment: manuscripts: 16 pages, 3 figures, 1 table. SI: 14 pages, 5 figures, 9 tables

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals., Comment: 14 pages, 8 figures

Ferroelectricity in atomically thin bilayer structures has been recently predicted1 and measured[2-4] in two-dimensional (2D) materials with hexagonal non-centrosymmetric unit-cells. Interestingly, the crystal symmetry translates lateral shifts between parallel 2D layers to a change of sign in their out-of-plane electric polarization, a mechanism referred to as "Slide-Tronics"[4]. These observations, however, have been restricted to switching between only two polarization states under low charge carrier densities[5-12], strongly limiting the practical application of the revealed phenomena[13]. To overcome these issues, one needs to explore the nature of the polarization that arises in multi-layered van der Waals (vdW) stacks, how it is governed by intra- and inter-layer charge redistribution, and to which extent it survives the introduction of mobile charge carriers, all of which are presently unknown14. To explore these questions, we conduct surface potential measurements of parallel WSe2 and MoS2 multi-layers with aligned and anti-aligned configurations of the polar interfaces. We find evenly spaced, nearly decoupled potential steps, indicating highly confined interfacial electric fields, which provide means to design multi-state "ladder ferroelectrics". Furthermore, we find that the internal polarization remains significant upon electrostatic doping of a mobile charge carrier density as high as 1013 cm-2, with substantial in-plane conductivity. Using first-principles calculations based on density functional theory (DFT), we trace the extra charge redistribution in real and momentum space and identify an eventual doping-induced depolarization mechanism.

The dependence of ab initio many-body perturbation theory within the $GW$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $GW$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely PBE, PBE0, and HSE, as well as to self-consistent $GW$ schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $GW$ (QS$GW$) and eigenvalue self-consistent $GW$ (ev$GW$). We also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

We provide a critical overview of the theory of the chirality-induced spin selectivity (CISS) effect, i.e., phenomena in which the chirality of molecular species imparts significant spin selectivity to various electron processes. Based on discussions in a recently held workshop, and further work published since, we review the status of CISS effects - in electron transmission, electron transport, and chemical reactions. For each, we provide a detailed discussion of the state-of-the-art in theoretical understanding and identify remaining challenges and research opportunities.

Accurate prediction of fundamental band gaps of crystalline solid state systems entirely within density functional theory is a long standing challenge. Here, we present a simple and inexpensive method that achieves this by means of non-empirical optimal tuning of the parameters of a screened range-separated hybrid functional. The tuning involves the enforcement of an ansatz that generalizes the ionization potential theorem to the removal of an electron in an occupied state described by a localized Wannier function in a modestly sized supercell calculation. The method is benchmarked against experiment for a set of systems ranging from narrow band gap semiconductors to large band gap insulators, spanning a range of fundamental band gaps from 0.2 to 14.2 eV and is found to yield quantitative accuracy across the board, with a mean absolute error of $\sim$0.1 eV and a maximal error of $\sim$0.2 eV., Comment: 10 pages, 2 figures

The importance of many-body dispersion effects in layered materials subjected to high external loads is evaluated. State-of-the-art many-body dispersion density functional theory calculations performed for graphite, hexagonal boron nitride, and their hetero-structures were used to fit the parameters of a classical registry-dependent interlayer potential. Using the latter, we performed extensive equilibrium molecular dynamics simulations and studied the mechanical response of homogeneous and heterogeneous bulk models under hydrostatic pressures up to 30 GPa. Comparison with experimental data demonstrates that the reliability of the many-body dispersion model extends deep into the sub-equilibrium regime. Friction simulations demonstrate the importance of many-body dispersion effects for the accurate description of the tribological properties of layered materials interfaces under high pressure., Comment: 40 pages, 18 figures

We consider the Br vacancy in CsPbBr3 as a prototype for the impact of structural dynamics on defect energetics in halide perovskites (HaPs). Using first-principles molecular dynamics based on density functional theory, we find that the static picture of defect energetics breaks down; the energy of the Br vacancy level is found to be intrinsically dynamic, oscillating by as much as 1 eV on the ps time scale at room temperature. These significant energy fluctuations are correlated with the distance between the neighboring Pb atoms across the vacancy and with the electrostatic potential at these Pb atomic sites. We expect this unusually strong coupling of structural dynamics and defect energetics to bear important implications for both experimental and theoretical analysis of defect characteristics in HaPs. It may also hold significant ramifications for carrier transport and defect tolerance in this class of photovoltaic materials., Comment: 5 figures, 1 table

Electronic coupling matrix elements are important to the theoretical description of electron transfer processes. However, they are notoriously difficult to obtain accurately from time- dependent density functional theory (TDDFT). Here, we use the HAB11 benchmark dataset of coupling matrix elements to assess whether TDDFT using optimally-tuned range-separated hybrid functionals, already known to be successful for the description of charge transfer excitation energies, also allows for an improved accuracy in the prediction of coupling matrix elements. We find that this approach outperforms all previous TDDFT calculations, based on semi-local, hybrid, or non-tuned range-separated hybrid functionals, with a remaining average deviation as low as ~12%. We discuss potential sources for the remaining error., Comment: Mol. Phys., in press (Michael Baer festschrift)

By studying the lowest excitations of an exactly solvable one-dimensional molecular model, we show that components of Kohn-Sham ensembles can be used to describe charge transfers. Furthermore, we compute the approximate excitation energies obtained by using thee exact ensemble densities in the recently formulated ensemble Hartree-exchange theory [Gould and Pittalis, Phys. Rev. Lett. 119, 243001 (2017)]. Remarkably, our results show that triplet excitations are accurately reproduced across a dissociation curve in all cases tested, even in systems where ground state energies are poor due to strong static correlations. Singlet excitations exhibit larger deviations from exact results but are still reproduced semi-quantitatively.

Hybrid organic-inorganic perovskites (HOIPs) have recently emerged as highly promising solution-processable materials for photovoltaic (PV) and other optoelectronic devices. HOIPs represent a broad family of materials with properties highly tuneable by the ions that make up the perovskite structure as well as their multiple combinations. Interestingly, recent high-efficiency PV devices using HOIPs with substantially improved long-term stability have used combinations of different ionic compositions. The structural dynamics of these systems are unique for semiconducting materials and are currently argued to be central to HOIPs stability and charge-transport properties. Here, we studied the impact of ionic composition on phonon speeds of HOIPs from Brillouin spectroscopy experiments and density functional theory calculations for FAPbBr$_3$, MAPbBr$_3$, MAPbCl$_3$, and the mixed halide MAPbBr$_{1.25}$Cl$_{1.75}$. Our results show that the acoustic phonon speeds can be strongly modified by ionic composition, which we explain by analysing the lead-halide sublattice in detail. The vibrational properties of HOIPs are therefore tuneable by using targeted ionic compositions in the perovskite structure. This tuning can be rationalized with non-trivial effects, for example, considering the influence of the shape and dipole moment of organic cations. This has an important implication to further improvements in the stability and charge-transport properties of these systems.

The alignment of the frontier orbital energies of an adsorbed molecule with the substrate Fermi level at metal-organic interfaces is a fundamental observable of significant practical importance in nanoscience and beyond. Typical density functional theory calculations, especially those using local and semi-local functionals, often underestimate level alignment leading to inaccurate electronic structure and charge transport properties. In this work, we develop a new fully self-consistent predictive scheme to accurately compute level alignment at certain classes of complex heterogeneous molecule-metal interfaces based on optimally-tuned range-separated hybrid functionals. Starting from a highly accurate description of the gas-phase electronic structure, our method by construction captures important nonlocal surface polarization effects via tuning of the long-range screened exchange in a range-separated hybrid in a non-empirical and system-specific manner. We implement this functional in a plane-wave code and apply it to several physisorbed and chemisorbed molecule-metal interface systems. Our results are in quantitative agreement with experiments, both the level alignment and work function changes. Our approach constitutes a new practical scheme for accurate and efficient calculations of the electronic structure of molecule-metal interfaces., Comment: 15 pages, 8 figures

Lead-halide perovskites are promising materials for opto-electronic applications. Recent reports indicated that their mechanical and electronic properties are strongly affected by the lattice vibrations. Herein we report far-infrared spectroscopy measurements of CH$_{3}$NH$_{3}$Pb(I/Br/Cl)$_{3}$ thin films and single crystals at room temperature and a detailed quantitative analysis of the spectra. We find strong broadening and anharmonicity of the lattice vibrations for all three halide perovskites, which indicates dynamic disorder of the lead-halide cage at room temperature. We determine the frequencies of the transversal and longitudinal optical phonons, and use them to calculate, via appropriate models, the static dielectric constants, polaron masses, electron-phonon coupling constants, and upper limits for the phonon-scattering limited charge carrier mobilities. Within the limitations of the model used, we can place an upper limit of 200$\,$cm$^{2}$V$^{-1}$s$^{-1}$ for the room temperature charge carrier mobility in MAPbI$_{3}$ single crystals. Our findings are important for the basic understanding of charge transport processes and mechanical properties in metal halide perovskites.

We observe large, reversible, bias driven changes in the vibrational energies of PCBM, based on simultaneous transport and surface-enhanced Raman spectroscopy (SERS) measurements on PCBM-gold junctions. A combination of linear and quadratic shifts in vibrational energies with voltage is analyzed and compared with similar measurements involving C60-gold junctions. A theoretical model based on density functional theory (DFT) calculations suggests that both a vibrational Stark effect and bias-induced charging of the junction contribute to the shifts in vibrational energies. In the PCBM case, a linear vibrational Stark effect is observed due to the permanent electric dipole moment of PCBM. The vibrational Stark shifts shown here for PCBM junctions are comparable to or larger than the charging effects that dominate in C60 junctions., Comment: 16 pages, 4 figures

Hybrid lead-halide perovskites have emerged as an excellent class of photovoltaic materials. Recent reports suggest that the organic molecular cation is responsible for local polar fluctuations that inhibit carrier recombination. We combine low frequency Raman scattering with first-principles molecular dynamics (MD) to study the fundamental nature of these local polar fluctuations. Our observations of a strong central peak in both hybrid (CH$_3$NH$_3$PbBr$_3$) and all-inorganic (CsPbBr$_3$) lead-halide perovskites show that anharmonic, local polar fluctuations are intrinsic to the general lead-halide perovskite structure, and not unique to the dipolar organic cation. MD simulations show that head-to-head Cs motion coupled to Br face expansion, on a few hundred femtosecond time scale, drives the local polar fluctuations in CsPbBr$_3$.

Molecular crystals are a prototypical class of van der Waals (vdW) bound organic materials with excited state properties relevant for optoelectronics applications. Predicting the structure and excited state properties of molecular crystals presents a challenge for electronic structure theory, as standard approximations to density functional theory (DFT) do not capture long range vdW dispersion interactions and do not yield excited state properties. In this work, we use a combination of DFT including vdW forces) using both non local correlation functionals and pair wise correction methods (together with many body perturbation theory (MBPT) to study the geometry and excited states, respectively, of the entire series of oligoacene crystals, from benzene to hexacene. We find that vdW methods can predict lattice constants within 1 percent of the experimental measurements, on par with the previously reported accuracy of pairwise approximations for the same systems. We further find that excitation energies are sensitive to geometry, but if optimized geometries are used MBPT can yield excited state properties within a few tenths of an eV from experiment. We elucidate trends in MBPT computed charged and neutral excitation energies across the acene series and discuss the role of common approximations used in MBPT.

We present an experimental and theoretical viewpoint on the electronic carrier mobilities of typical hybrid organic-inorganic perovskites (HOIPs). While these mobilities are often quoted as high, a review of them shows that although otherwise the semiconducting properties of HOIPs are impressively good, mobilities of HOIPs used in most solar cells are actually not that high. This is especially apparent if they are compared to those of inorganic semiconductors used in other high efficiency solar cells. We critically examine possible causes and focus on electron-lattice coupling mechanisms that are active at room temperature, and can lead to carrier scattering. From this, we propose scattering due to acoustic phonons or polarons as possible causes, but also point out the difficulties with each of these in view of additional experimental and theoretical findings in the literature. Further research in this direction will contribute to making HOIP solar cells even more efficient than they already are.

We investigate periodical oscillations in the conductance of suspended Au and Pt atomic chains during elongation under mechanical stress. Analysis of conductance and shot noise measurements reveals that the oscillations are mainly related to variations in a specific conduction channel as the chain undergoes transitions between zigzag and linear atomic configurations. The calculated local electronic structure shows that the oscillations originate from varying degrees of hybridization between the atomic orbitals along the chain as a function of the zigzag angle. These variations are highly dependent on the directionality and symmetry of the relevant orbitals, in agreement with the order-of-magnitude difference between the Pt and Au oscillation amplitudes observed in experiment. Our results demonstrate that the sensitivity of conductance to structural variations can be controlled by designing atomic-scale conductors in view of the directional interactions between atomic orbitals.

One-electron self-interaction and an incorrect asymptotic behavior of the Kohn-Sham exchange-correlation potential are among the most prominent limitations of many present-day density functionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct long-range potential. This is here shown explicitly for local hybrid functionals. Furthermore, carefully studying the ratio of the von Weizs\"acker kinetic energy density to the (positive) Kohn-Sham kinetic energy density, $\tau_\mathrm{W}/\tau$, reveals that this ratio, which frequently serves as an iso-orbital indicator and is used to eliminate one-electron self-interaction effects in meta-generalized-gradient approximations and local hybrid functionals, can fail to approach its expected value in the vicinity of orbital nodal planes. This perspective article suggests that the nature and consequences of one-electron self-interaction and some of the strategies for its correction need to be reconsidered., Comment: 12 pages, 8 figures

There are several approximations to the exchange-correlation functional in density-functional theory that accurately predict total energy-related properties of many-electron systems, such as binding energies, bond lengths, and crystal structures. Other approximations are designed to describe potential-related processes, such as charge transfer and photoemission. However, the development of a functional which can serve the two purposes simultaneously is a long-standing challenge. Trying to address it, we employ in the current work the ensemble generalization procedure proposed in Phys. Rev. Lett. 110, 126403 (2013). Focusing on the prediction of the ionization potential via the highest occupied Kohn-Sham eigenvalue, we examine a variety of exchange-correlation approximations: the local spin-density approximation, semi-local generalized gradient approximations, and global and local hybrid functionals. Results for a test set of 26 diatomic molecules and single atoms are presented. We find that the aforementioned ensemble generalization systematically improves the prediction of the ionization potential, for various systems and exchange-correlation functionals, without compromising the accuracy of total energy-related properties. We specifically examine hybrid functionals. These depend on a parameter controlling the ratio of semi-local to non-local functional components. The ionization potential obtained with ensemble-generalized functionals is found to depend only weakly on the parameter value, contrary to common experience with non-generalized hybrids, thus eliminating one aspect of the so-called `parameter dilemma' of hybrid functionals.

Many approximations within density-functional theory spuriously predict that a many-electron system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated diatomic molecules, known to exhibit this problem when studied within standard approaches, including the local spin-density approximation (LSDA). By employing our recently proposed [E. Kraisler and L. Kronik, Phys. Rev. Lett. 110, 126403 (2013)] ensemble-generalization we find that asymptotic fractional dissociation is eliminated in all systems examined, even if the underlying exchange-correlation (xc) is still the LSDA. Furthermore, as a result of the ensemble generalization procedure, the Kohn-Sham potential develops a spatial step between the dissociated atoms, reflecting the emergence of the derivative discontinuity in the xc energy functional. This step, predicted in the past for the exact Kohn-Sham potential and observed in some of its more advanced approximate forms, is a desired feature that prevents any fractional charge transfer between the system's fragments. It is usually believed that simple xc approximations such as the LSDA cannot develop this step. Our findings show, however, that ensemble generalization to fractional electron densities automatically introduces the desired step even to the most simple approximate xc functionals and correctly predicts asymptotic integer dissociation., Comment: includes a supplementary material

The fundamental gap is a central quantity in the electronic structure of matter. Unfortunately, the fundamental gap is not generally equal to the Kohn-Sham gap of density functional theory (DFT), even in principle. The two gaps differ precisely by the derivative discontinuity, namely, an abrupt change in slope of the exchange-correlation (xc) energy as a function of electron number, expected across an integer-electron point. Popular approximate functionals are thought to be devoid of a derivative discontinuity, strongly compromising their performance for prediction of spectroscopic properties. Here we show that, in fact, all exchange-correlation functionals possess a derivative discontinuity, which arises naturally from the application of ensemble considerations within DFT, without any empiricism. This derivative discontinuity can be expressed in closed form using only quantities obtained in the course of a standard DFT calculation of the neutral system. For small, finite systems, addition of this derivative discontinuity indeed results in a greatly improved prediction for the fundamental gap, even when based on the most simple approximate exchange-correlation density functional - the local density approximation (LDA). For solids, the same scheme is exact in principle, but when applied to LDA it results in a vanishing derivative discontinuity correction. This failure is shown to be directly related to the failure of LDA in predicting fundamental gaps from total energy differences in extended systems., Comment: 11 pages, 3 figures

The energy positions of frontier orbitals in organic electronic materials are often studied experimentally by (inverse) photoemission spectroscopy and theoretically within density functional theory. However, standard exchange-correlation functionals often result in too small fundamental gaps, may lead to wrong orbital energy ordering, and do not capture polarization-induced gap renormalization. Here, we examine these issues and a strategy for overcoming them by studying the gas phase and bulk electronic structure of the organic molecule quinacridone (5Q), a promising material with many interesting properties for organic devices. Experimentally, we perform angle-resolved photoemission spectroscopy (ARUPS) on thin films of the crystalline $\beta$-phase of 5Q. Theoretically, we employ an optimally-tuned range-separated hybrid functional (OT-RSH) within density functional theory. For the gas phase molecule, our OT-RSH result for the ionization potential (IP) represents a substantial improvement over the semi-local PBE and the PBE0 hybrid functional results, producing an IP in quantitative agreement with experiment. For the bulk crystal, we take into account the correct screening in the bulk, using the recently developed OT-SRSH approach, while retaining the optimally-tuned parameters for the range-separation and the short-range Fock exchange. This leads to a band gap narrowing due to polarization effects and results in a valence band spectrum in excellent agreement with experimental ARUPS data, with respect to both peak positions and heights. Finally, full-frequency $G_0W_0$ results based on a hybrid functional starting point are shown to agree with the OT-SRSH approach, improving substantially on the PBE-starting point.