# Input file for LDA+DMFT # This set of parameters can be edited during run. At each DMFT step, these parameters are dynamically updated. # The commands are in standard Python syntax solver="OCA" # Impurity solver recompute_mu=True # Recompute the chemical potential at each step DCs="fixn" # Double counting scheme. finish=1 # Number of charge self-consistent steps base="DMFTbase" # In which base is the greens function computed max_iterations=1000 # Number of charge self-consistent steps use_tetra=False # Use tetrahedron method to determin the chemical potential Niter=10 # Number of iterations in the inner DMFT loop runLDA=True # Run LDA at each step of the outside loop runIMP=True # Run impurity solver at each iteration in the inside loop DoubleCounting="impurity" # Is double counting computed from impurity density matrix or charge density matrix cut_ab=[-1000.0, 1000.0] # The eigenvalues which are outside this interval are ignored orthog=True # The DMFT base is made exactly orthonormal LowerBound=-30.0 # The lower limit when computing charge gbroad=0.01 # Broadening of the hybridization function gamma=[0.01, 0.01] # Broadening of the self-energy for non-correlated and correlated orbitals gammag=[0.01, 0.01] # Broadening of the self-energy for non-correlated and correlated orbitals sdmu=0.1 # Used to bracket the chemical potential mix_mu=1.0 # Mixing of the chemical potential max_metropolis_steps=50000 # Maximum number of metropolis steps in sorting eigenvalues. UpdateAtom=False # Does the impurity cix file needs to be recomputed at each step tetrahedronOm=-1 symmetrizeQ=False # Impurity problem number 0 iparams0={"exe" : ["nca" , "# Name of the executable"], "U" : [5.0 , "# Coulomb repulsion (F0)"], "T" : [0.00025 , "# Temperature"], "nf0" : [0.9 , "# Double counting parameter"], "nc" : [[0, 1, 2] , "# Impurity occupancies"], "Ncentral" : [[1] , "# Central occupancies for OCA diagrams evaluation"], "alpha" : [0.5 , "# Mixing for bath spectral function"], "max_steps" : [30 , "# Maximum number of impurity steps"], "max_diff" : [0.001 , "# Maximum difference between steps"], "followPeak" : [-1 , "# A mode to determin lambda0"], "Q" : [1.0 , "# A parameter to determin lambda0"], "StartLambda" : [-20.0 , "# Where to start looking for zero to determin lambda0"], "dLambda" : [0.1 , "# Step in searching for the lambda"], "EndLambda" : [1.0 , "# Where to stop searching for the lambda0"], "cutAc" : [[-10.0, 10.0] , "# Only window [La,Lb] of baths spectra is taken into account"], "Gh" : [0.5 , "# Parameter to improve the high frequency self-energy"], "epsilon" : [[[-10, -2.0], [2.0, 10.0]], "# Parameter to improve the high frequency self-energy"], "Th" : [0.5 , "# Parameter to improve the high frequency self-energy"], "lorentz" : [1 , "# Weather to subtract lorentz from diverging spectral functions and treat it analytically"], "SearchLorentz" : [2.5 , "# How far from zero to search fro Lorentz"], "LorentzMaxRatio" : [1.0 , "# How far from zero to search for Lorentz"], "FirstLorentz" : [0 , "# First pseudoparticle which could be augmented with lorentz"], "LastLorentz" : [10000 , "# Last pseudoparticle which could be augmented with lorentz"], "CmpDiff" : [-1 , "# When calculating Difference, only first CmpDiff particles should be taken into account (-1,all)"], "para" : [0 , "# Non paramagnetic mode"], }