[ALPS-users] Conserve quantum number for composite unit cell

Matthias Troyer troyer at phys.ethz.ch
Wed Jan 18 21:47:41 CET 2017


I would suggest to use one quantum number for both Sz and N instead of two or three different quantum numbers. 

> On 18 Jan 2017, at 12:24, Axel Gagge <axel.gagge at fysik.su.se> wrote:
> 
> With extra quantum number, fixed n.o. excitations: runs in 124 s, E0, trunc_error = -1.0000000000001179, 0.00082901673330108058
> With no extra quantum number, no fixed n.o. excitations (grand canonical...): runs in 24 s, E0, trunc_error = -0.99999999999997691, -1.3941131613197575e-16
> 
> 
> I will proceed with the grand canonical and this is certainly enough, but is there some way around this slow-down? It would be nice to be able to conserve any charge which correspond to an abelian symmetry.
> 
> Thanks,
> /Axel
> 
> Från: Comp-phys-alps-users <comp-phys-alps-users-bounces at lists.phys.ethz.ch> för Matthias Troyer <troyer at phys.ethz.ch>
> Skickat: den 18 januari 2017 18:47
> Till: comp-phys-alps-users at lists.phys.ethz.ch
> Ämne: Re: [ALPS-users] Conserve quantum number for composite unit cell
>  
> Have you tried using a grandcanonical instead of fixing the quantum number?
> 
>> On 18 Jan 2017, at 09:15, Axel Gagge <axel.gagge at fysik.su.se> wrote:
>> 
>> ​It wasn't very hard to define such a quantum number, but it seems like it actually makes the simulation run much slower. I suppose this is because another quantum number needs to be stored which leads to a larger basis size again? On the other hand, I can think of no way to properly define the spin operators without using the quantum numbers S, Sz... ​
>> 
>> Many thanks,
>> Axel
>> 
>> Från: Comp-phys-alps-users <comp-phys-alps-users-bounces at lists.phys.ethz.ch> för Matthias Troyer <troyer at phys.ethz.ch>
>> Skickat: den 17 januari 2017 19:42
>> Till: comp-phys-alps-users at lists.phys.ethz.ch
>> Ämne: Re: [ALPS-users] Conserve quantum number for composite unit cell
>>  
>> You need to define a new combined quantum number, which is changed by the training and lowering operators for both types of particles.
>> 
>> 
>>> On 17 Jan 2017, at 08:54, Michele Dolfi <dolfim at phys.ethz.ch> wrote:
>>> 
>>> I hope somebody else could give you a better answer, but I don’t think this is possible out of the box, and I’m not 100% sure it is at all possible in the way you think about it.
>>> 
>>> What I would try is away to play around with a new model definition where you make a composite quantum number explicit.
>>> The "fermion Hubbard” and “alternative fermion Hubbard” could be of inspiration for how to conserve N and/or Sz instead of Nup and/or Ndown.
>>> 
>>> 
>>> Best,
>>> Michele
>>> 
>>> 
>>> --
>>> ETH Zurich
>>> Dr. Michele Dolfi
>>> Institute for Theoretical Physics
>>> HIT G 32.4
>>> Wolfgang-Pauli-Str. 27
>>> 8093 Zurich
>>> Switzerland
>>> 
>>> dolfim at phys.ethz.ch
>>> www.itp.phys.ethz.ch
>>> 
>>> +41 44 633 78 56 phone
>>> +41 44 633 11 15 fax 
>>> 
>>> 
>>> 
>>> 
>>> 
>>> 
>>>> On Jan 17, 2017, at 4:10 PM, Axel Gagge <axel.gagge at fysik.su.se> wrote:
>>>> 
>>>> ​Hello, I just wanted to ask if it is possible to preserve a combination of quantum numbers such as the total number of excitations:
>>>> 
>>>> N = \sum_{bosonic sites} n_j + 1/2 \sum_{spinful sites} s_z
>>>> 
>>>> and what the syntax for this would be? I guess that one option would be to define a new quantum number and let Splus raise it by 1/2 but I just thought that I would check this with you first!
>>>> 
>>>> Otherwise, thank you for an excellent library! I have had no problems so far with the new release.
>>>> 
>>>> Thanks,
>>>> /Axel Gagge
>>>> First year PhD Student, Stockholm University
>>>> 
>>>> 
>>>> ----
>>>> Comp-phys-alps-users Mailing List for the ALPS Project
>>>> http://alps.comp-phys.org/
>>>> 
>>>> List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users
>>>> Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
>>>> 
>>>> Unsubscribe by writing a mail to comp-phys-alps-users-leave at lists.phys.ethz.ch.
>>> 
>>> 
>>> 
>>> ----
>>> Comp-phys-alps-users Mailing List for the ALPS Project
>>> http://alps.comp-phys.org/
>>> 
>>> List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users
>>> Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
>>> 
>>> Unsubscribe by writing a mail to comp-phys-alps-users-leave at lists.phys.ethz.ch.
>> 
>> <my_models.xml><my_lattices.xml><Test of iDMRG.ipynb>
>> 
>> ----
>> Comp-phys-alps-users Mailing List for the ALPS Project
>> http://alps.comp-phys.org/
>> 
>> List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users
>> Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
>> 
>> Unsubscribe by writing a mail to comp-phys-alps-users-leave at lists.phys.ethz.ch.
> 
> 
> 
> ----
> Comp-phys-alps-users Mailing List for the ALPS Project
> http://alps.comp-phys.org/
> 
> List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users
> Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
> 
> Unsubscribe by writing a mail to comp-phys-alps-users-leave at lists.phys.ethz.ch.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.phys.ethz.ch/pipermail/comp-phys-alps-users/attachments/20170118/548bc123/attachment.html>


More information about the Comp-phys-alps-users mailing list