Band structure interpolation using optimized local orbitals from linear-scaling density functional theory

© 2018 American Physical Society. Several approaches to linear-scaling density functional theory (LS-DFT) that seek to achieve accuracy equivalent to plane-wave methods do so by optimizing in situ a set of local orbitals in terms of which the density matrix can be accurately expressed. These local orbitals, which can also accurately represent the canonical Kohn-Sham orbitals, qualitatively resemble the maximally localized Wannier functions employed in band structure interpolation. As LS-DFT methods are increasingly being used in real-world applications demanding accurate band structures, it is natural to question the extent to which these optimized local orbitals can provide sufficient accuracy. In this paper, we present and compare, in principle and in practice, two methods for obtaining band structures. We apply these to a (10, 0) carbon nanotube as an example. By comparing with the results from a traditional plane-wave pseudopotential calculation, the optimized local orbitals are found to provide an excellent description of the occupied bands and some low-lying unoccupied bands, with consistent agreement across the Brillouin zone. However free-electron-like states derived from weakly bound states independent of the σ and π orbitals can only be found if additional local orbitals are included

Unravelling the roles of size, ligands and pressure in the piezochromic properties of CdS nanocrystals

Understanding the effects of pressure-induced deformations on the optoelectronic properties of nanomaterials is important not only from the fundamental point of view, but also for po- tential applications such as stress sensors and electromechanical devices. Here we describe the novel insights into these piezochromic ef- fects gained from using a linear-scaling den- sity functional theory framework and an elec- tronic enthalpy scheme, which allow us to ac- curately characterize the electronic structure of CdS nanocrystals with a zincblende-like core of experimentally relevant size. In particu- lar we focus on unravelling the complex inter- play of size and surface (phenyl) ligands with pressure. We show that pressure-induced de- formations are not simple isotropic scaling of the original structures and that the change in HOMO-LUMO gap with pressure results from two competing factors: (i) a bulk-like linear in- crease due to compression, which is offset by (ii) distortions/disorder and, to a lesser ex- tent, orbital hybridization induced by ligands affecting the frontier orbitals. Moreover, we observe that the main peak in the optical ab- sorption spectra is systematically red-shifted or blue-shifted, as pressure is increased up to 5 GPa, depending on the presence or absence of phenyl ligands. These heavily hybridize the frontier orbitals, causing a reduction in over- lap and oscillator strength, so that at zero pres- sure the lowest energy transition involves deeper hole orbitals than in the case of hydrogen- capped nanocrystals; the application of pressure induces greater delocalisation over the whole nanocrystals bringing the frontier hole orbitals into play and resulting in an unexpected red shift for the phenyl-capped nanocrystals, in part caused by distortions. In response to a growing interest in relatively small nanocrystals that can be difficult to accurately characterize with ex- perimental techniques, this work exemplifies the detailed understanding of structure-property re- lationships under pressure that can be obtained for realistic nanocrystals with state-of-the-art first principles methods and used for the charac- terization and design of devices based on these and similar nanomaterials

Implicit and explicit host effects on excitons in pentacene derivatives

An ab initio study of the effects of implicit and explicit hosts on the excited state properties of pentacene and its nitrogen-based derivatives has been performed using ground state density func- tional theory (DFT), time-dependent DFT and ∆SCF. We observe a significant solvatochromic redshift in the excitation energy of the lowest singlet state (S 1 ) of pentacene from inclusion in a p -terphenyl host compared to vacuum; for an explicit host consisting of six nearest neighbour p -terphenyls, we obtain a redshift of 65 meV while a conductor-like polarisable continuum model (CPCM) yields a 78 meV redshift. Comparison is made between the excitonic properties of pen- tacene and four of its nitrogen-based analogues, 1,8-, 2,9-, 5,12-, and 6,13-diazapentacene with the latter found to be the most distinct due to local distortions in the ground state electronic struc- ture. We observe that a CPCM is insufficient to fully understand the impact of the host due to the presence of a mild charge-transfer (CT) coupling between the chromophore and neighbouring p -terphenyls, a phenomenon which can only be captured using an explicit model. The strength of this CT interaction increases as the nitrogens are brought closer to the central acene ring of pentacene

Linear-scaling time-dependent density-functional theory in the linear response formalism

We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on a system of carbon nanotubes

On the Growth of Al_2 O_3 Scales

Understanding the growth of Al2O3 scales requires knowledge of the details of the chemical reactions at the scale–gas and scale–metal interfaces, which in turn requires specifying how the creation/annihilation of O and Al vacancies occurs at these interfaces. The availability of the necessary electrons and holes to allow for such creation/annihilation is a crucial aspect of the scaling reaction. The electronic band structure of polycrystalline Al2O3 thus plays a decisive role in scale formation and is considered in detail, including the implications of a density functional theory (DFT) calculation of the band structure of a Σ7 View the MathML source bicrystal boundary, for which the atomic structure of the boundary was known from an independent DFT energy-minimization calculation and comparisons with an atomic-resolution transmission electron micrograph of the same boundary. DFT calculations of the formation energy of O and Al vacancies in bulk Al2O3 in various charge states as a function of the Fermi energy suggested that electronic conduction in Al2O3 scales most likely involves excitation of both electrons and holes, which are localized on singly charged O vacancies, View the MathML source and doubly charged Al vacancies, View the MathML source, respectively. We also consider the variation of the Fermi level across the scale and bending (“tilting”) of the conduction band minimum and valence band maximum due to the electric field developed during the scaling reaction. The band structure calculations suggest a new mechanism for the “reactive element” effect—a consequence of segregation of Y, Hf, etc., to grain boundaries in Al2O3 scales, which results in improved oxidation resistance—namely, that the effect is due to the modification of the near-band edge grain-boundary defect states rather than any blocking of diffusion pathways, as previously postulated. Secondly, Al2O3 scale formation is dominated by grain boundary as opposed to lattice diffusion, and there is unambiguous evidence for both O and Al countercurrent transport in Al2O3 scale-forming alloys. We postulate that such transport is mediated by migration of grain boundary disconnections containing charged jogs, rather than by jumping of isolated point defects in random high-angle grain boundaries

Simulation of electron energy loss spectra of nanomaterials with linear-scaling density functional theory.

Experimental techniques for electron energy loss spectroscopy (EELS) combine high energy resolution with high spatial resolution. They are therefore powerful tools for investigating the local electronic structure of complex systems such as nanostructures, interfaces and even individual defects. Interpretation of experimental electron energy loss spectra is often challenging and can require theoretical modelling of candidate structures, which themselves may be large and complex, beyond the capabilities of traditional cubic-scaling density functional theory. In this work, we present functionality to compute electron energy loss spectra within the onetep linear-scaling density functional theory code. We first demonstrate that simulated spectra agree with those computed using conventional plane wave pseudopotential methods to a high degree of precision. The ability of onetep to tackle large problems is then exploited to investigate convergence of spectra with respect to supercell size. Finally, we apply the novel functionality to a study of the electron energy loss spectra of defects on the (1 0 1) surface of an anatase slab and determine concentrations of defects which might be experimentally detectable.This work was performed using the Darwin Supercomputer of the University of Cambridge High Performance Computing Service (http://www.hpc.cam.ac.uk/), provided by Dell Inc. using Strategic Research Infrastructure Funding from the Higher Education Funding Council for England and funding from the Science and Technology Facilities Council. E.W.T. was supported by the EPSRC Cambridge NanoDTC, EP/G037221/1. L.E.R. was supported, in part, by DOE Office of Science User Facility under Contract DE-AC02-06CH11357. N.D.M.H. Acknowledges the support of the Winton Programme for the Physics of Sustainability. We wish to acknowledge the use of the EPSRC's Chemical Database Service at Daresbury

Tuning the role of charge-transfer states in intramolecular singlet exciton fission through side-group engineering

Understanding the mechanism of singlet exciton fission, in which a singlet exciton separates into a pair of triplet excitons, is crucial to the development of new chromophores for efficient fission-sensitized solar cells. The challenge of controlling molecular packing and energy levels in the solid state precludes clear determination of the singlet fission pathway. Here, we circumvent this difficulty by utilizing covalent dimers of pentacene with two types of side groups. We report rapid and efficient intramolecular singlet fission in both molecules, in one case via a virtual charge-transfer state and in the other via a distinct charge-transfer intermediate. The singlet fission pathway is governed by the energy gap between singlet and charge-transfer states, which change dynamically with molecular geometry but are primarily set by the side group. These results clearly establish the role of charge-transfer states in singlet fission and highlight the importance of solubilizing groups to optimize excited-state photophysics.S.L. thanks AGS(O) Scholarship support from A*STAR Singapore. J.W. acknowledges financial support from MOE Tier 3 grant (MOE2014-T3-1-004). This work was supported by the Engineering and Physical Sciences Research Council, U.K. (Grant numbers EP/M005143/1 and EP/G060738/1). D.H.P.T. and N.D.M.H. acknowledge the Winton Programme for the Physics of Sustainability. K.C. and J.M.H. acknowledge support from a Rutherford Discovery Fellowship to J.M.H

Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer's, Parkinson's and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth

Research Data Supporting "Impact of carbonyl formation on cobalt ripening over titania surface"

Input and output for the simulations performed to obtain the results of the publication "Impact of carbonyl formation on cobalt ripening over titania surface". All calculations were carried out in VASP 5.3.3 with the Vasp Transition State Tools (VTST-Tools) maintained by the Henkelman group at UT Austin.Shell Global Solutions International B.V