Поляризационная модовая дисперсия: DGD и закон корня

Идеальное одномодовое волокно круглое и изотропное, и обе поляризации света бегут по нему с одной скоростью. Реальное - слегка эллиптично, имеет внутренние напряжения, изгибы и неоднородности, поэтому две взаимно перпендикулярные поляризации видят чуть разный показатель преломления и приходят на выход с задержкой друг относительно друга. Это и есть поляризационная модовая дисперсия (PMD): импульс «расщепляется» по поляризации и размывается, ограничивая скорость и дальность волоконных линий. Ниже - откуда берётся дифференциальная групповая задержка, почему она накапливается как корень из длины, и как это считать на типовых задачах. Если хочется сразу прикинуть DGD своей линии, воспользуйтесь калькулятором сразу под этим абзацем.

Двулучепреломление: откуда берутся две оси

Корень PMD - двулучепреломление (birefringence) волокна. Из-за неидеально круглой сердцевины и анизотропных механических напряжений показатель преломления зависит от направления поляризации. В сечении волокна выделяются две ортогональные оси: «быстрая» с меньшим показателем $n_x$ и «медленная» с большим $n_y$. Мера различия - групповое двулучепреломление $\Delta n = n_y - n_x$.

Рассчитать для своей задачи

Свет, поляризованный вдоль быстрой оси, обгоняет свет вдоль медленной. Эти две ортогональные поляризации называют главными состояниями поляризации (Principal States of Polarization, PSP) - на коротком отрезке они совпадают с осями двулучепреломления. Произвольно поляризованный входной импульс раскладывается по этим двум осям, и каждая компонента бежит со своей групповой скоростью.

Дифференциальная групповая задержка (DGD)

Ключевая величина PMD - дифференциальная групповая задержка $\Delta\tau$ (Differential Group Delay, DGD): разность времён прихода двух поляризационных компонент. Для короткого однородного участка волокна длины $L$ с двулучепреломлением она равна

$\Delta\tau = \frac{\Delta n_g}{c}\, L,$

где $\Delta n_g$ - групповое двулучепреломление, $c$ - скорость света. Здесь задержка растёт линейно с длиной - потому что ось двулучепреломления на коротком отрезке не меняется. На входе импульс расщепляется на две компоненты, и на выходе они разделены во времени на $\Delta\tau$: приёмник видит «двоящийся» или размытый символ.

Удобнее выражать DGD через двулучепреломление в единицах задержки на длину: $\Delta\tau / L$ измеряют в пс/км для однородного волокна. Например, при $\Delta n_g = 3\cdot10^{-7}$ задержка набегает примерно на 1 пс/км - мелочь на метрах, но заметная величина на сотнях километров.

Почему DGD растёт как корень из длины

В реальной линии ось двулучепреломления не постоянна: каждые несколько метров - десятков метров (длина корреляции $h$) она случайно поворачивается из-за изменения напряжений, температуры и изгибов. Это режим случайной связи мод (random mode coupling). Компоненты импульса то опережают, то отстают, и задержка накапливается не как сумма, а как случайное блуждание - статистически, как корень из числа участков.

Поэтому для длинной линии ($L \gg h$) средняя DGD растёт не линейно, а как корень из длины:

$\langle\Delta\tau\rangle = D_{\text{PMD}} \sqrt{L},$

где $D_{\text{PMD}}$ - PMD-коэффициент волокна, измеряемый в пс/$\sqrt{\text{км}}$. Это центральная формула инженерных расчётов PMD.

Рассчитать для своей задачи

Переход от линейного $\propto L$ к корневому $\propto \sqrt{L}$ происходит на масштабе длины корреляции $h$: на коротком отрезке (несколько метров) DGD ещё линейна, на длинной магистрали - уже корневая. Современные волокна имеют $D_{\text{PMD}} \approx 0{,}05{-}0{,}2$ пс/$\sqrt{\text{км}}$, старые (проложенные в 1980-х) - до $1{-}2$ пс/$\sqrt{\text{км}}$ и хуже.

Статистика Максвелла и мгновенная DGD

Важная тонкость: $D_{\text{PMD}}\sqrt{L}$ - это среднее, а мгновенная DGD случайна и зависит от длины волны, температуры и времени. При сильной связи мод мгновенная DGD подчиняется распределению Максвелла (как модуль скорости в кинетической теории газов):

$p(\Delta\tau) = \sqrt{\frac{2}{\pi}}\, \frac{\Delta\tau^2}{q^3}\, \exp\!\left(-\frac{\Delta\tau^2}{2q^2}\right),$

где параметр $q$ связан со средним как $\langle\Delta\tau\rangle = 2q\sqrt{2/\pi}$. Практическое следствие - у максвелловского распределения длинный хвост: мгновенная DGD изредка достигает 3-3,5 средних значений. Поэтому инженеры закладывают запас: допустимое среднее берут как $1/10$ длительности символа, чтобы вероятность выброса до критической DGD была около $10^{-5}$ и ниже.

Как PMD ограничивает линию связи

PMD размывает импульс, и при DGD порядка длительности бита соседние символы накладываются - растёт ошибка. Грубая оценка предельной дальности при заданной скорости $B$: допустимая DGD $\approx 0{,}1/B$, откуда

$L_{\max} \approx \left(\frac{0{,}1}{B\, D_{\text{PMD}}}\right)^2.$

Квадрат здесь - прямое следствие корневого закона: вдвое снизив $D_{\text{PMD}}$, выигрываем в дальности вчетверо. Для 10 Гбит/с и $D_{\text{PMD}} = 0{,}5$ пс/$\sqrt{\text{км}}$ предел - порядка 1600 км; при переходе на 40 Гбит/с он падает в 16 раз. Именно PMD, а не хроматическая дисперсия (которую компенсируют), стала жёстким барьером для апгрейда старых линий на высокие скорости.

В отличие от хроматической дисперсии, PMD нельзя скомпенсировать раз и навсегда фиксированным элементом - она дрейфует во времени. Поэтому применяют адаптивные PMD-компенсаторы с обратной связью либо цифровую компенсацию в DSP когерентного приёмника. Родственный эффект вращения поляризации в магнитном поле разбирается в заметке про эффект Фарадея.

PMD первого и второго порядка

То, что описано выше, - это PMD первого порядка: фиксированная DGD между двумя главными состояниями поляризации на данной длине волны. Но сами PSP и величина DGD зависят от частоты, и эта частотная зависимость - PMD второго порядка. Её раскладывают на два вклада:

• Поляризационно-зависимая хроматическая дисперсия - производная DGD по частоте $d\Delta\tau/d\omega$ добавляется к обычной хроматической дисперсии (со знаком, зависящим от поляризации).
• Деполяризация PSP - поворот главных состояний поляризации с частотой; широкий спектр импульса видит «размазанные» оси, и даже выровненный по PSP сигнал частично расщепляется.

Для широкополосных систем (40-100 Гбит/с и выше) второй порядок уже нельзя игнорировать: он растёт как $L$ (а не $\sqrt{L}$) и при больших длинах может доминировать над первым. Полное описание дают через вектор PMD $\vec{\Omega}(\omega)$, модуль которого равен DGD, а производная $d\vec{\Omega}/d\omega$ задаёт оба вклада второго порядка.

Как измеряют PMD

Стандартные методы измерения DGD и PMD-коэффициента:

• Метод фиксированного анализатора (wavelength scanning) - сканируют длину волны и считают число экстремумов в спектре прошедшей мощности; их плотность пропорциональна средней DGD.
• Интерферометрический метод (GINTY/TINTY) - по автокорреляционной функции поля восстанавливают распределение задержек; быстр и устойчив, годится для полевых измерений.
• Метод матрицы Джонса (JME) и анализ Мюллера - измеряют выходную поляризацию на серии длин волн и извлекают вектор PMD; дают и DGD, и ориентацию PSP.

Все методы дают статистическую оценку, усреднённую по длинам волн и времени, потому что мгновенная DGD флуктуирует.

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

• Путать линейный и корневой закон. На коротком однородном участке DGD $\propto L$, на длинной линии со связью мод $\langle\Delta\tau\rangle \propto \sqrt{L}$. Применять линейную формулу к магистрали - переоценить DGD в разы.
• Считать DGD постоянной. Мгновенная DGD случайна (Максвелл), $D_{\text{PMD}}\sqrt{L}$ - лишь среднее. Расчёт по среднему без запаса на хвост распределения занижает риск выбросов.
• Путать размерности PMD-коэффициента. $D_{\text{PMD}}$ имеет единицы пс/$\sqrt{\text{км}}$ (для длинной линии), а не пс/км. Под корнем - длина, поэтому и единица «нестандартная».
• Смешивать PMD с хроматической дисперсией. Хроматическая дисперсия зависит от длины волны и компенсируется фиксированно; PMD связана с поляризацией, дрейфует и требует адаптивной компенсации.
• Игнорировать связь мод при малой длине. Если длина линии меньше длины корреляции, статистика Максвелла неприменима, и DGD ведёт себя детерминированно-линейно.

FAQ

Чем PMD отличается от хроматической дисперсии? Хроматическая дисперсия - это зависимость групповой скорости от длины волны (разные спектральные компоненты импульса расходятся). PMD - зависимость скорости от поляризации (две ортогональные поляризации расходятся). Хроматическую дисперсию компенсируют фиксированными элементами, PMD дрейфует во времени и требует адаптивной или цифровой компенсации.

Почему DGD растёт как корень из длины, а не линейно? Из-за случайной связи мод: ось двулучепреломления случайно поворачивается вдоль линии, и задержка накапливается как случайное блуждание, а не как простая сумма. Случайное блуждание даёт зависимость $\propto\sqrt{L}$. Линейный рост сохраняется только на отрезке короче длины корреляции, где ось двулучепреломления постоянна.

Что такое PMD-коэффициент и в чём он измеряется? PMD-коэффициент $D_{\text{PMD}}$ - характеристика волокна, связывающая среднюю DGD с длиной по формуле $\langle\Delta\tau\rangle = D_{\text{PMD}}\sqrt{L}$. Измеряется в пс/$\sqrt{\text{км}}$. У современных волокон он около 0,05-0,2 пс/$\sqrt{\text{км}}$, у старых - на порядок больше.

Коротко

Поляризационная модовая дисперсия возникает из двулучепреломления волокна: две ортогональные поляризации бегут с разными групповыми скоростями и приходят с дифференциальной групповой задержкой (DGD) $\Delta\tau$. На коротком однородном участке DGD линейна по длине, но в реальной линии из-за случайной связи мод средняя задержка растёт как корень из длины: $\langle\Delta\tau\rangle = D_{\text{PMD}}\sqrt{L}$, где $D_{\text{PMD}}$ - PMD-коэффициент в пс/$\sqrt{\text{км}}$. Мгновенная DGD случайна и подчиняется распределению Максвелла с длинным хвостом, поэтому в проектировании закладывают запас. PMD ограничивает дальность высокоскоростных линий квадратично по $D_{\text{PMD}}$ и, в отличие от хроматической дисперсии, требует адаптивной компенсации.