NOTE: these are archival examples, and provided files may not be compatible with the present API of the W-SLDA Toolkit!
- Example 1: gas confined in a tube
- Example 2: vortex solution within BdG
- Example 3: mass imbalanced gas in a harmonic trap
Example 1: gas confined in a tube
The solution of cold atomic gas in an external potential in the form of a tube. Gas with is in BCS regime with ak_F=-0.9. In the calculation, we assumed translation symmetry along the z direction, and st-wslda-2d was used. In the computation, double arithmetic is utilized.
The graph below shows the density distribution (left) and the absolute value of delta (right) for the converged solution.
Example 2: vortex solution within BdG
The solution represents a vortex confined in a tube. The conditions are the same as for Example 1. To speed up convergence, we start from the state provided in Example 1. In the computation, double complex arithmetic is utilized.
The graph below shows the density distribution (left) and the absolute value of delta (right) for the converged solution. By arrows, currents are plotted.

Example 3: mass imbalanced gas in a harmonic trap
This example is motivated by work arXiv:1909.03424. Namely, let us consider a gas of:
- component a:
{}^{161}\textrm{Dy}, - component b:
{}^{40}\textrm{K},
confined in a harmonic trap:
V_{a,b}(x)=\dfrac{m_a \omega_a^2 x^2}{2}
where trapping frequencies of both components are different (in the example we use, according to arXiv:1909.03424 \omega_a/\omega_b=120/430). In addition N_{\textrm{Dy}} / N_{\textrm{K}} = 20000/8000. In the calculations BdG functional is used, with the scattering length exceeding other length scales (in the example a=100). In the calculation, st-wslda-1d was used.
The graph below shows density distribution along x axis of {}^{161}\textrm{Dy} (red) and {}^{40}\textrm{K} (blue). The phase separation phenomenon is visible.

