Tight Binding Hamiltonian for Graphene

Lockoman
Messages
1
Reaction score
0
Hello, I am trying to write a program that will automate the creation of a tight binding Hamiltonian matrix for armchair cut graphene. However, I have almost no experience coding and would need some help to get started.

This would be assuming that the energy between nearest neighbor carbon atoms is t and everything else zero.

The only real info I have gathered myself is that for an armchair configuration, the number of atoms along a straight x or y line will be even. (as opposed to zigzag which would be odd).

Does anyone have any tips for me or examples of this being done?

Thank you very much!

P.S. I have a copy of Maple and would ideally like to use that as opposed to matlab, mathematica, etc
 
Physics news on Phys.org
. if possible.The best way to get started with this is to look into existing implementations of tight binding Hamiltonian matrixes for armchair cut graphene, as this will give you a good starting point. You can then modify the code to suit your needs. A few good resources to check out include: - The SourceForge project "Graphene-Tight-Binding" which provides a MATLAB implementation of a tight binding Hamiltonian matrix for armchair cut graphene. - The paper "Tight-Binding Model for Graphene: A First Principles Study" by M.-H. Bae et al. which provides a detailed description of the tight binding model and how it can be used to generate a Hamiltonian matrix for armchair cut graphene. - The paper "Band Structure and Electronic Properties of Armchair Graphene Nanoribbons" by L. K. Ang et al. which uses a tight binding approach to study the band structure of armchair cut graphene nanoribbons.Hopefully these resources will help you get started with creating your own Hamiltonian matrix for armchair cut graphene!
 
##|\Psi|^2=\frac{1}{\sqrt{\pi b^2}}\exp(\frac{-(x-x_0)^2}{b^2}).## ##\braket{x}=\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dx\,x\exp(-\frac{(x-x_0)^2}{b^2}).## ##y=x-x_0 \quad x=y+x_0 \quad dy=dx.## The boundaries remain infinite, I believe. ##\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dy(y+x_0)\exp(\frac{-y^2}{b^2}).## ##\frac{2}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,y\exp(\frac{-y^2}{b^2})+\frac{2x_0}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,\exp(-\frac{y^2}{b^2}).## I then resolved the two...
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top