st-wslda-examples.md

+[[_TOC_]]
+
 # 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](uploads/b0f9c475f5dbd44640310b713ecb5ea5/tube_predefines.h)
@@ -17,4 +19,20 @@ The solution representing a vortex confined in a tube. The conditions are the sa
 * [output](uploads/633a07650caee1345b0a1fbdd9513d51/vortex.out)
 
 The graph below shows density distribution (left) and the absolute value of delta (right) for the converged solution. By arrows currents are plotted. 
-![example-st-2](uploads/f896b92eeaa489872600109615ef066b/example-st-2.png)
\ No newline at end of file
+![example-st-2](uploads/f896b92eeaa489872600109615ef066b/example-st-2.png)
+
+# Example 3: Mass imbalanced gas in a harmonic trap  
+This example is motivated by work [arXiv:1909.03424](https://arxiv.org/abs/1909.03424).
+Namely, let us consider gas of:  
+* component *a*: $`{}^{161}\textrm{Dy}`$,
+* component *b*: $`{}^{40}\textrm{K}`$,  
+
+confined in harmonic trap:  
+```math
+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](https://arxiv.org/abs/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](uploads/8bed3548acc9dd04a2ccf54a35687ac7/DyK_predefines.h) 
+* [problem-definition.h](uploads/475df575833b5d4e92acd77b30992930/DyK_problem-definition.h)
+* [logger.h](uploads/8ac2b1262a442fa45092727564919f94/DyK_logger.h)
+* [input.txt](uploads/287f965385d1da19e7d84f007ec26a0e/DyK_input.txt)
\ No newline at end of file