Fifth tutorial on QSGW+DMFT

The density loop

The standard way (valid for both QSGW and LDA) to close the external loop is by updating the density taking fully into account the dynamical contribution of DMFT, unlike what we have done in the previous tutorial.

Self-consistent update of the electronic density

Once you obtained a converged Sig.out.brd from the DMFT loop, you can can use lmfdmft to sum over all Matsubara’s frequencies and update the density until it is self-consistent with the impurity self-energy and the QSGW potential.

As a result of the last command, you have an updated density file rst.ni. You can monitor the difference with respect to the previous density by

grep 'RMS DQ=' log

To achieve self-consistency you shall repeat the last lmfdmft calculation until RSM DE< 2.0e-4, or any other tolerance you may wish although it’s not easy to get better than 1.0e-4. To do that you may prefer to insert the last command in a loop like