Hi Kuang-Shing Chen,
first you need to change to the Hubbard model density of states (for 1d, 2d, 3d, …). You can do that, for example, by using the loop for the general density of states and specifying a DOS file that has the DOS for your lattice.
The parameter name is "DOSFILE" (see class FSDOSHilbertTransformer).
To adjust the filling: there is a parameter 'MU' corresponding to the chemical potential. If you change mu you change the filling. The ALPS code does not have a self-consistency for constant particle number (but that's a small thing that we could add if you really need it).
Emanuel
On Dec 13, 2011, at 4:21 PM, Kuangshing Chen wrote:
Thank you Emanuel,
Now I understand the code deeper. If I want to calculate the DMFT Hubbard Model, how do I change to filling?
Regards,
Kuang-Shing Chen
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
--
Kuang-Shing Chen