اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدAll-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
71
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
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.
قيم البحث
اقرأ أيضاً
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 re
sults 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 o
ccupied 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 comput
ationally 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.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
