INPUT_PARAMETERS
#Parameters (1.System)
suffix                         ABACUS #the name of main output directory
ntype                          6 #atom species number
calculation                    scf #test; scf; relax; nscf; get_wf; get_pchg
esolver_type                   ksdft #the energy solver: ksdft, sdft, ofdft, tddft, lj, dp, ks-lr, lr
symmetry                       0 #the control of symmetry
symmetry_prec                  1e-06 #accuracy for symmetry
symmetry_autoclose             1 #whether to close symmetry automatically when error occurs in symmetry analysis
cal_stress                     1 #calculate the stress or not
cal_force                      1 #if calculate the force at the end of the electronic iteration
kpar                           1 #devide all processors into kpar groups and k points will be distributed among
bndpar                         1 #devide all processors into bndpar groups and bands will be distributed among each group
latname                        none #the name of lattice name
ecutwfc                        100 #energy cutoff for wave functions
ecutrho                        400 #energy cutoff for charge density and potential
nx                             0 #number of points along x axis for FFT grid
ny                             0 #number of points along y axis for FFT grid
nz                             0 #number of points along z axis for FFT grid
ndx                            0 #number of points along x axis for FFT smooth grid
ndy                            0 #number of points along y axis for FFT smooth grid
ndz                            0 #number of points along z axis for FFT smooth grid
cell_factor                    1.2 #used in the construction of the pseudopotential tables
erf_ecut                       0 #the value of the constant energy cutoff
erf_height                     0 #the height of the energy step for reciprocal vectors
erf_sigma                      0.1 #the width of the energy step for reciprocal vectors
fft_mode                       0 #mode of FFTW
init_wfc                       atomic #start wave functions are from 'atomic', 'atomic+random', 'random' or
pw_seed                        0 #random seed for initializing wave functions
init_chg                       atomic #start charge is from 'atomic' or file
dm_to_rho                      0 #reads dmr in npz format and calculates electron density
chg_extrap                     atomic #atomic; first-order; second-order; dm:coefficients of SIA
init_vel                       0 #read velocity from STRU or not
stru_file                      STRU #the filename of file containing atom positions
kpoint_file                    KPT #the name of file containing k points
pseudo_dir                      #the directory containing pseudo files
orbital_dir                     #the directory containing orbital files
read_file_dir                  OUT.ABACUS/ #directory of files for reading
restart_load                   0 #restart from disk
wannier_card                   none #input card for wannier functions
mem_saver                      0 #Only for nscf calculations. if set to 1, then a memory saving technique will be used for many k point calculations.
diago_proc                     16 #the number of procs used to do diagonalization
nbspline                       -1 #the order of B-spline basis
kspacing                       0.1 0.1 0.1  #unit in 1/bohr, should be > 0, default is 0 which means read KPT file
min_dist_coef                  0.2 #factor related to the allowed minimum distance between two atoms
device                         cpu #the computing device for ABACUS
precision                      double #the computing precision for ABACUS

#Parameters (2.Electronic structure)
ks_solver                      genelpa #cg; dav; lapack; genelpa; elpa; scalapack_gvx; cusolver
basis_type                     lcao #PW; LCAO in pw; LCAO
use_paw                        0 #whether to use PAW in pw calculation
nbands                         0 #number of bands
nbands_mul                     1.2 #value to multiply the number of bands
nelec                          0 #input number of electrons
nelec_delta                    0 #change in the number of total electrons
nupdown                        0 #the difference number of electrons between spin-up and spin-down
dft_functional                 GGA_X_RPBE+GGA_C_PBE #exchange correlation functional
xc_temperature                 0 #temperature for finite temperature functionals
pseudo_rcut                    15 #default #exchange correlation functional
pseudo_mesh                    0 #0: use our own mesh to do radial renormalization; 1: use mesh as in QE
nspin                          2 #1: single spin; 2: up and down spin; 4: noncollinear spin
pw_diag_nmax                   50 #max iteration number for cg
pw_diag_thr                    0.01 #threshold for eigenvalues is cg electron iterations
diago_smooth_ethr              0 #smooth ethr for iter methods
use_k_continuity               0 #whether to use k-point continuity for initializing wave functions
pw_diag_ndim                   4 #dimension of workspace for Davidson diagonalization
diago_cg_prec                  1 #diago_cg_prec
smearing_method                gaussian #type of smearing_method: gauss; fd; fixed; mp; mp2; mv
smearing_sigma                 0.01 #energy range for smearing
mixing_type                    broyden #plain; pulay; broyden
mixing_beta                    0.2 #mixing parameter: 0 means no new charge
mixing_ndim                    8 #mixing dimension in pulay or broyden
mixing_restart                 0 #threshold to restart mixing during SCF
mixing_gg0                     1 #mixing parameter in kerker
mixing_beta_mag                0.8 #mixing parameter for magnetic density
mixing_gg0_mag                 0 #mixing parameter in kerker
mixing_gg0_min                 0.1 #the minimum kerker coefficient
mixing_angle                   -10 #angle mixing parameter for non-colinear calculations
mixing_tau                     0 #whether to mix tau in mGGA calculation
mixing_dftu                    0 #whether to mix locale in DFT+U calculation
mixing_dmr                     0 #whether to mix real-space density matrix
gamma_only                     0 #Only for localized orbitals set and gamma point. If set to 1, a fast algorithm is used
scf_nmax                       100 #number of electron iterations
scf_thr                        1e-07 #charge density error
scf_ene_thr                    -1 #total energy error threshold
scf_os_stop                    0 #whether to stop scf when oscillation is detected
scf_os_thr                     -0.01 #charge density threshold for oscillation
scf_os_ndim                    8 #number of old iterations used for oscillation detection
sc_os_ndim                     5 #number of old iterations used for oscillation detection, for Spin-Constrained DFT
scf_thr_type                   2 #type of the criterion of scf_thr, 1: reci drho for pw, 2: real drho for lcao
lspinorb                       0 #consider the spin-orbit interaction
noncolin                       0 #using non-collinear-spin
soc_lambda                     1 #The fraction of averaged SOC pseudopotential is given by (1-soc_lambda)

#Parameters (3.LCAO)
nb2d                           0 #matrix 2d division
lmaxmax                        2 #maximum of l channels used
lcao_ecut                      100 #energy cutoff for LCAO
lcao_dk                        0.01 #delta k for 1D integration in LCAO
lcao_dr                        0.01 #delta r for 1D integration in LCAO
lcao_rmax                      30 #max R for 1D two-center integration table
search_radius                  -1 #input search radius (Bohr)
search_pbc                     1 #input periodic boundary condition
bx                             0 #division of an element grid in FFT grid along x
by                             0 #division of an element grid in FFT grid along y
bz                             0 #division of an element grid in FFT grid along z
elpa_num_thread                -1 #Number of threads need to use in elpa
num_stream                     4 #the nstream in compute the LCAO with CUDA
bessel_nao_ecut                100.000000 #energy cutoff for spherical bessel functions(Ry)
bessel_nao_tolerence           1e-12 #tolerence for spherical bessel root
bessel_nao_rcut                 #radial cutoff for spherical bessel functions(a.u.)
bessel_nao_smooth              1 #spherical bessel smooth or not
bessel_nao_sigma               0.1 #spherical bessel smearing_sigma

#Parameters (4.Relaxation)
relax_method                   cg #cg; bfgs; sd; cg; cg_bfgs;
relax_new                      1 #whether to use the new relaxation method
relax                          0 #allow relaxation along the specific direction
relax_scale_force              0.5 #controls the size of the first CG step if relax_new is true
relax_nmax                     1 #number of ion iteration steps
relax_cg_thr                   0.5 #threshold for switching from cg to bfgs, unit: eV/Angstrom
force_thr                      0.001 #force threshold, unit: Ry/Bohr
force_thr_ev                   0.0257112 #force threshold, unit: eV/Angstrom
force_thr_ev2                  0 #force invalid threshold, unit: eV/Angstrom
stress_thr                     0.5 #stress threshold
press1                         0 #target pressure, unit: KBar
press2                         0 #target pressure, unit: KBar
press3                         0 #target pressure, unit: KBar
relax_bfgs_w1                  0.01 #wolfe condition 1 for bfgs
relax_bfgs_w2                  0.5 #wolfe condition 2 for bfgs
relax_bfgs_rmax                0.2 #maximal trust radius, unit: Bohr
relax_bfgs_rmin                1e-05 #minimal trust radius, unit: Bohr
relax_bfgs_init                0.5 #initial trust radius, unit: Bohr
fixed_axes                     None #which axes are fixed
fixed_ibrav                    0 #whether to preseve lattice type during relaxation
fixed_atoms                    0 #whether to preseve direct coordinates of atoms during relaxation

#Parameters (5.Molecular dynamics)
md_type                        nvt #choose ensemble
md_thermostat                  nhc #choose thermostat
md_nstep                       10 #md steps
md_dt                          1 #time step
md_tchain                      1 #number of Nose-Hoover chains
md_tfirst                      -1 #temperature first
md_tlast                       -1 #temperature last
md_dumpfreq                    1 #The period to dump MD information
md_restartfreq                 5 #The period to output MD restart information
md_seed                        -1 #random seed for MD
md_prec_level                  0 #precision level for vc-md
ref_cell_factor                1 #construct a reference cell bigger than the initial cell
md_restart                     0 #whether restart
lj_rule                        2 #combination rules used to construct the parameter matrix for LJ potential
lj_eshift                      0 #whether to use energy shift for LJ potential
lj_rcut                         #cutoff radius of LJ potential
lj_epsilon                      #the value of epsilon for LJ potential
lj_sigma                        #the value of sigma for LJ potential
pot_file                       graph.pb #the filename of potential files for CMD such as DP
dp_rescaling                   1 #rescaling factor for dp potential
dp_fparam                       #the frame parameter for dp potential
dp_aparam                       #the atomic parameter for dp potential
msst_direction                 2 #the direction of shock wave
msst_vel                       0 #the velocity of shock wave
msst_vis                       0 #artificial viscosity
msst_tscale                    0.01 #reduction in initial temperature
msst_qmass                     1 #mass of thermostat
md_tfreq                       0 #oscillation frequency, used to determine qmass of NHC
md_damp                        1 #damping parameter (time units) used to add force in Langevin method
md_nraise                      1 #parameters used when md_type=nvt
cal_syns                       0 #calculate asynchronous overlap matrix to output for Hefei-NAMD
dmax                           0.01 #maximum displacement of all atoms in one step (bohr)
md_tolerance                   100 #tolerance for velocity rescaling (K)
md_pmode                       iso #NPT ensemble mode: iso, aniso, tri
md_pcouple                     xyz #whether couple different components: xyz, xy, yz, xz, none
md_pchain                      1 #num of thermostats coupled with barostat
md_pfirst                      -1 #initial target pressure
md_plast                       -1 #final target pressure
md_pfreq                       0 #oscillation frequency, used to determine qmass of thermostats coupled with barostat
dump_force                     1 #output atomic forces into the file MD_dump or not
dump_vel                       1 #output atomic velocities into the file MD_dump or not
dump_virial                    1 #output lattice virial into the file MD_dump or not

#Parameters (6.orbital free density functional theory)
of_kinetic                     wt #kinetic energy functional, such as tf, vw, wt
of_method                      tn #optimization method used in OFDFT, including cg1, cg2, tn (default)
of_conv                        energy #the convergence criterion, potential, energy (default), or both
of_tole                        1e-06 #tolerance of the energy change (in Ry) for determining the convergence, default=2e-6 Ry
of_tolp                        1e-05 #tolerance of potential for determining the convergence, default=1e-5 in a.u.
of_tf_weight                   1 #weight of TF KEDF
of_vw_weight                   1 #weight of vW KEDF
of_wt_alpha                    0.833333 #parameter alpha of WT KEDF
of_wt_beta                     0.833333 #parameter beta of WT KEDF
of_wt_rho0                     0 #the average density of system, used in WT KEDF, in Bohr^-3
of_hold_rho0                   0 #If set to 1, the rho0 will be fixed even if the volume of system has changed, it will be set to 1 automaticly if of_wt_rho0 is not zero
of_lkt_a                       1.3 #parameter a of LKT KEDF
of_full_pw                     1 #If set to 1, ecut will be ignored when collect planewaves, so that all planewaves will be used
of_full_pw_dim                 0 #If of_full_pw = true, dimention of FFT is testricted to be (0) either odd or even; (1) odd only; (2) even only
of_read_kernel                 0 #If set to 1, the kernel of WT KEDF will be filled from file of_kernel_file, not from formula. Only usable for WT KEDF
of_kernel_file                 WTkernel.txt #The name of WT kernel file.

#Parameters (7.Stochastic DFT)
method_sto                     2 #1: slow and save memory, 2: fast and waste memory
npart_sto                      1 #Reduce memory when calculating Stochastic DOS
nbands_sto                     256 #number of stochstic orbitals
nche_sto                       100 #Chebyshev expansion orders
emin_sto                       0 #trial energy to guess the lower bound of eigen energies of the Hamitonian operator
emax_sto                       0 #trial energy to guess the upper bound of eigen energies of the Hamitonian operator
seed_sto                       0 #the random seed to generate stochastic orbitals
initsto_ecut                   0 #maximum ecut to init stochastic bands
initsto_freq                   0 #frequency to generate new stochastic orbitals when running md

#Parameters (8.DeepKS)
deepks_out_labels              0 #>0 compute descriptor for deepks
deepks_scf                     0 #>0 add V_delta to Hamiltonian
deepks_equiv                   0 #whether to use equivariant version of DeePKS
deepks_bandgap                 0 #>0 for bandgap label
deepks_v_delta                 0 #>0 for v_delta label. when output, 1 for v_delta_precalc, 2 for phialpha and grad_evdm ( can save memory )
deepks_out_unittest            0 #if set 1, prints intermediate quantities that shall be used for making unit test
deepks_model                   None #file dir of traced pytorch model: 'model.ptg
bessel_descriptor_lmax         2 #lmax used in generating spherical bessel functions
bessel_descriptor_ecut         100.000000 #energy cutoff for spherical bessel functions(Ry)
bessel_descriptor_tolerence    1e-12 #tolerence for spherical bessel root
bessel_descriptor_rcut         6 #radial cutoff for spherical bessel functions(a.u.)
bessel_descriptor_smooth       1 #spherical bessel smooth or not
bessel_descriptor_sigma        0.1 #sphereical bessel smearing_sigma

#Parameters (9.rt-tddft)
td_force_dt                    0.02 #time of force change
td_vext                        0 #add extern potential or not
td_vext_dire                   1 #extern potential direction
init_vecpot_file               0 #init vector potential through file or not
td_print_eij                   -1 #print eij or not
td_edm                         0 #the method to calculate the energy density matrix
td_propagator                  0 #method of propagator
td_stype                       0 #type of electric field in space domain
td_ttype                       0 #type of electric field in time domain
td_tstart                      1 # number of steps where electric field starts
td_tend                        1000 #number of steps where electric field ends
td_lcut1                       0.05 #cut1 of interval in length gauge
td_lcut2                       0.95 #cut2 of interval in length gauge
td_gauss_freq                  22.13 #frequency (freq) of Gauss type electric field
td_gauss_phase                 0.0 #phase of Gauss type electric field
td_gauss_sigma                 30.0 #sigma of Gauss type electric field
td_gauss_t0                    100.0 #step number of time center (t0) of Gauss type electric field
td_gauss_amp                   0.25 #amplitude of Gauss type electric field
td_trape_freq                  1.60 #frequency of Trapezoid type electric field
td_trape_phase                 0.0 #phase of Trapezoid type electric field
td_trape_t1                    1875.0 #t1 of Trapezoid type electric field
td_trape_t2                    5625.0 #t2 of Trapezoid type electric field
td_trape_t3                    7500.0 #t3 of Trapezoid type electric field
td_trape_amp                   2.74 #amplitude of Trapezoid type electric field
td_trigo_freq1                 1.164656 #frequency 1 of Trigonometric type electric field
td_trigo_freq2                 0.029116 #frequency 2 of Trigonometric type electric field
td_trigo_phase1                0.0 #phase 1 of Trigonometric type electric field
td_trigo_phase2                0.0 #phase 2 of Trigonometric type electric field
td_trigo_amp                   2.74 #amplitude of Trigonometric type electric field
td_heavi_t0                    100.0 #t0 of Heaviside type electric field
td_heavi_amp                   1.0 #amplitude of Heaviside type electric field
ocp                            0 #change occupation or not
ocp_set                         #set occupation

#Parameters (10.lr-tddft)
lr_nstates                     1 #the number of 2-particle states to be solved
nocc                           0 #the number of occupied orbitals to form the 2-particle basis ( <= nelec/2)
nvirt                          1 #the number of virtual orbitals to form the 2-particle basis (nocc + nvirt <= nbands)
xc_kernel                      LDA #exchange correlation (XC) kernel for LR-TDDFT
lr_init_xc_kernel              default  #The method to initalize the xc kernel
lr_solver                      dav #the eigensolver for LR-TDDFT
lr_thr                         0.01 #convergence threshold of the LR-TDDFT eigensolver
out_wfc_lr                     0 #whether to output the eigenvectors (excitation amplitudes) in the particle-hole basis
lr_unrestricted                0 #Whether to use unrestricted construction for LR-TDDFT
abs_wavelen_range              0 0  #the range of wavelength(nm) to output the absorption spectrum 
abs_gauge                      length #whether to use length or velocity gauge to calculate the absorption spectrum in LR-TDDFT
abs_broadening                 0.01 #the broadening (eta) for LR-TDDFT absorption spectrum

#Parameters (11.Output)
out_stru                       0 #output the structure files after each ion step
out_freq_elec                  100 #the frequency of electronic iter to output charge density and wavefunction 
out_freq_ion                   0 #the frequency ( >= 0 ) of ionic step to output charge density and wavefunction. 0: output only when ion steps are finished
out_chg                        0 3  #> 0 output charge density for selected electron steps, second parameter controls the precision, default is 3.
out_pot                        0 #output realspace potential
out_wfc_pw                     0 #output wave functions
out_wfc_r                      0 #output wave functions in realspace
printe                         100 #Print out energy for each band for every printe steps
out_band                       0 8  #output energy and band structure (with precision 8)
out_dos                        0 #output energy and dos
out_mul                        0 #mulliken charge or not
out_proj_band                  0 #output projected band structure
out_level                      ie #ie(for electrons); i(for ions);
out_dm                         0 #>0 output density matrix
out_dm1                        0 #>0 output density matrix (multi-k points)
out_bandgap                    0 #if true, print out bandgap
out_mat_hs                     0 8  #output H and S matrix (with precision 8)
out_mat_tk                     0 8  #output T(k)
out_mat_hs2                    0 #output H(R) and S(R) matrix
out_mat_dh                     0 #output of derivative of H(R) matrix
out_mat_xc                     0 #output exchange-correlation matrix in KS-orbital representation
out_eband_terms                0 #output the band energy terms separately
out_hr_npz                     0 #output hr(I0,JR) submatrices in npz format
out_dm_npz                     0 #output dmr(I0,JR) submatrices in npz format
out_interval                   1 #interval for printing H(R) and S(R) matrix during MD
out_app_flag                   1 #whether output r(R), H(R), S(R), T(R), and dH(R) matrices in an append manner during MD
out_ndigits                    8 #the length of decimal part of output data
out_mat_t                      0 #output T(R) matrix
out_element_info               0 #output (projected) wavefunction of each element
out_mat_r                      0 #output r(R) matrix
out_wfc_lcao                   0 #ouput LCAO wave functions, 0, no output 1: text, 2: binary
out_dipole                     0 #output dipole or not
out_efield                     0 #output dipole or not
out_current                    0 #output current or not
out_current_k                  0 #output current for each k
out_vecpot                     0 #output TDDFT vector potential or not
restart_save                   0 #print to disk every step for restart
rpa                            0 #true:generate output files used in rpa calculation; false:(default)
nbands_istate                  5 #number of bands around Fermi level for get_wf and get_pchg calulation
bands_to_print                  #specify the bands to be calculated for the partial (band-decomposed) charge densities
out_pchg                        #specify the bands to be calculated for the partial (band-decomposed) charge densities
out_wfc_norm                    #specify the bands to be calculated for the norm of wavefunctions
out_wfc_re_im                   #specify the bands to be calculated for the real and imaginary parts of wavefunctions
if_separate_k                  0 #specify whether to write the partial charge densities for all k-points to individual files or merge them
out_elf                        0 3  #> 0 output electron localization function (ELF) for selected electron steps, second parameter controls the precision, default is 3.

#Parameters (12.Postprocess)
dos_emin_ev                    -15 #minimal range for dos
dos_emax_ev                    15 #maximal range for dos
dos_edelta_ev                  0.01 #delta energy for dos
dos_scale                      0.01 #scale dos range by
dos_sigma                      0.07 #gauss b coefficeinet(default=0.07)
dos_nche                       100 #orders of Chebyshev expansions for dos
cal_cond                       0 #calculate electronic conductivities
cond_che_thr                   1e-08 #control the error of Chebyshev expansions for conductivities
cond_dw                        0.1 #frequency interval for conductivities
cond_wcut                      10 #cutoff frequency (omega) for conductivities
cond_dt                        0.02 #t interval to integrate Onsager coefficiencies
cond_dtbatch                   0 #exp(iH*dt*cond_dtbatch) is expanded with Chebyshev expansion
cond_smear                     1 #Smearing method for conductivities
cond_fwhm                      0.4 #FWHM for conductivities
cond_nonlocal                  1 #Nonlocal effects for conductivities
berry_phase                    0 #calculate berry phase or not
gdir                           3 #calculate the polarization in the direction of the lattice vector
towannier90                    0 #use wannier90 code interface or not
nnkpfile                       seedname.nnkp #the wannier90 code nnkp file name
wannier_spin                   up #calculate spin in wannier90 code interface
wannier_method                 1 #different implementation methods under Lcao basis set
out_wannier_mmn                1 #output .mmn file or not
out_wannier_amn                1 #output .amn file or not
out_wannier_unk                0 #output UNK. file or not
out_wannier_eig                1 #output .eig file or not
out_wannier_wvfn_formatted     1 #output UNK. file in text format or in binary format

#Parameters (13.Model)
efield_flag                    0 #add electric field
dip_cor_flag                   0 #dipole correction
efield_dir                     2 #the direction of the electric field or dipole correction
efield_pos_max                 -1 #position of the maximum of the saw-like potential along crystal axis efield_dir
efield_pos_dec                 -1 #zone in the unit cell where the saw-like potential decreases
efield_amp                     0 #amplitude of the electric field
gate_flag                      0 #compensating charge or not
zgate                          0.5 #position of charged plate
block                          0 #add a block potential or not
block_down                     0.45 #low bound of the block
block_up                       0.55 #high bound of the block
block_height                   0.1 #height of the block
imp_sol                        0 #calculate implicit solvation correction or not
eb_k                           80 #the relative permittivity of the bulk solvent
tau                            1.0798e-05 #the effective surface tension parameter
sigma_k                        0.6 #the width of the diffuse cavity
nc_k                           0.00037 #the cut-off charge density

#Parameters (14.vdW Correction)
vdw_method                     d3_bj #the method of calculating vdw (none ; d2 ; d3_0 ; d3_bj
vdw_s6                         default #scale parameter of d2/d3_0/d3_bj
vdw_s8                         default #scale parameter of d3_0/d3_bj
vdw_a1                         default #damping parameter of d3_0/d3_bj
vdw_a2                         default #damping parameter of d3_bj
vdw_d                          20 #damping parameter of d2
vdw_abc                        0 #third-order term?
vdw_c6_file                    default #filename of C6
vdw_c6_unit                    Jnm6/mol #unit of C6, Jnm6/mol or eVA6
vdw_r0_file                    default #filename of R0
vdw_r0_unit                    A #unit of R0, A or Bohr
vdw_cutoff_type                radius #expression model of periodic structure, radius or period
vdw_cutoff_radius              95 #radius cutoff for periodic structure
vdw_radius_unit                Bohr #unit of radius cutoff for periodic structure
vdw_cn_thr                     40 #radius cutoff for cn
vdw_cn_thr_unit                Bohr #unit of cn_thr, Bohr or Angstrom
vdw_cutoff_period              3 3 3 #periods of periodic structure

#Parameters (15.exx)
exx_hybrid_alpha               0 #fraction of Fock exchange in hybrid functionals
exx_hse_omega                  0.11 #range-separation parameter in HSE functional
exx_separate_loop              1 #if 1, a two-step method is employed, else it will start with a GGA-Loop, and then Hybrid-Loop
exx_hybrid_step                100 #the maximal electronic iteration number in the evaluation of Fock exchange
exx_mixing_beta                1 #mixing_beta for outer-loop when exx_separate_loop=1
exx_lambda                     0.3 #used to compensate for divergence points at G=0 in the evaluation of Fock exchange using lcao_in_pw method
exx_real_number                0 #exx calculated in real or complex
exx_pca_threshold              0.0001 #threshold to screen on-site ABFs in exx
exx_c_threshold                0.0001 #threshold to screen C matrix in exx
exx_v_threshold                0.1 #threshold to screen C matrix in exx
exx_dm_threshold               0.0001 #threshold to screen density matrix in exx
exx_schwarz_threshold          0 #threshold to screen exx using Cauchy-Schwartz inequality
exx_cauchy_threshold           1e-07 #threshold to screen exx using Cauchy-Schwartz inequality
exx_c_grad_threshold           0.0001 #threshold to screen nabla C matrix in exx
exx_v_grad_threshold           0.1 #threshold to screen nabla V matrix in exx
exx_c_grad_r_threshold         0.0001 #threshold to screen nabla C matrix in exx
exx_v_grad_r_threshold         0.1 #threshold to screen nabla V matrix in exx
exx_cauchy_force_threshold     1e-07 #threshold to screen exx force using Cauchy-Schwartz inequality
exx_cauchy_stress_threshold    1e-07 #threshold to screen exx stress using Cauchy-Schwartz inequality
exx_ccp_rmesh_times            1 #how many times larger the radial mesh required for calculating Columb potential is to that of atomic orbitals
exx_distribute_type            htime #exx_distribute_type
exx_opt_orb_lmax               0 #the maximum l of the spherical Bessel functions for opt ABFs
exx_opt_orb_ecut               0 #the cut-off of plane wave expansion for opt ABFs
exx_opt_orb_tolerence          0 #the threshold when solving for the zeros of spherical Bessel functions for opt ABFs
exx_symmetry_realspace         0 #whether to reduce real-space sector in Hexx calculation
rpa_ccp_rmesh_times            10 #how many times larger the radial mesh required for calculating Columb potential is to that of atomic orbitals
out_ri_cv                      0 #Whether to output the coefficient tensor C and ABFs-representation Coulomb matrix V

#Parameters (16.dft+u)
dft_plus_u                     0 #DFT+U correction method
dft_plus_dmft                  0 #true:DFT+DMFT; false: standard DFT calcullation(default)
yukawa_lambda                  -1 #default:0.0
yukawa_potential               0 #default: false
uramping                       -1 #increasing U values during SCF
omc                            0 #the mode of occupation matrix control
onsite_radius                  3 #radius of the sphere for onsite projection (Bohr)
hubbard_u                      0 0 0 0 0 0  #Hubbard Coulomb interaction parameter U(ev)
orbital_corr                   -1 -1 -1 -1 -1 -1  #which correlated orbitals need corrected ; d:2 ,f:3, do not need correction:-1

#Parameters (17.non-collinear spin-constrained DFT)
sc_mag_switch                  0 #switch to control spin-constrained DFT
decay_grad_switch              0 #switch to control gradient break condition
sc_thr                         1e-06 #Convergence criterion of spin-constrained iteration (RMS) in uB
nsc                            100 #Maximal number of spin-constrained iteration
nsc_min                        2 #Minimum number of spin-constrained iteration
sc_scf_nmin                    2 #Minimum number of outer scf loop before initializing lambda loop
alpha_trial                    0.01 #Initial trial step size for lambda in eV/uB^2
sccut                          3 #Maximal step size for lambda in eV/uB
sc_drop_thr                    0.001 #Convergence criterion ratio of lambda iteration in Spin-constrained DFT
sc_scf_thr                     0.001 #Density error threshold for inner loop of spin-constrained SCF

#Parameters (18.Quasiatomic Orbital analysis)
qo_switch                      0 #switch to control quasiatomic orbital analysis
qo_basis                       szv #type of QO basis function: hydrogen: hydrogen-like basis, pswfc: read basis from pseudopotential
qo_thr                         1e-06 #accuracy for evaluating cutoff radius of QO basis function
qo_strategy                    all all all all all all  #strategy to generate generate radial orbitals
qo_screening_coeff             0.1 0.1 0.1 0.1 0.1 0.1  #rescale the shape of radial orbitals

#Parameters (19.PEXSI)
pexsi_npole                    40 #Number of poles in expansion
pexsi_inertia                  1 #Whether inertia counting is used at the very beginning of PEXSI process
pexsi_nmax                     80 #Maximum number of PEXSI iterations after each inertia counting procedure
pexsi_comm                     1 #Whether to construct PSelInv communication pattern
pexsi_storage                  1 #Storage space used by the Selected Inversion algorithm for symmetric matrices
pexsi_ordering                 0 #Ordering strategy for factorization and selected inversion
pexsi_row_ordering             1 #Row permutation strategy for factorization and selected inversion, 0: NoRowPerm, 1: LargeDiag
pexsi_nproc                    1 #Number of processors for parmetis
pexsi_symm                     1 #Matrix symmetry
pexsi_trans                    0 #Whether to transpose
pexsi_method                   1 #pole expansion method, 1: Cauchy Contour Integral, 2: Moussa optimized method
pexsi_nproc_pole               1 #Number of processes used by each pole
pexsi_temp                     0.015 #Temperature, in the same unit as H
pexsi_gap                      0 #Spectral gap
pexsi_delta_e                  20 #An upper bound for the spectral radius of S^{-1} H
pexsi_mu_lower                 -10 #Initial guess of lower bound for mu
pexsi_mu_upper                 10 #Initial guess of upper bound for mu
pexsi_mu                       0 #Initial guess for mu (for the solver)
pexsi_mu_thr                   0.05 #Stopping criterion in terms of the chemical potential for the inertia counting procedure
pexsi_mu_expand                0.3 #If the chemical potential is not in the initial interval, the interval is expanded by muInertiaExpansion
pexsi_mu_guard                 0.2 #Safe guard criterion in terms of the chemical potential to reinvoke the inertia counting procedure
pexsi_elec_thr                 0.001 #Stopping criterion of the PEXSI iteration in terms of the number of electrons compared to numElectronExact
pexsi_zero_thr                 1e-10 #if the absolute value of matrix element is less than ZERO_Limit, it will be considered as 0

#Parameters (20.Test)
out_alllog                     0 #output information for each processor, when parallel
nurse                          0 #for coders
t_in_h                         1 #calculate the kinetic energy or not
vl_in_h                        1 #calculate the local potential or not
vnl_in_h                       1 #calculate the nonlocal potential or not
vh_in_h                        1 #calculate the hartree potential or not
vion_in_h                      1 #calculate the local ionic potential or not
test_force                     0 #test the force
test_stress                    0 #test the stress
test_skip_ewald                0 #whether to skip ewald
ri_hartree_benchmark           none #whether to use the RI approximation for the Hartree term in LR-TDDFT for benchmark (with FHI-aims/ABACUS read-in style)
aims_nbasis                     #the number of basis functions for each atom type used in FHI-aims (for benchmark)
rdmft                          0 #whether to perform rdmft calculation, default is false
rdmft_power_alpha              0.656 #the alpha parameter of power-functional, g(occ_number) = occ_number^alpha used in exx-type functionals such as muller and power
