[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".
ops=q_d q_d^2 pair_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

options=SCHUR
