• 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 of the form of a tube. Gas with is in BCS regime with ak_F=-0.9. In the calculation, we assumed translation symmetry along z direction and st-wslda-2d was used. In the computation, double arithmetic is utilized.

• predefines.h
• problem-definition.h
• input.txt
• output

The graph below shows density distribution (left) and the absolute value of delta (right) for the converged solution.

Example 2: vortex solution within BdG

The solution representing a vortex confined in a tube. The conditions are the same as for Example 1. To speed up the convergence process we start from the state provided by Example 1. In the computation, double complex arithmetic is utilized.

• predefines.h
• problem-definition.h
• input.txt
• output

The graph below shows 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 gas of:

• component a: {}^{161}\textrm{Dy},
• component b: {}^{40}\textrm{K},

confined in harmonic trap:

V_{a,b}(x)=\dfrac{m_a \omega_a^2 x^2}{2}

where traping 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.

• predefines.h
• problem-definition.h
• logger.h
• input.txt
• DyK.out

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.