Kronig-Penney Model: Plotting Energy Bands in Excel

  • Thread starter Thread starter Fez07
  • Start date Start date
  • Tags Tags
    Model
Fez07
Messages
7
Reaction score
0
Hi there,
I want to plot a graph of the right hand side of the equation VS the energy, E, in eV for a periodic array of square wells each of depth 300eV, width 0.05nm, and period 0.1nm using SI units and MICROSOFT EXCEL spreadsheet to determine the allowed energy bands and band gaps:

cos(kp)=F=cos(α(p-w))cosh(βw)+((β squared - α squared)/2αβ))sin(α(p-w))sinh(βw)


α = (1/h bar)* square root(2mE)

β = (1/h bar)* square root(2m(V-E))

w = width of barrier separating two adjacent potential wells
p = period
m = mass of an electron
V = well depth

Solutions are only obtained if F is between -1 and 1.

I have no idea how to do this so any help would be much appreciated.

Thanks in advance.
 
Physics news on Phys.org
First, you need to set up your spreadsheet with the variable quantities and constants. We insist that you attempt a solution before we help you.
 
Thanks for replying. I set up the EXCEL spread sheet with columns of:

V = 300*1.6*10^-19 (energy in J)
p = 0.1*10^-9 (period in m)
m = 9.1094*10^-31 (mass of the electron in kg)
h_bar = 1.05*10^-34
w = 0.05*10^-9 (width of separating barrier in m)

E = alpha^2*h_bar^2 / (2*m)
β = (1/h_bar)* square root(2m(V-E))

and the RHS of the equation.

But i don't know what to do from here.
Please help me.

Thanks
 
So if you've set up your columns, plot 3 lines, \pm 1 and the RHS vs E. The areas where the values are greater than 1 are solutions where there are bands. Later I'll plot up your solutions and see what I get.
 
Hi again,
I still can't figure it out. I'm getting a "Circular Reference" error between the alpha and E values..?

Thanks

Fez07
 
Hi there, it's me again:

I've figured out what i was doing wrong and gotten the solution.

Thanks for the help. I couldn't have done it otherwise.

Two thumbs up.

Fez07
 
##|\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