... | ... | @@ -58,17 +58,17 @@ Role of $`f_{\textrm{reg.}}`$ is to exclude from computation regions of small de |
|
|
/**
|
|
|
* Meaningful only in case of ASLDA.
|
|
|
* Parameters defining stabilization procedure of ASLDA functional.
|
|
|
* For regions with density smaller than ASLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY
|
|
|
* For regions with density smaller than SLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY
|
|
|
* contribution from current term j^2/2n is assumed to be zero.
|
|
|
* For regions with density above ASLDA_STABILIZATION_RETAIN_ABOVE_DENSITY
|
|
|
* For regions with density above SLDA_STABILIZATION_RETAIN_ABOVE_DENSITY
|
|
|
* the contribution is assumed to be intact by stabilization procedure.
|
|
|
* */
|
|
|
#define ASLDA_STABILIZATION_RETAIN_ABOVE_DENSITY 1.0e-5
|
|
|
#define ASLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY 1.0e-7
|
|
|
#define SLDA_STABILIZATION_RETAIN_ABOVE_DENSITY 1.0e-5
|
|
|
#define SLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY 1.0e-7
|
|
|
```
|
|
|
There meaning is presented on the figure below.
|
|
|
![stabilization](uploads/740d294dbfb30a203cabec0b140fa326/stabilization.png)
|
|
|
Smooth transition between densities `ASLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY` and `ASLDA_STABILIZATION_RETAIN_ABOVE_DENSITY` is introduced in order to avoid discontinuities for quantities, that may lead to divergences when computing derivatives.
|
|
|
Smooth transition between densities `SLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY` and `SLDA_STABILIZATION_RETAIN_ABOVE_DENSITY` is introduced in order to avoid discontinuities for quantities, that may lead to divergences when computing derivatives.
|
|
|
|
|
|
## Chemical potential
|
|
|
The spin-symmetric and uniform unitary Fermi is expected to be scale-invariant where the following relations for the total energy and the chemical potential satisfy: $`E=\xi E_\textrm{ffg}`$ and $`\mu=\xi \varepsilon_{\textrm{F}}`$, where $`\xi\approx 0.4`$ is the Bertch parameter. However, in numerical realization, the relation for the chemical potential is satisfied only in the limit of very low densities. It is demonstrated in the table below, where results for the uniform solution obtained on lattice `128x128x128` with `DX=1` is presented, and $`k_{\textrm{F}}=\sqrt{2\varepsilon_{\textrm{F}}}=(6\pi^2 n_{\uparrow})^{1/3}`$.
|
... | ... | |