Search

In this work, we use density functional theory to investigate the electronic structure of poly(3,4-ethylenedioxythiophene) (PEDOT) oligomers with co-located AlCl4− anions, a promising combination for energy storage. The 1980s bipolaron model remains the dominant interpretation of the electronic structure of PEDOT despite recent theoretical progress that has provided new definitions of bipolarons and polarons. By considering the influence of oligomer length, oxidation or anion concentration and spin state, we find no evidence for many of the assertions of the 1980s bipolaron model and so further contribute to a new understanding. No self-localisation of positive charges in PEDOT is found, as predicted by the bipolaron model at the hybrid functional level. Instead, our results show distortions that exhibit a single or a double peak in bond length alternations and charge density. Either can occur at different oxidation or anion concentrations. Rather than representing bipolarons or polaron pairs in the original model, these are electron distributions driven by a range of factors. Distortions can span an arbitrary number of nearby anions. We also contribute a novel conductivity hypothesis. Conductivity in conducting polymers has been observed to reduce at anion concentrations above 0.5. We show that at high anion concentrations, the energy of the localised, non-bonding anionic orbitals approaches that of the system HOMO due to Coulombic repulsion between anions. We hypothesize that with nucleic motion in the macropolymer, these orbitals will interfere with the hopping of charge carriers between sites of similar energy, lowering conductivity. [ABSTRACT FROM AUTHOR]

Density functional theory (DFT) is often used for simulating extended materials such as infinite crystals or surfaces, under periodic boundary conditions (PBCs). In such calculations, when the simulation cell has non-zero charge, electrical neutrality has to be imposed, and this is often done via a uniform background charge of opposite sign ("jellium"). This artificial neutralization does not occur in reality, where a different mechanism is followed as in the example of a charged electrode in electrolyte solution, where the surrounding electrolyte screens the local charge at the interface. The neutralizing effect of the surrounding electrolyte can be incorporated within a hybrid quantum–continuum model based on a modified Poisson–Boltzmann equation, where the concentrations of electrolyte ions are modified to achieve electroneutrality. Among the infinite possible ways of modifying the electrolyte charge, we propose here a physically optimal solution, which minimizes the deviation of concentrations of electrolyte ions from those in open boundary conditions (OBCs). This principle of correspondence of PBCs with OBCs leads to the correct concentration profiles of electrolyte ions, and electroneutrality within the simulation cell and in the bulk electrolyte is maintained simultaneously, as observed in experiments. This approach, which we call the Neutralization by Electrolyte Concentration Shift (NECS), is implemented in our electrolyte model in the Order-N Electronic Total Energy Package (ONETEP) linear-scaling DFT code, which makes use of a bespoke highly parallel Poisson–Boltzmann solver, DL_MG. We further propose another neutralization scheme ("accessible jellium"), which is a simplification of NECS. We demonstrate and compare the different neutralization schemes on several examples. [ABSTRACT FROM AUTHOR]

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Finally, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms. [ABSTRACT FROM AUTHOR]

First principles electronic structure calculations are typically performed in terms of molecular orbitals (or bands), providing a straightforward theoretical avenue for approximations of increasing sophistication, but do not usually provide any qualitative chemical information about the system. We can derive such information via post-processing using natural bond orbital (NBO) analysis, which produces a chemical picture of bonding in terms of localized Lewis-type bond and lone pair orbitals that we can use to understand molecular structure and interactions. We present NBO analysis of large-scale calculations with the ONETEP linear-scaling density functional theory package, which we have interfaced with the NBO 5 analysis program. In ONETEP calculations involving thousands of atoms, one is typically interested in particular regions of a nanosystem whilst accounting for long-range electronic effects from the entire system. We show that by transforming the Non-orthogonal Generalized Wannier Functions of ONETEP to natural atomic orbitals, NBO analysis can be performed within a localized region in such a way that ensures the results are identical to an analysis on the full system. We demonstrate the capabilities of this approach by performing illustrative studies of large proteins-namely, investigating changes in charge transfer between the heme group of myoglobin and its ligands with increasing system size and between a protein and its explicit solvent, estimating the contribution of electronic delocalization to the stabilization of hydrogen bonds in the binding pocket of a drug-receptor complex, and observing, in situ, the n → π* hyperconjugative interactions between carbonyl groups that stabilize protein backbones. © 2012 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

A new method for the constrained optimization of non-orthogonal, spatially localized orbitals using direct energy minimization techniques, in the context of electronic structure calculations, is presented. The variationally localized search direction (VLSD) method, as it was named, ensures that strict localization constraints are imposed upon the search direction vectors exactly, analytically and in a fully variational fashion. In contrast, the truncated search direction (TSD) method, of standard use in many electronic structure approaches with localization constraints, relies on the approximation that the truncated search direction vectors of the unconstrained problem resemble the exact search direction vectors of the constrained problem. With the TSD method, in order to maintain the localization constraints, a part of the pre-calculated information that is stored in the search direction vectors has to be deleted via an ad hoc, non-variational truncation step. The results on an extensive set of test molecules show that, in general, calculations with the VLSD method require less iterations to converge than with the TSD method for any size of the localization region. It was found that in calculations on certain systems where the TSD method is forced to delete a very large amount of information, the VLSD method is capable of achieving convergence in up to three times less iterations. Validation tests show that structural and electronic properties calculated with either method are accurate and in agreement with other electronic structure approaches. [ABSTRACT FROM AUTHOR]