Download Common Exam - 2010 Department of Physics University of Utah August 21, 2010

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts

Mathematical descriptions of the electromagnetic field wikipedia , lookup

Routhian mechanics wikipedia , lookup

Computational electromagnetics wikipedia , lookup

Relativistic quantum mechanics wikipedia , lookup

Transcript
Common Exam - 2010
Department of Physics
University of Utah
August 21, 2010
Examination booklets have been provided for recording your work and your solutions.
Please note that there is a separate booklet for each numbered question (i.e., use
booklet #1 for problem #1, etc.).
To receive full credit, not only should the correct solutions be given, but a sufficient
number of steps should be given so that a faculty grader can follow your reasoning.
Define all algebraic symbols that you introduce. If you are short of time it may be helpful
to give a clear outline of the steps you intended to complete to reach a solution. In some
of the questions with multiple parts you will need the answer to an earlier part in order to
work a later part. If you fail to solve the earlier part you may represent its answer with an
algebraic symbol and proceed to give an algebraic answer to the later part. This is a
closed book exam: No notes, books, or other records should be consulted. YOU MAY
ONLY USE THE CALCULATORS PROVIDED. The total of 250 points is divided
equally among the ten questions of the examination.
All work done on scratch paper should be NEATLY transferred to answer booklets.
SESSION 1
COMMON EXAM DATA SHEET
Physical constants and units
e = −1.602 × 10−19 C = −4.803 × 10−10 esu
1eV = 1.602 × 10−19 J
c = 299792458m · s−1 ≈ 3 × 108 m · s−1
h = 6.626 × 10−34 J · s = 6.626 × 10−27 erg · s = 4.135 × 10−15 eV · s
hc = 1.240 × 104 eV · Å
h̄ = 1.054 × 10−34 J · s = 1.054 × 10−27 erg · s = 6.582 × 10−16 eV · s
k = 1.380 × 10−23 J · K −1 = 1.380 × 10−16 erg · K −1 = 8.617 × 10−5 eV · K −1
G = 6.674 × 10−11 N · m2 · kg −2
g = 9.80 m · s−2
NA = 6.022 × 1023
!0 = 8.854 × 10−12 F · m−1
µ0 = 4π × 10−7 H · m−1
MElectron = 9.109 × 10−31 kg = 5.4858 × 10−4 amu = 511keV /c2
MP roton = 1.673 × 10−27 kg = 1.007276amu = 938.3M eV /c2
MN eutron = 1.675 × 10−27 kg = 1.008665amu = 939.5M eV /c2
1mile = 1609 m
1m = 3.28 f t
Some useful integrals
a, b, n > 0
!
!
√ x
dx
a+bx
√1
dx
x2 a+bx
! (a+bx)n/2
x
=−
√
√
−2(2a−bx) a+bx
2
3b
a+bx
ax
x
!
=
x
!
!
√
a − bx2 −
√ 1
dx
a2 +x2
! √a2 +x2
dx
!
!
=
√
b
2a
√
"
= ln ""x +
!
√1
dx
x a+bx
! (a+bx)(n−2)/2
x
√
"
a2 + x2 ""
"
"
= arcsin xa
1
dx
(x2 ±a2 )3/2
=
ñx
a2 x2 ±a2
1
dx
(a2 −x2 )3/2
=
√x
a2 a2 −x2
√
sin2 xdx = 12 (x − sin x cos x)
2
e−x dx =
! ∞ n−1 −x
x e dx
0
! ∞ n+1 −a2 x2
x e
dx
=
√
π
= Γ(n)
1
2an+2
#
Γ
"
a2 +x2 "
"
x
cos2 xdx = 12 (x + sin x cos x)
!∞
dx
"√ √
"
" a+ a−bx2 "
"
x
a ln "
a2 + x2 − a ln " a+
√ 1
dx
a2 −x2
−∞
0
−
dx = n2 (a + bx)n/2 + a
! √a−bx2
dx
!
=
n+2
2
$
Maxwell’s equations (SI unit system)
∇·D
=
ρ
∇·B
=
0
∇ × E = − ∂B
∂t
∇×H
=
J+
∂D
∂t
Maxwell’s equations (Gaussian unit system)
∇·D
=
4πρ
∇·B
=
0
∇ × E = − 1c ∂B
∂t
∇×H
=
4π
1 ∂D
J+
c
c ∂t
Simplest material equations
D = $·E
B = µ·H
J = σ·E
Problem 1 – General physics
Consider a cylindrical bowl of radius r = 10 cm and depth ! = 5 mm. The bowl sits on a horizontal
surface with its axis vertical, and is filled with black ink at an initial temperature T = 300 K. The dish is
then exposed to sunlight. The sun can be considered an ideal Black Body and at the time of measurement
the Sun is ψ = 70◦ above the horizon. Only η = 75% of the light striking the top of the atmosphere gets
through to the surface of the earth, and the sun has an angular radius θ = 0.26◦ . After ∆t = 60 s seconds,
you observe that the temperature of the ink has increased by ∆T = 1.35 K. The heat capacity of the ink is
Cp = 4185.5 J · kg −1 · K −1 and its density is ρ = 1000 kg · m−3 . The sides and bottom of the bowl are coated
with a perfectly reflecting material and the top surface of the ink radiates as a ideal Black Body and there is
no non-radiative heat transfer. The Stefan-Boltzmann constant, σ = 5.67 × 10−8 W m−2 K −4 .
(a) [8 pts.] Write a differential equation for the ink temperature as a function of time (Hint: θ ≈ RSun /DSun ).
(b) [7 pts.] From the provided data, estimate the surface temperature TSun of the Sun.
(c) [4 pts.] Estimate the ink temperature at which thermal equilibrium is reached. (Note: if you were not
able to answer the previous question, use TSun = 6000K as the Sun temperature). Give two reasons why
such a high temperature cannot be reached in practice with the setup described.
(d) [6 pts.] At night, the same setup is prepared with the ink at an initial temperature T1 = 320 K. How
much time does it take for the water to cool to T2 = 300 K?
Problem 1 – General physics – sample solution
Consider a cylindrical bowl of radius r = 10 cm and depth ! = 5 mm. The bowl sits on a horizontal
surface with its axis vertical, and is filled with black ink at an initial temperature T = 300 K. The dish is
then exposed to sunlight. The sun can be considered an ideal Black Body and at the time of measurement
the Sun is ψ = 70◦ above the horizon. Only η = 75% of the light striking the top of the atmosphere gets
through to the surface of the earth, and the sun has an angular radius θ = 0.26◦ . After ∆t = 60 s seconds,
you observe that the temperature of the ink has increased by ∆T = 1.35 K. The heat capacity of the ink is
Cp = 4185.5 J · kg −1 · K −1 and its density is ρ = 1000 kg · m−3 . The sides and bottom of the bowl are coated
with a perfectly reflecting material and the top surface of the ink radiates as a ideal Black Body and there is
no non-radiative heat transfer. The Stefan-Boltzmann constant, σ = 5.67 × 10−8 W m−2 K −4 .
(a) [8 pts.] Write a differential equation for the ink temperature as a function of time (Hint: θ ≈ RSun /DSun ).
Let Q be the heat contained in the ink. Then
dQ
dt
= πr2 !ρCp dT
= −πr2 σT 4 + η
dt
2
4
4πRSun
σTSun
2
4πDSun
· πr2 sinψ
where we have introduced RSun and DSun the radius and distance of the Sun. After simplification and
identifying θ ≈ RSun /DSun it comes
dT
dt
=
σ
4
(ηθ2 sinψTSun
#ρCp
− T 4)
(b) [7 pts.] From the provided data, estimate the surface temperature TSun of the Sun.
From the data we have, we see that, at first,
dT
dt
=
1.35
60
= 0.0225 Ks−1 . Solving the above equation for
TSun we get
p dT
1
TSun = [ ηθ2 sinψ
( #ρC
+ T 4 )]1/4 . The angle θ = 0.26◦ = 4.538 × 10−3 rad so numerically we get
σ dt
TSun = [ 0.75×(4.538×101−3 rad)2 ×sin70◦ ( 5×10
−3 m×1000kg·m−3 ×4185.5J·kg −1 K −1
5.67×10−8 W ·m−2 ·K −4
× 0.0225K · s−1 + (300K)4 )]1/4 =
5798 K
(c) [4 pts.] Estimate the temperature at which thermal equilibrium is reached. (Note: if you were not able
to answer the previous question, use TSun = 6000K as the Sun temperature). Give two reasons why such
a high temperature cannot be reached in practice with the setup described.
Thermal equilibrium corresponds to
dT
dt
= 0 and using the differential equation, T = (ηθ2 sinψ)1/4 TSun .
Numerically T = (0.75 × (4.538 × 10−3 )2 × sin70◦ )1/4 TSun = 6.172 × 10−2 × TSun . So for TSun = 5798K
we get T = 358K and for TSun = 6000K we get T = 370K.
These temperatures are close to the boiling water temperature under standard atmospheric conditions.
In practice such high temperatures can not be reached in the setup of this exercise because of heat loss
by thermal conduction and because of evaporation, however solar water heater are designed to prevent
water from boiling in the Sun exposed conduits.
(d) [6 pts.] At night, the same setup is prepared with the ink at an initial temperature T1 = 320 K. How
much time does it take for the water to cool to T2 = 300 K?
At night, the heat source disappears from our differential equation which becomes
−A · T 4 with A =
T1
σ
.
"ρCp
−[ 3T1 3 ]T = A∆t so ∆t =
2
Numerically, ∆t =
This differential equation can be integrated as −
"ρCp
(T2−3
3σ
− T1−3 ).
5×10−3 m×103 kg·m−3 ×4185.5J·kg −1 K −1
((300K)−3
3×5.67×10−8 W ·m−2 ·K −4
! T1 dT
T2
T4
dT
dt
σ
= − "ρC
T4 =
p
= −
! ∆t
0
− (320K)−3 ) = 802 s = 13.3 min
A · dt or
Problem 2 – Electrodynamics
(a) [12 pts.] Neglecting edge effects, find the capacitance per unit area of a parallel plate capacitor with
infinite metallic plates in the planes x = 0 and x = d. The capacitor is filled with a media with
permittivity:
!(x) = !0
(x + d)
.
d
(1)
(b) [13 pts.] Neglecting edge effects, find the total capacitance of a parallel plate capacitor with metallic
plates in the planes x = 0 and x = d. The sizes of the plates in z and y directions are Z and Y
respectively. The capacitor is filled with a media with permittivity:
!(y, z) = !0 (1 + sin2 (2πzm/Z) sin2 (2πyn/Y )),
where n and m are integers.
(2)
Problem 2 – Electrodynamics – sample solution
(a) [12 pts.] Neglecting edge effects, find the capacitance per unit area of a parallel plate capacitor with
infinite metallic plates in the planes x = 0 and x = d. The capacitor is filled with a media with
permittivity:
!(x) = !0
(x + d)
.
d
(1)
" and displacement D
" obey the following equations within the capacitor:
Electric field E
" = 0
divD
(2)
" = 0
curlE
(3)
" = !(x)E
"
D
(4)
There is also the boundary condition that may be written in a form
−
!
" · d"s = V,
E
(5)
where V is the potential difference between the two metallic plates and the integral may be taken along
any contour connecting these plates.
" and D
" have only x-components, Dx = D0 ,
This system of equations is satisfied if we assume that both E
where D0 is x-independent, and Ex ≡ E(x) = D0 /!(x).
As follows from the Gauss theorem the surface charge on metallic plates σ = ±D0 and the capacitance
per unit area is
C=
On the other hand,
V =
! d
0
D0
.
V
E(x)dx = D0
(6)
! d
0
dx
!(x)
(7)
Thus,
C=
!0
.
d ln(2)
(8)
(b) [13 pts.] Neglecting edge effects, find the total capacitance of a parallel plate capacitor with metallic
plates in the planes x = 0 and x = d. The sizes of the plates in z and y directions are Z and Y
respectively. The capacitor is filled with a media with permittivity:
!(y, z) = !0 (1 + sin2 (2πzm/Z) sin2 (2πyn/Y )),
(9)
where n and m are integers.
!
Due to the fact that E · dl = 0, the tangential component of the electric field is zero inside the capacitor.
Since ! is independent of x , the normal component of the electric field is constant: E = V /d. On the
other hand, Dn (z, y) = σ(z, y) = E!(z, y). Therefore, the total charge is
Q=
" Y " Z
0
0
#
$
%
$
%
" Y
" Z
V !0
2πzm
2 2πyn
σ dydz =
YZ +
sin
dy
sin2
dz
d
Y
Z
0
0
Since, Q = Ct V , where Ct is the total capacitance we obtain
Ct =
5 !0 Y Z
.
4 d
&
Problem 3 – Classical Mechanics
A uniform cylindrical log of mass M = 5 kg and radius R = 10 cm is initially at rest on the ground when
it is hit by a bullet of mass m = 10 g, moving horizontally (and also perpendicularly to the log’s axis O) with
a velocity v = 1000 m/s as shown in the picture. The bullet lodges very close to the axis O.
(a) [10 pts.] Assuming that no slipping occurs between the log and the ground, find the frequency ν of the
log’s rotation after the impact.
(b) [7 pts.] What fraction of the initial mechanical energy is dissipated into heat during the collision?
(c) [8 pts.] What is the total impulse
!
F dt provided by the friction force between the log and the ground
during the collision (as compared with the initial momentum of the bullet)?
Problem 3 – Classical Mechanics – sample solution
A uniform cylindrical log of mass M = 5 kg and radius R = 10 cm is initially at rest on the ground when
it is hit by a bullet of mass m = 10 g, moving horizontally (and also perpendicularly to the log’s axis O) with
a velocity v = 1000 m/s as shown in the picture. The bullet lodges very close to the axis O.
(a) [10 pts.] Assuming that no slipping occurs between the log and the ground, find the frequency ν of the
log’s rotation after the impact.
The momentum of inertia of a uniform cylinder with respect to its axis is IO = M R2 /2 and with respect
to the point on its rim, IG = IO + M R2 = 3M R2 /2. Angular momentum conservation about point G
gives,
3
mvR = M R2 ω + mRω,
2
which gives the angular frequency,
ω=
mv
.
(3M/2 + m)R
The frequency is
ν=
ω
mv
mv
=
≈
≈ 2.1 Hz.
2π
2π(3M/2 + m)R
3πM R
(b) [7 pts.] What fraction of the initial mechanical energy is dissipated into heat during the collision?
The kinetic energy after the collision is
(IG + mR2 )
ω2
mv 2
m
mv 2 2m
mv 2
=
≈
≈ 0.0013
.
2
2 3M/2 + m
2 3M
2
About 99.9% of the initial kinetic energy is converted into heat.
(c) [8 pts.] What is the total impulse
!
F dt provided by the friction force between the log and the ground
during the collision (as compared with the initial momentum of the bullet)?
One solution is to equate the impulse of the friction force to the decrease of the total momentum.
Alternatively, it is related to the angular momentum increase with respect to point O:
R
"
F dt = IO ω =
M R2
mv
mvR
≈
.
2 (3M/2 + m)R
3
Thus the total impulse provided by the friction force in the collision is close to 1/3 of the initial momentum
mv.
Problem 4 – Quantum Mechanics
A particle of mass m moves in a one-dimensional repulsive delta-function potential U (x) = αδ(x), with
α > 0.
(a) [10 pts.] Derive the boundary conditions for the wave function ψ(x) and its derivative dψ/dx at the
origin. (Hint: integrate the Schrödinger equation over x in the interval (−d, d), and take the limit
d → +0 after that.)
(b) [10 pts.] Find the reflection R and transmission T coefficients for the particle with energy E > 0
incoming on the delta-function barrier from the left.
(c) [5 pts.] Discuss how your results for R and T apply to the attractive potential with α < 0. Also discuss
E → ∞ limit of your solution: is it consistent with classical expectations?
!!""!#
x
FIG. 1: The scattering problem.
Problem 4 – Quantum Mechanics – sample solution
A particle of mass m moves in a one-dimensional repulsive delta-function potential U (x) = αδ(x), with
α > 0.
(a) [10 pts.] Derive the boundary conditions for the wave function ψ(x) and its derivative dψ/dx at the
origin. (Hint: integrate the Schrödinger equation over x in the interval (−d, d), and take the limit
d → +0 after that.)
The Schrödinger equation reads
−
h̄2 d2 ψ(x)
+ αδ(x)ψ(x) = Eψ(x).
2m dx2
(1)
Integrate both sides over x from some −d to d, and take the limit d → +0. Then
dψ(+d)/dx − dψ(−d)/dx
! +d
−d
→
! +d
−d
dx d2 ψ(x)/dx2 =
ψ " (+0) − ψ " (−0), where prime denotes spatial derivative, and
dx δ(x)ψ(x) = ψ(0) irrespective of d, while the right-hand side
as d → 0. Thus we obtain
ψ " (+0) − ψ " (−0) =
! +d
−d
dxEψ(x) → 0 approaches zero
2mα
ψ(0).
h̄2
(2)
The wave function itself is always continuous, ψ(+0) = ψ(−0).
(b) [10 pts.] Find the reflection R and transmission T coefficients for the particle with energy E > 0
incoming on the delta-function barrier from the left.
For x #= 0 Eq.(1) reads, for E > 0,
d2 ψ(x)
2mE
=
−
Eψ(x),
dx2
h̄2
so that its general solution is ψ(x) = a1 eikx + b1 e−ikx , with k =
(3)
"
2mE/h̄2 .
We now set up the scattering problem:
ψ(x) = eikx + Ce−ikx for x < 0 , ψ(x) = Beikx for x > 0.
(4)
Transmission coefficient D = |B|2 while reflection R = |C|2 . Using derived boundary conditions we
obtain two equations for unknown B, C
1 + C = B ; ikB − ik(1 − C) =
2mα
(1 + C).
h̄2
(5)
Solving for C we get
C=
Hence
!
R = |C|2 = 1 +
! ikh̄2
"−1
−1
.
(6)
h̄4 k 2 "−1
1
=
.
2
2
2
mα
1 + 2h̄ E/(mα2 )
(7)
mα
Clearly D = 1 − R.
(c) [5 pts.] Discuss how your results for R and T apply to the attractive potential with α < 0. Also discuss
E → ∞ limit of your solution: is it consistent with classical expectations?
The result, Eq.(7), is not sensitive to the sign of α: hence the well (α < 0) reflect the particle as well as
the barrier (α > 0) does!
In the limit E → ∞: R(E) →
mα2
2h̄2 E
→ 0 while D(E) → 1 −
mα2
2h̄2 E
→ 1. Thus at high energy reflection
disappears and all particles pass through, which is quite reasonable classically.
!!""!#
x
FIG. 1: The scattering problem.
Problem 5 – Thermodynamics
The equilibrium energy of a photon gas in a cavity with a volume V and temperature T is given by the
equation E = αV T 4 . Useful equations for this problem are dE = T dS − P dV and dF = −SdT − pdV .
(a) [9 pts.] Find the entropy of the photon gas.
(b) [9 pts.] Find the pressure of the photon gas under the same conditions.
(c) [7 pts.] The volume of the cavity is expanded adiabatically up to V2 . What is the equilibrium temperature
of the photon gas? Neglect the heat capacity of the cavity.
Problem 5 – Thermodynamics – sample solution
The equilibrium energy of a photon gas in a cavity with a volume V and temperature T is given by the
equation E = αV T 4 . Useful equations for this problem are dE = T dS − P dV and dF = −SdT − pdV .
(a) [9 pts.] Find the entropy of the photon gas.
If the volume of the cavity is kept constant, from the thermodynamic identity dE = T dS − P dV we
have T dS = dE = 4V αT 3 dT . From the third law of thermodynamics, the entropy of a system at zero
temperature is zero. Then, after integration, we have
S=
! S
0
dS =
! T
0
4
4V αT 2 dT = αV T 3 .
3
(1)
(b) [9 pts.] Find the pressure of the photon gas under the same conditions.
The free energy of the photon gas is
1
F = E − T S = − αV T 4 .
3
(2)
Then, the pressure of the photon gas is
"
∂F
P =−
∂V
#
T
1
= αT 4
3
(3)
(c) [7 pts.] The volume of the cavity is expanded adiabatically up to V2 . What is the equilibrium temperature
of the photon gas? Neglect the heat capacity of the cavity.
In an adiabatic expansion of the photon gas the entropy is conserved. That means V1 T13 = V2 T23 or
T2 = T1
$
V1
V2
% 13
.
(4)
Common Exam - 2010
Department of Physics
University of Utah
August 21, 2010
Examination booklets have been provided for recording your work and your solutions.
Please note that there is a separate booklet for each numbered question (i.e., use
booklet #1 for problem #1, etc.).
To receive full credit, not only should the correct solutions be given, but a sufficient
number of steps should be given so that a faculty grader can follow your reasoning.
Define all algebraic symbols that you introduce. If you are short of time it may be helpful
to give a clear outline of the steps you intended to complete to reach a solution. In some
of the questions with multiple parts you will need the answer to an earlier part in order to
work a later part. If you fail to solve the earlier part you may represent its answer with an
algebraic symbol and proceed to give an algebraic answer to the later part. This is a
closed book exam: No notes, books, or other records should be consulted. YOU MAY
ONLY USE THE CALCULATORS PROVIDED. The total of 250 points is divided
equally among the ten questions of the examination.
All work done on scratch paper should be NEATLY transferred to answer booklets.
SESSION 2
COMMON EXAM DATA SHEET
Physical constants and units
e = −1.602 × 10−19 C = −4.803 × 10−10 esu
1eV = 1.602 × 10−19 J
c = 299792458m · s−1 ≈ 3 × 108 m · s−1
h = 6.626 × 10−34 J · s = 6.626 × 10−27 erg · s = 4.135 × 10−15 eV · s
hc = 1.240 × 104 eV · Å
h̄ = 1.054 × 10−34 J · s = 1.054 × 10−27 erg · s = 6.582 × 10−16 eV · s
k = 1.380 × 10−23 J · K −1 = 1.380 × 10−16 erg · K −1 = 8.617 × 10−5 eV · K −1
G = 6.674 × 10−11 N · m2 · kg −2
g = 9.80 m · s−2
NA = 6.022 × 1023
!0 = 8.854 × 10−12 F · m−1
µ0 = 4π × 10−7 H · m−1
MElectron = 9.109 × 10−31 kg = 5.4858 × 10−4 amu = 511keV /c2
MP roton = 1.673 × 10−27 kg = 1.007276amu = 938.3M eV /c2
MN eutron = 1.675 × 10−27 kg = 1.008665amu = 939.5M eV /c2
1mile = 1609 m
1m = 3.28 f t
Some useful integrals
a, b, n > 0
!
!
√ x
dx
a+bx
√1
dx
x2 a+bx
! (a+bx)n/2
x
=−
√
√
−2(2a−bx) a+bx
2
3b
a+bx
ax
x
!
=
x
!
!
√
a − bx2 −
√ 1
dx
a2 +x2
! √a2 +x2
dx
!
!
=
√
b
2a
√
"
= ln ""x +
!
√1
dx
x a+bx
! (a+bx)(n−2)/2
x
√
"
a2 + x2 ""
"
"
= arcsin xa
1
dx
(x2 ±a2 )3/2
=
ñx
a2 x2 ±a2
1
dx
(a2 −x2 )3/2
=
√x
a2 a2 −x2
√
sin2 xdx = 12 (x − sin x cos x)
2
e−x dx =
! ∞ n−1 −x
x e dx
0
! ∞ n+1 −a2 x2
x e
dx
=
√
π
= Γ(n)
1
2an+2
#
Γ
"
a2 +x2 "
"
x
cos2 xdx = 12 (x + sin x cos x)
!∞
dx
"√ √
"
" a+ a−bx2 "
"
x
a ln "
a2 + x2 − a ln " a+
√ 1
dx
a2 −x2
−∞
0
−
dx = n2 (a + bx)n/2 + a
! √a−bx2
dx
!
=
n+2
2
$
Maxwell’s equations (SI unit system)
∇·D
=
ρ
∇·B
=
0
∇ × E = − ∂B
∂t
∇×H
=
J+
∂D
∂t
Maxwell’s equations (Gaussian unit system)
∇·D
=
4πρ
∇·B
=
0
∇ × E = − 1c ∂B
∂t
∇×H
=
4π
1 ∂D
J+
c
c ∂t
Simplest material equations
D = $·E
B = µ·H
J = σ·E
Problem 6 – Classical Mechanics
Two bodies with masses m1 and m2 are connected by a spring with a spring constant k. The equilibrium
distance between the two bodies is L.
(a) [8 pts.] Find the frequency of oscillation of the system.
Body 2 is now charged with charge +Q and body 1 is neutral and non-polarizable. The system is initially
at rest with the distance between the two bodies being L. At time t = 0, an electric field E is turned
on. The direction of the field is along the axis of the system as shown in the figure.
(b) [7 pts.] Write down Newton’s equations of motion for the first and the second bodies.
(c) [10 pts.] Find the distance between two bodies as a function of time.
Problem 6 – Classical Mechanics – sample solution
Two bodies with masses m1 and m2 are connected by a spring with a spring constant k. The equilibrium
distance between the two bodies is L.
(a) [8 pts.] Find the frequency of oscillation of the system.
The position of the center of mass of the system can be found from the equations m1 L1 = m2 L2 and
L = L1 + L2 , where L1 and L2 are the distances from the center of mass to the first and the second body,
respectively. We find that L2 =
m2 +m1
L.
m1
Since the center of mass is immobile we may describe the oscillation of the second body as if it is
1
connected to a spring with a length L2 and a spring constant k2 = k LL2 = k m2m+m
. The frequency of
1
oscillations of the system is
!
m2 + m1
ω= k
m1 · m2
"1
2
.
(1)
Body 2 is now charged with charge +Q and body 1 is neutral and non-polarizable. The system is initially
at rest with the distance between the two bodies being L. At time t = 0, an electric field E is turned
on. The direction of the field is along the axis of the system as shown in the figure.
(b) [7 pts.] Write down Newton’s equations of motion for the first and the second bodies.
m2 x¨2 = eQ − k(x2 − x1 − L),
(2)
m1 x¨1 = k(x2 − x1 − L),
(3)
where x2 and x1 are the coordinates of the second and the first body, respectively.
(c) [10 pts.] Find the distance between two bodies as a function of time.
Let S be the deviation of the distance between bodies 1 and 2 from their equilibrium separation L,
S = x2 − x1 − L. Then from eq. (2) and (3) we find an equation of motion for S.
S̈ = −ω 2 S + β,
where ω is given by eq. (1) and β =
(4)
eQ
.
m2
The equation (4) is an inhomogeneous ordinary differential equation. Its solution can be found as a sum
of a general solution of a homogeneous equation and any solution of an inhomogeneous equation. For
example, it can be chosen as
S = A cos(ωt) + B sin(ωt) +
β
ω2
(5)
From initial conditions S(t = 0) = 0 and Ṡ(t = 0) = 0, we find that B = 0 and A = − ωβ2 .
We have then finally
x2 (t) − x1 (t) = L −
eQ
eQ
cos(ωt) +
.
2
m2 ω
m2 ω 2
(6)
Problem 7 – Electrodynamics
Consider some properties of inhomogeneous plane waves that are also called evanescent waves. These waves
appear in the theory of total internal reflection and in many other examples of near field optics.
Suppose that the electric field of the evanescent wave has only a z component. It has the form
Ez = E0 exp(iky y − κx − iωt),
(1)
where κ, ky are real. The wave propagates in a dielectric medium with electric permittivity # and magnetic
permeability µ.
$ as given by Eq.(1) obeys the wave equation of electrodynamics and find ω(κ, ky ).
(a) [3 pts.] Show that E
(b) [2 pts.] Find the region of (κ, ky ), where the frequency ω is real.
$
(c) [10 pts.] Using the Maxwell equations, find all components of the magnetic field H.
(d) [10 pts.] Prove that if A(t) = A0 exp −iωt and B(t) = B0 exp −iωt, then < ReAReB >= Re(A∗ · B)/2,
where < ... > means time averaging. Using this relation, find all components of the time averaged
$ =< [ReE
$ × ReH]
$ >.
Poynting vector S
Problem 7 – Electrodynamics – sample solution
Consider some properties of inhomogeneous plane waves that are also called evanescent waves. These waves
appear in the theory of total internal reflection and in many other examples of near field optics.
Suppose that the electric field of the evanescent wave has only a z component. It has the form
Ez = E0 exp(iky y − κx − iωt),
(1)
where κ, ky are real. The wave propagates in a dielectric medium with electric permittivity # and magnetic
permeability µ.
$ as given by Eq.(2) obeys the wave equation of electrodynamics and find ω(κ, ky ).
(a) [3 pts.] Show that E
Substituting
Ez = E0 exp(iky y − κx − iωt),
(2)
(∂ 2 /c2 ∂ 2 t − ∇2 )Ez = 0
(3)
into equation
one gets ω 2 = c2 (ky2 − κ2 ).
(b) [2 pts.] Find the region of (κ, ky ), where the frequency ω is real.
The frequency ω is real if ky2 ≥ κ2 .
$
(c) [10 pts.] Using the Maxwell equations, find all components of the magnetic field H.
Using Maxwell equation
$ = −∂ B/∂t
$
curlE
(4)
$ = iω B,
$
i[$k × E]
(5)
one finds
$ = µH.
$ Then
where $k has projections (iκ, ky , 0) and B
ky
Ez
ωµ
iκEz
Hy = −
ωµ
Hx =
(6)
(d) [10 pts.] Prove that if A(t) = A0 exp −iωt and B(t) = B0 exp −iωt, then < ReAReB >= Re(A∗ · B)/2,
where < ... > means time averaging. Using this relation, find all components of the time averaged
" =< [ReE
" × ReH]
" >.
Poynting vector S
To prove the relations one should write ReA = (A + A∗ )/2 and similar for B. Then take into account
that < exp(±2iωt) >= 0
" are
The components of the Poynting vector S
Sx = −Hy Ez
(7)
Sy = Hx Ez
(8)
Sz = 0
(9)
After time averaging one gets
< Sx >= −Re(Hy Ez∗ )/2 =
κ
Re(i|Ez |2 ) = 0.
2ωµ
(10)
and
< Sy >=
ky
ky
|Ez |2 =
|E0 |2 exp (−2κx)
2ωµ
2ωµ
Thus, the energy of the evanescent wave propagates only in the direction with real momentum.
(11)
Problem 8 – Statistical Mechanics
A classical gas of atoms, each of mass m, is maintained at temperature T inside an enclosure. The gas is
described by Maxwell-Boltzmann distribution,
f (!v )d3 v =
!
m "3/2
m!v 2 3
exp[−
]d v.
2πkB T
2kB T
(1)
The atoms emit light which passes through a window of the enclosure and can then be observed as a spectral
line in a spectroscope. A stationary atom would emit light at a sharply defined wavelength λ0 . But, because of
the Doppler effect, the observed wavelength λ of the light from an atom with the component vs of the velocity
along the line of sight is instead approximately given by
λ = λ0 (1 +
vs
),
c
(2)
where c is the speed of light. Note that components of the velocity orthogonal to the line of sight do not affect
the wavelength of the radiation.
As a result, the light arriving at the spectroscope has an intensity in the wavelength range (λ, λ + dλ)
characterized by some intensity distribution I(λ)dλ.
(a) [15 pts.] Derive the distribution function I(λ)dλ of the observed light for a gas of N atoms.
(b) [5 pts.] Observe that Eq.(2) implies that |vs | " c. What restriction does this impose on the temperature
T of the gas?
(c) [5 pts.] Using the result of (b), normalize the distribution by requiring
#∞
0
I(λ)dλ = N I0 . (Here I0 is
the intensity of radiation detected by the spectroscope from a stationary atom.) (Hint: This question
requires making an approximation which relies on the result of question (b).)
Problem 8 – Statistical Mechanics – sample solution
A classical gas of atoms, each of mass m, is maintained at temperature T inside an enclosure. The gas is
described by Maxwell-Boltzmann distribution,
f (!v )d3 v =
!
m "3/2
m!v 2 3
exp[−
]d v.
2πkB T
2kB T
(1)
The atoms emit light which passes through a window of the enclosure and can then be observed as a spectral
line in a spectroscope. A stationary atom would emit light at a sharply defined wavelength λ0 . But, because of
the Doppler effect, the observed wavelength λ of the light from an atom with the component vs of the velocity
along the line of sight is instead approximately given by
λ = λ0 (1 +
vs
),
c
(2)
where c is the speed of light. Note that components of the velocity orthogonal to the line of sight do not affect
the wavelength of the radiation.
As a result, the light arriving at the spectroscope has an intensity in the wavelength range (λ, λ + dλ)
characterized by some intensity distribution I(λ)dλ.
(a) [15 pts.] Derive the distribution function I(λ)dλ of the observed light for a gas of N atoms.
We only need the distribution of the vs component of the velocity as the transverse ones do not affect the
wavelength of the emitted light. Hence, choosing for a moment vs to be along the x-axis, we integrate
over the two orthogonal to it velocity components, vy and vz , in the full range (−∞, ∞) of allowed values.
Note that !v
2
= vs2 + vy2 + vz2 and d3 v = dvs dvy dvz . Then, the number of particles with the s-component
of their velocity in the interval (vs , vs + dvs ) is obtained
dN (vs ) = N
!
m "1/2
mvs2
exp[−
]dvs .
2πkB T
2kB T
(3)
Using (2) we convert this into the number of atoms dN (λ) emitting light in the interval (λ, λ + dλ):
dN (λ) = N
c ! m "1/2
(λ − λ0 )2
exp[−
]dλ,
λ0 2πkB T
δ2
(4)
where δ 2 = 2kB T λ20 /(mc2 ).
The distribution in question follows as
I(λ)dλ = αdN (λ)
(5)
where the unknown normalization constant α is to be fixed by the normalization condition. The distribution is centered around λ0 and its width is given by δ. Observe that the distribution is very sharp:
λ0 /δ =
!
mc2 /(2kB T ) ! 1.
(b) [5 pts.] Observe that Eq.(2) implies that |vs | " c. What restriction does this impose on the temperature
T of the gas?
The average velocity of particles is determined by the temperature alone: v̄ ∼
velocity of the s-component is also given by this expression, v̄s ∼
!
!
kB T /m. Thus average
kB T /m. Hence, vs " c translates
into kB T " mc2 : the thermal energy of the particles in much smaller than their rest energy. The motion
is non-relativistic.
Observe that this condition implies the sharpness of the derived light distribution: its width is much
smaller than the average wavelength, δ/λ0 =
!
2kB T /(mc2 ) " 1.
(c) [5 pts.] Using the result of (b), normalize the distribution by requiring
"∞
0
I(λ)dλ = N I0 . (Here I0 is
the intensity of radiation detected by the spectroscope from a stationary atom.) (Hint: This question
requires making an approximation which relies on the result of question (b).)
Normalization requires α
"∞
0
dN (λ) = N I0 , which leads to
αN # ∞
αN # ∞ −(λ−λ0 )2 /δ2
2
N I0 = √
e
dλ = √
e−x dx.
πδ 0
π −λ0 /δ
(6)
Here we substituted λ = λ0 + δx. Since, for the non-relativistic gas λ0 /δ ! 1, we can replace the
lower limit of the integration by negative infinity, −∞, which results in the standard gaussian integral
√
" ∞ −x2
dx = π. Thus we fix the normalization constant as α = I0 .
−∞ e
Problem 9 – Quantum mechanics
Consider a single 1/2 spin that initially is aligned with the positive direction of the y-axis. At time t = 0, a
magnetic field along the positive direction of the z-axis is suddenly switched on. Interaction of the spin with
this field is described by the Hamiltonian,
Ĥ = −bσ̂z ,
b>0
(a) [9 pts.] Find the wavefunction of the spin at a moment of time t > 0.
(b) [4 pts.] Describe qualitatively the spin dynamics for t > 0.
(c) [5 pts.] Will there be an occasion when the spin will be aligned with the positive direction of the x-axis?
If this is the case, when will it happen?
(d) [4 pts.] What about the spin being aligned with the negative direction of the x-axis? When will it
happen?
(e) [3 pts.] What about the spin being aligned with the the z-axis? When will it happen if it ever does?
For reference, Pauli matrices:


0 1
σ̂x = 
,
1 0


 0 −i 
σ̂y = 
,
i 0


1 0 
σ̂z = 
.
0 −1
Problem 9 – Quantum mechanics – sample solution
Consider a single 1/2 spin that initially is aligned with the positive direction of the y-axis. At time t = 0, a
magnetic field along the positive direction of the z-axis is suddenly switched on. Interaction of the spin with
this field is described by the Hamiltonian,
Ĥ = −bσ̂z ,
b>0
(a) [9 pts.] Find the wavefunction of the spin at a moment of time t > 0.
The eigenstates of the Hamiltonian correspond to two states with sz = ±h̄/2:


1
χsz =+h̄/2 =   e−ibt/h̄ ,
0


0
χsz =−h̄/2 =   eibt/h̄ .
1
(1)
At t < 0 the wave function corresponds to sy = +h̄/2. This eigenfunction of σ̂y is simply


1 1
χsy =+h̄/2 = √  
2 i
(2)
The linear combination of the functions (1) that coincides with (2) at t = 0 is




1 1
i 0
χ(t) = √   e−ibt/h̄ + √   eibt/h̄ .
2 0
2 1
(3)
(b) [4 pts.] Describe qualitatively the spin dynamics for t > 0.
The motion of spin is precession around the z-axis with the angular frequency b/h̄.
(c) [5 pts.] Will there be an occasion when the spin will be aligned with the positive direction of the x-axis?
If this is the case, when will it happen?
Spin will be aligned with the positive direction of the x-axis when its wavefunction becomes


1
χsx =+h̄/2 = c   .
1
From (4) we obtain that it happens when simultaneously
eibt/h̄ = c,
ie−ibt/h̄ = c,
→ e2ibt/h̄ = i,
which happens at
t=
h̄
h̄n
+
,
4b
b
n = 0, 1, 2, ...
(4)
(d) [4 pts.] What about the spin being aligned with the negative direction of the x-axis? When will it
happen?
Similarly, spin is aligned along the negative direction of the x-axis at times shifted by half-period compared to the above,
t=
3h̄ h̄n
+
,
4b
b
n = 0, 1, 2, ...
(e) [3 pts.] What about the spin being aligned with the the z-axis? When will it happen if it ever does?
As long as dissipation of energy is absent, the spin is never going to be along the z-axis. In practice
however, energy dissipation, although it may be small, is always present and at t − ∞ the spin is going
to get aligned with the positive direction of the z-axis, which is its minumum energy (ground) state.
Problem 10 – Modern physics
Consider photons γ with frequency ν = 160 GHz and density n = 2 × 108 m−3 propagating in the negative
x-direction and protons propagating in the opposite direction. The protons interact with the photons as
p + γ → ∆(1232) → p + π 0 with a total cross section σ = 300µb. The mean lifetime of the ∆(1232) resonance
is τ∆ = 5 × 10−25 s. (m∆ = 1232 M eV /c2 , mp = 938 M eV /c2 and mπ0 = 135 M eV /c2 , 1 b = 10−28 m2 ).
(a) [6 pts.] What is the minimum energy protons must have to produce ∆(1232) resonances when interacting
with the photons? You may assume the proton to be ultra-relativistic and verify with your result that
this assumption is well justified.
(b) [6 pts.] A proton of energy E = 1021 eV interacts with a photon, producing a ∆(1232). What is the
mean distance traveled by the ∆(1232) before it decays?
(c) [7 pts.] From the point of view of a proton of energy Ep = 1021 eV , what is the wavelength of the
photons? How does this wavelength compare to the physical size of a proton?
(d) [6 pts.] In the laboratory rest frame, what is the average distance traveled by E = 1021 eV protons
without interacting? Express your result in light years (a light year is the distance travelled at the speed
of light in one year).
Problem 10 – Modern physics – sample solution
Consider photons γ with frequency ν = 160 GHz and density n = 2 × 108 m−3 propagating in the negative
x-direction and protons propagating in the opposite direction. The protons interact with the photons as
p + γ → ∆(1232) → p + π 0 with a total cross section σ = 300µb. The mean lifetime of the ∆(1232) resonance
is τ∆ = 5 × 10−25 s. (m∆ = 1232 M eV /c2 , mp = 938 M eV /c2 and mπ0 = 135 M eV /c2 , 1 b = 10−28 m2 ).
(a) [6 pts.] What is the minimum energy protons must have to produce ∆(1232) resonances when interacting
with the photons? You may assume the proton to be ultra-relativistic and verify with your result that
this assumption is well justified.
For the reaction to be possible, the mass invariant must be (Ep + Eγ )2 − (cPp − cPγ )2 > m2∆ c4 or
Ep2 + Eγ2 + 2Ep Eγ − c2 Pp2 − c2 Pγ2 + 2Pp Pγ c2 > m2∆ c4 which simplifies to
m2p c4 + 2Ep Eγ + 2Pp cEγ > m2∆ c4 and considering the proton to be ultra-relativistic, Ep ≈ cPp so:
Ep >
m2∆ c4 −m2p c4
4hν
and, numerically: Ep >
(1232 M eV )2 −(938 M eV )2
4×6.63×10−34 J·s×160×109 Hz/1.602×10−19 J·eV −1
= 2.4 × 1020 eV
Since we find Ep >> mp c2 for this reaction to be possible, the ultra-relativistic approximation is well
justified.
(b) [6 pts.] A proton of energy E = 1021 eV interacts with a photon, producing a ∆(1232). What is the
mean distance traveled by the ∆(1232) before it decays?
The lifetime of the resonance is going to be affected by time dilation so we need to find the Lorentz
factor of ∆(1232). We can do this by using conservation of energy: Ep + Eγ = γm∆ c2 so γ =
1021 eV
1232×106 eV
Ep +Eγ
m∆ c2
≈
= 8.117 × 1011 . Since γ >> 1, the ∆(1232) can be considered to be going at the speed of light
and the distance it travels is D = γ · τ∆ · c = 8.117 × 1011 × 5−25 s × 3 × 108 m · s−1 = 0.12 mm.
(c) [7 pts.] From the point of view of a proton of energy Ep = 1021 eV , what is the wavelength of the
photons? How does this wavelength compare to the physical size of a proton?
The Lorentz factor of the protons is γ =
Ep
mp c2
=
1021 eV
938×106 eV
= 1.066 × 1012 . The energy of the photons in
the laboratory reference frame is Eγ = h · ν. The energy of the photons in the proton rest frame can be
obtained by a Lorentz transformation: Eγ# = γEγ + βγEγ ≈ 2γEγ as since γ >> 1 we have β ≈ 1. So
hc/λ# = 2γhν and λ# =
Numerically λ# =
c
.
2γν
3×108 m·s−1
1.066×1012 ×160×109 Hz
as the physical size of a proton.
= 8.8 × 10−16 m = 0.88 f m. This wavelength is of the same order
(d) [6 pts.] In the laboratory rest frame, what is the average distance traveled by E = 1021 eV protons
without interacting? Express your result in light years (a light year is the distance travelled at the speed
of light in one year).
The average distance D traveled without interacting is the speed v of the protons multiplying the mean
time before interaction so D =
D ≈
1
.
2σ·n
Numerically, D =
v
.
σ·n(c+v)
We have seen the protons are ultra-relativistic so v ≈ c and
1
2×300×10−6 ×10−28 m2 ×2×108 m−3
= 8.3 × 1022 m. A light year corresponds to
365 × 24 × 3600 s × 3 × 108 m · s−1 = 9.46 × 1015 m and so the mean distance traveled by protons without
interacting is 8.8 × 106 ly.