Мультипольное разложение — это математический ряд, представляющий функцию , которая зависит от углов — обычно двух углов, используемых в сферической системе координат (полярный и азимутальный углы) для трехмерного евклидова пространства , . Подобно рядам Тейлора , мультипольные разложения полезны, поскольку часто для обеспечения хорошего приближения исходной функции требуются только первые несколько членов. Раскрываемая функция может быть действительной или комплексной и определяется либо на , либо, реже, на для некоторых других .
Мультипольные разложения часто используются при изучении электромагнитных и гравитационных полей , где поля в удаленных точках задаются в терминах источников в малой области. Мультипольное разложение с углами часто сочетается с разложением по радиусу . Такое сочетание дает разложение, описывающее функцию во всем трехмерном пространстве. [1]
Мультипольное расширение выражается как сумма членов с постепенно усложняющимися угловыми характеристиками ( моментами ). Первый член (нулевого порядка) называется монопольным моментом, второй член (первого порядка) называется дипольным моментом, третий (второго порядка) — квадрупольным моментом, четвертый (третьего порядка) — октупольным моментом и т. д. Учитывая ограничение греческих числовых префиксов , члены более высокого порядка обычно именуются путем добавления «-полюс» к числу полюсов — например, 32-полюсный (редко дотриаконтаполь или триаконтадиполь) и 64-полюсный (редко тетрагексаконтаполь или гексаконтатетраполь). [2] [3] [4] Мультипольный момент обычно включает в себя мощности (или обратные мощности) расстояния до начала координат, а также некоторую угловую зависимость.
В принципе, мультипольное разложение дает точное описание потенциала и, как правило, сходится при двух условиях: (1) если источники (например, заряды) локализованы близко к началу координат, а точка, в которой наблюдается потенциал, находится далеко от начала координат; или (2) наоборот, т. е. если источники расположены далеко от начала координат, а потенциал наблюдается близко к началу координат. В первом (более распространенном) случае коэффициенты разложения в ряд называются внешними мультипольными моментами или просто мультипольными моментами , тогда как во втором случае они называются внутренними мультипольными моментами .
Чаще всего ряд записывается как сумма сферических гармоник . Таким образом, мы могли бы записать функцию как сумму , где — стандартные сферические гармоники, а — постоянные коэффициенты, зависящие от функции. Член представляет монополь; представляет диполь; и так далее. Эквивалентно ряд также часто записывается [5] как , где представляют компоненты единичного вектора в направлении, заданном углами и , а индексы неявно суммируются . Здесь член — монополь; — набор из трех чисел, представляющих диполь; и так далее.
В приведенных выше разложениях коэффициенты могут быть действительными или комплексными . Однако, если функция, выраженная как многополюсное разложение, действительна, коэффициенты должны удовлетворять определенным свойствам. В сферическом гармоническом разложении мы должны иметь В многовекторном разложении каждый коэффициент должен быть действительным:
Хотя разложения скалярных функций являются наиболее распространенным применением мультипольных разложений, их также можно обобщить для описания тензоров произвольного ранга. [6] Это находит применение в мультипольных разложениях векторного потенциала в электромагнетизме или метрических возмущениях при описании гравитационных волн .
Для описания функций трех измерений, вдали от начала координат, коэффициенты мультипольного разложения можно записать как функции расстояния до начала координат, —чаще всего, как ряд Лорана по степеням . Например, для описания электромагнитного потенциала, , от источника в небольшой области вблизи начала координат, коэффициенты можно записать как:
Мультипольные разложения широко используются в задачах, связанных с гравитационными полями систем масс , электрическими и магнитными полями зарядов и распределений токов, а также распространением электромагнитных волн . Классическим примером является вычисление внешних мультипольных моментов атомных ядер из их энергий взаимодействия с внутренними мультиполями электронных орбиталей. Мультипольные моменты ядер сообщают о распределении зарядов внутри ядра и, таким образом, о форме ядра. Усечение мультипольного разложения до его первого ненулевого члена часто полезно для теоретических расчетов.
Мультипольные разложения также полезны в численном моделировании и составляют основу быстрого мультипольного метода Грингарда и Рохлина , общей методики эффективного вычисления энергий и сил в системах взаимодействующих частиц . Основная идея заключается в разложении частиц на группы; частицы внутри группы взаимодействуют нормально (т. е. полным потенциалом), тогда как энергии и силы между группами частиц вычисляются из их мультипольных моментов. Эффективность быстрого мультипольного метода в целом аналогична эффективности суммирования Эвальда , но выше, если частицы сгруппированы, т. е. система имеет большие флуктуации плотности.
Рассмотрим дискретное распределение зарядов, состоящее из N точечных зарядов q i с векторами положения r i . Мы предполагаем, что заряды сгруппированы вокруг начала координат, так что для всех i : r i < r max , где r max имеет некоторое конечное значение. Потенциал V ( R ) , обусловленный распределением зарядов, в точке R вне распределения зарядов, т. е. | R | > r max , может быть разложен по степеням 1/ R . В литературе можно найти два способа сделать это разложение: первый - это ряд Тейлора в декартовых координатах x , y , и z , тогда как второй - в терминах сферических гармоник , которые зависят от сферических полярных координат . Декартов подход имеет то преимущество, что не требуется никаких предварительных знаний о функциях Лежандра, сферических гармониках и т. д. Его недостаток в том, что выводы довольно громоздки (фактически большая их часть представляет собой неявный повторный вывод разложения Лежандра 1 / | r − R | , которое было сделано Лежандром раз и навсегда в 1780-х годах). Также трудно дать замкнутое выражение для общего члена многополюсного разложения — обычно приводятся только первые несколько членов с последующим многоточием.
Предположим для удобства v ( r ) = v (− r ) . Разложение Тейлора v ( r − R ) вокруг начала координат r = 0 можно записать как с коэффициентами Тейлора Если v ( r − R ) удовлетворяет уравнению Лапласа , то согласно приведенному выше разложению имеем и разложение можно переписать в терминах компонентов бесследового декартова тензора второго ранга : где δ αβ — символ Кронекера , а r 2 ≡ | r | 2 . Удаление следа является обычным делом, поскольку оно выводит инвариантный относительно вращения r 2 из тензора второго ранга.
Рассмотрим теперь следующую форму v ( r − R ) : Тогда путем прямого дифференцирования следует, что Определим монополь, диполь и (бесследовый) квадруполь как, соответственно, и мы получаем в итоге первые несколько членов мультипольного разложения полного потенциала, который является суммой кулоновских потенциалов отдельных зарядов: [7] : 137–138
Это расширение потенциала дискретного распределения заряда очень похоже на то, что в реальных твердых гармониках, приведенных ниже. Главное отличие в том, что настоящее расширение находится в терминах линейно зависимых величин, для
Примечание: Если распределение зарядов состоит из двух зарядов противоположного знака, находящихся на бесконечно малом расстоянии d друг от друга, так что d / R ≫ ( d / R ) 2 , легко показать, что доминирующим членом в разложении является электрическое дипольное потенциальное поле .
Потенциал V ( R ) в точке R вне распределения заряда, т.е. | R | > r max , может быть разложен с помощью разложения Лапласа : где — нерегулярная сплошная гармоника (определяемая ниже как сферическая гармоническая функция, деленная на ), а — регулярная сплошная гармоника (сферическая гармоника, умноженная на r ℓ ). Мы определяем сферический мультипольный момент распределения заряда следующим образом. Обратите внимание, что мультипольный момент определяется исключительно распределением заряда (положениями и величинами N зарядов).
Сферическая гармоника зависит от единичного вектора . (Единичный вектор определяется двумя сферическими полярными углами.) Таким образом, по определению, нерегулярные сплошные гармоники можно записать в виде так, что мультипольное разложение поля V ( R ) в точке R вне распределения заряда задается выражением
Это расширение является полностью общим, поскольку оно дает замкнутую форму для всех членов, а не только для первых нескольких. Оно показывает, что сферические мультипольные моменты появляются как коэффициенты в 1/ R -расширении потенциала.
Интересно рассмотреть первые несколько членов в действительной форме, которые являются единственными членами, обычно встречающимися в учебниках для студентов. Поскольку слагаемое суммы m инвариантно относительно унитарного преобразования обоих множителей одновременно и поскольку преобразование комплексных сферических гармоник в действительную форму осуществляется унитарным преобразованием , мы можем просто заменить действительные нерегулярные сплошные гармоники и действительные мультипольные моменты. Член ℓ = 0 становится Это фактически снова закон Кулона . Для члена ℓ = 1 мы вводим Тогда Этот член идентичен найденному в декартовой форме.
Чтобы записать член ℓ = 2 , нам нужно ввести сокращенные обозначения для пяти действительных компонентов квадрупольного момента и действительных сферических гармоник. Обозначения такого типа можно найти в литературе. Очевидно, что действительные обозначения очень скоро становятся неудобными, что демонстрирует полезность комплексных обозначений.
Рассмотрим два набора точечных зарядов, один набор { q i }, сгруппированный вокруг точки A , и один набор { q j }, сгруппированный вокруг точки B. Рассмотрим, например, две молекулы и вспомним, что молекула по определению состоит из электронов (отрицательных точечных зарядов) и ядер (положительных точечных зарядов). Полная энергия электростатического взаимодействия U AB между двумя распределениями равна Эту энергию можно разложить в степенной ряд по обратному расстоянию между A и B. Это расширение известно как мультипольное расширение U AB .
Чтобы вывести это мультипольное разложение, мы записываем r XY = r Y − r X , что является вектором, направленным от X к Y . Обратите внимание, что
Мы предполагаем, что два распределения не перекрываются:
При этом условии мы можем применить разложение Лапласа в следующей форме
, где и являются нерегулярными и регулярными сплошными гармониками соответственно. Перевод регулярной сплошной гармоники дает конечное разложение,
где величина между заостренными скобками является коэффициентом Клебша–Гордана . Далее мы использовали
Использование определения сферических мультиполей Q м
ли покрытие диапазонов суммирования в несколько ином порядке (что допускается только для бесконечного диапазона L ) дает в конечном итоге
Это мультипольное разложение энергии взаимодействия двух неперекрывающихся распределений зарядов, находящихся на расстоянии R AB друг от друга. Поскольку это разложение явно находится в степенях 1 / R AB . Функция Y m l является нормированной сферической гармоникой .
Все атомы и молекулы (кроме атомов в S -состоянии ) имеют один или несколько неисчезающих постоянных мультипольных моментов. В литературе можно найти различные определения, но следующее определение в сферической форме имеет то преимущество, что оно содержится в одном общем уравнении. Поскольку оно находится в сложной форме, оно имеет еще одно преимущество, что его легче манипулировать в расчетах, чем его реальный аналог.
Рассмотрим молекулу, состоящую из N частиц (электронов и ядер) с зарядами eZ i . (У электронов значение Z равно −1, тогда как у ядер это атомный номер ). Частица i имеет сферические полярные координаты r i , θ i , и φ i и декартовы координаты x i , y i , и z i . (Комплексный) электростатический мультипольный оператор равен , где — регулярная сплошная гармоническая функция в нормализации Рака (также известной как полунормализация Шмидта). Если молекула имеет полную нормализованную волновую функцию Ψ (зависящую от координат электронов и ядер), то мультипольный момент порядка молекулы задается математическим ожиданием (ожидаемым) значением : Если молекула имеет определенную точечную групповую симметрию , то это отражается в волновой функции: Ψ преобразуется в соответствии с определенным неприводимым представлением λ группы ( «Ψ имеет тип симметрии λ»). Это имеет следствием то, что правила отбора справедливы для ожидаемого значения оператора мультиполя, или, другими словами, что ожидаемое значение может исчезать из-за симметрии. Хорошо известным примером этого является тот факт, что молекулы с центром инверсии не несут диполь (ожидаемые значения исчезают при m = −1, 0, 1) . Для молекулы без симметрии никакие правила отбора не действуют, и такая молекула будет иметь неисчезающие мультиполи любого порядка (она будет нести диполь и одновременно квадруполь, октуполь, гексадекаполь и т. д.).
Низшие явные формы регулярных твердых гармоник (с фазой Кондона-Шортли ) дают: (полный заряд молекулы). (Комплексные) дипольные компоненты:
Обратите внимание, что с помощью простой линейной комбинации можно преобразовать сложные мультипольные операторы в действительные. Действительные мультипольные операторы имеют тип косинуса или синуса . Вот несколько самых низких из них:
Определение комплексного молекулярного мультипольного момента, данное выше, является комплексно сопряженным определением, данным в этой статье , которое следует определению стандартного учебника по классической электродинамике Джексона, [7] : 137 за исключением нормализации. Более того, в классическом определении Джексона эквивалент N -частичного квантово-механического ожидаемого значения является интегралом по одночастичному распределению заряда. Помните, что в случае одночастичной квантово-механической системы ожидаемое значение есть не что иное, как интеграл по распределению заряда (модуль волновой функции в квадрате), так что определение этой статьи является квантово-механическим N -частичным обобщением определения Джексона.
Определение в этой статье согласуется, среди прочего, с определением Фано и Рака [8], а также Бринка и Сэтчлера [9] .
Существует много типов мультипольных моментов, поскольку существует много типов потенциалов и много способов аппроксимации потенциала путем разложения в ряд , в зависимости от координат и симметрии распределения заряда. Наиболее распространенные разложения включают:
Примерами 1/ R -потенциалов являются электрический потенциал , магнитный потенциал и гравитационный потенциал точечных источников. Примером ln R- потенциала является электрический потенциал бесконечного линейного заряда.
Мультипольные моменты в математике и математической физике образуют ортогональный базис для разложения функции, основанный на реакции поля на точечные источники, которые расположены бесконечно близко друг к другу. Их можно рассматривать как расположенные в различных геометрических формах, или, в смысле теории распределения , как производные по направлению .
Мультипольные разложения связаны с базовой вращательной симметрией физических законов и связанных с ними дифференциальных уравнений . Даже если исходные члены (такие как массы, заряды или токи) могут быть несимметричными, их можно разложить в терминах неприводимых представлений группы вращательной симметрии , что приводит к сферическим гармоникам и связанным с ними наборам ортогональных функций. Используется метод разделения переменных для извлечения соответствующих решений для радиальных зависимостей.
На практике многие поля можно хорошо аппроксимировать конечным числом мультипольных моментов (хотя для точной реконструкции поля может потребоваться бесконечное число). Типичным применением является аппроксимация поля локализованного распределения заряда его монопольными и дипольными членами. Задачи, решенные один раз для заданного порядка мультипольного момента, можно линейно объединить для создания окончательного приближенного решения для заданного источника.