A Radically Modern Approach to Introductory Physics

Second Edition

Volume 1: Fundamental Principles

David J. Raymond

The New Mexico Tech Press

Socorro, New Mexico, USA

A Radically Modern Approach to Introductory Physics Second Edition

Volume 1: Fundamental Principles

David J. Raymond

Copyright © 2011, 2016 David J. Raymond

Second Edition

Original state, 17 May 2016

Content of this book available under the Creative Commons Attribution-Noncommercial-Share-Alike License. See http://creativecommons.org/licenses/by-nc-sa/3.0/ for details.

Publisher’s Cataloguing-in-Publication Data

Raymond, David J.

A radically modern approach to introductory physics, second edition:

Volume 1: Fundamental principles / David J. Raymond.

xii, 227 p.: ill. ; 28 cm

Preface -- Waves in one dimension -- Waves in two and three

dimensions -- Geometrical optics -- Special relativity -- Applications

of special relativity -- Acceleration and general relativity --

Matter waves -- Geometrical optics and Newton’s Laws -- Symmetry

and bound states -- Dynamics of multiple particles -- Rotational

dynamics -- Harmonic oscillator -- Constants, units and conversions

ISBN 978-1-938159-06-0 (print) -- 978-1-938159-07-7 (ebook)

1. Physics. 2. Waves. 3. Wave mechanics. 4. Mathematical

physics. 5. Relativity (Physics). I. Title.

QC23.R268 2016 v.1 (print) -- QC23.R268 2016eb v.1 (ebook)

A radically modern approach to introductory physics, second edition:

Volume 1: Fundamental principles / David J. Raymond.

xii, 227 p.: ill. ; 28 cm

Preface -- Waves in one dimension -- Waves in two and three

dimensions -- Geometrical optics -- Special relativity -- Applications

of special relativity -- Acceleration and general relativity --

Matter waves -- Geometrical optics and Newton’s Laws -- Symmetry

and bound states -- Dynamics of multiple particles -- Rotational

dynamics -- Harmonic oscillator -- Constants, units and conversions

ISBN 978-1-938159-06-0 (print) -- 978-1-938159-07-7 (ebook)

1. Physics. 2. Waves. 3. Wave mechanics. 4. Mathematical

physics. 5. Relativity (Physics). I. Title.

QC23.R268 2016 v.1 (print) -- QC23.R268 2016eb v.1 (ebook)

OCLC Number: 952471268 (print) – 952469817 (ebook)

Published by The New Mexico Tech Press, a New Mexico nonprofit corporation

This copy printed by CreateSpace, Charleston, SC

The New Mexico Tech Press

Socorro, New Mexico, USA

http://press.nmt.edu

To my wife Georgia and my daughters Maria and Elizabeth.

Preface

Preface to Second Edition

1 Waves in One Dimension

1.1 Transverse and Longitudinal Waves

1.2 Sine Waves

1.3 Types of Waves

1.3.1 Ocean Surface Waves

1.3.2 Sound Waves

1.3.3 Light

1.4 Superposition Principle

1.5 Beats

1.6 Interferometers

1.6.1 The Michelson Interferometer

1.7 Thin Films

1.8 Math Review — Derivatives

1.9 Group Velocity

1.9.1 Derivation of Group Velocity Formula

1.9.2 Examples

1.10 Problems

2 Waves in Two and Three Dimensions

2.1 Math Tutorial — Vectors

2.2 Plane Waves

2.3 Superposition of Plane Waves

2.3.1 Two Waves of Identical Wavelength

2.3.2 Two Waves of Differing Wavelength

2.3.3 Many Waves with the Same Wavelength

2.4 Diffraction Through a Single Slit

2.5 Two Slits

2.6 Diffraction Gratings

2.7 Problems

3 Geometrical Optics

3.1 Reflection and Refraction

3.2 Total Internal Reflection

3.3 Anisotropic Media

3.4 Thin Lens Equation and Optical Instruments

3.5 Fermat’s Principle

3.6 Problems

4 Special Relativity

4.1 Galilean Spacetime Thinking

4.2 Spacetime Thinking in Special Relativity

4.3 Postulates of Special Relativity

4.3.1 Simultaneity

4.3.2 Spacetime Pythagorean Theorem

4.4 Time Dilation

4.5 Lorentz Contraction

4.6 Twin Paradox

4.7 Problems

5 Applications of Special Relativity

5.1 Waves in Spacetime

5.2 Math Tutorial – Four-Vectors

5.3 Principle of Relativity Applied

5.4 Characteristics of Relativistic Waves

5.5 The Doppler Shift

5.6 Addition of Velocities

5.7 Problems

6 Acceleration and General Relativity

6.1 Acceleration

6.2 Circular Motion

6.3 Acceleration, Force, and Mass

6.4 Acceleration in Special Relativity

6.5 Accelerated Reference Frames

6.6 Gravitational Red Shift

6.7 Event Horizons

6.8 Problems

7 Matter Waves

7.1 Bragg’s Law

7.2 X-Ray Diffraction Techniques

7.2.1 Single Crystal

7.2.2 Powder Target

7.3 Meaning of Quantum Wave Function

7.4 Sense and Nonsense in Quantum

Mechanics

7.5 Mass, Momentum, and Energy

7.5.1 Planck, Einstein, and de Broglie

7.5.2 Wave and Particle Quantities

7.5.3 Non-Relativistic Limits

7.5.4 An Experimental Test

7.6 Heisenberg Uncertainty Principle

7.7 Problems

8 Geometrical Optics and Newton’s Laws

8.1 Fundamental Principles of Dynamics

8.1.1 Pre-Newtonian Dynamics

8.1.2 Newtonian Dynamics

8.1.3 Quantum Dynamics

8.2 Potential Energy

8.2.1 Gravity as a Conservative Force

8.3 Work and Power

8.4 Mechanics and Geometrical Optics

8.5 Math Tutorial – Partial Derivatives

8.6 Motion in Two and Three Dimensions

8.7 Kinetic and Total Momentum

8.8 Problems

9 Symmetry and Bound States

9.1 Math Tutorial — Complex Waves

9.2 Symmetry and Quantum Mechanics

9.2.1 Free Particle

9.2.2 Symmetry and Definiteness

9.2.3 Invariance

9.2.4 Compatible Variables

9.2.5 Compatibility and Conservation

9.2.6 New Symmetries and Variables

9.3 Confined Matter Waves

9.3.1 Particle in a Box

9.3.2 Barrier Penetration

9.3.3 Orbital Angular Momentum

9.3.4 Spin Angular Momentum

9.4 Problems

10 Dynamics of Multiple Particles

10.1 Momentum and Newton’s Second Law

10.2 Newton’s Third Law

10.3 Conservation of Momentum

10.4 Collisions

10.4.1 Elastic Collisions

10.4.2 Inelastic Collisions

10.5 Rockets and Conveyor Belts

10.6 Problems

11 Rotational Dynamics

11.1 Math Tutorial — Cross Product

11.2 Torque and Angular Momentum

11.3 Two Particles

11.4 The Uneven Dumbbell

11.5 Many Particles

11.6 Rigid Bodies

11.7 Statics

11.8 Problems

12 Harmonic Oscillator

12.1 Energy Analysis

12.2 Analysis Using Newton’s Laws

12.3 Forced Oscillator

12.4 Complex Exponential Solutions

12.5 Quantum Mechanical Harmonic Oscillator

12.6 Problems

A Constants, Units, and Conversions

A.1 SI Units

A.1.1 Derived Units

A.1.2 SI Multipliers

A.1.3 CGS or Centimeter-Gram-Second Units

A.1.4 Miscellaneous Conversions

A.2 Advice on Calculations

A.2.1 Substituting Numbers

A.2.2 Significant Digits

A.2.3 Changing Units

A.3 Constants of Nature

A.4 Properties of Stable Particles

A.5 Properties of Solar System Objects

Preface to Second Edition

1 Waves in One Dimension

1.1 Transverse and Longitudinal Waves

1.2 Sine Waves

1.3 Types of Waves

1.3.1 Ocean Surface Waves

1.3.2 Sound Waves

1.3.3 Light

1.4 Superposition Principle

1.5 Beats

1.6 Interferometers

1.6.1 The Michelson Interferometer

1.7 Thin Films

1.8 Math Review — Derivatives

1.9 Group Velocity

1.9.1 Derivation of Group Velocity Formula

1.9.2 Examples

1.10 Problems

2 Waves in Two and Three Dimensions

2.1 Math Tutorial — Vectors

2.2 Plane Waves

2.3 Superposition of Plane Waves

2.3.1 Two Waves of Identical Wavelength

2.3.2 Two Waves of Differing Wavelength

2.3.3 Many Waves with the Same Wavelength

2.4 Diffraction Through a Single Slit

2.5 Two Slits

2.6 Diffraction Gratings

2.7 Problems

3 Geometrical Optics

3.1 Reflection and Refraction

3.2 Total Internal Reflection

3.3 Anisotropic Media

3.4 Thin Lens Equation and Optical Instruments

3.5 Fermat’s Principle

3.6 Problems

4 Special Relativity

4.1 Galilean Spacetime Thinking

4.2 Spacetime Thinking in Special Relativity

4.3 Postulates of Special Relativity

4.3.1 Simultaneity

4.3.2 Spacetime Pythagorean Theorem

4.4 Time Dilation

4.5 Lorentz Contraction

4.6 Twin Paradox

4.7 Problems

5 Applications of Special Relativity

5.1 Waves in Spacetime

5.2 Math Tutorial – Four-Vectors

5.3 Principle of Relativity Applied

5.4 Characteristics of Relativistic Waves

5.5 The Doppler Shift

5.6 Addition of Velocities

5.7 Problems

6 Acceleration and General Relativity

6.1 Acceleration

6.2 Circular Motion

6.3 Acceleration, Force, and Mass

6.4 Acceleration in Special Relativity

6.5 Accelerated Reference Frames

6.6 Gravitational Red Shift

6.7 Event Horizons

6.8 Problems

7 Matter Waves

7.1 Bragg’s Law

7.2 X-Ray Diffraction Techniques

7.2.1 Single Crystal

7.2.2 Powder Target

7.3 Meaning of Quantum Wave Function

7.4 Sense and Nonsense in Quantum

Mechanics

7.5 Mass, Momentum, and Energy

7.5.1 Planck, Einstein, and de Broglie

7.5.2 Wave and Particle Quantities

7.5.3 Non-Relativistic Limits

7.5.4 An Experimental Test

7.6 Heisenberg Uncertainty Principle

7.7 Problems

8 Geometrical Optics and Newton’s Laws

8.1 Fundamental Principles of Dynamics

8.1.1 Pre-Newtonian Dynamics

8.1.2 Newtonian Dynamics

8.1.3 Quantum Dynamics

8.2 Potential Energy

8.2.1 Gravity as a Conservative Force

8.3 Work and Power

8.4 Mechanics and Geometrical Optics

8.5 Math Tutorial – Partial Derivatives

8.6 Motion in Two and Three Dimensions

8.7 Kinetic and Total Momentum

8.8 Problems

9 Symmetry and Bound States

9.1 Math Tutorial — Complex Waves

9.2 Symmetry and Quantum Mechanics

9.2.1 Free Particle

9.2.2 Symmetry and Definiteness

9.2.3 Invariance

9.2.4 Compatible Variables

9.2.5 Compatibility and Conservation

9.2.6 New Symmetries and Variables

9.3 Confined Matter Waves

9.3.1 Particle in a Box

9.3.2 Barrier Penetration

9.3.3 Orbital Angular Momentum

9.3.4 Spin Angular Momentum

9.4 Problems

10 Dynamics of Multiple Particles

10.1 Momentum and Newton’s Second Law

10.2 Newton’s Third Law

10.3 Conservation of Momentum

10.4 Collisions

10.4.1 Elastic Collisions

10.4.2 Inelastic Collisions

10.5 Rockets and Conveyor Belts

10.6 Problems

11 Rotational Dynamics

11.1 Math Tutorial — Cross Product

11.2 Torque and Angular Momentum

11.3 Two Particles

11.4 The Uneven Dumbbell

11.5 Many Particles

11.6 Rigid Bodies

11.7 Statics

11.8 Problems

12 Harmonic Oscillator

12.1 Energy Analysis

12.2 Analysis Using Newton’s Laws

12.3 Forced Oscillator

12.4 Complex Exponential Solutions

12.5 Quantum Mechanical Harmonic Oscillator

12.6 Problems

A Constants, Units, and Conversions

A.1 SI Units

A.1.1 Derived Units

A.1.2 SI Multipliers

A.1.3 CGS or Centimeter-Gram-Second Units

A.1.4 Miscellaneous Conversions

A.2 Advice on Calculations

A.2.1 Substituting Numbers

A.2.2 Significant Digits

A.2.3 Changing Units

A.3 Constants of Nature

A.4 Properties of Stable Particles

A.5 Properties of Solar System Objects

The idea for a “radically modern” introductory physics course arose out of frustration in the physics department at New Mexico Tech with the standard two-semester treatment of the subject. It is basically impossible to incorporate a significant amount of “modern physics” (meaning post-19th century!) in that format. It seemed to us that largely skipping the “interesting stuff” that has transpired since the days of Einstein and Bohr was like teaching biology without any reference to DNA. We felt at the time (and still feel) that an introductory physics course for non-majors should make an attempt to cover the great accomplishments of physics in the 20th century, since they form such an important part of our scientific culture.

It would, of course, be easy to pander to students – teach them superficially about the things they find interesting, while skipping the “hard stuff”. However, I am convinced that they would ultimately find such an approach as unsatisfying as would the educated physicist. What was needed was a unifying vision which allowed the presentation of all of physics from a modern point of view.

The idea for this course came from reading Louis de Broglie’s Nobel Prize
address.^{1}
De Broglie’s work is a masterpiece based on the principles of optics and special
relativity, which qualitatively foresees the path taken by Schrödinger and others in
the development of quantum mechanics. It thus dawned on me that perhaps
optics and waves together with relativity could form a better foundation for all of
physics, providing a more interesting way into the subject than classical
mechanics.

Whether this is so or not is still a matter of debate, but it is indisputable that such a path is much more fascinating to most college students interested in pursuing studies in physics — especially those who have been through the usual high school treatment of classical mechanics. I am also convinced that the development of physics in these terms, though not historical, is at least as rigorous and coherent as the classical approach.

After 15 years of gradual development, it is clear that the course failed in its original purpose, as a replacement for the standard, one-year introductory physics course with calculus. The material is way too challenging, given the level of interest of the typical non-physics student. However, the course has found a niche at the sophomore level for physics majors (and occasional non-majors with a special interest in physics) to explore some of the ideas that drew them to physics in the first place. It was placed at the sophomore level because we found that having some background in both calculus and introductory college-level physics is advantageous for most students. However, we allow incoming freshmen into the course if they have an appropriate high school background in physics and math.

The course is tightly structured, and it contains little or nothing that can be omitted. However, it is designed to fit into two semesters or three quarters. In broad outline form, the structure is as follows:

- Optics and waves occur first on the menu. The idea of group velocity is central to the entire course, and is introduced in the first chapter. This is a difficult topic, but repeated reviews through the year cause it eventually to sink in. Interference and diffraction are done in a reasonably conventional manner. Geometrical optics is introduced, not only for its practical importance, but also because classical mechanics is later introduced as the geometrical optics limit of quantum mechanics.
- Relativity is treated totally in terms of space-time diagrams – the Lorentz transformations seem to me to be quite confusing to students at this level (“Does gamma go upstairs or downstairs?”), and all desired results can be obtained by using the “space-time Pythagorean theorem” instead, with much better effect.
- Relativity plus waves leads to a dispersion relation for free matter waves. Optics in a region of variable refractive index provides a powerful analogy for the quantum mechanics of a particle subject to potential energy. The group velocity of waves is equated to the particle velocity, leading to the classical limit and Newton’s equations. The basic topics of classical mechanics are then done in a more or less conventional, though abbreviated fashion.
- Gravity is treated conventionally, except that Gauss’s law is introduced for the gravitational field. This is useful in and of itself, but also provides a preview of its deployment in electromagnetism. The repetition is useful pedagogically.
- Electromagnetism is treated in a highly unconventional way, though the endpoint is Maxwell’s equations in their usual integral form. The connection to relativity is exploited rather than buried. In particular, the seemingly simple question of how potential energy can be extended to the relativistic context gives rise to the idea of potential momentum. The potential energy and potential momentum together form a four-vector which is closely related to the scalar and vector potential of electromagnetism. The Aharonov-Bohm effect is easily explained using the idea of potential momentum in one dimension, while extension to three dimensions results in a version of Snell’s law valid for matter waves, from which the Lorentz force law is derived.
- The generation of electromagnetic fields comes from Coulomb’s law plus relativity (I borrowed from my graduate advisor Mel Schwartz’s text on electromagnetism here), with the scalar and vector potential being used to produce a much more straightforward treatment than is possible with electric and magnetic fields. Electromagnetic radiation is a lot simpler in terms of the potential fields as well.
- Resistors, capacitors, and inductors are treated for their practical value, but also because their consideration leads to an understanding of energy in electromagnetic fields.
- At this point the book shifts to a more qualitative (but non-trivial) treatment of atoms, atomic nuclei, the standard model of elementary particles, and techniques for observing the very small. Ideas from optics, waves, and relativity reappear here. The Bohr model of the hydrogen atom is not presented for the simple reason that it gets the angular momentum of the electron wrong!
- The final section of the course deals with heat and statistical mechanics. Only at this point do non-conservative forces appear in the context of classical mechanics. Counting as a way to compute the entropy is introduced, and is applied to the Einstein model of a collection of harmonic oscillators (conceptualized as a “brick”), and in a limited way to an ideal gas. The second law of thermodynamics follows. The book ends with a fairly conventional treatment of heat engines.

A few words about how I have taught the course at New Mexico Tech are in order. As with our standard course, each week contains three lecture hours and a two-hour recitation. The recitation is the key to making the course accessible to the students. I generally have small groups of students working on assigned homework problems during recitation while I wander around giving hints. After all groups have completed their work, a representative from each group explains their problem to the class. The students are then required to write up the problems on their own and hand them in at a later date. The problems are the key to student learning, and associating course credit with the successful solution of these problems insures virtually 100% attendance in recitation.

In addition, chapter reading summaries are required, with the students urged to ask questions about material in the text that gave them difficulties. Significant lecture time is taken up answering these questions. Students tend to do the summaries, as they also count for their grade. The summaries and the questions posed by the students have been quite helpful to me in indicating parts of the text which need clarification.

The writing style of the text is quite terse. This partially reflects its origin in a set of lecture notes, but it also focuses the students’ attention on what is really important. Given this structure, a knowledgeable instructor able to offer one-on-one time with students (as in our recitation sections) is essential for student success. The text is most likely to be useful in a sophomore-level course introducing physics majors to the broad world of physics viewed from a modern perspective.

I freely acknowledge stealing ideas from Edwin Taylor, John Archibald Wheeler, Thomas Moore, Robert Mills, Bruce Sherwood, and many other creative physicists, and I owe a great debt to them. The physics department at New Mexico Tech has been quite supportive of my efforts over the years relative to this course, for which I am exceedingly grateful. Finally, my humble thanks go out to the students who have enthusiastically (or on occasion unenthusiastically) responded to this course. It is much, much better as a result of their input.

My colleagues Alan Blyth, David Westpfahl, Ken Eack, and Sharon Sessions were brave enough to teach this course at various stages of its development, and I welcome the feedback I have received from them. Their experience shows that even seasoned physics teachers require time and effort to come to grips with the content of this textbook!

The reviews of Allan Stavely and Paul Arendt in conjunction with the publication of this book by the New Mexico Tech Press have been enormously helpful, and I am very thankful for their support and enthusiasm. Penny Bencomo and Keegan Livoti taught me a lot about written English with their copy editing.

David J. Raymond

Aside from numerous corrections, clarifications, and minor enhancements, the main additions to this edition include the following:

- The treatment of the harmonic oscillator is expanded, with complex exponential solutions, damped and forced oscillators, and several new problems.
- The discussion of Faraday’s law is enhanced, with a distinction being made between EMFs generated by changing magnetic fluxes through fixed loops and those generated by magnetic forces on charge moving through translating or deforming loops.
- A section on angular momentum and magnetic moments is added to the chapter on atoms. The Stern-Gerlach experiment is treated in a problem.

As in the first edition, I am greatful for the reviews of Paul Arendt and Allan Stavely, who always manage to catch things that I have overlooked.

David J. Raymond

Physics Department

New Mexico Tech

Socorro, NM, USA

djraymondnm@gmail.com

Waves in One Dimension

The wave is a universal phenomenon which occurs in a multitude of physical contexts. The purpose of this section is to describe the kinematics of waves, i. e., to provide tools for describing the form and motion of all waves irrespective of their underlying physical mechanisms.

Many examples of waves are well known to you. You undoubtedly know about ocean waves and have probably played with a stretched slinky toy, producing undulations which move rapidly along the slinky. Other examples of waves are sound, vibrations in solids, and light.

In this chapter we learn first about the basic properties of waves and introduce a special type of wave called the sine wave. Examples of waves seen in the real world are presented. We then learn about the superposition principle, which allows us to construct complex wave patterns by superimposing sine waves. Using these ideas, we discuss the related ideas of beats and interferometry. Finally, the ideas of wave packets and group velocity are introduced.

With the exception of light, waves are undulations in a material medium. For instance, ocean waves are (nearly) vertical undulations in the position of water parcels. The oscillations in neighboring parcels are phased such that a pattern moves across the ocean surface. Waves on a slinky are either transverse, in that the motion of the material of the slinky is perpendicular to the orientation of the slinky, or they are longitudinal, with material motion in the direction of the stretched slinky. (See figure 1.1.) Some media support only longitudinal waves, others support only transverse waves, while yet others support both types. Light waves are purely transverse, while sound waves are purely longitudinal. Ocean waves are a peculiar mixture of transverse and longitudinal, with parcels of water moving in elliptical trajectories as waves pass.

Light is a form of electromagnetic radiation. The undulations in an electromagnetic wave occur in the electric and magnetic fields. These oscillations are perpendicular to the direction of motion of the wave (in a vacuum), which is why we call light a transverse wave.

A particularly simple kind of wave, the sine wave, is illustrated in figure 1.2. This has the mathematical form

| (1.1) |

where h is the displacement (which can be either longitudinal or transverse), h_{0} is
the maximum displacement, also called the amplitude of the wave, and λ is the
wavelength. The oscillatory behavior of the wave is assumed to carry on to infinity
in both positive and negative x directions. Notice that the wavelength is the
distance through which the sine function completes one full cycle. The crest and
the trough of a wave are the locations of the maximum and minimum
displacements, as seen in figure 1.2.

So far we have only considered a sine wave as it appears at a particular time. All interesting waves move with time. The movement of a sine wave to the right a distance d may be accounted for by replacing x in the above formula by x - d. If this movement occurs in time t, then the wave moves at velocity c = d∕t. Solving this for d and substituting yields a formula for the displacement of a sine wave as a function of both distance x and time t:

| (1.2) |

The time for a wave to move one wavelength is called the period of the wave: T = λ∕c. Thus, we can also write

| (1.3) |

Physicists actually like to write the equation for a sine wave in a slightly simpler form. Defining the wavenumber as k = 2π∕λ and the angular frequency as ω = 2π∕T, we write

| (1.4) |

We normally think of the frequency of oscillatory motion as the number of cycles
completed per second. This is called the rotational frequency, and is given by
f = 1∕T. It is related to the angular frequency by ω = 2πf. The rotational
frequency is usually easier to measure than the angular frequency, but the angular
frequency tends to be used more often in theoretical discussions. As shown above,
converting between the two is not difficult. Rotational frequency is measured in
units of hertz, abbreviated Hz; 1 Hz = 1 cycle s^{-1}. Angular frequency also has the
dimensions of inverse time, e. g., radian s^{-1}, but the term “hertz” is generally
reserved only for rotational frequency.

The argument of the sine function is by definition an angle. We refer to this angle as the phase of the wave, ϕ = kx - ωt. The difference in the phase of a wave at fixed time over a distance of one wavelength is 2π, as is the difference in phase at fixed position over a time interval of one wave period.

Since angles are dimensionless, we normally don’t include this in the units for frequency. However, it sometimes clarifies things to refer to the dimensions of rotational frequency as “rotations per second” or angular frequency as “radians per second”.

As previously noted, we call h_{0}, the maximum displacement of the wave, the
amplitude. Often we are interested in the intensity of a wave, which is
proportional to the square of the amplitude. The intensity is often related to the
amount of energy being carried by a wave.

The wave speed we have defined above, c = λ∕T, is actually called the phase speed. Since λ = 2π∕k and T = 2π∕ω, we can write the phase speed in terms of the angular frequency and the wavenumber:

| (1.5) |

In order to make the above material more concrete, we now examine the characteristics of various types of waves which may be observed in the real world.

These waves are manifested as undulations of the ocean surface as seen in figure 1.3. The speed of ocean waves is given by the formula

| (1.6) |

where g = 9.8 m s^{-2} is the earth’s gravitational force per unit mass,
H is the depth of the ocean, and the hyperbolic tangent is defined
as^{1}

| (1.7) |

The equation for the speed of ocean waves comes from the theory for oscillations of a fluid surface in a gravitional field.

As figure 1.4 shows, for |x|≪ 1, we can approximate the hyperbolic tangent by tanh(x) ≈ x, while for |x|≫ 1 it is +1 for x > 0 and -1 for x < 0. This leads to two limits: Since x = kH, the shallow water limit, which occurs when kH ≪ 1, yields a wave speed of

| (1.8) |

while the deep water limit, which occurs when kH ≫ 1, yields

| (1.9) |

Notice that the speed of shallow water waves depends only on the depth of the water and on g. In other words, all shallow water waves move at the same speed. On the other hand, deep water waves of longer wavelength (and hence smaller wavenumber) move more rapidly than those with shorter wavelength. Waves for which the wave speed varies with wavelength are called dispersive. Thus, deep water waves are dispersive, while shallow water waves are non-dispersive.

For water waves with wavelengths of a few centimeters or less, surface tension becomes important to the dynamics of the waves. In the deep water case, the wave speed at short wavelengths is given by the formula

| (1.10) |

where the constant A is related to an effect called surface tension. For an
air-water interface near room temperature, A ≈ 74 cm^{3} s^{-2}.

Sound is a longitudinal compression-expansion wave in a fluid. The wave speed for sound in an ideal gas is

| (1.11) |

where γ and R are constants and T_{abs} is the absolute temperature. The
absolute temperature is measured in Kelvins and is numerically given
by

| (1.12) |

where T_{C} is the temperature in Celsius degrees. The angular frequency of sound
waves is thus given by

| (1.13) |

The speed of sound in air at normal temperatures is about 340 m s^{-1}.

Light moves in a vacuum at a speed of c_{vac} = 3 × 10^{8} m s^{-1}. In transparent
materials it moves at a speed less than c_{vac} by a factor n which is called the
refractive index of the material:

| (1.14) |

Often the refractive index takes the form

| (1.15) |

where k is the wavenumber and k_{R} and A are positive constants characteristic of
the material. The angular frequency of light in a transparent medium is
thus

| (1.16) |

It is found empirically that as long as the amplitudes of waves in most media are small, two waves in the same physical location don’t interact with each other. Thus, for example, two waves moving in the opposite direction simply pass through each other without their shapes or amplitudes being changed. When collocated, the total wave displacement is just the sum of the displacements of the individual waves. This is called the superposition principle. At sufficiently large amplitude the superposition principle often breaks down — interacting waves may scatter off of each other, lose amplitude, or change their form.

Interference is a consequence of the superposition principle. When two or more waves are superimposed, the net wave displacement is just the algebraic sum of the displacements of the individual waves. Since these displacements can be positive or negative, the net displacement can either be greater or less than the individual wave displacements. The former case, which occurs when both displacements are of the same sign, is called constructive interference, while destructive interference occurs when they are of opposite sign.

Let us see what happens when we superimpose two sine waves with
different wavenumbers. Figure 1.5 shows the superposition of two waves with
wavenumbers k_{1} = 4 and k_{2} = 5. Notice that the result is a wave with
about the same wavelength as the two initial waves, but which varies
in amplitude depending on whether the two sine waves are interfering
constructively or destructively. We say that the waves are in phase if they are
interfering constructively, and they are out of phase if they are interfering
destructively.

What happens when the wavenumbers of the two sine waves are changed?
Figure 1.6 shows the result when k_{1} = 10 and k_{2} = 11. Notice that though
the wavelength of the resultant wave is decreased, the locations where
the amplitude is maximum have the same separation in x as in figure
1.5.

If we superimpose waves with k_{1} = 10 and k_{2} = 12, as is shown in figure 1.7,
we see that the x spacing of the regions of maximum amplitude has decreased by
a factor of two. Thus, while the wavenumber of the resultant wave seems to be
related to something like the average of the wavenumbers of the component waves,
the spacing between regions of maximum wave amplitude appears to go inversely
with the difference of the wavenumbers of the component waves. In other words, if
k_{1} and k_{2} are close together, the amplitude maxima are far apart and vice
versa.

We can symbolically represent the sine waves that make up figures 1.5, 1.6, and 1.7 by a plot such as that shown in figure 1.8. The amplitudes and wavenumbers of each of the sine waves are indicated by vertical lines in this figure.

The regions of large wave amplitude are called wave packets. Wave packets will
play a central role in what is to follow, so it is important that we acquire a good
understanding of them. The wave packets produced by only two sine waves are
not well separated along the x-axis. However, if we superimpose many waves,
we can produce an isolated wave packet. For example, figure 1.9 shows
the results of superimposing 20 sine waves with wavenumbers k = 0.4m,
m = 1, 2,…, 20, where the amplitudes of the waves are largest for wavenumbers
near k = 4. In particular, we assume that the amplitude of each sine
wave is proportional to exp[-(k - k_{0})^{2}∕Δk^{2}], where k_{
0} = 4 defines the
maximum of the distribution of wavenumbers and Δk = 1 defines the
half-width of this distribution. The amplitudes of each of the sine waves
making up the wave packet in figure 1.9 are shown schematically in figure
1.10.

The quantity Δk controls the distribution of the sine waves being
superimposed — only those waves with a wavenumber k within approximately Δk
of the central wavenumber k_{0} of the wave packet, i. e., for 3 ≤ k ≤ 5 in
this case, contribute significantly to the sum. If Δk is changed to 2, so
that wavenumbers in the range 2 ≤ k ≤ 6 contribute significantly, the
wavepacket becomes narrower, as is shown in figures 1.11 and 1.12. Δk is
called the wavenumber spread of the wave packet, and it evidently plays
a role similar to the difference in wavenumbers in the superposition of
two sine waves — the larger the wavenumber spread, the smaller the
physical size of the wave packet. Furthermore, the wavenumber of the
oscillations within the wave packet is given approximately by the central
wavenumber.

We can better understand how wave packets work by mathematically
analyzing the simple case of the superposition of two sine waves. Let us define
k_{0} = (k_{1} + k_{2})∕2 where k_{1} and k_{2} are the wavenumbers of the component waves.
Furthermore let us set Δk = (k_{2} - k_{1})∕2. The quantities k_{0} and Δk are
graphically illustrated in figure 1.8. We can write k_{1} = k_{0} - Δk and k_{2} = k_{0} + Δk
and use the trigonometric identity sin(a + b) = sin(a) cos(b) + cos(a) sin(b) to find

| (1.18) |

Thus, as we observed in the earlier examples, the length of the wave packet Δx is inversely related to the spread of the wavenumbers Δk (which in this case is just the difference between the two wavenumbers) of the component waves. This relationship is central to the uncertainty principle of quantum mechanics.

Suppose two sound waves of different frequency but equal amplitude impinge on your ear at the same time. The displacement perceived by your ear is the superposition of these two waves, with time dependence

| (1.19) |

where we have used the above math trick, and where ω_{0} = (ω_{1} + ω_{2})∕2 and
Δω = (ω_{2} - ω_{1})∕2. What you actually hear is a tone with angular frequency ω_{0}
which fades in and out with period

| (1.20) |

| (1.21) |

Note how beats are the time analog of wave packets — the mathematics are the same except that frequency replaces wavenumber and time replaces space.

An interferometer is a device which splits a beam of light (or other wave) into two sub-beams, shifts the phase of one sub-beam with respect to the other, and then superimposes the sub-beams so that they interfere constructively or destructively, depending on the magnitude of the phase shift between them. In this section we study the Michelson interferometer and interferometric effects in thin films.

The American physicist Albert Michelson invented the optical interferometer illustrated in figure 1.13. The incoming beam is split into two beams by the half-silvered mirror. Each sub-beam reflects off of another mirror which returns it to the half-silvered mirror, where the two sub-beams recombine as shown. One of the reflecting mirrors is movable by a sensitive micrometer device, allowing the path length of the corresponding sub-beam, and hence the phase relationship between the two sub-beams, to be altered. As figure 1.13 shows, the difference in path length between the two sub-beams is 2x because the horizontal sub-beam traverses the path twice. Thus, constructive interference occurs when this path difference is an integral number of wavelengths, i. e.,

| (1.22) |

where λ is the wavelength of the wave and m is an integer. Note that m is the number of wavelengths that fits evenly into the distance 2x.

One of the most revealing examples of interference occurs when light interacts with a thin film of transparent material such as a soap bubble. Figure 1.14 shows how a plane wave normally incident on the film is partially reflected by the front and rear surfaces. The waves reflected off the front and rear surfaces of the film interfere with each other. The interference can be either constructive or destructive depending on the phase difference between the two reflected waves.

If the wavelength of the incoming wave is λ, one would naively expect constructive interference to occur between the A and B beams if 2d were an integral multiple of λ.

Two factors complicate this picture. First, the wavelength inside the film is not λ, but λ∕n, where n is the index of refraction of the film. Constructive interference would then occur if 2d = mλ∕n. Second, it turns out that an additional phase shift of half a wavelength occurs upon reflection when the wave is incident on material with a higher index of refraction than the medium in which the incident beam is immersed. This phase shift doesn’t occur when light is reflected from a region with lower index of refraction than felt by the incident beam. Thus beam B doesn’t acquire any additional phase shift upon reflection. As a consequence, constructive interference actually occurs when

| (1.23) |

while destructive interference results when

| (1.24) |

When we look at a soap bubble, we see bands of colors reflected back from a
light source. What is the origin of these bands? Light from ordinary sources is
generally a mixture of wavelengths ranging from roughly λ = 4.5 × 10^{-7} m (violet
light) to λ = 6.5 × 10^{-7} m (red light). In between violet and red we also
have blue, green, and yellow light, in that order. Because of the different
wavelengths associated with different colors, it is clear that for a mixed light
source we will have some colors interfering constructively while others
interfere destructively. Those undergoing constructive interference will be
visible in reflection, while those undergoing destructive interference will
not.

Another factor enters as well. If the light is not normally incident on the film, the difference in the distances traveled between beams reflected off of the front and rear faces of the film will not be just twice the thickness of the film. To understand this case quantitatively, we need the concept of refraction, which will be developed later in the context of geometrical optics. However, it should be clear that different wavelengths will undergo constructive interference for different angles of incidence of the incoming light. Different portions of the thin film will in general be viewed at different angles, and will therefore exhibit different colors under reflection, resulting in the colorful patterns normally seen in soap bubbles.

This section provides a quick review of the idea of the derivative. Often we are interested in the slope of a line tangent to a function y(x) at some value of x. This slope is called the derivative and is denoted dy∕dx. Since a tangent line to the function can be defined at any point x, the derivative itself is a function of x:

| (1.25) |

As figure 1.15 illustrates, the slope of the tangent line at some point on the function may be approximated by the slope of a line connecting two points, A and B, set a finite distance apart on the curve:

| (1.26) |

As B is moved closer to A, the approximation becomes better. In the limit when B moves infinitely close to A, it is exact.

Derivatives of some common functions are now given. In each case a is a constant.

| (1.27) |

| (1.28) |

| (1.29) |

| (1.30) |

| (1.31) |

| (1.32) |

| (1.33) |

| (1.34) |

| (1.35) |

The product and chain rules are used to compute the derivatives of complex functions. For instance,

and

We now ask the following question: How fast do wave packets move? Surprisingly, we often find that wave packets move at a speed very different from the phase speed, c = ω∕k, of the wave composing the wave packet.

We shall find that the speed of motion of wave packets, referred to as the group velocity, is given by

| (1.36) |

The derivative of ω(k) with respect to k is first computed and then evaluated at
k = k_{0}, the central wavenumber of the wave packet of interest.

The relationship between the angular frequency and the wavenumber for a wave, ω = ω(k), depends on the type of wave being considered. Whatever this relationship turns out to be in a particular case, it is called the dispersion relation for the type of wave in question.

As an example of a group velocity calculation, suppose we want to find the
velocity of deep ocean wave packets for a central wavelength of λ_{0} = 60 m. This
corresponds to a central wavenumber of k_{0} = 2π∕λ_{0} ≈ 0.1 m^{-1}. The phase
speed of deep ocean waves is c = (g∕k)^{1∕2}. However, since c ≡ ω∕k, we
find the frequency of deep ocean waves to be ω = (gk)^{1∕2}. The group
velocity is therefore u ≡ dω∕dk = (g∕k)^{1∕2}∕2 = c∕2. For the specified central
wavenumber, we find that u ≈ (9.8 m s^{-2}∕0.1 m^{-1})^{1∕2}∕2 ≈ 5 m s^{-1}. By
contrast, the phase speed of deep ocean waves with this wavelength is
c ≈ 10 m s^{-1}.

Dispersive waves are waves in which the phase speed varies with wavenumber. It is easy to show that dispersive waves have unequal phase and group velocities, while these velocities are equal for non-dispersive waves.

We now derive equation (1.36). It is easiest to do this for the simplest wave packets, namely those constructed out of the superposition of just two sine waves. We will proceed by adding two waves with full space and time dependence:

| (1.37) |

After algebraic and trigonometric manipulations familiar from earlier sections, we find

| (1.38) |

where as before we have k_{0} = (k_{1} + k_{2})∕2, ω_{0} = (ω_{1} + ω_{2})∕2, Δk = (k_{2} - k_{1})∕2,
and Δω = (ω_{2} - ω_{1})∕2.

Again think of this as a sine wave of frequency ω_{0} and wavenumber
k_{0} modulated by a cosine function. In this case the modulation pattern
moves with a speed so as to keep the argument of the cosine function
constant:

| (1.39) |

Differentiating this with respect to t while holding Δk and Δω constant yields

| (1.40) |

In the limit in which the deltas become very small, this reduces to the derivative

| (1.41) |

which is the desired result.

We now illustrate some examples of phase speed and group velocity by showing the displacement resulting from the superposition of two sine waves, as given by equation (1.38), in the x-t plane. This is an example of a spacetime diagram, of which we will see many examples later on.

The upper panel of figure 1.16 shows a non-dispersive case in which the phase speed equals the group velocity. The white and black regions indicate respectively strong wave crests and troughs (i. e., regions of large positive and negative displacements), with grays indicating a displacement near zero. Regions with large displacements indicate the location of wave packets. The positions of waves and wave packets at any given time may therefore be determined by drawing a horizontal line across the graph at the desired time and examining the variations in wave displacement along this line. The lower panel of this figure shows the wave displacement as a function of x at time t = 0 as an aid to interpretation of the upper panel.

Notice that as time increases, the crests move to the right. This corresponds to the motion of the waves within the wave packets. Note also that the wave packets, i. e., the broad regions of large positive and negative amplitudes, move to the right with increasing time as well.

Since velocity is distance moved Δx divided by elapsed time Δt, the slope of a line in figure 1.16, Δt∕Δx, is one over the velocity of whatever that line represents. The slopes of lines representing crests are the same as the slopes of lines representing wave packets in this case, which indicates that the two move at the same velocity. Since the speed of movement of wave crests is the phase speed and the speed of movement of wave packets is the group velocity, the two velocities are equal and the non-dispersive nature of this case is confirmed.

Figure 1.17 shows a dispersive wave in which the group velocity is twice the phase speed, while figure 1.18 shows a case in which the group velocity is actually opposite in sign to the phase speed. See if you can confirm that the phase and group velocities seen in each figure correspond to the values for these quantities calculated from the specified frequencies and wavenumbers.

- Measure your pulse rate. Compute the ordinary frequency of your heart beat in cycles per second. Compute the angular frequency in radians per second. Compute the period.
- An important wavelength for radio waves in radio astronomy is 21 cm.
(This comes from neutral hydrogen.) Compute the wavenumber of this
wave. Compute the ordinary and angular frequencies. (The speed of
light is 3 × 10
^{8}m s^{-1}.) - Sketch the resultant wave obtained from superimposing the waves A = sin(2x) and B = sin(3x). By using the trigonometric identity given in equation (1.17), obtain a formula for A+B in terms of sin(5x∕2) and cos(x∕2). Does the wave obtained from sketching this formula agree with your earlier sketch?
- Two sine waves with wavelengths λ
_{1}and λ_{2}are superimposed, making wave packets of length L. If we wish to make L larger, should we make λ_{1}and λ_{2}closer together or farther apart? Explain your reasoning. - By examining figure 1.9 versus figure 1.10 and then figure 1.11 versus figure 1.12, determine whether equation (1.18) works at least in an approximate sense for isolated wave packets.
- The frequencies of the chromatic scale in music are given by
(1.42) where f

_{0}is a constant equal to the frequency of the lowest note in the scale.- Compute f
_{1}through f_{11}if f_{0}= 440 Hz (the “A” note). - Using the above results, what is the beat frequency between the “A” (i = 0) and “B” (i = 2) notes? (The frequencies are given here in cycles per second rather than radians per second.)
- Which pair of the above frequencies f
_{0}- f_{11}yields the smallest beat frequency? Explain your reasoning.

- Compute f
- Large ships in general cannot move faster than the phase speed of surface
waves with a wavelength equal to twice the ship’s length. This is because
most of the propulsive force goes into making big waves under these
conditions rather than accelerating the ship.
- How fast can a 300 m long ship move in very deep water?
- As the ship moves into shallow water, does its maximum speed increase or decrease? Explain.

- Given the formula for refractive index of light quoted in this section, for what range of k does the phase speed of light in a transparent material take on real values which exceed the speed of light in a vacuum?
- A police radar works by splitting a beam of microwaves, part of which is
reflected back to the radar from your car where it is made to interfere
with the other part which travels a fixed path, as shown in figure
1.19.
- If the wavelength of the microwaves is λ, how far do you have to travel in your car for the interference between the two beams to go from constructive to destructive to constructive?
- If you are traveling toward the radar at speed v = 30 m s
^{-1}, use the above result to determine the number of times per second constructive interference peaks will occur. Assume that λ = 3 cm.

- Suppose you know the wavelength of light passing through a Michelson interferometer with high accuracy. Describe how you could use the interferometer to measure the length of a small piece of material.
- A Fabry-Perot interferometer (see figure 1.20) consists of two parallel half-silvered mirrors placed a distance d from each other as shown. The beam passing straight through interferes with the beam which reflects once off of both of the mirrored surfaces as shown. For wavelength λ, what values of d result in constructive interference?
- A Fabry-Perot interferometer has spacing d = 2 cm between the glass
plates, causing the direct and doubly reflected beams to interfere (see
figure 1.20). As air is pumped out of the gap between the plates, the
beams go through 23 cycles of constructive-destructive-constructive
interference. If the wavelength of the light in the interfering beams is
5 × 10
^{-7}m, determine the index of refraction of the air initially in the interferometer. - Measurements on a certain kind of wave reveal that the angular
frequency of the wave varies with wavenumber as shown in the following
table:
ω (s ^{-1})k (m ^{-1})5 1 20 2 45 3 80 4 125 5 - Compute the phase speed of the wave for k = 3 m
^{-1}and for k = 4 m^{-1}. - Estimate the group velocity for k = 3.5 m
^{-1}using a finite difference approximation to the derivative.

- Compute the phase speed of the wave for k = 3 m
- Suppose some type of wave has the (admittedly weird) dispersion relation
shown in figure 1.21.
- For what values of k is the phase speed of the wave positive?
- For what values of k is the group velocity positive?

- Group velocities of various waves.
- Compute the group velocity for shallow water waves. Compare it with the phase speed of shallow water waves. (Hint: You first need to derive a formula for ω(k) from c(k).)
- Repeat the above problem for deep water waves.
- Repeat for sound waves. What does this case have in common with shallow water waves?

Waves in Two and Three Dimensions

In this chapter we extend the ideas of the previous chapter to the case of waves in more than one dimension. The extension of the sine wave to higher dimensions is the plane wave. Wave packets in two and three dimensions arise when plane waves moving in different directions are superimposed.

Diffraction results from the disruption of a wave which is impingent upon an object. Those parts of the wave front hitting the object are scattered, modified, or destroyed. The resulting diffraction pattern comes from the subsequent interference of the various pieces of the modified wave. A knowledge of diffraction is necessary to understand the behavior and limitations of optical instruments such as telescopes.

Diffraction and interference in two and three dimensions can be manipulated to produce useful devices such as the diffraction grating.

Before we can proceed further we need to explore the idea of a vector. A vector is a quantity which expresses both magnitude and direction. Graphically we represent a vector as an arrow. In typeset notation a vector is represented by a boldface character, while in handwriting an arrow is drawn over the character representing the vector.

Figure 2.1 shows some examples of displacement vectors, i. e., vectors which represent the displacement of one object from another, and introduces the idea of vector addition. The tail of vector B is collocated with the head of vector A, and the vector which stretches from the tail of A to the head of B is the sum of A and B, called C in figure 2.1.

The quantities A_{x}, A_{y}, etc., represent the Cartesian components of the vectors
in figure 2.1. A vector can be represented either by its Cartesian components,
which are just the projections of the vector onto the Cartesian coordinate axes, or
by its direction and magnitude. The direction of a vector in two dimensions is
generally represented by the counterclockwise angle of the vector relative to the x
axis, as shown in figure 2.2. Conversion from one form to the other is given by the
equations

| (2.1) |

| (2.2) |

where A is the magnitude of the vector. A vector magnitude is sometimes represented by absolute value notation: A ≡|A|.

Notice that the inverse tangent gives a result which is ambiguous relative to
adding or subtracting integer multiples of π. Thus the quadrant in which the
angle lies must be resolved by independently examining the signs of A_{x} and A_{y}
and choosing the appropriate value of θ.

To add two vectors, A and B, it is easiest to convert them to Cartesian component form. The components of the sum C = A + B are then just the sums of the components:

| (2.3) |

Subtraction of vectors is done similarly, e. g., if A = C - B, then

| (2.4) |

A unit vector is a vector of unit length. One can always construct a unit vector from an ordinary (non-zero) vector by dividing the vector by its length: n = A∕|A|. This division operation is carried out by dividing each of the vector components by the number in the denominator. Alternatively, if the vector is expressed in terms of length and direction, the magnitude of the vector is divided by the denominator and the direction is unchanged.

Unit vectors can be used to define a Cartesian coordinate system.
Conventionally, i, j, and k indicate the x, y, and z axes of such a system. Note
that i, j, and k are mutually perpendicular. Any vector can be represented in
terms of unit vectors and its Cartesian components: A = A_{x}i + A_{y}j + A_{z}k. An
alternate way to represent a vector is as a list of components: A = (A_{x},A_{y},A_{z}).
We tend to use the latter representation since it is somewhat more economical
notation.

There are two ways to multiply two vectors, yielding respectively what are known as the dot product and the cross product. The cross product yields another vector while the dot product yields a number. Here we will discuss only the dot product. The cross product will be presented later when it is needed.

Given vectors A and B, the dot product of the two is defined as

| (2.5) |

where θ is the angle between the two vectors. In two dimensions an alternate expression for the dot product exists in terms of the Cartesian components of the vectors:

| (2.6) |

It is easy to show that this is equivalent to the cosine form of the dot
product when the x axis lies along one of the vectors, as in figure 2.3.
Notice in particular that A_{x} = |A| cos θ, while B_{x} = |B| and B_{y} = 0. Thus,
A ⋅ B = |A| cos θ|B| in this case, which is identical to the form given in equation
(2.5).

All that remains to be proven for equation (2.6) to hold in general is
to show that it yields the same answer regardless of how the Cartesian
coordinate system is oriented relative to the vectors. To do this, we must
show that A_{x}B_{x} + A_{y}B_{y} = A_{x}′B_{x}′ + A_{y}′B_{y}′, where the primes indicate
components in a coordinate system rotated from the original coordinate
system.

Figure 2.4 shows the vector R resolved in two coordinate systems rotated with
respect to each other. From this figure it is clear that X′ = a + b. Focusing on the
shaded triangles, we see that a = X cos θ and b = Y sin θ. Thus, we find
X′ = X cos θ + Y sin θ. Similar reasoning shows that Y ′ = -X sin θ + Y cos θ.
Substituting these and using the trigonometric identity cos ^{2}θ + sin ^{2}θ = 1 results
in

A numerical quantity that doesn’t depend on which coordinate system is being used is called a scalar. The dot product of two vectors is a scalar. However, the components of a vector, taken individually, are not scalars, since the components change as the coordinate system changes. Since the laws of physics cannot depend on the choice of coordinate system being used, we insist that physical laws be expressed in terms of scalars and vectors, but not in terms of the components of vectors.

In three dimensions the cosine form of the dot product remains the same, while the component form is

| (2.8) |

A plane wave in two or three dimensions is like a sine wave in one dimension except that crests and troughs aren’t points, but form lines (2-D) or planes (3-D) perpendicular to the direction of wave propagation. Figure 2.5 shows a plane sine wave in two dimensions. The large arrow is a vector called the wave vector, which defines (1) the direction of wave propagation by its orientation perpendicular to the wave fronts, and (2) the wavenumber by its length. We can think of a wave front as a line along the crest of the wave. The equation for the displacement associated with a plane sine wave (of unit amplitude) in three dimensions at some instant in time is

| (2.9) |

Since wave fronts are lines or surfaces of constant phase, the equation defining a wave front is simply k ⋅ x = const.

In the two dimensional case we simply set k_{z} = 0. Therefore, a wave
front, or line of constant phase ϕ in two dimensions is defined by the
equation

| (2.10) |

This can be easily solved for y to obtain the slope and intercept of the wave front in two dimensions.

As for one dimensional waves, the time evolution of the wave is obtained by adding a term -ωt to the phase of the wave. In three dimensions the wave displacement as a function of both space and time is given by

| (2.11) |

The frequency depends in general on all three components of the wave vector. The
form of this function, ω = ω(k_{x},k_{y},k_{z}), which as in the one dimensional case is
called the dispersion relation, contains information about the physical behavior of
the wave.

Some examples of dispersion relations for waves in two dimensions are as follows:

- Light waves in a vacuum in two dimensions obey
(2.12) where c is the speed of light in a vacuum.

- Deep water ocean waves in two dimensions obey
(2.13) where g is the strength of the Earth’s gravitational field as before.

- Certain kinds of atmospheric waves confined to a vertical x-z plane called
gravity waves (not to be confused with the gravitational waves of general
relativity)
^{1}obey(2.14) where N is a constant with the dimensions of inverse time called the Brunt-Väisälä frequency.

Contour plots of these dispersion relations are plotted in the upper
panels of figure 2.6. These plots are to be interpreted like topographic
maps, where the lines represent contours of constant elevation. In the case
of figure 2.6, constant values of frequency are represented instead. For
simplicity, the actual values of frequency are not labeled on the contour plots,
but are represented in the graphs in the lower panels. This is possible
because frequency depends only on wave vector magnitude (k_{x}^{2} + k_{
y}^{2})^{1∕2}
for the first two examples, and only on wave vector direction θ for the
third.

We now study wave packets in two dimensions by asking what the superposition of two plane sine waves looks like. If the two waves have different wavenumbers, but their wave vectors point in the same direction, the results are identical to those presented in the previous chapter, except that the wave packets are indefinitely elongated without change in form in the direction perpendicular to the wave vector. The wave packets produced in this case move in the direction of the wave vectors and thus appear to a stationary observer like a series of passing pulses with broad lateral extent.

Superimposing two plane waves which have the same frequency results in a stationary wave packet through which the individual wave fronts pass. This wave packet is also elongated indefinitely in some direction, but the direction of elongation depends on the dispersion relation for the waves being considered. These wave packets are in the form of steady beams, which guide the individual phase waves in some direction, but don’t themselves change with time. By superimposing multiple plane waves, all with the same frequency, one can actually produce a single stationary beam, just as one can produce an isolated pulse by superimposing multiple waves with wave vectors pointing in the same direction.

If the frequency of a wave depends on the magnitude of the wave vector, but not on its direction, the wave’s dispersion relation is called isotropic; otherwise it is anisotropic. In the isotropic case, two waves have the same frequency only if the lengths of their wave vectors, and hence their wavelengths, are the same. The first two examples in figure 2.6 satisfy this condition, while the last example is anisotropic.

We now use the language of vectors to investigate the superposition of two
plane waves with wave vectors k_{1} and k_{2}:

| (2.15) |

Applying the trigonometric identity for the sine of the sum of two angles (as we have done previously), equation (2.15) can be reduced to

| (2.16) |

where

| (2.17) |

This is in the form of a sine wave moving in the k_{0} direction with phase speed
c_{phase} = ω∕|k_{0}| and wavenumber |k_{0}|, modulated in the Δk direction by a cosine
function. The lines of destructive interference are normal to Δk. The distance w
between lines of destructive interference is the distance between successive zeros of
the cosine function in equation (2.16), implying that |Δk|w = π, which leads
to

| (2.18) |

Thus, the smaller |Δk|, the greater is the beam diameter.

In this section we investigate the beams produced by superimposing isotropic
waves of the same frequency. Figure 2.7 illustrates what happens in such a
superposition. Vectors k_{1} and k_{2} of equal length give rise to a mean wave vector
k_{0} and half the difference, Δk. As illustrated, the lines of constructive and
destructive interference are perpendicular to Δk. Figure 2.8 shows a concrete
example of the beams produced by superposition of two plane waves of equal
wavelength oriented as in figure 2.7. The beams are aligned vertically, since Δk is
horizontal, with the lines of destructive interference separating the beams located
near x = ±16. The transverse width of the beams of ≈ 32 satisfies equation (2.18)
with |Δk| = 0.1. Each beam is made up of vertically propagating phase
waves, with the crests and troughs indicated by the regions of white and
black.

In the third example of figure 2.6, the frequency of the wave depends only on the direction of the wave vector, independent of its magnitude, which is the reverse of the case for an isotropic dispersion relation. In this highly anisotropic case, different plane waves with the same frequency have wave vectors which point in the same direction, but have different lengths.

More generally, one might have waves for which the frequency depends on both the direction and magnitude of the wave vector. In this case, two different plane waves with the same frequency would typically have wave vectors which differ both in direction and magnitude. Such an example is illustrated in figures 2.9 and 2.10.

Figure 2.11 summarizes what we have learned about adding plane waves with the same frequency. In general, the beam orientation (and the lines of constructive interference) are not perpendicular to the wave fronts. This only occurs when the wave frequency is independent of wave vector direction.

As with wave packets in one dimension, we can add together more than two waves to produce an isolated wave packet. We will confine our attention here to the case of an isotropic dispersion relation in which all the wave vectors for a given frequency are of the same length.

Figure 2.12 shows an example of this in which wave vectors of the same
wavelength but different directions are added together. Defining α_{i} as the angle of
the ith wave vector clockwise from the vertical, as illustrated in figure
2.12, we could write the superposition of these waves at time t = 0 as

| (2.20) |

The exponential function decreases rapidly as its argument becomes more
negative, and for practical purposes, only wave vectors with |α_{i}|≤ α_{max}
contribute significantly to the sum. We call α_{max} the spreading angle.

Figure 2.13 shows what h(x,y) looks like when α_{max} = 0.8 radians and k = 1.
Notice that for y = 0 the wave amplitude is only large for a small region in the
range -4 < x < 4. However, for y > 0 the wave spreads into a broad, semicircular
pattern.

Figure 2.14 shows the computed pattern of h(x,y) when the spreading angle
α_{max} = 0.2 radians. The wave amplitude is large for a much broader range of x
at y = 0 in this case, roughly -12 < x < 12. On the other hand, the
subsequent spread of the wave is much smaller than in the case of figure
2.13.

We conclude that a superposition of plane waves with wave vectors spread narrowly about a central wave vector which points in the y direction (as in figure 2.14) produces a beam which is initially broad in x but for which the breadth increases only slightly with increasing y. However, a superposition of plane waves with wave vectors spread more broadly (as in figure 2.13) produces a beam which is initially narrow in x but which rapidly increases in width as y increases.

The relationship between the spreading angle α_{max} and the initial breadth of
the beam is made more understandable by comparison with the results for the
two-wave superposition discussed at the beginning of this section. As
indicated by equation (2.18), large values of k_{x}, and hence α, are associated
with small wave packet dimensions in the x direction and vice versa.
The superposition of two waves doesn’t capture the subsequent spread
of the beam which occurs when many waves are superimposed, but it
does lead to a rough quantitative relationship between α_{max} (which is
just tan ^{-1}(k_{
x}∕k_{y}) in the two wave case) and the initial breadth of the
beam. If we invoke the small angle approximation for α = α_{max} so that
α_{max} = tan ^{-1}(k_{
x}∕k_{y}) ≈ k_{x}∕k_{y} ≈ k_{x}∕k, then k_{x} ≈ kα_{max} and equation (2.18) can
be written w = π∕k_{x} ≈ π∕(kα_{max}) = λ∕(2α_{max}). Thus, we can find the
approximate spreading angle from the wavelength of the wave λ and the initial
breadth of the beam w:

| (2.21) |

How does all of this apply to the passage of waves through a slit? Imagine a plane wave of wavelength λ impingent on a barrier with a slit. The barrier transforms the plane wave with infinite extent in the lateral direction into a beam with initial transverse dimensions equal to the width of the slit. The subsequent development of the beam is illustrated in figures 2.13 and 2.14, and schematically in figure 2.15. In particular, if the slit width is comparable to the wavelength, the beam spreads broadly as in figure 2.13. If the slit width is large compared to the wavelength, the beam doesn’t spread as much, as figure 2.14 illustrates. Equation (2.21) gives us an approximate quantitative result for the spreading angle if w is interpreted as the width of the slit.

One use of the above equation is in determining the maximum angular
resolution of optical instruments such as telescopes. The primary lens or mirror
can be thought of as a rather large “slit”. Light from a distant point source is
essentially in the form of a plane wave when it arrives at the telescope.
However, the light passed by the telescope is no longer a plane wave, but is a
beam with a tendency to spread. The spreading angle α_{max} is given by
equation (2.21), and the telescope cannot resolve objects with an angular
separation less than α_{max}. Replacing w with the diameter of the lens or
mirror in equation (2.21) thus yields the telescope’s angular resolution. For
instance, a moderate sized telescope with aperture 1 m observing red
light with λ ≈ 6 × 10^{-7} m has a maximum angular resolution of about
3 × 10^{-7} radians.

Let us now imagine a plane sine wave normally impingent on a screen with two narrow slits spaced by a distance d, as shown in figure 2.16. Since the slits are narrow relative to the wavelength of the wave impingent on them, the spreading angle of the beams is large and the diffraction pattern from each slit individually is a cylindrical wave spreading out in all directions, as illustrated in figure 2.13. The cylindrical waves from the two slits interfere, resulting in oscillations in wave intensity at the screen on the right side of figure 2.16.

Constructive interference occurs when the difference in the paths traveled by
the two waves from their originating slits to the screen, L_{2} - L_{1}, is an integer
multiple m of the wavelength λ: L_{2} - L_{1} = mλ. If L_{0} ≫ d, the lines L_{1}
and L_{2} are nearly parallel, which means that the narrow end of the dark
triangle in figure 2.16 has an opening angle of θ. Thus, the path difference
between the beams from the two slits is L_{2} - L_{1} = d sin θ. Substitution of
this into the above equation shows that constructive interference occurs
when

| (2.22) |

Destructive interference occurs when m is an integer plus 1∕2. The integer m is called the interference order and is the number of wavelengths by which the two paths differ.

Since the angular spacing Δθ of interference peaks in the two slit case depends on the wavelength of the incident wave, the two slit system can be used as a crude device to distinguish between the wavelengths of different components of a non-sinusoidal wave impingent on the slits. However, if more slits are added, maintaining a uniform spacing d between slits, we obtain a more sophisticated device for distinguishing beam components. This is called a diffraction grating.

Figures 2.17-2.19 show the intensity of the diffraction pattern as a function of position x on the display screen (see figure 2.16) for gratings with 2, 4, and 16 slits respectively, with the same slit spacing. Notice how the interference peaks remain in the same place but increase in sharpness as the number of slits increases.

The width of the peaks is actually related to the overall width of the grating, w = nd,
where n is the number of slits. Thinking of this width as the dimension of a large
single slit, the single slit equation, α_{max} = λ∕(2w), tells us the angular width of the
peaks.^{2}

Whereas the angular width of the interference peaks is governed by the single
slit equation, their angular positions are governed by the two slit equation. Let us
assume for simplicity that |θ|≪ 1 so that we can make the small angle
approximation to the two slit equation, mλ = d sin θ ≈ dθ, and ask the following
question: How different do two wavelengths differing by Δλ have to be in order
that the interference peaks from the two waves not overlap? In order for the
peaks to be distinguishable, they should be separated in θ by an angle
Δθ = mΔλ∕d, which is greater than the angular width of each peak,
α_{max}:

| (2.23) |

Substituting in the above expressions for Δθ and α_{max} and solving for Δλ, we get
Δλ > λ∕(2mn), where λ is the average of the two wavelengths and n = w∕d is the
number of slits in the diffraction grating. Thus, the fractional difference
between wavelengths which can be distinguished by a diffraction grating
depends solely on the interference order m and the number of slits n in the
grating:

| (2.24) |

- Point A is at the origin. Point B is 3 m distant from A at 30
^{∘}counterclockwise from the x axis. Point C is 2 m from point A at 100^{∘}counterclockwise from the x axis.- Obtain the Cartesian components of the vector D
_{1}which goes from A to B and the vector D_{2}which goes from A to C. - Find the Cartesian components of the vector D
_{3}which goes from B to C. - Find the direction and magnitude of D
_{3}.

- Obtain the Cartesian components of the vector D
- For the vectors in the previous problem, find D
_{1}⋅ D_{2}using both the cosine form of the dot product and the Cartesian form. Check to see if the two answers are the same. - Show graphically or otherwise that |A + B|≠|A| + |B| except when the vectors A and B are parallel.
- A wave in the x-y plane is defined by h = h
_{0}sin(k ⋅ x) where k = (1, 2) cm^{-1}.- On a piece of graph paper draw x and y axes and then plot a line passing through the origin which is parallel to the vector k.
- On the same graph plot the line defined by k⋅x = k
_{x}x + k_{y}y = 0, k ⋅ x = π, and k ⋅ x = 2π. Check to see if these lines are perpendicular to k.

- A plane wave in two dimensions in the x - y plane moves in the direction
45
^{∘}counterclockwise from the x-axis as shown in figure 2.20. Determine how fast the intersection between a wave front and the x-axis moves to the right in terms of the phase speed c of the wave. Hint: What is the distance between wave fronts along the x-axis compared to the wavelength? - Two deep plane ocean waves with the same frequency ω are moving
approximately to the east. However, one wave is oriented a small angle β
north of east and the other is oriented β south of east.
- Determine the orientation of lines of constructive interference between these two waves.
- Determine the spacing between lines of constructive interference.

- An example of waves with a dispersion relation in which the frequency is a
function of both wave vector magnitude and direction is shown graphically
in figure 2.21.
- What is the phase speed of the waves for each of the three wave vectors? Hint: You may wish to obtain the length of each wave vector graphically.
- For each of the wave vectors, what is the orientation of the wave fronts?
- For each of the illustrated wave vectors, sketch two other wave vectors whose average value is approximately the illustrated vector, and whose tips lie on the same frequency contour line. Determine the orientation of lines of constructive interference produced by the superimposing pairs of plane waves for which each of the vector pairs are the wave vectors.

- Two gravity waves have the same frequency, but slightly different
wavelengths.
- If one wave has an orientation angle θ = π∕4 radians, what is the orientation angle of the other? (See figure 2.6.)
- Determine the orientation of lines of constructive interference between these two waves.

- A plane wave impinges on a single slit, spreading out a half-angle α after the
slit. If the whole apparatus is submerged in a liquid with index of refraction
n = 1.5, how does the spreading angle of the light change? (Hint: Recall that
the index of refraction in a transparent medium is the ratio of the speed of
light in a vacuum to the speed in the medium. Furthermore, when light
goes from a vacuum to a transparent medium, the light frequency
doesn’t change. Therefore, how does the wavelength of the light
change?)

- Determine the diameter of the telescope needed to resolve a planet
2 × 10
^{8}km from a star which is 6 light years from the earth. (Assume blue light which has a wavelength λ ≈ 4 × 10^{-7}m = 400 nm. Also, don’t worry about the great difference in brightness between the two for the purposes of this problem.) - A laser beam from a laser on the earth is bounced back to the earth by a
corner reflector on the moon.
- Engineers find that the returned signal is stronger if the laser beam is initially spread out by the beam expander shown in figure 2.22. Explain why this is so.
- The beam has a diameter of 1 m leaving the earth. How broad is
it when it reaches the moon, which is 4 × 10
^{5}km away? Assume the wavelength of the light to be 5 × 10^{-7}m. - How broad would the laser beam be at the moon if it weren’t initially passed through the beam expander? Assume its initial diameter to be 1 cm.

- Suppose that a plane wave impinges on two slits in a barrier at an angle, such that the phase of the wave at one slit lags the phase at the other slit by half a wavelength. How does the resulting interference pattern change from the case in which there is no lag?
- Suppose that a thin piece of glass of index of refraction n = 1.33 is placed in
front of one slit of a two slit diffraction setup.
- How thick does the glass have to be to slow down the incoming
wave so that it lags the wave going through the other slit by
a phase difference of π? Take the wavelength of the light to be
λ = 6 × 10
^{-7}m. - For the above situation, describe qualitatively how the diffraction pattern changes from the case in which there is no glass in front of one of the slits. Explain your results.

- How thick does the glass have to be to slow down the incoming
wave so that it lags the wave going through the other slit by
a phase difference of π? Take the wavelength of the light to be
λ = 6 × 10
- A light source produces two wavelengths, λ
_{1}= 400 nm (blue) and λ_{2}= 600 nm (red).- Qualitatively sketch the two slit diffraction pattern from this source. Sketch the pattern for each wavelength separately.
- Qualitatively sketch the 16 slit diffraction pattern from this source, where the slit spacing is the same as in the two slit case.

- A light source produces two wavelengths, λ
_{1}= 631 nm and λ_{2}= 635 nm. What is the minimum number of slits needed in a grating spectrometer to resolve the two wavelengths? (Assume that you are looking at the first order diffraction peak.) Sketch the diffraction peak from each wavelength and indicate how narrow the peaks must be to resolve them.

Geometrical Optics

As was shown previously, when a plane wave is impingent on an aperture which has dimensions much greater than the wavelength of the wave, diffraction effects are minimal and a segment of the plane wave passes through the aperture essentially unaltered. This plane wave segment can be thought of as a wave packet, called a beam or ray, consisting of a superposition of wave vectors very close in direction and magnitude to the central wave vector of the wave packet. In most cases the ray simply moves in the direction defined by the central wave vector, i. e., normal to the orientation of the wave fronts. However, this is not true when the medium through which the light propagates is optically anisotropic, i. e., light traveling in different directions moves at different phase speeds. An example of such a medium is a calcite crystal. In the anisotropic case, the orientation of the ray can be determined once the dispersion relation for the waves in question is known, by using the techniques developed in the previous chapter.

If light moves through some apparatus in which all apertures are much greater in dimension than the wavelength of light, then we can use the above rule to follow rays of light through the apparatus. This is called the geometrical optics approximation.

Most of what we need to know about geometrical optics can be summarized in two rules, the laws of reflection and refraction. These rules may both be inferred by considering what happens when a plane wave segment impinges on a flat surface. If the surface is polished metal, the wave is reflected, whereas if the surface is an interface between two transparent media with differing indices of refraction, the wave is partially reflected and partially refracted. Reflection means that the wave is turned back into the half-space from which it came, while refraction means that it passes through the interface, acquiring a different direction of motion from that which it had before reaching the interface.

Figure 3.1 shows the wave vector and wave front of a wave being reflected from
a plane mirror. The angles of incidence, θ_{I}, and reflection, θ_{R}, are defined to be
the angles between the incoming and outgoing wave vectors respectively and the
line normal to the mirror. The law of reflection states that θ_{R} = θ_{I}. This is a
consequence of the need for the incoming and outgoing wave fronts to be in phase
with each other all along the mirror surface. This plus the equality of
the incoming and outgoing wavelengths is sufficient to insure the above
result.

Refraction, as illustrated in figure 3.2, is slightly more complicated. Since
n_{R} > n_{I}, the speed of light in the right-hand medium is less than in the left-hand
medium. (Recall that the speed of light in a medium with refractive index n is
c_{medium} = c_{vac}∕n.) The frequency of the wave packet doesn’t change as it passes
through the interface, so the wavelength of the light on the right side is less than
the wavelength on the left side.

Let us examine the triangle ABC in figure 3.2. The side AC is equal to the side
BC times sin(θ_{I}). However, AC is also equal to 2λ_{I}, or twice the wavelength of the
wave to the left of the interface. Similar reasoning shows that 2λ_{R}, twice the
wavelength to the right of the interface, equals BC times sin(θ_{R}). Since
the interval BC is common to triangles ABC and DBC, we easily see
that

| (3.1) |

Since λ_{I} = c_{I}T = c_{vac}T∕n_{I} and λ_{R} = c_{R}T = c_{vac}T∕n_{R} where c_{I} and c_{R} are the
wave speeds to the left and right of the interface, c_{vac} is the speed of light in a
vacuum, and T is the (common) period, we can easily recast the above equation in
the form

| (3.2) |

This is called Snell’s law, and it governs how a ray of light bends as it
passes through a discontinuity in the index of refraction. The angle θ_{I}
is called the incident angle and θ_{R} is called the refracted angle. Notice
that these angles are measured from the normal to the surface, not the
tangent.

When light passes from a medium of lesser index of refraction to one with greater
index of refraction, Snell’s law indicates that the ray bends toward the normal to
the interface. The reverse occurs when the passage is in the other direction.
In this latter circumstance a special situation arises when Snell’s law
predicts a value for the sine of the refracted angle greater than one. This is
physically untenable. What actually happens is that the incident wave
is reflected from the interface. This phenomenon is called total internal
reflection. The minimum incident angle for which total internal reflection
occurs is obtained by substituting θ_{R} = π∕2 into equation (3.2), resulting
in

| (3.3) |

Notice that Snell’s law makes the implicit assumption that rays of light move in the direction of the light’s wave vector, i. e., normal to the wave fronts. As the analysis in the previous chapter makes clear, this is valid only when the optical medium is isotropic, i. e., the wave frequency depends only on the magnitude of the wave vector, not on its direction.

Certain kinds of crystals, such as those made of calcite, are not isotropic — the speed of light in such crystals, and hence the wave frequency, depends on the orientation of the wave vector. As an example, the angular frequency in an anisotropic medium might take the form

| (3.4) |

where c_{1} is the speed of light for waves in which k_{y} = k_{x}, and c_{2} is its speed when
k_{y} = -k_{x}.

Figure 3.3 shows an example in which a ray hits a calcite crystal oriented so that constant frequency contours are as specified in equation (3.4). The wave vector is oriented normal to the surface of the crystal, so that wave fronts are parallel to this surface. Upon entering the crystal, the wave front orientation must stay the same to preserve phase continuity at the surface. However, due to the anisotropy of the dispersion relation for light in the crystal, the ray direction changes as shown in the right panel. This behavior is clearly inconsistent with the usual version of Snell’s law!

It is possible to extend Snell’s law to the anisotropic case. However, we will not present this here. The following discussions of optical instruments will always assume that isotropic optical media are used.

Given the laws of reflection and refraction, one can see in principle how the passage of light through an optical instrument could be traced. For each of a number of initial rays, the change in the direction of the ray at each mirror surface or refractive index interface can be calculated. Between these points, the ray traces out a straight line.

Though simple in conception, this procedure can be quite complex in practice. However, the procedure simplifies if a number of approximations, collectively called the thin lens approximation, are valid. We begin with the calculation of the bending of a ray of light as it passes through a prism, as illustrated in figure 3.4.

The pieces of information needed to find θ, the angle through which the ray is deflected, are as follows: the geometry of the triangle defined by the entry and exit points of the ray and the upper vertex of the prism leads to

| (3.5) |

which simplifies to

| (3.6) |

Snell’s law at the entrance and exit points of the ray tell us that

| (3.7) |

where n is the index of refraction of the prism. (The index of refraction of the surroundings is assumed to be unity.) One can also infer that

| (3.8) |

This comes from the fact that the the sum of the internal angles of the shaded quadrangle in figure 3.4 is

| (3.9) |

Combining equations (3.6), (3.7), and (3.8) allows the ray deflection θ to be
determined in terms of θ_{1} and α, but the resulting expression is very
messy. However, great simplification occurs if the following conditions are
met:

- The angle α ≪ 1.
- The angles θ
_{1}, θ_{2}, θ_{3}, θ_{4}≪ 1.

With these approximations it is easy to show that

| (3.10) |

Generally speaking, lenses and mirrors in optical instruments have curved rather than flat surfaces. However, we can still use the laws for reflection and refraction by plane surfaces as long as the segment of the surface on which the wave packet impinges is not curved very much on the scale of the wave packet dimensions. This condition is easy to satisfy with light impinging on ordinary optical instruments. In this case, the deflection of a ray of light is given by equation (3.10) if α is defined as the intersection of the tangent lines to the entry and exit points of the ray, as illustrated in figure 3.5.

A positive lens is thicker in the center than at the edges. The angle α between
the tangent lines to the two surfaces of the lens at a distance r from the
central axis takes the form α = Cr, where C is a constant. The deflection
angle of a beam hitting the lens a distance r from the center is therefore
θ = Cr(n - 1), as indicated in figure 3.5. The angles o and i sum to the
deflection angle: o + i = θ = Cr(n - 1). However, to the extent that the
small angle approximation holds, o = r∕d_{o} and i = r∕d_{i} where d_{o} is the
distance to the object and d_{i} is the distance to the image of the object.
Putting these equations together and cancelling the r results in the thin lens
equation:

| (3.11) |

The quantity f is called the focal length of the lens. Notice that f = d_{i} if the
object is very far from the lens, i. e., if d_{o} is extremely large.

Figure 3.6 shows how a positive lens makes an image. The image is produced
by all of the light from each point on the object falling on a corresponding point
in the image. If the arrow on the left is an illuminated object, an image of the
arrow will appear at right if the light coming from the lens is allowed to fall on a
piece of paper or a ground glass screen. The size of the object S_{o} and the size of
the image S_{i} are related by simple geometry to the distances of the object and the
image from the lens:

| (3.12) |

Notice that a positive lens inverts the image.

An image will be produced to the right of the lens only if d_{o} > f. If d_{o} < f, the
lens is unable to converge the rays from the image to a point, as is seen in figure
3.7. However, in this case the backward extension of the rays converge at a point
called a virtual image, which in the case of a positive lens is always farther away
from the lens than the object. The image is called virtual because it does
not appear on a ground glass screen placed at this point. Unlike the real
image seen in figure 3.6, the virtual image is not inverted. The thin lens
equation still applies if the distance from the lens to the image is taken to be
negative.

A negative lens is thinner in the center than at the edges and produces only virtual images. As seen in figure 3.8, the virtual image produced by a negative lens is closer to the lens than is the object. Again, the thin lens equation is still valid, but both the distance from the image to the lens and the focal length must be taken as negative. Only the distance to the object remains positive.

Curved mirrors also produce images in a manner similar to a lens, as shown in figure 3.9. A concave mirror, as seen in this figure, works in analogy to a positive lens, producing a real or a virtual image depending on whether the object is farther from or closer to the mirror than the mirror’s focal length. A convex mirror acts like a negative lens, always producing a virtual image. The thin lens equation works in both cases as long as the angles are small.

An alternate approach to geometrical optics can be developed from Fermat’s principle. This principle states (in its simplest form) that light waves of a given frequency traverse the path between two points which takes the least time. The most obvious example of this is the passage of light through a homogeneous medium in which the speed of light doesn’t change with position. In this case the shortest time corresponds to the shortest distance between the points, which, as we all know, is a straight line. Thus, Fermat’s principle is consistent with light traveling in a straight line in a homogeneous medium.

Fermat’s principle can also be used to derive the laws of reflection and refraction. For instance, figure 3.10 shows a candidate ray for reflection in which the angles of incidence and reflection are not equal. The time required for the light to go from point A to point B is

| (3.13) |

where c is the speed of light. We find the minimum time by differentiating t with respect to y and setting the result to zero, with the result that

| (3.14) |

However, we note that the left side of this equation is simply sin θ_{I}, while the
right side is sin θ_{R}, so that the minimum time condition reduces to sin θ_{I} = sin θ_{R}
or θ_{I} = θ_{R}, which is the law of reflection.

A similar analysis may be done to derive Snell’s law of refraction. The speed of light in a medium with refractive index n is c∕n, where c is its speed in a vacuum. Thus, the time required for light to go some distance in such a medium is n times the time light takes to go the same distance in a vacuum. Referring to figure 3.11, the time required for light to go from A to B becomes

| (3.15) |

This results in the condition

| (3.16) |

where θ_{R} is now the refracted angle. We recognize this result as Snell’s
law.

Notice that the reflection case illustrates a point about Fermat’s principle: The minimum time may actually be a local rather than a global minimum — after all, in figure 3.10, the global minimum distance from A to B is still just a straight line between the two points! In fact, light starting from point A will reach point B by both routes — the direct route and the reflected route.

It turns out that trajectories allowed by Fermat’s principle don’t strictly have to be minimum time trajectories. They can also be maximum time trajectories, as illustrated in figure 3.12. In this case light emitted at point O can be reflected back to point O from four points on the mirror, A, B, C, and D. The trajectories O-A-O and O-C-O are minimum time trajectories while O-B-O and O-D-O are maximum time trajectories.

Fermat’s principle seems rather mysterious. However, the American physicist Richard Feynman made sense out of it by invoking an even more fundamental principle, as we now see.

If a light ray originates at point O in figure 3.12, reflects off of the ellipsoidal mirror surface at point A, and returns to point O, the elapsed time isn’t much different from that experienced by a ray which reflects off the mirror a slight distance from point A and returns to O. This is because at point A the beam from point O is perpendicular to the tangent to the surface of the mirror at point A. In contrast, the time experienced by a ray going from point O to point E and back would differ by a much greater amount than the time experienced by a ray reflecting off the mirror a slight distance from point E. This is because the tangent to the mirror surface at point E is not perpendicular to the beam from point O.

Technically, the change in the round trip time varies linearly with the deviation in the reflection point from point E, but quadratically with the deviation from point A. If this deviation is small in the first place, then the change in the round trip time will be much smaller for the quadratic case than for the linear case.

It seems odd that we would speak of a beam reflecting back to point O if it hit the mirror at any point except A, B, C, or D, due to the requirements of the law of reflection. However, recall that the law of reflection itself depends on Fermat’s principle, so we cannot assume the validity of that law in this investigation.

Feynman postulated that light rays explore all possible paths from one point to another, but that the only paths realized in nature are those for which light taking closely neighboring paths experiences nearly the same elapsed time (or more generally, traverses nearly the same number of wavelengths) as the original path. If this is true, then neighboring rays interfere constructively with each other, resulting in a much brighter beam than would occur in the absence of this constructive interference. Thus, the round-trip paths O-A-O, O-B-O, O-C-O, and O-D-O in figure 3.12 actually occur, but not O-E-O. Feynman explains Fermat’s principle by invoking constructive and destructive interference!

Figure 3.13 illustrates a rather peculiar situation. Notice that all the rays from point O which intercept the lens end up at point I. This would seem to contradict Fermat’s principle, in that only the minimum (or maximum) time trajectories should occur. However, a calculation shows that all the illustrated trajectories in this particular case take the same time. Thus, the light cannot choose one trajectory over another using Fermat’s principle and all of the trajectories are equally favored. Note that this inference applies not to just any set of trajectories, but only those going from an object point to the corresponding image point.

- The index of refraction varies as shown in figure 3.14:
- Given θ
_{1}, use Snell’s law to find θ_{2}. - Given θ
_{2}, use Snell’s law to find θ_{3}. - From the above results, find θ
_{3}, given θ_{1}. Do n_{2}or θ_{2}matter?

- Given θ
- A 45
^{∘}-45^{∘}-90^{∘}prism is used to totally reflect light through 90^{∘}as shown in figure 3.15. What is the minimum index of refraction of the prism needed for this to work?

- Show graphically which way the wave vector must point inside the calcite
crystal of figure 3.3 for a light ray to be horizontally oriented. Sketch the
orientation of the wave fronts in this case.

- The human eye is a lens which focuses images on a screen called the retina. Suppose that the normal focal length of this lens is 4 cm and that this focuses images from far away objects on the retina. Let us assume that the eye is able to focus on nearby objects by changing the shape of the lens, and thus its focal length. (The lens-retina distance remains the same.) If an object is 20 cm from the eye, what must the altered focal length of the eye be in order for the image of this object to be in focus on the retina?
- An amoeba 0.01 cm in diameter has its image projected on a screen as
shown in figure 3.18 by a positive lens of diameter 0.1 cm.
- How big is the image of the amoeba?
- What is the focal length of the lens?
- What is the minimum distance between features on the image of
the amoeba which can be resolved? Assume that the wavelength
of light used is 5 × 10
^{-7}m. (Hint: What is the spreading angle of a beam of light passing through an opening the size of the lens?)

- The great refractor telescope of Yerkes Observatory in Wisconsin (see figure
3.19) has primary lens D = 1.02 m in diameter with a focal length of
L = 19.4 m. Use the small angle approximation in all calculations and
assume that the light has wavelength 5 × 10
^{-7}m.- Jupiter has a diameter of 1.5 × 10
^{5}km and an average distance from the earth of 8 × 10^{8}km. How big is the image of Jupiter (in cm) at the focal plane of the primary lens? - Given perfect atmospheric “seeing” conditions, how far apart must two features be on Jupiter (in km) for the Yerkes telescope to be able to resolve them?
- What should the focal length l of the secondary lens or eyepiece be
for Jupiter to subtend the same angle as the moon subtends to the
naked eye? The moon’s diameter is 3.5 × 10
^{3}km and its distance from the earth is 3.8 × 10^{5}km. Hint: Imagine that a translucent sheet of ground glass is placed at the focal plane so that the image is seen projected on this ground glass, which scatters light over a broad range of angles. The eyepiece can then be thought of as a magnifying glass with which you can examine the image on the ground glass. Using this artiface, you need consider only light rays that pass through the center of each lens.

- Jupiter has a diameter of 1.5 × 10
- Show that a concave mirror that focuses incoming rays parallel to the optical axis of the mirror to a point on the optical axis, as illustrated in figure 3.16, is parabolic in shape. Hint: Since rays following different paths all move from the distant source to the focal point of the mirror, Fermat’s principle implies that all of these rays take the same time to do so (why is this?), and therefore all traverse the same distance.
- Use Fermat’s principle to explain qualitatively why a ray of light follows the solid rather than the dashed line through the wedge of glass shown in figure 3.17.
- Test your knowledge of Fermat’s principle by using equation (3.15) to derive Snell’s law.

Special Relativity

Albert Einstein invented the special and general theories of relativity early in the 20th century, though many other people contributed to the intellectual climate which made these discoveries possible. The special theory of relativity arose out of a conflict between the ideas of mechanics as developed by Galileo and Newton, and the theory of electromagnetism. For this reason relativity is often discussed in textbooks after electromagnetism is developed. However, special relativity is actually a valid extension to the Galilean world view which is needed when objects move at very high speeds, and it is only coincidentally related to electromagnetism. For this reason we discuss relativity before electromagnetism.

The only fact from electromagnetism that we need is introduced now: There is
a maximum speed at which objects can travel. This is coincidentally equal to the
speed of light in a vacuum, c = 3 × 10^{8} m s^{-1}. Furthermore, a measurement of the
speed of a particular light beam yields the same answer regardless of the speed
of the light source or the speed at which the measuring instrument is
moving.

This rather bizarre experimental result is in contrast to what occurs in
Galilean relativity. If two cars pass a pedestrian standing on a curb, one at
20 m s^{-1} and the other at 50 m s^{-1}, the faster car appears to be moving at
30 m s^{-1} relative to the slower car. However, if a light beam moving at
3 × 10^{8} m s^{-1} passes an interstellar spaceship moving at 2 × 10^{8} m s^{-1}, then the
light beam appears to occupants of the spaceship to be moving at 3 × 10^{8} m s^{-1},
not 1 × 10^{8} m s^{-1}. Furthermore, if the spaceship beams a light signal forward
to its (stationary) destination planet, then the resulting beam appears
to be moving at 3 × 10^{8} m s^{-1} to instruments at the destination, not
5 × 10^{8} m s^{-1}.

The fact that we are talking about light beams is only for convenience. Any other means of sending a signal at the maximum allowed speed would result in the same behavior. We therefore cannot seek the answer to this apparent paradox in the special properties of light. Instead we have to look to the basic nature of space and time.

In order to gain an understanding of both Galilean and Einsteinian relativity it is important to begin thinking of space and time as being different dimensions of a four-dimensional space called spacetime. Actually, since we can’t visualize four dimensions very well, it is easiest to start with only one space dimension and the time dimension. Figure 4.1 shows a graph with time plotted on the vertical axis and the one space dimension plotted on the horizontal axis. An event is something that occurs at a particular time and a particular point in space. (“Julius X. wrecks his car in Lemitar, NM on 21 June at 6:17 PM.”) A world line is a plot of the position of some object as a function of time on a spacetime diagram, although it is conventional to put time on the vertical axis. Thus, a world line is really a line in spacetime, while an event is a point in spacetime. A horizontal line parallel to the position axis is a line of simultaneity in Galilean relativity — i. e., all events on this line occur simultaneously.

In a spacetime diagram the slope of a world line has a special meaning. Notice that a vertical world line means that the object it represents does not move — the velocity is zero. If the object moves to the right, then the world line tilts to the right, and the faster it moves, the more the world line tilts. Quantitatively, we say that

| (4.1) |

in Galilean relativity. Notice that this works for negative slopes and velocities as well as positive ones. If the object changes its velocity with time, then the world line is curved, and the instantaneous velocity at any time is the inverse of the slope of the tangent to the world line at that time.

The hardest thing to realize about spacetime diagrams is that they represent the past, present, and future all in one diagram. Thus, spacetime diagrams don’t change with time — the evolution of physical systems is represented by looking at successive horizontal slices in the diagram at successive times. Spacetime diagrams represent evolution, but they don’t evolve themselves.

The principle of relativity states that the laws of physics are the same
in all inertial reference frames. An inertial reference frame is one that
is not accelerated. Reference frames attached to a car at rest and to a
car moving steadily down the freeway at 30 m s^{-1} are both inertial. A
reference frame attached to a car accelerating away from a stop light is not
inertial.

The principle of relativity is an educated guess or hypothesis based on extensive experience. If the principle of relativity weren’t true, we would have to do all our calculations in some preferred reference frame. This would be very annoying. However, the more fundamental problem is that we have no idea what the velocity of this preferred frame might be. Does it move with the earth? That would be very earth-centric. How about the velocity of the center of our galaxy or the mean velocity of all the galaxies? Rather than face the issue of a preferred reference frame, physicists have chosen to stick with the principle of relativity.

If an object is moving to the left at velocity v relative to a particular reference
frame, it appears to be moving at a velocity v′ = v - U relative to another
reference frame which itself is moving at velocity U. This is the Galilean velocity
transformation law, and it is based on everyday experience. If you are
traveling 30 m s^{-1} down the freeway and another car passes you doing
40 m s^{-1}, then the other car moves past you at 10 m s^{-1} relative to your
car.

Figure 4.2 shows how the world line of an object is represented differently in the unprimed (x, t) and primed (x′, t′) reference frames. The difference between the velocity of the object and the velocity of the primed frame (i. e., the difference in the inverses of the slopes of the corresponding world lines) is the same in both reference frames in this Galilean case. This illustrates the difference between a physical law independent of reference frame (the difference between velocities in Galilean relativity) and the different motion of the object in the two different reference frames.

In special relativity we find that space and time “mix” in a way that they don’t in Galilean relativity. This suggests that space and time are different aspects of the same “thing”, which we call spacetime.

If time and position are simply different dimensions of the same abstract space, then they should have the same units. The easiest way to arrange this is to multiply time by the maximum speed, c, resulting in the kind of spacetime diagram shown in figure 4.3. Notice that world lines of light have slope ±1 when the time axis is scaled this way. Furthermore, the relationship between speed and the slope of a world line must be revised to read

| (4.2) |

Notice that it is physically possible for an object to have a world line which connects event O at the origin and the events A and D in figure 4.3, since the slope of the resulting world line would exceed unity, and thus represent a velocity less than the speed of light. Events which can be connected by a world line are called timelike relative to each other. On the other hand, event O cannot be connected to events B and C by a world line, since this would imply a velocity greater than the speed of light. Events which cannot be connected by a world line are called spacelike relative to each other. Notice the terminology in figure 4.3: Event A is in the past of event O, while event D is in the future. Events B and C are elsewhere relative to event O.

As we learned previously, the principle of relativity states that the laws of physics are the same in all inertial reference frames. The principle of relativity applies to Einsteinian relativity just as it applies to Galilean relativity.

Notice that the constancy of the speed of light in all reference frames is consistent with the principle of relativity. However, as noted above, it is inconsistent with our notions as to how velocities add, or alternatively, how we think the world should look from reference frames moving at different speeds. We have called the classical way of understanding the view from different reference frames Galilean relativity. The new way that reconciles the behavior of objects moving at very high speeds is called Einsteinian relativity. Einstein’s great contribution was to discover the laws that tell us how the world looks from reference frames moving at high speeds relative to each other. These laws constitute a geometry of spacetime, and from them all of special relativity can be derived.

All of the observed facts about spacetime can be derived from two postulates:

- Whether two events are simultaneous depends on the reference frame from which they are viewed.
- Spacetime obeys a modified Pythagorean theorem, which gives the distance,
I, in spacetime or spacetime interval as
(4.3) where X, T, and I are defined in figure 4.4.

Let us discuss these postulates in turn.

The classical way of thinking about simultaneity is so ingrained in our everyday habits that we have a great deal of difficulty adjusting to what special relativity has to say about this subject. Indeed, understanding how relativity changes this concept is the single most difficult part of the theory — once you understand this, you are well on your way to mastering relativity!

Before tackling simultaneity, let us first think about collocation. Two events (such as A and E in figure 4.5) are collocated if they have the same x value. However, collocation is a concept that depends on the reference frame. For instance, George is driving from Boston to Washington, with the line passing through events A and D being his world line. Just as he passes New York he sneezes (event A in figure 4.5). As he drives by Baltimore, he sneezes again (event D). In the reference frame of the earth, these two sneezes are not collocated, since they are separated by many kilometers. However, in the reference frame of George’s car, they occur in the same place — assuming that George hasn’t left the driver’s seat!

Notice that any two events separated by a timelike interval are collocated in some reference frame. The speed of the reference frame is given by equation (4.2), where the slope is simply the slope of the world line connecting the two events.

In Galilean relativity, if two events are simultaneous, we consider them to be simultaneous in all reference frames. For instance, if two clocks, one in New York and one in Los Angeles, strike the hour at the same time in the earth reference frame, then in Galilean relativity these events also appear to be simultaneous to instruments in the space shuttle as it flies over the United States. However, if the space shuttle is moving from west to east, i. e., from Los Angeles toward New York, careful measurements will show that the clock in New York strikes the hour before the clock in Los Angeles! Thus, the Galilean point of view is not accurate.

Just as collocation depends on one’s reference frame, this result shows that simultaneity also depends on the reference frame. Figure 4.5 shows how this works. In figure 4.5 events A and B are simultaneous in the rest or unprimed reference frame. However, in the primed reference frame, events A and C are simultaneous, and event B occurs at an earlier time. If A and B correspond to the clocks striking in Los Angeles and New York respectively, then it is clear that B must occur at an earlier time in the primed frame if indeed A and C are simultaneous in that frame.

The tilted line passing through events A and C in figure 4.5 is called the line of simultaneity for the primed reference frame. Its slope is related to the speed, U, of the reference frame by

| (4.4) |

Notice that this is the inverse of the slope of the world line attached to the primed
reference frame. There is thus a symmetry between the world line and the line of
simultaneity of a moving reference frame — as the reference frame moves faster to
the right, these two lines close like the blades of a pair of scissors on the 45^{∘}
line.

In Galilean relativity it is fairly obvious what we mean by two events being simultaneous — it all boils down to coordinating portable clocks which are sitting next to each other, and then moving them to the desired locations. Two events separated in space are simultaneous if they occur at the same time on clocks located near each event, assuming that the clocks have been coordinated in the above manner.

In Einsteinian relativity this doesn’t work, because the very act of moving the clocks changes the rate at which the clocks run. Thus, it is more difficult to determine whether two distant events are simultaneous.

An alternate way of experimentally determining simultaneity is shown in figure 4.6. Since we know from observation that light travels at the same speed in all reference frames, the pulses of light emitted by the light sources in figure 4.6 will reach the two equidistant observers simultaneously in both cases. The line passing through these two events, A and B, defines a line of simultaneity for both stationary and moving observers. For the stationary observers this line is horizontal, as in Galilean relativity. For the moving observers the light has to travel farther in the rest frame to reach the observer receding from the light source, and it therefore takes longer in this frame. Thus, event B in the right panel of figure 4.6 occurs later than event A in the stationary reference frame and the line of simultaneity is tilted. We see that the postulate that light moves at the same speed in all reference frames leads inevitably to the dependence of simultaneity on reference frame.

The Pythagorean theorem of spacetime differs from the usual Pythagorean theorem in two ways. First, the vertical side of the triangle is multiplied by c. This is a trivial scale factor that gives time the same units as space. Second, the right side of equation (4.3) has a minus sign rather than a plus sign. This highlights a fundamental difference between spacetime and the ordinary xyz space in which we live. Spacetime is said to have a non-Euclidean (but not curved) geometry — in other words, the normal rules of geometry that we learn in high school don’t always work for spacetime!

The main consequence of the minus sign in equation (4.3) is that I^{2} can be
negative and therefore I can be imaginary. Furthermore, in the special case where
X = ±cT, we actually have I = 0 even though X, T≠0 — i. e., the “distance”
between two well-separated events can be zero. Clearly, spacetime has some weird
properties!

The quantity I is usually called an interval in spacetime. Generally speaking,
if I^{2} is positive, the interval is called spacelike, while for a negative I^{2}, the interval
is called timelike.

A concept related to the spacetime interval is the proper time τ. The proper time between the two events A and C in figure 4.4 is defined by the equation

| (4.5) |

Notice that I and τ are related by

| (4.6) |

so the spacetime interval and the proper time are not independent concepts. However, I has the dimensions of length and is real when the events defining the interval are spacelike relative to each other, whereas τ has the dimensions of time and is real when the events are timelike relative to each other. Both equation (4.3) and equation (4.5) express the spacetime Pythagorean theorem.

If two events defining the end points of an interval have the same t value, then the interval is the ordinary space distance between the two events. On the other hand, if they have the same x value, then the proper time is just the time interval between the events. If the interval between two events is spacelike, but the events are not simultaneous in the initial reference frame, they can always be made simultaneous by choosing a reference frame in which the events lie on the same line of simultaneity. Thus, the meaning of the interval in that case is just the distance between the events in the new reference frame. Similarly, for events separated by a timelike interval, the proper time is just the time between two events in a reference frame in which the two events are collocated.

Stationary and moving clocks run at different rates in relativity. This is illustrated in figure 4.7. The triangle ABC in the left panel of figure 4.7 can be used to illustrate this point. Suppose that the line passing through the events A and C in this figure is the world line of a stationary observer. At zero time another observer moving with velocity V passes the stationary observer. The moving observer’s world line passes through events A and B.

We assume that events B and C are simultaneous in the rest frame, so ABC is a right triangle. Application of the spacetime Pythagorean theorem thus yields

| (4.7) |

Since the second observer is moving at velocity V , the slope of his world line is

| (4.8) |

where the right side of the above equation is the slope calculated as the rise of the world line cT over the run X between events A and B. Eliminating X between the above two equations results in a relationship between T and T′:

| (4.9) |

where

| (4.10) |

The quantity γ occurs so often in relativistic calculations that we give it this special symbol. Note by its definition that γ ≥ 1.

Equation (4.9) tells us that the time elapsed for the moving observer is less than that for the stationary observer, which means that the clock of the moving observer runs more slowly. This is called the time dilation effect.

Let us view this situation from the reference frame of the moving observer. In this frame the moving observer becomes stationary and the stationary observer moves in the opposite direction, as illustrated in the right panel of figure 4.7. By symmetric arguments, one infers that the clock of the initially stationary observer who is now moving to the left runs more slowly in this reference frame than the clock of the initially moving observer. One might conclude that this contradicts the previous results. However, examination of the right panel of figure 4.7 shows that this is not so. The interval cT is still greater than the interval cT′, because such intervals are relativistically invariant quantities. However, events B and C are no longer simultaneous, so one cannot use these results to infer anything about the rate at which the two clocks run in this frame. Thus, the relative nature of the concept of simultaneity saves us from an incipient paradox, and we see that the relative rates at which clocks run depends on the reference frame in which these rates are observed.

A similar argument can be made to show how the postulates of relativity result in the Lorentz contraction. Figure 4.8 compares the length X′ of a moving object measured in its own reference frame (left panel) with its length X as measured in a stationary reference frame (right panel). The length of a moving object is measured by simultaneously measuring the positions of the front and the rear of the object and subtracting these two numbers. Events A and C correspond to these position measurements for the stationary reference frame since they are respectively on the rear and front world lines of the object. Thus, the interval AC, which is equal to X, is the length of the object as measured in the stationary frame.

In the left panel, X is the hypotenuse of a right triangle. Therefore, by the Pythagorean theorem of spacetime, we have

| (4.11) |

Now, the line passing through A and C in the left panel is the line of simultaneity of the stationary reference frame. The slope of this line is -V∕c, where V is the speed of the object relative to the stationary reference frame. Geometrically in figure 4.8, the slope of this line is -cT′∕X′, so we find by equating these two expressions for the slope that

| (4.12) |

Finally, eliminating T′ between (4.11) and (4.12) results in

| (4.13) |

This says that the length of a moving object as measured in a stationary reference frame (X) is less than the actual length of the object as measured in its own reference frame (X′). This reduction in length is called the Lorentz contraction.

Note that the Lorentz contraction only occurs in the direction of motion. The dimensions of a moving object perpendicular to the motion remain unchanged.

An interesting application of time dilation is the so-called twin paradox, which
turns out not to be a paradox at all. Two twins are initially the same
age. One twin travels to a distant star on an interstellar spaceship which
moves at speed V , which is close to the speed of light. Upon reaching the
star, the traveling twin immediately turns around and heads home. When
reaching home, the traveling twin has aged less than the twin that stayed
home. This is easily explained by the time dilation effect, which shows
that the proper time elapsed along the world line of the traveling twin is
(1 - V ^{2}∕c^{2})^{1∕2} times the proper time elapsed along the world line of the other
twin.

The “paradox” part of the twin paradox arises from making the symmetric argument in which one assumes the reference frame of the traveling twin to be stationary. The frame of the earth-bound twin must then travel in the sense opposite that of the erstwhile traveling twin, which means that the earth-bound twin must age less rather than more. This substitution is justified on the basis of the principle of relativity, which states that the laws of physics must be the same in all inertial reference frames.

However, the above argument is fallacious, because the reference frame of the traveling twin is not inertial throughout the entire trip, since at various points the spaceship has to accelerate and decelerate. Thus, the principle of relativity cannot be used to assert the equivalence of the traveler’s reference frame to the stationary frame.

Of particular importance is the period of deceleration and acceleration near the distant star. During this interval, the line of simulaneity of the traveling twin rotates, as illustrated in figure 4.9, such that the twin staying at home rapidly ages in the reference frame of the traveler. Thus, even though the acceleration of the traveling twin may occupy only a negligible segment of the twin’s world line, the overall effect is not negligible. In fact, the shorter and more intense the period of acceleration, the more rapidly the earth-bound twin ages in the traveling frame during this interval!

- Sketch your personal world line on a spacetime diagram for the last 24 hours, labeling by time and location special events such as meals, physics classes, etc. Relate the slope of the world line at various times to how fast you were walking, riding in a car, etc.
- Spacetime conversions:
- What is the distance from New York to Los Angeles in seconds? From here to the moon? From here to the sun?
- What is one nanosecond in meters? One second? One day? One year?

- Three events have the following spacetime coordinates: A is at
(x,ct) = (2 m, 1 m); B is at (x,ct) = (-2 m, 0 m); C is at (x,ct) = (0 m, 3 m).
- A world line for an object passes through events B and C. How fast and in which direction is the object moving?
- A line of simultaneity for a coordinate system passes through events A and B. How fast and in which direction is the coordinate system moving?
- What is the invariant interval between events A and B? B and C? A and C?
- Can a signal from event B reach event A? Can it reach event C? Explain.

Hint: Draw a spacetime diagram with all the events plotted before trying to answer the above questions.

- In the following problem be sure to indicate the slope of all pertinent lines
drawn.
- In a spacetime diagram, sketch a line of simultaneity for a reference frame moving to the left at V = c∕2, where c is the speed of light.
- Sketch the world line of an object which is initially stationary, but which accelerates to a velocity of v = c∕3.

- If the slopes of the world lines of the observers in the right panel of figure 4.6 are both 1∕β, find the slope of their line of simultaneity, AB.
- Suppose that an interstellar spaceship goes a distance X = 100 light years
relative to the rest frame in T′ = 10 years of its own time.
- Draw a spacetime diagram in the rest frame showing X, T′, and the time T needed for this journey relative to the rest frame.
- Compute T, using your spacetime diagram as an aid.
- Compute the speed of the spaceship.

- If an airline pilot flies 80 hr per month (in the rest frame) at 300 m s
^{-1}for 30 years, how much younger will she be than her twin brother (who handles baggage) when she retires? Hint: Use (1 + ϵ)^{x}≈ 1 + xϵ for small ϵ. - A mu particle normally lives about 2 × 10
^{-6}sec before it decays. However, muons created by cosmic rays 20 km up in the atmosphere reach the Earth’s surface. How fast must they be going? - The Stanford Linear Accelerator accelerates electrons to a speed such that
the 3 km long accelerator appears to be 8 cm long to the electron,
due to the Lorentz contraction. How much less than the speed of
light is the electron traveling? Hint: It is best to first develop an
approximation for the relationship between γ = (1 - v
^{2}∕c^{2})^{-1∕2}and the difference between c and v for a particle moving close to the speed of light. - How fast do you have to go to reach the center of our galaxy in your expected lifetime? At this speed, what does this distance appear to be? (We are about 30000 light years from the galactic center.)
- Two identical spaceships pass each other going in the opposite direction at
the same speed.
- Sketch a spacetime diagram showing the world lines of the front and rear of each spaceship as well as lines of simultaneity for each spaceship.
- Indicate an interval on the diagram corresponding to the rightward-moving spaceship’s length in its own reference frame.
- Indicate an interval corresponding to the leftward-moving spaceship’s length in the reference frame of the rightward-moving spaceship.
- Indicate an interval equal to the length of either spaceship in the rest frame.

- George and Sally are twins initially separated by a distance d and at rest. In
the rest frame they are initially the same age. At time t = 0 both George
and Sally get in their spaceships and head to the right at velocity U. Both
move a distance d to the right and decelerate to a halt. (See figure
4.10.)
- When both are moving, how far away is Sally according to George?
- How much older or younger is Sally relative to George while both are moving?
- How much older or younger is Sally relative to George after both stop?

Hint: Draw the triangle ABC in a reference frame moving with George and Sally.

Applications of Special Relativity

In this chapter we continue the study of special relativity. Three important applications of the ideas developed in the previous chapter are made here. First, we show how to describe waves in the context of spacetime. We then see how waves which have no preferred reference frame (such as that of a medium supporting them) are constrained by special relativity to have a dispersion relation of a particular form. This dispersion relation turns out to be that of the relativistic matter waves of quantum mechanics. Second, we investigate the Doppler shift phenomenon, in which the frequency of a wave takes on different values in different coordinate systems. Third, we show how to add velocities in a relativistically consistent manner.

A new mathematical idea is presented in the context of relativistic waves, namely the spacetime vector or four-vector. Writing the laws of physics totally in terms of relativistic scalars and four-vectors insures that they will be valid in all inertial reference frames.

We now look at the characteristics of waves in spacetime. Recall that a sine wave moving to the right in one space dimension can be represented by

| (5.1) |

where A_{0} is the (constant) amplitude of the wave, k is the wavenumber, and ω is
the angular frequency, and that the quantity ϕ = kx - ωt is called the phase of
the wave. For a plane wave in three space dimensions, the wave is represented in a
similar way,

| (5.2) |

where x is now the position vector and k is the wave vector. The magnitude of the wave vector, |k| = k is just the wavenumber of the wave and the direction of this vector indicates the direction the wave is moving. The phase of the wave in this case is ϕ = k ⋅ x - ωt.

In the case of a one-dimensional wave moving to the right ϕ = kx-ωt. A wave front has constant phase ϕ, so solving this equation for t and multiplying by c, the speed of light in a vacuum, gives us an equation for the world line of a wave front:

| (5.3) |

The slope of the world line in a spacetime diagram is the coefficient of x, or c∕u_{p},
where u_{p} = ω∕k is the phase speed. The world lines of the wave fronts of a wave
are illustrated in figure 5.1.

Also shown in figure 5.1 is a spacetime vector or four-vector which represents the frequency and wavenumber of the wave, which we refer to as the wave four-vector. It is called a four-vector because it has 3 spacelike components and one timelike component when there are 3 space dimensions. In the case shown, there is only a single space dimension. The spacelike component of the wave four-vector is just k (or k when there are 3 space dimensions), while the timelike component is ω∕c. The c is in the denominator to give the timelike component the same dimensions as the spacelike component. From figure 5.1 it is clear that the slope of the line representing the four-vector is ω∕ck, which is just the inverse of the slope of the wave fronts.

Let us define some terminology. We indicate a four-vector by underlining and write the components in the following way: k = (k,ω∕c), where k is the wave four-vector, k is its spacelike component, and ω∕c is its timelike component. For three space dimensions, where we have a wave vector rather than just a wavenumber, we write k = (k,ω∕c).

Another example of a four-vector is simply the position vector in spacetime, x = (x,ct), or x = (x,ct) in three space dimensions. The c multiplies the timelike component in this case, because that is what is needed to give it the same dimensions as the spacelike component.

In three dimensions we define a vector as a quantity with magnitude and
direction. Extending this to spacetime, a four-vector is a quantity with magnitude
and direction in spacetime. Implicit in this definition is the notion that the
vector’s magnitude is a quantity independent of coordinate system or reference
frame. We have seen that the invariant interval in spacetime from the origin to the
point (x,ct) is I = (x^{2} - c^{2}t^{2})^{1∕2}, so it makes sense to identify this as the
magnitude of the position vector. This leads to a way of defining a dot product of
four-vectors. Given two four-vectors A = (A,A_{t}) and B = (B,B_{t}), the dot
product is

| (5.4) |

This is consistent with the definition of invariant interval if we set A = B = x,
since then x ⋅ x = x^{2} - c^{2}t^{2} = I^{2}.

In the odd geometry of spacetime it is not obvious what “perpendicular” means. We therefore define two four-vectors A and B to be perpendicular if their dot product is zero, A ⋅ B = 0, in analogy with ordinary vectors.

The dot product of two four-vectors is a scalar result, i. e., its value is
independent of coordinate system. This can be used to our advantage on occasion.
For instance, consider the dot product of a four-vector A which resolves into
(A_{x},A_{t}) in the unprimed frame. Let us further suppose that the spacelike
component is zero in some primed frame, so that the components in this frame are
(0,A_{t}′). The fact that the dot product is independent of coordinate system means
that

| (5.5) |

This constitutes an extension of the spacetime Pythagorean theorem to four-vectors other than the position four-vector. Thus, for instance, the wavenumber for some wave may be zero in the primed frame, which means that the wavenumber and frequency in the unprimed frame are related to the frequency in the primed frame by

| (5.6) |

Returning to the phase of a wave, we immediately see that

| (5.7) |

Thus, a compact way to rewrite equation (5.2) is

| (5.8) |

Since x is known to be a four-vector and since the phase of a wave is known to be a scalar independent of reference frame, it follows that k is also a four-vector rather than just a set of numbers. Thus, the square of the length of the wave four-vector must also be a scalar independent of reference frame:

| (5.9) |

Let us review precisely what this means. As figure 5.2 shows, we can resolve a
position four-vector into components in two different reference frames,
x = (X,cT) = (X′,cT′). However, even though X≠X′ and T≠T′, the vector
lengths computed from these two sets of components are necessarily the same:
x ⋅ x = X^{2} - c^{2}T^{2} = X′^{2} - c^{2}T′^{2}.

Applying this to the wave four-vector, we infer that

| (5.10) |

where the unprimed and primed values of k and ω refer to the components of the wave four-vector in two different reference frames.

Up to now, this argument applies to any wave. However, waves can be divided into two categories, those for which a “special” reference frame exists, and those for which there is no such special frame. As an example of the former, sound waves look simplest in the reference frame in which the gas carrying the sound is stationary. The same is true of light propagating through a material medium with an index of refraction not equal to unity. In both cases the speed of the wave is the same in all directions only in the frame in which the material medium is stationary.

Suppose we have a machine that produces a wave with wavenumber k and frequency ω in its own rest frame. If we observe the wave from a moving reference frame, the wavenumber and frequency will be different, say, k′ and ω′. However, these quantities will be related by equation (5.10).

Up to this point the argument applies to any wave whether a special reference frame exists or not; the observed changes in wavenumber and frequency have nothing to do with the wave itself, but are just consequences of how we have chosen to observe it. However, if there is no special reference frame for the type of wave under consideration, then the same result can be obtained by keeping the observer stationary and moving the wave-producing machine in the opposite direction. By moving it at various speeds, any desired value of k′ can be obtained in the initial reference frame (as opposed to some other frame), and the resulting value of ω′ can be computed using equation (5.10).

This is actually an amazing result. We have shown on the basis of the principle of relativity that any wave type for which no special reference frame exists can be made to take on a full range of frequencies and wavenumbers in any given reference frame, and furthermore that these frequencies and wavenumbers obey

| (5.11) |

Equation (5.11) comes from solving equation (5.10) for ω^{2} and the constant μ^{2}
equals the constant in equation (5.10) times -c^{2}. Equation (5.11) relates
frequency to wavenumber and therefore is the dispersion relation for such
waves. We call waves which have no special reference frame and therefore
necessarily obey equation (5.11) relativistic waves. The only difference in the
dispersion relations between different types of relativistic waves is the
value of the constant μ. The meaning of this constant will become clear
later.

Light in a vacuum is an example of a wave for which no special reference frame exists. For light, μ = 0, and we have (taking the positive root) ω = ck. This simply states what we know already, namely that the phase speed of light in a vacuum is c.

If μ≠0, waves of this type are dispersive. The phase speed is

| (5.12) |

This phase speed always exceeds c, which at first seems like an unphysical conclusion. However, the group velocity of the wave is

| (5.13) |

which is always less than c. Since wave packets and hence signals propagate at the group velocity, waves of this type are physically reasonable even though the phase speed exceeds the speed of light.

Another interesting property of such waves is that the wave four-vector is
parallel to the world line of a wave packet in spacetime. This is easily shown by
the following argument. As figure 5.1 shows, the spacelike component of a wave
four-vector is k, while the timelike component is ω∕c. The slope of the four-vector
on a spacetime diagram is therefore ω∕kc. However, the slope of the world line of
a wave packet moving with group velocity u_{g} is c∕u_{g} = ω∕(kc), which is the same
as the slope of the k four-vector.

Note that when k = 0 we have ω = μ. In this case the group velocity of the wave is zero. For this reason we call μ the rest frequency of the wave.

You have probably heard how the pitch of a train horn changes as it passes you. When the train is approaching, the pitch or frequency is higher than when it is moving away from you. This is called the Doppler shift. A similar, but distinct shift occurs if you are moving past a source of sound. If a stationary whistle is blowing, the pitch perceived from a moving car is higher while moving toward the source than when moving away. The first case thus has a moving source, while the second case has a moving observer.

In this section we compute the Doppler shift as it applies to light moving through a vacuum. Figure 5.3 shows the geometry for computing the time between wave fronts of light for a stationary and a moving reference frame. The time in the stationary frame is just T. Since the world lines of the wave fronts have a slope of unity, the sides of the shaded triangle have the same value, X. If the observer is moving at speed U, the slope of the observer’s world line is c∕U, which means that c∕U = (cT + X)∕X. Solving this for X yields X = UT∕(1 - U∕c), which can then be used to compute T′ = T + X∕c = T∕(1 - U∕c). This formula as it stands leads to the classical Doppler shift for a moving observer. However, with relativistic velocities, one additional factor needs to be taken into account: The observer experiences time dilation since he or she is moving. The actual time measured by the observer between wave fronts is actually

| (5.14) |

where the last step uses 1 - U^{2}∕c^{2} = (1 - U∕c)(1 + U∕c). From this we infer the
relativistic Doppler shift formula for light in a vacuum:

| (5.15) |

where the frequency measured by the moving observer is ω′ = 2π∕τ and the frequency observed in the stationary frame is ω = 2π∕T.

We could go on to determine the Doppler shift resulting from a moving source. However, by the principle of relativity, the laws of physics should be the same in the reference frame in which the observer is stationary and the source is moving. Furthermore, the speed of light is still c in that frame. Therefore, the problem of a stationary observer and a moving source is conceptually the same as the problem of a moving observer and a stationary source when the wave is moving at speed c. This is unlike the case for, say, sound waves, where the stationary observer and the stationary source yield different formulas for the Doppler shift.

Figure 5.4 shows the world line of a moving object from the point of view of two different reference frames, with the primed frame (left panel) moving to the right at speed U relative to the unprimed frame (right panel). The goal is to calculate the velocity of the object relative to the unprimed frame, v, assuming its velocity in the primed frame, v′ is known. The classical result is simply

| (5.16) |

However, this is inconsistent with the speed of light being constant in all reference frames, since if we substitute c for v′, this formula predicts that the speed of light in the unprimed frame is U + c.

We can use the geometry of figure 5.4 to come up with the correct relativistic formula. From the right panel of this figure we infer that

| (5.17) |

This follows from the fact that the slope of the world line of the object in this frame is c∕v. The slope is calculated as the ratio of the rise, c(T + ΔT), to the run, X + ΔX.

From the left panel of figure 5.4 we similarly see that

| (5.18) |

However, we can apply our calculations of Lorentz contraction and time dilation
from the previous chapter to triangles ABD and OAE in the right panel.
The slope of AB is U∕c because AB is horizontal in the left panel, so
X′ = ΔX(1 - U^{2}∕c^{2})^{1∕2}. Similarly, the slope of OA is c∕U since OA is vertical in
the left panel, and T′ = T(1 - U^{2}∕c^{2})^{1∕2}. Substituting these formulas into the
equation for v′∕c yields

| (5.19) |

Again using what we know about the triangles ABD and OAE, we see that

| (5.20) |

Finally, we calculate ΔT∕T by noticing that

| (5.21) |

Substituting equations (5.19), (5.20), and (5.21) into equation (5.17) and simplifying yields the relativistic velocity addition formula:

| (5.22) |

Notice how this equation behaves in various limits. If |Uv′|≪ c^{2}, the
denominator of equation (5.22) is nearly unity, and the special relativistic formula
reduces to the classical case. On the other hand, if v′ = c, then equation (5.22)
reduces to v = c. In other words, if the object in question is moving at the speed
of light in one reference frame, it is moving at the speed of light in all
reference frames, i. e., for all possible values of U. Thus, we have found a
velocity addition formula that 1) reduces to the classical formula for low
velocities and 2) gives the observed results for very high velocities as
well.

Equation (5.22) is valid even if v′ is negative, i. e., if the object is moving to the right less rapidly than the primed reference frame, or even if it is moving to the left in the unprimed frame.

- Sketch the wave fronts and the k four-vector in a spacetime diagram
for the case where ω∕k = 2c. Label your axes and space the wave fronts
correctly for the case k = 4π m
^{-1}. - If the four-vector k = (0, 1 nm
^{-1}) in the rest frame, find the space and time components of k in a frame moving to the left at speed c∕2. - Let’s examine the four-vector u = (u
_{g},c)∕(1 - β^{2})^{1∕2}where β = u_{ g}∕c, u_{g}being the velocity of some object. The four-vector u is called the four-velocity.- Show that u is parallel to the world line of the object.
- Show that u ⋅ u = -c
^{2}. - If u
_{g}is the group velocity of a relativistic wave packet, show that k = (μ∕c^{2})u, where k is the central wave four-vector of the wave packet.

- Find the Doppler shift for a moving source of light from figure 5.5, roughly following the procedure used in the text to find the shift for a moving observer. (Assume that the source moves to the left at speed U.) Is the result the same as for the moving observer, as demanded by the principle of relativity?
- Find the Doppler shift for a stationary source of sound (c ≈ 380 m s
^{-1}in a stationary atmosphere) as measured by a moving observer. Follow the procedure as for light in section 5.5 except note that the proper time τ is virtually the same as the time T′ in figure 5.3 since speeds are much less than the speed of light. - Suppose you shine a laser with frequency ω and wavenumber k on a mirror
moving toward you at speed v, as seen in figure 5.6. What are the frequency
ω′ and wavenumber k′ of the reflected beam? Hint: Find the frequency of
the incident beam in the reference frame of the mirror. The frequency of
the reflected beam will be the same as that of the incident beam
in this frame. Then transform back to the reference frame of the
laser.

- Suppose the moving twin in the twin paradox has a powerful telescope so that she can watch her twin brother back on earth during the entire trip. Describe how the earthbound twin appears to age to the travelling twin compared to her own rate of aging. Use a spacetime diagram to illustrate your argument and consider separately the outbound and return legs. Remember that light travels at the speed of light! Hint: Does the concept of Doppler shift help here?
- Find the velocity of an object with respect to the rest frame if it is
moving at a velocity of 0.1c with respect to another frame which itself
is moving in the same direction at 0.1c relative to the rest frame
using
- the Galilean formula and
- the formula of special relativity.

Determine the fractional error made in using the Galilean formula.

- Each stage of a high performance 3 stage rocket can accelerate to a speed of 0.9c from rest. If the rocket starts from rest, how fast does the final stage eventually go?
- An interstellar spaceship is going from Earth to Sirius with speed U = 0.8c relative to the rest frame. It passes a spaceship which is going from Sirius to Earth at a speed of 0.95c in the reference frame of the first spaceship. What is the velocity (direction and speed) of the second spaceship in the rest frame?

Acceleration and General Relativity

General relativity is Einstein’s extension of special relativity to include gravity. An important aspect of general relativity is that spacetime is no longer necessarily flat, but in fact may be curved under the influence of mass. Understanding curved spacetime is an advanced topic which is not easily accessible at the level of this text. However, it turns out that some insight into general relativistic phenomena may be obtained by investigating the effects of acceleration in the flat (but non-Euclidean) space of special relativity.

The central assumption of general relativity is the equivalence principle, which states that gravity is a force which arises from being in an accelerated reference frame. To understand this we must first investigate the concept of acceleration. We then see how this leads to phenomena such as the gravitational red shift, event horizons, and black holes. We also introduce in a preliminary way the notions of force and mass.

Imagine that you are in a powerful luxury car stopped at a stoplight. As you sit there, gravity pushes you into the comfortable leather seat. The light turns green and you “floor it”. The car accelerates and an additional force pushes you into the seat back. You round a curve, and yet another force pushes you toward the outside of the curve. (But the well designed seat and seat belt keep you from feeling discomfort!)

Let us examine the idea of acceleration more closely. Considering first acceleration in one dimension, figure 6.1 shows the position of an object as a function of time, x(t). The velocity is simply the time rate of change of the position:

| (6.1) |

The acceleration is the time rate of change of velocity:

| (6.2) |

In figure 6.1, only the segment OA has zero velocity. Velocity is increasing in AB, and the acceleration is positive there. Velocity is constant in BC, which means that the acceleration is zero. Velocity is decreasing in CD, and the acceleration is negative. Finally, in DE, the velocity is negative and the acceleration is zero.

In two or three dimensions, position x, velocity v, and acceleration a are all vectors, so that the velocity is

| (6.3) |

while the acceleration is

| (6.4) |

Thus, over some short time interval Δt, the changes in x and v can be written

| (6.5) |

These are vector equations, so the subtractions implied by the “delta” operations must be done vectorially. An example where the vector nature of these quantities is important is motion in a circle at constant speed, which is discussed in the next section.

Imagine an object constrained by an attached string to move in a circle at constant speed, as shown in the left panel of figure 6.2. We now demonstrate that the acceleration of the object is toward the center of the circle. The acceleration in this special case is called the centripetal acceleration.

Figure 6.3 shows the position of the object at two times spaced by the time interval Δt. The position vector of the object relative to the center of the circle rotates through an angle Δθ during this interval, so the angular rate of revolution of the object about the center is ω = Δθ∕Δt. The magnitude of the velocity of the object is v, so the object moves a distance vΔt during the time interval. To the extent that this distance is small compared to the radius r of the circle, the angle Δθ = vΔt∕r. Solving for v and using ω = Δθ∕Δt, we see that

| (6.6) |

The direction of the velocity vector changes over this interval, even though the magnitude v stays the same. Figure 6.3 shows that this change in direction implies an acceleration a which is directed toward the center of the circle, as noted above. The magnitude of the vectoral change in velocity in the time interval Δt is aΔt. Since the angle between the initial and final velocities is the same as the angle Δθ between the initial and final radius vectors, we see from the geometry of the triangle in figure 6.3 that aΔt∕v = Δθ. Solving for a results in

| (6.7) |

Combining equations (6.6) and (6.7) yields the equation for centripetal acceleration:

| (6.8) |

The second form is obtained by eliminating ω from the first form using equation (6.6).

We have a good intuitive feel for the concepts of force and mass because they are very much a part of our everyday experience. We think of force as how hard we push on something. Mass is the resistance of an object to acceleration if it is otherwise free to move. Thus, pushing on a bicycle on a smooth, level road causes it to accelerate more readily than pushing on a car. We say that the car has greater mass. We can summarize this relationship with Newton’s second law

| (6.9) |

where F is the total force on an object, m is its mass, and a is the acceleration resulting from the force.

Three provisos apply to equation (6.9). First, it only makes sense in unmodified form when the velocity of the object is much less than the speed of light. For relativistic velocities it is best to write this equation in a slightly different form which we introduce later. Second, the force must be the total force, including all frictional and other incidental forces which might otherwise be neglected by an uncritical observer. Third, it only works in a reference frame which itself is unaccelerated, i.e., an inertial reference frame. We deal below with accelerated reference frames.

As noted above, acceleration is just the time rate of change of velocity. We use the above results to determine how acceleration transforms from one reference frame to another. Figure 6.4 shows the world line of an accelerated reference frame, with a time-varying velocity U(t) relative to the unprimed inertial rest frame. Defining ΔU = U(T) - U(0) as the change in the velocity of the accelerated frame (relative to the unprimed frame) between events A and C, we can relate this to the change of velocity, ΔU′, of the accelerated frame relative to an inertial frame moving with the initial velocity, U(0). Applying the equation for the relativistic addition of velocities, we find

| (6.10) |

We now note that the mean acceleration of the reference frame between events A and C in the rest frame is just a = ΔU∕T, whereas the mean acceleration in the primed frame between the same two events is a′ = ΔU′∕T′. From equation (6.10) we find that

| (6.11) |

and the acceleration of the primed reference frame as it appears in the unprimed frame is

| (6.12) |

Since we are interested in the instantaneous rather than the average
acceleration, we let T become small. This has three consequences. First, ΔU
and ΔU′ become small, which means that the term U(0)ΔU′∕c^{2} in the
denominator of equation (6.12) can be ignored compared to 1. This means
that

| (6.13) |

with the approximation becoming exact as T → 0. Second, the “triangle” with the
curved side in figure 6.4 becomes a true triangle, with the result that
T′ = T[1 - U(0)^{2}∕c^{2}]^{1∕2}. The acceleration of the primed frame with respect to an
inertial frame moving at speed U(0) can therefore be written

| (6.14) |

Third, we can replace U(0) with U, since the velocity of the accelerated frame doesn’t change very much over a short time interval.

Dividing equation (6.13) by equation (6.14) results in a relationship between the two accelerations:

| (6.15) |

which shows that the acceleration of a rapidly moving object, a, as observed
from the rest frame, is less than its acceleration relative to an inertial
reference frame in which the object is nearly stationary, a′, by the factor
(1 - U^{2}∕c^{2})^{3∕2}. We call this latter acceleration the intrinsic acceleration. This
difference in observed acceleration between the two inertial reference frames is
purely the result of the geometry of spacetime, but it has interesting
consequences.

Identifying a with dU∕dt, we can integrate the acceleration equation assuming that the intrinsic acceleration a′ is constant and that the velocity U = 0 at time t = 0. We get the following result (verify this by differentiating with respect to time):

| (6.16) |

which may be solved for U∕c:

| (6.17) |

This is plotted in figure 6.5. Classically, the velocity U would reach the speed of light when a′t∕c = 1. However, as figure 6.5 shows, the rate at which the velocity increases with time slows as the object moves faster, such that U approaches c asymptotically, but never reaches it.

The results for this section are valid only for acceleration components in the direction of motion. The components perpendicular to this direction behave differently and are treated in more advanced texts.

Referring back to the forces being felt by the occupant of a car, it is clear that the forces associated with accelerations are directed opposite the accelerations and proportional to their magnitudes. For instance, when accelerating away from a stoplight, the acceleration is forward and the perceived force is backward. When turning a corner, the acceleration is toward the corner while the perceived force is away from the corner. Such forces are called inertial forces.

The origin of these forces can be understood by determining how acceleration changes when one observes it from a reference frame which is itself accelerated. Suppose that the primed reference frame is accelerating to the right with acceleration A relative to the unprimed frame. The position x′ in the primed frame can be related to the position x in the unprimed frame by

| (6.18) |

where X is the position of the origin of the primed frame in the unprimed frame. Taking the second time derivative, we see that

| (6.19) |

where a = d^{2}x∕dt^{2} is the acceleration in the unprimed frame and a′ = d^{2}x′∕dt^{2} is
the acceleration according to an observer in the primed frame.

We now substitute this into equation (6.9) and move the term involving A to the left side:

| (6.20) |

This shows that Newton’s second law represented by equation (6.9) is not valid in an accelerated reference frame, because the total force F and the acceleration a′ in this frame don’t balance as they do in the unaccelerated frame; the additional term -mA messes up this balance.

We can fix this problem by considering -mA to be a type of force, in which case we can include it as a part of the total force F. This is the inertial force which we mentioned above. Thus, to summarize, we can make Newton’s second law work when objects are observed from accelerated reference frames if we include as part of the total force an inertial force which is equal to -mA, A being the acceleration of the reference frame of the observer and m the mass of the object being observed.

The right panel of figure 6.2 shows the inertial force observed in the reference frame of an object moving in circular motion at constant speed. In the case of circular motion the inertial force is called the centrifugal force. It points away from the center of the circle and just balances the tension in the string. This makes the total force on the object zero in its own reference frame, which is necessary since the object cannot move (or accelerate) in this frame.

General relativity says that gravity is nothing more than an inertial force. This was called the equivalence principle by Einstein. Since the gravitational force on the Earth points downward, it follows that we must be constantly accelerating upward as we stand on the surface of the Earth! The obvious problem with this interpretation of gravity is that we don’t appear to be moving away from the center of the Earth, which would seem to be a natural consequence of such an acceleration. However, relativity has a surprise in store for us here.

It follows from the above considerations that something can be learned about general relativity by examining the properties of accelerated reference frames. In particular, we can gain insight into the above apparent paradox. Equation (6.17) shows that the velocity of an object undergoing constant intrinsic acceleration a (note that we have dropped the “prime” from a to simplify the notation) is

| (6.21) |

where t is the time and c is the speed of light. A function x(t) which satisfies equation (6.21) is

| (6.22) |

(Verify this by differentiating it.) The interval OB in figure 6.6 is of length
x(0) = c^{2}∕A.

The slanted line OA is a line of simultaneity associated with the unaccelerated world line tangent to the accelerated world line at point A. This line of simultaneity goes through the origin, as is shown in figure 6.6. To demonstrate this, multiply equations (6.21) and (6.22) together and solve for v∕c:

| (6.23) |

From figure 6.6 we see that ct∕x is the slope of the line OA, where (x, ct) are the coordinates of event A. Equation (6.23) shows that this line is indeed the desired line of simultaneity, since its slope is the inverse of the slope of the world line, c∕v. Since there is nothing special about the event A, we infer that all lines of simultaneity associated with the accelerated world line pass through the origin.

We now inquire about the length of the invariant interval OA in figure 6.6.
Recalling that I^{2} = x^{2} -c^{2}t^{2} and using equation (6.22), we find that the length of
OA is

| (6.24) |

which is the same as the length of the interval OB. By extension, all events on the accelerated world line are the same invariant interval from the origin. Recalling that the interval along a line of simultaneity is the distance in the associated reference frame, we reach the astonishing conclusion that even though the object associated with the curved world line in figure 6.6 is accelerating away from the origin, it always remains the same distance (in its own frame) from the origin.

The analogy between this problem and the apparent paradox in which one remains a fixed distance from the center of the earth while accelerating away from it is not perfect. In particular, the earth case depends on the existence of the earth’s mass.

Light emitted at a lower level in a gravitational field has its frequency reduced as it travels to a higher level. This phenomenon is called the gravitational red shift. Figure 6.26 shows why this happens. Since experiencing a gravitational force is equivalent to being in an accelerated reference frame, we can use the tools of special relativity to view the process of light emission and absorption from the point of view of the unaccelerated or inertial frame. In this reference frame the observer of the light is accelerating to the right, as indicated by the curved world line in figure 6.26, which is equivalent to a gravitational force to the left. The light is emitted at point A with frequency ω by a source which is stationary at this instant. At this instant the observer is also stationary in this frame. However, by the time the light gets to the observer, he or she has a velocity to the right which means that the observer measures a Doppler shifted frequency ω′ for the light. Since the observer is moving away from the source, ω′ < ω, as indicated above.

The relativistic Doppler shift is given by

| (6.25) |

so we need to compute U∕c. The line of simultaneity for the observer at point B
goes through the origin, and is thus given by line segment OB in figure 6.26. The
slope of this line is U∕c, where U is the velocity of the observer at point B.
From the figure we see that this slope is also given by the ratio X′∕X.
Equating these, eliminating X in favor of L = (X^{2} - X′^{2})^{1∕2} = c^{2}∕g,
which is the actual invariant distance of the observer from the origin, and
substituting into equation (6.25) results in our gravitational red shift
formula:

| (6.26) |

If X′ = 0, then there is no redshift, because the source is collocated with the observer. On the other hand, if the source is located at the origin, so X′ = X, the Doppler shifted frequency is zero. In addition, the light never gets to the observer, since the world line is asymptotic to the light world line passing through the origin. If the source is at a higher level in the gravitational field than the observer, so that X′ < 0, then the frequency is shifted to a higher value, i. e., it becomes a “blue shift”.

Equation (6.26) works for more complex geometries than that associated with
an accelerated reference frame, e.g., for the gravitational field g associated with a
star, as long as |X′|≪ c^{2}∕g. In this case L is no longer the distance to the center
of the star but remains equal to c^{2}∕g.

The 45^{∘} diagonal line passing through the origin in figure 6.6 is called
the event horizon for the accelerated observer in this figure. Notice that
light from the “twilight zone” above and to the left of the event horizon
cannot reach the accelerated observer. However, the reverse is not true
— a light signal emitted to the left by the observer can cross the event
horizon into the “twilight zone”. The event horizon thus has a peculiar
one-way character — it passes signals from right to left, but not from left to
right.

- An object moves as described in figure 6.8, which shows its position x as a
function of time t.
- Is the velocity positive, negative, or zero at each of the points A, B, C, D, E, and F?
- Is the acceleration positive, negative, or zero at each of the points A, B, C, D, E, and F?

- An object is moving counterclockwise at constant speed around the circle
shown in figure 6.9 due to the fact that it is attached by a string to the
center of the circle at point O.
- Sketch the object’s velocity vectors at points A, B, and C.
- Sketch the object’s acceleration vectors at points A, B, and C.
- If the string breaks at point A, sketch the subsequent trajectory followed by the object.

- How fast are you going after accelerating from rest with intrinsic
acceleration a = 10 m s
^{-2}after the given times measured in the rest frame:- 10 y?
- 100000 y?

Express your answer as the speed of light minus your actual speed. Hint: You may have a numerical problem on the second part, which you should try to resolve using the approximation (1 + ϵ)

^{x}≈ 1 + xϵ, which is valid for |ϵ|≪ 1. - An object’s world line is defined by x(t) = (d
^{2}+ c^{2}t^{2})^{1∕2}where d is a constant and c is the speed of light.- Find the object’s velocity as a function of time.
- Using the above result, find the slope of the tangent to the world line as a function of time.
- Find where the line of simultaneity corresponding to each tangent world line crosses the x axis.

- A car accelerates in the positive x direction at 3 m s
^{-2}.- What is the net force on a 100 kg man in the car as viewed from an inertial reference frame?
- What is the inertial force experienced by this man in the reference frame of the car?
- What is the net force experienced by the man in the car’s (accelerated) reference frame?

- A person is sitting in a comfortable chair in her home in Bogotá, Colombia,
which is essentially on the equator.
- What would the rotational period of the earth have to be to make this person weightless?
- What is her acceleration according to the equivalence principle in the earth frame in this situation?

- At time t = 0 a Zork (a creature from the planet Zorkheim) accelerating to
the right at a = 10
^{3}m s^{-2}in a spaceship accidently drops its stopwatch from the spaceship just when its velocity is zero.- Describe qualitatively how the hands of the watch appear to move to the Zork as it observes the watch through a powerful telescope.
- After a very long time what does the watch read?

Hint: Draw a spacetime diagram with the world lines of the spaceship and the watch. Then send light rays from the watch to the spaceship.

- Using a spacetime diagram, show why signals from events on the hidden side of the event horizon from an accelerating spaceship cannot reach the spaceship.
- Approximate equation (6.26) to first order in X′ for the case in which X′≪ L.
- Imagine two identical clocks, one on top of the volcano Chimborazo in
Ecuador (6300 m above sea level), the other in the Ecuadorian city of
Guayaquil (at sea level).
- From the perspective of Chimborazo, does the clock in Guayaquil appear to be running faster or slower than the Chimborazo clock? Explain.
- Compute the fractional frequency difference (ω - ω′)∕ω in this case, where ω is the freqency of the Guayaquil clock as observed in Guayaquil (and the frequency of the Chimborazo clock on Chimborazo) and ω′ is the frequency of the Guayaquil clock as observed from Chimborazo. You may wish to use the results of the previous problem.

Matter Waves

We begin our study of quantum mechanics by discussing the diffraction undergone by X-rays and electrons when they interact with a crystal. X-rays are a form of electromagnetic radiation with wavelengths comparable to the distances between atoms. Scattering from atoms in a regular crystalline structure results in an interference pattern which is in many ways similar to the pattern from a diffraction grating. We first develop Bragg’s law for diffraction of X-rays from a crystal. Two practical techniques for doing X-ray diffraction are then described.

It turns out that electrons have wave-like properties and also undergo Bragg diffraction by crystals. Bragg diffraction thus provides a crucial bridge between the worlds of waves and particles. With this bridge we introduce the classical ideas of momentum and energy by relating them to the wave vector and frequency of a wave. The properties of waves also give rise to the Heisenberg uncertainty principle.

Table 7.1 shows a table of the Nobel prizes associated with the ideas presented in this chapter. This gives us a feel for the chronology of these discoveries and indicates how important they were to the development of physics in the early 20th century.

Year | Recipient | Contribution |

1901 | W. K. Röntgen | Discovery of X-rays |

1906 | J. J. Thomson | Discovery of electron |

1914 | M. von Laue | X-ray diffraction in crystals |

1915 | W. and L. Bragg | X-ray analysis of crystal structure |

1918 | M. Planck | Energy quantization |

1921 | A. Einstein | Photoelectric effect |

1922 | N. Bohr | Structure of atoms |

1929 | L.-V. de Broglie | Wave nature of electrons |

1932 | W. Heisenberg | Quantum mechanics |

1933 | Schrödinger and Dirac | Atomic theory |

1937 | Davisson and Thomson | Electron diffraction in crystals |

____________________________

Figure 7.1 schematically illustrates interference between waves scattering from two adjacent rows of atoms in a crystal. The net effect of scattering from a single row is equivalent to partial reflection from a mirror imagined to be aligned with the row. Thus, the angle of “reflection” equals the angle of incidence for each row. Interference then occurs between the beams reflecting off different rows of atoms in the crystal.

For the two adjacent rows shown in figure 7.1, the path difference between beams is 2h = 2d sin θ. For constructive interference this must be an integer number of wavelengths, mλ, where the integer m is called the order of interference. The result is Bragg’s law of diffraction:

| (7.1) |

If only two rows are involved, the transition from constructive to destructive interference as θ changes is gradual. However, if interference from many rows occurs, then the constructive interference peaks become very sharp with mostly destructive interference in between. This sharpening of the peaks as the number of rows increases is very similar to the sharpening of the diffraction peaks from a diffraction grating as the number of slits increases.

Two types of targets are used in Bragg diffraction experiments: single crystals and powder targets.

In a single crystal setup, an X-ray detector is mounted as shown in figure 7.2. A mechanical device keeps the detector oriented so that the angle of incidence equals the angle of reflection for the desired crystal plane. Peaks in the X-ray detection rate are sought as the angle θ is varied.

The advantage of this type of apparatus is that diffraction peaks from only the selected crystal plane are observed.

The powder in a powder target is really a conglomeration of many tiny crystals randomly oriented. Thus, for each possible Bragg diffraction angle there are crystals oriented correctly for Bragg diffraction to take place. The detector is usually a photographic plate or an equivalent electronic device as illustrated in figure 7.3. For each Bragg diffraction angle one sees a ring on the plate concentric with the axis of the incident X-ray beam.

The advantage of this type of system is that no knowledge is needed of the crystal plane orientations. Furthermore, a single large crystal is not required. However, all possible Bragg scattering angles are seen at once, which can lead to confusion in the interpretation of the results.

Bragg diffraction illustrates the most difficult thing to understand about quantum mechanics, namely that particles can have wave-like properties and waves can have particle-like properties.

The variation of X-ray intensity with angle seen in a Bragg diffraction apparatus is very difficult to explain in any terms other than wave interference. Yet, X-rays are typically detected by a device such as a Geiger counter which produces a pulse of electricity for each X-ray particle, or photon, which hits it. Thus, X-rays sometimes act like particles and sometimes like waves.

Light isn’t alone in having both particle and wave properties. Davisson and Germer and later G. P. Thomson (son of J. J. Thomson, the discoverer of the electron) showed that electrons also can act like waves. They did this by demonstrating that electrons undergo Bragg diffraction in crystals, much as X-rays do.

Most physicists (including Albert Einstein) have found quantum mechanics to be extremely bizarre, so if you feel the same way, you are in good company! However, there is a useful interpretation of quantum mechanics which at least allows us to get on with using it to solve problems, even though it may not satisfy our intuitive reservations about the theory.

The displacement of the matter wave associated with a particle is usually
called the wave function, ψ. It is not at all clear what ψ is a displacement of,
but its use is straightforward. The absolute square of the wave function,
|ψ(x,t)|^{2}, is proportional to the probability of finding the associated particle
at position x and time t. The absolute square is taken because under
many circumstances the wave function is actually complex, i. e., it has
both real and imaginary parts. The reasons for this will be discussed
later.

Due to the interpretation of the wave function, quantum mechanics is a probabilistic theory. It does not tell us with certainty what happens to a particular particle. Instead, it tells us the probability for detecting the particle in any given location. If many experiments are done, with one particle per experiment, the numbers of experiments with particles being detected in the various possible locations are in proportion to the quantum mechanical probabilities.

Mechanics

The essential mystery of quantum mechanics becomes clearer when discussed in the conceptually simpler context of two slit interference. If light and electrons can have both particle and wave properties, then one might ask through which of the two slits the particle passed. However, in physics a question simply doesn’t make sense if it cannot be answered by experiment.

One can indeed perform an experiment to determine which slit the X-ray photon or electron passes through in the two slit experiment. However, by the very act of making this measurement, the form of the associated wave is altered. In particular, since the absolute square of the wave displacement represents the probability of finding the particle, once the particle has definitely been found passing through one or the other of the slits, the wave function collapses into a very small wave packet located at the observed position of the particle. Thus, the wave displacement becomes zero at the slit it didn’t pass through. However, the interference pattern results from the superposition of waves emanating from two slits. If no wave comes from one of the slits (because the wave displacement is zero there), then there can be no interference pattern!

We can now make the inverse argument. If there is an interference pattern, then we know that the wave displacement is non-zero at both slits. From the probability interpretation of the wave displacement, we conclude that we cannot say, even in principle, through which slit the particle passed. It is not just that we don’t know the answer to this question; there is simply no experiment which can give us an answer without destroying the interference pattern. In other words, the question “Which slit did the particle pass through?” is a nonsensical question in the case where an interference pattern is actually produced.

The American physicist Richard Feynman noticed that the above behavior can
be interpreted as violating the normal laws of probability. These laws
say that the probability of an event is the sum of the probabilities of
alternate independent ways for that event to occur. For instance, the
probability for a particle to reach point A on the detection screen of a
two slit setup is just the probability P_{1} for the particle to reach point
A after going through slit 1, plus the probability P_{2} for the particle to
reach point A after going through slit 2. Thus, if P_{1} = P_{2} = 0.1, then the
probability for the particle to reach point A irrespective of which slit
it went through should be P_{total} = P_{1} + P_{2} = 0.2. However, if point A
happens to be a point of destructive interference, then we know that
P_{total} = 0.

Feynman proposed that the above rule stating that alternate independent probabilities add, is simply incorrect. In its place Feynman asserted that probability amplitudes add instead, where the probability amplitude in this case is just the wave function associated with the particle. The total probability for a process is obtained by adding the alternate probability amplitudes together and taking the absolute square of the sum.

Feynman’s view is a particularly compact expression of the so-called Copenhagen interpretation of quantum mechanics which evolved from the ideas of Niels Bohr, Werner Heisenberg, Max Born, and others in the 1920s. It dispenses with the wave-particle duality and other philosophical baggage by saying “Particles are real in that we can observe them, but the only theory we have is about probability amplitudes for particles.” This interpretation of quantum mechanics may be weird, but it appears to be self-consistent and in agreement with experiment.

In this section we relate the classical ideas of mass, momentum and energy to what we have done so far. Historically, these connections were first made by Max Planck and Louis de Broglie with help from Albert Einstein. Bragg diffraction of electrons is invoked as an experimental test of the Planck and de Broglie relations.

Technically, we don’t need the ideas of mass, momentum, and energy to do physics – the notions of wavenumber, frequency, and group velocity are sufficient to describe and explain all observed phenomena. However, mass, momentum, and energy are so firmly embedded in physics that one couldn’t talk to other physicists without an understanding of these quantities!

Max Planck was the first to develop a theory explaining the energy density of electromagnetic radiation in a box at a fixed temperature. Albert Einstein extended Planck’s ideas by postulating that the energy of electromagnetic radiation is quantized into chunks called photons. The energy E of a photon is related to the frequency of the electromagnetic radiation by the equation

| (7.2) |

where f is the rotational frequency of the associated electromagnetic wave and ω
is its angular frequency. The constant h = 6.63 × 10^{-34} kg m^{2} s^{-1} is called
Planck’s constant. The related constant ℏ = h∕2π = 1.06 × 10^{-34} kg m^{2} s^{-1} is also
referred to as Planck’s constant, but to avoid confusion with the original constant,
we will generally refer to it as “h bar”.

Notice that a new physical dimension has appeared, namely mass, with the unit kilogram, abbreviated “kg”. The physical meaning of mass is much like our intuitive understanding of the concept, i. e., as a measure of the resistance of an object to its velocity being changed. The precise scientific meaning will emerge shortly.

Einstein showed that Planck’s idea could be used to explain the emission of
electrons which occurs when light impinges on the surface of a metal. This
emission, which is called the photoelectric effect, can only occur when electrons
are supplied with a certain minimum energy E_{B} required to break them loose
from the metal. Experiment shows that this emission occurs only when the
frequency of the light exceeds a certain minimum value. This value turns out to
equal ω_{min} = E_{B}∕ℏ, which suggests that electrons gain energy by absorbing
a single photon. If the photon energy, ℏω, exceeds E_{B}, then electrons
are emitted, otherwise they are not. It is much more difficult to explain
the photoelectric effect from the classical theory of light. The value of
E_{B}, called the binding energy or work function, is different for different
metals.

Louis de Broglie proposed that Planck’s energy-frequency relationship be extended to all kinds of particles. In addition he hypothesized that the momentum Π of the particle and the wave vector k of the corresponding wave were similarly related:

| (7.3) |

Note that this can also be written in scalar form in terms of the wavelength as Π = h∕λ. (We use Π rather than the more common p for momentum, because as we shall see, there are two different kinds of momentum, one related to the wavenumber, the other related to the velocity of a particle. In many cases they are equal, but there are certain important situations in which they are not.)

De Broglie’s hypothesis was inspired by the fact that wave frequency and wavenumber are components of the same four-vector according to the theory of relativity, and are therefore closely related to each other. Thus, if the energy of a particle is related to the frequency of the corresponding wave, then there ought to be some similar quantity which is correspondingly related to the wavenumber. It turns out that the momentum is the appropriate quantity. The physical meaning of momentum will become clear as we proceed.

We will also find that the rest frequency, μ, of a particle is related to its mass, m:

| (7.4) |

The quantity E_{rest} is called the rest energy of the particle.

From our perspective, energy, momentum, and rest energy are just scaled versions of frequency, wave vector, and rest frequency, with a scaling factor ℏ. We can therefore define a four-momentum as a scaled version of the wave four-vector:

| (7.5) |

The spacelike component of Π is just Π, while the timelike part is E∕c.

Planck, Einstein, and de Broglie had extensive backgrounds in classical mechanics, in which the concepts of energy, momentum, and mass have precise meaning. In this text we do not presuppose such a background. Perhaps the best strategy at this point is to think of these quantities as scaled versions of frequency, wavenumber, and rest frequency, where the scale factor is ℏ. The significance of these quantities to classical mechanics will emerge bit by bit.

Let us now recapitulate what we know about relativistic waves, and how this knowledge translates into knowledge about the mass, energy, and momentum of particles. In the following equations, the form on the left is expressed in wave terms, i. e., in terms of frequency, wavenumber, and rest frequency. The form on the right is the identical equation expressed in terms of energy, momentum, and mass. Since the latter variables are just scaled forms of the former, the two forms of each equation are equivalent.

We begin with the dispersion relation for relativistic waves:

| (7.6) |

Calculation of the group velocity, u_{g} = dω∕dk, from the dispersion relation
yields

| (7.7) |

These two sets of equations represent what we know about relativistic waves, and what this knowledge tells us about the relationships between the mass, energy, and momentum of relativistic particles. When in doubt, refer back to these equations, as they work in all cases, including for particles with zero mass!

It is useful to turn equations (7.6) and (7.7) around so as to express the frequency as a function of rest frequency and group velocity,

| (7.8) |

and the wavenumber as a similar function of these quantities:

| (7.9) |

Note that equations (7.8) and (7.9) work only for particles with non-zero mass! For zero mass particles both the numerators and denominators of equations (7.8) and (7.9) are zero, making these equations undefined, and you need to use equations (7.6) and (7.7) with m = 0 and μ = 0 instead.

The quantity ω - μ indicates how much the frequency exceeds the rest
frequency. Notice that if ω = μ, then from equation (7.6) k = 0. Thus,
positive values of ω_{k} ≡ ω - μ indicate |k| > 0, which means that the
particle is moving according to equation (7.7). Let us call ω_{k} the kinetic
frequency:

| (7.10) |

We call K the kinetic energy for similar reasons. Again, equation (7.10) only works for particles with non-zero mass. For zero mass particles the kinetic energy equals the total energy.

Note that the results of this section are valid only for free particles, i. e., particles to which no force is applied. Force in classical and quantum mechanics is treated in the next chapter.

When the mass is non-zero and the group velocity is much less than the speed of
light, it is useful to compute approximate forms of the above equations valid in
this limit. Using the approximation (1 + ϵ)^{x} ≈ 1 + xϵ, we find that the dispersion
relation becomes

| (7.11) |

and the group velocity equation takes the approximate form

| (7.12) |

The non-relativistic limits for equations (7.8) and (7.9) become

| (7.13) |

| (7.14) |

while the approximate kinetic energy equation is

| (7.15) |

Just a reminder — the equations in this section are not valid for massless particles!

How can we test the above predictions against experiment? The key point is to be
able to relate the wave aspects to the particle aspects of a quantum mechanical
wave-particle. Equation (7.9), or equation (7.14) in the non-relativistic case,
relates a particle’s wavenumber k to its velocity u_{g}. Both of these quantities can
be measured in a Bragg’s law experiment with electrons. In this experiment
electrons are fired at a crystal with known atomic dimensions at a known speed,
which we identify with the group velocity u_{g}. The Bragg angle which yields
constructive interference can be used to calculate the wavelength of the
corresponding electron wave, and hence the wavenumber and momentum. If
the momentum is plotted against group velocity in the non-relativistic
case, a straight line should be found, the slope of which is the particle’s
mass. In the fully relativistic case one needs to plot momentum versus
u_{g}∕(1 -u_{g}^{2}∕c^{2})^{1∕2}. Again, a straight line indicates agreement with the theory and
the slope of the line is the particle’s mass. This particular experiment
is difficult to do, but the corresponding theories verify in many other
experiments.

Classically, we consider the location of a particle to be a knowable piece of information. In quantum mechanics the position of a particle is well known if the wave packet representing it is small in size. However, quantum mechanics imposes a price on accurately knowing the position of a particle in terms of the future predictability of its position. This is because a small wave packet, which corresponds to accurate knowledge of the corresponding particle’s position, implies the superposition of plane waves corresponding to a broad distribution of wavenumbers. This translates into a large uncertainty in the wavenumber, and hence the momentum of the particle. In contrast, a broad wave packet corresponds to a narrower distribution of wavenumbers, and correspondingly less uncertainty in the momentum.

Referring back to chapters 1 and 2, recall that both the longitudinal (along the
direction of motion) and transverse (normal to the direction of motion)
dimensions of a wave packet, Δx_{L} and Δx_{T }, can be related to the spread of
longitudinal and transverse wavenumbers, Δk_{L} and Δk_{T }:

| (7.16) |

| (7.17) |

We have omitted numerical constants which are of order unity in these approximate relations so as to show their essential similarity.

The above equations can be interpreted in the following way. Since the
absolute square of the wave function represents the probability of finding a
particle, Δx_{L} and Δx_{T } represent the uncertainty in the particle’s position.
Similarly, Δk_{L} and Δk_{T } represent the uncertainty in the particle’s longitudinal
and transverse wave vector components. This latter uncertainty leads to
uncertainty in the particle’s future motion — larger or smaller longitudinal k
results respectively in larger or smaller particle speed, while uncertainty in the
transverse wavenumber results in uncertainty in the particle’s direction of motion.
Thus uncertainties in any component of k result in uncertainties in the
corresponding component of the particle’s velocity, and hence in its future
position.

The equations (7.16) and (7.17) show that uncertainty in the present and future positions of a particle are complimentary. If the present position is accurately known due to the small size of the associated wave packet, then the future position is not very predictable, because the wave packet disperses rapidly. On the other hand, a broad-scale initial wave packet means that the present position is poorly known, but the uncertainty in position, poor as it is, doesn’t rapidly increase with time, since the wave packet has a small uncertainty in wave vector and thus disperses slowly. This is a statement of the Heisenberg uncertainty principle.

The uncertainty principle also applies between frequency and time:

| (7.18) |

This shows up in the beat frequency equation 1∕T_{beat} = Δf = Δω∕2π. The beat
period T_{beat} may be thought of as the size of a “wave packet in time”. The beat
frequency equation may be rewritten as ΔωT_{beat} = 2π, which is the same as
equation (7.18) if the factor of 2π is ignored and T_{beat} is identified with
Δt.

The above forms of the uncertainty principle are not relativistically invariant. A useful invariant form may be obtained by transforming to the coordinate system in which a particle is stationary. In this reference frame the time t becomes the proper time τ associated with the particle. Furthermore, the frequency ω becomes the rest frequency μ. The uncertainty principle thus becomes

| (7.19) |

in this reference frame. However, since Δμ and Δτ are relativistic invariants, this expression of the uncertainty principle is valid in any reference frame.

It is more common to express the uncertainty principle in terms of the mass, momentum, and energy by multiplying equations (7.16) - (7.19) by ℏ. Lumping the momentum equations, we find

| (7.20) |

| (7.21) |

| (7.22) |

Classical mechanics is the realm of quantum mechanics in which the dimensions of the system of interest are much larger than the wavelengths of the waves corresponding to the particles constituting the system. In this case the uncertainties induced by the uncertainty principle are unimportant. This limit is analogous to the geometrical optics limit for light. Thus, we can say that classical mechanics is the geometrical optics limit of quantum mechanics.

- An electron with wavelength λ = 1.2 × 10
^{-10}m undergoes Bragg diffraction from a single crystal with atomic plane spacing of d = 2 × 10^{-10}m.- Calculate the Bragg angles (all of them!) for which constructive interference occurs.
- Calculate the speed of the electron.

- Suppose that electrons impinge on two slits in a plate, resulting in a two slit
diffraction pattern on a screen on the other side of the plate. The
probability for an electron to pass through either one of the slits and
reach point A on the screen is P, assuming that the other slit is
blocked.
- If there are two slits open and A is a point of constructive interference, what is the probability of an electron reaching A? Hint: Remember that amplitudes, not probabilities add.
- If there are two slits open and A is a point of destructive interference, what is the probability of an electron reaching A?
- If there are two slits open, what is the probability for an electron to reach point A according to the conventional rule that probabilities add? (This is the result one would expect if, for instance, the particles were machine gun bullets and the slits were, say, 5 cm apart.)
- If the slit separation is very much greater than the electron wavelength, how does this affect the spacing of regions of constructive and destructive interference? Explain how the results of parts (a) and (b) become approximately consistent with those of part (c) in this case.

- Compute the (angular) rest frequency of an electron and a neutron. (Look up their masses.)
- How does the dispersion relation for relativistic waves simplify if the rest frequency (and hence the particle mass) is zero? What is the group velocity in this case?
- X-rays are photons with frequencies about 2000 times the frequencies of ordinary light photons. From this information and what you know about light, infer the approximate velocity of electrons which have Bragg diffraction properties similar to X-rays. Are the electrons relativistic or non-relativistic?
- Electrons with velocity v = 0.6c are diffracted with a 0.2 radian half-angle of diffraction when they hit an object. What is the approximate size of the object? Hint: Diffraction of a wave by an object of a certain size is quite similar to diffraction by a hole in a screen of the same size.
- Work out an approximate formula for the kinetic energy of a particle as a
function of mass m and velocity u
_{g}which is valid when u_{g}^{2}≪ c^{2}. Hint: Use the approximation (1 + ϵ)^{x}≈ 1 + xϵ, which is valid for |ϵ|≪ 1. As u_{ g}∕c becomes larger, how does this approximate formula deviate from the exact formula? - Work out an approximate formula for the momentum of a particle as a
function of m and u
_{g}in the case where u_{g}^{2}≪ c^{2}. You may wish to use the approximation mentioned in the previous problem. - If a photon is localized to within a distance Δx, what is the uncertainty in the photon energy?
- If an electron is localized to within a distance Δx, what is the uncertainty in
the electron kinetic energy? Hint: As long as ΔΠ ≪ Π, ΔΠ
^{2}≈ 2ΠΔΠ. To see why, compute dΠ^{2}∕dΠ. - A grocer dumps some pinto beans onto a scale, estimates their mass as 2 kg,
and then dumps them off after 5 s. What is the quantum mechanical
uncertainty in this measurement? Assume this occurs in Quantum World
where the speed of light is 10 m s
^{-1}(speed of a fast buggy) and Planck’s constant ℏ = 1 kg m^{2}s^{-1}. - Mary’s physics text (mass 0.3 kg) has to be kept on a leash (length 0.5 m)
to prevent it from wandering away from her in Quantum World
(ℏ = 1 kg m
^{2}s^{-1}).- If the leash suddenly breaks, what is the maximum speed at which the book is likely to move away from its initial location?
- In order to reduce this speed, should Mary make the new leash shorter or longer than the old one? Explain.

- A proton (mass M = 1.7 × 10
^{-27}kg) is confined to an atomic nucleus of diameter D = 2 × 10^{-15}m.- What is the uncertainty in the proton’s momentum?
- Roughly what kinetic energy might you expect the proton to have?

Planck’s constant is ℏ = 1.06 × 10

^{-34}kg m^{2}s^{-1}. You may use the non-relativistic equation for the energy.

Geometrical Optics and Newton’s Laws

The question that motivates us to study physics is “What makes things go?” The answers we conceive to this question constitute the subject of dynamics. This is in contrast to the question we have primarily addressed so far, namely “How do things go?” The latter question is about kinematics. Extensive preparation in the kinematics of waves and particles in relativistic spacetime is needed to intelligently address dynamics. This preparation is now complete.

In this chapter we outline three different dynamical principles based respectively on pre-Newtonian, Newtonian, and quantum mechanical thinking. We first discuss the Newtonian mechanics of conservative forces in one dimension. Certain ancillary concepts in mechanics such as work and power are introduced at this stage. We then show that Newtonian and quantum mechanics are consistent with each other in the realm in which they overlap, i. e., in the geometrical optics limit of quantum mechanics. For simplicity, this relationship is first developed in one dimension in the non-relativistic limit. Higher dimensions require the introduction of partial derivatives, and the relativistic case will be considered later.

Roughly speaking, there have been three eras of physics, characterized by three different answers to the question of what makes things go.

Aristotle expounded a view of dynamics which agrees closely with our everyday experience of the world. Objects only move when a force is exerted upon them. As soon as the force goes away, the object stops moving. The act of pushing a box across the floor illustrates this principle — the box certainly doesn’t move by itself!

In contrast to earthly behavior, the motions of celestial objects seem effortless. No obvious forces act to keep the planets in motion around the sun. In fact, it appears that celestial objects simply coast along at constant velocity unless something acts on them. The Newtonian view of dynamics — objects change their velocity rather than their position when a force is exerted on them — is expressed by Newton’s second law:

| (8.1) |

where F is the force exerted on a body, m is its mass, and a is its acceleration. Newton’s first law, which states that an object remains at rest or in uniform motion unless a force acts on it, is actually a special case of Newton’s second law which applies when F = 0.

It is no wonder that the first successes of Newtonian mechanics were in the celestial realm, namely in the predictions of planetary orbits. It took Newton’s genius to realize that the same principles which guided the planets also applied to the earthly realm as well. In the Newtonian view, the tendency of objects to stop when we stop pushing on them is simply a consequence of frictional forces opposing the motion. Friction, which is so important on the earth, is negligible for planetary motions, which is why Newtonian dynamics is more obviously valid for celestial bodies.

Note that the principle of relativity is closely related to Newtonian physics and is incompatible with pre-Newtonian views. After all, two reference frames moving relative to each other cannot be equivalent in the pre-Newtonian view, because objects with nothing pushing on them can only come to rest in one of the two reference frames. Newton’s second law obeys the principle of relativity because the acceleration of an object is the same when viewed from two different reference frames moving at a constant velocity with respect to each other.

Einstein’s relativity is often viewed as a repudiation of Newton, but this is far from the truth — Newtonian physics makes the theory of relativity possible through its invention of the principle of relativity. Compared with the differences between pre-Newtonian and Newtonian dynamics, the changes needed to go from Newtonian to Einsteinian physics constitute minor tinkering.

In quantum mechanics, particles are represented by matter waves, with the absolute square of the wave displacement yielding the probability of finding the particle. The behavior of particles thus follows from the reflection, refraction, diffraction, and interference of the associated waves. The connection with Newtonian dynamics comes from tracing the trajectories of matter wave packets. Changes in the speed and direction of motion of these packets correspond to the accelerations of classical mechanics. When wavelengths are small compared to the natural length scale of the problem at hand, the wave packets can be made small, thus pinpointing the position of the associated particle, without generating excessive uncertainty in the particle’s momentum. This is the geometrical optics limit of quantum mechanics.

We now address Newtonian mechanics in the case where the force on a particle is conservative. A conservative force is one that can be derived from a so-called potential energy U. We assume that the potential energy of the particle depends only on its position. The force is obtained from the potential energy by the equation

| (8.2) |

Using this equation we write Newton’s second law as

| (8.3) |

We then notice that the acceleration can be written in terms of the x derivative
along the particle’s trajectory of v^{2}∕2:

| (8.4) |

The last step in the above derivation can be verified by applying the product rule:
dv^{2}∕dt = d(vv)∕dt = v(dv∕dt) + (dv∕dt)v = 2v(dv∕dt). Putting equations (8.3)
and (8.4) together, we find that d(mv^{2}∕2 + U)∕dt = 0, which implies that
mv^{2}∕2 + U is constant. We call this constant the total energy E and the quantity
K = mv^{2}∕2 the kinetic energy. We thus have the principle of conservation of
energy for conservative forces:

| (8.5) |

Recall that in quantum mechanics the momentum is related to the group
velocity u_{g} by

| (8.6) |

in the nonrelativistic case. Equating the group velocity to v and eliminating it in the kinetic energy results in an alternate expression for this quantity:

| (8.7) |

Since the total energy E is constant or conserved, increases in the potential energy coincide with decreases in the kinetic energy and vice versa, as is illustrated in figure 8.1. In classical mechanics the kinetic energy cannot be negative, since it is the product of half the mass and the square of the velocity, both of which are positive. Thus, a particle with total energy E and potential energy U is forbidden to venture into regions in which the kinetic energy K = E - U is less than zero.

The points at which the kinetic energy is zero are called turning points. This is because a particle decreases in speed as it approaches a turning point, stops there for an instant, and reverses direction. Note also that a particle with a given total energy always has the same speed at some point x regardless of whether it approaches this point from the left or the right:

| (8.8) |

An example of a conservative force is gravity. An object of mass m near the surface of the earth has the gravitational potential energy

| (8.9) |

where z is the height of the object above some reference point such as the earth’s
surface and g = 9.8 m s^{-2} is the local value of the gravitational field near the
surface. Notice that the gravitational potential energy increases upward. The
speed of the object in this case is |u_{g}| = [2(E - mgz)∕m]^{1∕2}. If |u_{
g}| is known to
equal the constant value u_{0} at elevation z = 0, then equations (8.8) and (8.9) tell
us that u_{0} = (2E∕m)^{1∕2} and |u_{
g}| = (u_{0}^{2} - 2gz)^{1∕2}.

There are certain types of questions which energy conservation cannot directly
answer. For instance, if an object is released at elevation h with zero velocity at
t = 0, at what time will it reach z = 0 under the influence of gravity? In
such cases it is often easiest to return to Newton’s second law. Since the
force on the object is F = -dU∕dz = -mg in this case, we find that the
acceleration is a = F∕m = -mg∕m = -g. However, a = du∕dt = d^{2}z∕dt^{2},
so

| (8.10) |

where C_{1} and C_{2} are constants to be determined by the initial conditions. These
results can be verified by differentiating to see if the original acceleration is
recovered. Since u = 0 and z = h at t = 0, we have C_{1} = 0 and C_{2} = h.
With these results it is easy to show that the object reaches z = 0 when
t = (2h∕g)^{1∕2}.

When a force is exerted on an object, energy is transferred to the object. The amount of energy transferred is called the work done on the object. However, energy is only transferred if the object moves. The work W done is

| (8.11) |

where the distance moved by the object is Δx and the force exerted on it is F. Notice that work can either be positive or negative. The work is positive if the object being acted upon moves in the same direction as the force, with negative work occurring if the object moves opposite to the force.

Equation (8.11) assumes that the force remains constant over the full displacement Δx. If it is not, then it is necessary to break up the displacement into a number of smaller displacements, over each of which the force can be assumed to be constant. The total work is then the sum of the works associated with each small displacement.

If more than one force acts on an object, the works due to the different forces each add or subtract energy, depending on whether they are positive or negative. The total work is the sum of these individual works.

There are two special cases in which the work done on an object is related to other quantities. If F is the total force acting on the object, then W = FΔx = maΔx by Newton’s second law. However, a = dv∕dt where v is the velocity of the object, and Δx = (Δx∕Δt)Δt ≈ vΔt, where Δt is the time required by the object to move through distance Δx. The approximation becomes exact when Δx and Δt become very small. Putting all of this together results in

| (8.12) |

where K is the kinetic energy of the object. Thus, when F is the only
force, W = W_{total} is the total work on the object, and this equals the
change in kinetic energy of the object. This is called the work-energy
theorem, and it demonstrates that work really is a transfer of energy to an
object.

The other special case occurs when the force is conservative, but is not necessarily the total force acting on the object. In this case

| (8.13) |

where ΔU is the change in the potential energy of the object associated with the force of interest.

The power associated with a force is simply the amount of work done by the force divided by the time interval Δt over which it is done. It is therefore the energy per unit time transferred to the object by the force of interest. From equation (8.11) we see that the power is

| (8.14) |

where v is the velocity at which the object is moving. The total power is just the sum of the powers associated with each force. It equals the time rate of change of kinetic energy of the object:

| (8.15) |

Louis de Broglie^{1}
made an analogy between matter waves and light waves, pointing out that wave
packets of light change their velocity as the result of spatial variations in
the index of refraction of the medium in which they are travelling. This
behavior comes about because the dispersion relation for light traveling
through a medium with index of refraction n is ω = kc∕n, so that the group
velocity, u_{g} = dω∕dk = c∕n. Thus, when n increases, u_{g} decreases, and vice
versa.^{2}

In this section we pursue de Broglie’s analogy to see if we can come up
with a theory of matter waves which gives the same results as classical
mechanics in the geometrical optics limit of these waves. The dispersion
relation for free matter waves is ω = (k^{2}c^{2} + μ^{2})^{1∕2}. In the non-relativistic
limit k^{2}c^{2} ≪ μ^{2}. As done previously, we use (1 + ϵ)^{n} ≈ 1 + nϵ for small ϵ.
In the non-relativistic limit, the dispersion relation for free waves thus
becomes

| (8.16) |

The above equation can be transformed into the total energy equation for a
free, non-relativistic particle, E = mc^{2} + K, where mc^{2} is the rest energy and K is
the kinetic energy, by multiplying by ℏ. We convert the free particle energy
equation into the equation for a particle subject to a conservative force by adding
the potential energy U the right side. The analogous change to equation (8.16)
is to add S = U∕ℏ to the right side, resulting in a modified dispersion
relation:

| (8.17) |

(Since the rest energy is just a constant, we have absorbed it into S.) This gives us the dispersion relation for one-dimensional matter waves subject to a spatially varying potential energy. The quantity S, which we see is just a scaled potential energy, plays a role for matter waves which is analogous to the role played by a spatially variable index of refraction for light waves.

Let us now imagine that all parts of the wave governed by this dispersion relation oscillate in phase. The only way this can happen is if ω is constant, i. e., it takes on the same value in all parts of the wave.

If ω is constant, the only way S can vary with x in equation (8.17) is if the wavenumber varies in a compensating way. Thus, constant frequency and spatially varying S together imply that k = k(x). Solving equation (8.17) for k yields

| (8.18) |

Since ω is constant, the wavenumber becomes smaller and the wavelength larger as the wave moves into a region of increased S.

In the geometrical optics limit, we assume that S doesn’t change much over one wavelength so that the wave remains reasonably sinusoidal in shape with approximately constant wavenumber over a few wavelengths. However, over distances of many wavelengths the wavenumber and amplitude of the wave are allowed to vary considerably.

The group velocity calculated from the dispersion relation given by equation (8.17) is

| (8.19) |

where k is eliminated in the last step with the help of equation (8.18). The
resulting equation tells us how the group velocity varies as a matter wave
traverses a region of slowly varying S. Thus, as S increases, u_{g} decreases and vice
versa.

We can now calculate the acceleration of a wave packet resulting from the spatial
variation in S. We assume that x(t) represents the position of the wave packet, so
that u_{g} = dx∕dt. Using the chain rule du_{g}∕dt = (du_{g}∕dx)(dx∕dt) = (du_{g}∕dx)u_{g}, we
find

| (8.20) |

The group velocity is eliminated in favor of S by squaring equation (8.19) and substituting the result into equation (8.20).

Recalling that U = ℏS, equation (8.20) becomes

| (8.21) |

which is just Newton’s second law! Thus, the geometrical optics approach to particle motion is completely equivalent to the classical mechanics of a particle moving under the influence of a conservative force, at least in one dimension. We therefore have two ways of solving for the motion of a particle subject to a potential energy U(x). We can apply the principles of classical mechanics to get the force and the acceleration of the particle, from which we can derive the motion. Alternatively, we can apply the principles of geometrical optics to compute the spatially variable velocity of the wave packet using equation (8.19). The results are completely equivalent, though the methods are conceptually very different.

In order to understand the generalization of Newtonian mechanics to two and three dimensions, we first need to understand a new type of derivative called the partial derivative. The partial derivative is used in functions of more than one variable. It is just like an ordinary derivative, except that when taking the derivative of the function with respect to one of the variables, the other variables are held constant. As an example, let us consider the function

| (8.22) |

where A, B, and C are constants. The partial derivative of f with respect to x is

| (8.23) |

and the partial derivative with respect to y is

| (8.24) |

That’s it! Note that a special symbol “∂” is used in place of the normal “d” for the partial derivative. This is sometimes called a “curly d”.

When a matter wave moves through a region of variable potential energy in one dimension, only the wavenumber changes. In two or three dimensions, the wave vector can change in both direction and magnitude. This complicates the calculation of particle movement. However, we already have an example of how to handle this situation, namely, the refraction of light. In that case Snell’s law tells us how the direction of the wave vector changes, while the dispersion relation combined with the constancy of the frequency gives us information about the change in the magnitude of the wave vector. For matter waves a similar procedure works, though the details are different, because we seek the consequences of a change in potential energy rather than a change in the index of refraction.

Figure 8.2 illustrates the refraction of matter waves at a discontinuity in the
potential energy. Let us suppose that the discontinuity occurs at x = 0. If the
matter wave to the left of the discontinuity is ψ_{1} = sin(k_{1x}x + k_{1y}y - ω_{1}t) and to
the right is ψ_{2} = sin(k_{2x}x + k_{2y}y - ω_{2}t), then the wavefronts of the waves
will match across the discontinuity for all time only if ω_{1} = ω_{2} ≡ ω and
k_{1y} = k_{2y} ≡ k_{y}. We are already familiar with the first condition from the
one-dimensional problem, so the only new ingredient is the constancy of the y
component of the wave vector.

In two dimensions the momentum is a vector: Π = mu when |u|≪ c,
where u is the particle velocity. Furthermore, the kinetic energy is
K = m|u|^{2}∕2 = |Π|^{2}∕(2m) = (Π_{
x}^{2} + Π_{
y}^{2})∕(2m). The relationship between kinetic,
potential, and total energy is unchanged from the one-dimensional case, so we
have

| (8.25) |

The de Broglie relationship tells us that Π = ℏk, so the constancy of k_{y} across the
discontinuity in U tells us that

| (8.26) |

there.

Let us now approximate a continuously variable U(x) by a series of steps of constant U oriented normal to the x axis. The above analysis can be applied at the jumps or discontinuities in U between steps, as illustrated in figure 8.3, with the result that equations (8.25) and (8.26) are valid across all discontinuities. If we now let the step width go to zero, these equations then become valid for U continuously variable in x.

An example from classical mechanics of a problem of this type is a ball rolling down an inclined ramp with an initial velocity component across the ramp, as illustrated in figure 8.4. The potential energy decreases in the down ramp direction, resulting in a force down the ramp. This accelerates the ball in that direction, but leaves the component of momentum across the ramp unchanged.

Using the procedure which we invoked before, we find the force components
associated with U(x) in the x and y directions to be F_{x} = -dU∕dx and F_{y} = 0.
This generalizes to

| (8.27) |

in the three-dimensional case where the orientation of constant U surfaces is arbitrary. It is also valid when U(x,y,z) is not limited to a simple ramp form, but takes on a completely arbitrary structure.

The definitions of work and power are slightly different in two and three dimensions. In particular, work is defined as

| (8.28) |

where Δx is now a vector displacement. The vector character of this expression yields an additional possibility over the one dimensional case, where the work is either positive or negative depending on the direction of Δx relative to F. If the force and the displacement of the object on which the force is acting are perpendicular to each other, the work done by the force is actually zero, even though the force and the displacement both have non-zero magnitudes. The power exhibits a similar change:

| (8.29) |

Thus the power is zero if an object’s velocity is normal to the force being exerted on it.

As in the one-dimensional case, the total work done on a particle equals the change in the particle’s kinetic energy. In addition, the work done by a conservative force equals minus the change in the associated potential energy.

Energy conservation by itself is somewhat less useful for solving problems in
two and three dimensions than it is in one dimension. This is because
knowing the kinetic energy at some point tells us only the magnitude of the
velocity, not its direction. If conservation of energy fails to give us the
information we need, then we must revert to Newton’s second law, as we
did in the one-dimensional case. For instance, if an object of mass m
has initial velocity u_{0} = (u_{0}, 0) at location (x,z) = (0,h) and has the
gravitational potential energy U = mgz, then the force on the object is
F = (0,-mg). The acceleration is therefore a = F∕m = (0,-g). Since
a = du∕dt = d^{2}x∕dt^{2} where x = (x,z) is the object’s position, we find
that

| (8.30) |

where C_{1}, C_{2}, C_{3}, and C_{4} are constants to be evaluated so that the solution
reduces to the initial conditions at t = 0. The specified initial conditions
tell us that C_{1} = u_{0}, C_{2} = 0, C_{3} = 0, and C_{4} = h in this case. From
these results we can infer the position and velocity of the object at any
time.

If you have previously taken a physics course then you have probably
noticed that a rather odd symbol is used for momentum, namely Π, rather
than the more commonly employed p. The reason for this peculiar usage
is that there are actually two kinds of momentum, kinetic momentum
and total momentum, just as there are two kinds of energy, kinetic and
total.^{3}
The symbol Π represents total momentum while p represents kinetic momentum.
Normally we don’t need to distinguish between the two quantities, as they are
generally equal to each other. However, we will find later in the course that it is
crucial to make this distinction in the case of charged particles in a magnetic field.
As a general rule, the total momentum is related to a particle’s wave vector via
the de Broglie relation, Π = ℏk, while the kinetic momentum is related to a
particle’s velocity, p = mu∕(1 - u^{2}∕c^{2})^{1∕2}.

- Suppose the dispersion relation for a matter wave under certain
conditions is ω = μ + (k - a)
^{2}c^{2}∕(2μ) where k is the wavenumber of the wave, μ = mc^{2}∕ℏ, m is the associated particle’s mass, a is a constant, c is the speed of light, and ℏ is Planck’s constant divided by 2π.- Use this disperson relation and the Planck and de Broglie relations to determine the relationship between energy E, momentum Π, and mass m.
- Compute the group velocity of the wave and use this to determine how the group velocity depends on mass and momentum in this case.

- A matter wave function associated with a particle of definite (constant)
total energy E takes the form shown in figure 8.5. Make a sketch showing
how the kinetic, potential, and total energies of the particle vary with
x.

- Compute ∂∕∂x and ∂∕∂y of the following functions. Other symbols are
constants.
- f(x,y) = ax
^{2}+ by^{3} - f(x,y) = ax
^{2}y^{2} - f(x,y) = (x + a)∕(y + b)

- f(x,y) = ax
- Given a potential energy for a particle of mass M of the form U(x) = Ax
^{3}-Bx where A and B are positive constants:- Find the force on the particle.
- Find the values of x where the force is zero.
- Sketch U(x) versus x and graphically compare the slope of U(x) to the force computed above. Do the two qualitatively match?
- If the total energy of the particle is zero, where are its turning points?
- What is the particle’s speed as a function of position assuming that the total energy E is known?

- Given a potential energy function U(x,y) = A(x
^{2}+ y^{2}) where A is a positive constant:- Sketch lines of constant U in the x-y plane.
- Compute the components of force as a function of x and y and draw sample force vectors in the x-y plane on the same plot used above. Do the force vectors point “uphill” or “downhill”?

- Do the same as in the previous question for the potential energy function U(x,y) = Axy.
- Suppose that the components of the force vector in the x-y plane are
F = (2Axy
^{3}, 3Ax^{2}y^{2}) where A is a constant. See if you can find a potential energy function U(x,y) which gives rise to this force. - You are standing on top of a cliff of height H with a rock of mass
M.
- If you throw the rock horizontally outward at speed u
_{0}, what will its speed be when it hits the ground below? - If you throw the rock upward at 45
^{∘}to the horizontal at speed u_{ 0}, what will its speed be when it hits the ground?

Hint: Can you use conservation of energy to solve this problem? Ignore air friction.

- If you throw the rock horizontally outward at speed u
- A car of mass 1200 kg initially moving 30 m s
^{-1}brakes to a stop.- What is the net work done on the car due to all the forces acting on it during the indicated period?
- Describe the motion of the car relative to an inertial reference frame initially moving with the car.
- In the above reference frame, what is the net work done on the car during the indicated period?

Is work a relativistically invariant quantity?

- A soccer player kicks a soccer ball, which is caught by the goal keeper as
shown in figure 8.6. At various points forces exerted by gravity, air friction,
the foot of the offensive player, and the hands of the goal keeper act on the
ball.
- List the forces acting on the soccer ball at each of the points A, B, C, D, and E.
- State whether the instantaneous power being applied to the soccer ball due to each of the forces listed above is positive, negative, or zero at each of the labeled points.

- A cannon located at (x,z) = (0, 0) shoots a cannon ball upward at an angle
of θ from the horizontal at initial speed u
_{0}. Hint: In order to solve this problem you must first obtain the x and z components of acceleration from Newton’s second law. Second, you must find the velocity components as a function of time from the components of acceleration. Third, you must find x and z as a function of time from the the components of velocity. Only then should you attempt to answer the questions below.- How long does it take the cannon ball to reach its peak altitude?
- How high does the cannon ball go?
- At what value of x does the cannon ball hit the ground (z = 0)?
- Determine what value of θ yields the maximum range.

Symmetry and Bound States

When quantum mechanics was first invented, the dynamical principles used were the same as those underlying classical mechanics. The initial development of the field thus proceeded largely by imposing quantum laws on classical variables such as position, momentum, and energy. However, as quantum mechanics advanced, it became clear that there were many situations in which no classical analogs existed for new types of quantum mechanical systems, especially those which arose in the study of elementary particles. To understand these systems it was necessary to seek guidance from novel sources. One of the most important of these sources was the idea of symmetry, and in particular the relationship between symmetry and conserved variables. This type of relationship was first developed in the early 20th century by the German mathematician Emmy Nöther in the context of classical mechanics. However, her idea is easier to express and use in quantum mechanics than it is in classical mechanics. Emmy Nöther showed that there is a relationship between the symmetries of a system and conserved dynamical variables. This idea is naturally called Nöther’s theorem.

In classical mechanics a particle is bound or confined to a particular region if its total energy exceeds the potential energy only in this region. In quantum mechanics, matter waves can also be confined for the same reason, though the confinement is often less perfect than in the classical case. In this chapter we examine the consequences of wave confinement. We first look at the so-called “particle in a box” in one spatial dimension. We find that confined particles can take on only discrete energy values. When confinement isn’t perfect we see how a quantum mechanical particle can leak through a potential energy barrier which is classically impenetrable. Movement of a particle on a circular ring leads us to another form of confinement and the introduction of angular momentum. This brings us finally to a discussion of the intrinsic or spin angular momentum of elementary particles.

Until now we have represented quantum mechanical plane waves by sine and cosine functions, just as with other types of waves. However, plane matter waves cannot be truly represented by sines and cosines. We need instead mathematical functions in which the wave displacement is complex rather than real. This requires the introduction of a bit of new mathematics, which we tackle first. Using our new mathematical tool, we are then able to explore two crucially important ideas in quantum mechanics; (1) the relationship between symmetry and conservation laws, and (2) the dynamics of spatially confined waves.

A complex number z is the sum of a real number and an imaginary number.
An imaginary number is just a real number multiplied by i ≡ (-1)^{1∕2}. Thus, we
can write z = a + ib for any complex z, where a and b are real. The quantities a
and b are the real and imaginary parts of z, sometimes written Re(z) and
Im(z).

Quantum mechanics requires wave functions to be complex, i. e., to possess real and imaginary parts. Plane waves in quantum mechanics actually take the form ψ = exp[i(kx - ωt)] rather than, say, cos(kx - ωt). The reason for this is the need to distinguish between waves with positive and negative frequencies. If we replace k and ω with -k and -ω in the cosine form, we get cos(-kx + ωt) = cos[-(kx - ωt)] = cos(kx - ωt). In other words, changing the sign of k and ω results in no change in a wave expressed as a cosine function. The two quantum mechanical states, one with wavenumber and frequency k and ω and the other with -k and -ω, yield indistinguishable wave functions and therefore would represent physically indistinguishable states. The cosine form is thus insufficiently flexible to represent quantum mechanical waves. On the other hand, if we replace k and ω with their negatives in the complex exponential form of a plane wave we get ψ = exp[-i(kx - ωt)], which is different from exp[i(kx - ωt)]. These two wave functions are distinguishable and thus correspond to distinct physical states.

It is not immediately obvious that a complex exponential function provides the oscillatory behavior needed to represent a plane wave. However, the complex exponential can be expressed in terms of sines and cosines using Euler’s equation:

| (9.1) |

If we define r = (a^{2} + b^{2})^{1∕2} and ϕ = tan ^{-1}(b∕a), then an alternate way of
expressing a complex number is z = r exp(iϕ), which by Euler’s equation equals
r cos(ϕ) + ir sin(ϕ). Comparison shows that a = r cos(ϕ) and b = r sin(ϕ). Thus, a
complex number can be thought of as a point in the a-b plane with Cartesian
coordinates a and b and polar coordinates r and ϕ. The a-b plane is called the
complex plane.

We now see how the complex wave function represents an oscillation. If ψ = exp[i(kx-ωt)], the complex function ψ(x,t) moves round and round the unit circle in the complex plane as x and t change, as illustrated in figure 9.1. This contrasts with the back and forth oscillation along the horizontal axis of the complex plane represented by cos(kx - ωt).

We will not present a formal proof of Euler’s equation — you will eventually see it in your calculus course. However, it may be helpful to note that the ϕ derivatives of exp(iϕ) and cos(ϕ) + i sin(ϕ) have the same behavior:

| (9.2) |

We indicate the complex conjugate of a complex number z by a superscripted
asterisk, i. e., z^{*}. It is obtained by replacing i by -i. Thus, (a + ib)^{*} = a-ib. The
absolute square of a complex number is the number times its complex
conjugate:

| (9.4) |

Notice that the absolute square of a complex exponential function is one:

| (9.5) |

In quantum mechanics the absolute square of the wave function at any point
expresses the relative probability of finding the associated particle at that point.
Thus, the probability of finding a particle represented by a plane wave is uniform
in space. Contrast this with the relative probability associated with a
sine wave: | sin(kx - ωt)|^{2} = sin ^{2}(kx - ωt). This varies from zero to one,
depending on the phase of the wave. The “waviness” in a complex exponential
plane wave resides in the phase rather than in the magnitude of the wave
function.

One more piece of mathematics is needed. The complex conjugate of Euler’s equation is

| (9.6) |

Taking the sum and the difference of this with the original Euler’s equation results in the expression of the sine and cosine in terms of complex exponentials:

| (9.7) |

We aren’t used to having complex numbers show up in physical theories and it is hard to imagine how we would measure such a number. However, everything observable comes from taking the absolute square of a wave function, so we deal only with real numbers in experiments.

The idea of symmetry plays a huge role in physics. We have already used symmetry arguments in the theory of relativity — applying the principle of relativity to obtain the dispersion relation for relativistic matter waves is just such an argument. In this section we begin to explore how symmetry can be used to increase our understanding of quantum mechanics.

For our first example we take the case of a free particle in quantum mechanics, i. e., a particle subject to no force. The wave function for a free particle of definite momentum Π and energy E is given by

| (9.8) |

For this wave function |ψ|^{2} = 1 everywhere, so the probability of finding the
particle anywhere in space and time is uniform. This contrasts with the
probability distribution which arises if we assume a free particle to have the wave
function ψ = cos[(Πx - Et)∕ℏ]. In this case |ψ|^{2} = cos ^{2}[(Πx - Et)∕ℏ], which
varies with position and time, and is inconsistent with a uniform probability
distribution.

Quantum mechanics is a probabilistic theory, in the sense that the predictions it makes tell us, for instance, the probability of finding a particle somewhere in space. If we know nothing about a particle’s previous history, and if there are no physical constraints that would make it more likely for a particle to be at one point along the x axis than any another, then the probability distribution must be P(x) = constant.

This is an example of a symmetry argument. Expressed more formally, it states that if the above conditions apply, then the probability distribution ought to be subject to the condition P(x + D) = P(x) for any constant value of D. The only possible P(x) in this case is P = constant. In the language of physics, if there is nothing that gives the particle a higher probability of being at one point rather than another, then the probability is independent of position and the system is invariant under displacement in the x direction.

The above argument doesn’t suffice for quantum mechanics, since as we have learned, the fundamental quantity describing a particle is not the probability distribution, but the wave function ψ(x). Thus, the wave function rather than the probability distribution ought to be the quantity which is invariant under displacement, i. e., ψ(x + D) = ψ(x).

This condition turns out to be too restrictive, because it implies that ψ(x) = constant, whereas we know that a one-dimensional plane wave, which describes a particle with a uniform probability of being found anywhere along the x axis, has the form ψ(x) = exp(ikx). (For simplicity we temporarily ignore the time dependence.) If we make the substitution x → x + D in a plane wave, we get exp[ik(x + D)] = exp(ikx) exp(ikD). The wave function is thus technically not invariant under displacement, in that the displaced wave function is multiplied by the factor exp(ikD). However, the probability distribution of the displaced wave function still equals one everywhere, so there is no change in what we observe. Thus, in determining invariance under displacement, we are allowed to ignore changes in the wave function which consist only of multiplying it by a complex constant with an absolute value of one. Such a multiplicative constant is called a phase factor.

It is easy to convince oneself by trial and error or by more sophisticated means that the only form of wave function ψ(x) which satisfies the condition ψ(x + D) = ψ(x) × (phase factor) is ψ(x) = A exp(ikx) where A is a (possibly complex) constant. This is just in the form of a complex exponential plane wave with wavenumber k. Thus, not only is the complex exponential wave function invariant under displacements in the manner defined above, it is the only wave function which is invariant to displacements. Furthermore, the phase factor which appears for a displacement D of such a plane wave takes the form exp(iC) = exp(ikD), where k is the wavenumber of the plane wave.

As an experiment, let us see if a wave packet is invariant under displacement. Let’s define a wave packet consisting of two plane waves:

| (9.9) |

Making the substitution x → x + D in this case results in

The impossibility of writing ψ(x + D) = ψ(x) × (phase factor) lends plausibility to the assertion that a single complex exponential is the only possible form of the wave function that is invariant under displacement. Notice that the wave packet does not have definite wavenumber, and hence,
momentum. In particular, the wave packet is a sum of complex exponentials with
wavenumbers k_{1} and k_{2}, which means that the associated particle can have either
momentum Π_{1} = ℏk_{1} or Π_{2} = ℏk_{2}. This makes sense from the point of view of
the uncertainty principle – for a single plane wave the uncertainty in
position is complete and the uncertainty in momentum is zero. For a wave
packet the uncertainty in position is reduced and the uncertainty in the
momentum is non-zero. As we have seen, this idea can be carried further: A
definite value of momentum must be associated with a completely indefinite
probability distribution in position, i. e., with P = constant. This corresponds
to a wave function which has the form of a complex exponential plane
wave. However, such a plane wave is invariant under displacement D,
except for the multiplicative phase factor exp(ikD), which has no physical
consequences since it disappears when the probability distribution is obtained.
Thus, we see that invariance under displacement of the wave function and
a definite value of the momentum are linked, in that each implies the
other:

| (9.11) |

The idea of potential energy was introduced in the previous chapter. In particular, we found that if the total energy is constant, then the momentum cannot be constant in the presence of spatially varying potential energy. This means that the wavenumber, and hence the wavelength of the oscillations in the wave function also vary with position. The spatial inhomogeneity of the potential energy gives rise to spatial inhomogeneity in the wave function, and hence an indefinite momentum.

The above argument can be extended to other variables besides momentum. In particular since the time dependence of a complex exponential plane wave is exp(-iωt) = exp(-iEt∕ℏ), where E is the total energy, we have by analogy with the above argument that

| (9.12) |

Thus, invariance of the wave function under a displacement in time implies a definite value of the energy of the associated particle.

In the previous chapter we assumed that the frequency (and hence the energy) was definite and constant for a particle passing through a region of variable potential energy. We now see that this assumption is justified only if the potential energy doesn’t change with time. This is because a time-varying potential energy eliminates the possibility of invariance under time shift.

We have seen a few examples of invariance in quantum mechanics. It is now time to define this concept more precisely. A quantum mechanical wave function is said to be invariant under some transformation if the transformed wave function is observationally indistinguishable from the original.

In the above examples, the transformation is accomplished by replacing x by x + D in the case of displacement in space and similarly by replacing t by t + T for displacement in time. However, the idea of a transformation is much more general; other examples will be discussed as they arise.

The idea of “observationally indistinguishable” can be tricky. For example, if some transformation results in a new wave function which is the old wave function times a constant phase factor, then the new wave function is observationally indistinguishable from the old one. This is because physical measurements capture phase differences between different parts of wave functions (think of how interferometers work), but not absolute phases. The constant phase factor disappears in this difference calculation. However, if the multiplicative phase factor created by some transformation is a function of position, then the phase difference between different parts of a wave function changes as a result of the transformation. The wave function is not invariant under this transformation.

We already know that definite values of certain pairs of variables cannot be obtained simultaneously in quantum mechanics. For instance, the indefiniteness of position and momentum are related by the uncertainty principle — a definite value of position implies an indefinite value of the momentum and vice versa. If definite values of two variables can be simultaneously obtained, then we call these variables compatible. If not, the variables are incompatible.

If the wave function of a particle is invariant under the displacements associated with both variables, then the variables are compatible. For instance, the complex exponential plane wave associated with a free particle is invariant under displacements in both space and time. Since momentum is associated with space displacements and energy with time displacements, the momentum and energy are compatible variables for a free particle.

Variables which are compatible with the energy have a special status. The wave function which corresponds to a definite value of such a variable is invariant to displacements in time. In other words, the wave function doesn’t change under this displacement except for a trivial phase factor. Thus, if the wave function is also invariant to some other transformation at a particular time, it is invariant to that transformation for all time. The variable associated with that transformation therefore retains its definite value for all time — i. e., it is conserved.

For example, the plane wave implies a definite value of energy, and is thus invariant under time displacements. At time t = 0, it is also invariant under x displacements, which corresponds to the fact that it represents a particle with a known value of momentum. However, since momentum and energy are compatible for a free particle, the wave function will represent the same value of momentum at all other times. In other words, if the momentum is definite at t = 0, it will be definite at all later times, and furthermore will have the same value. This is how the conservation of momentum (and by extension, the conservation of any other variable compatible with energy) is expressed in quantum mechanics.

In modern quantum physics, the discovery of new symmetries leads to new dynamical variables. In the problems we show how that comes about for the symmetries of parity (x →-x), time reversal t →-t), and charge conjugation (the interchange of particles with antiparticles). One of the key examples of this was the development of the quark theory of matter, which came from the observation that the interchange of certain groups of elementary particles left the universe approximately unchanged, meaning that the universe was (approximately) symmetric under these interchanges.

Confinement of a wave to a limited spatial region results in rather peculiar behavior — the wave can only fit comfortably into the confined region if the wave frequency, and hence the associated particle energy, takes on a limited set of possible values. This is the origin of the famous quantization of energy, from which the “quantum” in quantum mechanics comes. We will explore two types of confinement, position confinement due to a potential energy well, and rotational confinement due to the fact that rotation of an object through 2π radians returns the object to its original orientation.

We now imagine how a particle confined to a region 0 ≤ x ≤ a on the x axis
must behave. As with the displacement of a guitar string, the wave function must
be zero at x = 0 and a, i. e., at the ends of the guitar string. A single complex
exponential plane wave cannot satisfy this condition, since | exp[i(kx - ωt)]|^{2} = 1
everywhere. However, a superposition (with a minus sign) of leftward and
rightward traveling waves creates a standing wave, in which the the wave
function separates into a function of space alone times a function of time
alone.

| (9.13) |

Notice that the time dependence is still a complex exponential, which means that
|ψ|^{2} is independent of time. This insures that the probability of finding the
particle somewhere in the box remains constant with time. It also means that the
wave packet corresponds to a definite energy E = ℏω.

Because we took a difference rather than a sum of plane waves, the condition ψ = 0 is already satisfied at x = 0. To satisfy it at x = a, we must have ka = nπ, where n = 1, 2, 3,…. Thus, the absolute value of the wavenumber must take on the discrete values

| (9.14) |

(The wavenumbers of the two plane waves equal plus or minus this absolute value
respectively.) This implies that the absolute value of the particle momentum is
Π_{n} = ℏk_{n} = nπℏ∕a, which in turn means that the energy of the particle must
be

| (9.15) |

where m is the particle mass. In the non-relativistic limit this becomes

| (9.16) |

where we have dropped the rest energy mc^{2} since it is a constant offset.
In the ultra-relativistic case where we can ignore the particle mass, we
find

| (9.17) |

The shapes of the wave functions for the first three values of n for the particle in the box are illustrated in figure 9.2.

In both limits the energy takes on only a certain set of possible values.
This is called energy quantization and the integer n is called the energy
quantum number. In the non-relativistic limit the energy is proportional
to n^{2}, while in the ultra-relativistic case the energy is proportional to
n.

We can graphically represent the allowed energy levels for the particle in a box by an energy level diagram. Such a diagram is shown in figure 9.3 for the non-relativistic case.

One aspect of this problem deserves a closer look. Equation (9.13) shows that
the wave function for this problem is a superposition of two plane waves
corresponding to momenta Π_{1} = +ℏk and Π_{2} = -ℏk and is therefore a kind of
wave packet. Thus, the wave function is not invariant under displacement and
does not correspond to a definite value of the momentum — the momentum’s
absolute value is definite, but its sign is not. Following Feynman’s prescription,
equation (9.13) tells us that the amplitude for the particle in the box to have
momentum +ℏk is exp[i(kx - ωt)], while the amplitude for it to have
momentum -ℏk is - exp[i(-kx - ωt)]. The absolute square of the sum of these
amplitudes gives us the relative probability of finding the particle at position
x:

| (9.18) |

Which of the two possible values of the momentum the particle takes on is unknowable, just as it is impossible in principle to know which slit a particle passes through in two slit interference. If an experiment is done to measure the momentum, then the wave function is irreversibly changed, just as the interference pattern in the two slit problem is destroyed if the slit through which the particle passes is unambiguously determined.

Unlike the situation in classical mechanics, quantum mechanics allows the kinetic
energy K to be negative. This makes the momentum Π (equal to (2mK)^{1∕2} in the
nonrelativistic case) imaginary, which in turn gives rise to an imaginary
wavenumber.

Let us investigate the nature of a wave with an imaginary wavenumber. Let us assume that k = iκ in a complex exponential plane wave, where κ is real:

| (9.19) |

The wave function doesn’t oscillate in space when K = E - U < 0, but grows or decays exponentially with x, depending on the sign of κ.

For a particle moving to the right, with positive k in the allowed region, κ turns out to be positive, and the solution decays to the right. Thus, a particle impingent on a potential energy barrier from the left (i. e., while moving to the right) will have its wave amplitude decay in the classically forbidden region, as illustrated in figure 9.4. If this decay is very rapid, then the result is almost indistinguishable from the classical result — the particle cannot penetrate into the forbidden region to any great extent. However, if the decay is slow, then there is a reasonable chance of finding the particle in the forbidden region. If the forbidden region is finite in extent, then the wave amplitude will be small, but non-zero at its right boundary, implying that the particle has a finite chance of completely passing through the classical forbidden region. This process is called barrier penetration.

The probability for a particle to penetrate a barrier is the absolute square of the amplitude after the barrier divided by the square of the amplitude before the barrier. Thus, in the case of the wave function illustrated in equation (9.19), the probability of penetration is

| (9.20) |

where d is the thickness of the barrier.

The rate of exponential decay with x in the forbidden region is related to how negative K is in this region. Since

| (9.21) |

we find that

| (9.22) |

where the potential energy barrier is B ≡-K = U - E. The smaller B is, the smaller is κ, resulting in less rapid decay of the wave function with x. This corresponds to stronger barrier penetration. (Note that the way B is defined, it is positive in forbidden regions.)

If the energy barrier is very high, then the exponential decay of the wave function is very rapid. In this case the wave function goes nearly to zero at the boundary between the allowed and forbidden regions. This is why we specify the wave function to be zero at the walls for the particle in a box. These walls act in effect as infinitely high potential barriers.

Barrier penetration is important in a number of natural phenomena. Certain types of radioactive decay and the fissioning of heavy nuclei are governed by this process.

Another type of bound state motion occurs when a particle is constrained to move in a circle. (Imagine a bead sliding on a circular loop of wire, as illustrated in figure 9.5.) We can define x in this case as the path length around the wire and relate it to the angle θ: x = Rθ. For a plane wave we have

| (9.23) |

This plane wave differs from the normal plane wave for motion along a Cartesian axis in that we must have ψ(θ) = ψ(θ + 2π). This can only happen if the circumference of the loop, 2πR, is an integral number of wavelengths, i. e., if 2πR∕λ = m where m is an integer. However, since 2π∕λ = k, this condition becomes kR = m.

Since Π = ℏk, the above condition can be written Π_{m}R = mℏ. The
quantity

| (9.24) |

is called the angular momentum, leading to our final result,

| (9.25) |

We see that the angular momentum can only take on values which are integer multiples of ℏ. This represents the quantization of angular momentum, and m in this case is called the angular momentum quantum number. Note that this quantum number differs from the energy quantum number for the particle in the box in that zero and negative values are allowed.

The energy of our bead on a loop of wire can be expressed in terms of the angular momentum:

| (9.26) |

This means that angular momentum and energy are compatible variables in this case, which further means that angular momentum is a conserved variable. Just as definite values of linear momentum are related to invariance under translations, definite values of angular momentum are related to invariance under rotations. Thus, we have

| (9.27) |

for angular momentum.

We need to briefly address the issue of angular momentum in three dimensions. Angular momentum is actually a vector oriented perpendicular to the wire loop in the example we are discussing. The direction of the vector is defined using a variation on the right-hand rule: Curl your fingers in the direction of motion of the bead around the loop (using your right hand!). The orientation of the angular momentum vector is defined by the direction in which your thumb points. This tells you, for instance, that the angular momentum in figure 9.5 points out of the page.

In quantum mechanics it is only possible to measure simultaneously the square of the length of the angular momentum vector and one component of this vector. Two different components of angular momentum cannot be simultaneously measured because of the uncertainty principle. However, the length of the angular momentum vector may be measured simultaneously with one component. Thus, in quantum mechanics, the angular momentum is completely specified if the length and one component of the angular momentum vector are known.

Figure 9.6 illustrates the angular momentum vector associated with a bead
moving on a wire loop which is tilted from the horizontal. One component (taken
to be the z component) is shown as well. For reasons we cannot explore here, the
square of the length of the angular momentum vector L^{2} is quantized with the
following values:

| (9.28) |

One component (say, the z component) of angular momentum is quantized just
like angular momentum in the two-dimensional case, except that l acts as an
upper bound on the possible values of |m|. In other words, if the square of the
length of the angular momentum vector is ℏ^{2}l(l + 1), then the z component can
take on the values

| (9.29) |

The quantity l is called the angular momentum quantum number, while m is called the orientation or magnetic quantum number, the latter for historical reasons.

The type of angular momentum discussed above is associated with the movement of particles in orbits. However, it turns out that even stationary particles can possess angular momentum. This is called spin angular momentum. The spin quantum number s plays a role analogous to l for spin angular momentum, i. e., the square of the spin angular momentum vector of a particle is

| (9.30) |

The spin orientation quantum number m_{s} is similarly related to s:

| (9.31) |

The spin angular momentum for an elementary particle is absolutely conserved, i. e., it can never change. Thus, the value of s is an intrinsic property of a particle. The major difference between spin and orbital angular momentum is that the spin quantum number can take on more values, i. e., s = 0, 1∕2, 1, 3∕2, 2, 5∕2,….

Particles with integer spin values s = 0, 1, 2,… are called bosons after the Indian physicist Satyendra Nath Bose. Particles with half-integer spin values s = 1∕2, 3∕2, 5∕2,… are called fermions after the Italian physicist Enrico Fermi. As we shall see later in the course, bosons and fermions play very different roles in the universe.

- Suppose that a particle is represented by the wave function
ψ = sin(kx - ωt) + sin(-kx - ωt).
- Use trigonometry to simplify this wave function.
- Compute the x and t dependence of the probability of finding the particle by squaring the wave function.
- Explain what this result says about the time dependence of the probability of finding the particle. Does this make sense?

- Repeat the above problem for a particle represented by the wave function ψ = exp[i(kx - ωt)] + exp[i(-kx - ωt)].
- Determine if the wavefunction ψ(x) = exp(iCx
^{2}) is invariant under displacement in the sense that the displaced wave function differs from the original wave function by just a phase factor. - Just as invariance under the substitution x → x + D is associated with
momentum, invariance under the substitution x →-x is associated with a
quantum mechanical variable called parity, denoted P. However, unlike
momentum, which can take on any numerical value, parity can take on only
two possible values, ±1. The parity of a wave function ψ(x) is +1 if
ψ(-x) = ψ(x), while the parity is -1 if ψ(-x) = -ψ(x). If ψ(x)
satisfies neither of these conditions, then it has no definite value of
parity.
- What is the parity of ψ = sin(kx)? Of ψ = cos(kx)? The quantity k is a constant.
- Is ψ(x) = cos(kx) invariant under the substitution x = x + D for all possible values of D? Does this wave function have a definite value of the momentum?
- Show that a wave function with a definite value of the momentum does not have a definite value of parity. Are momentum and parity compatible variables?

- Realizing that cos(kx - ωt) can be written in terms of complex exponential functions, give a physical interpretation of the meaning of the above cosine wave function. In particular, what are the possible values of the associated particle’s momentum and energy?
- The time reversal operation T makes the substitution t →-t. Similar to parity, time reversal can only take on values ±1. Is symmetry of a wave function under time reversal, i. e., ψ(-t) = ψ(t), consistent with a definite value of the energy? Hint: Any wave function corresponding to a definite value of energy E must have the form ψ = A exp(-iEt∕ℏ) where A is not a function of time t. (Why?)
- The operation C takes the complex conjugate of the wave function, i. e., it makes the substitution i →-i. In modern quantum mechanics this corresponds to interchanging particles and antiparticles, and is called charge conjugation. What does the combined operation CPT do to a complex plane wave, i. e., one with definite wave vector and frequency?
- Make an energy level diagram for the case of a massless particle in a box.
- Compare |Π| for the ground state of a non-relativistic particle in a box of size a with ΔΠ obtained from the uncertainty principle in this situation. Hint: What should you take for Δx?
- Imagine that a billiard table has an infinitely high rim around it. For this
problem assume that ℏ = 1 kg m
^{2}s^{-1}.- If the table is 1.5 m long and if the mass of a billiard ball is M = 0.5 kg, what is the billiard ball’s lowest or ground state energy? Hint: Even though the billiard table is two dimensional, treat this as a one-dimensional problem. Also, treat the problem nonrelativistically and ignore the contribution of the rest energy to the total energy.
- The energy required to lift the ball over a rim of height H against
gravity is U = MgH where g = 9.8 m s
^{-2}. What rim height makes the gravitational potential energy equal to the ground state energy of the billiard ball calculated above? - If the rim is actually twice as high as calculated above but is only 0.1 m thick, determine the probability of the ball penetrating the rim.

- The real part of the wave function of a particle with positive energy E
passing through a region of negative potential energy is shown in figure
9.7.
- If the total energy is definitely E, what is the dependence of this wave function on time?
- Is the wave function invariant under displacement in space in this case? Why or why not?
- Does this wave function correspond to a definite value of momentum? Why or why not?
- Is the momentum compatible with the energy in this case? Why or why not?

- Assuming again that ℏ = 1 kg m
^{2}s^{-1}, what are the possible speeds of a toy train of mass 3 kg running around a circular track of radius 0.8 m? - If a particle of zero mass sliding around a circular loop of radius R can take
on angular momenta L
_{m}= mℏ where m is an integer, what are the possible kinetic energies of the particle? Hint: Remember that L = ΠR.

Dynamics of Multiple Particles

So far we have considered only the dynamics of a single particle subject to an externally imposed potential energy. The particle has no way of influencing this external agent. In the real world particles interact with each other. In this chapter we learn how this happens.

We first rewrite Newton’s second law in terms of momentum. This is useful in the subsequent consideration of Newton’s third law, which leads to the principle of the conservation of momentum. Collisions between particles and the behavior of rockets and conveyor belts are then studied as applications of the conservation laws to more than one particle.

Up to this point we have stated Newton’s second law in its conventional form, F = ma. However, in the non-relativistic case ma = mdu∕dt = d(mu)∕dt, so we can also write Newton’s second law as

| (10.1) |

where p = mu is the non-relativistic kinetic momentum. This form of Newton’s
second law is actually closer to Newton’s original statement of the law. It also has
the advantage that it is correct even in the relativistic case when the relativistic
definition of kinetic momentum, p = mu∕(1 - u^{2}∕c^{2})^{1∕2} (as defined earlier), is
substituted.

Newton’s third law states that if particle A exerts a force F on particle B, then particle B exerts a force -F on particle A. Newton’s third law makes it possible to apply Newton’s second law to systems of particles without considering the detailed interactions between particles within the system. For instance, if we (arbitrarily) define the system in figure 10.1 to be the particles A and B inside the dashed lines, then we can divide the forces acting on these particles into internal and external parts,

| (10.2) |

| (10.3) |

Adding these equations together results in the net force F_{net} being equal
to

| (10.4) |

However, the internal interactions in this case are A acting on B and B acting on
A. These forces are equal in magnitude but opposite in direction, so they cancel
out, leaving us with the net force equal to the sum of the external parts,
F_{net} = F_{A-external} + F_{B-external}. The external forces in figure 10.1 are the force of
C on A and the force of C on B. Defining the total kinetic momentum of the
system as the sum of the A and B momenta, p_{tot} = p_{A} + p_{B}, the above equation
becomes

| (10.5) |

which looks just like Newton’s second law for a single particle, except that it now applies to the system of particles (A and B in the present case) as a whole. This argument easily generalizes to any number of particles inside and outside the system. Thus, for instance, even though a soccer ball consists of billions of atoms, we are sure that the forces between atoms within the soccer ball cancel out, and the trajectory of the ball as a whole is determined solely by external forces such as gravity, wind drag, friction with the ground, and the kicks of soccer players.

Remember that for two forces to be a third law pair, they have to be acting on different particles. Furthermore, if one member of the pair is the force of particle A acting on particle B, then the other must be the force of particle B acting on particle A. A counterexample would be gravity and the upward normal force acting on a mass sitting on a table; these forces are equal and opposite in the stationary case, but act on the same object, and therefore do not constitute a third law pair. However, the upward normal force of the table on the mass and the downward normal force of the mass on the table would be a third law pair.

If all external forces on a system are zero, then equation (10.5) reduces to

| (10.6) |

A system of particles with no external forces acting on it is called isolated. Newton’s third law thus tells us that the kinetic momentum of an isolated system doesn’t change with time. This law is called the conservation of momentum.

Let us now consider the situation in which two particles collide with each other. There can be several outcomes to this collision, of which we will study two:

- The two particles collide elastically, in essence bouncing off of each other.
- The two particles stick together, resulting in the production of a single particle, or a single particle breaks apart into two particles. These are inelastic processes.

In both of the above cases energy and momentum are conserved. We assume that the forces acting between the particles are short range, so that except in the instant of collision, we need not worry about potential energy or potential momentum — all energy is in the form of rest plus kinetic energy except in this short interval, and all momenta are kinetic momenta.

Because of the principle of relativity, we are free to consider collisions in any convenient reference frame. We can then transform the results to any reference frame we please. Generally speaking, the most convenient reference frame to consider is the one in which the total momentum of the two particles is zero. For the sake of simplicity we only consider collisions in one dimension.

Suppose a particle with mass m_{1} and initial velocity u_{1} in the center of
momentum frame, i. e., the reference frame in which the total momentum is zero,
collides elastically with another particle of mass m_{2} with initial velocity u_{2}. The
momenta of the two particles are

| (10.7) |

In the center of momentum frame we must have

| (10.8) |

Figure 10.2 shows what happens when these two particles collide. The first
particle acquires momentum p_{1}′ while the second acquires momentum p_{2}′. The
conservation of momentum tells us that the total momentum after the collision is
the same as before the collision, namely zero, so

| (10.9) |

In the center of momentum frame we know that |p_{1}| = |p_{2}| and we know that
the two momentum vectors point in opposite directions. Similarly, |p_{1}′| = |p_{2}′|.
However, we as yet don’t know how p_{1}′ is related to p_{1}. Conservation of
energy,

| (10.10) |

gives us this information. Notice that if p_{1}′ = -p_{1}, then
E_{1}′^{2} = p_{
1}′^{2}c^{2} + m_{
1}^{2}c^{4} = p_{
1}^{2}c^{2} + m_{
1}^{2}c^{4} = E_{
1}^{2}. Assuming positive energies, we
therefore have E_{1}′ = E_{1}. If p_{2}′ = -p_{2}, then we can similiarly infer that E_{2}′ = E_{2}.
If these conditions are satisfied, then so is equation (10.10). Therefore, a complete
solution to the problem is

| (10.11) |

and

| (10.12) |

In other words, the particles just exchange momenta.

The left panel of figure 10.2 shows what happens in a collision when the
masses of the two colliding particles are equal. If m_{1} = m_{2}, then the incoming and
outgoing velocities of the two particles are the same, as indicated by the inverse
slopes of the world lines. On the other hand, if m_{1} > m_{2}, then the velocity of
particle 2 is greater than the velocity of particle 1, as is illustrated in the right
panel of figure 10.2.

Suppose we wish to view the results of an elastic collision in a reference frame
in which particle 2 is initially stationary. All we have to do is to transform the
velocities into a reference frame moving with the initial velocity of particle 2, as
illustrated in figure 10.3. We do this by relativistically adding U = -u_{2} to each
velocity. (Note that the velocity U of the moving frame is positive since u_{2}
is negative.) Using the relativistic velocity translation formula, we find
that

| (10.13) |

where u_{1}, u_{1}′, u_{2}, and u_{2}′ indicate velocities in the original, center of momentum
reference frame and v_{1}, v_{1}′, etc., indicate velocities in the transformed
frame.

In the special case where the masses of the two particles are equal to each
other, we have v_{1} = 2U∕(1 + U^{2}∕c^{2}), v_{
1}′ = 0, and v_{2}′ = 2U∕(1 + U^{2}∕c^{2}) = v_{
1}.
Thus, when the masses are equal, the particles simply exchange velocities.

If the velocities are nonrelativistic, then the simpler Galilean transformation law v = u + U can be used in place of the relativistic equations invoked above.

An inelastic collision is one in which the particles coming out of the collision are not the same as the particles going into it. Inelastic collisions conserve both total momentum and energy just as elastic collisions do. However, unlike elastic collisions, inelastic collisions generally do not conserve the total kinetic energy of the particles, as some rest energy is generally created or destroyed.

Figure 10.4 shows the fundamental building blocks of inelastic collisions. We can consider even the most complex inelastic collisions to be made up of composites of only two processes, the creation of one particle from two, and the disintegration of one particle into two.

Let us consider each of these in the center of momentum frame. In both cases
the single particle must be stationary in this frame since it carries the total
momentum of the system, which has to be zero. By conservation of momentum, if
particle 1 in the left panel of figure 10.4 has momentum p, then the momentum of
particle 2 is -p. If the two particles have masses m_{1} and m_{2}, then their energies
are E_{1} = (p^{2}c^{2} + m_{
1}^{2}c^{4})^{1∕2} and E_{
2} = (p^{2}c^{2} + m_{
2}^{2}c^{4})^{1∕2}. The energy of particle
3 is therefore E_{3} = E_{1} + E_{2}, and since it is at rest, all of its energy is
in the form of “mc^{2}” or rest energy, and so the mass of this particle is

The last line in the above equation shows that m_{3} > m_{1} + m_{2} because it is in
the form m_{1}A + m_{2}B where both A and B are greater than one. Thus, rest energy
is created in the amount ΔE_{rest} = (m_{3} - m_{1} - m_{2})c^{2}.

Actually, it is easy to calculate the mass of particle 3 in the above case from
any reference frame as long as the momenta and energies of particles 1 and 2 are
known in this frame. By conservation of energy and momentum, E_{3} = E_{1} + E_{2}
and p_{3} = p_{1} + p_{2}. Furthermore, E_{3}^{2} = p_{
3}^{2}c^{2} + m_{
3}^{2}c^{4}, so we can solve for
m_{3}:

| (10.15) |

The right panel of figure 10.4 shows the process of particle decay. This is just the inverse of the particle creation process, and all of the analysis we have done for creation is valid for particle decay except that rest energy is converted to kinetic energy rather than vice versa.

Normally when we define a system to which Newton’s second law is to be applied, the system is closed in the sense that mass cannot enter or exit the system. However, sometimes it is convenient to work with open systems for which this is not true. The classic example is the rocket, where exhaust gases leave the system, thus decreasing the mass of the rocket with time.

Open systems can be analyzed if momentum is considered to be a quantity which is accounted for much as money is accounted for in a bank account. The bank account can change in three ways: money can be deposited in the account, it can be withdrawn from the account, and the amount can grow or shrink as a consequence of interest payments or fees. Similarly, the amount of momentum in a system can change as the result of mass entering the system, mass leaving the system, and forces acting on the system. The time rate of change of momentum in a system is therefore

| (10.16) |

where F is the net force on the system, (dp∕dt)_{in} is the momentum per unit time
added by mass entering the system, and (dp∕dt)_{out} is the amount lost
per unit time by mass exiting the system. In the non-relativistic case,
(dp∕dt)_{in} = u_{in}(dm∕dt)_{in} and (dp∕dt)_{out} = u_{out}(dm∕dt)_{out}, where (dm∕dt)_{in} is the
mass entering the system per unit time with velocity u_{in} and (dm∕dt)_{out} is the
mass per unit time exiting the system with velocity u_{out}.

For non-relativistic velocities, the momentum of the system can be written as p = mu so that

| (10.17) |

To complete the analysis, we need an accounting of the mass entering and leaving the system:

| (10.18) |

Let us see how to apply this to a rocket for which all velocities are
non-relativistic. As figure 10.5 indicates, a rocket spews out a stream of exhaust
gas. The system is defined by the dashed box and includes the rocket and the part
of the exhaust gas inside the box. The reaction to the momentum carried off in
this stream of gas is what causes the rocket to accelerate. We note that
(dm∕dt)_{in} = 0 since no mass is entering the system, and (dm∕dt)_{out} = R, the rate
at which mass is ejected by the rocket in the form of exhaust gas. The
rocket is assumed to be moving to the right at speed V and the gas is
ejected at a speed u_{x} relative to the rocket, which means that its actual
velocity after ejection is V - u_{x}. We call u_{x} the exhaust velocity. Notice
that V - u_{x} may be either positive or negative, depending on how big V
is.

Equating the mass of the rocket to the system mass, we find that R = -dm∕dt.
The momentum balance equation (10.16) becomes dp∕dt = -(V - u_{x})R.
The force on the rocket is actually zero, so the force term does not enter
the momentum balance equation. This is non-intuitive, because we are
used to acceleration being the result of a force. However, nothing,
including the ejected gas, is actually pushing on the system, so we
must indeed conclude that there is no force — all of the change in the
system’s momentum arises from the ejection of gas with the opposite
momentum.^{1}

Finally, we see that dp∕dt = (dm∕dt)V + m(dV∕dt) = -RV + m(dV∕dt).
Equating this to the results of the momentum balance calculation gives us
-RV + m(dV∕dt) = -(V - u_{x})R. Solving for the acceleration dV∕dt results
in

| (10.19) |

Thus, the acceleration of the rocket depends on the exhaust velocity of the ejected gas, the rate at which the gas is being ejected, and the mass of the rocket.

Figure 10.6 illustrates another type of open system problem. A hopper
dumps sand on a conveyor belt at a rate of R kilograms per second. The
conveyor belt is moving to the right at (non-relativistic) speed V and
the sand is dumped off at the end. What force F is needed to keep the
conveyor belt moving at a constant speed, assuming that the conveyor belt
mechanism itself is frictionless? In this case (dm∕dt)_{in} = (dm∕dt)_{out} = R.
Furthermore, since the system outlined by the dashed line is in a steady state,
dp∕dt = 0.

The key to understanding this problem is that the sand enters the system with zero horizontal velocity, but exits the system with the horizontal velocity of the conveyor belt, V . The momentum balance equation is thus

| (10.20) |

and the force is

| (10.21) |

This force serves to accelerate the sand up to the velocity of the conveyor belt.

- Imagine a block of mass M resting on a plate under the influence of gravity,
as shown in figure 10.7.
- Determine the force of the plate on the block, N
_{b}, and the force of the block on the plate, N_{p}. - State which of the three forces, Mg, N
_{b}, and N_{p}, form a Newton’s third law pair.

- Determine the force of the plate on the block, N
- Repeat the previous problem assuming that the block and the plate are in an elevator accelerating upward with acceleration a.
- Straighten out the misunderstanding of Newton’s third law implicit in the question “If the force of the horse on the cart equals the force of the cart on the horse, why does anything ever go anywhere”? Examine in particular the conditions under which the horse-cart system accelerates.
- A pusher boat (mass M) on the Mississippi is pushing two barges (each
mass m) at a steady speed as shown in figure 10.8. Each barge is subject to
a drag force by the water of F
_{B}. Consider only horizontal force components in the following.- What is the total horizontal force of the water on the barge-boat system? Explain.
- What is the direction and magnitude of the force of the pusher boat on barge 1? Explain.

- A train with an engine of mass M and 2 freight cars, each of mass
m, is accelerating to the right with acceleration a on a horizontal
track as shown in figure 10.9. Assume that the two freight cars roll
with negligible friction. Consider only horizontal force components
below.
- Find the direction and magnitude of the force of the rails on the engine and specify the system to which Newton’s second law is applied.
- Find the direction and magnitude of the force of the engine on the first car and specify the system to which Newton’s second law is applied.
- Find the direction and magnitude of the force of the first car on the second car and specify the system to which Newton’s second law is applied.
- Find the direction and magnitude of the force of the second car on the first car and specify the law used to obtain this force.

- A car and trailer are descending a hill as shown in figure 10.10. Assume that
the trailer rolls without friction and that air friction can be ignored.
Consider only forces parallel to the road surface.
- Compute the force of the road on the car if the car-trailer system shown in figure 10.10 is moving down the hill at constant speed.
- Compute the force of the trailer on the car in the above conditions.
- If the driver takes his foot off the brake and lets the car coast frictionlessly, recompute the force of the trailer on the car.

- Consider a one-dimensional elastic collision between particles of masses m
_{1}and m_{2}. If particle 2 is initially stationary, what range of values must m_{1}∕m_{2}have for the initial particle to rebound backwards along its initial track after the collision? (Do this problem non-relativistically.) - A stationary pion (mass M) decays into a muon (mass m < M) and a
neutrino (massless).
- What is the (fully relativistic) momentum of the muon after the decay?
- What is the energy of the neutrino?

- In an elastic collision viewed in the center of momentum frame, the energy of each particle is conserved individually. Is this true for the same process viewed from a reference frame in which one of the particles is initially stationary?
- A space probe approaches a planet in the -x direction, curves around it
under the influence of the planet’s powerful gravity (a conservative force)
and recedes from the planet in the +x direction, as seen in figure 10.11. The
planet is moving in the +x direction at speed V , while the space
probe is initially moving in the -x direction at speed u
_{1}. What is its speed u_{2}in the +x direction after this close approach to the planet? Treat this problem nonrelativistically. Hint: First transform to the center of mass frame in which the planet is essentially stationary. Work out the interaction between the probe and the planet in this frame. Then transform back to the original reference frame. Assume that the mass of the probe is negligible compared to that of the planet.

- Two asteroids, each with mass 10
^{10}kg and initial speed 10^{5}m s^{-1}, collide head on. The whole mess congeals into one large mass. How much rest mass (rest energy divided by c^{2}) is created? - Two equal objects, both with mass m, collide and stick together. Before the
collision, one mass is stationary and the other is moving at speed v. In the
following, assume that velocities are fully relativistic.
- Compute the total momentum and energy (including rest energy) of the two masses before the collision.
- Compute the mass M of the combined system after the collision, taking the conversion of energy into mass into account.

- Explain qualitatively why a fireman needs to push forward on a firehose to keep it stationary. Hint: The water is flowing faster after it comes out of the nozzle of the hose than before.
- Solve equation (10.19) for V as a function of m, assuming that
V = 0 and m = m
_{0}at t = 0. Hint: Since R = -dm∕dt, we have R∕m = -d ln(m)∕dt. - Bottles are filled with soft drink at a bottling plant as shown in figure 10.12.
The bottles sit on a scale which is used to determine when to shut off the
flow of soft drink. If the desired mass of the bottle plus soft drink after
filling is M, what weight should the scale read when the bottle is full? The
rate at which mass is being added to the bottle is R and its velocity entering
the bottle is V .

- An interstellar space probe has frontal area A, initial mass M
_{0}, and initial velocity V_{0}, which is non-relativistic. The tenuous gas between the stars has mass density ρ. These gas molecules stick to the probe when they hit it. Find the probe’s acceleration. Hint: In a frame of reference in which the gas is stationary, does the momentum of the space probe change with time? Does its mass? - A light beam with power J hits a plate which is oriented normally to
the beam. Compute the force required to hold the plate in place
if
- the plate completely absorbs the light, and
- the plate completely reflects the light.

Hint: Photons are massless, so the momentum of a photon with energy E is E∕c. Thus, the momentum per unit time hitting the plate is J∕c.

- Find the acceleration of a rocket when the exhaust “gas” is actually a laser beam of power J. Assume that the rocket moves at non-relativistic velocities and that the decrease in mass due to the loss of energy in the laser beam is negligible.

Rotational Dynamics

We have already seen the quantum mechanical treatment of angular momentum and rotational dynamics. In this section we study these subjects in a classical, non-relativistic context. We first define the concepts of torque and angular momentum in order to understand the orbital motion of a single particle. Next we examine two particles in arbitrary motion and learn how kinetic energy and angular momentum are partitioned between orbital and internal components. Two particles fixed to the ends of a light rod constitute a dumbbell, which serves as a prototype for the rotation of rigid bodies. We then see how what we learned for two particles extends to an arbitrary number of particles. Finally, we explore the physics of structures in static equilibrium.

Before we begin, we need to extend our knowledge of vectors to the cross product.

There are two ways to multiply two vectors together, the dot product and the cross product. We have already studied the dot product of two vectors, which results in a scalar or single number.

The cross product of two vectors results in a third vector, and is written symbolically as follows:

| (11.1) |

As illustrated in figure 11.1, the cross product of two vectors is perpendicular to the plane defined by these vectors. However, this doesn’t tell us whether the resulting vector in figure 11.1 points upward out of the plane or downward. This ambiguity is resolved using the right-hand rule:

- Point the uncurled fingers of your right hand along the direction of the first vector A.
- Rotate your arm until you can curl your fingers in the direction of the second vector B.
- Your stretched out thumb now points in the direction of the cross product vector C.

The magnitude of the cross product is given by

| (11.2) |

where |A| and |B| are the magnitudes of A and B, and θ is the angle between these two vectors. Note that the magnitude of the cross product is zero when the vectors are parallel or anti-parallel, and maximum when they are perpendicular. This contrasts with the dot product, which is maximum for parallel vectors and zero for perpendicular vectors.

Notice that the cross product does not commute, i. e., the order of the vectors is important. In particular, it is easy to show using the right-hand rule that

| (11.3) |

An alternate way to compute the cross product is most useful when the two
vectors are expressed in terms of components, i. e., A = (A_{x},A_{y}.A_{z}) and
B = (B_{x},B_{y},B_{z}):

Torque is the action of a force F on a mass M which induces it to revolve about some point, called the origin. It is defined

| (11.5) |

where r is the position of the mass relative to the origin, as illustrated in figure 11.2.

Notice that the torque is zero in a number of circumstances. If the force points directly toward or away from the origin, the cross product is zero, resulting in zero torque, even though the force is non-zero. Likewise, if r = 0, the torque is zero. Thus, a force acting at the origin produces no torque. Both of these limits make sense intuitively, since neither induces the mass to revolve around the origin.

The angular momentum of a mass M relative to a point O is

| (11.6) |

where p is the ordinary kinetic momentum of the
mass.^{1}
The angular momentum is zero if the motion of the object is directly towards or
away from the origin, or if it is located at the origin.

If we take the cross product of the position vector and Newton’s second law, we obtain an equation that relates torque and angular momentum:

| (11.7) |

The second term on the right side of the above equation is zero because dr∕dt equals the velocity of the mass, which is parallel to its momentum and the cross product of two parallel vectors is zero. This equation can therefore be written

| (11.8) |

It is the rotational version of Newton’s second law.

For both torque and angular momentum the location of the origin is arbitrary, and is generally chosen for maximum convenience. However, it is necessary to choose the same origin for both the torque and the angular momentum.

For the case of a central force, i. e., one which acts along the line of centers between two objects (such as gravity), there often exists a particularly convenient choice of origin. Imagine a planet revolving around the sun, as illustrated in figure 11.3. If the origin is placed at the center of the sun (which is assumed not to move under the influence of the planet’s gravity), then the torque exerted on the planet by the sun’s gravity is zero, which means that the angular momentum of the planet about the center of the sun is constant in time. No other choice of origin would yield this convenient result.

We already know about two fundamental conservation laws — those of energy and linear momentum. We believe that angular momentum is similarly conserved in isolated systems. In other words, particles can exchange angular momentum between themselves, but the vector sum of the angular momentum of all the particles in a system isolated from outside influences must remain constant.

Conservation of angular momentum is not an automatic consequence of the
conservation of linear momentum, even though the governing equation (11.8) for
angular momentum is derived from Newton’s second law. As an example, figure
11.4 shows a hypothetical situation in which the force F_{21} of M_{1} on M_{2} is equal
in magnitude but opposite in sign to the force F_{12} of M_{2} on M_{1}, i. e., Newton’s
third law holds, and the sum of the momenta of the two masses is conserved.
However, because the forces are non-central, the angular momentum of
the masses is not conserved. This scenario is impossible if the forces are
central.

Suppose we wish to apply Newton’s second law to two particles considered
together as a single system. As we showed previously, only external forces act on
the total momentum, p_{total} = p_{1} + p_{2}, of the two particles:

| (11.9) |

Let’s write the total non-relativistic momentum of the two particles in a special way:

| (11.10) |

where M_{total} = M_{1} + M_{2}. The quantity V_{cm} is the velocity of the center of mass
and can be expressed as the time derivative of the position of the center of mass,
R_{cm},

| (11.11) |

where

| (11.12) |

We now see how the kinetic energy and the angular momentum of the two
particles may be split into two parts, one having to do with the motion
of the center of mass of the two particles, the other having to do with
the motion of the two particles relative their center of mass. Figure 11.5
shows graphically how the vectors r′_{1} = r_{1} - R_{cm} and r′_{2} = r_{2} - R_{cm}
are defined. These vectors represent the positions of the two particles
relative to the center of mass. Substitution into equation (11.12) shows
that M_{1}r′_{1} + M_{2}r′_{2} = 0. This leads to the conclusion that M_{1}d_{1} = M_{2}d_{2}
in figure 11.5. We also define the velocity of each mass relative to the
center of mass as v′_{1} = dr′_{1}∕dt and v′_{2} = dr′_{2}∕dt, and we therefore have
M_{1}v′_{1} + M_{2}v′_{2} = 0.

The total kinetic energy is just the sum of the kinetic energies of the two
particles, K = M_{1}v_{1}^{2}∕2 + M_{
2}v_{2}^{2}∕2, where v_{
1} and v_{2} are the magnitudes of the
corresponding velocity vectors. Substitution of v_{1} = V_{cm} + v′_{1} etc., into the
kinetic energy formula and rearranging yields

| (11.13) |

Terms like V_{cm} ⋅ v′_{1} cancel out because M_{1}v′_{1} + M_{2}v′_{2} = 0.

The first term on the right side of equation (11.13) in square brackets is the kinetic energy the two particles would have if all of the mass were concentrated at the center of mass. The second term is the kinetic energy computed relative to the motion of the center of mass. The first is called the translational kinetic energy of the system while the second is called the internal kinetic energy.

The angular momentum of the system is just the sum of the angular momenta
of the two particles: L_{total} = M_{1}r_{1} × v_{1} + M_{2}r_{2} × v_{2}. By reasoning similar to the
case of kinetic energy, we can rewrite this as

| (11.14) |

The first term in square brackets on the right is called the orbital angular momentum while the second term is called the spin angular momentum. The former is the angular momentum the system would have if all the mass were concentrated at the center of mass, while the latter is the angular momentum of motion about the center of mass.

Interestingly, the idea of center of mass and the corresponding split of kinetic
energy and angular momentum into orbital and spin parts has no useful
relativistic generalization. This is due to the factor of γ ≡ (1 - v^{2}∕c^{2})^{-1∕2} in the
relativistic definition of momentum, p = mvγ, which means that

| (11.15) |

So far we have put no restrictions on the movements of the two particles. An interesting special case occurs when the particles are connected by a lightweight, rigid rod, giving us a dumbbell. In order to further simplify things, we assume that the rod is connected rigidly to a fixed axle at the center of mass of the two particles, as shown in figure 11.6. The masses constituting the ends of the dumbbell are therefore free to revolve in circles about the axle, but they are prevented from executing any other motion. The key effect of this constraint is that both masses rotate about the axle with the same angular frequency ω.

If the particles are respectively distances d_{1} and d_{2} from the axle, then their
speeds are v_{1} = d_{1}ω and v_{2} = d_{2}ω. Thus the kinetic energy of the rotating
dumbbell is

| (11.16) |

where I = M_{1}d_{1}^{2} + M_{
2}d_{2}^{2} is called the moment of inertia. Similarly, the
magnitude of the spin angular momentum, which is a vector parallel to the axle,
is

| (11.17) |

Finally, Newton’s second law for rotation becomes

| (11.18) |

where τ is the component of torque along the rotation axis.

Note that the rightmost expression in equation (11.18) assumes that I is
constant, which only is true if d_{1} and d_{2} are constant – i. e., the dumbbell must
truly be rigid.

The generalization from two particles to many particles is quite easy in principle. If a subscripted i indicates the value of a quantity for the ith particle, then the center of mass is given by

| (11.19) |

where

| (11.20) |

Furthermore, if we define r′_{i} = r_{i} - R_{cm}, etc., then the kinetic energy is
just

| (11.21) |

and the angular momentum is

| (11.22) |

In other words, both the kinetic energy and the angular momentum can be separated into two parts: one part is related to the overall motion of the system and the other is due to motions of system components relative to the center of mass, just as for the case of the dumbbell.

For a rigid body rotating about a fixed axle, the moment of inertia is

| (11.23) |

where d_{i} is the perpendicular distance of the ith particle from the axle. Equations
(11.16)-(11.18) are valid for a rigid body consisting of many particles.
Furthermore, the moment of inertia is constant in this case, so it can be taken out
of the time derivative:

| (11.24) |

The quantity α = dω∕dt is called the angular acceleration.

The sum in the equation for the moment of inertia can be converted to an integral for a continuous distribution of mass. We shall not pursue this here, but simply quote the results for a number of solid objects of uniform density:

- For rotation of a sphere of mass M and radius R about an axis piercing
its center: I = 2MR
^{2}∕5. - For rotation of a cylinder of mass M and radius R about its axis of
symmetry: I = MR
^{2}∕2. - For rotation of a thin rod of mass M and length L about an axis
perpendicular to the rod passing through its center: I = ML
^{2}∕12. - For rotation of an annulus of mass M, inner radius R
_{a}, and outer radius R_{b}about its axis of symmetry: I = M(R_{a}^{2}+ R_{ b}^{2})∕2.

If a rigid body is initially at rest, it will remain at rest if and only if the sum of
all the forces and the sum of all the torques acting on the body are zero. As an
example, a mass balance with arms of differing length is shown in figure 11.7. The
balance beam is subject to three forces pointing upward or downward, the tension
T in the string from which the beam is suspended and the weights M_{1}g and M_{2}g
exerted on the beam by the two suspended masses. The parameter g is the local
gravitational field and the balance beam itself is assumed to have negligible
mass. Taking upward as positive, the force condition for static equilibrium
is

| (11.25) |

Defining a counterclockwise torque to be positive, the torque balance computed about the pivot point in figure 11.7 is

| (11.26) |

where d_{1} and d_{2} are the lengths of the beam arms.

The first of the above equations shows that the tension in the string must be

| (11.27) |

while the second shows that

| (11.28) |

Thus, the tension in the string is just equal to the weight of the masses attached to the balance beam, while the ratio of the two masses equals the inverse ratio of the associated beam arm lengths.

- Show using the component form of the cross product given by equation (11.4) that A × B = -B × A.
- A mass M is sliding on a frictionless table, but is attached to a string which
passes through a hole in the center of the table as shown in figure 11.8. The
string is gradually drawn in so the mass traces out a spiral pattern as shown
in figure 11.8. The initial distance of the mass from the hole in the table is
R and its initial tangential velocity is v. After the string is drawn in, the
mass is a distance R′ from the hole and its tangential velocity is
v′.
- Given R, v, and R′, find v′.
- Compute the change in the kinetic energy of the mass in going from radius R to radius R′.
- If the above change is non-zero, determine where the extra energy
came from.

- A car of mass 1000 kg is heading north on a road at 30 m s
^{-1}which passes 2 km east of the center of town.- Compute the angular momentum of the car about the center of town when the car is directly east of the town.
- Compute the angular momentum of the car about the center of town when it is 3 km north of the above point.

- The apparatus illustrated in figure 11.9 is used to raise a bucket of mass M
out of a well.
- What force F must be exerted to keep the bucket from falling back into the well?
- If the bucket is slowly raised a distance d, what work is done on the bucket by the rope attached to it?
- What work is done by the force F on the handle in the above case?

- Derive equations (11.13) and (11.14).
- A mass M is held up by the structure shown in figure 11.10. The support
beam has negligible mass. Find the tension T in the diagonal wire. Hint:
Compute the net torque on the support beam about point A due to the
tension T and the weight of the mass M.

- A system consists of two stars, one of mass M moving with velocity
v
_{1}= (0,v, 0) at position r_{1}= (d, 0, 0), the other of mass 2M with zero velocity at the origin.- Find the center of mass position and velocity of the system of two stars.
- Find the spin angular momentum of the system.
- Find the internal kinetic energy of the system.

- A solid disk is rolling down a ramp tilted an angle θ from the horizontal. Compute the acceleration of the disk down the ramp and compare it with the acceleration of a block sliding down the ramp without friction.
- A mass M is suspended from the ceiling by two strings as shown in figure 11.11. Find the tensions in the strings.
- A man of mass M is a distance D up a ladder of length L which makes an angle θ with respect to the vertical wall as shown in figure 11.12. Take the mass of the ladder to be negligible. Find the force F needed to keep the ladder from sliding if the wall and floor are frictionless and therefore can only exert normal forces A and B on the ladder.

Harmonic Oscillator

Figure 12.1 illustrates the prototypical harmonic oscillator, the mass-spring system. A mass M is attached to one end of a spring. The other end of the spring is attached to something rigid such as a wall. The spring exerts a restoring force F = -kx on the mass when it is stretched by an amount x, i. e., it acts to return the mass to its initial position. This is called Hooke’s law and k is called the spring constant.

The potential energy of the mass-spring system is

| (12.1) |

which may be verified by noting that the Hooke’s law force is derived from this
potential energy: F = -d(kx^{2}∕2)∕dx = -kx. This is shown in figure 12.2. Since a
potential energy exists, the total energy E = K + U is conserved, i. e., is
constant in time. If the total energy is known, this provides a useful tool for
determining how the kinetic energy varies with the position x of the mass M:
K(x) = E - U(x). Since the kinetic energy is expressed (non-relativistically) in
terms of the velocity u as K = Mu^{2}∕2, the velocity at any point on the graph in
figure 12.2 is

| (12.2) |

Given all this, it is fairly evident how the mass moves. From Hooke’s law, the mass is always accelerating toward the equilibrium position, x = 0. However, at any point the velocity can be either to the left or the right. At the points where U(x) = E, the kinetic energy is zero. This occurs at the turning points

| (12.3) |

If the mass is moving to the left, it slows down as it approaches the left turning point. It stops when it reaches this point and begins to move to the right. It accelerates until it passes the equilibrium position and then begins to decelerate, stopping at the right turning point, accelerating toward the left, etc. The mass thus oscillates between the left and right turning points. (Note that equations (12.2) and (12.3) are only true for the harmonic oscillator.)

How does the period of the oscillation depend on the total energy of the
system? Notice that from equation (12.2) the maximum speed of the mass (i. e.,
the speed at x = 0) is equal to u_{max} = (2E∕M)^{1∕2}. The average speed must be
some fraction of this maximum value. Let us guess here that it is half the
maximum speed:

| (12.4) |

However, the distance d the mass has to travel for one full oscillation is twice the
distance between turning points, or d = 4(2E∕k)^{1∕2}. Therefore, the period of
oscillation must be approximately

| (12.5) |

The acceleration of the mass at any time is given by Newton’s second law:

| (12.6) |

An equation of this type is known as a differential equation since it involves a derivative of the dependent variable x. Equations of this type are generally more difficult to solve than algebraic equations, as there are no universal techniques for solving all forms of such equations. In fact, it is fair to say that the solutions of most differential equations were originally obtained by guessing!

We already have the basis on which to make an intelligent guess for the
solution to equation (12.6) since we know that the mass oscillates back and forth
with a period that is independent of the amplitude of the oscillation. A function
which might fill the bill is the cosine function. Let us try substituting
x = cos(ωt), where ω is a constant, into this equation. The second derivative of x
with respect to t is -ω^{2} cos(ωt), so performing this substitution results
in

| (12.7) |

Notice that the cosine function cancels out, leaving us with -ω^{2} = -k∕M. The
guess thus works if we set

| (12.8) |

The constant ω is the angular oscillation frequency for the oscillator, from
which we infer the period of oscillation to be T = 2π(M∕k)^{1∕2}. This agrees with
the earlier approximate result of equation (12.5), except that the approximation
has a numerical factor of 8 rather than 2π ≈ 6. Thus, the earlier guess is only off
by about 30%!

It is easy to show that x = B cos(ωt) is also a solution of equation (12.6),
where B is any constant and ω = (k∕M)^{1∕2}. This confirms that the oscillation
frequency and period are independent of amplitude. Furthermore, the sine
function is equally valid as a solution: x = A sin(ωt), where A is another constant.
In fact, the most general possible solution is just a combination of these two,
i. e.,

| (12.9) |

The values of A and B depend on the position and velocity of the mass at time t = 0. The right side of equation (12.9) shows an alternate way of writing the general harmonic oscillator solution that uses a cosine function with a phase factor ϕ.

If we wiggle the left end of the spring by the amount d = d_{0} cos(ω_{F }t), as in
figure 12.3, rather than rigidly fixing it as in figure 12.1, we have a forced
harmonic oscillator. The constant d_{0} is the amplitude of the imposed wiggling
motion. The forcing frequency ω_{F } is not necessarily equal to the natural or
resonant frequency ω = (k∕M)^{1∕2} of the mass-spring system. Very different
behavior occurs depending on whether ω_{F } is less than, equal to, or greater than
ω.

Given the above wiggling, the force of the spring on the mass becomes
F = -k(x - d) = -k[x - d_{0} cos(ω_{F }t)] since the length of the spring is the
difference between the positions of the left and right ends. Proceeding as for the
unforced mass-spring system, we arrive at the differential equation

| (12.10) |

The solution to this equation turns out to be the sum of a forced part in which x
is proportional to cos(ω_{F }t) and a free part which is the same as the solution to
the unforced equation (12.9). We are primarily interested in the forced part of the
solution, so let us set x = x_{0} cos(ω_{F }t) and substitute this into equation
(12.10):

| (12.11) |

Again the cosine factor cancels and we are left with an algebraic equation for x_{0},
the amplitude of the oscillatory motion of the mass.

Solving for the ratio of the oscillation amplitude of the mass to the amplitude
of the wiggling motion, x_{0}∕d_{0}, we find

| (12.12) |

where we have recognized that k∕M = ω^{2}, the square of the frequency of the free
oscillation. This function is plotted in figure 12.4.

Notice that if ω_{F } < ω, the motion of the mass is in phase with the
wiggling motion and the amplitude of the mass oscillation is greater than
the amplitude of the wiggling. As the forcing frequency approaches the
natural frequency of the oscillator, the response of the mass grows in
amplitude. When the forcing is at the resonant frequency, the response
is technically infinite, though practical limits on the amplitude of the
oscillation will intervene in this case — for instance, the spring cannot
stretch or shrink an infinite amount. In many cases friction will act to
limit the response of the mass to forcing near the resonant frequency.
When the forcing frequency is greater than the natural frequency, the
mass actually moves in the opposite direction of the wiggling motion —
i. e., the response is out of phase with the forcing. The amplitude of the
response decreases as the forcing frequency increases above the resonant
frequency.

Forced and free harmonic oscillators form an important part of many physical systems. For instance, any elastic material body such as a bridge or an airplane wing has harmonic oscillatory modes. A common engineering problem is to ensure that such modes are damped by friction or some other physical mechanism when there is a possibility of exitation of these modes by naturally occurring processes. A number of disasters can be traced to a failure to properly account for oscillatory forcing in engineered structures.

Complex exponential functions of the form x = exp(±iωt) also constitute solutions to the free harmonic oscillator governed by equation (12.6). This makes sense, as the complex exponential is the sum of sines and cosines. However, for the frictionless harmonic oscillator, the exponential solutions provide no particular advantage over sines and cosines. Furthermore, oscillator displacements are real, not complex quantities.

The superposition principle solves the problem of complex versus real solutions. For an equation like (12.6) which has real coefficients, if exp(iωt) is a solution, then so is exp(-iωt), so the superposition of these two solutions is also a solution. Furthermore

| (12.13) |

This shows a shortcut for getting the physical part of a complex exponential solution to equations like the harmonic oscillator equation; simply take the real part.

Complex exponential solutions come into their own for more complicated equations. For instance, suppose the force on the mass in the mass-spring system takes the form

| (12.14) |

The term containing b represents a frictional damping effect on the harmonic oscillator and the governing differential equation becomes

| (12.15) |

Trying the exponential function exp(σt) in this equation results in

| (12.16) |

where we have set

| (12.17) |

The quantity ω ≡ (ω_{0}^{2} - β^{2})^{1∕2} is the actual frequency of oscillation of the
damped oscillator, which one can see is less than the oscillation frequency ω_{0} that
occurs with the damping turned off. The physical solution to the damped
oscillator is thus

| (12.18) |

as long as ω_{0}^{2} > β^{2}. Notice that this solution is in the form of an oscillation
cos(ωt) multiplied by a decaying exponential exp(-βt). This confirms that the b
term decreases the amplitude of the oscillation with time.

The quantum mechanical harmonic oscillator shares the characteristic of other quantum mechanical bound state problems in that the total energy can take on only discrete values. Calculation of these values is too difficult for this book, but the problem is sufficiently important to warrant reporting the results here. The energies accessible to a quantum mechanical mass-spring system are given by the formula

| (12.19) |

In other words, the energy difference between successive quantum mechanical
energy levels in this case is constant and equals the classical resonant frequency
for the oscillator, ω = (k∕M)^{1∕2}, times ℏ.

- An oscillator (non-harmonic) has the potential energy function U(x) =
Cx
^{4}, where C is a constant. How does the oscillation frequency depend on energy? Explain your reasoning. - Show that C cos(ωt - ϕ) is an alternate way of writing A sin(ωt) + B cos(ωt) by finding the values of A and B in terms of the constants C and ϕ. Hint: Expand cos(ωt-ϕ) by using the trigonometric identity for the cosine of the sum of two angles.
- If a mass-spring harmonic oscillator has displacement x = 0 and velocity dx∕dt = V at time t = 0, determine the values of A and B as well as those of C and ϕ in equation (12.9).
- A mass M is suspended against gravity by a spring of spring constant k.
The unstretched length of the spring is x
_{0}and under the influence of gravity the spring is stretched to a resting length x_{1}> x_{0}.- Compute the length of the spring x
_{1}in the steady, resting case. - Set up the equation of motion for the mass moving under the influence of the two forces, gravity and spring. Solve the equation for the frequency of the oscillation and the position of the spring as a function of time x(t). Does the oscillation frequency change from the case without gravity?

- Compute the length of the spring x
- Determine the two real solutions to the damped harmonic oscillator problem
in the case in which ω
_{0}^{2}< β^{2}.

- Consider the pendulum in figure 12.5. The mass M moves along an arc
with x denoting the distance along the arc from the equilibrium
point.
- Find the component of the gravitational force tangent to the arc (and thus in the direction of motion of the mass) as a function of the angle θ. Use the small angle approximation on sin(θ) to simplify this answer.
- Get the force in terms of x rather than θ. (Recall that θ = x∕L.)
- Use Newton’s second law for motion in the x direction (i. e., along the arc followed by the mass) to get the equation of motion for the mass.
- Solve the equation of motion using the solution to the mass-spring problem as a guide.

- Forced damped oscillator:
- Add a damping term to the forced harmonic oscillator equation
(12.10) and solve for the forced part of the solution using complex
exponential methods. Hint: Change the cosine on the right side
of this equation to exp(iω
_{F }t) to convert the equation to complex form and then try the solution x = x_{0}exp(iω_{F }t) where x_{0}will depend on ω_{F }. Also, write the equation in terms of β = b∕(2M) and ω^{2}= k∕M. - Find the physical part of this solution by taking the real part
of x(t). Hint: While taking the real part of x, it may be helpful
to recall that the inverse of any complex number can be written
1∕(a + ib) = (a - ib)∕(a
^{2}+ b^{2}). - Determine how x
_{0}differs from that in the undamped case when ω_{F }is near the resonant frequency of the unforced oscillator. In particular, show how the phase and amplitude of the oscillation change as the forcing frequency changes from less than to greater than the resonant frequency.

- Add a damping term to the forced harmonic oscillator equation
(12.10) and solve for the forced part of the solution using complex
exponential methods. Hint: Change the cosine on the right side
of this equation to exp(iω
- A massless particle is confined to a box of length a. (Think of a photon between two mirrors.) Treating the particle classically, compute the period of one round trip from one end of the box to the other and back again. From this compute an angular frequency for the oscillation of this particle in the box. Does this frequency depend on the particle’s energy?
- Compute the ground state energy E
_{ground}of a massless particle in a box of length a using quantum mechanics. Compare E_{ground}∕ℏ with the angular frequency computed in the previous problem.

Constants, Units, and Conversions

This appendix contains various useful constants and conversion factors as well as information on the International System of Units.

“SI” is the French abbreviation for the International System of Units, the system used universally in science. See http://physics.nist.gov/cuu/Units/ for the last word on this subject. This treatment is derived from the National Institute of Science and Technology (NIST) website.

The most fundamental units of measure are length (meters; m), mass (kilograms; kg), time (seconds; s), electric current (ampere; A), temperature (kelvin; K), amount of a substance (mole; mol), and the luminous intensity (candela; cd). The candela is a rather specialized unit related to the perceived brightness of a light source by a “standard” human eye. As such, it is rather anthropocentric and hardly seems to merit the designation “fundamental”. The mole is also less fundamental than the other units, as it is simply a convenient way to refer to a multiple of Avogadro’s number of atoms or molecules.

Fundamental units can be combined to form derived units with special names. Some of these derived units are listed below.

Fundamental and derived SI units can have multipliers expressed as prefixes, e. g., 1 km = 1000 m. The NIST website points out a minor irregularity with the fundamental unit of mass, the kilogram. This already has the multiplier “kilo” prefixed to the unit “gram”. In this case 1000 kg is written 1 Mg, not 1 kkg, etc. SI multipliers are listed below as well.

Name | Abbrev. | Units | Meaning |

hertz | Hz | s^{-1} | frequency (cycles/sec) |

(unnamed) | s^{-1} | angular frequency (radians/sec) | |

newton | N | kg m s^{-2} | force |

pascal | Pa | N m^{-2} | pressure |

joule | J | N m | energy |

watt | W | J s^{-1} | power |

coulomb | C | A s | electric charge |

volt | V | N m C^{-1} | scalar potential |

(unnamed) | N s C^{-1} | vector potential | |

(unnamed) | V m^{-1} | electric field | |

tesla | T | N s C^{-1} m^{-1} | magnetic field |

(unnamed) | V m | electric flux | |

weber | Wb | T m^{2} | magnetic flux |

volt | V | V | electric circulation (EMF) |

(unnamed) | T m | magnetic circulation | |

farad | F | C V^{-1} | capacitance |

ohm | Ω | V A^{-1} | resistance |

henry | H | V s^{2} C^{-1} | inductance |

Multiplier | Name | Prefix |

10^{24} | yotta | Y |

10^{21} | zetta | Z |

10^{18} | exa | E |

10^{15} | peta | P |

10^{12} | tera | T |

10^{9} | giga | G |

10^{6} | mega | M |

10^{3} | kilo | k |

10^{2} | hecto | h |

10^{1} | deka | da |

10^{-1} | deci | d |

10^{-2} | centi | c |

10^{-3} | milli | m |

10^{-6} | micro | μ |

10^{-9} | nano | n |

10^{-12} | pico | p |

10^{-15} | femto | f |

10^{-18} | atto | a |

10^{-21} | zepto | z |

10^{-24} | yocto | y |

An older system of scientific units is the CGS system. This system is
still used widely in certain areas of physics. The fundamental units of
length, mass, and time are as implied by the title given above. The most
common CGS derived units are those for force (1 dyne = 10^{-5} N) and energy
(1 erg = 10^{-7} J).

Electromagnetism is expressed in several different ways in CGS units. Electromagnetic quantities in CGS not only have different units than in SI, they also have different physical dimensions, with different versions differing among themselves. The most common variant of CGS electromagnetic units is called “Gaussian” units. This variant is advocated by some physicists, though many others consider the whole subject of CGS electromagnetic units to be a terrible mess! SI units for electromagnetism are used in this text and CGS units will not be discussed further here.

1 lb = 4.448 N |

1 ft = 0.3048 m |

1 mph = 0.4470 m s^{-1} |

1 eV = 1.60 × 10^{-19} J |

1 mol = 6.022 × 10^{23} molecules |

(One mole of carbon-12 atoms has a mass of 12 g.) |

1 gauss = 10^{-4} T (CGS unit of magnetic field) |

1 millibar = 1 mb = 100 Pa (Old unit of pressure) |

When faced with solving an algebraic equation to obtain a numerical answer, solve the equation symbolically first and then substitute numbers. For example, given the equation

| (A.1) |

where a = 2 and b = 8, first solve for x,

| (A.2) |

and then substitute the numerical values:

| (A.3) |

This procedure is far better than substituting numbers first,

| (A.4) |

and then solving for x. Solving first and then substituting has two advantages: (1) It is easier to make algebraic manipulations with symbols than it is with numbers. (2) If you decide later that numerical values should be different, then the entire solution procedure doesn’t have to be repeated, only the substitutions at the end.

In numerical calculations, keep only one additional digit beyond those present in the least accurate input number. For instance, if you are taking the square root of 3.4, your calculator might tell you that the answer is 1.843908891. The answer you write down should be 1.84. Keeping all ten digits of the calculator’s answer gives a false sense of the accuracy of the result.

Round the result up if the digit following the last significant digit is 5 or greater and round it down if it is less than 5. Thus, the square root of 4.1, which the calculator tells us is 2.049390153, should be represented as 2.05 rather than 2.04.

It is easy to make mistakes when changing the units of a quantity. Adopting a systematic approach to changing units greatly reduces the chance of error. We illustrate a systematic approach to this problem with an example in which we change the units of acceleration from meters per second squared to kilometers per minute squared:

Symbol | Value | Meaning |

h | 6.63 × 10^{-34} J s | Planck’s constant |

ℏ | 1.06 × 10^{-34} J s | h∕(2π) |

c | 2.998 × 10^{8} m s^{-1} | speed of light |

G | 6.67 × 10^{-11} m^{3} s^{-2} kg^{-1} | universal gravitational constant |

k_{B} | 1.38 × 10^{-23} J K^{-1} | Boltzmann’s constant |

σ | 5.67 × 10^{-8} W m^{-2} K^{-4} | Stefan-Boltzmann constant |

K | 3.67 × 10^{11} s^{-1} K^{-1} | thermal frequency constant |

ϵ_{0} | 8.85 × 10^{-12} C^{2} N^{-1} m^{-2} | permittivity of free space |

μ_{0} | 4π × 10^{-7} N s^{2} C^{-2} | permeability of free space (= 1∕(ϵ_{
0}c^{2})). |

Symbol | Value | Meaning |

e | 1.60 × 10^{-19} C | fundamental unit of charge |

m_{e} | 9.11 × 10^{-31} kg = 0.511 MeV | mass of electron |

m_{p} | 1.672648 × 10^{-27} kg = 938.280 MeV | mass of proton |

m_{n} | 1.674954 × 10^{-27} kg = 939.573 MeV | mass of neutron |

Symbol | Value | Meaning |

M_{e} | 5.98 × 10^{24} kg | mass of earth |

M_{m} | 7.36 × 10^{22} kg | mass of moon |

M_{s} | 1.99 × 10^{30} kg | mass of sun |

R_{e} | 6.37 × 10^{6} m | radius of earth |

R_{m} | 1.74 × 10^{6} m | radius of moon |

R_{s} | 6.96 × 10^{8} m | radius of sun |

D_{m} | 3.82 × 10^{8} m | earth-moon distance |

D_{s} | 1.50 × 10^{11} m | earth-sun distance |

g | 9.81 m s^{-2} | earth’s surface gravity |

Absolute temperature, 6

Acceleration, 103, 104, 107, 137

centripetal, 106, 107

in relativity, 109

intrinsic, 109, 112

of wave packets, 137

Angle of incidence, 55

Angle of reflection, 55

Angular momentum, 153, 166, 189

conservation of, 193

definition, 192

location of origin, 192

orbital, 165

spin, 153, 168, 195

vector, 167

quantized length, 167

angular momentum

orbital, 195

Aristotle, 136

Atom, 119

Acceleration, 103, 104, 107, 137

centripetal, 106, 107

in relativity, 109

intrinsic, 109, 112

of wave packets, 137

Angle of incidence, 55

Angle of reflection, 55

Angular momentum, 153, 166, 189

conservation of, 193

definition, 192

location of origin, 192

orbital, 165

spin, 153, 168, 195

vector, 167

quantized length, 167

angular momentum

orbital, 195

Aristotle, 136

Atom, 119

Barrier penetration, 152, 163, 165

and nuclear fission, 165

and radioactive decay, 165

Beam, 14, 36–40, 43, 54, 67

diameter, 38

orientation, 39, 41

spread, 43

stationary, 36

width, 38

Beats, 14

beat frequency, 14

Black hole, 103

Bose, Satyendra Nath, 169

Boson, 169

Bragg diffraction, 119, 123

by a powder, 121

by a single crystal, 121

of electrons, 119, 123, 125

of X-rays, 119

Bragg’s law, 119, 121

test of quantum mechanics, 130

Calcite crystal, 54, 57, 58

Center of mass, 194

no relativistic generalization, 196

position, 194

velocity, 194

Centrifugal force, 111

Charge conjugation, 160

Circular motion, 105

Classical mechanics, 132

and geometrical optics, 141

Clocks in relativity, 82

Closed system, 180

Collisions, 176

elastic, 176

inelastic, 176, 179

Collocation, 78

reference frame dependence, 78

Compatible variables, 159, 167

and conservation, 160

Complex number, 153

absolute square, 155

complex conjugate, 155

graphical representation, 154

real and imaginary parts, 153

Confinement, 152

and potential energy, 161

in classical mechanics, 152

in quantum mechanics, 152

Contour plot, 35, 36

Conveyor belt, 182

Coordinate system, 30

Cartesian, 30, 31

choice of, 32

rotated, 31, 32

Crystal, 119

scattering of X-rays, 119

Davisson, Clinton, 123

De Broglie, Louis, 125

Nobel Prize address, 141

Derivative, 17

partial, 135, 143

Diffraction, 28, 119, 137

pattern, 28, 46

Diffraction grating, 28, 46, 119, 121

resolution, 48

Dispersion relation, 19, 34

anisotropic, 54, 57, 58

dispersive, 6, 20, 22, 23

isotropic, 36, 40, 57

non-dispersive, 6, 20

relativistic wave, 90, 95

Doppler shift, 90, 96, 114

light, 96

sound, 97

Double slit, 45, 124

Dynamics, 135

Newtonian, 135, 136

and geometrical optics, 143

pre-Newtonian, 135, 136

quantum mechanical, 135

Einstein, Albert, 72, 77, 103, 111, 123, 125, 137

Electromagnetic radiation, 119

Electromagnetism, 72

Electron, 119

as matter wave, 119, 123

particle and wave properties, 124

Energy, 119, 125

collisions, 176

conservation of, 138

and conservative forces, 138

kinetic, 138

non-relativistic, 129

of free particle, 156

potential, 137

total, 138

Energy level diagram, 162

Equivalence principle, 103, 111

Euler’s equation, 154

inverse, 155

Event, 73

Event horizon, 103, 115

Exhaust velocity, 181

Fermat’s principle, 63, 65–67

maximum time, 65

minimum time, 65

Fermi, Enrico, 169

Fermion, 169

Feynman, Richard, 66, 67, 163

view of quantum mechanics, 124

Force, 103, 107

central, 192

angular momentum conservation, 194

conservative, 135, 137, 143

inertial, 110

internal and external, 174

short range in collisions, 176

Four-momentum, 127

Four-vector, 90, 92

direction, 92

dot product, 92

magnitude, 92

position, 92

spacelike component, 92

timelike component, 92

wave, 92, 93, 96

Free particle, 156

in quantum mechanics, 156

wave function for, 156

Frequency, 3, 39, 90, 125

angular, 3, 4, 34

of wave, 3

rotational, 3

Galilei, Galileo, 72

Geiger counter, 123

Geometrical optics, 17, 54, 63

and classical mechanics, 132, 135

Germer, Lester, 123

Gravitational field, 113, 139

Gravitational red shift, 103, 113

Gravity, 103

as a conservative force, 139

as inertial force, 111

Group velocity, 19, 21, 125

and wave packets, 19, 96

of matter waves, 142

Harmonic oscillator, 204

angular frequency, 207

differential equation, 206

energy analysis, 205

forced, 207

mass-spring system, 204

Newton’s laws analysis, 206

quantum, 211

quantum frequencies, 211

resonance, 207

Heisenberg uncertainty principle, 119, 130

and wave packets, 130

position-momentum, 132

proper time-mass, 132

time-energy, 132

Hooke’s law, 204

spring constant, 204

Image, 61

real, 62, 63

virtual, 62, 63

Inclined ramp, 145

Index of refraction, 7, 16, 55, 56, 59

spatial variations, 141

Inertial force, 111

Interference, 7, 28, 137

constructive, 8, 16, 37, 46, 67, 121

destructive, 8, 16, 37, 46, 67

order, 46, 121

Interferometer, 14

Michelson, 14

optical, 14

thin film, 14, 15

Interval in spacetime, 81

meaning, 82

spacelike, 81

timelike, 81

Invariance, 159

and definiteness, 158

displacement in time, 158

under displacement, 156

under rotation, 167

up to a phase factor, 157

Kinematics, 135

of waves, 1

Kinetic energy, 128, 138, 142

inelastic collision, 179

internal, 195

non-relativistic, 129

translational, 195

Kinetic frequency, 128

Law of reflection, 64

Law of refraction, 64

Lens, 60

focal length, 61, 63

negative, 62

positive, 60, 61, 63

Light, 2

and photoelectric effect, 126

and soap bubbles, 17

and thin films, 15

dispersion relation, 34, 141

Doppler shift, 96

electromagnetic radiation, 2

in telescopes, 45

no special reference frame, 95

particle and wave properties, 123, 124

slope of world line, 76

speed in matter, 7

speed in vacuum, 7, 72, 91

transverse wave, 2

Line of simultaneity, 74, 79, 112

slope, 79, 85

Lorentz contraction, 84, 85

Mass, 103, 107, 125

and rest frequency, 127

in inelastic collision, 180

physical meaning, 126

Matter wave, 119

dispersion relation, 141

displacement of, 123

example of relativistic wave, 90

refraction, 145

theory of, 141

two and three dimensions, 144

Michelson, Albert, 14

Mirror, 55

concave, 63

convex, 63

curved, 60, 63

ellipsoidal, 66

focal length, 63

plane, 55

Moment of inertia, 197, 198

for solid bodies, 198

Momentum, 119, 125, 126

and group velocity, 138

and Newton’s second law, 173

collisions, 176

conservation of, 175

kinetic, 148

and particle velocity, 148

non-relativistic, 129

of free particle, 156

total (or canonical), 148

and wave vector, 126, 148

Nöther’s theorem, 152

Nöther, Emmy, 152

Newton’s first law, 136

Newton’s second law, 107, 136

accelerated reference frame, 111

and accelerating wave packets, 143

conservative force, 137

for open systems, 180

relativistic form, 173

rotational version, 192

Newton’s third law, 174

and Newton’s second law, 174

Newton, Isaac, 72

Nobel prize winners, 119, 120

Non-relativistic limit, 129

Open system, 180

conveyor belt, 182

mass equation, 181

Newton’s second law, 181

rocket, 180, 181

Parity, 160

Particle in a box, 152, 161

boundary condition, 165

non-relativistic, 162

ultra-relativistic, 162

Phase shift, 16

Phase speed, 4, 21, 23

Photoelectric effect, 126

Photon, 123, 126

Planck’s constant, 126, 142

Planck, Max, 125

Plane wave, 28, 33, 36

Potential energy, 137, 138, 142, 158

barrier, 165

gravitational, 139

Power, 135, 141

in three dimensions, 147

total, 141

Principle of relativity, 74, 77, 176

and Newtonian dynamics, 136

Prism, 59

Probability, 124

in quantum mechanics, 125

Probability amplitude, 125

Proper time, 81

meaning, 82

relation to spacetime interval, 81

Pulse, 36

isolated, 36

Quantization, 162

angular momentum, 166

energy, 162

Quantum mechanics, 14, 90, 119, 123, 132, 137

bizarreness of, 123

effect of measurement, 124

geometrical optics limit, 135

particles and waves, 129

probabilistic theory, 123

sense and nonsense, 124

two-slit interference, 124

uncertainty principle, 130

Quantum number, 162

angular momentum, 166, 168

energy, 162

orientation, 168

spin, 168

spin orientation, 168

Quark theory, 160

Ray, 54, 56, 58

Reference frame, 74, 75

accelerated, 86, 103, 107, 110, 112, 113

center of momentum, 176–179

for collisions, 176

inertial, 74, 77, 107

preferred, 74

lack of, 74

Reflection, 16, 137

in Bragg diffraction, 119

law of, 54

Refraction, 17, 137

law of (see Snell’s law), 54

Relativistic wave, 95, 127

dispersion relation, 95, 128

group velocity, 95, 128

phase speed, 95

Relativity, 72

Einsteinian, 77, 137

Galilean, 72, 74, 75, 77

general, 72, 103, 111

postulates, 78

special, 72, 75, 90

Rest energy, 127, 142

Rest frequency, 127

Rigid body, 198

Rocket, 181

acceleration, 182

exhaust velocity, 181

Scalar, 32

in spacetime, 93

Scattering, 119

Simultaneity, 78

Einsteinian relativity, 80

Galilean relativity, 78

reference frame dependence, 81

Sine wave, 2, 3, 8, 10, 28, 33, 36, 37

superposition, 12

Single slit, 44

Snell’s law, 56, 57, 59, 64, 65

and matter waves, 145

anisotropic case, 58

Spacetime, 73, 75

curved, 103

geometry of, 77

non-Euclidean geometry, 81

Spacetime diagram, 21, 74

elsewhere, 77

evolution, 74

future, 77

Spacetime interval, 78

spacelike, 76

timelike, 76

Spacetime Pythagorean theorem, 81, 82, 85

Speed of light, 72, 112

constancy, 77

in vacuum, 72

Spreading angle, 42–45

Statics, 198

Superposition principle, 7, 36

Surface tension, 6

Symmetry, 152

and conserved variables, 152

in quantum mechanics, 156

Telescope, 44

angular resolution, 44, 45

Thin lens approximation, 59

Thin lens equation, 61, 63

Thomson, G. P., 123

Thomson, J. J., 123

Time dilation, 83, 85

Time reversal, 160

Torque, 189

definition, 191

location of origin, 192

Total energy, 138, 142

conserved, 138

Total internal reflection, 57

Turning point, 139, 205

Twin paradox, 85

Uncertainty principle, 14

Vector, 28

absolute value, 30

addition, 30

components, 30, 31

cross product, 31, 189

component form, 190

does not commute, 190

magnitude, 190

right-hand rule, 190

direction, 30

displacement, 28

dot product, 31, 189

component form, 31, 33

cosine form, 31, 33

head, 30

magnitude, 30

tail, 30

unit vector, 30

Velocity, 104, 106

constant intrinsic acceleration, 112

Velocity addition

Galilean, 74

relativistic, 90, 99, 108

Wave, 1

longitudinal, 1

transverse, 1

Wave amplitude, 2, 7, 10, 42, 90

Wave crest, 3, 21

speed of movement, 23

Wave displacement, 2, 14, 21, 34, 42, 43, 137

Wave front, 33, 41, 54, 91

world line, 91

Wave function, 123, 125

collapse, 124

is complex, 123, 153

probability as absolute square, 123

Wave packet, 10, 23, 28, 36, 43, 54

acceleration of, 143

and group velocity, 19

and Newtonian dynamics, 137

isolated, 10, 40

of light, 141

speed of movement, 23

Wave period, 3

Wave phase, 4, 8, 15, 34, 91

Wave trough, 3, 21

Wave types, 4

gravitational waves, 34

gravity waves, 34, 35

light waves, 1, 7, 34, 35

ocean waves, 4, 35

deep water, 5, 34

shallow water, 5

sound waves, 1, 6

vibrations in solids, 1

Wave vector, 33, 34, 145

central, 42, 54

direction, 36

magnitude, 36

perpendicular, 92

Wave-particle duality, 125

Wavelength, 3, 15

Wavenumber, 3, 4, 10, 90, 125

central, 12, 19

difference, 10

imaginary, 164

spread, 12, 13

Work, 135, 140

by conservative force, 140

in three dimensions, 147

Work-energy theorem, 140

World line, 73, 75, 112

accelerated object, 113

slope, 74–76, 82

X-ray, 119, 121