|
|
|
# Documentation - Wiki:
|
|
|
|
# https://gitlab.fizyka.pw.edu.pl/wtools/wslda/-/wikis/home
|
|
|
|
|
|
|
|
# -------------- USER DEFINED PARAMETERS --------------
|
|
|
|
# See: Wiki -> User defined parameters
|
|
|
|
# Data flow: [Read params from input file] -> [execute process_params( )] -> [pass params to functions]
|
|
|
|
params[1] = 0.80 # Thomas-Fermi radius for x direction, in units of LX/2
|
|
|
|
params[2] = 0.90 # Smoothing parameter x1 for harmonic_oscillator_smooth_edges, in units of L/2
|
|
|
|
params[3] = 0.97 # Smoothing parameter x2 for harmonic_oscillator_smooth_edges, in units of L/2
|
|
|
|
params[4] = 0.24 # height of the barrier, in eF units
|
|
|
|
params[5] = 5.0 # width of barrier, in kF^1 units
|
|
|
|
params[6] = 0.0 # no tilt
|
|
|
|
|
|
|
|
# ------------------- INITIALIZATION ------------------
|
|
|
|
# See: Wiki -> Initialization of the solver
|
|
|
|
inittype 1 # 0 - create uniform solution and evolve it,
|
|
|
|
# 10 - read uniform solution from file `inprefix`/uniform.solution
|
|
|
|
#
|
|
|
|
# 1 - start from st-wslda-1d solution, inprefix points to folder with s1dpca binary files
|
|
|
|
# 2 - start from st-wslda-2d solution, inprefix points to folder with s2dpca binary files
|
|
|
|
# 3 - start from st-wslda-3d solution, inprefix points to folder with s3dpca binary files
|
|
|
|
#
|
|
|
|
# 5 - start from td-wslda checkpoint, inprefix points to folder with checkpoint binary files
|
|
|
|
|
|
|
|
# ------------------- INPUT/OUTPUT ------------------
|
|
|
|
outprefix jj-r2 # all output files will start with this prefix
|
|
|
|
inprefix ../st-jj-example/jj-r1 # prefix of calculation you need to use as input (when inittype!=0)
|
|
|
|
overwrite 1 # 1 - allow for overwriting existing files, 0 - do not allow for overwriting existing files
|
|
|
|
iogroups 8 # number of IO groups used for parallel reading of data
|
|
|
|
# it must be the same as the value used for static calculations.
|
|
|
|
checkpoint 1 # do checkpoint at the end of calculation: 1-yes, 0-no
|
|
|
|
checkperiod 0 # periodicity of the checkpoint, meaningful only if checkpoint=1
|
|
|
|
# if checkperiod>0 then after each checkperiod measurements, the checkpoint will be conducted
|
|
|
|
# if checkperiod=0 the option is deactivated (default)
|
|
|
|
walltime 11.9 # after this time force the checkpoint [hr]
|
|
|
|
|
|
|
|
# list of variables for writing
|
|
|
|
# all = rho delta j nu tau V V_ext delta_ext velocity_ext alpha A
|
|
|
|
# default = rho delta j
|
|
|
|
writevar default V_ext
|
|
|
|
writeprec d # precision for data writing: d-double(default), f-float.
|
|
|
|
|
|
|
|
# ----------------- SCATTERING LENGTH ---------------
|
|
|
|
sclgth -10000.0 # scattering length in units of lattice spacing
|
|
|
|
# meaningful only for FUNCTIONAL=BDG,SLDAE
|
|
|
|
# in case of FUNCTIONAL=(A)SLDA it is set automatically to infinity
|
|
|
|
#akF -1.0 # scattering length times Fermi wave vector
|
|
|
|
# this tag modifies the value of tag sclgth via process_params(...) function
|
|
|
|
# by default is inactive (akF=0)
|
|
|
|
|
|
|
|
# ------------------- TIME EVOLUTION ----------------
|
|
|
|
dt 0.0035 # integration time step, in units of eF^-1
|
|
|
|
measurements 200 # number of requested measurements
|
|
|
|
timesteps 1430 # number of time steps between each measurement
|
|
|
|
# measurements are written with time resolution: timesteps*dt
|
|
|
|
# the total trajectory length is: measurements*timesteps*dt
|
|
|
|
selfstart 1 # use Taylor expansion for first 5 steps?
|
|
|
|
# See: Wiki -> Starting ABM algorithm
|
|
|
|
|
|
|
|
# ------------------ GPUS THREADING SET --------------
|
|
|
|
# nthreads 512 # number of parallel threads executing by each GPU, default=512
|
|
|
|
# batch 1000000 # maximal number of wave-functions processed in parallel by each card, default=1000000
|
|
|
|
# gpuspernode 1 # number of GPUs per node (resource set), default=GPUS_PER_NODE (in machine.h)
|
|
|
|
|
|
|
|
# -------------- HIGH K-WAVES FILTER ----------------
|
|
|
|
# See: Wiki -> Stabilization of the time-dependent code
|
|
|
|
hkf_mode 1 # 0 - no noise filtering (default)
|
|
|
|
# 1 - noise filtering for mean-fields only
|
|
|
|
# 2 - noise filtering of mean-fields and pairing field
|
|
|
|
hkf_mu 0.9 # mu parameter of the Fermi-Dirac (filtering) function, in Ec units, default=0.9
|
|
|
|
hkf_T 0.02 # T parameter of the Fermi-Dirac (filtering) function, in Ec units, default=0.02
|
|
|
|
|
|
|
|
# ---------------- CONSERVATION LAWS -----------------
|
|
|
|
# See: Wiki -> Monitoring of conservation laws
|
|
|
|
# # -> Energy conservation monitoring
|
|
|
|
Econservation_start 0 # Start to monitor energy conservation from time*eF, default=1e12 (infinity)
|
|
|
|
Econservation_stop 10000 # Stop to monitor energy conservation at time*eF, default=1e12 (infinity)
|
|
|
|
Econservation_tol 0.05 # if |[E(t)-E(start_t)]/E(start_t)|>tol then the code will terminate, default=0.05
|
|
|
|
# # -> Particle number conservation monitoring
|
|
|
|
Nconservation_start 0 # Start to monitor total particle number conservation from time*eF, default=0
|
|
|
|
Nconservation_stop 10000 # Stop to monitor total particle number conservation at time*eF, default=1e12 (infinity)
|
|
|
|
Nconservation_tol 0.05 # if |[N(t)-N(start_t)]/N(start_t)|>tol then the code will terminate, default=0.05
|
|
|
|
|
|
|
|
# --------------- QUANTUM FRICTION ------------------
|
|
|
|
# See: Wiki -> Quantum friction
|
|
|
|
# qfalpha 0.0 # alpha parameter for quantum friction term (cooling in the normal channel), default=0
|
|
|
|
# qfbeta 0.0 # beta parameter for quantum friction term (cooling in the pairing channel), default=0
|
|
|
|
# qfgamma 0.0 # particle control parameter (cooling in the pairing channel), default=0
|
|
|
|
# qfNreq 0.0 # requested particle number, meaningful only if qfgamma>0
|
|
|
|
# qfstart 0.0 # start time for evolving with quantum friction [eF], default=0
|
|
|
|
# qfstop 0.0 # stop time for evolving with quantum friction [eF], default=infinity
|
|
|
|
# qfswitch 0.0 # the friction will be activated and deactivated gradually over this period of time [eF], default=1
|
|
|
|
|
|
|
|
# ---------------- SUBSET TRACKING ------------------
|
|
|
|
# See: Wiki -> Tracking of selected states
|
|
|
|
# observables (densities, currents) arising from state in the energy interval En in [subsetMinEn,subsetMaxEn]
|
|
|
|
# will be computed, if subsetMinEn=subsetMaxEn=0 the functionality is disables
|
|
|
|
# subsetMinEn 0.0 # in eF units, deafault=0.0
|
|
|
|
# subsetMaxEn 0.0 # in eF units, deafault=0.0
|
|
|
|
# subsetShiftDmu 1 # if 1 then apply extra shift of quasiparticle energies by (mu_a-mu_b)/2, default=0
|
|
|
|
|
|
|
|
# ----------- IN CASE OF: inittype=0 ----------------
|
|
|
|
# See: Wiki -> Initialization of the solver
|
|
|
|
# init0Na 28.0 # requested number of particles of type a, default=Na
|
|
|
|
# init0Nb 28.0 # requested number of particles of type b, default=Nb
|
|
|
|
init0na 0.01689 # ...alternatively, you can request the desired density for particles of type a
|
|
|
|
init0nb 0.01689 # ...alternatively, you can request the desired density for particles of type b
|
|
|
|
init0muchange 0.1 # Change rate of chemical potential
|
|
|
|
init0Tstart 0.0 # Start temperature, in units of eF, default 0.2
|
|
|
|
init0Tstop 0.0 # Stop temperature, in units of eF, default 0.05
|
|
|
|
init0DeltaT 0.01 # change of temperature in units of eF, default 0.01
|
|
|
|
init0eps 1.0e-9 # epsilon for convergence, default 1.0e-6
|
|
|
|
init0scmix 0.5 # mixing parameter in self-consistent process, default 0.25
|
|
|
|
init0maxiter 2000 # maximum number of iterations, default 100000
|
|
|
|
# init0kc 3.1415 # momentum cut-off used for uniform solver, default=pi/max(DX,DY,DZ)
|
|
|
|
init0debug 1 # debug level, default=0 (no debug info),
|
|
|
|
# init0save 1 # save solution |
