Physics 252 Pledged Homework

 

Due Friday April 17, Noon

1. Read the notes below on constructing a "leapfrog" spreadsheet to solve Schrödinger's equation for a potential, and test it out by finding the first few eigenvalues (bound state energies) for the simple harmonic oscillator potential.

Solving f"(x) = (V(x) -E)f(x) by the leapfrog method.

Take a set of points x0, x1, x2, x3, ... as usual equally spaced dx apart, so xi+1 = xi + dx.

We shall find f(x) and its second derivative f"(x) at this series of points, using the value of the first derivative f(x) at the points midway between this set of points, for example at the point x1+1/2 which is at x1 + 0.5dx.

The reason for this procedure is as follows: if we know f'(x) at x1+1/2, and f(x) at x1, then to a very good approximation f(x2) = f(x1) + dx*f'(x1+1/2), since this is using the slope of the function midway between x1 and x2 as a straight line approximation to the function in that interval.

We actually find f'(x) at these midpoints by using the value of f"(x) at the original set of points. It's really the same trick: f'(x2+1/2) = f'(x1+1/2) + dx* f"(x2). Again, this approximates a function with a straight line having the slope in the middle of the interval.

The method is implemented as follows:

Take five columns of the spreadsheet and label them as below:

A

B

C

D

E

position

V(x)-E

wavefn

2nd deriv

deriv at 1/2

x j

v(x j)-E

f(x j)

f"(xj)

f'(xj+1/2)

(You don't need to use A, B, C, D, E -- I just did that for convenience here)

The first two columns are filled in as in the previous spreadsheet.

The initial entry (which we take to be row 3) in the third column is f(x0), in D3 (say) =B3*C3 (using f" = (V-E)f), in E3 =f'(x0) + 0.5*dx*D3.

Going on to the fourth row, in C4 put =C3 + dx*E3 ... and you figure out the rest.

Once you have set up this spreadsheet, and got it to run, use it to find the first three energy levels of a harmonic oscillator. Compare your results with:
(a) the exact answer,
(b) the results using the simpler spreadsheet. State approximately what level of accuracy the two spreadsheets give for these problems. How is this improved if you double the number of rows used in the spreadsheet?

2. Square Wells

(a)Use the new spreadsheet, or the one from last week's homework if you prefer, to find the lowest four energy levels in a one-dimensional potential that consists of two square wells, each of depth 20, and width 3.14159, a distance S apart. Find how the energy difference between the two lowest energy levels varies with S: find it for S = 0.1, 1, 5, 10. How does the energy in these lowest states relate to the energy levels of a particle in one of these wells by itself? Plot some of the wavefunctions.

(b)An ionized hydrogen molecule looks to the single electron orbiting the two protons like a double well. Do you think the wavefunction in the lowest state is odd or even about the central point? Explain.

(c)Now consider a potential in one dimension which has three such wells. Take the separation S=1. Find the lowest energy wavefunction, and discuss the energy difference to the next wavefunctions. What happens as S is varied? This model is a first attempt at understanding electronic energy levels in metals.

3. Muonic atoms

Recall that the muon is a particle having the same charge as an electron, but with a mass 207 times greater. The muon has a lifetime of only two microseconds, but that's long enough, if trapped in matter, for it to drop through a succession of Bohr-type orbits around a nucleus, emitting photons corresponding to the energy differences between the orbits.

(a) Suppose a muon is orbiting a hydrogen nucleus, and drops from n = 3 to n = 2 then to n = 1. What would be the frequency of the photons emitted? What kind of photons are these? (Take the nucleus to have infinite mass.)

(b) Describe how your answer would change if you took into account the fact that the nucleus has a finite mass. Is this more important for this muonic atom than for a regular hydrogen atom?

(c) Will the presence of the muon affect the atom's electron significantly?

(d) Now consider a muon captured by a Uranium nucleus, of charge 92 and nuclear radius 7.10-15 meters. Find the radii of the lowest three Bohr orbits. Don't just use the formula in this part - derive from first principles!