********------ 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': [95, '# 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': [30.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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76833117740405, 26.148093450005931, 26.571612684419389, 26.57465030538404, 26.141288080907881] 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/5417.1.wp02m/machines -port 20417 -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.04system 1:39.16elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1737minor)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/5417.1.wp02m/machines -port 20417 -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.0758495725 Found DMFT-mu= 11.0758351594 ... 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.0758351594 EF= 11.0758351594 Charge 0 = 36.0000000309 EF= 11.0757351594 Charge 0 = 35.9997180565 EF= 11.0759351594 Charge 0 = 36.0002820572 EF= 11.0758351392 Charge 0 = 35.999999974 EF= 11.0758351484 Charge 0 = 36.0 EF= 11.0758351484 Charge 0 = 36.0 New chemical potential found at -1.09482769517e-08 -> Chemical potential becomes 11.0758351484 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.76833117740405, 26.148093450005931, 26.571612684419389, 26.57465030538404, 26.141288080907881] 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/5417.1.wp02m/machines -port 20417 -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.06system 0:41.63elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+104outputs (0major+1738minor)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.04260202 -26.85688297 -26.99338334 -26.99520039 -26.85554516] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21637619085 /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.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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': [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': [95, '# 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': [30.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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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.00user 0.02system 1:36.97elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1704minor)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/5417.1.wp02m/machines -port 20417 -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.0757135156 Found DMFT-mu= 11.075775382 ... 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.075775382 EF= 11.075775382 Charge 0 = 35.9999734948 EF= 11.175775382 Charge 0 = 36.3216219853 (mu,dens)= 11.175775382 0.321621985275 (mu0,mu1), (dens0,dens1) 11.075775382 11.175775382 -2.65051642785e-05 0.321621985275 EF= 11.075775382 Charge 0 = 35.9999734948 EF= 11.175775382 Charge 0 = 36.3216219853 EF= 11.0757836225 Charge 0 = 35.9999966988 EF= 11.0757847948 Charge 0 = 36.0 EF= 11.0757847948 Charge 0 = 36.0 New chemical potential found at 9.41274291577e-06 -> Chemical potential becomes 11.0757847948 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.767515839597909, 26.147673830869952, 26.572430100089139, 26.574572016746849, 26.141068306373189] 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/5417.1.wp02m/machines -port 20417 -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.03system 0:37.08elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1753minor)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.04131718 -26.85599064 -26.99217224 -26.99398925 -26.85465324] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21543268395 /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.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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': [95, '# 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': [30.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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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.04system 1:38.82elapsed 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/5417.1.wp02m/machines -port 20417 -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.0738087187 Found DMFT-mu= 11.073836754 ... 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.073836754 EF= 11.073836754 Charge 0 = 35.9998525194 EF= 11.173836754 Charge 0 = 36.3212764384 (mu,dens)= 11.173836754 0.321276438443 (mu0,mu1), (dens0,dens1) 11.073836754 11.173836754 -0.000147480550822 0.321276438443 EF= 11.073836754 Charge 0 = 35.9998525194 EF= 11.173836754 Charge 0 = 36.3212764384 EF= 11.0738826375 Charge 0 = 35.99998162 EF= 11.0738891695 Charge 0 = 35.9999999998 EF= 11.0738891696 Charge 0 = 36.0 EF= 11.0738891696 Charge 0 = 36.0 New chemical potential found at 5.24156179438e-05 -> Chemical potential becomes 11.0738891696 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.766058352882361, 26.14570738278681, 26.57021963302148, 26.57331916720241, 26.13976916361413] 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/5417.1.wp02m/machines -port 20417 -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.04system 0:42.86elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1706minor)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.03940224 -26.85361144 -26.99014317 -26.99196011 -26.85227417] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21535755253 /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.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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 95 -t 800 case=FeSi, nom=95, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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': [110, '# 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': [30.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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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 2:07.07elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+120outputs (0major+1753minor)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/5417.1.wp02m/machines -port 20417 -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.0735366049 Found DMFT-mu= 11.073586206 ... 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.073586206 EF= 11.073586206 Charge 0 = 35.9999909763 EF= 11.173586206 Charge 0 = 36.3227407005 (mu,dens)= 11.173586206 0.322740700512 (mu0,mu1), (dens0,dens1) 11.073586206 11.173586206 -9.02374393519e-06 0.322740700512 EF= 11.073586206 Charge 0 = 35.9999909763 EF= 11.173586206 Charge 0 = 36.3227407005 EF= 11.0735890019 Charge 0 = 35.9999988773 EF= 11.0735893992 Charge 0 = 36.0 EF= 11.0735893992 Charge 0 = 36.0 New chemical potential found at 3.19317979702e-06 -> Chemical potential becomes 11.0735893992 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765291310157849, 26.144792956871441, 26.56862212881871, 26.571569590340161, 26.13827217682168] 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/5417.1.wp02m/machines -port 20417 -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.02system 0:38.19elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1752minor)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.03895047 -26.85201441 -26.98941105 -26.99122793 -26.8506775 ] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.2157224072 /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.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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': [110, '# 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': [30.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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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.01user 0.02system 1:40.27elapsed 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/5417.1.wp02m/machines -port 20417 -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.0753053449 Found DMFT-mu= 11.0753703028 ... 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.0753703028 EF= 11.0753703028 Charge 0 = 36.0000560082 EF= 10.9753703028 Charge 0 = 35.708139893 (mu,dens)= 10.9753703028 -0.291860107037 (mu0,mu1), (dens0,dens1) 11.0753703028 10.9753703028 5.60082094623e-05 -0.291860107037 EF= 11.0753703028 Charge 0 = 36.0000560082 EF= 10.9753703028 Charge 0 = 35.708139893 EF= 11.0753511164 Charge 0 = 36.0000018278 EF= 11.0753504691 Charge 0 = 36.0 EF= 11.0753504691 Charge 0 = 36.0 EF= 11.0753504691 Charge 0 = 36.0 New chemical potential found at -1.98336805965e-05 -> Chemical potential becomes 11.0753504691 /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.764260587662061, 26.14474392808939, 26.569174638843482, 26.571457477291879, 26.137555559906509] 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/5417.1.wp02m/machines -port 20417 -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.00user 0.05system 0:43.41elapsed 0%CPU (0avgtext+0avgdata 0maxresident)k 0inputs+112outputs (0major+1721minor)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.03967934 -26.8533164 -26.99024041 -26.99205723 -26.85197972] Edc= [ 25.75 25.75 25.75 25.75 25.75] ncorr= 6.21490635644 /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.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] Edc= [25.75, 25.75, 25.75, 25.75, 25.75] /home/users/jtomczak/dmft_w2k_111/bin/ssplit.py -n 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] 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/5417.1.wp02m/machines -port 20417 -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 110 -t 800 case=FeSi, nom=110, ntail=800, insig=sig.inp, outsig=sig.inp[1-n] s_oo= [25.765227027191809, 26.144689244224459, 26.56912277785052, 26.571817631221911, 26.138490082426859] 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/5417.1.wp02m/machines -port 20417 -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