IvanovMaksyuta
.pdfЛекція 11. Задача двох тіл. Розсіювання частинок
§ 1 Багаточастинкові задачі з ньютонівською взаємодією
Розглядаючи задачу Кеплера, ми вважали одне із взаємодіючих тіл нерухомим. Однак здебільшого це не так, зокрема, коли кілька тіл взаємодіють між собою й кожне з них рухається під дією інших. Задача про рух таких тіл називається багаточастинковою задачею. У переважній більшості випадків її точні розв'язки невідомі. Переходимо до найпростіших багаточастинкових задач.
1. Задача двох тіл. Розглянемо дві матеріальні точки з масами m1 і m2 , які взаємодіють одна з одною. У цьому випадку згід-
но з третім законом Ньютона потенціальна енергія залежить лише від відстані між ними, тобто від абсолютної величини різниці їх радіус-векторів r1 і r2 (рис. 1). Тому функція Лагранжа
такої системи записується у вигляді
|
G2 |
|
G2 |
|
|
|||
L(rG1, rG1, rG2 , rG2 ) = |
m1r1 |
+ |
m2r2 |
−U ( |
|
rG1 − rG2 |
|
) . |
|
|
|||||||
2 |
2 |
|
|
|||||
|
|
Вводимо тепер радіус-вектор центра інерції системи двох тіл
RG = m1rG1 + m2rG2 m1 + m2
інерції, тобто переходимо в систему координат ( x′, y′, z′) . Щодо
системи центра інерції (ц-система) функція Лагранжа і радіусвектор центра інерції записуються у вигляді
G G G G |
m rG′2 |
|
m rG′2 |
|
|
G G |
|
G |
|
|
m rG′+ m rG′ |
|
|||
′ ′ ′ |
′ |
1 1 |
+ |
2 2 |
|
|
′ |
′ |
|
|
′ |
= |
1 1 |
2 2 |
= 0 . |
|
|
|
|||||||||||||
L(r1, r1, r2 |
, r2 ) = |
2 |
2 |
−U ( |
|
r1 |
− r2 |
|
), R |
|
m |
+ m |
|||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
Вводячи вектор взаємної відстані між матеріальними точками |
|||||||
rG = rG1 − rG2 = rG1′− rG2′ і враховуючи |
m1r1′+ m2rG2′ = 0 , знаходимо раді- |
||||||
ус-вектори rG1′ і r2′: |
|
|
|
|
|
|
|
G |
m2 |
G |
G |
|
m1 |
G |
|
r1′= |
|
r, |
r2′ = − |
|
|
r . |
|
m1 + m2 |
m1 |
+ m2 |
|||||
|
|
|
|
||||
|
|
111 |
|
|
|
З урахуванням цих обчислень функція Лагранжа тепер набуває такого вигляду:
L(rG, rG) = μrG2 −U (r) , 2
де μ = m1m2 (m1 + m2 ) – зведена маса. Отже, одержана функція Лагранжа формально описує рух деякої фіктивної точки з масою μ і радіус-вектором r у зовнішньому центральному полі U (r) . Цю фіктивну точку називають μ -точкою. Рух μ -точки репрезентує відносний рух реальних матеріальних точок. Знайшовши радіус-вектор r , можна знайти радіус-вектори r1′ і r2′ обох матеріальних точок у ц-системі за формулами
G |
G |
μ G |
G |
G |
μ G |
||
r1 |
= R + |
|
r , |
r2 |
= R − |
|
r |
|
|
||||||
|
|
m1 |
|
|
m2 |
у лабораторній системі координат ( x, y, z) (л-система). Із формул rG1′= μrGm1 , r2′ = −μrGm2 видно, що радіус-вектори частинки в ц-системі направлені протилежно, а їхні довжини обернено пропорційні масам ( r1′ r2′ = m2 m1 ), тобто радіус-вектор частин-
ки з більшою масою менший, і навпаки.
Якщо траєкторія руху μ -точки є еліпсом, то і реальні частин-
ки рухаються еліптичними орбітами, причому обидва еліпси мають спільний фокус, розташований y центрі інерції (рис. 2).
Наприклад, для Сонця і Землі m1 m2 ≈ 3 10−6 , тобто Сонце відхи-
ляється від центра мас на відстань 450 км, що дуже мало з астрономічної точки зору. Тому можна вважати, що Сонце практично розташоване у фокусі.
Перший закон Кеплера можна сформулювати таким чином: кожна планета Сонячної системи рухається по еліпсу, в одному із фокусів якого розташоване Сонце. В узагальненій формі (на випадок параболи і гіперболи) цей закон формулюється так: у ньютонівському полі матеріальна точка рухається конічним перетином.
112
z |
|
|
z′ |
|
|
|
|
m1 |
r1′ |
|
y′ |
m1 |
|
|
|
|
|
|||
G |
G |
G |
r ′ |
r2′ |
||
r1 |
R |
|
r′ |
m |
1 |
|
|
|
|
2 |
2 |
|
m2 |
|
x′ |
r2 |
|
y |
|
|
|
|
|
|
|||
x |
|
|
|
|
Рис. 2. Схематичне зображення |
|
Рис. 1. Координатні системи |
еліптичних траєкторій |
|||||
двох взаємодіючих матеріальних |
||||||
до задачі двох тіл |
тіл з різними масами |
Закон збереження моменту імпульсу L = mρ2ϕ = const допускає геометричну інтерпретацію, а саме, вираз ρ2dϕ2 являє со-
бою площу сектора (рис. 3), що утворюється двома нескінченно близькими радіус-векторами та елементом дуги траєкторії. По-
значаючи її як df , запишемомомент частинки у вигляді L = 2mf , де похідну f називають сек-
торною швидкістю. Тому збе-
реження моменту вказує на постійність секторної швидкості. Другий закон Кеплера можна сформулювати так: за рівні проміжки часу радіус-вектор Землі описує рівні площі.
Інтегруючи тепер рівність
ρdϕ
ρ
dϕ
Рис. 3. Площа, яка описується радіус-вектором за час dt
L = 2mf за часом від нуля до
періоду T , одержуємо 2mf = TL . Для еліпса площа f = πab і, враховуючи, що
a = |
|
p |
= |
α |
, |
b = |
p |
= |
L |
, |
||||
|
− e2 |
2 |
|
E |
|
|
||||||||
1 |
|
|
|
|
1 − e2 |
|
2m |
E |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
знаходимо T = πα m 2 E 3 , а з урахуванням знову виразу для ве-
ликої півосі a і того, що α = GMm , одержуємо T = 2πa3/2 GM .
113
Звідси можна сформулювати третій закон Кеплера: відношення куба великої півосі еліпса планети до квадрата її періоду обертання навколо Сонця є величиною сталою (вона не залежить від маси планети m , тобто однакова для всіх планет Сонячної системи).
2. Задача трьох тіл. Загальний розв'язок задачі трьох тіл, який виражається через відомі функції, знайти не вдається. У кінці ХІХ ст. А. Пуанкаре і Х. Брунс показали, що такого розв'язку у вигляді алгебраїчних або однозначних трансцендентних функцій від координат швидкостей взаємодіючих тіл не існує. Знайдено загальний розв'язок цієї задачі у вигляді рядів, але ці ряди збігаються дуже повільно. Були знайдені деякі частинні розв'язки задачі трьох тіл. У 1772 р. Ж. Л. Лагранж показав, що рівняння руху будуть справджуватися, якщо до двох тіл, що обертаються, додати третє, помістивши його в одну із п'яти точок, які називаються точками лібрації (лат. libra – терези). Колінеарні точки лібрації L1, L2 , L3 , які зображені на рис. 4, а, є точками нестійкої рівнова-
ги. Трикутні точки лібрації L4 , L5 утворюють з двома тілами рів-
носторонні трикутники (рис. 4, б); рівноважне положення третього тіла в цих точках стійке.
У трикутних точках лібрації системи Сонце–Юпітер розташовані дві групи астероїдів, які були названі – Греки і Троянці. У трикутних точках лібрації системи Земля–Місяць польський астроном К. Кордилевський у 1961 р. виявив скупчення метеорного пилу.
L4
С |
|
С |
L3 |
L1 |
L2 |
|
|
L5 |
а) |
|
б) |
Рис. 4. Схематичне зображення колінеарних (а) і трикутних (б) точок лібрації
114
§ 2. Пружні зіткнення частинок
Зіткнення двох частинок називається пружним, якщо воно не супроводжується зміною їх внутрішнього стану. Найпростіше зіткнення розглядаються в ц-системі. Згідно з r1′=μrGm1 ,
rG2′ = −μrGm2 швидкості частинок до зіткнення в ц-системі пов'яза-
ні з їх швидкостями v10 |
і v20 |
у л-системі співвідношеннями |
||||||||
G′ |
|
|
μ |
G |
G′ |
|
μ |
G |
, |
|
= m |
v0 |
= − m |
v0 |
|||||||
v10 |
, v20 |
|||||||||
|
|
1 |
|
|
2 |
|
|
де vG0 = vG10 −vG20 . У силу закону збереження імпульсу імпульси
обох частинок залишаються після зіткнення рівними за величиною і протилежними за напрямком, а в силу закону збереження енергії залишаються незмінними й їх абсолютні величини.
Таким чином, результат зіткнення зводиться в ц-системі до повороту швидкостей обох частинок, які залишаються взаємно протилежнимиG і незмінними за величиною. Якщо позначити через n одиничний вектор у напрямку швидкості частинки з
масою m1 після зіткнення, то швидкості обох частинок після зіткнення будуть
vG1′ = |
μ |
vnG, |
vG2′ = − |
μ |
vnG . |
|
|
||||
|
m1 |
|
m2 |
Для швидкостей цих частинок вл-системі одержуємо такі вирази:
vG1 = vG1′ +VG
vG2 = vG2′ +VG
|
μ |
G |
|
|
m1vG10 + m2vG20 |
|
|
||||
= |
|
|
vn |
+ |
|
|
|
, |
|
||
m |
|
|
m + m |
|
|||||||
1 |
|
|
|
|
|
1 |
2 |
|
|
||
|
|
μ |
|
G |
|
|
m1vG10 + m2vG20 |
|
|||
= − |
|
vn |
+ |
|
|
|
. |
||||
m |
m |
+ m |
|
||||||||
2 |
|
|
|
|
1 |
2 |
|
|
Одержані результати можна інтерпретувати геометрично. При цьому зручно перейти від швидкостей до імпульсів. Помно-
жаючи швидкості v1 і |
v2 |
на |
m1 |
і |
m2 , |
дістанемо: |
|||
pG1 = μvnG + μ( pG10 + pG |
20 ) m2 , |
pG |
2 |
= −μvnG |
+ μ( pG10 + pG |
20 ) |
m1 . Побу- |
||
|
|
|
|
JJG |
G |
JJG |
|
G |
G |
дуємо коло з радіусом OC = μv ( OC =μv , |
AO = μ( p10 |
+ p20 ) m2 , |
|||||||
|
|
|
115 |
|
|
|
|
|
JJJG |
+ pG |
20 ) |
m1 ) (рис. 5). Якщо одиничний вектор n |
|
|||||
OB = μ( pG10 |
на- |
||||||||
|
|
|
JJG |
JJG |
JJG |
|
|
|
p1 |
правлений уздовж OC , то вектори AC і CB дають імпульси |
|||||||||
і pG2 , відповідно. За заданими |
p10 і pG20 радіус кола і положення |
||||||||
точок A і B залишаються незмінними, |
а точка C може мати |
||||||||
довільне положення на колі. |
Розглянемо |
|
тепер детальні- |
||||||
|
|
|
С |
|
|||||
|
|
|
ше випадок, коли одна із час- |
||||||
|
pG1 |
|
pG |
||||||
|
G |
тинок (нехай це буде частинка |
|||||||
|
|
2 |
з масою |
m ) до зіткнення пе- |
|||||
|
|
n |
|
||||||
|
|
|
|
|
2 |
|
|
|
|
A |
O |
B |
ребувала у стані спокою. Тоді |
||||||
|
|
|
m2 |
|
|
||||
|
|
|
|
довжина |
OB = |
|
p10 = μv |
||
|
|
|
|
|
m1 + m2 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
збігається з радіусом, тобто |
|||||
Рис. 5. Геометрична ілюстрація |
точка B лежить на колі, а век- |
||||||||
процесу пружного зіткнення |
JJJG |
JJJG |
|
JJJG |
|
|
|||
тор AB = AO +OB збігається з |
|||||||||
частинок у ц-системі |
імпульсом p10 |
першої частинки |
|||||||
|
|
|
|
||||||
до розсіювання. При цьому точка A лежить усередині круга, |
|||||||||
якщо m1 < m2 , і – зовні, якщо m1 > m2 (рис. 6). Кути θ1 і θ2 |
яв- |
ляють собою кути відхилення частинки після зіткнення щодо напрямку удару (напрямок p10 ). Центральний кут χ являє собою кут повороту першої частинки в системі центра інерції.
|
|
|
С |
|
|
G |
С |
|
pG |
|
|
|
|
|
|
|
|
G |
|
|
p1 |
G |
|
|
1 |
|
p2 |
|
|
|
p2 |
|
|
|
|
|
|
|
|
|
θ1 |
χ |
θ2 |
θ1max |
θ1 |
χ |
θ2 |
A |
O |
|
B |
|
A |
O |
B |
m1 < m2 |
m1 > m2 |
Рис. 6. Діаграми, що ілюструють процес пружного зіткнення частинок у ц-системі, одна з яких, що має масу m2 ,
перебувала до зіткнення в стані спокою
116
Із рис. 6 видно, що кути θ1 і θ2 можна виразити через кут χ за такими формулами:
tg θ1 |
= |
|
m2 sin χ |
, θ2 = |
π − χ |
. |
|
m1 |
+ m2 cos χ |
|
2 |
||||
|
|
|
|
|
|||
За m1 < m2 швидкість частинки з масою m1 |
після зіткнення має |
довільний напрямок. Якщо ж m1 > m2 , то, як видно із рис. 6, кут відхилення цієї частинки не перевищує кут θ1max , який обчислюється за формулою
θ |
= arcsin |
m2 |
. |
|
1max |
|
m |
|
|
|
|
|
1 |
|
За m1 = m2 маємо θ1 = χ2, θ2 = (π−χ)2 , тобто частинки розлітаються під прямим кутом одна щодо іншої.
§ 3. Розсіювання частинок
Процеси розсіювання двох частинок зображені на рис. 7. Найкоротша відстань між асимптотами АВ і СD називається
прицільною відстанню. Цю відстань можна обчислити таким чином: b = rlim→∞{r sin (rGm,vG)} .
|
C |
|
C |
|
B |
χ |
B |
χ |
|
rG1 (t) |
||||
|
|
r2 (t) |
||
|
rG2 (t) |
|
rG1 (t) |
|
|
b D |
|
b D |
|
|
A |
|
A |
|
|
а) |
|
б) |
Рис. 7. Схематичне зображення процесів розсіювання двох частинок: а) взаємне притягування; б) взаємне відштовхування
117
Будемо розглядати еквівалентну задачу про відхилення однієї частинки з ефективною масою μ у полі U (r) нерухомого сило-
вого центра, розташованого в центрі інерції частинок. Траєкторія частинки в центральному полі симетрична щодо прямої, проведеної в найближчу до центра точку орбіти ( OA на рис. 8). Тому обидві асимптоти орбіти перетинають вказану пряму під однаковими кутами. Якщо позначити ці кути як ϕ0 , то кут відхилення
частинки під час її прольоту повз центр, як видно із рис. 8, дорівнюватиме χ = π− 2ϕ0 . Кут же ϕ0 визначається інтегралом
|
|
L |
dρ |
|
|
|
|
∞ |
ρ2 |
|
|
|
|||
|
|
|
|
||||
ϕ0 = ∫ρmin |
|
|
|
|
. |
||
|
|
|
2 |
||||
|
2μ E −U (ρ) − |
L |
|
||||
ρ2 |
|||||||
|
|
|
|
У випадку інфінітного руху зручно ввести замість E і L інші величини – швидкість v∞ частинки на нескінченності й прицільну відстань b. Енергія і момент імпульсу виражаються через ці величини відпо-
відно E = μv2 |
2 і L =μbv |
, акут ϕ |
0 |
тодідаєтьсяінтегралом |
||||||||||
∞ |
|
|
∞ |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
b |
dρ |
|
|
|
|
||
|
ϕ0 |
= |
∞ |
|
|
ρ2 |
|
|
. |
|||||
|
|
|
|
|
|
|||||||||
|
∫ρmin |
|
b2 |
|
|
|
||||||||
|
|
|
− |
− |
2U |
|||||||||
|
|
|
|
1 |
|
|
|
|
|
|
||||
|
|
|
|
ρ2 |
μv2 |
|||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
χA
b
О ϕ0
Рис. 8. Траєкторія частинки при розсіюванні під дією відштовхувальної сили
118
Нехай пучок однакових частинок з однаковою швидкістю v∞ падає на центр розсіювання, як показано на рис. 9. Різні
частинки в пучку мають різні прицільні відстані, і тому розсіюються під різними кутами. Позначимо через dN кількість частинок, що розсіюються в одиницю часу на кути, які лежать в інтервалі від χ до χ+ dχ. Ця кількість не може бути зруч-
ною для характеристики процесу розсіювання, оскільки вона залежить від густини пучка. Тому зручно ввести відношення
dσ = dNn , де n – кількість частинок, які проходять в одиницю
часу через одиницю площі поперечного перерізу пучка. Це відношення має розмірність площі і називається ефективним перерізом розсіювання. Воно повністю визначається виглядом поля розсіювання.
dχ
A A′
|
b b + db |
dσ |
O |
|
Рис. 9. Розсіювання пучка елементарних частинок під дією центральної сили
Як видно із рис. 9, зв'язок між χ і b є взаємно однознач-
ний; це так, оскільки кут розсіювання є монотонною функцією прицільної відстані. У такому випадку розсіюються в інтервал кутів між χ і χ+ dχ лише ті частинки, які летять з
прицільними відстанями між b(χ) і b(χ) + db(χ) . Кількість
119
таких частинок дорівнює добутку n на площу кільця між колами з радіусами b і b +db , тобто dN = 2πbdbn . Тому ефективний переріз дорівнює dσ = 2πbdb . Щоб знайти залежність ефективного перерізу від кута розсіювання, достатньо переписати цей вираз у вигляді
dσ = 2πb(χ) |
|
db(χ) |
|
dχ . |
|
|
|||
|
dχ |
|
Тут записуємо абсолютне значення похідної dbdχ , оскільки вона може набувати і від'ємного значення. Якщо віднести dσ до елементу тілесного кута dο = 2πsin χdχ , то дістанемо
dσ = b(χ) db dο. sin χ dχ
Для визначення ефективного перерізу залежно від кута розсіювання θ в л-системі необхідно виразити кут χ через
кут θ згідно із формулами tg θ1 = m2 sin χ(m1 + m2 cos χ)
і θ2 = (π−χ)2 .
§ 4. Формула Резерфорда
Одним із найважливіших застосувань отриманих вище формул є розсіювання заряджених частинок у кулонівському полі
U (ρ) = αρ .
Інтеграл для знаходження кута ϕ0 у цьому конкретному випадку записується в такому вигляді:
|
|
|
|
b |
dρ |
|
|
|
|
|
∞ |
|
|
|
|
|
|
, |
|||
|
|
ρ2 |
|
|
|
|||||
ϕ0 = ∫ρ |
|
|
|
|
|
|
|
|
||
min |
|
b2 |
|
|
|
|
||||
|
− |
− |
|
2α |
||||||
|
1 |
|
|
|
|
|
|
|||
|
ρ2 |
|
μv2 |
ρ |
||||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
∞ |
|
|
|
|
120 |
|
|
|
|
|
|