#EXAMPLE# Bethe Lattice 1-band Hubbard model, single-site DMFT, Bethe lattice ##################################### # To load Sigma(w) # put the correct path or '' to avoid loading # It must be a correct path from the directory # where the lisa script is executed ##################################### Reload_Sigma= False ##################################### # General parameters ##################################### # Maximum number of iterations of the main loop N_Iter_SelfCons_Max = 5 # Temperature Beta = 10.0 # 1/Temperature is_SUPRA = False ##################################### # Scheme ##################################### SCHEME = "DMFT_DOS" # The number of orbitals is given by the size of this array Half_Bandwidths=[sqrt(2)] AF_Computation = True ##################################### # Forcing terms ##################################### Forcing_SC = [0.0] Forcing_AF = [0.0] Forcing_F = [0.0] Forcing_DDW = [0.0] Number_Iterations_With_Forcing =1 ##################################### # Should we compute at fixed density ? ##################################### Fix_Density = True Chemical_potential_init = 0.0 + 1.5 # Not used if fixdoping=False but required for compilation doping = 0.1 Density_Required = 1 - doping Precision_Density = 0.01 Delta_Chemical_Potential = 0.1 # MUST be non 0 or it crashes the code !! Max_loop_mu = 10 Nbre_Iter_wait_adjust_dop = 0 ##################################### # Interactions # Must be a matrix ..... ##################################### import utils U_interact = 3.0 # The Coulomb interaction U_matrix_anti = utils.ScalarMatrix(1,U_interact) U_matrix_para = utils.ScalarMatrix(1,0.0) ##################################### # Choose the solver ##################################### Solver = "QMC_Hirsh_Fye" QMC_Number_Time_Slices = 16 #int(1.0*Beta* 2.0 ) QMC_N_cycles_MAX = 100000 #32*1000 QMC_Compute_Susceptibility = True ##################################### # Where do we put the data ? ##################################### Data_Directory = 'ess-Beta%3.2f-U%3.2f-dop%3.2f'%(Beta,U_interact,doping) Data_OverWrite = False