Генерация второй гармоники: нелинейная оптика и SHG

Генерация второй гармоники - самое наглядное явление нелинейной оптики: лазер на длине волны $\lambda_\omega$ проходит через кристалл и часть света выходит уже на $\lambda_{2\omega} = \lambda_\omega / 2$. Так из инфракрасного Nd:YAG (1064 нм) получают зелёную линию 532 нм, а из Ti:Sapphire (800 нм) - фиолетовую 400 нм. За эффект отвечает квадратичная восприимчивость среды $\chi^{(2)}$, и работает он только при выполнении двух условий: у кристалла нет центра инверсии и для двух волн выполнен фазовый синхронизм $\Delta k = 0$.

Откуда берётся вторая гармоника: нелинейная поляризация

В обычной (линейной) оптике поляризация среды пропорциональна полю: $P = \varepsilon_0 \chi^{(1)} E$. При больших интенсивностях разложение продолжается:

$$P = \varepsilon_0 \left( \chi^{(1)} E + \chi^{(2)} E^2 + \chi^{(3)} E^3 + \ldots \right).$$

Квадратичный член $P^{(2)} = \varepsilon_0 \chi^{(2)} E^2$ при гармонической накачке $E = E_0 \cos(\omega t)$ даёт

$$P^{(2)} = \tfrac{1}{2}\varepsilon_0 \chi^{(2)} E_0^2 \big[ 1 + \cos(2\omega t) \big].$$

Первое слагаемое - постоянная поляризация (оптическое выпрямление), второе колеблется на удвоенной частоте $2\omega$ и излучает вторую гармонику. Это и есть SHG.

Почему нужен кристалл без центра инверсии

Если среда центросимметрична (стекло, газ, жидкость с хаотически ориентированными молекулами), то операция $E \to -E$ должна давать $P \to -P$. Но $E^2$ при смене знака поля не меняется, значит $P^{(2)} = 0$ - иначе нарушается симметрия. Поэтому в стекле SHG в объёме запрещена и наблюдается только на поверхности (где центр инверсии нарушен) или в нелинейности третьего порядка $\chi^{(3)}$.

Все классические SHG-кристаллы - нецентросимметричные: $\beta$-BaB$_2$O$_4$ (BBO), KH$_2$PO$_4$ (KDP), LiNbO$_3$, LiB$_3$O$_5$ (LBO), KTP. У каждого свой тензор $d_{ijk} = \tfrac{1}{2}\chi^{(2)}_{ijk}$ - эффективная нелинейность зависит от того, как поляризованы волны накачки и гармоники относительно кристаллографических осей.

Фазовый синхронизм: главное препятствие

Накачка и гармоника бегут по кристаллу с разными показателями преломления, потому что $n(\omega) \neq n(2\omega)$ (нормальная дисперсия). Волновые векторы:

$$k_\omega = \frac{\omega}{c} n(\omega), \qquad k_{2\omega} = \frac{2\omega}{c} n(2\omega).$$

Гармоники, рождённые в разных точках кристалла, складываются в фазе только если $\Delta k = k_{2\omega} - 2 k_\omega = 0$, то есть $n(2\omega) = n(\omega)$. Без специальных мер $\Delta k \neq 0$, и поток гармоники сначала растёт, потом убывает с периодом когерентной длины $L_c = \pi / |\Delta k|$ (для типичной видимой накачки $L_c \sim 10$ мкм).

Формула интенсивности SHG

В режиме слабого истощения накачки интенсивность гармоники на выходе кристалла длиной $L$:

$$I_{2\omega} \propto |\chi^{(2)}|^2\, I_\omega^2\, L^2\, \text{sinc}^2\!\left(\frac{\Delta k L}{2}\right).$$

Три ключевых вывода:

• Квадратично растёт с накачкой $I_\omega$ - характерная подпись нелинейного процесса.
• Квадратично растёт с длиной $L$, но только при идеальном синхронизме.
• Множитель $\text{sinc}^2(\Delta k L / 2)$ режет всё, что не попадает в узкое окно по $\Delta k$.

Из последнего получается угловая, спектральная и температурная ширина синхронизма - диапазоны, в которых ещё можно эффективно генерировать 2ω.

Типы синхронизма I и II

Чтобы скомпенсировать дисперсию, используют двулучепреломление в одноосном кристалле: показатель необыкновенной волны зависит от угла $\theta$ между волновым вектором и оптической осью.

• Тип I (oo→e): обе волны накачки - обыкновенные, гармоника - необыкновенная. Подбирают такой угол $\theta$, при котором $n_e(2\omega, \theta) = n_o(\omega)$. Это самый частый сценарий для BBO/KDP.
• Тип II (oe→e): одна волна накачки обыкновенная, другая - необыкновенная, гармоника - необыкновенная. У типа II шире спектральная и угловая ширина, удобнее для широкополосных импульсов.

В обоих случаях есть угловой walk-off: необыкновенная волна уходит в сторону под небольшим углом $\rho$ к волновому вектору. Это ограничивает полезную длину кристалла даже при идеальном фазовом согласовании.

Кристаллы: BBO, KDP, LBO, LiNbO₃

Каждый кристалл оптимален в своём окне.

• BBO - самый универсальный нелинейный кристалл для видимого диапазона: $d_{\text{eff}} \approx 2.0$ пм/В, прозрачен 190–3500 нм, высокая лучевая прочность. Стандарт для второй и третьей гармоник от Ti:Sapphire и Nd:YAG.
• KDP - исторический рабочий лоток, $d_{\text{eff}} \approx 0.4$ пм/В, зато выращивается метровыми блоками. Используется на установках лазерного термоядерного синтеза (NIF, ЛФИ).
• LBO - среднее звено: $d_{\text{eff}} \approx 0.8$ пм/В, очень малый walk-off, удобен для непрерывной и пико-секундной накачки.
• LiNbO$_3$ - большой коэффициент $d_{33} \approx 30$ пм/В, но в чистом виде нет фазового синхронизма для типичной видимой накачки - используется в форме PPLN (см. ниже).

Квази-фазовый синхронизм и PPLN

Когда обычный фазовый синхронизм недоступен (например, для самого большого тензорного коэффициента $d_{33}$), применяют квази-фазовый синхронизм (QPM): каждые $\Lambda / 2 = L_c$ направление спонтанной поляризации в кристалле инвертируют - получается периодически поляризованный кристалл, обычно периодически поляризованный ниобат лития (PPLN).

Условие QPM:

$$\Delta k_{\text{QPM}} = k_{2\omega} - 2 k_\omega - \frac{2\pi}{\Lambda} = 0,$$

где $\Lambda$ - период доменной структуры (обычно 5–30 мкм). Период подбирается под выбранную длину волны и температуру; перестройку гармоники делают либо температурой, либо набором кристаллов с разными $\Lambda$ на одной подложке (fan-out PPLN). QPM даёт эффективность в 10–100 раз выше угловой схемы при тех же длинах.

Где встречается SHG

• Зелёные лазеры 532 нм в указках, накачке Ti:Sapphire, медицине, голографии - все они работают через SHG в KTP/LBO/PPLN из 1064 нм.
• УФ-источники: третья и четвёртая гармоники Nd:YAG (355 и 266 нм) - это каскад SHG + сумма частот в BBO.
• Сверхбыстрая диагностика: автокорреляторы фемтосекундных импульсов измеряют длительность через SHG в тонком BBO.
• SHG-микроскопия: коллаген, миозин, микротрубочки имеют упорядоченную нецентросимметричную структуру - они светятся на 2ω без меток и контраста, в отличие от изотропной воды.
• Метрология длин волн: удвоение и сложение частот - основа оптических часов и привязки фемтосекундных гребёнок.

Типичные задачи в курсе нелинейной оптики

1. Найти длину когерентности $L_c = \pi / |\Delta k|$ для заданной пары $\lambda_\omega, \lambda_{2\omega}$ и показателей $n(\omega), n(2\omega)$.
2. Найти угол синхронизма $\theta_s$ типа I в одноосном кристалле из уравнения $n_e(2\omega, \theta_s) = n_o(\omega)$.
3. Оценить эффективность $\eta = I_{2\omega} / I_\omega$ в режиме слабой накачки и сравнить с экспериментальным значением.
4. Подобрать период $\Lambda$ для PPLN под заданную длину волны и температуру.
5. Оценить угловую и спектральную ширину синхронизма из условия $\Delta k L / 2 = \pi / 2$.

Частые ошибки

• Путают $\chi^{(2)}$ и $d_{\text{eff}}$: $d_{ijk} = \tfrac{1}{2}\chi^{(2)}_{ijk}$, и в формулах фигурирует свёртка тензора с поляризациями - это и есть эффективная нелинейность для выбранной геометрии.
• Забывают про walk-off: даже при $\Delta k = 0$ полезная длина кристалла ограничена расходимостью пучка $w_0 / \rho$.
• Считают, что SHG возможна в любом кристалле: в центросимметричных (CaF$_2$, NaCl, KBr) объёмный $\chi^{(2)} = 0$, гармоника может появиться только на интерфейсе.
• Игнорируют температуру: показатели преломления зависят от $T$, поэтому угол или период QPM подбирают вместе с термостабилизацией ($\pm 0.1$ °C для узкополосной 1 нм SHG).
• Подают слишком плотную накачку: при $I_\omega > 1$ ГВт/см² начинается истощение накачки, обратное преобразование и оптический пробой; формула $I_{2\omega} \propto I_\omega^2$ перестаёт работать.

FAQ

В чём разница между SHG и комбинационным рассеянием? SHG - параметрический процесс из $\chi^{(2)}$: энергия двух фотонов $\omega$ переходит в один фотон $2\omega$, среда возвращается в исходное состояние, фазовое согласование обязательно. Комбинационное рассеяние - неупругий процесс с участием колебаний молекул, частоты $\omega \pm \Omega_{\text{vib}}$, фазовое согласование не нужно, но сечение на 6–8 порядков меньше.

Можно ли получить SHG в стекле? В объёме оптического стекла нельзя - есть центр инверсии. На поверхности или интерфейсе SHG регистрируется и используется в поверхностной нелинейной спектроскопии (SFG). В оптоволокне можно индуцировать постоянную поляризацию термополингом - тогда стекло «учится» генерировать SHG.

Почему зелёные лазерные указки именно 532 нм, а не 510 или 550? Потому что доступный и дешёвый диодно-накачиваемый Nd:YVO$_4$ или Nd:YAG генерирует ровно 1064 нм, а $1064 / 2 = 532$ нм. Никакого «фундаментального» 532 нм - это просто вторая гармоника от инфракрасной линии неодима.

Коротко

Генерация второй гармоники описывается квадратичным членом нелинейной поляризации $P^{(2)} = \varepsilon_0 \chi^{(2)} E^2$ и работает только в кристаллах без центра инверсии. Чтобы гармоники, рождённые в разных точках, складывались конструктивно, нужен фазовый синхронизм $\Delta k = k_{2\omega} - 2 k_\omega = 0$ - его реализуют через двулучепреломление (тип I и II) в BBO, KDP, LBO или через периодическую инверсию доменов (QPM) в PPLN. Эффективность растёт как $I_\omega L^2$ и режется фактором $\text{sinc}^2(\Delta k L / 2)$ - на этом строятся и оценка выхода 2ω, и анализ угловой/спектральной/температурной ширины синхронизма.