Dear Joe, 

I am trying to investigate  a graphene using kwant and I would like to see if  any defect in the material (e.g. using del sys [a(0,0)] ) could make the  material becomes magnetic. Can this transformation be indicated by some changes in DOS ? I enclose the python codes.

Regards,

Dedi 


> Date: Tue, 29 Sep 2015 08:32:29 +0200
> From: joseph.weston@cea.fr
> To: setiabudidaya@hotmail.com
> CC: kwant-discuss@kwant-project.org
> Subject: Re: [Kwant] ldos.cumsum or ldos.cumprod
>
> Hi Dedi,
>
> > Hi there,
> > I tried the functions ldos.cunsum() and ldos.cumprod() and succeeded in obtaining the plot;but I do not know what the plots told me. Anyone can help ?
>
> Can you try and give some context as to what you are doing? It will be
> extremely difficult to help you if we do not know what it is that you
> are trying to do.
>
> Was there a previous message in this thread that I missed?
>
>
> > Below is the codes for the plot
> > data1 = []energies = np.linspace(1,4,101)for energy in energies: ldos = kwant.ldos(sys,energy) data1.append(ldos.cumsum())pyplot.figure()pyplot.plot(energies, data1)pyplot.xlabel("energy [t]")pyplot.ylabel("ldos")pyplot.show()
> > and for the other plot I just changed the function ldos.cumsum() with ldos.cumprod().
>
> Could you please submit the code either: in the body of the email but
> properly formatted (so that it can be easily cut/pasted and be a valid
> Python program), or put it as an attachment.
>
>
> Thanks,
>
> Joe