Functionals.md

@@ -11,8 +11,13 @@ The functional has the generic form:
     &+ [1-\alpha_{\uparrow}(n_{\uparrow},n_{\downarrow})]\dfrac{\bm{j}_{\uparrow}^2}{2n_{\uparrow}}
     + [1-\alpha_{\downarrow}(n_{\uparrow},n_{\downarrow})]\dfrac{\bm{j}_{\downarrow}^2}{2n_{\downarrow}}
 \end{aligned}
-```
-TODO
+```  
+The functional is fitted to quantum Monte Carlo data, and for spin symmetric and uniform system it provides:
+```math
+E/E_{\textrm{ffg}}=\xi=0.40(1),\qquad\Delta/\varepsilon_F=0.504(24).
+```  
+For more info about fitting procedure and explicit form of functional terms see: [arXiv:1008.3933](https://arxiv.org/abs/1008.3933).
+
 ## Stabilization of ASLDA functional
 In the case of calculations for trapped system term $`\frac{\bm{j}_{\sigma}^2}{2n_{\sigma}}`$ is source of numerical instabilities. Precisely, for small density regions we have:  
 * $`n_{\sigma}\rightarrow 0`$, 
@@ -38,4 +43,22 @@ Role of $`f_{\textrm{reg.}}`$ is to exclude from computation regions of small de
 ```
 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.
\ No newline at end of file
+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.
+
+# SLDA - superfluid local density approximation
+The term in ASLDA functional that depends on the currents $`\bm{j}_{\sigma}`$ introduces a significant cost to the computation. This term is responsible for maintaining Gallilean invariance of the ASLDA theory. The effective mass was found to be consistent with the bare mass to within 10% for a large range of polarizations,  Therefore, in many applications, one can set $`\alpha_{\sigma}=1`$ without losing qualitative features of ASLDA theory. Under this ausumptin we obtain SLDA functional:
+```math
+\mathcal{E}_{\textrm{edf}}(n,\nu,\ldots)  =  \frac{\tau_{\uparrow}}{2} + 
+    \frac{\tau_{\downarrow}}{2} 
++
+    D(n_{\uparrow},n_{\downarrow})
+    +
+    g(n_{\uparrow},n_{\downarrow})\nu^{\dagger}\nu 
+```  
+
+*Note 1*: SLDA functional exhibits much better convergence properties than ASLDA, i.e `st-wslda` codes typically converge in a significantly smaller number of iterations.  
+
+*Note 2*: SLDA functional provides by factor about 2x better performance in case of `td-wslda` codes than ASLDA variant.  
+
+Below we provide plot showing relative energy change (in %) between simulations of SLDA and ASLDA that compute the energy of N fermions (both even and odd particle numbers) in the spin-symmetric unitary Fermi gas trapped in an isotropic harmonic oscillator.  For more info see [supplemental material of this paper](https://arxiv.org/abs/1306.4266).  
+![slda_vs_aslda](uploads/073efe244d1ab67bed1314af5cf0fc21/slda_vs_aslda.png)