VERSION>=2021.09.18

Tracking the states

W-SLDA Toolkit allows to track observables arsing from selected states in time. Presently following quantities can be tracked:

• subset_rho_a:
  n_{\uparrow}^{\textrm{(subset)}}(r)= \sum_{E_{\textrm{min}}^{\textrm{(subset)}}<E_n<E_{\textrm{max}}^{\textrm{(subset)}}}|u_{n,\uparrow}(r)|^2 f_{\beta}(E_n)
• subset_rho_b:
  n_{\downarrow}^{\textrm{(subset)}}(r) = \sum_{E_{\textrm{min}}^{\textrm{(subset)}}<E_n<E_{\textrm{max}}^{\textrm{(subset)}}}|v_{n,\downarrow}(r)|^2 f_{\beta}(-E_n)
• subset_j_a_x, subset_j_a_y, subset_j_a_z:
  \vec{j}_{\uparrow}^{\textrm{(subset)}}(r) = -\sum_{E_{\textrm{min}}^{\textrm{(subset)}}<E_n<E_{\textrm{max}}^{\textrm{(subset)}}} \textrm{Im}[u_{n,\uparrow}(r)\nabla u_{n,\uparrow}^*(r)]f_{\beta}(E_n)
• subset_j_b_x, subset_j_b_y, subset_j_b_z:
  \vec{j}_{\downarrow}^{\textrm{(subset)}}(r) = \sum_{E_{\textrm{min}}^{\textrm{(subset)}}<E_n<E_{\textrm{max}}^{\textrm{(subset)}}} \textrm{Im}[v_{n,\downarrow}(r)\nabla v_{n,\downarrow}^*(r)]f_{\beta}(-E_n)

To control interval range E_{\textrm{min}}^{\textrm{(subset)}}<E_n<E_{\textrm{max}}^{\textrm{(subset)}} use input file variables:

subsetMinEn             0.0     # in eF units, deafault=0.0
subsetMaxEn             0.0     # in eF units, deafault=0.0

In the case of spin-imbalanced systems, it is also convenient to introduce a shift of quasiparticle energies when selecting the states: E_n \rightarrow E_n+\frac{\mu_{\uparrow}-\mu_{\downarrow}}{2}. This can be done automatically via enabling

subsetShiftDmu          1       # if 1 then apply extra shift of quasiparticle energies by (mu_a-mu_b)/2, default=0

Example

TODO: Andrea