On Tue, Dec 13, 2011 at 8:43 AM, Emanuel Gull <emanuel.gull@gmail.com> wrote:

Hi Kuang-Shing Chen,

> It sounds that you impose the particle-hole symmetry in G_s (the argument of it is beta -tau) and the minus sign is just the convention for G_s and there is another minus sign in the multiple_G in the function get_result() to cancel it. Is the particle-hole symmetry only true for the initial condition (no interaction)?

this is a special case for the self-consistency condition and, as you mention, it is only valid for Bethe lattice. We have only implemented it in the paramagnetic phase. You can find the implementation of it in F_selfconsistency_loop in the file selfconsistency.C

F is not the hybridization function, but something like it up to a minus sign, see Rev. Lett. 97, 076405 (2006) for its exact definition.

It's probably best if you start by running some of the tutorials for the hybridization expansion code. Tutorial 02 runs DMFT on a Bethe lattice with the hybridization expansion code in the antiferromagnetic case, I think this is exactly what you want.

Best,
Emanuel