********------ 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': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.01, '# 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': [145, '# 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': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [40.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': [1000, '# 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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76601711780151, 26.149617100312259, 26.588307549212441, 26.590736018984892, 26.143288269287481] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.05user 0.03system 1:49.91elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1722minor)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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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.0766659141 Found DMFT-mu= 11.0766098095 ... 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.0766098095 EF= 11.0766098095 Charge 0 = 36.0000000308 EF= 11.0765098095 Charge 0 = 35.9997365165 EF= 11.0767098095 Charge 0 = 36.0002636541 EF= 11.0766097771 Charge 0 = 35.9999999455 EF= 11.0766097978 Charge 0 = 36.0 EF= 11.0766097978 Charge 0 = 36.0 New chemical potential found at -1.16820881857e-08 -> Chemical potential becomes 11.0766097978 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76601711780151, 26.149617100312259, 26.588307549212441, 26.590736018984892, 26.143288269287481] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.03user 0.06system 0:44.62elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 2312inputs+104outputs (1major+1736minor)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.04834696 -26.85844388 -26.99855991 -27.000376 -26.85711201] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21770953906 /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.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 ********------ parameter change ------*********** iparams0 = {'Ntau': [1000, '# Number of imaginary time points (only for debugging)'], 'tsample': [10, '# How often to record measurements'], 'Ncout': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.01, '# 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': [145, '# 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': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [40.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': [1000, '# 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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.04user 0.09system 1:48.47elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1721minor)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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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.0773461986 Found DMFT-mu= 11.077290204 ... 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.077290204 EF= 11.077290204 Charge 0 = 36.0000252743 EF= 10.977290204 Charge 0 = 35.7239536517 (mu,dens)= 10.977290204 -0.276046348308 (mu0,mu1), (dens0,dens1) 11.077290204 10.977290204 2.5274312982e-05 -0.276046348308 EF= 11.077290204 Charge 0 = 36.0000252743 EF= 10.977290204 Charge 0 = 35.7239536517 EF= 11.077281049 Charge 0 = 36.0000010945 EF= 11.0772806347 Charge 0 = 36.0 EF= 11.0772806346 Charge 0 = 36.0 EF= 11.0772806346 Charge 0 = 36.0 New chemical potential found at -9.56937781638e-06 -> Chemical potential becomes 11.0772806346 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76631186438517, 26.15026914048692, 26.588214338154071, 26.590985289841552, 26.1426586426219] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.03user 0.08system 0:44.70elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 2512inputs+112outputs (1major+1736minor)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.04954981 -26.85866531 -26.99948518 -27.00130129 -26.85733348] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21762475222 /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.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.01, '# 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': [145, '# 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': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [40.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': [1000, '# 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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.02user 0.02system 1:47.39elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1721minor)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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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.0780264832 Found DMFT-mu= 11.0779769255 ... 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.0779769255 EF= 11.0779769255 Charge 0 = 35.9999983413 EF= 11.1779769255 Charge 0 = 36.331626557 (mu,dens)= 11.1779769255 0.331626556977 (mu0,mu1), (dens0,dens1) 11.0779769255 11.1779769255 -1.65867012214e-06 0.331626556977 EF= 11.0779769255 Charge 0 = 35.9999983413 EF= 11.1779769255 Charge 0 = 36.331626557 EF= 11.0779774257 Charge 0 = 35.9999996627 EF= 11.0779775533 Charge 0 = 36.0 EF= 11.0779775533 Charge 0 = 36.0 New chemical potential found at 6.27817147603e-07 -> Chemical potential becomes 11.0779775533 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76699796883403, 26.150975471000319, 26.589342786559079, 26.59170161932337, 26.143982169232888] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.04system 0:40.07elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 2464inputs+112outputs (1major+1751minor)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.05026437 -26.85946902 -27.00022471 -27.00204082 -26.85813718] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.2180258587 /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.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.01, '# 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': [145, '# 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': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [40.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': [1000, '# 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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.06user 0.06system 1:50.08elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1729minor)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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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.078434654 Found DMFT-mu= 11.0784475561 ... 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.0784475561 EF= 11.0784475561 Charge 0 = 36.0000382774 EF= 10.9784475561 Charge 0 = 35.7245992341 (mu,dens)= 10.9784475561 -0.27540076589 (mu0,mu1), (dens0,dens1) 11.0784475561 10.9784475561 3.82773957881e-05 -0.27540076589 EF= 11.0784475561 Charge 0 = 36.0000382774 EF= 10.9784475561 Charge 0 = 35.7245992341 EF= 11.0784336592 Charge 0 = 36.0000016167 EF= 11.0784330464 Charge 0 = 36.0000000001 EF= 11.0784330464 Charge 0 = 36.0 EF= 11.0784330464 Charge 0 = 36.0 New chemical potential found at -1.45097168449e-05 -> Chemical potential becomes 11.0784330464 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767493058784851, 26.151189445710578, 26.59001907691156, 26.592567904109931, 26.144920844793191] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.02user 0.05system 0:44.75elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 2472inputs+112outputs (1major+1720minor)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.05086213 -26.86063916 -27.00097992 -27.00279704 -26.8593001 ] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21830304674 /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.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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': [400000000, '# How often to print out info'], 'PChangeOrder': [0.90000000000000002, '# Ratio between trial steps: add-remove-a-kink / move-a-kink'], 'cx': [0.0, '# Spin-orbit'], 'sderiv': [0.01, '# 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': [145, '# 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': [400000000, '# Total number of Monte Carlo steps'], 'CoulombF': ["'Georges'", '# Full Coulomb repulsion'], 'beta': [40.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': [1000, '# 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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft0.def >> dmft0_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.03user 0.04system 1:48.81elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1730minor)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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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.0778904263 Found DMFT-mu= 11.0778477278 ... 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.0778477278 EF= 11.0778477278 Charge 0 = 36.0001890895 EF= 10.9778477278 Charge 0 = 35.7235344051 (mu,dens)= 10.9778477278 -0.276465594857 (mu0,mu1), (dens0,dens1) 11.0778477278 10.9778477278 0.000189089451545 -0.276465594857 EF= 11.0778477278 Charge 0 = 36.0001890895 EF= 10.9778477278 Charge 0 = 35.7235344051 EF= 11.0777793793 Charge 0 = 36.0000087668 EF= 11.0777760564 Charge 0 = 36.0000000015 EF= 11.0777760559 Charge 0 = 36.0 EF= 11.0777760559 Charge 0 = 36.0 New chemical potential found at -7.16719447027e-05 -> Chemical potential becomes 11.0777760559 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.768121864842151, 26.151146796844039, 26.589685148771149, 26.592817646521819, 26.145593552052489] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -env LD_LIBRARY_PATH /opt/intel/Compiler/11.1/072/mkl/lib/em64t -env OMP_NUM_THREADS 1 ./dmft dmft1.def >> dmft1_info.out DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END DMFT END 0.01user 0.04system 0:45.31elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 2472inputs+112outputs (1major+1737minor)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.05074665 -26.86024694 -27.00078665 -27.00260376 -26.85890772] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21765675301 /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.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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 145 -t 800 case=FeSi, nom=145, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76798702399693, 26.151333774546789, 26.590000720945959, 26.59227735564852, 26.145431872722462] 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 76 -machinefile /tmp/5342.1.wp02/machines -port 20342 -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