Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAll-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
71
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Within density-functional theory, perturbation theory~(PT) is the state-of-the-art formalism for assessing the response to homogeneous electric fields and the associated material properties, e.g., polarizabilities, dielectric constants, and Raman intensities. Here we derive a real-space formulation of PT and present an implementation within the all-electron, numeric atom-centered orbitals electronic structure code FHI-aims that allows for massively-parallel calculations. As demonstrated by extensive validation, this allows the rapid computation of accurate response properties of molecules and solids. As an application showcase, we present harmonic and anharmonic Raman spectra, the latter obtained by combining hundreds of thousands of PT calculations with textit{ab initio} molecular dynamics. By using the PBE exchange-correlation functional with many-body van der Waals corrections, we obtain spectra in good agreement with experiment especially with respect to lineshapes for the isolated paracetamol molecule and two polymorphs of the paracetamol crystal.
rate research
Read More
The methods of density-functional perturbation theory may be used to calculate various physical response properties of insulating crystals including elastic, dielectric, Born charge, and piezoelectric tensors. These and other important tensors may be defined as second derivatives of the total energy with respect to atomic-displacement, electric-field, or strain perturbations, or as mixed derivatives with respect to two of these perturbations. The resulting tensor quantities tend to be coupled in complex ways in polar crystals, giving rise to a variety of variant definitions. For example, it is generally necessary to distinguish between elastic tensors defined under different electrostatic boundary conditions, and between dielectric tensors defined under different elastic boundary conditions. Here, we describe an approach for computing all of these various response tensors in a unified and systematic fashion. Applications are presented for two materials, wurtzite ZnO and rhombohedral BaTiO3, at zero temperature.
We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections.
A real-space formalism for density-functional perturbation theory (DFPT) is derived and applied for the computation of harmonic vibrational properties in molecules and solids. The practical implementation using numeric atom-centered orbitals as basis functions is demonstrated exemplarily for the all-electron Fritz Haber Institute ab initio molecular simulations (FHI-aims) package. The convergence of the calculations with respect to numerical parameters is carefully investigated and a systematic comparison with finite-difference approaches is performed both for finite (molecules) and extended (periodic) systems. Finally, the scaling tests and scalability tests on massively parallel computer systems demonstrate the computational efficiency.
The bandstructure of gold is calculated using many-body perturbation theory (MBPT). Different approximations within the GW approach are considered. Standard single shot G0W0 corrections shift the unoccupied bands up by ~0.2 eV and the first sp-like occupied band down by ~0.4 eV, while leaving unchanged the 5d occupied bands. Beyond G0W0, quasiparticle self-consistency on the wavefunctions lowers the occupied 5d bands by 0.35 eV. Globally, many-body effects achieve an opening of the interband gap (5d-6sp gap) of 0.35 to 0.75 eV approaching the experimental results. Finally, the quasiparticle bandstructure is compared to the one obtained by the widely used HSE (Heyd, Scuseria, and Ernzerhof) hybrid functional.
We present an efficient implementation of the PBE0 hybrid functional within the full-potential linearized augmented-plane-wave (FLAPW) method. The Hartree-Fock exchange term, which is a central ingredient of hybrid functionals, gives rise to a computationally expensive nonlocal potential in the one-particle Schroedinger equation. The matrix elements of this exchange potential are calculated with the help of an auxiliary basis that is constructed from products of FLAPW basis functions. By representing the Coulomb interaction in this basis the nonlocal exchange term becomes a Brillouin-zone (BZ) sum over vector-matrix-vector products. We show that the Coulomb matrix can be made sparse by a suitable unitary transformation of the auxiliary basis, which accelerates the computation of the vector-matrix-vector products considerably. Additionally, we exploit spatial and time-reversal symmetry to identify the nonvanishing exchange matrix elements in advance and to restrict the k summations for the nonlocal potential to an irreducible set of k points. Favorable convergence of the self-consistent-field cycle is achieved by a nested density-only and density-matrix iteration scheme. We discuss the convergence with respect to the parameters of our numerical scheme and show results for a variety of semiconductors and insulators, including oxide materials, where the PBE0 hybrid functional improves the band gaps and the description of localized states in comparison with the PBE functional. Furthermore, we find that in contrast to conventional local exchange-correlation functionals ferromagnetic EuO is correctly predicted to be a semiconductor.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
