[param]
# Single impurity Anderson model
model=SIAM

# The only symmetry is the spin SU(2) symmetry
symtype=SPSU2

# BCS gap
bcsgap=1e-2

# SIAM model parameters
U=0.27
# We assume a flat band, i.e., Gamma(omega)=Gamma for all energies omega.
Gamma=0.01
# delta=epsilon+U/2
delta=0

# NRG discretization parameter. Lambda=4 is relatively large, but
# it turns out that the errors are still quite low with this choice,
# i.e., in the percent range.
Lambda=4

# Lowest temperature considered in the NRG iteration. 
# One must choose Tmin << bcsgap.
Tmin=1e-12

# NRG truncation parameters. Keep at least 100 states and at most
# 10000 states, with the goal of keeping states with the energies
# up to 9 omega_N, where omega_N is the characteristic energy scale
# at the N-th iteration. The values chosen here are sufficient for
# full convergence at Lambda=4.
keepmin=100
keepenergy=9.0
keep=10000

# NRG discretization scheme and the twist parameter z. The twist
# averaging is not necessary in this calculation.
discretization=Z
z=1

# Operators considered. Their expectation values will be stored
# in the output file "td".
# Here we add A_d = c^\dag_\sigma and Ac_d = c_-\sigma
ops=q_d q_d^2 pair_d A_d Ac_d

# We compute spectral functions <<c_\sigma;c_\sigma^\dag>> and
# <<c^\dag_-\sigma;c^\dag_\sigma>>.
specd=A_d-A_d Ac_d-A_d

# Saving of excited states in the output file "annotated.dat". 
# 100 states are stored after each iteration, energies are expressed 
# in absolute excitation energies, rather than rescaled in units
# of omega_N.
dumpannotated=100
dumpscaled=0
grouptol=1e-300

# Spectral function calculation settings:
# enable complete Fock space algorithm
cfs=true
# Number of digits of precision
prec_xy=12
# Broadening function kernel
smooth=new
# Logarithmic binning: number of bins per energy decade
bins=500
# Broadening parameter
alpha=0.8
# Save raw spectral data: necessary for determining the position of the sub-gap delta peaks
savebins=true
# Lorentzian broadening off since this is a T=0 calculation
omega0=1e-99

options=SCHUR
