|
*This page needs update*
|
|
*This page needs update*
|
|
|
|
|
|
In the W-SLDA toolkit it is possible to simulate "*strict 2d*" cases without any restriction for the z-direction. This can be done by simply choosing the size of the simulation box as Nx*Ny*1. In oder to do that one needs to modifiy the **predefines.h**:
|
|
In the W-SLDA toolkit, it is possible to simulate "*strict 2d*" cases without any restriction for the z-direction. This can be done by simply choosing the simulation box size as Nx*Ny*1. To do that, one needs to modify the **predefines.h**:
|
|
```c
|
|
```c
|
|
/**
|
|
/**
|
|
* Define lattice size and lattice spacing
|
|
* Define lattice size and lattice spacing
|
... | @@ -11,7 +11,7 @@ In the W-SLDA toolkit it is possible to simulate "*strict 2d*" cases without any |
... | @@ -11,7 +11,7 @@ In the W-SLDA toolkit it is possible to simulate "*strict 2d*" cases without any |
|
```
|
|
```
|
|
### Setting the Correct Fermi Momentum
|
|
### Setting the Correct Fermi Momentum
|
|
|
|
|
|
In the toolkit the default Fermi momentum is defined for the 3d case. If you are using strict 2d or 1d cases it is highly suggested to set the value of the Fermi momentum in the input file as **referencekF**.
|
|
The default Fermi momentum is defined in the toolkit for the 3d case. If you are using strict 2d or 1d cases, it is highly suggested to set the value of the Fermi momentum in the input file as **referencekF**.
|
|
The Fermi momentum for the 2d case:
|
|
The Fermi momentum for the 2d case:
|
|
* $`k_F = \sqrt{2\pi n}`$
|
|
* $`k_F = \sqrt{2\pi n}`$
|
|
|
|
|
... | @@ -22,9 +22,9 @@ where $`n`$ is the total density (spin up + spin down) of the gas. |
... | @@ -22,9 +22,9 @@ where $`n`$ is the total density (spin up + spin down) of the gas. |
|
|
|
|
|
# Choosing the Functional
|
|
# Choosing the Functional
|
|
|
|
|
|
The SLDA (and consequently ASLDA) functional has a self-energy term which is designed for 3d case. For the strcit 2d case we suggest **to use the BDG functional**. This can be done by choosing the BDG flag in the predefines.h file and requires the scattering length (*aBdG*) to be defined in the input file.
|
|
The SLDA (and consequently ASLDA) functional has a self-energy term, which is designed for 3d cases. For the strict 2d case, we suggest **using the BDG functional**. This can be done by choosing the BDG flag in the predefines.h file and requires the scattering length (*aBdG*) to be defined in the input file.
|
|
|
|
|
|
An alternative way to calculate the pairing interaction is to set a fixed coupling constant. This can be done by modifiying the pairing potential:
|
|
An alternative way to calculate the pairing interaction is to set a fixed coupling constant. This can be done by modifying the pairing potential:
|
|
|
|
|
|
```c
|
|
```c
|
|
void modify_potentials(int it, wslda_density h_densities, wslda_potential h_potentials, double *params, size_t extra_data_size, void *extra_data)
|
|
void modify_potentials(int it, wslda_density h_densities, wslda_potential h_potentials, double *params, size_t extra_data_size, void *extra_data)
|
... | @@ -53,4 +53,4 @@ params0 -2.30 # g_eff |
... | @@ -53,4 +53,4 @@ params0 -2.30 # g_eff |
|
|
|
|
|
```
|
|
```
|
|
|
|
|
|
By changing the value of this parameter one can set the magnitude of the pairing field to a desired value. |
|
By changing the value of this parameter, one can set the magnitude of the pairing field to a desired value. |
|
\ No newline at end of file |
|
\ No newline at end of file |