DFT fails to describe the ground state and/or the excited states such as many lanthanides, actinides or transition metals. Indeed, exchange correlation functionals are not (yet) able to describe the strong repulsive Coulomb interactions occurring among electrons in partly filled localized d or f orbitals.
A way to improve the description of strongly correlated systems is to explicitly include the strong repulsive Coulomb interactions in the Hamiltonian. Solving it in the static mean field approximation, gives the DFT+U method ([Anisimov1991],[Liechtenstein1995]), implemented in ABINIT [Amadon2008a]. The Dynamical Mean Field Theory [Georges1996] (DMFT), goes beyond, by solving exactly the local correlations for an atom in an effective field (i.e., an Anderson model). The effective field reproduces the effect of the surrounding correlated atoms and is thus self-consistently related to the solution of the Anderson model [Georges1996].
The combination of DFT with DMFT [Georges2004],[Kotliar2006] (usedmft= 1) relies on :
The practical solution of the DFT+DMFT scheme is usually presented as a double loop over, first, the local Green's function, and second the electronic local density [Amadon2012]. The number of iterations of the two loops are determined by dmft_iter and nstep. However, in the general case, the most efficient way to carry out fully consistent DFT+DMFT calculations is to keep only the loop governed by nstep, while dmft_iter=1 [Bieder2014], dmft_rslf= 1 (to read the self-energy file at each step of the DFT loop) and prtden= -1 (to be able to restart the calculation of each step of the DFT loop from the density file). Lastly, one linear and one logarithmic grid are used for Matsubara frequencies [Kotliar2006] determined by dmft_nwli and dmft_nwlo (Typical values are 10000 and 100, but convergence should be studied). More information can be obtained in the log file by setting pawprtvol=3.
The main output of the calculations are the imaginary time Green's function , from which spectral functions can be obtained using an external maximum entropy code [Bergeron2016], self-energies, from which quasiparticle renormalization weight can be extracted, the density matrix of correlated orbitals, and the internal energies [Amadon2006]. The electronic entropic contribution to the free energy can also be obtained using dmft_entropy and dmft_nlambda.
The efficient CTQMC code in ABINIT, which is the most time consuming part of DMFT, uses the hybridization expansion [Werner2006],[Gull2011] with a density-density multiorbital interaction [Gull2011]. Moreover, the hybridization function [Gull2011] is assumed to be diagonal in the orbital (or flavor) index. This is exact for cubic symmetry without spin orbit coupling but, in general, one should always check that the off-diagonal terms are much smaller than the diagonal ones. A link to the exact rotationally invariant interaction CTQMC code of the TRIQS library is also available using dmft_solv=7.
As the CTQMC solver uses a Fourier transform, the time grid dmftqmc_l in imaginary space must be chosen so that the Nyquist frequency, defined by π*dmftqmc_l*tsmear, is around 2 or 3 Ha.
A convergence study should be performed on this variable.
Moreover, the number of imaginary frequencies (dmft_nwlo) has to be set to at least twice the value of dmftqmc_l.
Typical numbers of steps for the thermalization (dmftqmc_therm) and for the Monte carlo runs (dmftqmc_n) are 106 and 109 respectively.
The random number generator can be initialized with the variable dmftqmc_seed.
Several other variables are available.
dmftctqmc_order gives a histogram of the perturbation orders during the simulation, dmftctqmc_gmove customizes the global move tries (mainly useful for systems with high/low spin configurations), and dmftctqmc_meas sets the frequency of measurement of quantities.
Go to the top
Compulsory input variables:
... dmft_iter [Dynamical Mean Fied Theory: number of ITERation]
... dmft_nwli [Dynamical Mean Fied Theory: Number of frequency omega (W) in the LInear mesh]
... dmft_nwlo [Dynamical Mean Fied Theory: Number of frequency omega (W) in the LOg mesh]
... dmftbandf [Dynamical Mean Field Theory: BAND: Final]
... dmftbandi [Dynamical Mean Field Theory: BAND: Initial]
... dmftqmc_l [Dynamical Mean Fied Theory: Quantum Monte Carlo time sLices]
... dmftqmc_n [Dynamical Mean Fied Theory: Quantum Monte Carlo Number of sweeps]
... dmftqmc_therm [Dynamical Mean Fied Theory: Quantum Monte Carlo THERMalization]
... usedmft [USE Dynamical Mean Field Theory]
Basic input variables:
... dmft_rslf [Dynamical Mean Fied Theory: Read SeLF energy]
... dmft_solv [Dynamical Mean Fied Theory: choice of SOLVer]
Useful input variables:
... dmft_dc [Dynamical Mean Fied Theory: Double Counting]
... dmft_mxsf [Dynamical Mean Fied Theory: MiXing parameter for the SelF energy]
... dmft_tollc [Dynamical Mean Fied Theory: TOLerance on Local Charge for convergence of the DMFT loop]
... dmftcheck [Dynamical Mean Fied Theory: CHECKs]
... dmftctqmc_check [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo CHECK]
... dmftctqmc_gmove [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo Global MOVEs]
... dmftctqmc_order [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo perturbation ORDER]
... dmftqmc_seed [Dynamical Mean Fied Theory: Quantum Monte Carlo SEED]
Input variables for experts:
... dmft_entropy [Dynamical Mean Fied Theory: ENTROPY]
... dmft_nlambda [Dynamical Mean Fied Theory: Number of LAMBDA points]
... dmft_read_occnd [Dynamical Mean Fied Theory: READ OCCupations (Non Diagonal)]
... dmft_t2g [Dynamical Mean Fied Theory: t2g orbitals]
... dmft_tolfreq [Dynamical Mean Fied Theory: TOLerance on DFT correlated electron occupation matrix for the definition of the FREQuency grid]
... dmftctqmc_basis [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo BASIS]
... dmftctqmc_correl [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo CORRELations]
... dmftctqmc_grnns [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo GReeNs NoiSe]
... dmftctqmc_meas [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo MEASurements]
... dmftctqmc_mov [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo MOVie]
... dmftctqmc_mrka [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo MARKov Analysis]
... dmftctqmc_triqs_nleg [Dynamical Mean Fied Theory: Continuous Time Quantum Monte Carlo perturbation of TRIQS, Number of LEGendre polynomials]
Go to the top
tests/v6/Input: t47.in
tests/v7/Input: t21.in
t27.in
t28.in
Go to the top