********------ parameter change ------*********** ntail = 800 ********------ parameter change ------*********** finish = 5 ********------ parameter change ------*********** recompute_mu = True ********------ parameter change ------*********** solver = CTQMC ********------ parameter change ------*********** wbroad = 0.01 ********------ parameter change ------*********** max_lda_iterations = 4 ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.02, '# Maximum derivative mismatch accepted for tail concatenation'], 'nOm': [100, 'num bos freq for vertex'], 'CleanUpdate': [100000, '# How often to make clean update'], 'minM': [1e-10, '# The smallest allowed value for the atomic trace'], 'SampleSusc': [0, '#bla'], 'nc': [[4, 5, 6, 7, 8], '# Impurity occupancies'], 'minD': [1e-10, '# The smallest allowed value for the determinant'], 'ChooseRandom': [0.80000000000000004, '# How often to choose kink in random bath and how often in those with less kinks'], 'Ncorrect': [-1, '# Which baths should not be corrected'], 'nomv': [100, 'num ferm freq for vertex'], 'nf0': [6.0, '# Double counting parameter'], 'nom': [320, '# Number of Matsubara frequency points sampled'], 'aom': [40, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [80.0, '# Inverse temperature'], 'U': [5.0, '# Coulomb repulsion (F0)'], 'GlobalFlip': [10000, '# How often to try a global flip'], 'Naver': [1000000000, '# How often to print out debug info'], 'Ud': [5.0, '# Coulomb repulsion for DC'], 'Nmax': [2500, '# Maximum perturbation order allowed'], 'exe': ['ctqmc', '# Name of the executable'], 'warmup': [100000, '# Warmup number of QMC steps'], 'SampleGtau': [1000, '# How often to update G(tau)']} ********------ parameter change ------*********** kbroad = 0.01 ********------ parameter change ------*********** fh_info = ***** Reading Structure file ***** Number of sorts, nat= 2 unit vectors (a,b,c) are 8.49054 8.49054 8.49054 angles betwen unit vectors are [ 90. 90. 90.] ---- over all types and atoms ---- ---------- atom type 0 ------------ iatnr= -1 Position inside the unit cell, pos= [0.1358, 0.1358, 0.1358] Number of atoms of this type, mult= 4 Symmetry of the atom, isplit= 4 *** Equivalent atoms *** iatnr= -1 pos= [0.36420000000000002, 0.86419999999999997, 0.63580000000000003] iatnr= -1 pos= [0.63580000000000003, 0.36420000000000002, 0.86419999999999997] iatnr= -1 pos= [0.86419999999999997, 0.63580000000000003, 0.36420000000000002] --------- Information about the current type ------------- aname= Fe1 Znuc= 26.0 Rmt= 2.26 RotLoc= [[0.70710680000000004, 0.40824830000000001, 0.57735029999999998], [-0.70710680000000004, 0.40824830000000001, 0.57735029999999998], [0.0, -0.81649660000000002, 0.57735029999999998]] ----------------------------------------------------------- ---------- atom type 1 ------------ iatnr= -2 Position inside the unit cell, pos= [0.84399999999999997, 0.84399999999999997, 0.84399999999999997] Number of atoms of this type, mult= 4 Symmetry of the atom, isplit= 4 *** Equivalent atoms *** iatnr= -2 pos= [0.65600000000000003, 0.156, 0.34399999999999997] iatnr= -2 pos= [0.34399999999999997, 0.65600000000000003, 0.156] iatnr= -2 pos= [0.156, 0.34399999999999997, 0.65600000000000003] --------- Information about the current type ------------- aname= Si2 Znuc= 14.0 Rmt= 2.0 RotLoc= [[0.70710680000000004, 0.40824830000000001, 0.57735029999999998], [-0.70710680000000004, 0.40824830000000001, 0.57735029999999998], [0.0, -0.81649660000000002, 0.57735029999999998]] ----------------------------------------------------------- Z= [26.0, 26.0, 26.0, 26.0, 14.0, 14.0, 14.0, 14.0] Znucs= {1: 26, 2: 26, 3: 26, 4: 26} iSiginds= {0: [[1, 0, 0, 0, 0], [0, 2, 0, 0, 0], [0, 0, 3, 0, 0], [0, 0, 0, 4, 0], [0, 0, 0, 0, 5]]} Impurity-lattice connection: icols_ind= {0: [1, 2, 3, 4]} --------- Preparing GF calculation --------- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py -d dmft0 >> dmft0_info.out 2>&1 --------- Preparing GF calculation --------- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py -d dmft1 >> dmft1_info.out 2>&1 --------- Preparing CHR calculation --------- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py -d -m 1 -x 1.0 -w 1.0 dmft2 >> ksum_info.out 2>&1 lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.7883401158095, 26.175248648238298, 26.63267300853521, 26.63573423927178, 26.166608033872791] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] --------- Computing eigenvalues --------- Running ---- ksum ---- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.02system 1:02.73elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2532minor)pagefaults 0swaps --------- Determine the chemical potential --------- Running ---- x_dmft.py ---- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py mu DmftEnvironment: mpi_prefix.dat exists -- running in parallel mode. /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 case= FeSi root= /home/users/jtomczak/dmft_w2k_111/bin cixs= [1, 2, 3, 4] NOE= 72.0 opening FeSi.scf2 Ef= 11.0913600613 Found DMFT-mu= 11.0914004984 ... Found matsubara and processing imaginary axis ... 1 Large matsubara mesh of size 1999 was created Looking for the chemical potential. Starting from old_mu= 11.0914004984 EF= 11.0914004984 Charge 0 = 36.0000000299 EF= 11.0913004984 Charge 0 = 35.9997687184 EF= 11.0915004984 Charge 0 = 36.0002318016 EF= 11.0914003861 Charge 0 = 35.9999997699 EF= 11.0914004855 Charge 0 = 36.0 EF= 11.0914004855 Charge 0 = 36.0 New chemical potential found at -1.291665086e-08 -> Chemical potential becomes 11.0914004855 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.7883401158095, 26.175248648238298, 26.63267300853521, 26.63573423927178, 26.166608033872791] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] Running ---- ksum ----- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.03system 0:25.47elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2531minor)pagefaults 0swaps /home/users/jtomczak/dmft_w2k_111/bin/sjoin.py -m 1.0 case=FeSi, indlt=FeSi.dlt?, outdlt=imp.?/Delta.imp, inE=FeSi.Eimp?, outE=imp.?/Eimp.inp, m_extn= icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} allcols= [1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5] noccur= [4 4 4 4 4] icix=1 reading from: FeSi.dlt1 cols= [1, 2, 3, 4, 5] icix=2 reading from: FeSi.dlt2 cols= [1, 2, 3, 4, 5] icix=3 reading from: FeSi.dlt3 cols= [1, 2, 3, 4, 5] icix=4 reading from: FeSi.dlt4 cols= [1, 2, 3, 4, 5] icix=0 saving into: imp.0/Delta.imp cols= [1, 2, 3, 4, 5] UpdateAtom= False Running ----- impurity solver ----- Taking care of high frequency Eimp= [-27.08142459 -26.88523195 -27.03085919 -27.03267736 -26.88388708] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22475392871 /home/users/jtomczak/dmft_w2k_111/bin/sgather.py case=FeSi, isig=imp.?/sig.out, osig=sig.inp icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} icix=0 reading from: imp.0/sig.out cols= [1, 2, 3, 4, 5] ------- DMFT part of the step # 0 done --------- Averaging over the following self-energies ['sig.inp.1.1'] saverage.py -o sig.inp sig.inp.1.1 files to average over: ['sig.inp.1.1'] output: sig.inp s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 0 0 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 1 0 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 2 0 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 3 0 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.02, '# Maximum derivative mismatch accepted for tail concatenation'], 'nOm': [100, 'num bos freq for vertex'], 'CleanUpdate': [100000, '# How often to make clean update'], 'minM': [1e-10, '# The smallest allowed value for the atomic trace'], 'SampleSusc': [0, '#bla'], 'nc': [[4, 5, 6, 7, 8], '# Impurity occupancies'], 'minD': [1e-10, '# The smallest allowed value for the determinant'], 'ChooseRandom': [0.80000000000000004, '# How often to choose kink in random bath and how often in those with less kinks'], 'Ncorrect': [-1, '# Which baths should not be corrected'], 'nomv': [100, 'num ferm freq for vertex'], 'nf0': [6.0, '# Double counting parameter'], 'nom': [320, '# Number of Matsubara frequency points sampled'], 'aom': [40, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [80.0, '# Inverse temperature'], 'U': [5.0, '# Coulomb repulsion (F0)'], 'GlobalFlip': [10000, '# How often to try a global flip'], 'Naver': [1000000000, '# How often to print out debug info'], 'Ud': [5.0, '# Coulomb repulsion for DC'], 'Nmax': [2500, '# Maximum perturbation order allowed'], 'exe': ['ctqmc', '# Name of the executable'], 'warmup': [100000, '# Warmup number of QMC steps'], 'SampleGtau': [1000, '# How often to update G(tau)']} /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] --------- Computing eigenvalues --------- Running ---- ksum ---- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.02user 0.03system 1:01.00elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2532minor)pagefaults 0swaps --------- Determine the chemical potential --------- Running ---- x_dmft.py ---- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py mu DmftEnvironment: mpi_prefix.dat exists -- running in parallel mode. /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 case= FeSi root= /home/users/jtomczak/dmft_w2k_111/bin cixs= [1, 2, 3, 4] NOE= 72.0 opening FeSi.scf2 Ef= 11.0910879475 Found DMFT-mu= 11.0910377882 ... Found matsubara and processing imaginary axis ... 1 Large matsubara mesh of size 1999 was created Looking for the chemical potential. Starting from old_mu= 11.0910377882 EF= 11.0910377882 Charge 0 = 35.9998360728 EF= 11.1910377882 Charge 0 = 36.3792289486 (mu,dens)= 11.1910377882 0.379228948615 (mu0,mu1), (dens0,dens1) 11.0910377882 11.1910377882 -0.00016392717275 0.379228948615 EF= 11.0910377882 Charge 0 = 35.9998360728 EF= 11.1910377882 Charge 0 = 36.3792289486 EF= 11.091080996 Charge 0 = 35.9999354015 EF= 11.0911090914 Charge 0 = 36.0000000339 EF= 11.0911090767 Charge 0 = 36.0 EF= 11.0911090767 Charge 0 = 36.0 EF= 11.0911090767 Charge 0 = 36.0 New chemical potential found at 7.12884765832e-05 -> Chemical potential becomes 11.0911090767 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788868207723208, 26.176618218393791, 26.63330434516277, 26.635854474293929, 26.167909358721779] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] Running ---- ksum ----- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.02system 0:25.86elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2532minor)pagefaults 0swaps /home/users/jtomczak/dmft_w2k_111/bin/sjoin.py -m 1.0 case=FeSi, indlt=FeSi.dlt?, outdlt=imp.?/Delta.imp, inE=FeSi.Eimp?, outE=imp.?/Eimp.inp, m_extn= icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} allcols= [1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5] noccur= [4 4 4 4 4] icix=1 reading from: FeSi.dlt1 cols= [1, 2, 3, 4, 5] icix=2 reading from: FeSi.dlt2 cols= [1, 2, 3, 4, 5] icix=3 reading from: FeSi.dlt3 cols= [1, 2, 3, 4, 5] icix=4 reading from: FeSi.dlt4 cols= [1, 2, 3, 4, 5] icix=0 saving into: imp.0/Delta.imp cols= [1, 2, 3, 4, 5] UpdateAtom= False Running ----- impurity solver ----- Taking care of high frequency Eimp= [-27.08092475 -26.88599668 -27.03066524 -27.03248351 -26.8846516 ] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22422717991 /home/users/jtomczak/dmft_w2k_111/bin/sgather.py case=FeSi, isig=imp.?/sig.out, osig=sig.inp icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} icix=0 reading from: imp.0/sig.out cols= [1, 2, 3, 4, 5] ------- DMFT part of the step # 0 done --------- Averaging over the following self-energies ['sig.inp.2.1'] saverage.py -o sig.inp sig.inp.2.1 files to average over: ['sig.inp.2.1'] output: sig.inp s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 0 1 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 1 1 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 2 1 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 3 1 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.02, '# Maximum derivative mismatch accepted for tail concatenation'], 'nOm': [100, 'num bos freq for vertex'], 'CleanUpdate': [100000, '# How often to make clean update'], 'minM': [1e-10, '# The smallest allowed value for the atomic trace'], 'SampleSusc': [0, '#bla'], 'nc': [[4, 5, 6, 7, 8], '# Impurity occupancies'], 'minD': [1e-10, '# The smallest allowed value for the determinant'], 'ChooseRandom': [0.80000000000000004, '# How often to choose kink in random bath and how often in those with less kinks'], 'Ncorrect': [-1, '# Which baths should not be corrected'], 'nomv': [100, 'num ferm freq for vertex'], 'nf0': [6.0, '# Double counting parameter'], 'nom': [320, '# Number of Matsubara frequency points sampled'], 'aom': [40, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [80.0, '# Inverse temperature'], 'U': [5.0, '# Coulomb repulsion (F0)'], 'GlobalFlip': [10000, '# How often to try a global flip'], 'Naver': [1000000000, '# How often to print out debug info'], 'Ud': [5.0, '# Coulomb repulsion for DC'], 'Nmax': [2500, '# Maximum perturbation order allowed'], 'exe': ['ctqmc', '# Name of the executable'], 'warmup': [100000, '# Warmup number of QMC steps'], 'SampleGtau': [1000, '# How often to update G(tau)']} /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] --------- Computing eigenvalues --------- Running ---- ksum ---- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.02system 1:02.25elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+96outputs (0major+2531minor)pagefaults 0swaps --------- Determine the chemical potential --------- Running ---- x_dmft.py ---- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py mu DmftEnvironment: mpi_prefix.dat exists -- running in parallel mode. /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 case= FeSi root= /home/users/jtomczak/dmft_w2k_111/bin cixs= [1, 2, 3, 4] NOE= 72.0 opening FeSi.scf2 Ef= 11.0908158337 Found DMFT-mu= 11.0908523145 ... Found matsubara and processing imaginary axis ... 1 Large matsubara mesh of size 1999 was created Looking for the chemical potential. Starting from old_mu= 11.0908523145 EF= 11.0908523145 Charge 0 = 35.9999735405 EF= 11.1908523145 Charge 0 = 36.3801675681 (mu,dens)= 11.1908523145 0.380167568058 (mu0,mu1), (dens0,dens1) 11.0908523145 11.1908523145 -2.6459515027e-05 0.380167568058 EF= 11.0908523145 Charge 0 = 35.9999735405 EF= 11.1908523145 Charge 0 = 36.3801675681 EF= 11.090859274 Charge 0 = 35.9999895215 EF= 11.0908638371 Charge 0 = 36.0000000009 EF= 11.0908638367 Charge 0 = 36.0 EF= 11.0908638367 Charge 0 = 36.0 New chemical potential found at 1.15222137378e-05 -> Chemical potential becomes 11.0908638367 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.789089298596849, 26.175651273342829, 26.63296176462984, 26.63493608014382, 26.16773735860475] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] Running ---- ksum ----- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.02user 0.03system 0:25.38elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2532minor)pagefaults 0swaps /home/users/jtomczak/dmft_w2k_111/bin/sjoin.py -m 1.0 case=FeSi, indlt=FeSi.dlt?, outdlt=imp.?/Delta.imp, inE=FeSi.Eimp?, outE=imp.?/Eimp.inp, m_extn= icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} allcols= [1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5] noccur= [4 4 4 4 4] icix=1 reading from: FeSi.dlt1 cols= [1, 2, 3, 4, 5] icix=2 reading from: FeSi.dlt2 cols= [1, 2, 3, 4, 5] icix=3 reading from: FeSi.dlt3 cols= [1, 2, 3, 4, 5] icix=4 reading from: FeSi.dlt4 cols= [1, 2, 3, 4, 5] icix=0 saving into: imp.0/Delta.imp cols= [1, 2, 3, 4, 5] UpdateAtom= False Running ----- impurity solver ----- Taking care of high frequency Eimp= [-27.08077493 -26.88542213 -27.03039207 -27.03221028 -26.88407714] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22488088397 /home/users/jtomczak/dmft_w2k_111/bin/sgather.py case=FeSi, isig=imp.?/sig.out, osig=sig.inp icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} icix=0 reading from: imp.0/sig.out cols= [1, 2, 3, 4, 5] ------- DMFT part of the step # 0 done --------- Averaging over the following self-energies ['sig.inp.3.1'] saverage.py -o sig.inp sig.inp.3.1 files to average over: ['sig.inp.3.1'] output: sig.inp s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 0 2 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 1 2 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 2 2 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 3 2 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.02, '# Maximum derivative mismatch accepted for tail concatenation'], 'nOm': [100, 'num bos freq for vertex'], 'CleanUpdate': [100000, '# How often to make clean update'], 'minM': [1e-10, '# The smallest allowed value for the atomic trace'], 'SampleSusc': [0, '#bla'], 'nc': [[4, 5, 6, 7, 8], '# Impurity occupancies'], 'minD': [1e-10, '# The smallest allowed value for the determinant'], 'ChooseRandom': [0.80000000000000004, '# How often to choose kink in random bath and how often in those with less kinks'], 'Ncorrect': [-1, '# Which baths should not be corrected'], 'nomv': [100, 'num ferm freq for vertex'], 'nf0': [6.0, '# Double counting parameter'], 'nom': [320, '# Number of Matsubara frequency points sampled'], 'aom': [40, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [80.0, '# Inverse temperature'], 'U': [5.0, '# Coulomb repulsion (F0)'], 'GlobalFlip': [10000, '# How often to try a global flip'], 'Naver': [1000000000, '# How often to print out debug info'], 'Ud': [5.0, '# Coulomb repulsion for DC'], 'Nmax': [2500, '# Maximum perturbation order allowed'], 'exe': ['ctqmc', '# Name of the executable'], 'warmup': [100000, '# Warmup number of QMC steps'], 'SampleGtau': [1000, '# How often to update G(tau)']} /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] --------- Computing eigenvalues --------- Running ---- ksum ---- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.02user 0.02system 1:01.96elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2531minor)pagefaults 0swaps --------- Determine the chemical potential --------- Running ---- x_dmft.py ---- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py mu DmftEnvironment: mpi_prefix.dat exists -- running in parallel mode. /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 case= FeSi root= /home/users/jtomczak/dmft_w2k_111/bin cixs= [1, 2, 3, 4] NOE= 72.0 opening FeSi.scf2 Ef= 11.0912240044 Found DMFT-mu= 11.0912651978 ... Found matsubara and processing imaginary axis ... 1 Large matsubara mesh of size 1999 was created Looking for the chemical potential. Starting from old_mu= 11.0912651978 EF= 11.0912651978 Charge 0 = 36.0000075894 EF= 10.9912651978 Charge 0 = 35.7748764464 (mu,dens)= 10.9912651978 -0.225123553588 (mu0,mu1), (dens0,dens1) 11.0912651978 10.9912651978 7.58944191404e-06 -0.225123553588 EF= 11.0912651978 Charge 0 = 36.0000075894 EF= 10.9912651978 Charge 0 = 35.7748764464 EF= 11.0912618267 Charge 0 = 35.9999998221 EF= 11.0912619039 Charge 0 = 36.0 EF= 11.0912619039 Charge 0 = 36.0 EF= 11.0912619039 Charge 0 = 36.0 New chemical potential found at -3.29391042797e-06 -> Chemical potential becomes 11.0912619039 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788059156609819, 26.175597028880709, 26.632881189693869, 26.63544179359635, 26.16762672831673] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] Running ---- ksum ----- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.02system 0:25.84elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2530minor)pagefaults 0swaps /home/users/jtomczak/dmft_w2k_111/bin/sjoin.py -m 1.0 case=FeSi, indlt=FeSi.dlt?, outdlt=imp.?/Delta.imp, inE=FeSi.Eimp?, outE=imp.?/Eimp.inp, m_extn= icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} allcols= [1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5] noccur= [4 4 4 4 4] icix=1 reading from: FeSi.dlt1 cols= [1, 2, 3, 4, 5] icix=2 reading from: FeSi.dlt2 cols= [1, 2, 3, 4, 5] icix=3 reading from: FeSi.dlt3 cols= [1, 2, 3, 4, 5] icix=4 reading from: FeSi.dlt4 cols= [1, 2, 3, 4, 5] icix=0 saving into: imp.0/Delta.imp cols= [1, 2, 3, 4, 5] UpdateAtom= False Running ----- impurity solver ----- Taking care of high frequency Eimp= [-27.08120888 -26.8857828 -27.03080016 -27.03261835 -26.8844378 ] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22505553294 /home/users/jtomczak/dmft_w2k_111/bin/sgather.py case=FeSi, isig=imp.?/sig.out, osig=sig.inp icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} icix=0 reading from: imp.0/sig.out cols= [1, 2, 3, 4, 5] ------- DMFT part of the step # 0 done --------- Averaging over the following self-energies ['sig.inp.4.1'] saverage.py -o sig.inp sig.inp.4.1 files to average over: ['sig.inp.4.1'] output: sig.inp s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 0 3 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 1 3 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 2 3 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 3 3 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.02, '# Maximum derivative mismatch accepted for tail concatenation'], 'nOm': [100, 'num bos freq for vertex'], 'CleanUpdate': [100000, '# How often to make clean update'], 'minM': [1e-10, '# The smallest allowed value for the atomic trace'], 'SampleSusc': [0, '#bla'], 'nc': [[4, 5, 6, 7, 8], '# Impurity occupancies'], 'minD': [1e-10, '# The smallest allowed value for the determinant'], 'ChooseRandom': [0.80000000000000004, '# How often to choose kink in random bath and how often in those with less kinks'], 'Ncorrect': [-1, '# Which baths should not be corrected'], 'nomv': [100, 'num ferm freq for vertex'], 'nf0': [6.0, '# Double counting parameter'], 'nom': [320, '# Number of Matsubara frequency points sampled'], 'aom': [40, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [80.0, '# Inverse temperature'], 'U': [5.0, '# Coulomb repulsion (F0)'], 'GlobalFlip': [10000, '# How often to try a global flip'], 'Naver': [1000000000, '# How often to print out debug info'], 'Ud': [5.0, '# Coulomb repulsion for DC'], 'Nmax': [2500, '# Maximum perturbation order allowed'], 'exe': ['ctqmc', '# Name of the executable'], 'warmup': [100000, '# Warmup number of QMC steps'], 'SampleGtau': [1000, '# How often to update G(tau)']} /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] --------- Computing eigenvalues --------- Running ---- ksum ---- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.03system 1:02.03elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+2532minor)pagefaults 0swaps --------- Determine the chemical potential --------- Running ---- x_dmft.py ---- /home/users/jtomczak/dmft_w2k_111/bin/x_dmft.py mu DmftEnvironment: mpi_prefix.dat exists -- running in parallel mode. /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 case= FeSi root= /home/users/jtomczak/dmft_w2k_111/bin cixs= [1, 2, 3, 4] NOE= 72.0 opening FeSi.scf2 Ef= 11.0914961183 Found DMFT-mu= 11.0914393404 ... Found matsubara and processing imaginary axis ... 1 Large matsubara mesh of size 1999 was created Looking for the chemical potential. Starting from old_mu= 11.0914393404 EF= 11.0914393404 Charge 0 = 35.9999887275 EF= 11.1914393404 Charge 0 = 36.3811681218 (mu,dens)= 11.1914393404 0.381168121751 (mu0,mu1), (dens0,dens1) 11.0914393404 11.1914393404 -1.12724507844e-05 0.381168121751 EF= 11.0914393404 Charge 0 = 35.9999887275 EF= 11.1914393404 Charge 0 = 36.3811681218 EF= 11.0914422977 Charge 0 = 35.9999955276 EF= 11.0914442427 Charge 0 = 36.0000000002 EF= 11.0914442426 Charge 0 = 36.0 EF= 11.0914442426 Charge 0 = 36.0 New chemical potential found at 4.90217991533e-06 -> Chemical potential becomes 11.0914442426 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.788078561633931, 26.175686391211279, 26.63337109289273, 26.635756541682809, 26.168126238981639] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] Running ---- ksum ----- /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.02user 0.02system 0:25.03elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+2533minor)pagefaults 0swaps /home/users/jtomczak/dmft_w2k_111/bin/sjoin.py -m 1.0 case=FeSi, indlt=FeSi.dlt?, outdlt=imp.?/Delta.imp, inE=FeSi.Eimp?, outE=imp.?/Eimp.inp, m_extn= icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} allcols= [1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5] noccur= [4 4 4 4 4] icix=1 reading from: FeSi.dlt1 cols= [1, 2, 3, 4, 5] icix=2 reading from: FeSi.dlt2 cols= [1, 2, 3, 4, 5] icix=3 reading from: FeSi.dlt3 cols= [1, 2, 3, 4, 5] icix=4 reading from: FeSi.dlt4 cols= [1, 2, 3, 4, 5] icix=0 saving into: imp.0/Delta.imp cols= [1, 2, 3, 4, 5] UpdateAtom= False Running ----- impurity solver ----- Taking care of high frequency Eimp= [-27.08138617 -26.88599894 -27.0309884 -27.03280659 -26.88465394] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22466648697 /home/users/jtomczak/dmft_w2k_111/bin/sgather.py case=FeSi, isig=imp.?/sig.out, osig=sig.inp icols= {0: [1, 2, 3, 4, 5]} cols= {1: [1, 2, 3, 4, 5], 2: [1, 2, 3, 4, 5], 3: [1, 2, 3, 4, 5], 4: [1, 2, 3, 4, 5]} icix=0 reading from: imp.0/sig.out cols= [1, 2, 3, 4, 5] ------- DMFT part of the step # 0 done --------- Averaging over the following self-energies ['sig.inp.5.1'] saverage.py -o sig.inp sig.inp.5.1 files to average over: ['sig.inp.5.1'] output: sig.inp s_oo= [25.78830986248644, 26.176661526309129, 26.633881717822039, 26.63615976747721, 26.167702628345332] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.78830986248644, 26.176661526309129, 26.633881717822039, 26.63615976747721, 26.167702628345332] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 0 4 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.78830986248644, 26.176661526309129, 26.633881717822039, 26.63615976747721, 26.167702628345332] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 1 4 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.78830986248644, 26.176661526309129, 26.633881717822039, 26.63615976747721, 26.167702628345332] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 2 4 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 320 -t 800 case=FeSi, nom=320, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.78830986248644, 26.176661526309129, 26.633881717822039, 26.63615976747721, 26.167702628345332] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] ..... Creating logarithmic mesh ..... Going over all correlated blocks icix= 1 cols= [1, 2, 3, 4, 5] icix= 2 cols= [1, 2, 3, 4, 5] icix= 3 cols= [1, 2, 3, 4, 5] icix= 4 cols= [1, 2, 3, 4, 5] dmft2: /usr/bin/time /opt/mpich2/intel/11.1/bin/mpiexec -n 144 -machinefile /tmp/8774.1.wp10/machines -port 23774 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft2 dmft2.def >> dmft2_info.out lcore: /home/users/jtomczak/wien2k/x -f FeSi lcore mixer: /home/users/jtomczak/wien2k/x -f FeSi mixer ------- LDA(charge) step 3 4 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c