Recent content by nuclearpasta

I'll post some snippets of what code in MATLAB would look like for a calculation like this. phase(i) = exp(1i*dot(k(:,index),R(:,i))) This gives you the phase you're going to multiply with the hopping parameter. You'd have to write a loop that runs from i = 1:3 in order to multiply each...

That's strange. When I do it I get that the two surfaces never touch each other, which is what is supposed to happen at the six corners of the hexagon.

I am encountering an issue when I plot attempting to plot 3d band structure of graphene in Mathematica. While the general shape and curvature looks fine, the cones at the K-points are not touching, which is an important qualitative electronic property of graphene. Since further down the line I...

If you're trying to plot this in something like MATLAB you can form the two block matrices and run a loop that adds the phase factors for the hopping parameters and then plot over those points if you know the region it spans (which in your case is just the Brillouin zone)