Комплексный анализ в физике: теорема о вычетах

Комплексный анализ в физике: теорема о вычетах

Магия комплексной плоскости

Многие интегралы по вещественной оси вычисляются изящно через замыкание контура в комплексной плоскости. Идея: замкнутый интеграл по контуру C от аналитической функции = 0 (теорема Коши). Но если внутри C есть особые точки (полюсы), интеграл равен сумме вычетов, умноженной на 2πi. Выбирая умный контур, мы «собираем» нужный вычет и получаем реальный интеграл.

В физике этот метод встречается повсюду: вычисление пропагаторов в квантовой теории поля, дисперсионные соотношения в оптике, функции Грина в статмеханике.

Теорема Коши и вычеты

Интеграл Коши: Для аналитической f в односвязной области D и контура C ⊂ D:

∮_C f(z) dz = 0 (теорема Коши)

Но если f имеет полюс при z₀ внутри C:

∮C f(z) dz = 2πi Res{z=z₀} f

Вычет — коэффициент при (z − z₀)^{−1} в разложении Лорана:

f(z) = ... + a₋₂/(z−z₀)² + a₋₁/(z−z₀) + a₀ + a₁(z−z₀) + ...

Res = a₋₁.

Для простого полюса: Res_{z=z₀} f(z) = lim_{z→z₀} (z−z₀)f(z). Для дроби f = g/h при h(z₀)=0, g(z₀)≠0: Res = g(z₀)/h'(z₀).

Теорема о вычетах: ∮C f(z) dz = 2πi Σ{полюсы внутри C} Res f

Вычисление реальных интегралов

Интегралы ∫_{−∞}^{+∞} f(x) dx: Замыкаем контур полуокружностью в верхней (или нижней) полуплоскости.

Лемма Жордана: для e^{iaz} f(z) при a > 0 и |f(z)| → 0 на полуокружности — интеграл по дуге → 0.

Результат: ∫{−∞}^{+∞} f(x) dx = 2πi × Σ{Im(z₀)>0} Res f(z₀)

Числовой пример 1: I = ∫_{−∞}^{+∞} dx/(x² + 1)

Полюсы f(z) = 1/(z²+1) = 1/((z+i)(z−i)) при z = ±i. В верхней полуплоскости: z = i.

Res_{z=i} = lim_{z→i} (z−i)/(z²+1) = 1/(2i)

I = 2πi × 1/(2i) = π ✓ (проверка: ∫dx/(x²+1) = arctan(x)|_{−∞}^{+∞} = π)

Числовой пример 2: I = ∫_{−∞}^{+∞} cos(ax)/(x² + b²) dx, a, b > 0

Записываем как Re[∫ e^{iax}/(x²+b²) dx]. Полюс в верхней полуплоскости: z = ib.

Res_{z=ib} e^{iaz}/(z²+b²) = e^{−ab}/(2ib)

I = Re[2πi × e^{−ab}/(2ib)] = Re[πe^{−ab}/b] = πe^{−ab}/b

Дисперсионные соотношения Крамерса–Кронига

Физическое обоснование: функция отклика χ(ω) системы аналитична в верхней полуплоскости — это следствие причинности (отклик после воздействия, не до). Применяя теорему Коши к χ(ω)/(ω' − ω):

Re χ(ω) = (2/π) P.V. ∫₀^∞ ω' Im χ(ω')/(ω'² − ω²) dω'

Im χ(ω) = −(2ω/π) P.V. ∫₀^∞ Re χ(ω')/(ω'² − ω²) dω'

P.V. — главное значение по Коши (интеграл в смысле пределов при обходе полюса).

Физические следствия: Реальная часть (дисперсия = отклонение волн) определяется мнимой (поглощение) и наоборот. Это фундаментальное ограничение: материал, не поглощающий волны ни на какой частоте, не изменяет и их скорость. Дисперсионные соотношения используются в оптике, акустике, ядерном рассеянии, для проверки причинности моделей.

Функции Грина через вычеты

Функция Грина для осциллятора с затуханием: G(ω) = 1/(m(ω₀² − ω² − iγω)).

Полюсы: ω_± = −iγ/2 ± √(ω₀² − γ²/4). При γ < 2ω₀: два полюса с Im(ω_±) < 0 (в нижней полуплоскости). Причинная G(t) = обратное преобразование Фурье, замыкание контура вниз при t > 0:

G(t) = θ(t) × (1/mω_d) sin(ω_d t) e^{−γt/2}, ω_d = √(ω₀² − γ²/4)

Физически: это и есть затухающий осциллятор, нарастающий от нуля при t = 0.

Приложения в квантовой теории поля

Пропагатор Фейнмана для свободного скалярного поля: G(k) = i/(k² − m² + iε). Полюс при k² = m² — «mass shell» (частицы на массовой поверхности). Сдвиг +iε (рецепт Фейнмана) выбирает причинный (ретардированный) пропагатор. Все диаграммы Фейнмана — это продукты таких пропагаторов, интегрированных по импульсам через теорему о вычетах.

Комплексный анализ в сигнальной обработке и квантовой механике

Комплексный анализ и дисперсионные соотношения — инструменты обработки сигналов и квантовой теории. В теории систем передаточная функция H(ω) (импеданс, фильтр, усилитель) аналитически продолжается в верхнюю полуплоскость: причинность системы (отклик идёт после воздействия) эквивалентна тому, что полюсы H(ω) лежат в нижней полуплоскости. Это и есть условие устойчивости фильтра. Дисперсионные соотношения Крамерса–Кронига связывают вещественную и мнимую части диэлектрической проницаемости ε(ω): знание поглощения на всех частотах (Im ε) полностью определяет дисперсию (Re ε). Это используется в спектроскопии материалов: из измеренного поглощения восстанавливается полный показатель преломления. В ядерной физике дисперсионные соотношения для амплитуды рассеяния связывают упругое и неупругое рассеяние. В квантовой хромодинамике правила сумм из дисперсионных соотношений позволяют определять структурные функции протона из данных рассеяния. Метод перевала (метод наискорейшего спуска) используется для асимптотики собственных значений и спектральных коэффициентов — в квантовой механике для вычисления ВКБ-амплитуд прохождения через барьер при туннелировании.

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

Задание: (а) Вычислите ∫₀^∞ sin(x)/x dx через контур с вырезом вокруг нуля. (б) Для G(ω) осциллятора с γ = ω₀/10: найдите полюсы в комплексной плоскости. Постройте причинную G(t). При каком t максимум? (в) Проверьте дисперсионное соотношение для простой модели χ(ω) = ω₀/(ω₀² − ω² − iεω) — вычислите Re χ и Im χ и проверьте формулу Крамерса–Кронига численно.