No Arabic abstract
We implemented the derivative of the free energy functional with respect to the atom displacements, so called force, within the combination of Density Functional Theory and the Embedded Dynamical Mean Field Theory. We show that in combination with the numerically exact quantum Monte Carlo (MC) impurity solver, the MC noise cancels to a great extend, so that the method can be used very efficiently for structural optimization of correlated electron materials. As an application of the method, we show how strengthening of the fluctuating moment in FeSe superconductor leads to a substantial increase of the anion height, and consequently to a very large effective mass, and also strong orbital differentiation.
Materials with correlated electrons often respond very strongly to external or internal influences, leading to instabilities and states of matter with broken symmetry. This behavior can be studied theoretically either by evaluating the linear response characteristics, or by simulating the ordered phases of the materials under investigation. We developed the necessary tools within the dynamical mean-field theory (DMFT) to search for electronic instabilities in materials close to spin-state crossovers and to analyze the properties of the corresponding ordered states. This investigation, motivated by the physics of LaCoO$_3$, led to a discovery of condensation of spinful excitons in the two-orbital Hubbard model with a surprisingly rich phase diagram. The results are reviewed in the first part of the article. Electronic correlations can also be the driving force behind structural transformations of materials. To be able to investigate correlation-induced phase instabilities we developed and implemented a formalism for the computation of total energies and forces within a fully charge self-consistent combination of density functional theory and DMFT. Applications of this scheme to the study of structural instabilities of selected correlated electron materials such as Fe and FeSe are reviewed in the second part of the paper.
We present a charge and self-energy self-consistent computational scheme for correlated systems based on the Korringa-Kohn-Rostoker (KKR) multiple scattering theory with the many-body effects described by the means of dynamical mean field theory (DMFT). The corresponding local multi-orbital and energy dependent self-energy is included into the set of radial differential equations for the single-site wave functions. The KKR Greens function is written in terms of the multiple scattering path operator, the later one being evaluated using the single-site solution for the $t$-matrix that in turn is determined by the wave functions. An appealing feature of this approach is that it allows to consider local quantum and disorder fluctuations on the same footing. Within the Coherent Potential Approximation (CPA) the correlated atoms are placed into a combined effective medium determined by the dynamical mean field theory (DMFT) self-consistency condition. Results of corresponding calculations for pure Fe, Ni and Fe$_{x}$Ni$_{1-x}$ alloys are presented.
The stationary functional of the all-electron density functional plus dynamical mean field theory (DFT+DMFT) formalism to perform free energy calculations and structural relaxations is implemented for the first time. Here, the first order error in the density leads to a much smaller, second order error in the free energy. The method is applied to several well known correlated materials; metallic SrVO$_3$, Mott insulating FeO, and elemental Cerium, to show that it predicts the lattice constants with very high accuracy. In Cerium, we show that our method predicts the iso-structural transition between the $alpha$ and $gamma$ phases, and resolve the long standing controversy in the driving mechanism of this transition.
We investigate the effect of charge self-consistency (CSC) in density functional theory plus dynamical mean-field theory (DFT+DMFT) calculations compared to simpler one-shot calculations for materials where interaction effects lead to a strong redistribution of electronic charges between different orbitals or between different sites. We focus on two systems close to a metal-insulator transition, for which the importance of CSC is currently not well understood. Specifically, we analyze the strain-related orbital polarization in the correlated metal CaVO$_3$ and the spontaneous electronic charge disproportionation in the rare-earth nickelate LuNiO$_3$. In both cases, we find that the CSC treatment reduces the charge redistribution compared to cheaper one-shot calculations. However, while the MIT in CaVO$_3$ is only slightly shifted due to the reduced orbital polarization, the effect of the site polarization on the MIT in LuNiO$_3$ is more subtle. Furthermore, we highlight the role of the double-counting correction in CSC calculations containing different inequivalent sites.
The most general way to describe localized atomic-like electronic states in strongly correlated compounds is to utilize Wannier functions. In the present paper we continue the development of widely-spread DFT+U method onto Wannier function basis set and propose the technique to calculate the Hubbard contribution to the forces. The technique was implemented as a part of plane-waves pseudopotential code Quantum-ESPRESSO and successfully tested on a charge transfer insulator NiO.