Уравнение Гамильтона – Якоби

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

В физике, уравнение Гамильтона – Якоби, названное в честь Уильяма Роуэна Гамильтона и Карла Густава Джейкоба Якоби, является альтернативной формулировкой классической механики, эквивалентные другим формулировкам, таким как законы движения Ньютона, лагранжева механика и гамильтонова механика. Уравнение Гамильтона – Якоби особенно полезно для определения сохраняющихся величин для механических систем, что может быть возможным, даже если сама механическая проблема не может быть решена полностью.

Уравнение Гамильтона – Якоби также является единственной формулировкой механики, в которой движение частицы может быть представлено как волна. В этом смысле он выполнил давнюю цель теоретической физики (датируемую, по крайней мере, Иоганном Бернулли в восемнадцатом веке) - найти аналогию между распространением света и движением частицы. Волновое уравнение, которому следуют механические системы, похоже, но не идентично уравнению Шредингера, как описано ниже; по этой причине уравнение Гамильтона – Якоби считается «наиболее близким подходом» классической механики к квантовой механике.

В математике уравнение Гамильтона – Якоби является необходимое условие, описывающее экстремальную геометрию в обобщениях задач из вариационного исчисления. Его можно рассматривать как частный случай уравнения Гамильтона – Якоби – Беллмана из динамического программирования.

Содержание
  • 1 Обозначение
  • 2 Математическая формулировка
  • 3 Сравнение с другими формулировки механики
  • 4 Вывод с использованием канонического преобразования
  • 5 Действие и функции Гамильтона
  • 6 Разделение переменных
    • 6.1 Примеры в различных системах координат
      • 6.1.1 Сферические координаты
      • 6.1.2 Эллиптические цилиндрические координаты
      • 6.1.3 Параболические цилиндрические координаты
  • 7 Волны и частицы
    • 7.1 Оптические волновые фронты и траектории
    • 7.2 Связь с уравнением Шредингера
  • 8 Приложения
    • 8.1 HJE в гравитационном поле
    • 8.2 HJE в электромагнитных полях
      • 8.2.1 Циркулярно поляризованная волна
      • 8.2.2 Монохроматическая линейно поляризованная плоская волна
      • 8.2.3 Электромагнитная волна с соленоидальным магнитным полем
  • 9 См. также
  • 10 Ссылки
  • 11 Дополнительная литература
Обозначение

Переменные, выделенные жирным шрифтом, такие как q {\ displaystyle \ mathbf {q}}\ mathbf {q} представляет список N {\ displaystyle N}N обобщенных координат,

q = (q 1, q 2,…, q N - 1, q N) {\ displaystyle \ mathbf {q} = (q_ {1}, q_ {2}, \ ldots, q_ {N-1}, q_ {N})}{\ displaystyle \ mathbf {q} = (q_ { 1}, q_ {2}, \ ldots, q_ {N-1}, q_ {N})}

точка над переменная или список обозначают производную по времени (см. нотация Ньютона ). Например,

q ˙ = d q d t. {\ displaystyle {\ dot {\ mathbf {q}}} = {\ frac {d \ mathbf {q}} {dt}}.}{\ displaystyle {\ dot {\ mathbf {q}}} = {\ frac {d \ mathbf {q}} {dt}}.}

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

p ⋅ q = ∑ k = 1 N pkqk. {\ displaystyle \ mathbf {p} \ cdot \ mathbf {q} = \ sum _ {k = 1} ^ {N} p_ {k} q_ {k}.}{\ displaystyle \ mathbf {p} \ cdot \ mathbf {q} = \ sum _ {k = 1} ^ {N} p_ {k} q_ {k}.}
Математическая формулировка

Учитывая гамильтониан H (q, p, t) {\ displaystyle H (q, p, t)}{\ displaystyle H (q, p, t)} механической системы (где q {\ displaystyle q}q , p {\ displaystyle p}p - координаты и импульсы системы, а t {\ displaystyle t}t - время) записывается уравнение Гамильтона – Якоби как нелинейное уравнение в частных производных для главной функции Гамильтона S (q, t) {\ displaystyle S (q, t)}{\ displaystyle S (q, t)} ,

- ∂ S ∂ t = H (q, ∂ S ∂ q, t). {\ displaystyle - {\ frac {\ partial S} {\ partial t}} = H \ left (q, {\ frac {\ partial S} {\ partial q}}, t \ right).}{\ displaystyle - {\ frac {\ partial S} {\ partial t}} = H \ left (q, {\ frac {\ partial S} {\ partial q}}, t \ right).}

Основная функция Гамильтона определяется как функция верхнего предела интеграла действия , взятого вдоль минимальной траектории действия системы,

S (q, t): = ∫ (q, t) L dt, {\ displaystyle S (q, t): = \ int ^ {(q, t)} {\ mathcal {L}} \, dt,}{\ displaystyle S (q, t): = \ int ^ {(q, t)} {\ mathcal {L}} \, dt,}

где L {\ displaystyle {\ mathcal { L}}}{\ mathcal { L}} - лагранжиан системы, и где траектория удовлетворяет уравнению Эйлера – Лагранжа системы,

ddt ∂ L ∂ q ˙ = ∂ L ∂ q. {\ displaystyle {\ frac {d} {dt}} {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} = {\ frac {\ partial {\ mathcal { L}}} {\ partial q}}.}{\ displaystyle {\ frac {d} {dt}} {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q} }}} = {\ frac {\ partial {\ mathcal {L}}} {\ partial q}}.}

Вычисление вариации S {\ displaystyle S}S относительно вариации координаты конечной точки,

δ S = ∫ (∂ L ∂ q δ q + ∂ L ∂ q ˙ δ q ˙) dt = ∫ (ddt ∂ L ∂ q ˙ δ q + ∂ L ∂ q ˙ ddt δ q) dt = ∫ ddt (∂ L ∂ q ˙ δ q) dt знак равно ∂ L ∂ q ˙ δ q знак равно п δ q, {\ displaystyle \ delta S = \ int \ left ({\ frac {\ partial {\ mathcal {L}}} {\ partial q}} \ delta q + {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta {\ dot {q}} \ right) dt = \ int \ left ({\ frac {d} {dt}} {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta q + {\ frac {\ partial {\ mathcal {L}} } {\ partial {\ dot {q}}}} {\ frac {d} {dt}} \ delta {q} \ right) dt = \ int {\ frac {d} {dt}} \ left ({\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta q \ right) dt = {\ frac {\ partial {\ mathcal {L}}} {\ partial { \ dot {q}}}} \ delta q = p \ delta q,}{\ displaystyle \ delta S = \ int \ left ({\ frac {\ partial {\ mathcal {L}}} {\ partial q}} \ delta q + {\ frac {\ partial { \ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta {\ dot {q}} \ right) dt = \ int \ left ({\ frac {d} {dt}} {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta q + {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q} }}} {\ frac {d} {dt}} \ delta {q} \ right) dt = \ int {\ frac {d} {dt}} \ left ({\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}}} \ delta q \ right) dt = {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta q = p \ delta q,}

приводит к

∂ S ∂ q = p. {\ displaystyle {\ frac {\ partial S} {\ partial q}} = p.}{\ displaystyle {\ frac {\ partial S} {\ partial q}} = p.}

Использование этого результата и вычисление вариации S {\ displaystyle S}S относительно изменение времени конечной точки приводит непосредственно к уравнению Гамильтона – Якоби,

δ S = L δ t + ∂ L ∂ q ˙ δ q = L δ t - ∂ L ∂ q ˙ q ˙ δ t = - ЧАС δ T, {\ displaystyle \ delta S = {\ mathcal {L}} \ delta t + {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta q = {\ mathcal {L}} \ delta t - {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}}} {\ dot {q}} \ delta t = -H \ delta t \ ;,}{\ displaystyle \ delta S = {\ mathcal { L}} \ delta t + {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} \ delta q = {\ mathcal {L}} \ delta t - {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} {\ dot {q}} \ delta t = -H \ delta t \ ;,}

или

∂ S ∂ t = - H, {\ displaystyle {\ frac {\ partial S} {\ partial t}} = - H,}{\ displaystyle {\ frac {\ partial S} {\ partial t}} = - H,}

где δ q = - q ˙ δ t {\ displaystyle \ delta q = - {\ dot {q}} \ delta t}{\ displaystyle \ delta q = - {\ dot {q }} \ delta t} - изменение траектории в старой конечной точке из-за временной сдвиг и где H = ∂ L ∂ q ˙ q ˙ - L {\ displaystyle H = {\ frac {\ partial {\ mathcal {L}}} {\ partial {\ dot {q}}}} {\ dot {q}} - {\ mathcal {L}}}{\ displaystyle H = {\ frac {\ partial {\ mathcal { L}}} { \ partial {\ dot {q}}}} {\ dot {q}} - {\ mathcal {L}}} - гамильтониан системы.

В качестве альтернативы, как описано ниже, уравнение Гамильтона – Якоби можно вывести из гамильтоновой механики, рассматривая S {\ displaystyle S}{\ displaystyle S} как производящая функция для канонического преобразования классического гамильтониана

H = H (q 1, q 2,…, q N; p 1, p 2,…, p N; t). {\ displaystyle H = H (q_ {1}, q_ {2}, \ ldots, q_ {N}; p_ {1}, p_ {2}, \ ldots, p_ {N}; t).}{\ displaystyle H = H (q_ {1}, q_ {2}, \ ldots, q_ {N}; p_ {1}, p_ {2}, \ ldots, p_ {N}; t). }

Сопряженные импульсы соответствуют первым производным S {\ displaystyle S}{\ displaystyle S} по обобщенным координатам

pk = ∂ S ∂ qk. {\ displaystyle p_ {k} = {\ frac {\ partial S} {\ partial q_ {k}}}.}p_k = \ frac {\ partial S } {\ partial q_k}.

В качестве решения уравнения Гамильтона – Якоби главная функция содержит N + 1 {\ displaystyle N + 1}N + 1 неопределенные константы, первые N {\ displaystyle N}N из них обозначаются как α 1, α 2,..., α N {\ displaystyle \ alpha _ {1}, \, \ alpha _ {2},..., \ alpha _ {N}}{\ displaystyle \ alpha _ {1}, \, \ alpha _ {2 },..., \ alpha _ {N}} , а последний получен в результате интеграции ∂ S ∂ t {\ displaystyle {\ frac {\ partial S} {\ partial t}}}{ \ displaystyle {\ frac {\ partial S} {\ partial t}}} .

Связь между p {\ displaystyle {\ textbf {p}}}{\ textbf {p }} и q {\ displaystyle {\ textbf {q}}}{\ displaystyle {\ textbf {q}}} затем описывает орбиту в фазовом пространстве в терминах этих констант движения. Кроме того, величины

β K = ∂ S ∂ α K, k = 1, 2,…, N {\ displaystyle \ beta _ {k} = {\ frac {\ partial S} {\ partial \ alpha _ { k}}}, \ quad k = 1,2, \ ldots, N}{\ displaystyle \ beta _ {k} = {\ frac { \ partial S} {\ par tial \ alpha _ {k}}}, \ quad k = 1,2, \ ldots, N}

также являются константами движения, и эти уравнения можно инвертировать, чтобы найти q {\ displaystyle {\ textbf {q}}}{\ displaystyle {\ textbf {q}}} как функция всех констант α {\ displaystyle \ alpha}\ альфа и β {\ displaystyle \ beta}\ beta и времени.

Сравнение с другими формулировками механики

HJE - это одно уравнение в частных производных первого порядка для функции обобщенных координат N {\ displaystyle N}N q 1, q 2,..., q N {\ displaystyle q_ {1}, \, q_ {2},..., q_ {N}}{\ displaystyle q_ {1}, \, q_ {2},..., q_ {N}} и время t {\ displaystyle t}{\ displaystyle t} . Обобщенные импульсы не отображаются, за исключением производных от S {\ displaystyle S}{\ displaystyle S} . Примечательно, что функция S {\ displaystyle S}{\ displaystyle S} эквивалентна классическому действию.

. Для сравнения, в эквивалентных уравнениях движения Эйлера – Лагранжа лагранжева механика, сопряженные импульсы тоже не появляются; однако эти уравнения представляют собой систему N {\ displaystyle N}N , как правило, уравнений второго порядка для временной эволюции обобщенных координат. Точно так же уравнения движения Гамильтона представляют собой другую систему из 2N уравнений первого порядка для временной эволюции обобщенных координат и их сопряженных импульсов p 1, p 2,..., p N {\ displaystyle p_ {1}, \, p_ {2},..., p_ {N}}{\ displaystyle p_ {1}, \, p_ {2},..., p_ {N}} .

Поскольку HJE является эквивалентным выражением интегральной задачи минимизации, такой как принцип Гамильтона, HJE может быть полезен в других задачах вариационного исчисления и, в более общем плане, в других разделах математики и физики, таких как динамические системы, симплектическая геометрия и квантовый хаос. Например, уравнения Гамильтона – Якоби можно использовать для определения геодезических на римановом многообразии, важной вариационной задачи в римановой геометрии.

Вывод с использованием канонического преобразования

Любое каноническое преобразование, включающее производящую функцию типа 2 G 2 (q, P, t) {\ displaystyle G_ {2} ({\ textbf {q}}, {\ textbf {P}}, t)}{\ displaystyle G_ {2} ({\ textbf {q}}, {\ textbf {P}}, t)} приводит к соотношениям

p = ∂ G 2 ∂ q, Q = ∂ G 2 ∂ П, К (Q, P, T) знак равно ЧАС (д, р, t) + ∂ G 2 ∂ T {\ Displaystyle \ mathbf {p} = {\ partial G_ {2} \ over \ partial \ mathbf {q} }, \ quad \ mathbf {Q} = {\ partial G_ {2} \ over \ partial \ mathbf {P}}, \ quad K (\ mathbf {Q}, \ mathbf {P}, t) = H (\ mathbf {q}, \ mathbf {p}, t) + {\ partial G_ {2} \ over \ partial t}}{\ displaystyle \ mathbf {p} = {\ partial G_ {2} \ over \ partial \ mathbf {q}}, \ quad \ mathbf {Q} = {\ partial G_ {2} \ over \ partial \ mathbf {P}}, \ quad K (\ mathbf {Q}, \ mathbf { P}, t) знак равно ЧАС (\ mathbf {q}, \ mathbf {p}, t) + {\ partial G_ {2} \ over \ partial t}}

и уравнения Гамильтона в терминах новых переменных P, Q {\ displaystyle \ mathbf {P}, \, \ mathbf {Q}}{\ displaystyle \ mathbf {P}, \, \ mathbf {Q}} и новый гамильтониан K {\ displaystyle K}К имеют одинаковую форму:

P ˙ = - ∂ К ∂ Q, Q ˙ = + ∂ K ∂ P. {\ displaystyle {\ dot {\ mathbf {P}}} = - {\ partial K \ over \ partial \ mathbf {Q}}, \ quad {\ dot {\ mathbf {Q}}} = + {\ partial K \ over \ partial \ mathbf {P}}.}{\ displaystyle {\ dot {\ mathbf {P}}} = - {\ partial K \ over \ partial \ mathbf {Q}}, \ quad {\ dot {\ mathbf {Q}}} = + {\ partial K \ over \ partial \ mathbf {P}}.}

Чтобы получить HJE, мы выбираем производящую функцию G 2 (q, P, t) {\ displaystyle G_ {2} ({\ textbf {q }}, {\ textbf {P}}, t)}{\ displaystyle G_ {2} ({\ textbf {q}}, {\ textbf {P}}, t)} таким образом, чтобы новый гамильтониан был K = 0 {\ displaystyle K = 0}K=0. Следовательно, все его производные также равны нулю, и преобразованные уравнения Гамильтона становятся тривиальными

P ˙ = Q ˙ = 0 {\ displaystyle {\ dot {\ mathbf {P}}} = {\ dot {\ mathbf {Q} }} = 0}\ dot {\ mathbf {P}} = \ dot {\ mathbf {Q}} = 0

, поэтому новые обобщенные координаты и импульсы являются константами движения. Поскольку они являются константами, в этом контексте новые обобщенные импульсы P {\ displaystyle {\ textbf {P}}}{ \ displaystyle {\ textbf {P}}} обычно обозначаются α 1, α 2,..., α N {\ displaystyle \ alpha _ {1}, \, \ alpha _ {2},..., \ alpha _ {N}}{\ displaystyle \ alpha _ {1}, \, \ alpha _ {2 },..., \ alpha _ {N}} , т.е. P m = α m { \ displaystyle P_ {m} = \ alpha _ {m}}{\ displaystyle P_ {m} = \ alpha _ {m}} и новые обобщенные координаты Q {\ displaystyle {\ textbf {Q}}}{\ displaystyle {\ textbf {Q}}} обычно обозначаются как β 1, β 2,..., β N {\ displaystyle \ beta _ {1}, \, \ beta _ {2},..., \ beta _ {N}}{\ displaystyle \ beta _ {1}, \, \ beta _ {2},..., \ beta _ {N}} , поэтому Q m = β m { \ displaystyle Q_ {m} = \ beta _ {m}}{\ displaystyle Q_ {m} = \ beta _ {m}} .

Установка производящей функции, равной главной функции Гамильтона, плюс произвольная константа A {\ displaystyle A}A :

G 2 (q, α, t) знак равно S (q, t) + A, {\ displaystyle G_ {2} (\ mathbf {q}, {\ boldsymbol {\ alpha}}, t) = S (\ mathbf {q}, t) + A,}{\ displaystyle G_ {2} (\ mathbf {q}, {\ boldsymbol {\ alpha}}, t) = S (\ mathbf {q}, t) + A,}

автоматически возникает HJE

p = ∂ G 2 ∂ q = ∂ S ∂ q → H (q, p, t) + ∂ G 2 ∂ t = 0 → H (q, ∂ S ∂ q, t) + ∂ S ∂ T = 0. {\ displaystyle \ mathbf {p} = {\ frac {\ partial G_ {2}} {\ partial \ mathbf {q}}} = {\ frac {\ partial S} {\ partial \ mathbf {q}}} \, \ rightarrow \, H (\ mathbf {q}, \ mathbf {p}, t) + {\ partial G_ {2} \ over \ partial t} = 0 \, \ rightarrow \, H \ left (\ mathbf {q}, {\ frac {\ partial S} {\ partial \ mathbf {q}}}, t \ right) + {\ partial S \ over \ partial t} = 0.}{\ displaystyle \ mathbf {p} = {\ frac {\ partial G_ {2 }} {\ partial \ mathbf {q}}} = {\ frac {\ partial S} {\ partial \ mathbf {q}}} \, \ rightarrow \, H (\ mathbf {q}, \ mathbf {p}, t) + {\ partial G_ {2} \ over \ partial t} = 0 \, \ rightarrow \, H \ left (\ mathbf {q}, {\ frac {\ partial S} {\ partial \ mathbf {q }}}, t \ right) + {\ partial S \ over \ partial t} = 0.}

Как только мы решили для S (q, α, t) {\ displaystyle S (\ mathbf {q}, {\ boldsymbol {\ alpha}}, t)}{\ displaystyle S (\ mathbf {q}, {\ boldsymbol {\ alpha}}, t)} , они также дают нам полезные уравнения

Q = β = ∂ S ∂ α, {\ displaystyle \ mathbf {Q} = {\ boldsymbol {\ beta}} = {\ partial S \ over \ partial {\ boldsymbol {\ alpha}}},}{\ displaystyle \ mathbf {Q} = {\ boldsymbol {\ beta}} = {\ частичный S \ over \ partial {\ boldsymbol {\ alpha}}},}

или записано в компонентах для ясности

Q m = β m = ∂ S (q, α, t) ∂ α m. {\ displaystyle Q_ {m} = \ beta _ {m} = {\ frac {\ partial S (\ mathbf {q}, {\ boldsymbol {\ alpha}}, t)} {\ partial \ alpha _ {m} }}.}{ \ Displaystyle Q _ {m} = \ beta _ {m} = {\ frac {\ partial S (\ mathbf {q}, {\ boldsymbol {\ alpha}}, t)} {\ partial \ alpha _ {m}}}. }

В идеале эти N уравнений можно инвертировать, чтобы найти исходные обобщенные координаты q {\ displaystyle {\ textbf {q}}}{\ textbf {q}} как функция констант α, β, {\ displaystyle {\ boldsymbol {\ alpha}}, \, {\ boldsymbol {\ beta}},}{\ displaystyle {\ boldsymbol {\ alpha}}, \, {\ boldsymbol {\ beta}},} и t {\ displaystyle t }t , таким образом решая исходную проблему.

Действие и функции Гамильтона

Основная функция Гамильтона S и классическая функция H тесно связаны с действием. Полная разность для S {\ displaystyle S}S составляет:

d S = ∑ i ∂ S ∂ qidqi + ∂ S ∂ tdt {\ displaystyle dS = \ sum _ {i} {\ frac {\ partial S} {\ partial q_ {i}}} dq_ {i} + {\ frac {\ partial S} {\ partial t}} dt}{\ displaystyle dS = \ sum _ {i} {\ frac {\ partial S} {\ partial q_ {i}}} dq_ {i} + {\ frac {\ partial S} {\ partial t}} dt}

, поэтому производная по времени от S равна

d S dt = ∑ i ∂ S ∂ qiq ˙ i + ∂ S ∂ t = ∑ ipiq ˙ i - H = L. {\ displaystyle {\ frac {dS} {dt}} = \ sum _ {i} {\ frac {\ partial S} {\ partial q_ {i}}} {\ dot {q}} _ {i} + { \ frac {\ partial S} {\ partial t}} = \ sum _ {i} p_ {i} {\ dot {q}} _ {i} -H = L.}{\ displaystyle {\ frac {dS} {dt}} = \ sum _ {i} {\ frac { \ partial S} {\ partial q_ {i}}} {\ dot {q}} _ {i} + {\ frac {\ partial S} {\ partial t}} = \ sum _ {i} p_ {i} {\ точка {q}} _ {i} -H = L.}

Следовательно,

S = ∫ L dt, {\ displaystyle S = \ int L \, dt,}{\ displaystyle S = \ int L \, dt,}

так что S на самом деле классическое действие плюс неопределенная константа.

Когда H явно не зависит от времени,

W = S + E t = S + H t = ∫ (L + H) dt = ∫ p ⋅ dq, {\ displaystyle W = S + Et = S + Ht = \ int (L + H) \, dt = \ int \ mathbf {p} \ cdot d \ mathbf {q},}{\ Displaystyle W = S + Et = S + Ht = \ int (L + H) \, dt = \ int \ mathbf {p} \ cdot d \ mathbf {q},}

в этом случае W то же самое, что и сокращенное действие.

Разделение переменных

HJE наиболее полезен, когда его можно решить с помощью аддитивного разделения переменных, которое напрямую определяет константы движения. Например, время t можно разделить, если гамильтониан не зависит явно от времени. В этом случае производная по времени ∂ S ∂ t {\ displaystyle {\ frac {\ partial S} {\ partial t}}}{ \ displaystyle {\ frac {\ partial S} {\ partial t}}} в HJE должна быть константой, обычно обозначаемой (- E {\ displaystyle -E}{\ displaystyle -E} ), что дает разделенное решение

S = W (q 1, q 2,…, q N) - E t {\ displaystyle S = W ( q_ {1}, q_ {2}, \ ldots, q_ {N}) - Et}{\ displaystyle S = W (q_ {1}, q_ {2}, \ ldots, q_ {N}) - Et}

где не зависящая от времени функция W (q) {\ displaystyle W ({\ textbf {q}}) }{\ displaystyle W ({\ textbf {q}})} иногда называют характеристической функцией Гамильтона . Тогда редуцированное уравнение Гамильтона – Якоби может быть записано

H (q, ∂ S ∂ q) = E. {\ displaystyle H \ left (\ mathbf {q}, {\ frac {\ partial S} {\ partial \ mathbf {q}}} \ right) = E.}{\ displaystyle H \ left (\ mathbf {q}, {\ frac {\ partial S} {\ partial \ mathbf {q}}} \ right) = E.}

Чтобы проиллюстрировать разделимость для других переменных, мы предполагаем что некоторая обобщенная координата qk {\ displaystyle q_ {k}}{\ displaystyle q_ {k}} и ее производная ∂ S ∂ qk {\ displaystyle {\ frac {\ partial S} {\ partial q_ {k}}}}{\ displaystyle {\ frac {\ partial S} {\ частичный q_ {k}}}} появляются вместе как одна функция

ψ (qk, ∂ S ∂ qk) {\ displaystyle \ psi \ left (q_ {k}, {\ frac {\ partial S} {\ partial q_ {k}}} \ right)}\ psi \ left (q_k, \ frac {\ partial S} {\ partial q_k} \ right)

в гамильтониане

H = H (q 1, q 2,…, qk - 1, qk + 1,…, q N; p 1, p 2,…, pk - 1, pk + 1,…, p N; ψ; t). {\ displaystyle H = H (q_ {1}, q_ {2}, \ ldots, q_ {k-1}, q_ {k + 1}, \ ldots, q_ {N}; p_ {1}, p_ {2) }, \ ldots, p_ {k-1}, p_ {k + 1}, \ ldots, p_ {N}; \ psi; t).}{\ стиль отображения H = H (q_ {1}, q_ {2}, \ ldots, q_ {k-1}, q_ {k + 1}, \ ldots, q_ {N}; p_ {1}, p_ {2}, \ ldots, p_ {k-1}, p_ {k + 1}, \ ldots, p_ {N}; \ psi; t).}

В этом случае функцию S можно разделить на две функции, одно зависит только от q k, а другое зависит только от оставшихся обобщенных координат

S = S k (qk) + S rem (q 1,…, qk - 1, qk + 1,…, q N, t). {\ Displaystyle S = S_ {k} (q_ {k}) + S _ {\ text {rem}} (q_ {1}, \ ldots, q_ {k-1}, q_ {k + 1}, \ ldots, q_ {N}, t).}{\ displaystyle S = S_ { k} (q_ {k}) + S _ {\ text {rem}} (q_ {1}, \ ldots, q_ {k-1}, q_ {k + 1}, \ ldots, q_ {N}, t).}

Подстановка этих формул в уравнение Гамильтона – Якоби показывает, что функция ψ должна быть константой (обозначенной здесь как Γ k {\ displaystyle \ Gamma _ {k}}\ Gamma_k ), что дает обыкновенное дифференциальное уравнение первого порядка для S k (qk), {\ displaystyle S_ {k} (q_ {k}),}{\ displaystyle S_ {k} (q_ {k}),}

ψ (qk, d S kdqk) = Γ k. {\ displaystyle \ psi \ left (q_ {k}, {\ frac {dS_ {k}} {dq_ {k}}} \ right) = \ Gamma _ {k}.}{\ displaystyle \ psi \ left (q_ {k}, {\ frac {dS_ {k}} {dq_ {k}) }} \ right) = \ Gamma _ {k}.}

В удачных случаях функция S {\ displaystyle S}S можно полностью разделить на N {\ displaystyle N}N функции S m (qm), {\ displaystyle S_ { m} (q_ {m}),}{\ displaystyle S_ {m} (q_ {m}),}

S = S 1 (q 1) + S 2 (q 2) + ⋯ + SN (q N) - E t. {\ displaystyle S = S_ {1} (q_ {1}) + S_ {2} (q_ {2}) + \ cdots + S_ {N} (q_ {N}) - Et.}S = S_1 (q_1) + S_2 (q_2) + \ cdots + S_N (q_N) -Et.

В таком В этом случае проблема переходит в N {\ displaystyle N}N обыкновенные дифференциальные уравнения.

Разделимость S зависит как от гамильтониана, так и от выбора обобщенных координат. Для ортогональных координат и гамильтонианов, которые не зависят от времени и являются квадратичными по обобщенным импульсам, S {\ displaystyle S}S будет полностью разделимым, если потенциальная энергия аддитивно отделима в каждой координате, где член потенциальной энергии для каждой координаты умножается на зависящий от координат множитель в соответствующем члене импульса гамильтониана (условия Штекеля ). Для иллюстрации в следующих разделах работают несколько примеров в ортогональных координатах.

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

Сферические координаты

В сферических координатах гамильтониан свободной частицы, движущейся в консервативном потенциале U, может быть записан

H = 1 2 m [пр 2 + p θ 2 r 2 + p ϕ 2 r 2 sin 2 ⁡ θ] + U (r, θ, ϕ). {\ displaystyle H = {\ frac {1} {2m}} \ left [p_ {r} ^ {2} + {\ frac {p _ {\ theta} ^ {2}} {r ^ {2}}} + {\ frac {p _ {\ phi} ^ {2}} {r ^ {2} \ sin ^ {2} \ theta}} \ right] + U (r, \ theta, \ phi).}H = \ frac {1} {2m} \ left [p_ {r} ^ {2} + \ frac {p _ {\ theta} ^ {2}} {r ^ {2}} + \ frac {p _ {\ phi} ^ {2}} {r ^ {2} \ sin ^ {2} \ theta} \ справа] + U (г, \ theta, \ phi).

Уравнение Гамильтона – Якоби полностью разделимо в этих координатах, если существуют функции U r (r), U θ (θ) и U ϕ (ϕ) такое, что U можно записать в аналогичном виде

U (r, θ, ϕ) = U r (r) + U θ (θ) r 2 + U ϕ (ϕ) r 2 sin 2 ⁡ θ. {\ Displaystyle U (г, \ theta, \ phi) = U_ {r} (r) + {\ frac {U _ {\ theta} (\ theta)} {r ^ {2}}} + {\ frac {U_ {\ phi} (\ phi)} {r ^ {2} \ sin ^ {2} \ theta}}.}U (r, \ theta, \ phi) = U_ {r} (r) + \ frac {U _ {\ theta} (\ theta)} {r ^ {2}} + \ frac {U _ {\ phi} (\ phi)} {r ^ {2} \ sin ^ {2} \ theta}.

Замена полностью разделенного решения

S = S r (r) + S θ ( θ) + S ϕ (ϕ) - E t {\ displaystyle S = S_ {r} (r) + S _ {\ theta} (\ theta) + S _ {\ phi} (\ phi) -Et}S = S_ {r} (г) + S _ {\ theta} (\ theta) + S _ {\ phi} (\ phi) - Et

в HJE дает

1 2 m (d S rdr) 2 + U r (r) + 1 2 mr 2 [(d S θ d θ) 2 + 2 m U θ (θ)] + 1 2 mr 2 sin 2 ⁡ θ [(d S ϕ d ϕ) 2 + 2 m U ϕ (ϕ)] = E. {\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {r}} {dr}} \ right) ^ {2} + U_ {r} (r) + {\ frac {1 } {2mr ^ {2}}} \ left [\ left ({\ frac {dS _ {\ theta}} {d \ theta}} \ right) ^ {2} + 2mU _ {\ theta} (\ theta) \ right ] + {\ frac {1} {2mr ^ {2} \ sin ^ {2} \ theta}} \ left [\ left ({\ frac {dS _ {\ phi}} {d \ phi}} \ right) ^ {2} + 2mU _ {\ phi} (\ phi) \ right] = E.}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {r}} {dr}} \ справа) ^ {2} + U_ {r} (r) + {\ frac {1} {2mr ^ {2}}} \ left [\ left ({\ frac {dS _ {\ theta}} {d \ theta} } \ right) ^ {2} + 2mU _ {\ theta} (\ theta) \ right] + {\ frac {1} {2mr ^ {2} \ sin ^ {2} \ theta}} \ left [\ left ( {\ frac {dS _ {\ phi}} {d \ phi}} \ right) ^ {2} + 2mU _ {\ phi} (\ phi) \ right] = E.}

Это уравнение может быть решено путем последовательного интегрирования обыкновенных дифференциальных уравнений, начиная с уравнения для ϕ {\ displaystyle \ phi}{\ displaystyle \ phi}

(d S ϕ d ϕ) 2 + 2 м U ϕ (ϕ) = Γ ϕ {\ displaystyle \ left ({\ frac {dS _ {\ phi}} {d \ phi}) } \ right) ^ {2} + 2mU _ {\ phi} (\ phi) = \ Gamma _ {\ phi}}{\ displaystyle \ left ({\ frac {dS _ {\ phi}} {d \ phi}} \ right) ^ {2} + 2mU _ {\ phi} (\ phi) = \ Гамма _ {\ phi}}

где Γ ϕ {\ displaystyle \ Gamma _ {\ phi}}{\ displaystyle \ Gamma _ {\ phi}} - константа движения, которая устраняет зависимость ϕ {\ displaystyle \ phi}{\ displaystyle \ phi} из уравнения Гамильтона – Якоби

1 2 m (d S rdr) 2 + U r (r) + 1 2 mr 2 [(d S θ d θ) 2 + 2 m U θ (θ) + Γ ϕ sin 2 ⁡ θ] = E. {\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {r}} {dr}} \ right) ^ {2} + U_ {r} (r) + {\ frac {1 } {2mr ^ {2}}} \ left [\ left ({\ frac {dS _ {\ theta}} {d \ theta}} \ right) ^ {2} + 2mU _ {\ theta} (\ theta) + { \ frac {\ Gamma _ {\ phi}} {\ sin ^ {2} \ theta}} \ right] = E.}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {r}} {dr}} \ right) ^ {2} + U_ {r} (r) + {\ frac {1 } {2mr ^ {2}}} \ left [\ left ({\ frac {dS _ {\ theta}} {d \ theta}} \ right) ^ {2} + 2mU _ {\ theta} (\ theta) + { \ frac {\ Gamma _ {\ phi}} {\ sin ^ {2} \ theta}} \ right] = E.}

Следующее обыкновенное дифференциальное уравнение включает θ { \ displaystyle \ theta}{\ displaystyle \ theta} обобщенная координата

(d S θ d θ) 2 + 2 m U θ (θ) + Γ ϕ sin 2 ⁡ θ = Γ θ {\ displaystyle \ left ({\ frac {dS_ {\ theta}} {d \ theta}} \ right) ^ {2} + 2mU _ {\ theta} (\ theta) + {\ frac {\ Gamma _ {\ phi}} {\ sin ^ {2} \ theta }} = \ Gamma _ {\ theta}}{\ displaystyle \ left ({\ frac {dS _ {\ theta}} {d \ theta}} \ right) ^ {2} + 2mU _ {\ theta} (\ theta) + {\ frac {\ Gamma _ {\ phi}} {\ sin ^ {2} \ theta}} = \ Gamma _ {\ theta}}

где Γ θ {\ displaystyle \ Gamma _ {\ theta}}{\ displaystyle \ Gamma _ {\ theta}} снова является константой движения который устраняет зависимость θ {\ displaystyle \ theta}{\ displaystyle \ theta} и сводит HJE к окончательному обыкновенному дифференциальному уравнению

1 2 m (d S rdr) 2 + U r (r) + Γ θ 2 mr 2 = E {\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {r}} {dr}} \ right) ^ {2} + U_ {r} (r) + {\ frac {\ Gamma _ {\ theta}} {2mr ^ {2}}} = E}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {r}} {dr}} \ right) ^ {2} + U_ {r} (r) + {\ frac {\ Gamma _ {\ theta}} {2mr ^ {2}}} = E}

, чья в Интеграция завершает решение для S {\ displaystyle S}{\ displaystyle S} .

Эллиптические цилиндрические координаты

Гамильтониан в эллиптических цилиндрических координатах можно записать

H = p μ 2 + п ν 2 2 ма 2 (зп 2 ⁡ μ + грех 2 ⁡ ν) + pz 2 2 m + U (μ, ν, z) {\ displaystyle H = {\ frac {p _ {\ mu} ^ {2} + p _ {\ nu} ^ {2}} {2ma ^ {2} \ left (\ sinh ^ {2} \ mu + \ sin ^ {2} \ nu \ right)}} + {\ frac {p_ {z} ^ {2}} {2m}} + U (\ mu, \ nu, z)}H = \ frac {p _ {\ mu} ^ {2} + p _ {\ nu} ^ {2}} {2ma ^ {2} \ left (\ sinh ^ {2} \ mu + \ sin ^ {2 } \ nu \ right)} + \ frac {p_ {z} ^ {2}} {2m} + U (\ mu, \ nu, z)

, где фокусы из эллипсов расположены в ± a {\ displaystyle \ pm a}\ pm a на оси x {\ displaystyle x}x . Уравнение Гамильтона – Якоби полностью разделимо в этих координатах при условии, что U {\ displaystyle U}{\ displaystyle U} имеет аналогичную форму

U (μ, ν, z) = U μ (μ) + U ν (ν) зп 2 ⁡ μ + грех 2 ⁡ ν + U Z (z) {\ displaystyle U (\ mu, \ nu, z) = {\ frac {U _ {\ mu} (\ mu) + U_ { \ nu} (\ nu)} {\ sinh ^ {2} \ mu + \ sin ^ {2} \ nu}} + U_ {z} (z)}U (\ mu, \ nu, z) = \ frac {U _ {\ mu} (\ mu) + U _ {\ nu } (\ nu)} {\ sinh ^ {2} \ mu + \ sin ^ {2} \ nu} + U_ {z} (z)

где: U μ (μ) {\ Displaystyle U _ {\ mu} (\ mu)}{\ displaystyle U _ {\ mu} (\ mu)} , U ν (ν) {\ displaystyle U _ {\ nu} (\ nu)}{\ displaystyle U _ {\ nu} (\ nu)} и U z (z) { \ displaystyle U_ {z} (z)}{\ displaystyle U_ {z} (z)} - произвольные функции. Замена полностью разделенного раствора

S = S μ (μ) + S ν (ν) + S z (z) - E t {\ displaystyle S = S _ {\ mu} (\ mu) + S _ {\ nu } (\ nu) + S_ {z} (z) -Et}S = S _ {\ mu} (\ mu) + S _ {\ nu} (\ nu) + S_ {z} (z) - Et в HJE дает
1 2 m (d S zdz) 2 + U z (z) + 1 2 ma 2 ( sinh 2 ⁡ μ + sin 2 ⁡ ν) [(d S μ d μ) 2 + (d S ν d ν) 2 + 2 ma 2 U μ (μ) + 2 ma 2 U ν (ν)] = E. {\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} (z) + {\ frac {1 } {2ma ^ {2} \ left (\ sinh ^ {2} \ mu + \ sin ^ {2} \ nu \ right)}} \ left [\ left ({\ frac {dS _ {\ mu}} {d \ mu}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ nu}} {d \ nu}} \ right) ^ {2} + 2ma ^ {2} U _ {\ mu} ( \ mu) + 2ma ^ {2} U _ {\ nu} (\ nu) \ right] = E.}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z }} {dz}} \ right) ^ {2} + U_ {z} (z) + {\ frac {1} {2ma ^ {2} \ left (\ sinh ^ {2} \ mu + \ sin ^ { 2} \ nu \ right)}} \ left [\ left ({\ frac {dS _ {\ mu}} {d \ mu}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ nu }} {d \ nu}} \ right) ^ {2} + 2ma ^ {2} U _ {\ mu} (\ mu) + 2ma ^ {2} U _ {\ nu} (\ nu) \ right] = E.}

Разделение первого обыкновенного дифференциального уравнения

1 2 m (d S zdz) 2 + U z (z) знак равно Γ z {\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} ( z) = \ Gamma _ {z}}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} (z) = \ Gamma _ {z}}

дает сокращенное уравнение Гамильтона – Якоби (после перестановки и умножения обеих частей на знаменатель)

(d S μ d μ) 2 + (d S ν d ν) 2 + 2 ma 2 U μ (μ) + 2 ma 2 U ν (ν) = 2 ma 2 (sinh 2 ⁡ μ + sin 2 ⁡ ν) (E - Γ z) {\ displaystyle \ left ({ \ frac {dS _ {\ mu}} {d \ mu}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ nu}} {d \ nu}} \ right) ^ {2} + 2ma ^ {2} U _ {\ mu} (\ mu) + 2ma ^ {2} U _ {\ nu} (\ nu) = 2ma ^ {2} \ left (\ sinh ^ {2} \ mu + \ sin ^ {2} \ nu \ right) \ left (E- \ Gamma _ {z} \ right)}{\ displaystyle \ left ({\ frac { dS _ {\ mu}} {d \ mu}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ nu}} {d \ nu}} \ right) ^ {2} + 2ma ^ { 2} U _ {\ mu} (\ mu) + 2ma ^ {2} U _ {\ nu} (\ nu) = 2ma ^ {2} \ left (\ sinh ^ {2} \ mu + \ sin ^ {2} \ Nu \ right) \ left (E- \ Gamma _ {z} \ right)}

который сам может быть разделен на tw o независимые обыкновенные дифференциальные уравнения

(d S μ d μ) 2 + 2 ma 2 U μ (μ) + 2 ma 2 (Γ z - E) sinh 2 ⁡ μ = Γ μ {\ displaystyle \ left ({\ frac {dS _ {\ mu}} {d \ mu}} \ right) ^ {2} + 2ma ^ {2} U _ {\ mu} (\ mu) + 2ma ^ {2} \ left (\ Gamma _ {z} -E \ right) \ sinh ^ {2} \ mu = \ Gamma _ {\ mu}}{\ displaystyle \ left ({\ frac {dS _ {\ mu}} {d \ mu}} \ right) ^ {2} + 2ma ^ {2} U _ {\ mu} (\ mu) + 2ma ^ {2} \ left (\ Gamma _ {z} -E \ right) \ sinh ^ {2} \ mu = \ Gamma _ {\ mu}}
(d S ν d ν) 2 + 2 ma 2 U ν (ν) + 2 ma 2 (Γ z - E) грех 2 ⁡ ν = Γ ν {\ displaystyle \ left ({\ frac {dS _ {\ nu}} {d \ nu}} \ right) ^ {2} + 2ma ^ {2} U_ { \ nu} (\ nu) + 2ma ^ {2} \ left (\ Gamma _ {z} -E \ right) \ sin ^ {2} \ nu = \ Gamma _ {\ nu}}{\ displaystyle \ left ({\ frac {dS_ {\ nu}} {d \ nu}} \ right) ^ {2} + 2ma ^ {2} U _ {\ nu} (\ nu) + 2ma ^ {2} \ left (\ Gamma _ {z} -E \ right) \ sin ^ {2} \ nu = \ Gamma _ {\ nu}}

что, когда решено, предоставить полное решение для S {\ displaystyle S}{\ displaystyle S} .

Параболических цилиндрических координат

Гамильтониан в параболических цилиндрических координатах может быть записан в виде

H = p σ 2 + p τ 2 2 m (σ 2 + τ 2) + pz 2 2 m + U (σ, τ, z). {\ displaystyle H = {\ frac {p _ {\ sigma} ^ {2} + p _ {\ tau} ^ {2}} {2m \ left (\ sigma ^ {2} + \ tau ^ {2} \ right) }} + {\ frac {p_ {z} ^ {2}} {2m}} + U (\ sigma, \ tau, z).}H = \ frac {p _ {\ sigma} ^ {2} + p _ {\ tau} ^ {2}} {2m \ left (\ sigma ^ {2} + \ tau ^ {2} \ right)} + \ frac {p_ { z} ^ {2}} {2m} + U (\ sigma, \ tau, z).

Уравнение Гамильтона – Якоби полностью разделимо в этих координатах при условии, что U {\ displaystyle U}{\ displaystyle U} имеет аналогичную форму

U (σ, τ, z) = U σ (σ) + U τ (τ) σ 2 + τ 2 + U z ( z) {\ Displaystyle U (\ sigma, \ tau, z) = {\ frac {U _ {\ sigma} (\ sigma) + U _ {\ tau} (\ tau)} {\ sigma ^ {2} + \ tau ^ {2}}} + U_ {z} (z)}U (\ sigma, \ tau, z) = \ frac {U _ {\ sigma} (\ sigma) + U _ {\ tau} (\ tau)} {\ sigma ^ {2} + \ tau ^ {2}} + U_ {z} (z)

где U σ (σ) {\ displaystyle U _ {\ sigma} (\ sigma)}{\ displaystyle U _ {\ sigma} (\ sigma)} , U τ (τ) {\ displaystyle U _ {\ tau} (\ tau)}{\ displaystyle U _ {\ tau} (\ tau)} и U z (z) {\ displaystyle U_ {z} (z)}{\ displaystyle U_ {z} (z)} - произвольные функции. Замена полностью разделенного решения

S = S σ (σ) + S τ (τ) + S z (z) - E t + constant {\ displaystyle S = S _ {\ sigma} (\ sigma) + S_ { \ tau} (\ tau) + S_ {z} (z) -Et + {\ text {constant}}}{\ displaystyle S = S _ {\ sigma} (\ sigma) + S _ {\ tau} (\ tau) + S_ {z} (z) -Et + {\ text {constant}}}

в HJE дает

1 2 m (d S zdz) 2 + U z (z) + 1 2 м (σ 2 + τ 2) [(d S σ d σ) 2 + (d S τ d τ) 2 + 2 м U σ (σ) + 2 м U τ (τ)] = E. {\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} (z) + {\ frac {1 } {2m \ left (\ sigma ^ {2} + \ tau ^ {2} \ right)}} \ left [\ left ({\ frac {dS _ {\ sigma}} {d \ sigma}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ tau}} {d \ tau}} \ right) ^ {2} + 2mU _ {\ sigma} (\ sigma) + 2mU _ {\ tau} (\ tau) \ right] = E.}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} (z) + {\ frac {1} {2m \ left ( \ sigma ^ {2} + \ tau ^ {2} \ right)}} \ left [\ left ({\ frac {dS _ {\ sigma}} {d \ sigma}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ tau}} {d \ tau}} \ right) ^ {2} + 2mU _ {\ sigma} (\ sigma) + 2mU _ {\ tau} (\ tau) \ right] = E. }

Разделение первого обыкновенного дифференциального уравнения

1 2 m (d S zdz) 2 + U z (z) = Γ z {\ displaystyle {\ frac {1} { 2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} (z) = \ Gamma _ {z}}{\ displaystyle {\ frac {1} {2m}} \ left ({\ frac {dS_ {z}} {dz}} \ right) ^ {2} + U_ {z} (z) = \ Gamma _ {z}}

возвращает сокращенное значение Гамильтона– Уравнение Якоби (после перестановки и умножения обеих частей на знаменатель)

(d S σ d σ) 2 + (d S τ d τ) 2 + 2 m U σ (σ) + 2 m U τ ( τ) знак равно 2 м (σ 2 + τ 2) (E - Γ z) {\ Displaystyle \ влево ({\ гидроразрыва {dS _ {\ sigma}} {d \ sigma}} \ вправо) ^ {2} + \ влево ({\ frac {dS _ {\ tau}} {d \ tau}} \ right) ^ {2} + 2mU _ {\ sigma} (\ sigma) + 2mU _ {\ tau} (\ tau) = 2m \ left (\ sigma ^ {2} + \ tau ^ {2} \ right) \ left (E- \ Gamma _ {z} \ right)}{\ displaystyle \ left ({\ frac {dS _ {\ sigma}} {d \ sigma}} \ right) ^ {2} + \ left ({\ frac {dS _ {\ tau}} {d \ tau}} \ right) ^ {2} + 2mU _ {\ sigma} (\ sigma) + 2mU _ {\ tau} (\ tau) = 2m \ left (\ sigma ^ {2} + \ tau ^ {2} \ right) \ left (E- \ Gamma _ {z} \ right)}

который сам может быть разделен на два независимых обычных разных ial уравнения

(d S σ d σ) 2 + 2 м U σ (σ) + 2 м σ 2 (Γ z - E) = Γ σ {\ displaystyle \ left ({\ frac {dS _ {\ sigma}) } {d \ sigma}} \ right) ^ {2} + 2mU _ {\ sigma} (\ sigma) + 2m \ sigma ^ {2} \ left (\ Gamma _ {z} -E \ right) = \ Gamma _ {\ sigma}}{\ displaystyle \ left ({\ frac {dS _ {\ sigma}} {d \ sigma}} \ right) ^ {2} + 2mU _ {\ sigma} (\ sigma) + 2m \ sigma ^ {2} \ слева (\ Gamma _ {z} -E \ right) = \ Gamma _ {\ sigma}}
(d S τ d τ) 2 + 2 м U τ (τ) + 2 м τ 2 (Γ z - E) = Γ τ {\ displaystyle \ left ({\ frac {dS_ { \ tau}} {d \ tau}} \ right) ^ {2} + 2mU _ {\ tau} (\ tau) + 2m \ tau ^ {2} \ left (\ Gamma _ {z} -E \ right) = \ Gamma _ {\ tau}}{\ displaystyle \ left ({\ frac {dS _ {\ tau}} {d \ tau}} \ справа) ^ {2} + 2mU _ {\ tau} (\ tau) + 2m \ tau ^ {2} \ left (\ Gamma _ {z} -E \ right) = \ Gamma _ {\ tau}}

, которые после решения предоставляют полное решение для S {\ displaystyle S}S .

Волны и частицы

Фронты и траектории оптических волн

HJE устанавливает двойственность между траекториями и волновыми фронтами. Например, в геометрической оптике свет можно рассматривать как «лучи» или как волны. Фронт волны можно определить как поверхность C t {\ textstyle {\ cal {C}} _ ​​{t}}{\ textstyle {\ cal {C}} _ ​​{t}} , на которой свет излучается в момент времени t = 0 {\ textstyle t = 0}{ \ textstyle t = 0} достиг в момент t {\ textstyle t}{\ textstyle t} . Световые лучи и волновые фронты двойственны: если один известен, другой можно вывести.

Точнее, геометрическая оптика - это вариационная задача, где «действие» - это время прохождения T {\ textstyle T}{\ textstyle T} по пути,

T = 1 c ∫ AB nds {\ displaystyle T = {\ frac {1} {c}} \ int _ {A} ^ {B} nds}{\ displaystyle T = {\ frac {1} {c}} \ int _ {A} ^ {B} nds} где n {\ textstyle n}{\ textstyle n} - показатель преломления и ds {\ textstyle ds}{\ textstyle ds} среды представляет собой бесконечно малую длину дуги. Из приведенной выше формулировки можно вычислить траектории лучей, используя формулировку Эйлера-Лагранжа; в качестве альтернативы можно вычислить волновые фронты, решив уравнение Гамильтона-Якоби. Знание одного ведет к знанию другого.

Вышеупомянутая двойственность является очень общей и применима ко всем системам, которые основаны на вариационном принципе: либо вычисляют траектории с помощью уравнений Эйлера-Лагранжа, либо волновые фронты с помощью уравнения Гамильтона-Якоби.

Фронт волны в момент времени t {\ textstyle t}{\ textstyle t} для системы, изначально равной q 0 {\ textstyle {\ textbf {q}} _ {0 }}{\ textstyle {\ textbf {q} } _ {0}} в момент времени t 0 {\ textstyle t_ {0}}{\ textstyle t_ {0}} , определяется как набор точек q {\ textstyle {\ textbf {q} }}{\ textstyle {\ textbf {q}}} такой, что S (q, t) = const {\ textstyle S ({\ textbf {q}}, t) = {\ text {const}}}{\ textstyle S ({\ textbf {q}}, t) = {\ text {const}}} . Если S (q, t) {\ textstyle S ({\ textbf {q}}, t)}{\ textstyle S ({\ textbf {q}}, t)} известен, моментально вычисляется импульс,

p = ∂ S ∂ q. {\ displaystyle {\ textbf {p}} = {\ frac {\ partial S} {\ partial {\ textbf {q}}}}.}{\ displaystyle {\ textbf {p}} = {\ frac {\ partial S} {\ partial {\ textbf {q}}}}.} Один раз p {\ textstyle {\ textbf {p} }}{\ textstyle {\ textbf {p}}} известно, касательные к траекториям q ˙ {\ textstyle {\ dot {\ textbf {q}}}}{\ textstyle {\ dot {\ textbf {q}}}} вычисляются путем решения уравнения ∂ L ∂ q ˙ знак равно п {\ Displaystyle {\ frac {\ partial {\ cal {L}}} {\ partial {\ dot {\ textbf {q}}}}} = {\ boldsymbol {p}}}{\ displaystyle {\ frac {\ partial {\ cal {L}}} {\ partial {\ dot {\ textbf {q}}}}} = {\ boldsymbol {p}}} для q ˙ {\ textstyle {\ dot {\ textbf {q}}}}{\ textstyle {\ dot {\ textbf {q}}}} , где L {\ textstyle {\ cal {L}}}{\ textstyle {\ cal {L}}} - лагранжиан. Затем траектории восстанавливаются на основе знания q ˙ {\ textstyle {\ dot {\ textbf {q}}}}{\ textstyle {\ dot {\ textbf {q}}}} .

Связи с уравнением Шредингера

изоповерхностей функции S (q, t) {\ displaystyle S ({\ textbf {q}}, t)}{\ displaystyle S ({\ textbf {q}}, t)} может быть определено в любое время t. Движение изоповерхности S {\ displaystyle S}S как функция времени определяется движениями частиц, начинающимися в точках q {\ displaystyle {\ textbf {q }}}{\ displaystyle {\ textbf {q}}} на изоповерхности. Движение такой изоповерхности можно представить как волну , движущуюся через q {\ displaystyle {\ textbf {q}}}{\ displaystyle {\ textbf {q}}} -пространство, хотя она не подчиняется волновое уравнение точно. Чтобы показать это, пусть S представляет собой фазу волны

ψ = ψ 0 ei S / ℏ {\ displaystyle \ psi = \ psi _ {0} e ^ {iS / \ hbar}}\ psi = \ psi_ {0} e ^ {iS / \ hbar}

где ℏ {\ displaystyle \ hbar}\ hbar - константа (постоянная Планка ), введенная, чтобы сделать экспоненциальный аргумент безразмерным; изменения амплитуды волны можно представить, если S {\ displaystyle S}S будет комплексным числом. We may then rewrite the Hamilton–Jacobi equation as

ℏ 2 2 m ∇ 2 ψ − U ψ = ℏ i ∂ ψ ∂ t {\displaystyle {\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi -U\psi ={\frac {\hbar }{i}}{\frac {\partial \psi }{\partial t}}}{\ displaystyle {\ frac {\ hbar ^ {2}} { 2m}} \ nabla ^ {2} \ psi -U \ psi = {\ frac {\ hbar} {i}} {\ frac {\ partial \ psi} {\ partial t}}}

which is the Schrödinger equation.

Conversely, starting with the Schrödinger equation and our ansatz for ψ {\displaystyle \psi }\ psi , we arrive at

1 2 m ( ∇ S) 2 + U + ∂ S ∂ t = i ℏ 2 m ∇ 2 S. {\displaystyle {\frac {1}{2m}}\left(\nabla S\right)^{2}+U+{\frac {\partial S}{\partial t}}={\frac {i\hbar }{2m}}\nabla ^{2}S.}\ frac {1} {2m} \ left (\ nabla S \ right) ^ {2} + U + \ frac {\ partial S} {\ partial t } = \ frac {i \ hbar} {2m} \ nabla ^ {2} S.

The classical limit (ℏ → 0 {\displaystyle \hbar \rightarrow 0}\ hbar \ rightarrow 0 ) of the Schrödinger equation above becomes identical to the following variant of the Hamilton–Jacobi equation,

1 2 m ( ∇ S) 2 + U + ∂ S ∂ t = 0. {\displaystyle {\frac {1}{2m}}\left(\nabla S\right)^{2}+U+{\frac {\partial S}{\partial t}}=0.}\ frac {1} {2m} \ left (\ nabla S \ right) ^ {2} + U + \ frac {\ partial S} {\ partial t} = 0.
Applications

HJE in a gravitational field

Using the energy–momentum relation in the form

g α β P α P β − ( m c) 2 = 0 {\displaystyle g^{\alpha \beta }P_{\alpha }P_{\beta }-(mc)^{2}=0}{\ displaystyle g ^ {\ alpha \ beta} P _ {\ alpha} P _ {\ beta} - (mc) ^ {2} = 0}

for a particle of rest mass m {\displaystyle m}m travelling in curved space, where g α β {\displaystyle g^{\alpha \beta }}г ^ {\ альфа \ бета} are the contravariant coordinates of the metric tensor (i.e., the inverse metric ) solved from the Einstein field equation s, and c {\displaystyle c}c is the speed of light. Setting the four-momentum P α {\displaystyle P_{\alpha }}P_ \ alpha equal to the four-gradient of the action S {\displaystyle S}S ,

P α = − ∂ S ∂ x α {\displaystyle P_{\alpha }=-{\frac {\partial S}{\partial x^{\alpha }}}}{\ displaystyle P _ {\ alpha} = - {\ frac {\ partial S} {\ partial x ^ {\ alpha}}}}

gives the Hamilton–Jacobi equation in the geometry determined by the metric g {\displaystyle g}g :

g α β ∂ S ∂ x α ∂ S ∂ x β − ( m c) 2 = 0, {\displaystyle g^{\alpha \beta }{\frac {\partial S}{\partial x^{\alpha }}}{\frac {\partial S}{\partial x^{\beta }}}-(mc)^{2}=0,}{\ displaystyle g ^ {\ alpha \ beta} {\ frac {\ partial S} {\ частичный x ^ {\ alpha}}} {\ frac {\ partial S} {\ partial x ^ {\ beta}}} - (mc) ^ {2} = 0,}

in other words, in a gravitational field.

HJE in electromagnetic fields

For a particle of rest mass m {\displaystyle m}m and electric charge e {\displaystyle e}e moving in electromagnetic field with four-potential A i = ( ϕ, A) {\displaystyle A_{i}=(\phi,\mathrm {A})}A_i = (\ phi, \ Alpha) in vacuum, the Hamilton–Jacobi equation in geometry determined by the metric tensor g i k = g i k {\di splaystyle g^{ik}=g_{ik}}g ^ {ik} = g_ {ik} has a form

g i k ( ∂ S ∂ x i + e c A i) ( ∂ S ∂ x k + e c A k) = m 2 c 2 {\displaystyle g^{ik}\left({\frac {\partial S}{\partial x^{i}}}+{\frac {e}{c}}A_{i}\right)\left({\frac {\partial S}{\partial x^{k}}}+{\frac {e}{c}}A_{k}\right)=m^{2}c^{2}}g ^ {ik} \ left (\ frac {\ partial S} {\ partial x ^ i} + \ frac {e} {c} A_i \ right) \ left (\ frac {\ partial S} {\ partial x ^ k} + \ frac {e} {c} A_k \ right) = m ^ 2 c ^ 2

and can be solved for the Hamilton principal action function S {\displaystyle S}S to obtain further solution for the particle trajectory and momentum:

x = − e c γ ∫ A z d ξ, {\displaystyle x=-{\frac {e}{c\gamma }}\int A_{z}\,d\xi,}{\ displaystyle x = - {\ frac {e} {c \ gamma }} \ int A_ {z} \, d \ xi,}
y = − e c γ ∫ A y d ξ, {\displaystyle y=-{\frac {e}{c\gamma }}\int A_{y}\,d\xi,}{\ displaystyle y = - {\ frac {e} {c \ gamma}} \ int A_ {y} \, d \ xi,}
z = − e 2 2 c 2 γ 2 ∫ ( A 2 − A 2 ¯) d ξ, {\displaystyle z=-{\frac {e^{2}}{2c^{2}\gamma ^{2}}}\int (\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}})\,d\xi,}{\ displaystyle z = - {\ frac {e ^ {2}} {2c ^ {2} \ gamma ^ {2}}} \ int (\ mathrm {A} ^ {2} - {\ overline {\ mathrm {A} ^ {2}}}) \, d \ xi,}
ξ = c t − e 2 2 γ 2 c 2 ∫ ( A 2 − A 2 ¯) d ξ, {\displaystyle \xi =ct-{\frac {e^{2}}{2\gamma ^{2}c^{2}}}\int (\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}})\,d\xi,}{\ displaystyle \ xi = ct - {\ frac {e ^ {2}} {2 \ gamma ^ {2} c ^ {2}}} \ int (\ mathrm {A} ^ {2} - { \ overline {\ mathrm {A} ^ {2}}}) \, d \ xi,}
p x = − e c A x {\displ aystyle p_{x}=-{\frac {e}{c}}A_{x}}p_x = - \ frac {e} {c} A_x , p y = − e c A y, {\displaystyle p_{y}=-{\frac {e}{c}}A_{y},}{\ displaystyle p_ {y} = - {\ frac {e} {c}} A_ {y},}
p z = e 2 2 γ c ( A 2 − A 2 ¯), {\displaystyle p_{z}={\frac {e^{2}}{2\gamma c}}(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}),}{\ displaystyle p_ {z} = {\ frac {e ^ {2}} {2 \ gamma c}} (\ mathrm {A} ^ {2} - {\ overline {\ mathrm {A} ^ {2}}}),}
E = c γ + e 2 2 γ c ( A 2 − A 2 ¯), {\displaystyle {\mathcal {E}}=c\gamma +{\frac {e^{2}}{2\gamma c}}(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}),}{\ displaystyle {\ mathcal {E}} = c \ gamma + {\ frac {e ^ {2}} {2 \ gamma c}} (\ mathrm {A} ^ {2} - {\ overline {\ mathrm {A} ^ {2 }}}),}

where ξ = c t − z {\displaystyle \xi =ct-z}\ xi = ct - z and γ 2 = m 2 c 2 + e 2 c 2 A ¯ 2 {\displaystyle \gamma ^{2}=m^{2}c^{2}+{\frac {e^{2}}{c^{2}}}{\overline {A}}^{2}}{\ displaystyle \ гамма ^ {2} = m ^ {2} c ^ {2} + {\ frac {e ^ {2}} {c ^ {2}}} {\ overline {A}} ^ {2}} with A ¯ {\displaystyle {\overline {\mathbf {A} }}}{\ overline {\ mathbf {A}}} the cycle average of the vector potential.

A circularly polarized wave

In the case of circular polarization,

E x = E 0 sin ⁡ ω ξ 1 {\displaystyle E_{x}=E_{0}\sin \omega \xi _{1}}E_x = E_0 \ sin \ omega \ xi_1 , E y = E 0 cos ⁡ ω ξ 1, {\displaystyle E_{y}=E_{0}\cos \omega \xi _{1},}{\ displaystyle E_ {y} = E_ {0} \ cos \ omega \ xi _ {1},}
A x = c E 0 ω cos ⁡ ω ξ 1 {\displaystyle A_{x}={\frac {cE_{0}}{\omega }}\cos \omega \xi _{1}}A_x = \ гидроразрыв {cE_0} {\ omega} \ cos \ omega \ xi_1 , A y = − c E 0 ω sin ⁡ ω ξ 1. {\displaystyle A_{y}=-{\frac {cE_{0}}{\omega }}\sin \omega \xi _{1}.}{\ displaystyle A_ {y} = - {\ frac {cE_ {0}} { \ omega}} \ sin \ omega \ xi _ {1}.}

Hence

x = - ec E 0 ω sin ⁡ ω ξ 1, {\ displaystyle x = - {\ frac {ecE_ {0}} {\ omega}} \ sin \ omega \ xi _ {1},}{\ displaystyle x = - {\ frac {ecE_ {0}} {\ omega}} \ sin \ omega \ xi _ {1},}
y = - ec E 0 ω cos ⁡ ω ξ 1, {\ displaystyle y = - {\ frac {ecE_ {0}} {\ omega}} \ cos \ omega \ xi _ {1},}{ \ displaystyle y = - {\ frac {ecE_ {0}} {\ omega}} \ cos \ omega \ xi _ {1},}
пикселей Знак равно - е E 0 ω соз ⁡ ω ξ 1, {\ displaystyle p_ {x} = - {\ frac {eE_ {0}} {\ omega}} \ cos \ omega \ xi _ {1},}{\ displaystyle p_ {x} = - {\ frac {eE_ {0}} {\ omega}} \ cos \ omega \ xi _ {1},}
py = e E 0 ω sin ⁡ ω ξ 1, {\ displaystyle p_ {y} = {\ frac {eE_ {0}} {\ omega}} \ sin \ omega \ xi _ {1},}{\ displaystyle p_ {y} = {\ frac {eE_ {0}} {\ omega}} \ sin \ omega \ xi _ {1},}

где ξ 1 = ξ / c {\ displaystyle \ xi _ {1} = \ xi / c}\ xi_1 = \ xi / c , подразумевая, что частица движется по круговой траектории с постоянным радиусом ec E 0 / γ ω 2 {\ displaystyle ecE_ {0} / \ gamma \ omega ^ {2}}ecE_0 / \ gamma \ omega ^ 2 и неизменное значение импульса e E 0 / ω 2 {\ displaystyle eE_ {0} / \ omega ^ {2}}eE_0 / \ omega ^ 2 , направленная вдоль вектора магнитного поля.

Монохроматическая линейно поляризованная плоская волна

Для плоской монохроматической линейно поляризованной волны с полем E {\ displaystyle E}E , направленным вдоль оси Y {\ displaystyle y}y

E y = E 0 cos ⁡ ω ξ 1, {\ displaystyle E_ {y} = E_ {0} \ cos \ omega \ xi _ {1},}{\ displaystyle E_ {y} = E_ {0} \ cos \ omega \ xi _ {1},}
A Y = - с E 0 ω грех ⁡ ω ξ 1, {\ displaystyle A_ {y} = - {\ frac {cE_ {0}} {\ omega}} \ sin \ omega \ xi _ {1},}{\ displaystyle A_ {y} = - {\ frac {cE_ {0}} {\ omega}} \ sin \ omega \ xi _ {1},}

отсюда

x = const, {\ displaystyle x = {\ text {const}},}{\ displaystyle x = { \ text {const}},}
y 0 = - ec E 0 γ ω 2, {\ displaystyle y_ {0} = - {\ frac { ecE_ {0}} {\ gamma \ omega ^ {2}}},}{\ displaystyle y_ {0} = - { \ frac {ecE_ {0}} {\ gamma \ omega ^ {2}}},}
y = y 0 cos ⁡ ω ξ 1 {\ displaystyle y = y_ {0} \ cos \ omega \ xi _ {1}}y = y_0 \ cos \ omega \ xi_1 , z знак равно C zy 0 грех ⁡ 2 ω ξ 1, {\ displaystyle z = C_ {z} y_ {0} \ sin 2 \ omega \ xi _ {1},}{\ displaystyle z = C_ {z} y_ {0} \ sin 2 \ omega \ xi _ {1},}
C z = e E 0 8 γ ω {\ displaystyle C_ {z} = {\ frac {eE_ {0}} {8 \ gamma \ omega}}}{\ displaystyle C_ {z} = {\ frac {eE_ {0}} {8 \ gamma \ omega }}} , γ 2 = m 2 c 2 + e 2 E 0 2 2 ω 2, { \ displaystyle \ gamma ^ {2} = m ^ {2} c ^ {2} + {\ frac {e ^ {2} E_ {0} ^ {2}} {2 \ omega ^ {2}}},}{\ displaystyle \ gamma ^ {2} = m ^ {2} c ^ {2} + {\ гидроразрыв {е ^ {2} E_ {0} ^ {2}} {2 \ omega ^ {2}}},}
px = 0, {\ displaystyle p_ {x} = 0,}{\ displaystyle p_ {x} = 0,}
py, 0 = e E 0 ω, {\ displaystyle p_ {y, 0} = {\ frac {eE _ {0}} {\ omega}},}{\ displaystyle p_ {y, 0} = {\ frac {eE_ {0}} {\ omega}},}
py = py, 0 грех ⁡ ω ξ 1, {\ displaystyle p_ {y} = p_ {y, 0} \ sin \ omega \ xi _ {1},}{\ displaystyle p_ {y} = p_ {y, 0} \ sin \ omega \ xi _ {1},}
pz = - 2 C zpy, 0 cos ⁡ 2 ω ξ 1 {\ displaystyle p_ {z} = - 2C_ {z} p_ {y, 0} \ cos 2 \ omega \ xi _ {1}}p_z = - 2C_z p_ {y, 0} \ cos 2 \ omega \ xi_1

подразумевает траекторию частицы в форме восьмерки с длинной осью, ориентированной вдоль вектора электрического поля E {\ displaystyle E}E .

Электромагнитная волна с соленоидальным магнитным полем

Для электромагнитной волны с аксиальным (соленоидальным) магнитным полем:

E = E ϕ = ω ρ 0 c B 0 cos ⁡ ω ξ 1, {\ displaystyle E = E _ {\ phi} = {\ frac {\ omega \ rho _ {0}} {c}} B_ {0} \ cos \ omega \ xi _ {1},}{\ displaystyle E = E _ {\ phi } = {\ frac {\ omega \ rho _ {0}} {c}} B_ {0} \ cos \ omega \ xi _ {1},}
A ϕ знак равно - ρ 0 В 0 грех ⁡ ω ξ 1 знак равно - L s π ρ 0 N s I 0 грех ⁡ ω ξ 1, {\ displaystyle A _ {\ phi} = - \ rho _ {0} B_ {0} \ sin \ omega \ xi _ {1} = - {\ frac {L_ {s}} {\ pi \ rho _ {0} N_ {s}}} I_ {0} \ sin \ omega \ xi _ {1}, }{\ displaystyle A _ {\ phi} = - \ rho _ {0} B_ {0} \ sin \ omega \ xi _ {1} = - {\ frac { L_ {s}} {\ pi \ rho _ {0} N_ {s}}} I_ {0} \ sin \ omega \ xi _ {1},}

следовательно,

x = константа, {\ displaystyle x = {\ text {constant}},}{\ displaystyle x = {\ text {constant}},}
y 0 = - e ρ 0 B 0 γ ω, {\ displaystyle y_ {0} = - {\ гидроразрыва {е \ rho _ {0} B_ {0}} {\ gamma \ omega}},}{\ displaystyle y_ {0} = - {\ frac {e \ rho _ {0} B_ {0}} {\ gamma \ omega}},}
y = y 0 cos ⁡ ω ξ 1, {\ displaystyle y = y_ {0} \ cos \ омега \ xi _ {1},}{\ displaystyle y = y_ {0} \ cos \ omega \ xi _ {1},}
z = C zy 0 sin ⁡ 2 ω ξ 1, {\ displaystyle z = C_ {z} y_ {0} \ sin 2 \ omega \ xi _ {1},}{\ displaystyle z = C_ {z} y_ {0} \ sin 2 \ omega \ xi _ {1},}
C z знак равно e ρ 0 B 0 8 c γ, {\ displaystyle C_ {z} = {\ frac {e \ rho _ {0} B_ {0}} {8c \ gamma}},}{\ displaystyle C_ {z} = {\ frac {e \ rho _ {0} B_ {0}} {8c \ gamma}},}
γ 2 знак равно м 2 с 2 + е 2 ρ 0 2 В 0 2 2 с 2, {\ displaystyle \ gamma ^ {2} = m ^ {2} c ^ {2} + {\ frac {e ^ {2} \ rho _ {0} ^ {2} B_ {0} ^ {2}} {2c ^ {2}}},}{\ displaystyle \ gamma ^ {2} = m ^ {2} c ^ {2} + {\ frac {e ^ {2} \ rho _ {0} ^ {2} B_ {0} ^ {2}} {2c ^ {2}}},}
px = 0, {\ displaystyle p_ {x} = 0,}{\ displaystyle p_ {x} = 0,}
py, 0 = e ρ 0 B 0 c, {\ displaystyle p_ {y, 0} = {\ frac {e \ rho _ {0} B_ {0}} {c}},}{\ displaystyle p_ {y, 0} = {\ frac {e \ rho _ {0} B_ {0}} {c}},}
py = py, 0 грех ⁡ ω ξ 1, {\ displaystyle p_ {y} = p_ {y, 0} \ sin \ omega \ xi _ {1},}{\ displaystyle p_ {y} = p_ {y, 0} \ sin \ omega \ xi _ {1},}
pz = - 2 C zpy, 0 cos ⁡ 2 ω ξ 1, {\ displaystyle p_ {z} = - 2C_ {z} p_ {y, 0} \ cos 2 \ omega \ xi _ {1},}{\ displaystyle p_ {z} = - 2C_ {z} p_ {y, 0} \ cos 2 \ omega \ xi _ {1},}

где B 0 {\ displaystyle B_ {0}}B_{0}- величина магнитного поля в соленоиде с эффективным радиусом ρ 0 {\ displaystyle \ rho _ {0}}\ rho _ {0} , индуктивность L s {\ displaystyle L_ {s}}L_s , количество витков N s {\ displaystyle N_ {s}}N_s , и величина электрического тока I 0 {\ displaystyle I_ {0}}I_ {0} через обмотки соленоида. Движение частицы происходит по траектории восьмерки в плоскости yz {\ displaystyle yz}yz , установленной перпендикулярно оси соленоида с произвольным азимутальным углом φ {\ displaystyle \ varphi}\ varphi из-за осевой симметрии соленоидального магнитного поля.

См. Также
  • значок Математический портал
  • значок Физический портал
Ссылки
Дополнительная литература
Последняя правка сделана 2021-05-22 11:54:31
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте