********------ 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': [200000000, '# 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': [45, '# Number of Matsubara frequency points sampled'], 'aom': [5, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [200000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [15.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': [2000, '# 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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780002958000171, 26.14221917925985, 26.517839851132958, 26.520344817223361, 26.135413335772078] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.02user 0.02system 1:10.06elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 2656inputs+128outputs (9major+924minor)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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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.0791149386 Found DMFT-mu= 11.0791263066 ... 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.0791263066 EF= 11.0791263066 Charge 0 = 35.9999999625 EF= 11.0790263066 Charge 0 = 35.9996874304 EF= 11.0792263066 Charge 0 = 36.0003124973 EF= 11.0791263181 Charge 0 = 35.9999999986 EF= 11.0791263186 Charge 0 = 36.0 EF= 11.0791263186 Charge 0 = 36.0 New chemical potential found at 1.2012979056e-08 -> Chemical potential becomes 11.0791263186 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780002958000171, 26.14221917925985, 26.517839851132958, 26.520344817223361, 26.135413335772078] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.00system 1:29.49elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+136outputs (0major+932minor)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.02669965 -26.85663684 -26.98009587 -26.98199827 -26.85532365] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21108646795 /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.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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': [200000000, '# 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': [45, '# Number of Matsubara frequency points sampled'], 'aom': [5, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [200000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [15.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': [2000, '# 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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.02system 1:09.86elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+136outputs (0major+932minor)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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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.0799312801 Found DMFT-mu= 11.0799143715 ... 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.0799143715 EF= 11.0799143715 Charge 0 = 36.0000715365 EF= 10.9799143715 Charge 0 = 35.6823352421 (mu,dens)= 10.9799143715 -0.317664757884 (mu0,mu1), (dens0,dens1) 11.0799143715 10.9799143715 7.15365251978e-05 -0.317664757884 EF= 11.0799143715 Charge 0 = 36.0000715365 EF= 10.9799143715 Charge 0 = 35.6823352421 EF= 11.0798918571 Charge 0 = 36.0000012666 EF= 11.0798914513 Charge 0 = 36.0 EF= 11.0798914513 Charge 0 = 36.0 EF= 11.0798914513 Charge 0 = 36.0 New chemical potential found at -2.29202426265e-05 -> Chemical potential becomes 11.0798914513 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780712407774189, 26.142549876232199, 26.518849377141809, 26.522040153706619, 26.137524341087889] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.03system 1:29.79elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+144outputs (0major+933minor)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.02715901 -26.85844115 -26.98087268 -26.9827752 -26.85712761] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21158111541 /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.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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': [200000000, '# 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': [45, '# Number of Matsubara frequency points sampled'], 'aom': [5, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [200000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [15.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': [2000, '# 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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.00user 0.01system 1:08.75elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+128outputs (0major+933minor)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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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.0806115647 Found DMFT-mu= 11.0805589242 ... 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.0805589242 EF= 11.0805589242 Charge 0 = 36.0000145335 EF= 10.9805589242 Charge 0 = 35.6825283596 (mu,dens)= 10.9805589242 -0.317471640372 (mu0,mu1), (dens0,dens1) 11.0805589242 10.9805589242 1.45335357189e-05 -0.317471640372 EF= 11.0805589242 Charge 0 = 36.0000145335 EF= 10.9805589242 Charge 0 = 35.6825283596 EF= 11.0805543465 Charge 0 = 36.0000002559 EF= 11.0805542644 Charge 0 = 36.0 EF= 11.0805542644 Charge 0 = 36.0 New chemical potential found at -4.65973899644e-06 -> Chemical potential becomes 11.0805542644 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.781151453551871, 26.142503009421269, 26.518903017035239, 26.521708234956591, 26.1369735326513] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.02system 1:29.50elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+144outputs (0major+933minor)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.02779568 -26.85904775 -26.98149745 -26.98339975 -26.85773471] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21157490251 /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.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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': [200000000, '# 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': [45, '# Number of Matsubara frequency points sampled'], 'aom': [5, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [200000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [15.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': [2000, '# 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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.02system 1:08.68elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+136outputs (0major+933minor)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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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.0791149386 Found DMFT-mu= 11.0791379781 ... 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.0791379781 EF= 11.0791379781 Charge 0 = 35.9999120575 EF= 11.1791379781 Charge 0 = 36.3186056722 (mu,dens)= 11.1791379781 0.31860567219 (mu0,mu1), (dens0,dens1) 11.0791379781 11.1791379781 -8.79425421729e-05 0.31860567219 EF= 11.0791379781 Charge 0 = 35.9999120575 EF= 11.1791379781 Charge 0 = 36.3186056722 EF= 11.0791655728 Charge 0 = 35.9999982098 EF= 11.0791661462 Charge 0 = 36.0 EF= 11.0791661462 Charge 0 = 36.0 EF= 11.0791661462 Charge 0 = 36.0 New chemical potential found at 2.81680975842e-05 -> Chemical potential becomes 11.0791661462 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.783465366729128, 26.14487872975926, 26.520275019642291, 26.52230186867375, 26.13821588968521] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.02system 1:31.32elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+144outputs (0major+932minor)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.02671608 -26.85835696 -26.98055941 -26.98245893 -26.85704648] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21120621549 /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.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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': [200000000, '# 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': [45, '# Number of Matsubara frequency points sampled'], 'aom': [5, '# Number of frequency points used to determin the value of sigma at nom'], 'J': [0.69999999999999996, '# Hunds'], 'M': [200000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [15.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': [2000, '# 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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.01system 1:09.13elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+128outputs (0major+933minor)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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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.0776183125 Found DMFT-mu= 11.0776636927 ... 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.0776636927 EF= 11.0776636927 Charge 0 = 35.9998300773 EF= 11.1776636927 Charge 0 = 36.3181866585 (mu,dens)= 11.1776636927 0.31818665848 (mu0,mu1), (dens0,dens1) 11.0776636927 11.1776636927 -0.000169922729015 0.31818665848 EF= 11.0776636927 Charge 0 = 35.9998300773 EF= 11.1776636927 Charge 0 = 36.3181866585 EF= 11.0777170676 Charge 0 = 35.9999966379 EF= 11.077718145 Charge 0 = 36.0 EF= 11.077718145 Charge 0 = 36.0 EF= 11.077718145 Charge 0 = 36.0 New chemical potential found at 5.4452365795e-05 -> Chemical potential becomes 11.077718145 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780136559852949, 26.141685625479859, 26.519180959476291, 26.522096718993499, 26.137418750671159] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 0.01user 0.02system 1:32.38elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+136outputs (0major+934minor)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.02619551 -26.85729082 -26.97982331 -26.98172272 -26.8559794 ] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21119360883 /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.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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 45 -t 800 case=FeSi, nom=45, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.780513820228531, 26.143334712542721, 26.519089332504169, 26.52128992580203, 26.13622663470515] 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 20 -machinefile /tmp/7597.1.wp10/machines -port 22597 -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