********------ 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 = 6 ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [300000000, '# 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.10000000000000001, '# 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': [300, '# Number of Matsubara frequency points sampled'], 'aom': [20, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [300000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [100.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': [3000, '# 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79952936436327, 26.182604927603808, 26.641883002626429, 26.64304948430868, 26.17406847336342] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.01user 0.00system 1:16.40elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+1085minor)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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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.0940811997 Found DMFT-mu= 11.0940649147 ... 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.0940649147 EF= 11.0940649147 Charge 0 = 35.9999837199 EF= 11.1940649147 Charge 0 = 36.3931471759 (mu,dens)= 11.1940649147 0.393147175931 (mu0,mu1), (dens0,dens1) 11.0940649147 11.1940649147 -1.6280149147e-05 0.393147175931 EF= 11.0940649147 Charge 0 = 35.9999837199 EF= 11.1940649147 Charge 0 = 36.3931471759 EF= 11.0940690555 Charge 0 = 35.9999928656 EF= 11.0940722856 Charge 0 = 36.0000000006 EF= 11.0940722853 Charge 0 = 36.0 EF= 11.0940722853 Charge 0 = 36.0 New chemical potential found at 7.37064014267e-06 -> Chemical potential becomes 11.0940722853 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79952936436327, 26.182604927603808, 26.641883002626429, 26.64304948430868, 26.17406847336342] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.01system 1:22.73elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1085minor)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.08822416 -26.88992262 -27.03723249 -27.03905094 -26.88857666] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.2265686994 /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.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 4 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 5 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': [300000000, '# 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.10000000000000001, '# 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': [300, '# Number of Matsubara frequency points sampled'], 'aom': [50, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [300000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [100.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': [3000, '# 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.00system 1:17.25elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+1085minor)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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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.0959859966 Found DMFT-mu= 11.0959989416 ... 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.0959989416 EF= 11.0959989416 Charge 0 = 36.0000257244 EF= 10.9959989416 Charge 0 = 35.7955864949 (mu,dens)= 10.9959989416 -0.204413505097 (mu0,mu1), (dens0,dens1) 11.0959989416 10.9959989416 2.57243825317e-05 -0.204413505097 EF= 11.0959989416 Charge 0 = 36.0000257244 EF= 10.9959989416 Charge 0 = 35.7955864949 EF= 11.0959863587 Charge 0 = 35.9999958482 EF= 11.0959881073 Charge 0 = 35.9999999994 EF= 11.0959881076 Charge 0 = 36.0 EF= 11.0959881076 Charge 0 = 36.0 New chemical potential found at -1.08340421399e-05 -> Chemical potential becomes 11.0959881076 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793710010322592, 26.183292073377501, 26.641862206509121, 26.645183981298221, 26.174133946931281] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.00system 1:23.30elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1085minor)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.08814418 -26.89310929 -27.03805933 -27.03987782 -26.89176321] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22664055428 /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.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 4 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 5 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': [300000000, '# 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.10000000000000001, '# 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': [300, '# Number of Matsubara frequency points sampled'], 'aom': [50, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [300000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [100.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': [3000, '# 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.00system 1:16.95elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1086minor)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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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.0944893705 Found DMFT-mu= 11.0944777879 ... 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.0944777879 EF= 11.0944777879 Charge 0 = 35.9998327194 EF= 11.1944777879 Charge 0 = 36.3957819852 (mu,dens)= 11.1944777879 0.395781985159 (mu0,mu1), (dens0,dens1) 11.0944777879 11.1944777879 -0.000167280594411 0.395781985159 EF= 11.0944777879 Charge 0 = 35.9998327194 EF= 11.1944777879 Charge 0 = 36.3957819852 EF= 11.0945200359 Charge 0 = 35.9999322336 EF= 11.0945488007 Charge 0 = 36.0000000526 EF= 11.0945487784 Charge 0 = 36.0 EF= 11.0945487784 Charge 0 = 36.0 EF= 11.0945487784 Charge 0 = 36.0 New chemical potential found at 7.09904952529e-05 -> Chemical potential becomes 11.0945487784 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795077189179359, 26.182055526844849, 26.641103300480822, 26.643389515349689, 26.174619343766491] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.01system 1:24.18elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1085minor)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.08657605 -26.89121821 -27.03633627 -27.03815469 -26.88987229] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22643571417 /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.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 4 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 5 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': [300000000, '# 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.10000000000000001, '# 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': [300, '# Number of Matsubara frequency points sampled'], 'aom': [50, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [300000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [100.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': [3000, '# 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.01system 1:10.68elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+1085minor)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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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.095305712 Found DMFT-mu= 11.0953025939 ... 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.0953025939 EF= 11.0953025939 Charge 0 = 35.9999917634 EF= 11.1953025939 Charge 0 = 36.398763795 (mu,dens)= 11.1953025939 0.398763795022 (mu0,mu1), (dens0,dens1) 11.0953025939 11.1953025939 -8.23662331584e-06 0.398763795022 EF= 11.0953025939 Charge 0 = 35.9999917634 EF= 11.1953025939 Charge 0 = 36.398763795 EF= 11.0953046594 Charge 0 = 35.9999966816 EF= 11.0953060531 Charge 0 = 36.0000000001 EF= 11.095306053 Charge 0 = 36.0 EF= 11.095306053 Charge 0 = 36.0 New chemical potential found at 3.45908160426e-06 -> Chemical potential becomes 11.095306053 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.795335473213171, 26.18255689709072, 26.640099784946731, 26.643481202644409, 26.173761195363099] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.01system 1:14.48elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1085minor)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.08733092 -26.89193505 -27.03707304 -27.03889148 -26.89058909] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22615083335 /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.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 4 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 300 -t 800 case=FeSi, nom=300, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 5 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': [300000000, '# 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.10000000000000001, '# 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': [325, '# Number of Matsubara frequency points sampled'], 'aom': [50, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [300000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [100.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': [3000, '# 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 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.00system 1:19.14elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+1085minor)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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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.0939451428 Found DMFT-mu= 11.0939001069 ... 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.0939001069 EF= 11.0939001069 Charge 0 = 36.0000129364 EF= 10.9939001069 Charge 0 = 35.7894999052 (mu,dens)= 10.9939001069 -0.210500094821 (mu0,mu1), (dens0,dens1) 11.0939001069 10.9939001069 1.29364163257e-05 -0.210500094821 EF= 11.0939001069 Charge 0 = 36.0000129364 EF= 10.9939001069 Charge 0 = 35.7894999052 EF= 11.0938939617 Charge 0 = 35.9999990786 EF= 11.0938943703 Charge 0 = 35.9999999999 EF= 11.0938943704 Charge 0 = 36.0 EF= 11.0938943704 Charge 0 = 36.0 New chemical potential found at -5.73656291003e-06 -> Chemical potential becomes 11.0938943704 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.79414878692468, 26.18133209829271, 26.640404096391439, 26.643736145304679, 26.173787464240931] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 0.00user 0.00system 1:16.54elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1085minor)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.08586965 -26.89002802 -27.03550638 -27.03732475 -26.88868224] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.22533368579 /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.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 4 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 325 -t 800 case=FeSi, nom=325, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.793195458533301, 26.180503368117009, 26.640130575810911, 26.642181988287408, 26.173165895063381] 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 32 -machinefile /tmp/5802.1.wp10/machines -port 20802 -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 5 4 done --------- lapw0: /home/users/jtomczak/wien2k/x -f FeSi lapw0 lapw1 -c: /home/users/jtomczak/wien2k/x -f FeSi lapw1 -c