Lecture 10
Lecture 10
Marat Siddikov
December 6, 2019
Outline Today’s plan:
Finalize parametric oscillations:
? Parametric resonance* (advanced
topic, see Landau, Section 27)
Motion of the system of particles.
Collisions. Scattering. Cross-
section
Technical notes
Please feel free to interrupt/ask
questions if needed
In case of connection problems, please
use your email for comments
I’ve uploaded the .PDF- file of this
lecture to Aula Virtual
Bibliography:
Symon→Marion→Landau, Fetter-
Walecka
Announcement
I’ve put a short Sample Problem to Aula Virtual. The goal
of the problem below is to let You practice with nonlinear and
forced oscillations which we have seen this week (and for me
your solutions would provide important feedback if we can pro-
ceed to other topics or should revisit “oscillations”). It is op-
tional assignment, not a regular homework (we cannot have
evaluations in December according to the letter of Vice Rec-
tor). If You have questions, please feel free to ask
The mathematical pendulum of length `
oscillates with small amplitude θmax
1 near its equlibrium. You need
to develop a perturbation theory with
4
. account of ∼ O θ term in the
lagrangian for the case of free oscilla-
tions near equilibrium and for forced os-
cillations (see the .pdf file in Aula Vir-
tual for more details)
Parametric resonance (example we know from childhood)
Up to now we’ve considered a system which If the child seats calmly, the amplitude de-
oscillates freely or under effects of known ex- creases due to dissipative forces (air resistance)
ternal force If the child moves in a certain way, the ampli-
There are systems where external force is not tude might increase despite dissipative forces.
given explicitly but instead the parameters of Formally, child’s motion changes parameters
the system change in a predetermined way of the system (center of mass position `eff (t),
air resistance)
Example from our childhood
Formally, child’s motion changes parameters
Parametric resonance
of the system (center of mass position `eff (t),
Up to now we’ve considered a system which air resistance)
oscillates freely or under effects of known
external force
There are systems where external force is d
m `2 (t)φ̇ + β(t) `eff (t) φ̇ +
not given explicitly but instead changes eff
dt | {z }
| {z }
parameters of the system Mz vφ
+ m g `eff (t) φ = 0
φ̈ + 2α(t)φ̇ + ω 2 (t) φ = 0
ξ̈ + Ω2 (t) ξ = 0 (1)
2 2 2
Ω (t) = ω (t) − α (t) − α̇(t)
Parametric resonance
ξ̈ + Ω2 (t) ξ = 0 (1)
Up to now we’ve considered a system which
oscillates freely or under effects of known Ω2 (t) = ω 2 (t) − α2 (t) − α̇(t)
external force
For arbitrary Ω(t) can integrate only numer-
There are systems where external force is
ically.
not given explicitly but instead changes
In many applications the t-dependence is
parameters of the system
weak,
φ̈ + 2α(t)φ̇ + ω 2 (t) φ = 0
Z t
φ(t) → exp − dτ α(τ ) ξ(t)
Parametric resonance Z
Up to now we’ve considered a system which ξE (t) + dt G (t − τ )fE (τ ),
oscillates freely or under effects of known
external force
There are systems where external force is fE (t) = −f (t)Ω20 Ae i Ω0 t + A∗ e −i Ω0 t
not given explicitly but instead changes
parameters of the system Case of periodic perturbation:
X
φ̈ + 2α(t)φ̇ + ω 2 (t) φ = 0 f (t) = f˜n e i ωn t ,
Z t n
φ(t) → exp − dτ α(τ ) ξ(t) 2πn
ωn = , fn∗ = f−n
T
ξ̈ + Ω2 (t) ξ = 0 (1) X iω t iΩ t
fE (t) = Ω20 f˜n e n Ae 0 + A∗ e −i Ω0 t ,
2 2 2
Ω (t) = ω (t) − α (t) − α̇(t) n
(n)
Near resonance frequencies it fails (all ξF
ξ(t) = Ae i Ω0 t + A∗ e −i Ω0 t + ξE (t), are large). Need to consider in detail what hap-
pens close to resonance:
ξ̈ + Ω2 (t) ξ = 0 (1)
a = a0 e st , a∗ = a0∗ e st
2 2 2
Ω (t) = ω (t) − α (t) − α̇(t)
Nontrivial solutions only if
Case of weak t-dependence:
s 2 + isω − Ω0 f0 Ω20
2 det ∗ 2 2 = 0,
Ω (t) = Ω20 (1 + f (t)) , ∀t, f (t) 1. f 0 Ω0 s − isω − Ω0
2
ω2
Need to consider in detail what happens s 2 + Ω20 − + s 2 ω 2 − |f0 |2 Ω40 = 0.
4
close to resonance:
ω2
q
2
:= ω − 2Ω0 , || 2Ω0 s1,2 = −Ω20 − ± |f0 |2 Ω40 + ω 2 Ω20 .
4
Which sign should we take ?
Seek for periodic solutions:
Parametric resonance
Up to now we’ve considered a system
which oscillates freely or under effects of ξF (t) = a(t)e i ωt/2 +a∗ (t)e −i ωt/2 , ȧ Ω a
known external force
There are systems where external force is
not given explicitly but instead changes a = a0 e st , a∗ = a0∗ e st
parameters of the system
Nontrivial solutions only if
Z t
s 2 + isω − Ω0 f0 Ω20
φ(t) → exp − dτ α(τ ) ξ(t)
det = 0,
f0∗ Ω20 2
s − isω − Ω0
ξ̈ + Ω2 (t) ξ = 0 (1)
ω2
q
2 2 2
Ω (t) = ω (t) − α (t) − α̇(t) s 2 = −Ω20 − + |f0 |2 Ω40 + ω 2 Ω20
4
Case of weak t-dependence:
|f0 |2 Ω20
≈ − 2 .
Ω2 (t) = Ω20 (1 + f (t)) , ∀t, f (t) 1.
4
Exponential increase of ξ(t) for || < |f0 |Ω0 /2
Need to consider in detail what happens Physical
amplitude increases
if e st grows faster
than exp − t dτ α(τ ) ;
R
close to resonance:
Fcase α = const:
:= ω − 2Ω0 , || 2Ω0
|f0 |2 Ω20
− α2 > 2
4
Seek for periodic solutions:
Parametric resonance
Up to now we’ve considered a system
which oscillates freely or under effects of ξF (t) = a(t)e i ωt/2 +a∗ (t)e −i ωt/2 , ȧ Ω a
known external force
There are systems where external force is
not given explicitly but instead changes a = a0 e st , a∗ = a0∗ e st
parameters of the system
Exponential increase of ξ(t) for || < |f0 |Ω0 /2
t Physical
amplitude increases if e st grows faster
Z
φ(t) → exp − dτ α(τ ) ξ(t) Rt
than exp − dτ α(τ ) ;
Fcase α = const:
ξ̈ + Ω2 (t) ξ = 0 (1)
2 2 2 |f0 |2 Ω20
Ω (t) = ω (t) − α (t) − α̇(t) − α2 > 2
4
Case of weak t-dependence:
Bibliography:
Symon→Marion→Landau, Fetter-
Walecka
Case of closed two-particle system masses
System of N particles m1 and m2 interacting with each other.:
Center of mass
P Equations of motion in terms of the vari-
~ CM := Pmk~rk .
R
ables
k mk
~r12 = ~r1 − ~r2 , ~ = m1~r1 + m2~r2 :
R
m1 + m2
~
d 2R
Mtot ~ (ext) = 0,
=F
dt 2
d 2~r12 ~ 12
µ12 =F
dt 2
Reduced mass
~ CM
dP ~ (ext)
=F
−1
dt m1 m2 1 1
µ12 = = +
m1 + m2 m1 m2
~ CM = P ~pk := Mtot R
P ~˙ CM is the total
k
momentum ofP the system Relative motion with reduced (“effective”)
~ (ext) =
F ~ (ext) is the total force mass
k Fk
acting on the system.
⇒If F ~ (ext) = 0 (the system is “closed”),
~
then P = const.
System of particles Angular momentum and the center of mass
~ CM (t) :
motion ~rk → ~rk − R
Angular momentum of the system
X X
~ →
M ~rk × ~pk + ~ CM × ~pk
R
X X
~ :=
M ~rk × ~pk = mk ~rk × ~vk k k
| {z } | {z }
k k ~ int ~ CM ×P
~ CM
M R
Total momentum
~˙ CM
X
~ CM =
P ~pk := Mtot R
k
We will call such scatterings elastic.
Inelastic scattering
Total angular momentum
We’ll call process inelastic if:
X X There are changes in internal state
~ :=
M ~rk × ~pk = mk ~rk × ~vk (=masses) of particles, e.g.
k k ?Particles stick together in 2 → 1
process
Energy of the system ?One of the particles disintegrates
into n pieces in 1 → n process
X m ~v 2 X X Mechanical energy is not conserved (i.e.
k
E = + Uk (~rk )+ Uik (~ri − ~rk ) might partially convert into other forms)
2
k k i <k ?Interaction is not potential.
?Can use conservation of momen-
tum and angular momentum to
reconstruct kinematics
2 → 2 elastic scattering
Assume that scattering is elas- Asymptotically (r12 → ∞) particles are free,
tic, no change of internal properties move along straight lines
(masses), no new particle production We call scattering angles θ1 , θ2 the angle
Interaction is potential, U = between the asymptotes of the particle before
U (~r1 − ~r2 ) , and assume and after interaction. Note that angles θ1 6= θ2
in general:
lim U(r ) = 0, r = |~r1 − ~r2 |
r →∞
sin θ1 = sin θ2 ⇒ θ1 = θ2
Remark on 2 → 2 inelastic scattering
In the inelastic scattering the inter- Scattering angles θ1 , θ2 the angle between
nal state of particle might change, so the asymptotes of the particle before and after
kinetic energy is not conserved (con- interaction.
verts to/from internal energy). Momentum conservation:
Interaction is potential, U = ~p1 I + ~p2 I = ~p1 F + ~p2 F
U (~r1 − ~r2 ) , and assume
(Kinetic) energy conservation
lim U(r ) = 0, r = |~r1 − ~r2 | ~p1 I 2 ~p2 I 2 ~p1 F 2 ~p2 F 2
r →∞ + +U = + +U+∆Einternal
2 m1 2 m2 2 m1 2 m2
dN = J dS dt
Scattering in the external field
Flux of the particles Differential cross-section
Number of particles which pass The number of particles scattered in the
through the unit area in transverse range of angles (θ, θ + d θ) and divided by
direction per unit of time the flux
dN = J dS dt dσ 1 dNscattered (θ, θ + d θ)
=
dθ J dθ
dσ 1 dNscattered 1 dNscattered
= =
dΩ J dΩ J sin θ dθ dφ
dΩ = sin θ dθdφ
dσ 1 dNscattered (E 0 , E 0 + dE 0 )
0
=
dE J d E0
Scattering in the external field
Flux of the particles Differential cross-section
Number of particles which pass dσ 1 dNscattered (θ, θ + d θ)
through the unit area in transverse =
dθ J dθ
direction per unit of time
dσ 1 dNscattered (E 0 , E 0 + dE 0 )
0
=
dN = J dS dt dE J d E0
θ = π − 2φ0 ,
∞ M
dr
Z
m r2
φ0 = r
M2
rmin 2
m
E − U(r ) − 2 m r2
where
m v∞2
E = , M = m v∞ ρ,
2
and thus
dσ dρ
= 2πρ(θ)
dθ dθ
Scattering in the Coulomb field