The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. However, the gap in the band structure of the exact multiplicative Kohn-Sham (KS) potential substantially underestimates the fundamental gap, a major limitation of KS densityfunctional theory. Here, we give a simple proof of a theorem: In generalized KS theory (GKS), the band gap of an extended system equals the fundamental gap for the approximate functional if the GKS potential operator is continuous and the density change is delocalized when an electron or hole is added. Our theorem explains how GKS band gaps from metageneralized gradient approximations (meta-GGAs) and hybrid functionals can be more realistic than those from GGAs or even from the exact KS potential. The theorem also follows from earlier work. The band edges in the GKS one-electron spectrum are also related to measurable energies. A linear chain of hydrogen molecules, solid aluminum arsenide, and solid argon provide numerical illustrations. T he most basic property of a periodic solid is its fundamental energy gap G, which vanishes for a metal but is positive for semiconductors and other insulators. G dominates many properties. As the unbound limit of an exciton series, G is an excitation energy of the neutral solid, but it is defined here as a difference of ground-state energies: If EðMÞ is the ground-state energy for a solid with a fixed number of nuclei and M electrons, and if M = N for electrical neutrality, thenis the difference between the first ionization energy IðNÞ and the first electron affinity AðNÞ of the neutral solid. Whereas I and A can be measured for a macroscopic solid, they can be computed directly (as ground-state energy differences) either by starting from finite clusters and extrapolating to infinite cluster size or (for I-A) by starting from a finite number of primitive unit cells, with periodic boundary condition on the surface of this finite collection, and extrapolating to an infinite number. Here we shall follow both approaches, which have been discussed in a recent study (1). (The energy to remove an electron to infinite separation cannot depend upon the crystal face through which it is removed, although the energy to remove an electron to a macroscopic separation, but much smaller than the dimensions of that face, may so depend. The gap is of course a bulk property.) Band-Gap Problem in Kohn-Sham Density-Functional TheoryKohn-Sham density-functional theory (2, 3) is a formally exact way to compute the ground-state energy and electron density of M interacting electrons in a multiplicative external potential. This theory sets up a fictitious system of noninteracting electrons with the same ground-state density as the real interacting system, found by solving self-consistent one-electron Schrödinger equations. These electrons move in a multiplicative effective Kohn-Sham (KS) potential, the sum of the external and Hartree potentials and the derivative of the density functional for...

A question of both fundamental as well as practical importance is the nature of one-dimensional carbon, in particular whether a one-dimensional carbon allotrope is polyynic or cumulenic, that is, whether bond-length alternation occurs or not. By combining the concept of aromaticity and antiaromaticity with the rule of Peierls distortion, the occurrence and magnitude of bond-length alternation in carbon chains with periodic boundary conditions and corresponding carbon rings as a function of the chain or ring length can be explained. The electronic properties of one-dimensional carbon depend crucially on the bond-length alternation. Whereas it is generally accepted that carbon chains in the limit of infinite length have a polyynic structure at the minimum of the potential energy surface with bond-length alternation, we show here that zero-point vibrations lead to an effective equalization of all carbon-carbon bond lengths and thus to a cumulenic structure.

We introduce new functionals for the Kohn–Sham correlation energy that are based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem and are named σ-functionals. Like in the well-established direct random phase approximation (dRPA), σ-functionals require as input exclusively eigenvalues σ of the frequency-dependent KS response function. In the new functionals, functions of σ replace the σ-dependent dRPA expression in the coupling-constant and frequency integrations contained in the ACFD theorem. We optimize σ-functionals with the help of reference sets for atomization, reaction, transition state, and non-covalent interaction energies. The optimized functionals are to be used in a post-self-consistent way using orbitals and eigenvalues from conventional Kohn–Sham calculations employing the exchange–correlation functional of Perdew, Burke, and Ernzerhof. The accuracy of the presented approach is much higher than that of dRPA methods and is comparable to that of high-level wave function methods. Reaction and transition state energies from σ-functionals exhibit accuracies close to 1 kcal/mol and thus approach chemical accuracy. For the 10 966 reactions of the W4-11RE reference set, the mean absolute deviation is 1.25 kcal/mol compared to 3.21 kcal/mol in the dRPA case. Non-covalent binding energies are accurate to a few tenths of a kcal/mol. The presented approach is highly efficient, and the post-self-consistent calculation of the total energy requires less computational time than a density-functional calculation with a hybrid functional and thus can be easily carried out routinely. σ-Functionals can be implemented in any existing dRPA code with negligible programming effort.

Recently, a new type of orbital-dependent functional for the Kohn–Sham (KS) correlation energy, σ-functionals, was introduced. Technically, σ-functionals are closely related to the well-known direct random phase approximation (dRPA). Within the dRPA, a function of the eigenvalues σ of the frequency-dependent KS response function is integrated over purely imaginary frequencies. In σ-functionals, this function is replaced by one that is optimized with respect to reference sets of atomization, reaction, transition state, and non-covalent interaction energies. The previously introduced σ-functional uses input orbitals and eigenvalues from KS calculations with the generalized gradient approximation (GGA) exchange–correlation functional of Perdew, Burke, and Ernzerhof (PBE). Here, σ-functionals using input orbitals and eigenvalues from the meta-GGA TPSS and the hybrid-functionals PBE0 and B3LYP are presented and tested. The number of reference sets taken into account in the optimization of the σ-functionals is larger than in the first PBE based σ-functional and includes sets with 3d-transition metal compounds. Therefore, also a reparameterized PBE based σ-functional is introduced. The σ-functionals based on PBE0 and B3LYP orbitals and eigenvalues reach chemical accuracy for main group chemistry. For the 10 966 reactions from the highly accurate W4-11RE reference set, the B3LYP based σ-functional exhibits a mean average deviation of 1.03 kcal/mol compared to 1.08 kcal/mol for the coupled cluster singles doubles perturbative triples method if the same valence quadruple zeta basis set is used. For 3d-transition metal chemistry, accuracies of about 2 kcal/mol are reached. The computational effort for the post-self-consistent evaluation of the σ-functional is lower than that of a preceding PBE0 or B3LYP calculation for typical systems.

Understanding and controlling exciton transport is a strategic way to enhance the optoelectronic properties of high-performance organic devices. In this article we study triplet exciton migration in crystalline poly(p-phenylene vinylene) polymer (PPV) using comprehensive electronic structure and quantum dynamical methods. We solve the coupled electron-nuclear dynamics for the triplet energy migrating between two neighboring Frenkel sites in J-and H-aggregate arrangements. From the two-site model we extract key parameters for use with a master-equation approach that allows us to treat nanosize systems where time-dependent Schrödinger equation becomes intractable. We calculate the transient exciton density evolution and determine the diffusion constants along the principal crystal axes of the PPV. The triplet diffusion is characterized by two distinctive components: fast intrachain, and slow interchain. At room temperature the interchain diffusion coefficients are found to be D a = 0.89 · 10 −2 cm 2 s −1 and D b = 1.49 · 10 −2 cm 2 s −1 along the respectiveā-andb-axes, and the intrachain is D c = 3.03 cm 2 s −1 along the fastc-axis. The exceptionally high exciton mobility along the π-conjugated backbone facilitates rapid triplet migration over long distances. Our results can be utilized in the design of efficient energy conversion and light-emitting devices with desired solid-state properties.J o u r n a l N a me , [ y e a r ] , [ v o l . ] , 1-10 | 1 arXiv:1907.07304v1 [cond-mat.mtrl-sci]