# 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