Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userLattice Dynamics Calculations based on Density-functional Perturbation Theory in Real Space
89
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We present a constrained density functional perturbation theory scheme for the calculation of structural and harmonic vibrational properties of insulators in the presence of an excited and thermalized electron-hole plasma. The method is ideal to tame ultrafast light induced structural transitions in the regime where the photocarriers thermalize faster than the lattice, the electron-hole recombination time is longer than the phonon period and the photocarrier concentration is large enough to be approximated by an electron-hole plasma. The complete derivation presented here includes total energy, forces and stress tensor, variable cell structural optimization, harmonic vibrational properties and the electron-phonon interaction. We discuss in detail the case of zone center optical phonons not conserving the number of electrons and inducing a Fermi shift in the photo-electron and hole distributions. We validate our implementation by comparing with finite differences in Te and VSe2. By calculating the evolution of the phonon spectrum of Te, Si and GaAs as a function of the fluence of the incoming laser light, we demonstrate that even at low fluences, corresponding to approximately 0.1 photocarriers per cell, the phonon spectrum is substantially modified with respect to the ground state one with new Kohn anomalies appearing and a substantial softening of zone center optical phonons. Our implementation can be efficiently used to detect reversible transient phases and irreversible structural transition induced by ultrafast light absorption.
We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $alpha$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $alpha$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant).
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.
Linear-response time-dependent density-functional theory (TDDFT) can describe excitonic features in the optical spectra of insulators and semiconductors, using exchange-correlation (xc) kernels behaving as $-1/k^{2}$ to leading order. We show how excitons can be modeled in real-time TDDFT, using an xc vector potential constructed from approximate, long-range corrected xc kernels. We demonstrate for various materials that this real-time approach is consistent with frequency-dependent linear response, gives access to femtosecond exciton dynamics following short-pulse excitations, and can be extended with some caution into the nonlinear regime.
The knowledge of effective masses is a key ingredient to analyze numerous properties of semiconductors, like carrier mobilities, (magneto-)transport properties, or band extrema characteristics yielding carrier densities and density of states. Currently, these masses are usually calculated using finite-difference estimation of density functional theory (DFT) electronic band curvatures. However, finite differences require an additional convergence study and are prone to numerical noise. Moreover, the concept of effective mass breaks down at degenerate band extrema. We assess the former limitation by developing a method that allows to obtain the Hessian of DFT bands directly, using density functional perturbation theory (DFPT). Then, we solve the latter issue by adapting the concept of `transport equivalent effective mass to the $vec{k} cdot hat{vec{p}}$ framework. The numerical noise inherent to finite-difference methods is thus eliminated, along with the associated convergence study. The resulting method is therefore more general, more robust and simpler to use, which makes it especially appropriate for high-throughput computing. After validating the developed techniques, we apply them to the study of silicon, graphane, and arsenic. The formalism is implemented into the ABINIT software and supports the norm-conserving pseudopotential approach, the projector augmented-wave method, and the inclusion of spin-orbit coupling. The derived expressions also apply to the ultrasoft pseudopotential method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
