Вектор Лапласа – Рунге – Ленца

редактировать
Вектор, используемый в астрономии

В классической механике, Вектор Лапласа – Рунге – Ленца (LRL) - это вектор , используемый в основном для описания формы и ориентации орбиты одного астрономического тела вокруг другого, например планеты. вращается вокруг звезды. Для двух тел, взаимодействующих с помощью ньютоновской гравитации, вектор LRL является константой движения, что означает, что он одинаков независимо от того, где он рассчитывается на орбите; эквивалентно, вектор LRL называется сохраняющимся. В более общем смысле, вектор LRL сохраняется во всех задачах, в которых два тела взаимодействуют посредством центральной силы, которая изменяется как обратный квадрат расстояния между ними; такие проблемы называются проблемами Кеплера.

атом водорода является проблемой Кеплера, поскольку он состоит из двух заряженных частиц, взаимодействующих по закону Кулона из электростатики, другой обратный квадрат центральная сила. Вектор LRL был существенным в первом квантовомеханическом выводе спектра атома водорода до разработки уравнения Шредингера. Однако сегодня такой подход используется редко.

В классической и квантовой механике сохраняющиеся величины обычно соответствуют симметрии системы. Сохранение вектора LRL соответствует необычной симметрии; задача Кеплера математически эквивалентна частице, свободно движущейся по поверхности четырехмерной (гипер-) сферы, так что вся проблема симметрична при определенных поворотах четырехмерного пространства. Эта более высокая симметрия является результатом двух свойств задачи Кеплера: вектор скорости всегда движется по идеальной окружности и для данной общей энергии все такие скоростные окружности пересекаются друг с другом в те же две точки.

Вектор Лапласа – Рунге – Ленца назван в честь Пьера-Симона де Лапласа, Карла Рунге и Вильгельма Ленца. Он также известен как вектор Лапласа, вектор Рунге – Ленца и вектор Ленца . По иронии судьбы, никто из этих ученых этого не обнаружил. Вектор LRL был повторно открыт несколько раз и также эквивалентен безразмерному вектору эксцентриситета из небесной механики. Были определены различные обобщения вектора LRL, которые включают эффекты специальной теории относительности, электромагнитных полей и даже различных типов центральных сил.

Содержание

  • 1 Контекст
  • 2 История повторных открытий
  • 3 Математическое определение
  • 4 Выведение орбит Кеплера
  • 5 Годографы кругового импульса
  • 6 Константы движения и суперинтегрируемость
  • 7 Эволюция при возмущенных потенциалах
  • 8 Скобки Пуассона
    • 8.1 Немасштабированные функции
    • 8.2 Масштабированные функции
    • 8.3 Инварианты Казимира и уровни энергии
  • 9 Квантовая механика атома водорода
  • 10 Сохранение и симметрия
  • 11 Вращательная симметрия в четырех измерениях
  • 12 Обобщения на другие потенциалы и относительность
  • 13 Доказательства того, что вектор Лапласа – Рунге – Ленца сохраняется в задачах Кеплера
    • 13.1 Прямое доказательство сохранения
    • 13.2 Уравнение Гамильтона – Якоби в параболических координатах
    • 13.3 Теорема Нётер
    • 13.4 Преобразование Ли n
  • 14 Альтернативные вычисления, символы и формулировки
  • 15 См. также
  • 16 Ссылки
  • 17 Дополнительная литература

Контекст

Отдельная частица, движущаяся под любым консервативным центральная сила имеет по крайней мере четыре константы движения, общую энергию E и три декартовых компонента углового импульс вектор Lотносительно начала координат. Орбита частицы ограничена плоскостью, определяемой начальным импульсом pчастицы (или, что эквивалентно, ее скоростью v) и вектором r между частицей и центром. силы (см. рисунок 1 ниже).

Как определено ниже (см. Математическое определение), вектор Лапласа – Рунге – Ленца (вектор LRL) A всегда лежит в плоскости движения для любого центральная сила. Однако A постоянно только для центральной силы, обратной квадрату. Однако для большинства центральных сил этот вектор Aне постоянен, но изменяется как по длине, так и по направлению; если центральная сила приблизительно равна закону обратных квадратов, вектор A приблизительно постоянен по длине, но медленно вращает свое направление. Обобщенный сохраняющийся вектор LRL A {\ displaystyle {\ mathcal {A}}}{\mathcal {A}}можно определить для всех центральных сил, но этот обобщенный вектор является сложной функцией положения и обычно не выражается в замкнутой форме.

Плоскость движения перпендикулярна вектору углового момента L, который постоянен; это может быть выражено математически векторным скалярным произведением уравнением r⋅ L= 0; аналогично, поскольку A лежит в этой плоскости, A⋅ L= 0.

Вектор LRL отличается от других сохраняемых величин следующим свойством. В то время как для типичных сохраняющихся величин существует соответствующая циклическая координата в трехмерном лагранжиане системы, такой координаты для вектора LRL не существует. Таким образом, сохранение вектора LRL должно быть получено напрямую, например, методом скобок Пуассона, как описано ниже. Такие сохраняющиеся величины называются «динамическими», в отличие от обычных «геометрических» законов сохранения, например, углового момента.

История повторного открытия

Вектор LRL A является константой движения важной проблемы Кеплера и полезен при описании астрономические орбиты, такие как движение планет. Тем не менее, он никогда не был хорошо известен среди физиков, возможно, потому, что он менее интуитивен, чем импульс и угловой момент. Следовательно, за последние три столетия он был независимо открыт повторно несколько раз.

Якоб Херманн был первым, кто показал, что A сохраняется для частного случая обратного квадрата центральную силу, и выяснил ее связь с эксцентриситетом орбитального эллипса . Работа Германа была обобщена до ее современной формы Иоганном Бернулли в 1710 году. В конце века Пьер-Симон де Лаплас заново открыл сохранение A, выводя его аналитически, а не геометрически. В середине девятнадцатого века Уильям Роуэн Гамильтон вывел эквивалентный вектор эксцентриситета, определенный ниже, используя его, чтобы показать, что вектор импульса p движется по окружности для движения под действием обратной квадрата центральной силы (рис. 3).

В начале двадцатого века Джозайя Уиллард Гиббс вывел тот же вектор с помощью векторного анализа. Вывод Гиббса был использован в качестве примера Карлом Рунге в популярном немецком учебнике векторов, на который ссылался Вильгельм Ленц в своей статье о (стар.) квантово-механическая обработка атома водорода . В 1926 году вектор был использован Вольфгангом Паули для получения спектра водорода водорода с использованием современной квантовой механики, но не Уравнение Шредингера ; после публикации Паули он стал известен в основном как вектор Рунге – Ленца.

Математическое определение

Для одиночной частицы, на которую действует обратный квадрат центральная сила, описываемая уравнением

F (r) = - kr 2 r ^, {\ displaystyle \ mathbf {F} (r) = {\ frac {-k} {r ^ {2}}} \ mathbf {\ hat {r}},}{\ displaystyle \ mathbf {F} (r) = {\ frac {-k} {r ^ {2}}} \ mathbf {\ hat {r}},}

LRL вектор A математически определяется формулой

Рисунок 1: Вектор LRL A (показан красным) в четырех точках (обозначенных 1, 2, 3 и 4) на эллиптическая орбита связанной точечной частицы, движущейся под действием центральной силы, обратной квадрату . Центр притяжения показан в виде небольшого черного кружка, из которого исходят векторы положения (также черные). Вектор углового момента L перпендикулярен орбите. Копланарные векторы p× Lи (mk / r) r показаны синим и зеленым цветом соответственно; эти переменные определены ниже. Вектор A постоянен по направлению и величине.

A = p × L - mkr ^, {\ displaystyle \ mathbf {A} = \ mathbf {p} \ times \ mathbf {L} - mk \ mathbf {\ hat {r}},}{\ displaystyle \ mathbf {A } = \ mathbf {p} \ times \ mathbf {L} -mk \ mathbf {\ hat {r}},}

где

. сила консервативна, общая энергия E представляет собой константу движения,

E = p 2 2 m - kr = 1 2 mv 2 - kr. {\ displaystyle E = {\ frac {p ^ {2}} {2m}} - {\ frac {k} {r}} = {\ frac {1} {2}} mv ^ {2} - {\ frac {k} {r}}.}{\displaystyle E={\frac {p^{2}}{2m}}-{\frac {k}{r}}={\frac {1}{2}}mv^{2}-{\frac {k}{r}}.}

Кроме того, предполагаемая сила является центральной силой, и, таким образом, вектор углового момента L также сохраняется и определяет плоскость, в которой частица путешествует. Вектор LRL A перпендикулярен вектору углового момента L, потому что оба p× Lи r перпендикулярны L . Отсюда следует, что A лежит в плоскости орбиты .

. Это определение вектора LRL A относится к единственной точечной частице массы м движется под действием фиксированной силы. Однако то же определение можно расширить до задач двух тел, таких как проблема Кеплера, приняв m как приведенную массу двух тел и r как вектор между двумя телами.

Также могут использоваться различные альтернативные формулировки для одной и той же постоянной движения. Наиболее распространенным является масштабирование с помощью mk для определения вектора эксцентриситета (который является вектором с тем же направлением, что и большая полуось, и модуль которого равен эксцентриситету коники):

e = A mk = 1 mk (p × L) - r ^. {\ displaystyle \ mathbf {e} = {\ frac {\ mathbf {A}} {mk}} = {\ frac {1} {mk}} (\ mathbf {p} \ times \ mathbf {L}) - \ mathbf {\ hat {r}}.}{\ displaystyle \ mathbf {e} = {\ frac {\ mathbf {A}} {mk}} = {\ frac {1} {mk}} (\ mathbf {p} \ times \ mathbf {L}) - \ mathbf {\ hat {r}}.}

Получение орбит Кеплера

Рисунок 2: Упрощенная версия рисунка 1, определяющая угол θ между A и r в одной точке орбиты.

Форма и ориентация орбит задачи Кеплера могут быть определены из вектора LRL следующим образом. Взяв скалярное произведение из A с вектором положения r, получаем уравнение

A ⋅ r = A r cos ⁡ θ = r ⋅ (p × L) - мкр, {\ displaystyle \ mathbf {A} \ cdot \ mathbf {r} = Ar \ cos \ theta = \ mathbf {r} \ cdot \ left (\ mathbf {p} \ times \ mathbf {L} \ справа) -mkr,}{\displaystyle \mathbf {A} \cdot \mathbf {r} =Ar\cos \theta =\mathbf {r} \cdot \left(\mathbf {p} \times \mathbf {L} \right)-mkr,}

где θ - это угол между r и A (рисунок 2). Переставляя тройное скалярное произведение

r ⋅ (p × L) = (r × p) ⋅ L = L ⋅ L = L 2 {\ displaystyle \ mathbf {r} \ cdot \ left (\ mathbf {p } \ times \ mathbf {L} \ right) = \ left (\ mathbf {r} \ times \ mathbf {p} \ right) \ cdot \ mathbf {L} = \ mathbf {L} \ cdot \ mathbf {L} = L ^ {2}}\mathbf{r} \cdot\left(\mathbf{p}\times \mathbf{L}\right) = \left(\mathbf{r} \times \mathbf{p}\right)\cdot\mathbf{L} = \mathbf{L}\cdot\mathbf{L}=L^2

и перестановка дает определяющую формулу для конического сечения , при условии, что A является константой, что имеет место для закона силы обратных квадратов,

1 r знак равно mk L 2 (1 + A mk cos ⁡ θ) {\ displaystyle {\ frac {1} {r}} = {\ frac {mk} {L ^ {2}}} \ left (1 + {\ frac { A} {mk}} \ cos \ theta \ right)}{\displaystyle {\frac {1}{r}}={\frac {mk}{L^{2}}}\left(1+{\frac {A}{mk}}\cos \theta \right)}

из эксцентриситета e,

e = A mk = | А | mk {\ displaystyle e = {\ frac {A} {mk}} = {\ frac {\ left | \ mathbf {A} \ right |} {mk}}}e = \frac{A}{mk} = \frac{\left|\mathbf{A}\right|}{m k}

и прямая кишка

| 2 ℓ | = 2 L 2 м k. {\ displaystyle \ left | 2 \ ell \ right | = {\ frac {2L ^ {2}} {mk}}.}{\ displaystyle \ left | 2 \ ell \ right | = {\ frac {2L ^ {2}} {mk}}.}

Большая полуось a конического сечения может быть определена с помощью прямой кишки и эксцентриситета

a (1 ± e 2) = ℓ = L 2 mk, {\ displaystyle a \ left (1 \ pm e ^ {2} \ right) = \ ell = {\ frac {L ^ {2}} { mk}},}{\ displaystyle a \ left (1 \ pm e ^ {2} \ right) = \ ell = {\ frac {L ^ {2}} {mk}},}

где знак минус относится к эллипсам, а знак плюс к гиперболам.

Взяв скалярное произведение A на себя, получаем уравнение с энергией E,

A 2 = m 2 k 2 + 2 m EL 2, {\ displaystyle A ^ {2} = m ^ {2} k ^ {2} + 2mEL ^ {2},}{\ displaystyle A ^ {2} = m ^ {2} k ^ {2} + 2mEL ^ {2},}

, которое можно переписать в терминах эксцентриситета,

e 2 - 1 = 2 L 2 mk 2 E. {\ displaystyle e ^ {2} -1 = {\ frac {2L ^ {2}} {mk ^ {2}}} E.}{\ displaystyle e ^ {2} -1 = {\ frac {2L ^ {2} } {mk ^ {2}}} E.}

Таким образом, если энергия E отрицательна (связанные орбиты), эксцентриситет меньше единицы, а орбита представляет собой эллипс . И наоборот, если энергия положительна (несвязанные орбиты, также называемые «рассеянными орбитами»), эксцентриситет больше единицы, а орбита представляет собой гиперболу. Наконец, если энергия точно равна нулю, эксцентриситет равен единице, а орбита представляет собой параболу. Во всех случаях направление A лежит вдоль оси симметрии конического сечения и указывает от центра силы к перицентру , точке наибольшего сближения.

Годографы кругового импульса

Рис. 3. Вектор импульса p (показан синим) движется по окружности, как частица движется по эллипсу. Четыре отмеченные точки соответствуют точкам на рисунке 1. Круг с центром на оси y в позиции A / L (показан пурпурным) с радиусом mk / L (показан зеленым). Угол η определяет эксцентриситет e эллиптической орбиты (cos η = e). По теореме о вписанном угле для окружностей, η также является углом между любой точкой окружности и двумя точками пересечения с осью p x, p x = ± p 0, которые зависят только от E, но не L.

Сохранение вектора LRL A и вектора углового момента L полезен для демонстрации того, что вектор импульса p движется по окружности под действием центральной силы, обратной квадрату.

Взяв скалярное произведение

mkr ^ = p × L - A {\ displaystyle mk {\ hat {\ mathbf {r}}} = \ mathbf {p} \ times \ mathbf {L} - \ mathbf {A}}{\ displaystyle mk {\ hat {\ mathbf {r}}} = \ mathbf {p} \ times \ mathbf {L} - \ mathbf {A}}

с собой дает

(mk) 2 = A 2 + p 2 L 2 + 2 L ⋅ (p × A). {\ displaystyle (mk) ^ {2} = A ^ {2} + p ^ {2} L ^ {2} +2 \ mathbf {L} \ cdot (\ mathbf {p} \ times \ mathbf {A}).}{\ displaystyle ( mk) ^ {2} = A ^ {2} + p ^ {2} L ^ {2} +2 \ mathbf {L} \ cdot (\ mathbf {p} \ times \ mathbf {A}).}

Далее выбирая L вдоль оси z и большую полуось в качестве оси x, получаем уравнение геометрического места для p,

px 2 + (py - AL) 2 = (mk L) 2. {\ displaystyle p_ {x} ^ {2} + \ left (p_ {y} - {\ frac {A} {L}} \ right) ^ {2} = \ left ({\ frac {mk} {L}) } \ right) ^ {2}.}{\ displaystyle p_ {x} ^ {2} + \ left (p_ { y} - {\ frac {A} {L}} \ right) ^ {2} = \ left ({\ frac {mk} {L}} \ right) ^ {2}.}

Другими словами, вектор импульса p ограничен кругом радиуса mk / L = L / ℓ с центром в (0, A / L). Эксцентриситет e соответствует косинусу угла η, показанного на рисунке 3.

В вырожденном пределе круговых орбит и, таким образом, исчезающем A, окружность центрируется в начале координат (0, 0). Для краткости также полезно ввести переменную p 0 = 2 m | E | {\ displaystyle p_ {0} = {\ sqrt {2m | E |}}}{\displaystyle p_{0}={\sqrt {2m|E|}}}.

Этот круговой годограф полезен для иллюстрации симметрии проблемы Кеплера.

Константы движения и суперинтегрируемость

Семь скалярных величин E, A и L (будучи векторами, последние две вносят по три сохраняющиеся величины каждая) связаны двумя уравнениями, A⋅ L= 0 и A = mk + 2 mEL, что дает пять независимых констант движения. (Поскольку величина A, следовательно, эксцентриситет e орбиты, может быть определена из полного углового момента L и энергии E, только направление A сохраняется независимо; кроме того, поскольку A должен быть перпендикулярен L, он вносит только одну дополнительную сохраняемую величину.)

Это согласуется с шестью начальными условиями (начальное положение частицы и векторы скорости , каждый с тремя компонентами), которые определяют орбиту частицы, поскольку начальное время не определяется постоянной движения. Таким образом, полученная 1-мерная орбита в 6-мерном фазовом пространстве полностью определена.

Механическая система с d степенями свободы может иметь не более 2d - 1 констант движения, так как есть 2d начальных условий, и начальное время не может быть определено константой движения. Система с более чем d константами движения называется суперинтегрируемой, а система с 2d - 1 константами называется максимально суперинтегрируемой. Поскольку решение уравнения Гамильтона – Якоби в одной системе координат может дать только d констант движения, суперинтегрируемые системы должны быть разделимы более чем в одной системе координат. Задача Кеплера максимально суперинтегрируема, поскольку имеет три степени свободы (d = 3) и пять независимых констант движения ; его уравнение Гамильтона – Якоби разделимо как в сферических координатах, так и в параболических координатах, как описано ниже.

Максимально суперинтегрируемые системы движутся по замкнутым одномерным орбитам в фазовое пространство, поскольку орбита является пересечением фазового пространства изоповерхностей их постоянных движения. Следовательно, орбиты перпендикулярны всем градиентам всех этих независимых изоповерхностей, пяти в этой конкретной задаче, и, следовательно, определяются обобщенными перекрестными произведениями всех этих градиентов. В результате все суперинтегрируемые системы автоматически описываются механикой Намбу, альтернативно и эквивалентно гамильтоновой механике.

Максимально суперинтегрируемые системы могут быть квантованы с использованием коммутационных соотношений, как показано ниже. Тем не менее, эквивалентным образом, они также квантованы в рамках Намбу, например, в этой классической проблеме Кеплера в квантовый атом водорода.

Эволюция при возмущенных потенциалах

Рис. 5: Постепенно прецессирующая эллиптическая орбита с эксцентриситетом e = 0,667. Такая прецессия возникает в задаче Кеплера, если притягивающая центральная сила немного отклоняется от закона обратных квадратов. Скорость прецессии может быть рассчитана с использованием формул в тексте.

Вектор Лапласа – Рунге – Ленца A сохраняется только для идеального обратного квадрата центрального сила. Однако в большинстве практических задач, таких как движение планет, взаимодействие потенциальной энергии между двумя телами не является в точности законом обратных квадратов, но может включать в себя дополнительную центральную силу, так называемую возмущение, описываемое потенциальной энергией h (r). В таких случаях вектор LRL медленно вращается в плоскости орбиты, что соответствует медленной апсидальной прецессии орбиты.

По предположению, возмущающий потенциал h (r) является консервативной центральной силой, что подразумевает, что полная энергия E и угловой момент вектор L сохранены. Таким образом, движение все еще находится в плоскости, перпендикулярной L, и величина A сохраняется из уравнения A = mk + 2mEL. Потенциал возмущения h (r) может быть функцией любого вида, но должен быть значительно слабее, чем основная сила обратного квадрата между двумя телами.

Скорость, с которой вектор LRL вращается, предоставляет информацию о возмущающем потенциале h (r). Используя каноническую теорию возмущений и координаты действие-угол, легко показать, что A вращается со скоростью

∂ ∂ L ⟨h (r)⟩ = ∂ ∂ L {1 T ∫ 0 T час (r) dt} знак равно ∂ ∂ L {м L 2 ∫ 0 2 π r 2 час (r) d θ}, {\ displaystyle {\ begin {align} {\ frac {\ partial } {\ partial L}} \ langle h (r) \ rangle = \ displaystyle {\ frac {\ partial} {\ partial L}} \ left \ {{\ frac {1} {T}} \ int _ { 0} ^ {T} час (г) \, dt \ right \} \\ [1em] = \ displaystyle {\ frac {\ partial} {\ partial L}} \ left \ {{\ frac {m} { L ^ {2}}} \ int _ {0} ^ {2 \ pi} r ^ {2} h (r) \, d \ theta \ right \}, \ end {align}}}{\ displaystyle {\ begin {align} { \ frac {\ partial} {\ partial L}} \ langle h (r) \ rangle = \ displaystyle {\ frac {\ partial} {\ partial L}} \ left \ {{\ frac {1} {T} } \ int _ {0} ^ {T} h (r) \, dt \ right \} \\ [1em] = \ displaystyle {\ frac {\ partial} {\ partial L}} \ left \ {{\ гидроразрыв {m} {L ^ {2}}} \ int _ {0} ^ {2 \ pi} r ^ {2} h (r) \, d \ theta \ right \}, \ end {align}}}

где T - период обращения, и тождество L dt = mr dθ использовалось для преобразования интеграла по времени в угловой интеграл (рисунок 5). Выражение в угловых скобках, ⟨h (r)⟩, представляет возмущающий потенциал, но усредненный за один полный период; то есть усредненное значение за один полный проход тела по орбите. Математически это среднее по времени соответствует следующей величине в фигурных скобках. Это усреднение помогает подавить колебания скорости вращения.

Этот подход использовался для проверки теории общей теории относительности Эйнштейна, которая добавляет небольшое эффективное обратнокубическое возмущение к нормальному ньютоновскому гравитационному потенциалу,

h (r) = k L 2 м 2 c 2 (1 r 3). {\ displaystyle h (r) = {\ frac {kL ^ {2}} {m ^ {2} c ^ {2}}} \ left ({\ frac {1} {r ^ {3}}} \ right).}{\ displaystyle h (r) = {\ frac {kL ^ {2}} {m ^ {2} c ^ {2}}} \ left ({\ frac {1} {r ^ {3}}} \ right).}

Подставляя эту функцию в интеграл и используя уравнение

1 r = mk L 2 (1 + A mk cos ⁡ θ) {\ displaystyle {\ frac {1} {r}} = {\ frac {mk} {L ^ {2}}} \ left (1 + {\ frac {A} {mk}} \ cos \ theta \ right)}{\displaystyle {\frac {1}{r}}={\frac {mk}{L^{2}}}\left(1+{\frac {A}{mk}}\cos \theta \right)}

, чтобы выразить r через θ, скорость прецессии перицентра , вызванная этим неньютоновским возмущением, вычисляется как

6 π k 2 TL 2 c 2, {\ displaystyle {\ frac {6 \ pi k ^ { 2}} {TL ^ {2} c ^ {2}}},}{\ displaystyle {\ frac {6 \ pi k ^ {2}} {TL ^ {2} c ^ {2}}},}

, что близко соответствует наблюдаемой аномальной прецессии Меркурия и бинарных пульсаров. Это согласие с экспериментом является убедительным доказательством общей теории относительности.

скобок Пуассона

Немасштабированных функций

Алгебраическая структура проблемы, как объясняется в последующих разделах, такова: SO (4) / 2 ~ SO (3) × SO (3). Три компонента L i вектора углового момента L имеют скобки Пуассона

{L i, L j} = ∑ s = 1 3 ε ijs L s, {\ displaystyle \ {L_ {i}, L_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} L_ {s},}{\ displaystyle \ {L_ {i}, L_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} L_ {s},}

где i = 1, 2, 3 и ϵ ijs - полностью антисимметричный тензор, то есть символ Леви-Чивита ; здесь используется индекс суммирования s, чтобы избежать путаницы с параметром силы k, определенным выше. Тогда, поскольку вектор LRL A преобразуется как вектор, мы имеем следующие отношения скобок Пуассона между A и L:

{A i, L j} = ∑ s = 1 3 ε ijs A s. {\ displaystyle \ {A_ {i}, L_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} A_ {s}.}{\ displaystyle \ {A_ {i}, L_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} A_ {s }.}

Наконец, скобки Пуассона между различными компонентами A следующие:

{A i, A j} = - 2 m ∑ s = 1 3 ε ijs L s H, {\ displaystyle \ {A_ {i}, A_ {j} \} = - 2m \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} L_ {s} H,}{\ displaystyle \ {A_ {i}, A_ {j} \} = - 2m \ sum _ {s = 1} ^ { 3} \ varepsilon _ {ijs} L_ {s} H,}

где H {\ displaystyle H}H- гамильтониан. Обратите внимание, что диапазон компонентов A и компонентов L не закрыт скобками Пуассона из-за множителя H {\ displaystyle H}Hв правой части этого последнего соотношения.

Наконец, поскольку и L, и A являются константами движения, мы имеем

{A i, H} = {L i, H} = 0. {\ displaystyle \ {A_ {i}, H \} = \ {L_ {i}, H \} = 0.}{\ displaystyle \ {A_ {i}, H \} = \ {L_ {i}, H \} = 0.}

Скобки Пуассона будут расширены до квантово-механического коммутационные отношения в следующем разделе и скобки в следующем разделе.

Масштабированные функции

Как отмечено ниже масштабированный вектор Лапласа – Рунге – Ленца D может быть определен в тех же единицах, что и угловой момент, путем деления A на p 0 = 2 м | H | {\ displaystyle p_ {0} = {\ sqrt {2m | H |}}}{\ displaystyle p_ {0} = {\ sqrt {2m | H |} }} . Поскольку D по-прежнему преобразуется как вектор, скобки Пуассона для D с вектором углового момента L могут быть записаны в аналогичном форма

{D i, L j} = ∑ s = 1 3 ε ijs D s. {\ displaystyle \ {D_ {i}, L_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} D_ {s}.}{\displaystyle \{D_{i}, L_{j}\}=\sum _{s=1}^{3}\varepsilon _{ijs}D_{s}.}

скобки Пуассона из D с самим собой зависят от знака H, то есть от того, является ли энергия отрицательной (создавая замкнутые эллиптические орбиты при обратном квадратная центральная сила) или положительная (создавая открытые гиперболические орбиты под действием центральной силы, обратной квадрату). Для отрицательных энергий, то есть для связанных систем, скобки Пуассона равны

{D i, D j} = ∑ s = 1 3 ε i j s L s. {\ displaystyle \ {D_ {i}, D_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} L_ {s}.}{\ displaystyle \ {D_ {i}, D_ {j} \} = \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} L_ {s}.}

Теперь мы можем оценить мотивацию для выбранного масштабирования D : При таком масштабировании гамильтониан больше не появляется в правой части предыдущего соотношения. Таким образом, диапазон трех компонентов L и трех компонентов D образует шестимерную алгебру Ли под скобкой Пуассона. Эта алгебра Ли изоморфна so (4), алгебре Ли 4-мерной группы вращений SO (4).

Напротив, для положительной энергии скобки Пуассона имеют противоположный знак,

{D i, D j} = - ∑ s = 1 3 ε ijs L s. {\ displaystyle \ {D_ {i}, D_ {j} \} = - \ sum _ {s = 1} ^ {3} \ varepsilon _ {ijs} L_ {s}.}{\displaystyle \{D_{i},D_{j}\}=-\sum _{s=1}^{3}\varepsilon _{ijs}L_{s}.}

В этом случае Алгебра Ли изоморфна so (3,1).

Различие между положительной и отрицательной энергиями возникает из-за того, что желаемое масштабирование - то, которое исключает гамильтониан из правой части скобок Пуассона между компонентами масштабированного вектора LRL - включает квадратный корень из гамильтониан. Чтобы получить функции с действительными значениями, мы должны затем взять абсолютное значение гамильтониана, который различает положительные значения (где | H | = H {\ displaystyle | H | = H}{\displaystyle |H|=H}) и отрицательные значения (где | H | = - H {\ displaystyle | H | = -H}{\ displaystyle | H | = -H} ).

Инварианты Казимира и уровни энергии

Инварианты Казимира для отрицательных энергий равны

C 1 = D ⋅ D + L ⋅ L = m k 2 2 | E |, {\ displaystyle C_ {1} = \ mathbf {D} \ cdot \ mathbf {D} + \ mathbf {L} \ cdot \ mathbf {L} = {\ frac {mk ^ {2}} {2 | E | }},}{\ displaystyle C_ {1} = \ mathbf {D} \ cdot \ mathbf { D} + \ mathbf {L} \ cdot \ mathbf {L} = {\ frac {mk ^ {2}} {2 | E |}},}
C 2 = D ⋅ L = 0, {\ displaystyle C_ {2} = \ mathbf {D} \ cdot \ mathbf {L} = 0,}{\ displaystyle C_ {2} = \ mathbf {D} \ cdot \ ma thbf {L} = 0,}

и имеют исчезающие скобки Пуассона со всеми компоненты D и L,

{C 1, L i} = {C 1, D i} = {C 2, L i} = {C 2, D i} = 0. {\ displaystyle \ {C_ {1}, L_ {i} \} = \ {C_ {1}, D_ {i} \} = \ {C_ {2}, L_ {i} \} = \ {C_ {2}, D_ {i} \} = 0.}{\ displaystyle \ {C_ {1}, L_ {i} \} = \ {C_ {1}, D_ {i} \} = \ {C_ {2}, L_ {i} \} = \ {C_ {2}, D_ {i} \} = 0.}

C2тривиально равно нулю, поскольку два вектора всегда перпендикулярны.

Однако другой инвариант, C 1, нетривиален и зависит только от m, k и E.При каноническом квантовании этот инвариант допускает уровни энергии водорода -подобные атомы должны быть получены с использованием только квантово-механических канонических коммутационных соотношений вместо обычного решения уравнения Шредингера. Этот вывод подробно обсуждается в следующем разделе.

Квантовая механика атома водорода

Рис. 6. Уровни энергии атома водорода, предсказанные из коммутационных соотношений углового момента и векторных операторов Лапласа – Рунге – Ленца; эти уровни энергии были проверены экспериментально.

Скобки Пуассона предоставляют простое руководство для квантования большинства классических систем : соотношение коммутации двух квантово-механических операторы задаются скобкой Пуассона соответствующих классических переменных, умноженной на iħ.

Путем выполнения этого квантования и вычисления собственных значений C 1 Оператор Казимира для задачи Кеплера, Вольфганг Паули смог получить уровни энергии водородоподобных атомов (рис. 6) и, таким образом, их спектр атомного излучения. Этот элегантный вывод 1926 года был получен до разработки уравнения Шредингера.

Тонкость квантово-механического оператора для вектора LRL A заключается в том, что операторы импульса и момента количества движения не коммутируют; следовательно, квантовый оператор перекрестное произведение из p и L должен быть определен тщательно. Обычно операторы для декартовых компонентов Asопределяются с использованием симметризованного (эрмитова) произведения,

A s = - mkr ^ s + 1 2 ∑ i = 1 3 ∑ j = 1 3 ε sij ( пи ℓ j + ℓ jpi), {\ displaystyle A_ {s} = - mk {\ hat {r}} _ {s} + {\ frac {1} {2}} \ sum _ {i = 1} ^ { 3} \ sum _ {j = 1} ^ {3} \ varepsilon _ {sij} (p_ {i} \ ell _ {j} + \ ell _ {j} p_ {i}),}{\ displaystyle A_ {s} = - mk {\ hat {r}} _ {s} + {\ frac {1} {2}} \ sum _ {i = 1} ^ {3} \ sum _ {j = 1} ^ {3} \ varepsilon _ {sij} (p_ {i} \ ell _ {j} + \ ell _ {j} p_ {i}),}

После этого сделано, можно показать, что квантовые операторы LRL удовлетворяют соотношениям коммутации, точно аналогичным отношениям скобок Пуассона из предыдущего раздела - просто заменив скобку Пуассона на 1 / (i ℏ) {\ displaystyle 1 / (i \ hbar)}{\ displaystyle 1 / (я \ hbar)} , умноженное на коммутатор.

Из этих операторов можно определить дополнительные лестничные операторы для L,

J 0 = A 3, {\ displaystyle J_ {0} = A_ {3},}{\ displaystyle J_ {0} = A_ {3},}
J ± 1 = ∓ 1 2 (A 1 ± i A 2). {\ displaystyle J _ {\ pm 1} = \ mp {\ frac {1} {\ sqrt {2}}} \ left (A_ {1} \ pm iA_ {2} \ right).}{ \ displaystyle J _ {\ pm 1} = \ mp {\ frac {1} {\ sqrt {2}}} \ left (A_ {1} \ pm iA_ {2} \ right).}

Эти дополнительные соединения разные собственные состояния L, поэтому разные спиновые мультиплеты между собой.

Нормализованный первый инвариантный оператор Казимира, квантовый аналог описанного выше, может быть определен аналогичным образом,

C 1 = - mk 2 2 ℏ 2 H - 1 - I, {\ displaystyle C_ {1} = - {\ frac {mk ^ {2}} {2 \ hbar ^ {2}}} H ^ {- 1} -I,}{\displaystyle C_{1}=-{\frac {mk^{2}}{2\hbar ^{2}}}H^{-1}-I,}

где H является обратным к гамильтонову оператору энергии, а I - тождественный оператор .

. Применяя эти лестничные операторы к собственным состояниям | ℓmn〉 полного углового момента, азимутального углового момента и операторов энергии, собственных значений первого оператора Казимира, C 1, как видно, квантованы, n - 1. Важно отметить, что из-за исчезновения C 2 они не зависят от квантовых чисел ℓ и m, что делает уровни энергии вырождаются.

Следовательно, уровни энергии задаются следующим образом:

E n = - mk 2 2 ℏ 2 n 2, {\ displaystyle E_ {n} = - {\ frac {mk ^ {2}} {2 \ hbar ^ {2} n ^ {2}}},}{\ displaystyle E_ {n} = - {\ frac {mk ^ {2}} {2 \ hbar ^ {2} n ^ {2}}},}

, что совпадает с формулой Ридберга для водородоподобных атомов (рисунок 6). Дополнительные операторы симметрии A соединили разные ℓ мультиплеты между собой для заданной энергии (и C 1), определяя n состояний на каждом уровне. Фактически они увеличили группу углового момента SO (3) до SO (4) /ℤ2~ SO (3) × SO (3).

Сохранение и симметрия

Сохранение вектора LRL соответствует тонкой симметрии системы. В классической механике симметрии - это непрерывные операции, которые отображают одну орбиту на другую без изменения энергии системы; в квантовой механике симметрии - это непрерывные операции, которые «смешивают» электронные орбитали с одинаковой энергией, то есть вырожденные уровни энергии. С такими симметриями обычно связывают сохраняющуюся величину. Например, каждая центральная сила симметрична относительно группы вращения SO (3), что приводит к сохранению углового момента L. Обычно полное вращение системы не влияет на энергию орбиты; квантово-механически, вращения смешивают сферические гармоники одного и того же квантового числа l без изменения энергии.

Рис. 7. Семейство годографов кругового импульса для заданной энергии E. Все круги проходят через одни и те же две точки ± p 0 = ± 2 м | E | {\ displaystyle \ pm p_ {0} = \ pm {\ sqrt {2m | E |}}}{\ displaystyle \ pm p_ {0} = \ pm {\ sqrt {2m | E |}}} на оси p x (см. рисунок 3). Это семейство годографов соответствует одному семейству аполлоновских кругов и σ-изоповерхностям биполярных координат.

. Симметрия для центральной силы обратного квадрата выше и тоньше. Своеобразная симметрия задачи Кеплера приводит к сохранению как вектора углового момента L, так и вектора LRL A (как определено выше) и квантово-механически, гарантирует, что уровни энергии водорода не зависят от квантовых чисел углового момента l и m. Однако симметрия более тонкая, потому что операция симметрии должна иметь место в многомерном пространстве ; такие симметрии часто называют «скрытыми симметриями».

Классически более высокая симметрия проблемы Кеплера допускает непрерывные изменения орбит, которые сохраняют энергию, но не угловой момент; Другими словами, орбиты с одинаковой энергией, но с различным угловым моментом (эксцентриситетом) могут непрерывно преобразовываться друг в друга. Квантово-механически это соответствует смешиванию орбиталей, которые различаются квантовыми числами l и m , например s (l = 0) и p (l = 1) атомными орбиталями. Такое смешивание не может быть выполнено с помощью обычных трехмерных перемещений или вращений, но эквивалентно вращению в более высоком измерении.

Для отрицательных энергий, т.е. для связанных систем, группа высшей симметрии SO (4), которая сохраняет длину четырехмерных векторов

| е | 2 = е 1 2 + е 2 2 + е 3 2 + е 4 2. {\ displaystyle | \ mathbf {e} | ^ {2} = e_ {1} ^ {2} + e_ {2} ^ {2} + e_ {3} ^ {2} + e_ {4} ^ {2}.}{\displaystyle |\mathbf {e} |^{2}=e_{1}^{2}+e_{2}^{2}+e_{3}^{2}+e_{4}^{2}.}

В 1935 году Владимир Фок показал, что квантово-механическая ограниченная проблема Кеплера эквивалентна проблеме свободной частицы, ограниченной трехмерной единичной сферой . в четырехмерном пространстве. В частности, Фок показал, что волновая функция Шредингера в импульсном пространстве для задачи Кеплера была стереографической проекцией сферических гармоник на сфера. Вращение сферы и повторное проецирование приводят к непрерывному отображению эллиптических орбит без изменения энергии; квантово-механически это соответствует смешиванию всех орбиталей с одинаковым квантовым числом энергии n. Валентин Баргманн впоследствии заметил, что скобки Пуассона для вектора углового момента L и масштабированный вектор LRL D образуют алгебру Ли для SO (4). Проще говоря, шесть величин D и L соответствуют шести сохраняющимся угловым моментам в четырех измерениях, связанных с шестью возможными простыми вращениями в этом пространстве (там шесть способов выбора двух осей из четырех). Этот вывод не означает, что наша вселенная представляет собой трехмерную сферу; это просто означает, что эта конкретная физическая проблема (задача двух тел для обратных квадратов центральных сил ) математически эквивалентна свободной частице на трехмерной сфере.

Для положительных энергий - т. Е. Для несвязанных, «рассеянных» систем - группа высшей симметрии SO (3,1), которая сохраняет длину Минковского 4-вектора

ds 2 = e 1 2 + e 2 2 + e 3 2 - e 4 2. {\ displaystyle ds ^ {2} = e_ {1} ^ {2} + e_ {2} ^ {2} + e_ {3} ^ {2} -e_ {4} ^ {2}.}{\ displaystyle ds ^ {2} = e_ {1} ^ {2} + e_ {2} ^ {2} + e_ {3} ^ {2} - e_ {4} ^ {2}.}

Оба случаи отрицательной и положительной энергии были рассмотрены Фоком и Баргманном и энциклопедически рассмотрены Бандером и Ициксоном.

Орбиты систем с центральной силой и орбиты проблемы Кеплера в частности - также симметричны под отражением. Следовательно, указанные выше группы SO (3), SO (4) и SO (3,1) не являются полными группами симметрии своих орбит; полные группы: O (3), O (4) и O (3,1), соответственно. Тем не менее, только связанные подгруппы, SO (3), SO (4) и SO (3,1) необходимы для демонстрации сохранения углового момента и векторы LRL; симметрия отражения не имеет отношения к сохранению, что может быть получено из алгебры Ли группы.

Вращательная симметрия в четырех измерениях

Рисунок 8: Годографы импульса на рисунке 7 соответствуют стереографическим проекциям больших кругов на трехмерной единичной сфере η. Все большие круги пересекают ось η x, которая перпендикулярна странице; проекция идет от северного полюса (единичный вектор w ) на плоскость η x−ηy, как показано здесь для пурпурного годографа пунктирными черными линиями. Большой круг на широте α соответствует эксцентриситету e = sin α. Цвета больших кругов, показанные здесь, соответствуют соответствующим им годографам на рисунке 7.

Связь между проблемой Кеплера и четырехмерной симметрией вращения SO (4) может быть легко визуализируется. Обозначим четырехмерные декартовы координаты (w, x, y, z), где (x, y, z) представляют собой декартовы координаты вектора нормального положения r. Трехмерный вектор импульса p связан с четырехмерным вектором η {\ displaystyle {\ boldsymbol {\ eta}}}{\ boldsymbol {\ eta}} на трехмерной единичной сфере.

η знак равно п 2 - p 0 2 p 2 + p 0 2 w ^ + 2 p 0 p 2 + p 0 2 p = mk - rp 0 2 mkw ^ + rp 0 mkp, {\ displaystyle {\ begin { выровнено} {\ boldsymbol {\ eta}} = \ displaystyle {\ frac {p ^ {2} -p_ {0} ^ {2}} {p ^ {2} + p_ {0} ^ {2}}} \ mathbf {\ hat {w}} + {\ frac {2p_ {0}} {p ^ {2} + p_ {0} ^ {2}}} \ mathbf {p} \\ [1em] = \ displaystyle {\ frac {mk-rp_ {0} ^ {2}} {mk}} \ mathbf {\ hat {w}} + {\ frac {rp_ {0}} {mk}} \ mathbf {p}, \ end {align}}}{\ displaystyle {\ begin {align} {\ boldsymbol {\ eta}} = \ displaystyle {\ frac {p ^ {2} -p_ {0} ^ {2}} {p ^ {2} + p_ {0} ^ {2}}} \ mathbf {\ hat {w}} + {\ frac {2p_ {0}} {p ^ {2} + p_ {0} ^ {2}}} \ mathbf {p} \\ [1em] = \ displaystyle {\ frac {mk -rp_ {0} ^ {2}} {mk}} \ mathbf {\ hat {w}} + {\ frac {rp_ {0}} {mk}} \ mathbf {p}, \ end {выравнивается}}}

где w ^ {\ displaystyle \ mathbf {\ hat {w}}}\ mathbf {\ hat {w}} - единичный вектор вдоль новой оси w. Преобразование p в η может быть однозначно инвертировано; например, компонент импульса x равен

px = p 0 η x 1 - η w, {\ displaystyle p_ {x} = p_ {0} {\ frac {\ eta _ {x}} {1- \ eta _ {w}}},}{\ displaystyle p_ {x} = p_ {0} {\ frac {\ eta _ {x}} {1- \ eta _ {w}}}, }

и аналогично для p y и p z. Другими словами, трехмерный вектор p является стереографической проекцией четырехмерного η {\ displaystyle {\ boldsymbol {\ eta}}}{\ boldsymbol {\ eta}} вектор, масштабированный на p 0 (рисунок 8).

Без потери общности, мы можем исключить нормальную вращательную симметрию, выбрав декартовы координаты так, чтобы ось z была выровнена с вектором углового момента L и Годографы импульса выровнены, как на рисунке 7, с центрами кружков на оси y. Поскольку движение является плоским, а p и L перпендикулярны, p z = η z = 0 и внимание может быть ограничено трехмерный вектор η {\ displaystyle {\ boldsymbol {\ eta}}}{\ boldsymbol {\ eta}} = (η w, η x, η y). Семейство аполлонических кругов годографов импульса (рис. 7) соответствует семейству больших кругов на трехмерном η {\ displaystyle {\ boldsymbol {\ eta} }}{\ boldsymbol {\ eta}} сфера, каждая из которых пересекает ось η x в двух фокусах η x = ± 1, соответствующих фокусам годографа импульса при p <391.>x = ± p 0. Эти большие круги связаны простым вращением вокруг оси η x (рисунок 8). Эта вращательная симметрия преобразует все орбиты с одинаковой энергией друг в друга; однако такой поворот ортогонален обычным трехмерным поворотам, поскольку он преобразует четвертое измерение η w. Эта более высокая симметрия характерна для задачи Кеплера и соответствует сохранению вектора LRL.

Элегантное переменные действие-угол решение проблемы Кеплера может быть получено путем исключения избыточных четырехмерных координат η {\ displaystyle {\ boldsymbol {\ eta}}}{\ boldsymbol {\ eta}} в пользу эллиптических цилиндрических координат (χ, ψ, φ)

η w = cn ⁡ χ cn ⁡ ψ, {\ displaystyle \ eta _ {w} = \ operatorname {cn} \ chi \ имя оператора {cn} \ psi,}{\ displaystyle \ eta _ {w} = \ operatorname {cn} \ chi \ operatorname {cn} \ psi,}
η x = sn χ dn ⁡ ψ cos ⁡ ϕ, {\ displaystyle \ eta _ {x} = \ operatorname {sn} \, \ chi \ operatorname {dn} \ psi \, \ соз \ фи,}{\ displaystyle \ eta _ {x} = \ OperatorName {sn} \, \ chi \ operatorname {dn} \ psi \, \ cos \ phi,}
η y = sn ⁡ χ dn ⁡ ψ sin ⁡ ϕ, {\ displaystyle \ eta _ {y} = \ operatorname {sn} \ chi \ operatorname {dn} \ psi \ sin \ фи,}{\displaystyle \eta _{y}=\operatorname {sn} \chi \operatorname {dn} \psi \sin \phi,}
η z = dn ⁡ χ sn ⁡ ψ, {\ displaystyle \ eta _ {z} = \ operatorname {dn} \ chi \ operatorname {sn} \ psi,}{\ отображает tyle \ eta _ {z} = \ operatorname {dn} \ chi \ operatorname {sn} \ psi,}

где sn, cn и dn - это эллиптические функции Якоби..

Обобщения для других потенциалов и теории относительности.

Вектор Лапласа – Рунге – Ленца также можно обобщить для определения сохраняющихся величин, применимых к другим ситуациям.

В присутствии однородного электрического поля E обобщенный вектор Лапласа – Рунге – Ленца A {\ displaystyle {\ mathcal {A}}}{\mathcal {A}}is

A = A + mq 2 [(r × E) × r], {\ displaystyle {\ mathcal {A}} = \ mathbf {A} + {\ frac {mq} {2}} \ left [\ left (\ mathbf { r} \ times \ mathbf {E} \ right) \ times \ mathbf {r} \ right],}{\displaystyle {\mathcal {A}}=\mathbf {A} +{\frac {mq}{2}}\left[\left(\mathbf {r} \times \mathbf {E} \right)\times \mathbf {r} \right],}

где q - заряд движущейся по орбите частицы. Хотя A {\ displaystyle {\ mathcal {A}}}{\mathcal {A}}не сохраняется, оно приводит к сохранению величины, а именно A ⋅ E {\ displaystyle {\ mathcal {A} } \ cdot \ mathbf {E}}{\mathcal {A}}\cdot {\mathbf {E}}.

Далее обобщая вектор Лапласа – Рунге – Ленца на другие потенциалы и специальную теорию относительности, наиболее общую форму можно записать как

A = (∂ ξ ∂ u) (п ​​× L) + [ξ - u (∂ ξ ∂ u)] L 2 r ^, {\ displaystyle {\ mathcal {A}} = \ left ({\ frac {\ partial \ xi} {\ частичный u}} \ right) \ left (\ mathbf {p} \ times \ mathbf {L} \ right) + \ left [\ xi -u \ left ({\ frac {\ partial \ xi} {\ partial u} } \ right) \ right] L ^ {2} \ mathbf {\ hat {r}},}{\ displaystyle {\ mathcal {A}} = \ left ({\ frac {\ partial \ xi} {\ partial u}} \ right) \ left (\ mathbf {p} \ times \ mathbf {L} \ right) + \ left [\ xi -u \ left ({\ frac {\ partial \ xi} {\ partial u}} \ right) \ right] L ^ {2} \ mathbf {\ hat {r}},}

где u = 1 / r (см. теорему Бертрана ) и ξ = cos θ, с угол θ, определяемый формулой

θ = L ∫ udum 2 c 2 (γ 2-1) - L 2 u 2, {\ displaystyle \ theta = L \ int ^ {u} {\ frac {du} {\ sqrt {m ^ {2} c ^ {2} (\ gamma ^ {2} -1) -L ^ {2} u ^ {2}}}},}{\displaystyle \theta =L\int ^{u}{\frac {du}{\sqrt {m^{2}c^{2}(\gamma ^{2}-1)-L^{2}u^{2}}}},}

и γ - фактор Лоренца. Как и раньше, мы можем получить сохраняющийся бинормальный вектор B, взяв векторное произведение с сохраняющимся вектором углового момента

B = L × A. {\ displaystyle {\ mathcal {B}} = \ mathbf {L} \ times {\ mathcal {A}}.}{\ displaystyle {\ mathcal {B}} = \ mathbf {L} \ times {\ mathcal {A}}.}

Эти два вектора можно аналогичным образом объединить в сохраняющийся диадический тензор W,

W = α A ⊗ A + β B ⊗ B. {\ displaystyle {\ mathcal {W}} = \ alpha {\ mathcal {A}} \ otimes {\ mathcal {A}} + \ beta \, {\ mathcal {B}} \ otimes {\ mathcal {B}}.}{\ displaystyle {\ mathcal {W}} = \ alpha {\ mathcal {A}} \ otimes {\ mathcal {A}} + \ beta \, {\ mathcal {B}} \ otimes {\ mathcal {B}}.}

На иллюстрации можно вычислить вектор LRL для нерелятивистского изотропного гармонического осциллятора. Поскольку сила центральная,

F (r) = - kr, {\ displaystyle \ mathbf {F} (r) = - k \ mathbf {r},}{\displaystyle \mathbf {F} (r)=-k\mathbf {r},}

, угловой момент вектор сохраняется, и движение лежит в плоскости.

Сохраняющийся диадический тензор можно записать в простой форме

W = 1 2 mp ⊗ p + k 2 r ⊗ r, {\ displaystyle {\ mathcal {W}} = {\ frac {1 } {2m}} \ mathbf {p} \ otimes \ mathbf {p} + {\ frac {k} {2}} \, \ mathbf {r} \ otimes \ mathbf {r},}{\ displaystyle {\ mathcal {W}} = {\ frac {1} {2m}} \ mathbf {p} \ otimes \ mathbf { p} + {\ frac {k} {2}} \, \ mathbf {r} \ otimes \ mathbf {r},}

хотя p и r не обязательно перпендикулярны.

Соответствующий вектор Рунге – Ленца более сложен,

A = 1 mr 2 ω 0 A - mr 2 E + L 2 {(p × L) + (mr ω 0 A - mr E) г ^}, {\ displaystyle {\ mathcal {A}} = {\ frac {1} {\ sqrt {mr ^ {2} \ omega _ {0} A-mr ^ {2} E + L ^ {2} }}} \ left \ {\ left (\ mathbf {p} \ times \ mathbf {L} \ right) + \ left (mr \ omega _ {0} A-mrE \ right) \ mathbf {\ hat {r} } \ right \},}{\ displaystyle {\ mathcal {A}} = {\ frac {1} {\ sqrt {mr ^ {2} \ omega _ {0} A-mr ^ {2} E + L ^ {2}}}} \ left \ {\ left (\ mathbf {p} \ times \ mathbf {L} \ right) + \ left (mr \ omega _ {0} A-mrE \ right) \ mathbf {\ hat {r}} \ right \},}

где

ω 0 = km {\ displaystyle \ omega _ {0} = {\ sqrt {\ frac {k} {m}}}}\ omega _ {0} = {\ sqrt {\ frac {k} {m}}}

- собственное колебание частота, и

A = (E 2 - ω 2 L 2) 1/2 / ω. {\ displaystyle A = (E ^ {2} - \ omega ^ {2} L ^ {2}) ^ {1/2} / \ omega.}{\displaystyle A=(E^{2}-\omega ^{2}L^{2})^{1/2}/\omega.}

Доказательства того, что вектор Лапласа – Рунге – Ленца сохраняется в Задачи Кеплера

Следующие аргументы показывают, что вектор LRL сохраняется при центральных силах, которые подчиняются закону обратных квадратов.

Прямое доказательство сохранения

Центральная сила F {\ displaystyle \ mathbf {F}}\ mathbf {F} , действующая на частицу, равна

F = dpdt = е (г) рр знак равно е (г) р ^ {\ Displaystyle \ mathbf {F} = {\ гидроразрыва {d \ mathbf {p}} {dt}} = f (г) {\ гидроразрыва {\ mathbf {г} } {r}} = f (r) \ mathbf {\ hat {r}}}\ mathbf {F} = \ frac {d \ mathbf {p}} {dt} = f (r) \ frac {\ mathbf {r}} {r} = f (r) \ mathbf {\ hat {r}}

для некоторой функции f (r) {\ displaystyle f (r)}f(r)радиуса р {\ displaystyle r}r. Поскольку угловой момент L = r × p {\ displaystyle \ mathbf {L} = \ mathbf {r} \ times \ mathbf {p}}\ mathbf {L} = \ mathbf {r} \ times \ mathbf {p} сохраняется при центральном силы, ddt L = 0 {\ displaystyle {\ frac {d} {dt}} \ mathbf {L} = 0}\ frac {d} {dt} \ mathbf {L} = 0 и

ddt (p × L) = dpdt × L знак равно е (г) р ^ × (г × mdrdt) = е (г) MR [г (г ⋅ drdt) - r 2 drdt], {\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) = {\ frac {d \ mathbf {p}} {dt}} \ times \ mathbf {L} = f (r) \ mathbf {\ hat {r}} \ times \ left (\ mathbf {r} \ times m {\ frac {d \ mathbf {r}} {dt}} \ right) = f (r) {\ frac {m} {r}} \ left [\ mathbf {r} \ left (\ mathbf {r} \ cdot {\ frac {d \ mathbf {r}} {dt}} \ right) -r ^ {2} {\ frac {d \ mathbf {r}} { dt}} \ right],}{\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) = {\ frac {d \ mathbf {p}} {dt}} \ times \ mathbf {L} = f (r) \ mathbf {\ hat {r}} \ times \ left (\ mathbf {r} \ times m {\ frac {d \ mathbf { r}} {dt}} \ right) = f (r) {\ frac {m} {r}} \ left [\ mathbf {r} \ left (\ mathbf {r} \ cdot {\ frac {d \ mathbf {r}} {dt}} \ right) -r ^ {2} {\ frac {d \ mathbf {r}} {dt}} \ right],}

где momentum p = mdrdt {\ displaystyle \ mathbf {p} = m {\ frac {d \ mathbf {r}} {dt} }}\ mathbf {p} = m \ frac {d \ mathbf {r}} {dt} и где тройное векторное произведение было упрощено с использованием формулы Лагранжа

r × (r × drdt) = r (r ⋅ drdt) - r 2 drdt. {\ displaystyle \ mathbf {r} \ times \ left (\ mathbf {r} \ times {\ frac {d \ mathbf {r}} {dt}} \ right) = \ mathbf {r} \ left (\ mathbf { r} \ cdot {\ frac {d \ mathbf {r}} {dt}} \ right) -r ^ {2} {\ frac {d \ mathbf {r}} {dt}}.}{\displaystyle \mathbf {r} \times \left(\mathbf {r} \times {\frac {d\mathbf {r} }{dt}}\right)=\mathbf {r} \left(\mathbf {r} \cdot {\frac {d\mathbf {r} } {dt}}\right)-r^{2}{\frac {d\mathbf {r} }{dt}}.}

Тождество

ddt (r ⋅ r) = 2 r ⋅ drdt = ddt (r 2) = 2 rdrdt {\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {r} \ cdot \ mathbf {r } \ right) = 2 \ mathbf {r} \ cdot {\ frac {d \ mathbf {r}} {dt}} = {\ frac {d} {dt}} (r ^ {2}) = 2r {\ frac {dr} {dt}}}{\displaystyle {\frac {d}{dt}}\left(\mathbf {r} \cdot \mathbf {r} \right)=2\mathbf {r} \cdot {\frac {d\mathbf {r} }{dt}}={\frac {d}{dt}}(r^{2})=2r{\frac {dr}{dt}}}

дает уравнение

ddt (p × L) = - mf (r) r 2 [1 rdrdt - rr 2 drdt] = - mf (r) r 2 ddt ( р-р). {\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) = - mf (r) r ^ {2} \ left [{\ frac {1 } {r}} {\ frac {d \ mathbf {r}} {dt}} - {\ frac {\ mathbf {r}} {r ^ {2}}} {\ frac {dr} {dt}} \ right] = - mf (r) r ^ {2} {\ frac {d} {dt}} \ left ({\ frac {\ mathbf {r}} {r}} \ right).}{\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) = - mf (r) r ^ {2} \ left [{\ frac {1} {r}} {\ frac {d \ mathbf {r}} {dt}} - {\ frac {\ mathbf {r}} {r ^ {2}}} {\ frac {dr} {dt}} \ right] = - mf (r) r ^ {2} {\ frac {d} {dt}} \ left ({\ frac {\ mathbf {r}} {r}) } \ right).}

Для частный случай центральной силы, обратной квадрату f (r) = - kr 2 {\ displaystyle f (r) = {\ frac {-k} {r ^ {2}}}}f (r) = \ frac {-k} {r ^ {2}} , это равно

ddt (p × L) = mkddt (rr) = ddt (mkr ^). {\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) = mk {\ frac {d} {dt}} \ left ({\ frac { \ mathbf {r}} {r}} \ right) = {\ frac {d} {dt}} (mk \ mathbf {\ hat {r}}).}{\ displaystyle {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) = mk {\ frac {d} {dt}} \ left ({\ frac {\ mathbf {r}} {r}} \ right) = {\ frac {d} {dt}} (mk \ mathbf {\ hat {r}}).}

Следовательно, A сохраняется для центральных сил, обратных квадрату

ddt A = ddt (p × L) - ddt (mkr ^) = 0. {\ displaystyle {\ frac {d} {dt}} \ mathbf {A} = {\ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) - {\ frac {d} {dt}} \ left (mk \ mathbf {\ hat {r}} \ right) = \ mathbf {0}.}{\ displaystyle {\ frac {d} {dt}} \ mathbf {A} = { \ frac {d} {dt}} \ left (\ mathbf {p} \ times \ mathbf {L} \ right) - {\ frac {d} {dt}} \ left (mk \ mathbf {\ hat {r}) } \ right) = \ mathbf {0}.}

Более короткое доказательство получается при использовании отношения углового момента к угловой скорости L = mr 2 ω {\ displaystyle \ mathbf {L} = mr ^ {2} {\ boldsymbol {\ omega}}}\ mathbf {L} = mr ^ 2 \ boldsymbol {\ omega} , что справедливо для частицы, движущейся в плоскости, перпендикулярной L {\ displaystyle \ mathbf {L}}\mathbf{L}. Если задать обратные квадраты центральных сил, производная по времени от p × L {\ displaystyle \ mathbf {p} \ times \ mathbf {L}}\mathbf{p} \times \mathbf{L}равна

ddtp × L = (- kr 2 r ^) × (mr 2 ω) знак равно mk ω × r ^ = mkddtr ^, {\ displaystyle {\ frac {d} {dt}} \ mathbf {p} \ times \ mathbf {L} = \ left ( {\ frac {-k} {r ^ {2}}} \ mathbf {\ hat {r}} \ right) \ times \ left (mr ^ {2} {\ boldsymbol {\ omega}} \ right) = mk \, {\ boldsymbol {\ omega}} \ times \ mathbf {\ hat {r}} = mk \, {\ frac {d} {dt}} \ mathbf {\ hat {r}},}{\ displaystyle {\ frac {d} {dt}} \ mathbf {p} \ times \ mathbf {L} = \ left ({\ frac {-k} {r ^ {2}}} \ mathbf {\ hat {r}} \ right) \ times \ left (mr ^ {2} {\ b oldsymbol {\ omega}} \ right) = mk \, {\ boldsymbol {\ omega}} \ times \ mathbf {\ hat {r}} = mk \, {\ frac {d} {dt}} \ mathbf {\ шляпа {r}},}

где последнее равенство выполняется, потому что единичный вектор может изменяться только вращением, и ω × r ^ {\ displaystyle {\ boldsymbol {\ omega}} \ times \ mathbf {\ hat {r}}}\boldsymbol{\omega}\times\mathbf{\hat{r}}- орбитальная скорость вращающегося вектора. Таким образом, A рассматривается как разность двух векторов с равными производными по времени.

Как описано в другом месте этой статьи, этот вектор LRL A является частным случаем общего сохраняемого вектора A {\ displaystyle {\ mathcal {A }}}{\mathcal {A}}, который может быть определен для всех центральных сил. Однако, поскольку большинство центральных сил не создают замкнутых орбит (см. теорему Бертрана ), аналогичный вектор A {\ displaystyle {\ mathcal {A}}}{\mathcal {A}}редко имеет простое определение и обычно является многозначной функцией угла θ между r и A {\ displaystyle {\ mathcal {A}}}{\mathcal {A}}.

уравнением Гамильтона – Якоби в параболические координаты

Постоянство вектора LRL также может быть получено из уравнения Гамильтона – Якоби в параболических координатах (ξ, η), которые определяются уравнения

ξ = r + x, {\ displaystyle \ xi = r + x,}{\ displaystyle \ xi = r + x,}
η = r - x, {\ displaystyle \ eta = rx,}{\displaystyle \ eta =r-x,}

где r представляет радиус в плоскости орбиты

r = x 2 + y 2. {\ displaystyle r = {\ sqrt {x ^ {2} + y ^ {2}}}.}{\ displaystyle r = {\ sqrt {x ^ { 2} + y ^ {2}}}.}

Инверсия этих координат:

x = 1 2 (ξ - η), {\ displaystyle x = {\ frac {1} {2}} (\ xi - \ eta),}{\ displaystyle x = {\ frac {1} {2}} (\ xi - \ eta),}
y = ξ η, {\ displaystyle y = {\ sqrt {\ xi \ eta}},}{\ displaystyle y = {\ sqrt {\ xi \ eta}},}

Разделение уравнение Гамильтона – Якоби в этих координатах дает два эквивалентных уравнения

2 ξ p ξ 2 - mk - m E ξ = - Γ, {\ displaystyle 2 \ xi p _ {\ xi} ^ {2} -mk-mE \ xi = - \ Gamma,}{\ displaystyle 2 \ xi p _ {\ xi} ^ {2} -mk-mE \ xi = - \ Gamma,}
2 η p η 2 - mk - m E η = Γ, {\ displaystyle 2 \ eta p _ {\ eta} ^ {2} -mk-mE \ eta = \ Gamma, }{\displaystyle 2\eta p_{\eta }^{2}-mk-mE\eta =\Gamma,}

где Γ - постоянная движения. Вычитание и повторное выражение с точки зрения декартовых импульсов p x и p y показывает, что Γ эквивалентно вектору LRL

Γ = py (xpy - ypx) - mkxr = А х. {\ displaystyle \ Gamma = p_ {y} (xp_ {y} -yp_ {x}) - mk {\ frac {x} {r}} = A_ {x}.}{\ displaystyle \ Gamma = p_ {y} (xp_ {y} -yp_ {x}) - mk {\ frac {x} {r}} = A_ {x}.}

Теорема Нётер

Связь между вращательной симметрией, описанной выше, и сохранением вектора LRL может быть сделана количественной посредством теоремы Нётер. Эта теорема, которая используется для нахождения постоянных движения, утверждает, что любое бесконечно малое изменение обобщенных координат физической системы

δ qi = ε gi (q, q ˙, t) {\ displaystyle \ delta q_ {i} = \ varepsilon g_ {i} (\ mathbf {q}, \ mathbf {\ dot {q}}, t)}{\ displaystyle \ delta q_ {i} = \ varepsilon g_ { i} (\ mathbf {q}, \ mathbf {\ dot {q}}, t)}

, который заставляет лагранжиан изменяться первым упорядочить по полной производной по времени

δ L = ε ddt G (q, t) {\ displaystyle \ delta L = \ varepsilon {\ frac {d} {dt}} G (\ mathbf {q}, t)}{\displaystyle \delta L=\varepsilon {\frac {d}{dt}}G(\mathbf {q},t)}

соответствует сохраняющейся величине Γ

Γ = - G + ∑ igi (∂ L ∂ q ˙ i). {\ displaystyle \ Gamma = -G + \ sum _ {i} g_ {i} \ left ({\ frac {\ partial L} {\ partial {\ dot {q}} _ {i}}} \ right).}{\displaystyle \Gamma =-G+\sum _{i}g_{i}\left({\frac {\partial L}{\partial {\dot {q}}_{i}}}\right).}

В частности, сохраняющаяся компонента вектора LRL A ​​s соответствует изменению координат

δ xi = ε 2 [2 пикселя - xips - δ равно (r ⋅ p)], { \ displaystyle \ delta x_ {i} = {\ frac {\ varepsilon} {2}} \ left [2p_ {i} x_ {s} -x_ {i} p_ {s} - \ delta _ {is} \ left ( \ mathbf {r} \ cdot \ mathbf {p} \ right) \ right],}{ \ displaystyle \ delta x_ {i} = {\ frac {\ varepsilon} {2}} \ left [2p_ {i} x_ {s} -x_ {i} p_ {s} - \ delta _ {is} \ left ( \ mathbf {r} \ cdot \ mathbf {p} \ right) \ right ],}

где i равно 1, 2 и 3, с x i и p i являющиеся i-ми компонентами векторов положения и импульса r и p, соответственно; как обычно, δ is представляет собой дельту Кронекера. Результирующее изменение первого порядка в лагранжиане составляет

δ L = 1 2 ε m k d d t (x s r). {\ displaystyle \ delta L = {\ frac {1} {2}} \ varepsilon mk {\ frac {d} {dt}} \ left ({\ frac {x_ {s}} {r}} \ right). }{\displaystyle \delta L={\frac {1}{2}}\varepsilon mk{\frac {d}{dt}}\left({\frac {x_{s}}{r}}\right).}

Подстановка в общую формулу для сохраняющейся величины Γ дает сохраняющуюся компоненту A s вектора LRL,

A s = [p 2 xs - ps (r ⋅ p)] - mk (xsr) = [p × (r × p)] s - mk (xsr). {\ Displaystyle A_ {s} = \ left [p ^ {2} x_ {s} -p_ {s} \ left (\ mathbf {r} \ cdot \ mathbf {p} \ right) \ right] -mk \ left ({\ frac {x_ {s}} {r}} \ right) = \ left [\ mathbf {p} \ times \ left (\ mathbf {r} \ times \ mathbf {p} \ right) \ right] _ {s} -mk \ left ({\ frac {x_ {s}} {r}} \ right).}{\ displaystyle A_ {s} = \ left [p ^ {2} x_ {s} -p_ {s} \ left (\ mathbf {r} \ cdot \ mathbf {p} \ right) \ right] -mk \ left ({\ frac {x_ {s}} {r}} \ right) = \ left [\ mathbf {p} \ times \ left (\ mathbf {r } \ times \ mathbf {p} \ right) \ right] _ {s} -mk \ left ({\ frac {x_ {s}} {r}} \ right).}

Преобразование Ли

Рис. 9: Преобразование Ли, из которого выполняется сохранение вектора LRL A выводится. При изменении масштабного параметра λ энергия и угловой момент изменяются, но эксцентриситет e, величина и направление A - нет.

Теорема Нётер вывод сохранения вектора LRL A элегантно, но имеет один недостаток: изменение координат δx i включает не только положение r, но и импульс p или, что то же самое, скорость v . Этот недостаток может быть устранен путем получения сохранения A с использованием подхода, впервые предложенного Sophus Lie. В частности, можно определить преобразование Ли, в котором координаты r и время t масштабируются с помощью различных степеней параметра λ (рисунок 9),

t → λ 3 t, r → λ 2 r, p → 1 λ p. {\ displaystyle t \ rightarrow \ lambda ^ {3} t, \ qquad \ mathbf {r} \ rightarrow \ lambda ^ {2} \ mathbf {r}, \ qquad \ mathbf {p} \ rightarrow {\ frac {1} {\ lambda}} \ mathbf {p}.}{\displaystyle t\rightarrow \lambda ^{3}t,\qquad \mathbf {r} \rightarrow \lambda ^{2 }\mathbf {r},\qquad \mathbf {p} \rightarrow {\frac {1}{\lambda }}\mathbf {p}.}

Это преобразование изменяет полный угловой момент L и энергию E,

L → λ L, E → 1 λ 2 E, {\ displaystyle L \ rightarrow \ lambda L, \ qquad E \ rightarrow {\ frac {1} {\ lambda ^ {2}}} E,}{\ displaystyle L \ rightarrow \ lambda L, \ qquad E \ rightarrow {\ frac {1} {\ lambda ^ {2}}} E, }

, но сохраняет их произведение EL. Следовательно, эксцентриситет e и величина A сохраняются, как можно видеть из уравнения для A

A 2 = m 2 k 2 e 2 = m 2 k 2 + 2 m E L 2. {\ displaystyle A ^ {2} = m ^ {2} k ^ {2} e ^ {2} = m ^ {2} k ^ {2} + 2mEL ^ {2}.}{\ displaystyle A ^ {2} = m ^ {2} k ^ {2} e ^ {2} = m ^ {2} k ^ {2} + 2mEL ^ {2}.}

Направление A также сохраняется, поскольку полуоси не изменяются глобальным масштабированием. Это преобразование также сохраняет третий закон Кеплера, а именно, что полуось a и период T образуют постоянную T / a.

Альтернативные вычисления, символы и формулировки

В отличие от импульса и углового момента векторов pи L общепринятого определения вектора Лапласа – Рунге – Ленца не существует; В научной литературе используется несколько различных коэффициентов масштабирования и символов. Наиболее распространенное определение дается выше, но другой распространенной альтернативой является деление на константу mk для получения безразмерного сохраненного вектора эксцентриситета

e = 1 mk (p × L) - r ^ знак равно mk (v × (r × v)) - r ^, {\ displaystyle \ mathbf {e} = {\ frac {1} {mk}} \ left (\ mathbf {p} \ times \ mathbf {L} \ справа) - \ mathbf {\ hat {r}} = {\ frac {m} {k}} \ left (\ mathbf {v} \ times \ left (\ mathbf {r} \ times \ mathbf {v} \ right) \ right) - \ mathbf {\ hat {r}},}{\ displaystyle \mathbf {e} ={\frac {1}{mk}}\left(\mathbf {p} \times \mathbf {L} \right)-\mathbf {\hat {r}} ={\frac {m }{k}}\left(\mathbf {v} \times \left(\mathbf {r} \times \mathbf {v} \right)\right)-\mathbf {\hat {r}},}

где v - вектор скорости. Этот масштабированный вектор e имеет то же направление, что и A, и его величина равна эксцентриситету орбиты и, таким образом, исчезает для круговых орбит.

Возможны и другие масштабированные версии, например, разделив A только на m

M = v × L - kr ^, {\ displaystyle \ mathbf {M} = \ mathbf {v} \ times \ mathbf {L} -k \ mathbf {\ hat {r}},}{\ displaystyle \ mathbf {M} = \ mathbf {v} \ times \ mathbf {L} -k \ mathbf {\ hat {r}},}

или на p 0

D = A p 0 = 1 2 m | E | {p × L - mkr ^}, {\ displaystyle \ mathbf {D} = {\ frac {\ mathbf {A}} {p_ {0}}} = {\ frac {1} {\ sqrt {2m | E | }}} \ left \ {\ mathbf {p} \ times \ mathbf {L} -mk \ mathbf {\ hat {r}} \ right \},}{\ displaystyle \ mathbf {D} = {\ frac {\ mathbf {A }} {p_ {0}}} = {\ frac {1} {\ sqrt {2m | E |}}} \ left \ {\ mathbf {p} \ times \ mathbf {L} -mk \ mathbf {\ hat {r}} \ right \},}

который имеет те же единицы измерения, что и angular импульс вектор L.

В редких случаях знак вектора LRL может быть изменен на противоположный, т. е. масштабироваться на -1. Другие общие символы для вектора LRL включают a, R, F, Jи V . Однако выбор масштаба и символа для вектора LRL не влияет на его сохранение.

Рисунок 4: Вектор углового момента L, вектор LRL A и Гамильтон вектор, бинормаль B, взаимно перпендикулярны; A и B точки вдоль большой и малой осей, соответственно, эллиптической орбиты задачи Кеплера.

Альтернативным сохраняющимся вектором является бинормаль вектор B, изученный Уильямом Роуэном Гамильтоном,

B = p - (mk L 2 r) (L × r), {\ displaystyle \ mathbf {B} = \ mathbf {p} - \ left ({\ frac {mk} {L ^ {2} r}} \ right) \ \ left (\ mathbf {L} \ times \ mathbf {r} \ right),}{\ displaystyle \ mathbf {B} = \ mathbf {p} - \ left ({\ frac {mk } {L ^ {2} r}} \ right) \ \ left (\ mathbf {L} \ times \ mathbf {r} \ right),}

который сохраняется и указывает по малой полуоси эллипса. (Он не определен для исчезающего эксцентриситета.)

Вектор LRL A= B× Lпредставляет собой векторное произведение из B и L (рисунок 4). На годографе импульса в соответствующем разделе выше легко увидеть, что B связывает начало импульсов с центром кругового годографа и имеет величину A / L. В перигелии он указывает в направлении импульса.

Вектор B обозначается как «бинормальный», поскольку он перпендикулярен как A, так и L . Подобно самому вектору LRL, бинормальный вектор может быть определен с различными масштабами и символами.

Два консервативных вектора, A и B, могут быть объединены, чтобы сформировать консервативный диадический тензор W,

W = α A ⊗ A + β B ⊗ В, {\ displaystyle \ mathbf {W} = \ alpha \ mathbf {A} \ otimes \ mathbf {A} + \ beta \, \ mathbf {B} \ otimes \ mathbf {B},}{\ displaystyle \ mathbf {W} = \ alpha \ mathbf {A} \ otimes \ mathbf {A} + \ beta \, \ mathbf {B} \ otimes \ mathbf {B},}

где α и β - произвольные константы масштабирования, а ⊗ {\ displaystyle \ otimes}\otimes представляет тензорное произведение (которое не связано с векторным перекрестным произведением , несмотря на их похожий символ). Написанное в явных компонентах, это уравнение имеет вид

W i j = α A i A j + β B i B j. {\ displaystyle W_ {ij} = \ alpha A_ {i} A_ {j} + \ beta B_ {i} B_ {j}.}{\ displaystyle W_ {ij} = \ alpha A_ {i} A_ {j } + \ beta B_ {i} B_ {j}.}

Будучи перпендикулярными друг другу, векторы A и B можно рассматривать как главные оси сохраняющегося тензора W, то есть его масштабированные собственные векторы. Wперпендикулярны L,

L ⋅ W знак равно α (L ⋅ A) A + β (L ⋅ B) B знак равно 0, {\ Displaystyle \ mathbf {L} \ cdot \ mathbf {W} = \ alpha \ left (\ mathbf {L} \ cdot \ mathbf {A} \ right) \ mathbf {A} + \ beta \ left (\ mathbf {L} \ cdot \ mathbf {B} \ right) \ mathbf {B} = 0,}{\displaystyle \mathbf {L} \cdot \mathbf {W} =\alpha \left(\mathbf {L} \cdot \mathbf {A} \right)\mathbf {A} +\beta \left(\mathbf {L} \cdot \mathbf {B} \right)\mathbf {B} =0,}

поскольку A и B также перпендикулярны L, L⋅ A= L⋅ B= 0.

Более конкретно, это уравнение в явных компонентах читается как

(L ⋅ W) J знак равно α (∑ я знак равно 1 3 L я A я) A J + β (∑ я = 1 3 L я В я) B J знак равно 0. {\ Displaystyle \ left (\ mathbf {L} \ cdot \ mathbf {W} \ right) _ {j} = \ alpha \ left (\ sum _ {i = 1} ^ {3} L_ {i} A_ {i} \ right) A_ {j} + \ beta \ left (\ sum _ {i = 1} ^ {3} L_ {i} B_ {i} \ right) B_ {j} = 0.}{\ displaystyle \ left (\ mathbf {L} \ cdot \ mathbf {W} \ right) _ {j} = \ alpha \ left (\ sum _ {i = 1} ^ {3} L_ {i} A_ { i} \ right) A_ {j} + \ beta \ left (\ sum _ {i = 1} ^ {3} L_ {i} B_ {i} \ right) B_ {j} = 0.}

См. также

Ссылки

Дополнительная литература

Последняя правка сделана 2021-05-26 13:25:38
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте