|
|
[[_TOC_]]
|
|
|
# Introduction
|
|
|
W-SLDA Toolkit utilizes a local pairing field $`\Delta(\bm{r})`$. In such case, renormalization
|
|
|
W-SLDA Toolkit utilizes a local pairing field $`\Delta(\bm{r})`$. In such a case, renormalization
|
|
|
procedure is required. There are two predefined regularization schemes that can be selected in `predefines.h` file:
|
|
|
```c
|
|
|
/**
|
|
|
* Scheme of pairing field renormalization procedure.
|
|
|
* SPHERICAL_CUTOFF: use spherical momentum space cutoff, in this case you need to set `ec` variable in input file (default).
|
|
|
* SPHERICAL_CUTOFF: use spherical momentum space cutoff, in this case, you need to set `ec` variable in the input file (default).
|
|
|
* CUBIC_CUTOFF: use cubic momentum space cutoff, in this case `ec` will be set to infinity automatically.
|
|
|
* */
|
|
|
#define REGULARIZATION_SCHEME SPHERICAL_CUTOFF
|
... | ... | @@ -50,8 +50,10 @@ The effective coupling constant is computed according to prescription: |
|
|
```math
|
|
|
\dfrac{1}{g_{\textrm{eff}}}=\dfrac{1}{g_0} - \dfrac{m}{2\alpha_+}\dfrac{K}{2\hbar^2\pi^2 dx},
|
|
|
```
|
|
|
where $`K=2.442 75`$ is a numerical constant. In this formula we assume that all states contribute to the densities. Physically it means that we take into account states up to maximal value of energy set by lattice, which is of the order $`E_c\approx 3\frac{\hbar^2\pi^2}{2mdx^2}`$ (assuming that $`dx=dy=dz`$).
|
|
|
where $`K=2.442 75`$ is a numerical constant. In this formula, we assume that all states contribute to the densities. Physically it means that we take into account states up to the maximal value of energy set by lattice, which is of the order $`E_c\approx 3\frac{\hbar^2\pi^2}{2mdx^2}`$ (assuming that $`dx=dy=dz`$).
|
|
|
*Note*: when working with this renormalization scheme value of tag `ec` will be ignored.
|
|
|
|
|
|
# Custom renormalization scheme
|
|
|
Static codes allow for defining your own renormalization scheme. You need to provide the formula in `void modify_potentials(...)` function. See [here](Strict 2D or 1D mode) for example.
|
|
|
|
|
|
# Regularization scheme and the energy conservation in td calculations |
|
|
\ No newline at end of file |