Deriving c

Screen Shot 2017-06-14 at 13.50.39Physics has its occasional eureka moment. James Clerk Maxwell had one such moment in 1862 when he discovered that the equations he had found describing electric and magnetic waves combined to describe light. The speed of light c, this article will demonstrate, can easily be found from Maxwell’s equations. As this derivation is an undergraduate degree exam question each step will be described qualitatively, but a brief glance at the vector operators guide in Resources is recommended.

One of the most profound 19th century discoveries was Faraday’s finding that electric currents passing through wires induce magnetic field loops, as shown below.Screen Shot 2017-06-14 at 15.16.15By 1862 Maxwell had shown that the gradient of the magnetic field H was directly proportional to the curl (or rotationality) of the electric field E. This essentially means that electricity and magnetism are one and the same, and a magnetic field flowing around a wire will equally induce a current through it, as described by:
Screen Shot 2017-06-14 at 15.20.08Where μ0 is the magnetic permeability constant, defined later. Taking the curl of both sides gives:Screen Shot 2017-06-14 at 15.20.14The curl of H (∇ x H) is defined by Ampére’s law as:Screen Shot 2017-06-14 at 15.29.30Where σ is the conductivity of the medium and ε₀ is the electric permittivity constant. Inserting this expression above gives:
Screen Shot 2017-06-14 at 15.29.36Using the following identity the Laplacian can be extracted. However, imposing the vacuum condition means that σ = 0 and the divergence of E (∇ . E) = 0.Screen Shot 2017-06-14 at 15.31.03Thus:Screen Shot 2017-06-14 at 15.31.11What has been derived is a wave equation, describing an electromagnetic wave, where μ0 = 1.26 x 10-6 N A⁻² and ε0 = 8.85 x10-12 F m⁻¹. Notice the similarity to the general equation for a wave:Screen Shot 2017-06-14 at 15.33.02In this equation v is the speed of the wave. The speed of an electromagnetic wave can therefore be found:Screen Shot 2017-06-14 at 15.33.41Screen Shot 2017-06-14 at 15.33.52This gives v = ~ 3 x10⁸ m s⁻¹, the speed of light. Thus, light can be treated as an electromagnetic wave. The century and a half since Maxwell first worked out this derivation has corroborated and built upon it.

Infinite Wells

Screen Shot 2017-06-12 at 17.57.14

The discovery of the wave nature of particles led Erwin Schrödinger to formulate his equation of the wavefunction in 1926, part of his quest to find a more reliable description of electrons in the atom. This article will lay out the mathematical basis of the wavefunction and its implications for quantum mechanics.

A bead on a string that is fixed at both ends is an example of an infinite potential well. It may move along one direction (call it x) between two limits (call these 0 and L). While a bead has no natural preference for one point on the string over any other, the following derivation will show that electrons do. In a one-dimensional well, their positions fixed between two limits, electrons ‘prefer’ to be in the middle and never occur at either edge. Let’s begin.

The general wavefunction Ψ (pronounced ‘sigh’) of an electron in a one-dimensional well is defined as:Screen Shot 2017-06-09 at 14.58.54A is the maximum height of the wave, k the wavenumber which relates to its wavelength, and x ranges from 0 to L as with the bead on a string.
Screen Shot 2017-06-09 at 15.07.33At this point it’s worth clarifying what the wavefunction is. Although it appears to be fundamental, its physical meaning is unknown, and is the subject of various interpretations which are beyond the remit of this article. Its square, however, is simple and measurable. For a well with one electron, the integral of Ψ² yields the probability of finding the electron in a particular region of the well.
Screen Shot 2017-06-12 at 15.36.07Probability is a familiar concept. In the case of a bead on the string, the probability of finding the bead off the string is zero and thus Ψ² = 0 at these points. Waves, unlike particles, are continuous in space, meaning they have no sudden breaks, and so if Ψ² = 0 just above L and just below 0 then the same must true at these points. So:Screen Shot 2017-06-12 at 15.36.14In case it isn’t apparent, this has shown why particles, whose positional probabilities relate to wave functions, can never appear at the edges of infinite wells. This has the consequence of allowing for k to be defined in simpler terms. If:Screen Shot 2017-06-12 at 15.50.05Then, since sin(x) only equals zero at points where x = n𝜋, then:Screen Shot 2017-06-12 at 15.53.12Thus k can be replaced to give:Screen Shot 2017-06-12 at 15.58.23At this point, it’s worth reflecting on n. n has integer values 1, 2, 3, 4… and its introduction is a consequence of treating the electron probability function as a wave, and by restricting both ends to zero. Since n relates to energy, which is shown later, this is the reason electrons are restricted to precise energy levels in atoms.

Finding A requires only the deduction that the probability of finding the particle in the entire well must be 1, so:Screen Shot 2017-06-12 at 16.17.22The above can be rewritten as:Screen Shot 2017-06-12 at 16.19.48Screen Shot 2017-06-12 at 16.47.41This follows a prolonged derivation, available here. Thus the full wavefunction is:
Screen Shot 2017-06-12 at 16.49.48For n = 1 and L = 1, the probability P of being in a region from x₁ to x₂ is given as:Screen Shot 2017-06-12 at 16.53.04Screen Shot 2017-06-12 at 16.55.41Below this wavefunction is graphed and regional probabilities are given.
figure_1~40% of the time the electron will be found in the middle fifth of the well, while ~10% of the time will it appear in the two fifths at both edges. Rather profoundly, the bead on a string is no different. It too obeys Schrödinger’s equation and has its own wavefunction just like the one above. The difference has to do with the associated energies involved, which as mentioned relate to n. Cranking up n yields the classical picture where there is no positional preference, as with beads on strings.figure₂figure_3In the n = 5 graph, there is equal probability of finding the electron in either fifth, while for the n = 20 graph, the classical picture is re-emerging where the electron can be virtually anywhere. As n is increased into the thousands and millions, the graphs are filled with yellow; the particle can be found anywhere.

The equation relating energy, mass and n is derived here, and a basic infinite well simulator, designed to produce the above graphs, can be found in Resources.

Dark Energy


In the 20th century as galaxies further afield were resolved astronomers found that the rate of cosmic expansion has increased over time; the Big Bang is, in fact, speeding up. This article will briefly cover the evidence leading Adam Riess et al. to this conclusion in their 1998 Nobel-winning paper, as well as its profound implications on the future of the universe.

The first method of finding the distance to another galaxy was provided by Henrietta Leavitt in 1912 when she, having tracked thousands of stars in the Magellanic cloud, found a relation between the period and luminosity of Cepheid variable stars. Distinguishable by the signficant He²⁺ content in their spectra, this transparent gas ionises at high temperatures to form the opaque He³⁺, making the star dimmer. A resulting expansion following this temperature increase causes the star to cool and contract, reverting to He²⁺, upon which this process repeats and the star oscillates. By measuring the period and relative brightness of Cepheid variables their distances, and those of the galaxies to which they belong, can be found. Another way of finding the distance to another galaxy is to measure the brightness of Type 1A supernovae. These are the result of white dwarfs gobbling up nearby stars, a process that the Pauli exclusion principle predicts must reach a critical instability at ~1.44 solar masses. Type 1A supernovae therefore share the same luminosity, making them a reliable standard candle, whose distance relates to their relative brightness. Generally, these and other methods are used in combination in order to improve precision and minimise distance uncertainties.

The velocity of a galaxy can be found by measuring the redshift of its spectrum. Elementary peaks are compared with their known laboratory wavelengths, and redshift z is defined as:Screen Shot 2017-06-06 at 14.28.06Where λ₀ is the known wavelength of a peak and Δλ is by how much it is shifted in a measured spectrum. The velocity relates to z by:
Screen Shot 2017-06-06 at 14.28.10The galaxies Hubble observed had, on average, a positive redshift, and the further away one was the redder its spectrum, and so the faster it was receding. More profoundly, Hubble showed that this would be the case at every point in space, meaning that the redshift was not purely a Doppler effect but is caused the expansion of space between the galaxies while the light is in transit. Hubble defined a rate of constant expansion H₀, which is presently defined as 67.8 km per second per megaparsec, and relates to velocity by:Screen Shot 2017-06-06 at 14.53.47Since c and H₀ are taken to be constants, the graph of this equation is simply a straight line. Observations agree with this on the relatively small scales Hubble observed in 1929, but as galaxies are mapped at ranges up to 12,000 megaparsecs, ~90% of the radius of the known universe, Hubble’s straight line bends, as shown below, indicating a much lower expansion rate in the past.
hubbleacc-2This is precisely the opposite of what is expected in a matter-dominated universe. In such a universe the force of gravity is sufficient to eventually slow down any expansion and cause a contraction. To overcome this force requires energy, and thus dark energy is introduced, with the word ‘dark’ conveying no more information than its being unknown. So the question of whether the universe will expand forever or eventually stop and collapse in on itself is no more than a question on the ratio of matter to dark energy. This cosmic density term, aptly given the letter Ω, is conveniently defined such that if Ω > 1 gravity prevails and the universe will collapse in on itself while if Ω < 1 dark energy prevails and the universe expands forever. Current measurements put a value on Ω at almost exactly 1, predicting that the universe will eventually slow down but never quite stop. Further research is underway to more precisely measure Ω and reduce what are considerable uncertainties.

For now the acceleration continues. It is predicted that within a few million years many galaxies now visible will have vanished from the night sky, since they will be receding from the Milky Way faster than the speed of light. In several billion years all distant galaxies will have vanished, leaving astronomers stranded in a much smaller universe, resembling the one in which they thought they lived at the start of the 20th century.

Riess, Adam et al. (1998). “Observational evidence from supernovae for an accelerating universe and a cosmological constant”. Astronomical Journal. 116 (3): 1009–38 [PDF]
Mazzali, P. A. et al. (2007). “A Common Explosion Mechanism for Type Ia Supernovae”. Science. 315 (5813): 825–828 [PDF]
Michael Richmond. (2000). Estimating Distances to Far-away Galaxies. [HTML]
Kirshner, R. P. (2003). “Hubble’s diagram and Cosmic expansion”. National Academy of Sciences. [HTML]
Peebles, P. J. E. (2003). “The cosmological constant and dark energy”. Reviews of Modern Physics. 75 (2): 559–606. [PDF]
Riess, Adam et al. (2004). “The Expanding Universe: From Slowdown to Speed Up”. Scientific American. 290 (2). [HTML]
De Bernardis, P. et al. (2000). “A flat Universe from high-resolution maps of the cosmic microwave background radiation”. Nature. 404 (6781): 955–9. [PDF]
Krauss, L. M.; Scherrer, R. J.  (2008). “The End of Cosmology?” Scientific American. 82. [PDF]

An Integral Quirk


Calculus revolutionised maths and physics in the 17th century, allowing Isaac Newton to form his equations of motion. This article will lay out the physical meaning of its primary components, differentiation and integration, and show their power by, from first principles, deriving Newton’s equations of motion.

Differentiation finds the slope of a curve. In others words, by how much y changes when x is increased by 1. For any function:Screen Shot 2017-06-01 at 20.57.34The slope (or derivative) is given by:Screen Shot 2017-06-01 at 20.57.41For f(x) = x, as plotted below, the gradient f'(x) = 1, meaning that y increases the same as x, which makes sense.Screen Shot 2017-06-01 at 21.06.34Differentiation is therefore the rate of change. The rate of change of position is velocity, and the rate of change of velocity is acceleration.

Integration, on the other hand, finds the curve from the slope. It is the opposite of differentiation. For any function:
Screen Shot 2017-06-02 at 22.19.34The curve (or integral) is given by:Screen Shot 2017-06-02 at 22.23.17.pngSo for f'(x) = 1, the integral f(x) = x + c. The c is an uncertainty constant. This arises because knowing the slope isn’t enough to know where the curve starts and ends, just as knowing the velocity of a car gives no information about its position. Yet the integral of velocity is position, and the integral of acceleration is velocity. This last one is where we begin.

Define acceleration as a. From first principles the velocity v will be the integral of a with respect to t. Performing this integral returns:Screen Shot 2017-06-02 at 23.00.37Screen Shot 2017-06-02 at 23.01.50.pngIn other words, the a gains a t term (as expected) and there is an uncertainty constant u. Now what is u? Well, if t = 0, we can see that v = u, and so u is the initial velocity, before time commenced. This is one of Newton’s equations of motion, but he has another.

If we integrate again, we will arrive at position:
Screen Shot 2017-06-02 at 23.09.38Screen Shot 2017-06-02 at 23.12.59Where, again, x₀ is the uncertainty constant and is equal to x at t = 0. Thus it is the initial position. By redefining xx₀, (i.e. distance between final and starting positions) as the displacement s, the second of Newton’s equations of motion is derived.Screen Shot 2017-06-02 at 23.16.01

Dark Matter


To approach the topic of dark matter is to approach the frontier of modern physics. This article presents briefly the evidence for this hypothetical substance, after a short overview on why it is widely considered to exist.¹

The first person to suggest the presence of dark matter was the Dutch astronomer Jacobus Kapteyn in 1922.² By this time telescopes had grown powerful enough to resolve galaxies, providing a testing ground for general relativity. Yet, while Einstein’s theory acquired overwhelming evidence closer to home, galaxies fell far short of its predictions. To account for this, Vera Rubin and Kent Ford published in 1980 a model for a new substance with the required mass and which did not emit light, that had to be six times more abundant in the universe than visible matter.³ The evidence leading to such a conclusion is summarised as follows:

  • First, the rotation of galaxies indicates a spread of mass throughout rather than a central concentration. Kepler’s law predicts that the average velocity of stars orbiting the centre of a galaxy is given by:Screen Shot 2017-05-31 at 23.20.46Where r is the distance from the centre. Instead of dropping, v is observed through optical spectroscopy to be constant.⁴ By equating the gravitational and centripetal forces the implications of this on galactic density are revealed:Screen Shot 2017-05-31 at 23.29.16Screen Shot 2017-05-31 at 23.26.21Screen Shot 2017-05-31 at 23.26.35Where G is the Newtonian constant, and M the mass of a galaxy whose derivative with respect to distance is dM/dr. In words, constant velocity requires almost-constant density meaning large amounts of matter must exist on the very edges of galaxies, far more than shows in optical imaging.
  • Second, and on an even grander scale, galaxy clusters appear to have many times more kinetic energy than their visible mass allows. For a galaxy of mass m moving in a cluster of mass M a distance r from the centre, its kinetic and potential energies will be:Screen Shot 2017-05-31 at 23.53.38Screen Shot 2017-05-31 at 23.53.44But for a stable system, the centripetal and gravitational forces must be equal, returning expressions for K and P:
    Screen Shot 2017-05-31 at 23.55.27Screen Shot 2017-05-31 at 23.55.33Screen Shot 2017-06-01 at 00.02.12.pngSo the potential should be equal to twice the kinetic energy, yet for observed clusters, the kinetic energy appears several times greater than the potential, implying an unseen source of extra gravitational energy.⁵
  • Third, gravitational lensing. Einstein showed that gravitational fields can deflect light, where the angle of deflection is:Screen Shot 2017-05-31 at 22.23.43Where d is the distance between the foreground and background galaxies, found by measuring their redshifts. The results again imply several times more mass in a galaxy cluster than is visible, as would be explained by dark matter.⁶
¹ Kroupa, P.; et al. (2010). “Local-Group tests of dark-matter Concordance Cosmology: Towards a new paradigm for structure formation”. p. 14. (PDF)
² Kapteyn, Jacobus Cornelius (1922). “First attempt at a theory of the arrangement and motion of the sidereal system”. p. 1. (PDF)
³ Rubin, Vera C.; Ford, W. Kent, Jr. (February 1970). “Rotation of the Andromeda Nebula from a Spectroscopic Survey of Emission Regions”. (PDF)
⁴ Corbelli, E. & Salucci, P. (2000). “The extended rotation curve and the dark matter halo of M33”. p. 1. (HTML)
⁵ Richmond, M (2004) “Using the virial theorem: the mass of a cluster of galaxies”. (HTML)
⁶ Taylor, A.N.; et al. (1998). “Gravitational Lens Magnification and the Mass of Abell 1689” (PDF)