Altitude of a satellite

A satellite is launched into a circular orbit around the earth so that it orbits the planet once every $T$ seconds. The altitude $h$ above the Earth’s surface that the satellite must have is:

$h = \left(\frac{GMT^2}{4\pi^2}\right)^{\frac{1}{3}} - R$

where $G$ is Newton’s gravitational constant, $M$ is the mass of the Earth, and $R$ is the radius.

a) Write a piece of code which calculates the altitude $h$ (in metres) for a given value of $T$ (in seconds).

b) Use this code to calculate the altitude of satellites that orbit the Earth once a day (a “geosynchronous” orbit) and once every 45 minutes. What can you conclude from this final calculation?

The emission lines of hydrogen

There is a simple and famous formula for calculating the wavelengths $\lambda$ of the emission lines of the hydrogen atom.

$\frac{1}{\lambda} = R\left(\frac{1}{m^2} - \frac{1}{n^2}\right)$

where R is the Rydberg constant $R = 1.097\times 10^{-2}\mathrm{nm}^{-1}$ and $m$ and $n$ are positive integers and $n>m$.

a) Write a piece of code to calculate the first five transitions in the Lyman series ($m=1$, transitions to the ground state, emission in the UV-range), Balmer series ($m=2$, transitions to the first excited state, emissions in the visible region) and Paschen series ($m=3$, emissions in the infra-red).

Calculating Planck’s constant

When light is shone on the surface of a metal, the photons in the light can excite (transfer energy to) electrons in the metal and, sometimes, eject them from the surface into the free space above. The energy of the ejected electron can be calculated by measuring the minimum voltage $V$ that stops the electron moving.

We know that the energy of an ejected electron is equal to the energy of the photon that excited it minus the workfunction $\phi$ (which is the energy needed to remove it from the surface) and that the energy of a single photon is $hf$ where $h$ is Planck’s constant and $f$ is the frequency of light. Mathematically this can be expressed as:

$eV = hf - \phi$,

where $e$ is the charge of the electron.

a) Read in the photoelectric measurement data from the file Planck.txt. The first column contains frequencies $f$ is hertz and the second column contains voltages $V$. Use this data to plot $V$ vs $f$. Think about the plot type - does a scatter plot of line plot make most sense?

The least-squares method is very commonly used for fitting a polynomial to a set of data. As it is so prevalent in physics and engineering, you are encouraged to watch this video which give an intuitive and mathematical description of the method.

b) Fit a straight line (polynomial of degree one) to the data using the least-squares method implemented in numpy.polyfit. Overlay this line on your data points.

c) Using the gradient of the fitted line calculate a value for Planck’s constant. Compare this to values you can find online.

The Madelung constant

The Madelung constant gives the total electric potential felt by an atom in a solid. It depends on the charge and position of other nearby atoms.

Consider the compound sodium chloride. These are arranged on a cubic lattice, with sodium having a positive charge ($+e$) and chlorine having a negative charge ($-e$). If each atom position is given by integers $(i,j,k)$ then the sodium atoms are at positions where $i+j+k$ is even and the chlorine atoms are at positions where $i+j+k$ is odd.

For an atom at $i=j=k=0$, the Madelung constant $M$ can be approximated by using the following formulae:

$V_\mathrm{total} = \sum_{i,j,k} V(i,j,k) = \frac{e}{4\pi\epsilon_0a}M$

$V(i,j,k) = \pm\frac{e}{4\pi\epsilon_0r}$

where $r$ is the distance from the origin to the atom at position $(i,j,k)$ and $a$ is the atom spacing. The summation runs from $i,j,k=-L$ to $i,j,k=L$ but not including $i,j,k=0$ (otherwise the expression would “blow up”).

a) Write an expression (in Markdown/LaTeX) for the distance $r$ in terms of $i$, $j$, $k$ and $a$.

b) Calculate the Madelung constant for sodium chloride using a large a value as L as you can (so the code runs in about a minute or less). How does it compare with published values?