solver.weiss.dat

Introduction

This file is used to store the imaginary-time Weiss's function $G_0(\tau)$. It will be output by the quantum impurity solvers when they are shut down.

Format

The solver.weiss.dat file contains norbs block. Each block is appended by two blank lines. The format of each block is as follows:

column 1: orbital index $i$, integer

column 2: imaginary-time index $j$, integer

column 3: imaginary-time point, $\tau$, double precision

column 4: imaginary-time Weiss's function, $G_0(\tau)$, double precision

column 5: error bar, $\delta G_0(\tau)$, double precision

Note

In the solver.weiss.dat file, we adopt the following orbital sequence: $1\uparrow$, $2\uparrow$, $3\uparrow$, $\cdots$, $1\downarrow$, $2\downarrow$, $3\downarrow$, $\cdots$ In other words, the spin up part is always before the spin down part.

Code

The corresponding Fortran code block for the writing of solver.weiss.dat file is as follows:

! open data file: solver.weiss.dat
open(mytmp, file='solver.weiss.dat', form='formatted', status='unknown')

! write it
do i=1,norbs
    do j=1,ntime
        write(mytmp,'(2i6,3f12.6)') i, j, tmesh(j), wtau(j,i,i), zero
    enddo ! over j={1,ntime} loop
    write(mytmp,*) ! write empty lines
    write(mytmp,*)
enddo ! over i={1,norbs} loop

! close data file
close(mytmp)

Note

The column for the error bar is always zero in this file.

In the iQIST software package, we provide a Python module to read the output files of quantum impurity solvers. You can use it to read the solver.weiss.dat file. See [script/u_reader.py] for more details.