Abstract:Within the context of first principles techniques we present a theoretical and computational framework to quickly determine, at finite momentum, the self-consistent (longitudinal) charge response to an external perturbation, that enters the determination of the scattering cross section of inelastic scattering processes such as EELS. We also determine the (tranverse) charge response computed in short-circuit condition. The all-order quasimomentum expansion of the tranverse charge response to an atomic displacement are the Born effective charges, quadrupoles, octupoles etc. We demonstrate that the transverse charge response can be related to the longitudinal one via a well-defined long-range dielectric function. Our advancements lead to an efficient use of perturbation theory. Due to its more favorable scaling, our method provides an interesting computational alternative to the use of the 2n+1 theorem, especially for semiconductors and metals with large unit cells. For semiconductors, we compute the piezoelectric properties of a large cell solid-solution of semiconducting hafniun oxide containing 96 atoms. We here show that the clamped ion piezoelectric response can be decomposed into real-space localized contributions that mostly depend on the chemical environment, paving the way for the use of machine-learning techniques in the material search for optimized piezoelectrics. We further apply our methodology to determine the density response of metals. Here, the leading terms of the charge expansion are related to the Fermi energy shift of the potential and by Born effective charges which do not sum to zero over the atoms. We apply our developments to the TEM-EELS spectroscopy of lithium intercalated graphites, where we find that the use of the atomic form-factor in the long-wavelength limit does not take into account for the anisotropy of the atomic chemical bonding.

Linear-Scaling Density Matrix Perturbation Treatment Of Electric Fields In Solids

MIXTURE EFFECTIVE PERMITTIVITY SIMULATIONS USING IMLMQRF METHOD ON PRECONDITIONED EFIE

A Stabilized and Well-Conditioned Marching-on-in Degree Method for Transient Analysis of Dielectric Objects

Efficient parallel linear scaling method to get the response density matrix in all-electron real-space density-functional perturbation theory

Self-consistent Density Functional Calculations without Diagonalization

Linear scaling calculation of band edge states and doped semiconductors

Generalized Drude-Lorentz Model Complying with the Singularity Expansion Method

Dynamic Work Distribution for PM Algorithm

First principles calculations of dynamical Born effective charges, quadrupoles and higher order terms from the charge response in large semiconducting and metallic systems

Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations

Ensemble Density-Functional Perturbation Theory: Spatial Dispersion in Metals

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

Linear Scaling Approach for Optical Excitations Using Maximally Localized Wannier Functions

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

Diagonalization without Diagonalization: A Direct Optimization Approach for Solid-State Density Functional Theory

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

Improved Dielectric Response of Solids: Combining the Bethe-Salpeter Equation with the Random Phase Approximation

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System

Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes

Calculating photonic Green's functions using a non-orthogonal finite difference time domain method

Current-density implementation for calculating flexoelectric coefficients