|
|
|
# Target
|
|
|
|
Simulate quasi-1D Josephson dynamics. Initialize the simulation of a unitary Fermi gas confined in an external potential
|
|
|
|
```math
|
|
|
|
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:
|
|
|
|
```bash
|
|
|
|
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:
|
|
|
|
```c
|
|
|
|
/**
|
|
|
|
* 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
|
|
|
|
// ...
|
|
|
|
``` |
|
|
|
\ No newline at end of file |