| ... | ... | @@ -2,6 +2,7 @@ |
|
|
|
title: KKRSusc Quickstart Guide
|
|
|
|
---
|
|
|
|
|
|
|
|
|
|
|
|
# Perform a KKRsusc calculation
|
|
|
|
|
|
|
|
In order to investigate spin-excitations via the KKRsusc program please
|
| ... | ... | @@ -58,7 +59,7 @@ To calculate the self-energy and renormalised spectrum of the system, a fitting |
|
|
|
- Use the file [emesh.dat](emesh.dat), that describes an energy contour from $-0.3$ to $1.3Ry$ to run one iteration using the host code with the run flag KKRSUSC like in step 2.
|
|
|
|
- Similarly to the susceptibility calculation in step 4, two runs of the code are needed.
|
|
|
|
- Set `lfit= T` to instruct the code to run the fitting procedure and set the numerator polynomial degree `numd` and denominator degree `dend` to proper values, with $dend=numd +1$, to match the $1/\omega$ behaviour of the susceptibility in high frequencies. A typical range for `numd` is between 10-30, with the best fits usually in the region of 17-21.
|
|
|
|
- Set the bias voltage range to match the frequency range defined by `omegamin, omegamax` and the corresponding parameters described in [newinpsusc.dat](newinpsusc.dat).
|
|
|
|
- Set the bias voltage range (`vbiasmin`, `vbiasmax`) to match the frequency range defined by `omegamin, omegamax` and the corresponding parameters described in [newinpsusc.dat](newinpsusc.dat).
|
|
|
|
- Run the mpi version of the code with `lhdio= T` and `lrestart= F` and `ldos= T` so that the DOS is calculated as well.
|
|
|
|
- Next perform the serial run with `lhdio= F` and `lrestart= T`
|
|
|
|
- After the calculation is done, plot the dos and compare it with the one you get from the impurity code. The goal here is to get a fit as loyal as possible.
|
| ... | ... | |
