Hi Fernando,
If you are using a recent version of ALPS, e.g. 2.2.0b3, you can use the new dmrg application `mps_optim`.
It works the same way as the `dmrg` application, but additionally allows you to measure correlation functions for i<j, which very often is good enough.
MEASURE_HALF_CORRELATIONS[cdag_up c_up]="cdag_up:c_up"
Best regards,
Michele
--
ETH Zurich
Michele Dolfi
Institute for Theoretical Physics
HIT G 32.4
Wolfgang-Pauli-Str. 27
8093 Zurich
Switzerland
+41 44 633 78 56 phone
+41 44 633 11 15 fax
Hi, Fernando,
Assuming that you use DMRG simulation since you specify 'MAXSTATES'.
(2014/06/12 0:10), Fernando Iemini wrote:
I've using the ''alternative fermion Hubbard'' model, with open boundary
conditions, and analyzing its single-particle correlations, <cdag_up:c_up>.
I'm not being able however to compute these observables for larger
lattices, e.g., L=~160. I'm getting a "segmentation fault" message, does
anybody know why is happening this error, maybe due to memory? If so, what
would be good parameters in order to avoid it?
Here is the script i'm using:
LATTICE="open chain lattice"
MODEL="alternative fermion Hubbard"
CONSERVED_QUANTUMNUMBERS="N,Sz,Nup,Ndown"
N_total=160
Sz_total=0
Nup_total=80
Ndown_total=80
mu=0
t=1
NUMBER_EIGENVALUES=1
L=160
MAXSTATES=100
MEASURE_CORRELATIONS[cdag_up c_up]="cdag_up:c_up"
{U = 9.000000e-01,V=5.500000e-01,TRUNCATION_ERROR=1e-7}
This is memory issue. The problem is that 'MEASURE_CORRELATIONS'
collects all pair sites. Observing ~10000 measurements costs too much
memory.
Not tested, but I believe that you can avoid this issue by tweaking
models.xml and lattices.xml. Change the name of some specific operator
and just calculating between them.
--
Ryo IGARASHI, Ph.D.
rigarash@issp.u-tokyo.ac.jp
tel: 04-7136-5343, fax: 04-7136-3441
OpenPGP fingerprint: BAD9 71E3 28F3 8952 5640 6A53 EC79 A280 6A19 2319