0 votes
asked by (250 points)

I want to know whether dmrg can realize hamiltonian contains:
b is the boson(cavity) and n is the number of electron at site i. b is interacted with all electron site. I don't know whether dmrg can calculate this .

1 Answer

0 votes
answered by (56.3k points)

Hi Weijie,
Yes, DMRG and ITensor can be used to study a Hamiltonian of this type (and any Hamiltonian which is a sum of local terms on a lattice).

But special care has to be taken with bosons, as you know, since one has to make their Hilbert space finite in some way. So you could use the Boson site set in ITensor and adjust the MaxOcc parameter to make sure your results are converged with the maximum boson occupancy per site.

Your other question about Hubbard-Holstein might be relevant for others reading this question, so I'm linking it here: http://itensor.org/support/1502/how-to-realize-hubbard-holstein-model?show=1502#q1502

commented by (250 points)
Hi miles
I find one example.The hamiltonian is H=He+Hph+He-ph
H_{e}=-t\sum_{i,\sigma}(c^{\dag}_{i\sigma}+h.c.)+U\sum_{i}n_{i\uparrow}n_{i\downarrow} \\
Hph=\Omega b^{dag}b \\
And I write a code. When I execute the code, the energy of every iteration is not always descending. It will be fluctuating, as -10,-11,-10,-9.... I try different parameters. However the result is same.
#include "itensor/all.h"
using namespace itensor;

using Holstein = MixedSiteSet<ElectronSite,BosonSite>;

    auto N = 20;

    auto U = 0.5;
    auto J = 1.0;
    auto g = 1.0;
    auto sigma = 1.0;

    auto sites = Holstein(N,{"ConserveNf=",true,

    auto ampo = AutoMPO(sites);
    for(int j=1;j  <= N-2; j+=2)
        ampo += -J,"Cdagup",j+2,"Cup",j;
        ampo += -J,"Cdagup",j,"Cup",j+2;
        ampo += -J,"Cdagdn",j+2,"Cdn",j;
        ampo += -J,"Cdagdn",j,"Cdn",j+2;
    ampo += -J,"Cdagup",1,"Cup",19;
    ampo += -J,"Cdagup",19,"Cup",1;
    ampo += -J,"Cdagdn",1,"Cdn",19;
    ampo += -J,"Cdagdn",19,"Cdn",1;
   for(int j=1;j < N; j += 2)
        ampo += U,"Nup",j,"Ndn",j;
    for(int j=1;j<N;j +=4)
    ampo += g,"Adag",2,"Nup",j;
    ampo += g,"Adag",2,"Ndn",j;
    for(int j=3;j<N;j +=4)
    ampo += -g,"Adag",2,"Nup",j;
    ampo += -g,"Adag",2,"Ndn",j;
    ampo += sigma,"Adag",2,"A",2;
    auto H = toMPO(ampo);

    auto sweeps = Sweeps(20);
    //very important to use noise for this model
    sweeps.noise() = 1E-6,1E-6,1E-8,1E-12;
    sweeps.maxdim() = 10,20,100,100,200;
    sweeps.cutoff() = 1E-10;

    auto state = InitState(sites);
    for(auto j : range1(N))
        if(j%4==1) state.set(j,"Up");
        else if(j%4==3) state.set(j,"Dn");
        else state.set(j,"0");
    auto psi0 = MPS(state);

    auto [energy,psi] = dmrg(H,psi0,sweeps,{"Quiet=",true});

    printfln("Ground State Energy = %.12f",energy);
commented by (56.3k points)
Hi, for cases like this I think it is best to compute local measurements (like of the density of particles) and print them out to a file, then plot the results to visualize what is going on. It could be that your calculation is getting stuck due to the details of the physics.

Welcome to ITensor Support Q&A, where you can ask questions and receive answers from other members of the community.

Formatting Tips:
  • To format code, indent by four spaces
  • To format inline LaTeX, surround it by @@ on both sides
  • To format LaTeX on its own line, surround it by $$ above and below
  • For LaTeX, it may be necessary to backslash-escape underscore characters to obtain proper formatting. So for example writing \sum\_i to represent a sum over i.
If you cannot register due to firewall issues (e.g. you cannot see the capcha box) please email Miles Stoudenmire to ask for an account.

To report ITensor bugs, please use the issue tracker.