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 want to test the effects of extra terms in the Hamiltonian, I need to first ensure that this is working as expected. Code is provided below:(adsbygoogle = window.adsbygoogle || []).push({});

In this case the variable a is set to 2.46. Would anyone know if there is something in the syntax I am doing incorrectly, or perhaps the dispersion is being calculated wrong? Any help is appreciated. Thanks!Code (Text):f[kx_, ky_,

t_] := -t E^(-I kx a) (1 + 2 E^(I (3 kx a)/2)*Cos[Sqrt[3]/2 ky a]);

GrapheneHam[kx_, ky_, t_] :=

ComplexExpand[{{0, f[kx, ky, t]}, {Conjugate[f[kx, ky, t]], 0}}];

Energies[kx_, ky_, t_] := Eigenvalues[GrapheneHam[kx,ky,t]];

Plot3D[Energies[kx,ky, 2.8], {kx, -Pi/a, Pi/a},{ky, -Pi/a, Pi/a}]

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Mathematica Plotting Graphene in Mathematica

Loading...

Similar Threads for Plotting Graphene Mathematica |
---|

Plot a stereographic hypercube in a Mac's Grapher.app? |

Mathematica Errors when plotting but not when evaluating |

Mathematica Color Contour Plots |

Mathematica Taking forever to plot |

Mathematica Cannot Plot This Function in Mathematica |

**Physics Forums | Science Articles, Homework Help, Discussion**