Зависимость диэлектрика от направления искажает форму сигнала в линии

Микрополосковая линия выглядит обманчиво простой: узкая металлическая дорожка над земляным слоем, разделённая тонкой пластиной диэлектрика. Но электромагнитная волна, бегущая вдоль неё, живёт сразу в двух средах одновременно, часть поля в диэлектрике, часть в воздухе над дорожкой. Из этого раздвоения рождается дисперсия, зависимость скорости волны от частоты, которая на высоких частотах искажает форму сигнала. А если подложка ещё и анизотропна, то есть её диэлектрические свойства зависят от направления, картина усложняется многократно, и расчёт линии превращается в тонкую задачу.

Сапфир, излюбленный материал сверхвысокочастотной техники за стабильность и низкие потери, как раз анизотропен. Его диэлектрическая проницаемость вдоль кристаллической оси и поперёк неё отличается заметно, и волна в микрополоске на сапфире чувствует разную среду в зависимости от ориентации поля. Игнорировать это нельзя: ошибка в проницаемости прямо переводится в ошибку фазовой скорости, а та в рассогласование и искажение. Разберём, как складывается эффективная проницаемость линии, как она плывёт с частотой и как учитывают анизотропию подложки через введённый специально параметр.

Раздвоение поля между диэлектриком и воздухом над дорожкой

Главная особенность микрополоска в том, что его волна не заперта целиком в одной среде. Силовые линии электрического поля идут от дорожки к земляному слою частью сквозь диэлектрик подложки, частью через воздух над дорожкой и по бокам от неё. Получается смешанная среда, и эффективная диэлектрическая проницаемость линии оказывается между проницаемостью подложки и единицей для воздуха. Этот промежуточный показатель называют эффективной проницаемостью, и именно он определяет фазовую скорость волны.

Доля поля в каждой среде зависит от геометрии. Чем шире дорожка относительно толщины подложки, тем плотнее поле прижато к диэлектрику под дорожкой, и тем ближе эффективная проницаемость к проницаемости подложки. Узкая дорожка, наоборот, выпускает большую долю поля в воздух, и эффективная проницаемость падает. Отношение ширины дорожки к толщине подложки, обозначаемое W/h, становится главным геометрическим параметром, и все универсальные кривые микрополоска строят именно как функции этого отношения.

Фазовую скорость волны эффективная проницаемость задаёт прямо: v = c / sqrt(eps_eff), где c скорость света, eps_eff эффективная проницаемость. Из неё же вытекает и характеристическое сопротивление линии, и длина волны в ней. Для подложки с проницаемостью около десяти эффективная проницаемость типичной линии ложится в районе шести-семи, то есть волна бежит примерно в два с половиной раза медленнее, чем в вакууме. Любое смещение эффективной проницаемости немедленно сдвигает и скорость, и сопротивление, и потому точность её расчёта критична для согласования.

Почему эффективная проницаемость растёт с частотой

На низких частотах эффективную проницаемость можно считать постоянной, по квазистатическому приближению, когда длина волны много больше размеров сечения линии. Но с ростом частоты картина меняется. Чем выше частота, тем сильнее поле концентрируется в диэлектрике под дорожкой, вытягиваясь из воздуха. Доля энергии в подложке растёт, и эффективная проницаемость поднимается, постепенно приближаясь к проницаемости самого диэлектрика. Это и есть дисперсия микрополоска: эффективная проницаемость становится функцией частоты.

Физически концентрация поля связана с тем, что на высоких частотах волна всё хуже распространяется как простая поперечная и приобретает продольные составляющие, прижимающие энергию к диэлектрику. Эффективная проницаемость eps_eff(f) монотонно растёт от квазистатического значения на нуле частоты до полной проницаемости подложки на бесконечной частоте, проходя характерную точку перегиба. Скорость этого подъёма зависит от геометрии и от частоты, на которой поле начинает заметно перераспределяться, и эту частоту иногда называют частотой дисперсии.

Масштаб эффекта вполне ощутим в реальных цифрах. Для сапфировой подложки измеренная проницаемость падает с примерно 9.32 на частоте 3.3 гигагерца до 8.62 на частоте 37.6 гигагерца, что даёт дисперсию около 7.5 процента по диапазону. Для подложки из карбида кремния эффект сопоставим. Семь процентов изменения проницаемости означают около трёх с половиной процентов изменения фазовой скорости, и для длинной линии или резонатора это уже существенный сдвиг фазы и частоты, который обязательно учитывают при проектировании.

Эмпирические формулы дисперсии и их сравнение с моделированием

Точно посчитать дисперсию микрополоска можно полноволновым численным методом, решая уравнения поля во всём сечении на каждой частоте. Но это долго и неудобно для инженерной работы, поэтому на практике пользуются эмпирическими формулами, которые хорошо приближают результат. Классический подход описывает рост эффективной проницаемости с частотой через выражение вида eps_eff(f) = eps_r - (eps_r - eps_eff0) / (1 + G * (f / f_p)^2), где eps_eff0 квазистатическое значение, eps_r проницаемость подложки, f_p характеристическая частота, G геометрический коэффициент.

Формул такого рода предложено несколько, и они отличаются точностью и диапазоном применимости. Известны выражения, носящие имена своих авторов, и каждое настроено на свой класс подложек и частот. Сравнение с точным моделированием показывает, что для типичной керамической подложки одни формулы ложатся почти вплотную к численному результату, другие заметно расходятся на краях диапазона. Для подложки из глинозёма с проницаемостью 9.9, шириной дорожки 1.1 миллиметра и толщиной 1.0 миллиметр расхождение лучших эмпирических формул с программным расчётом остаётся в пределах долей процента вплоть до миллиметровых волн.

Выбор формулы зависит от задачи. Для грубой прикидки на низких частотах годится простейшее квазистатическое значение без дисперсии. Для широкополосной линии или работы в миллиметровом диапазоне берут проверенную дисперсионную формулу с подтверждённой точностью до десятков гигагерц. Универсальные кривые, построенные как функции отношения W/h, дают сразу три величины: эффективную проницаемость, её низкочастотный предел и характеристическое сопротивление линии, что позволяет спроектировать линию без полноволнового расчёта на каждом шаге.

Анизотропия подложки и введение приведённой проницаемости

Сапфир и подобные кристаллические подложки добавляют принципиально новую сложность: их проницаемость зависит от направления поля. В одноосном кристалле проницаемость вдоль оптической оси отличается от проницаемости поперёк неё, и поле микрополоска, имеющее составляющие в разных направлениях, чувствует разную среду. Это означает, что нельзя просто подставить одно число проницаемости в формулы для изотропной подложки, нужна особая трактовка.

Изящное решение состоит в том, чтобы свести анизотропную подложку к эквивалентной изотропной с некоторой приведённой проницаемостью. Для этого вводят специальный параметр приведённой проницаемости, который сам зависит от геометрии линии, от отношения W/h. Идея в том, что для каждой ширины дорожки подбирается такое эффективное изотропное число, которое даёт ту же фазовую скорость, что и реальная анизотропная подложка. Это приведённое значение затем подставляют в обычные изотропные формулы и кривые, пользуясь всем накопленным аппаратом расчёта изотропных линий.

Зависимость приведённой проницаемости от отношения W/h рассчитывают численно, например конечно-разностным методом, один раз для данного кристалла и его ориентации, а результат сводят в универсальные кривые. Введённое таким образом изотропизированное значение проницаемости можно прямо использовать в коммерческих программах моделирования, не имеющих встроенной поддержки анизотропии. Любопытно, что для правильно ориентированного сапфира сама анизотропия после такого приведения вносит менее половины процента разброса, тогда как частотная дисперсия даёт упомянутые семь с половиной процентов, то есть дисперсия оказывается более сильным эффектом, чем анизотропия.

Ориентация кристалла и распределение проницаемости по срезам

Анизотропный кристалл можно резать под разными углами к оптической оси, и от среза прямо зависит, какую проницаемость почувствует линия. Сапфир выпускают в нескольких стандартных ориентациях, обозначаемых по кристаллографическим плоскостям, и для каждой выведены свои расчётные формулы приведённой проницаемости. Правильно выбранная ориентация позволяет минимизировать влияние анизотропии или, наоборот, использовать её осознанно, выстраивая поле линии вдоль выгодного направления кристалла.

Практическое значение ориентации велико. Для одного среза проницаемость, которую видит линия, оказывается ближе к большему значению тензора, для другого к меньшему, и фазовая скорость соответственно различается. Если ориентация подложки не учтена в расчёте, готовая линия даст неверную электрическую длину, и устройство выйдет за полосу. Поэтому для каждой ориентации сапфира выведены простые и корректные формулы расчёта изотропной эффективной проницаемости с учётом дисперсии, проверенные сравнением с полноволновым анализом на основе эквивалентного поверхностного импеданса.

К ориентации добавляется температурная зависимость. Проницаемость сапфира сама плывёт с температурой, и для прецизионных устройств этот дрейф учитывают наряду с дисперсией и анизотропией. Современные модели приведённой проницаемости включают температурную поправку, чтобы расчётная электрическая длина линии оставалась верной в рабочем диапазоне температур. Совокупность трёх факторов, ориентации среза, частотной дисперсии и температуры, и определяет полную картину поведения проницаемости, которую закладывают в проект сверхвысокочастотного устройства на анизотропной подложке.

Искажение импульса при распространении вдоль дисперсной линии

До сих пор речь шла о гармонической волне одной частоты, но реальный сигнал это спектр частот, и дисперсия искажает его форму. Поскольку эффективная проницаемость растёт с частотой, разные спектральные составляющие импульса бегут с разной скоростью: низкочастотные быстрее, высокочастотные медленнее. В результате импульс, вошедший в линию резким и коротким, на выходе расплывается, его фронты затягиваются, а форма размывается. Это прямое следствие зависимости фазовой скорости от частоты, и для цифровых и импульсных применений оно критично.

Расчёт искажения ведут через спектр. Импульс раскладывают преобразованием Фурье на гармоники, каждую гармонику проводят через линию с её собственной фазовой задержкой, зависящей от эффективной проницаемости на этой частоте, и затем собирают сигнал обратным преобразованием. Для гауссова импульса с полной шириной на половине высоты порядка десяти пикосекунд на сапфировой подложке с поперечной проницаемостью около 11 искажение формы при прохождении заметной длины линии хорошо предсказывается этим методом и подтверждается измерением.

Анизотропия добавляет к искажению свой вклад, потому что она смещает эффективную проницаемость, а вместе с ней и величину дисперсии. Учёт анизотропии, проводника и дисперсии вместе позволяет точно предсказать расплывание и цифрового прямоугольного, и гауссова импульса на анизотропной линии. Практический вывод прост: чем шире спектр сигнала и чем длиннее линия, тем сильнее проявляется дисперсионное искажение, и для широкополосных трактов на анизотропных подложках длину линий стараются минимизировать либо компенсировать дисперсию специальными корректирующими структурами.

Сведение всех факторов проницаемости в единый расчёт

Соберём картину поведения микрополоска на анизотропной подложке воедино. Эффективная проницаемость рождается из раздвоения поля между диэлектриком и воздухом и зависит от геометрии через отношение W/h. С частотой она растёт от квазистатического значения к проницаемости подложки, и этот рост описывают эмпирической формулой вида eps_eff(f) = eps_r - (eps_r - eps_eff0) / (1 + G * (f / f_p)^2). Анизотропию сводят к эквивалентной изотропной среде через приведённую проницаемость, зависящую от W/h и ориентации среза, а температурную зависимость учитывают отдельной поправкой.

Числовые ориентиры расставляют приоритеты. Частотная дисперсия сапфира даёт около 7.5 процента изменения проницаемости от единиц до десятков гигагерц, тогда как анизотропия после правильного приведения вносит менее половины процента, то есть дисперсия доминирует. Фазовая скорость v = c / sqrt(eps_eff) меняется примерно на половину дисперсии проницаемости, и для длинных линий это заметный сдвиг фазы. Искажение импульса считают спектральным методом, проводя каждую гармонику с её фазовой задержкой и собирая сигнал обратно.

Грамотный расчёт линии на анизотропной подложке это согласование всех этих факторов сразу. Геометрию W/h выбирают из требований к характеристическому сопротивлению и допустимой дисперсии. Ориентацию кристалла подбирают так, чтобы минимизировать или осознанно использовать анизотропию. Дисперсионную формулу берут с подтверждённой точностью до рабочей частоты. Температурную поправку закладывают для прецизионных устройств. Опытный разработчик сводит отношение ширины к толщине, ориентацию среза, дисперсионную модель и температурную зависимость в единый расчёт, где эффективная проницаемость и фазовая скорость предсказаны точно во всём рабочем диапазоне, согласование сохраняется, а искажение сигнала остаётся в допустимых пределах.