Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

IvanovMaksyuta

.pdf
Скачиваний:
35
Добавлен:
14.02.2016
Размер:
2.63 Mб
Скачать
і поміщаємо початок системи координат у центр

Лекція 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 rG2

 

m rG2

 

 

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 видно, що радіус-вектори частинки в ц-системі направлені протилежно, а їхні довжини обернено пропорційні масам ( r1r2′ = m2 m1 ), тобто радіус-вектор частин-

ки з більшою масою менший, і навпаки.

Якщо траєкторія руху μ -точки є еліпсом, то і реальні частин-

ки рухаються еліптичними орбітами, причому обидва еліпси мають спільний фокус, розташований y центрі інерції (рис. 2).

Наприклад, для Сонця і Землі m1 m2 3 106 , тобто Сонце відхи-

ляється від центра мас на відстань 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

 

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]