Determinantal QMC (BSS-QMC)
Documentation of code by Fakher Assaad
Usage
- Preparation
- Create
paramHUB
- Create
seeds
- Create
confin
- Running
-
./Hubb.out
- Analysis
- Create
Files_tau
- Run analysis:
bash annal.c
Example input files with documentation
paramHUB
4.0 5.00 0.0 40 20 10 100
4 1
RHUB, BETA, CHEM, LTROT, NWRAP, NBIN, NSWEEP
L, Ltau
Name |
Description |
Value |
RHUB |
Hubbard interaction U |
4.0 |
BETA |
inverse temperature 1/T |
5.0 |
CHEM |
chemical potential (relative to U/2) |
0.0 |
LTROT |
number of Trotter time slices |
40 |
NWRAP |
number of sweeps (?) between recalculations of G |
20 |
NBIN |
number of result bins (with NSWEEP sweeps each) |
10 |
NSWEEP |
number of sweeps (each site twice) |
100 |
L |
linear dimension of lattice |
4 |
Ltau |
flag for computation of ? (1=yes, 0=no) |
1 |
seeds
754837
just one integer number
confin
Two possibilities: (i) empty Ising spin configuration:
0
(ii) configuration copied from
confout
of previous run:
131358357
1
-1
-1
1
-1
1
1
-1
-1
1
-1
-1
-1
1
1
1
1
-1
-1
1
-1
1
-1
1
-1
-1
-1
-1
-1
1
1
1
1
-1
1
-1
-1
1
-1
-1
-1
1
-1
-1
-1
1
Files_tau
3 3
spin 2
gr 1
charge 2
First line:
nf
nf1
(
nf
is dummy,
nf1
is number of subsequent lines)
subsequent lines:
file
Nflag
(value of
Nflag
is written to file
Flag
before
cov_tau.out
is called).
NFlag=2
indicates symmetry in tau <-> beta-tau: (from
cov_tau.out
):
If (NFlag.eq.2) Then
NTDM = LTROT/2 + 1
else
NTDM = LTROT + 1
endif
Output files
- direct program output
-
enerJ
static observables
- ...
- results of analysis
-
grR0
local Green function
-
gr_n_m
(integers 1<=n,m<=L) k-dependent Green function
- similarly with
charge
and spin
instead of gr
--
NilsBluemer - 20 Mar 2011