Target

Simulate quasi-1D Josephson dynamics. Initialize the simulation of a unitary Fermi gas confined in an external potential

V^{(ext)}_{\sigma}(x) = \frac{1}{2}\omega_x^2 x^2 + V_0e^{-2x^2/w^2} + k x,

where the first term corresponds to harmonic confinement, the second describes the potential barrier separating the cloud into two reservoirs, and the last introduces a small linear tilt to generate a population imbalance between the left and right reservoirs.

Step 1: Creating working folders

We start by copying template folders for static (st) and time-dependent (td) simulations:

cp -r $WSLDA/st-project-template ./st-jj-example
cp -r $WSLDA/td-project-template ./td-jj-example

Step 2: Editing files for the static calculations

Edit predefines.h

Select the lattice size and functional, and inspect all other tags that can be relevant for your problem:

/**
 * Define lattice size and lattice spacing.
 * */
#define NX (256)
#define NY (16)
#define NZ (16)
// ...
#define FUNCTIONAL SLDAE
/**
 * Sets the effective mass to be equal 1.
 * */
#define SLDA_FORCE_A1
// ...