1. Введение: Инженерия свойств материалов с помощью геометрии вместо химии

На протяжении столетий базовой аксиомой материаловедения и классической оптики считалось утверждение, что электромагнитные свойства вещества — почему стекло пропускает свет, золото обладает металлическим блеском, а соломинка в воде кажется надломленной — определяются исключительно его химическим составом: типом атомов из периодической таблицы и характером химических связей. На рубеже XXI века на стыке электродинамики, физики конденсированного состояния и нанотехнологий произошла смена парадигмы — возникли метаматериалы. От греческой приставки meta- («сверх», «за пределами»), метаматериал представляет собой искусственно структурированный композит, чей макроскопический электромагнитный отклик определяется не атомной химией, а субволновой геометрией, топологией и пространственным расположением искусственных структурных ячеек («мета-атомов»).

Видимый свет охватывает длины волн от \(\lambda \approx 400\,\text{нм}\) до \(700\,\text{нм}\), тогда как межатомное расстояние в кристаллах составляет всего \(a_{\text{atom}} \sim 0.1\text{–}0.3\,\text{нм}\). Поскольку электромагнитное поле осциллирует на масштабах, в тысячи раз превышающих размеры отдельных атомов, световая волна воспринимает сплошную среду с усреднёнными макроскопическими параметрами: диэлектрической проницаемостью \(\varepsilon\) и магнитной проницаемостью \(\mu\), определяющими показатель преломления \(n = \sqrt{\varepsilon_r \mu_r}\). Если создать массив искусственных резонаторов с периодом решётки \(a \lesssim \lambda / 10\), падающая волна усредняет их локальные дипольные моменты и воспринимает структуру как однородную эффективную среду с заданными параметрами \(\varepsilon_{\text{eff}}(\omega)\) и \(\mu_{\text{eff}}(\omega)\).

Принципиально важно отличать метаматериалы от фотонных кристаллов и природной структурной окраски (опал, крылья бабочек Morpho). В фотонных кристаллах период решётки соизмерим с полуволной (\(a \sim \lambda / 2\)), а запрещённые зоны формируются за счёт брэгговской дифракции, где приближение эффективной среды (EMA) неприменимо. Метаматериалы же работают строго в глубоко субволновом режиме (\(a \ll \lambda\)), где локальные тензоры \(\varepsilon_{\text{eff}}\) и \(\mu_{\text{eff}}\) имеют строгий физический смысл.

2. Уравнения Максвелла и четыре квадранта плоскости \(\varepsilon\text{–}\mu\)

Для монохроматических плоских волн \(\exp(i\mathbf{k}\cdot\mathbf{r} - i\omega t)\) в однородной изотропной среде без источников роторные уравнения Максвелла имеют вид:

где \(\varepsilon = \varepsilon_0 \varepsilon_r\) и \(\mu = \mu_0 \mu_r\). Исключая \(\mathbf{H}\), получаем дисперсионное соотношение и волновой импеданс \(Z\):

Плоскость \((\varepsilon_r, \mu_r)\) делит все изотропные среды на четыре квадранта:

  flowchart TD
    Root["Электромагнитные параметры среды (ε_r, μ_r)"] --> Q1["Квадрант I: ε_r > 0, μ_r > 0<br/>[Дважды положительные среды (DPS) / Правосторонние]"]
    Root --> Q2["Квадрант II: ε_r < 0, μ_r > 0<br/>[Эпсилон-отрицательные среды (ENG) / Электрическая плазма]"]
    Root --> Q4["Квадрант IV: ε_r > 0, μ_r < 0<br/>[Мю-отрицательные среды (MNG) / Магнитные резонаторы]"]
    Root --> Q3["Квадрант III: ε_r < 0, μ_r < 0<br/>[Дважды отрицательные метаматериалы (DNG) / Левосторонние]"]

    Q1 --> Q1_P["n = +√(ε_r μ_r) > 0<br/>Бегущие волны (k сонаправлен с S)<br/>Примеры: Стекло, вода, воздух, диэлектрики"]
    Q2 --> Q2_P["n² = ε_r μ_r < 0 (k чисто мнимое)<br/>Эванесцентное затухание и металлическое отражение<br/>Примеры: Благородные металлы (Au, Ag) ниже ω_p"]
    Q4 --> Q4_P["n² = ε_r μ_r < 0 (k чисто мнимое)<br/>Эванесцентное затухание и магнитное экранирование<br/>Примеры: Ферриты вблизи ферромагнитного резонанса"]
    Q3 --> Q3_P["n = -√(ε_r μ_r) < 0<br/>Обратные бегущие волны (k антипараллелен S)<br/>Примеры: Искусственные метаматериалы SRR + проволочки"]

В природе Квадрант III пустует, поскольку электрические дипольные и плазменные резонансы лежат в оптическом и УФ-диапазонах (\(10^{14}\text{–}10^{16}\,\text{Гц}\)), а магнитные резонансы подавлены постоянной тонкой структуры \(\alpha \approx 1/137\) и затухают выше СВЧ- и ТГц-частот.

3. Мысленный эксперимент В. Г. Веселаго (1967–1968) и физика отрицательного преломления

В 1967 году (английский перевод в 1968 г.) советский физик Виктор Георгиевич Веселаго (ФИАН) впервые строго исследовал электродинамику среды с одновременными \(\varepsilon_r < 0\) и \(\mu_r < 0\).

Из \(\mathbf{k} \times \mathbf{E} = \omega \mu \mathbf{H}\) при \(\varepsilon < 0, \mu < 0\) следует, что тройка векторов \((\mathbf{E}, \mathbf{H}, \mathbf{k})\) образует левую тройку, почему Веселаго назвал такие среды «левыми» веществами (Left-Handed Materials, LHM). Однако поток электромагнитной энергии определяется вектором Пойнтинга \(\mathbf{S}\):

Так как в определение \(\mathbf{S}\) не входят \(\varepsilon\) и \(\mu\), векторы \((\mathbf{E}, \mathbf{H}, \mathbf{S})\) всегда образуют правую тройку! Следовательно, волновой вектор \(\mathbf{k}\) (фазовая скорость \(\mathbf{v}_p\)) и вектор Пойнтинга \(\mathbf{S}\) (групповая скорость \(\mathbf{v}_g\)) антипараллельны (\(\mathbf{k} \cdot \mathbf{S} < 0\), обратная волна). Комплексный анализ для пассивной причинной среды (\(\varepsilon'', \mu'' > 0\)) также однозначно требует выбора отрицательной ветви квадратного корня для выполнения условия затухания \(\text{Im}(n) > 0\):

Это приводит к трём обращённым оптическим эффектам:

1. Обращённый закон Снеллиуса и плоская линза Веселаго: Из непрерывности \(k_x\) и условия \(S_z > 0\) преломлённый луч отклоняется по ту же сторону от нормали, что и падающий: \( n_1 \sin\theta_1 = n_2 \sin\theta_2 = -|n_2|\sin|\theta_2| \) Плоскопараллельная пластина с \(n = -1\) (\(\varepsilon_r = \mu_r = -1\)) дважды фокусирует расходящиеся лучи точечного источника (внутри и снаружи пластины) без отражения от границ (\(Z = Z_0\sqrt{(-1)/(-1)} = Z_0\)).
2. Обратный эффект Доплера: Приближающийся к наблюдателю источник света испытывает красное смещение \(\omega' = \omega_0 (1 - |n|v/c)\).
3. Обратное излучение Вавилова — Черенкова: Быстрая заряженная частица (\(v > c/|n|\)) излучает конус ударной волны назад под тупым углом к траектории.

4. Прорыв сэра Джона Пендри (1996–1999): Разомкнутые кольцевые резонаторы (SRR)

В 1996–1999 годах сэр Джон Пендри (Sir John B. Pendry) из Имперского колледжа Лондона показал, как конструировать \(\varepsilon_{\text{eff}} < 0\) и \(\mu_{\text{eff}} < 0\) из немагнитных проводников:

• Решётка тонких металлических проволочек (1996): Пространственное разбавление и огромная магнитная самоиндукция придают электронам гигантскую эффективную массу \(m_{\text{eff}} = \frac{\mu_0 e^2 n_e \pi r^2}{2\pi}\ln(a/r) \gg m_e\), снижая плазменную частоту до ГГц-диапазона: \( \omega_{pe}^2 = \frac{2\pi c^2}{a^2 \ln(a/r)}, \qquad \varepsilon_{\text{eff}}(\omega) = 1 - \frac{\omega_{pe}^2}{\omega(\omega + i\gamma_e)} \)
• Разомкнутый кольцевой резонатор (Split-Ring Resonator, SRR, 1999): Два концентрических медных кольца с ёмкостными зазорами работают как субволновой планарный LC-контур (\(\omega_m = 1/\sqrt{LC}\)). Чуть выше \(\omega_m\) наведённый магнитный момент осциллирует в противофазе с внешним полем, создавая отрицательную магнитную проницаемость: \( \mu_{\text{eff}}(\omega) = 1 - \frac{F \omega^2}{\omega^2 - \omega_m^2 + i\gamma_m \omega} \)

В 2000–2001 годах группа Дэвида Смита (David R. Smith) в Калифорнийском университете в Сан-Диего (UCSD) объединила проволочки и SRR, впервые экспериментально измерив на частоте \(10.5\,\text{ГГц}\) отрицательный показатель преломления \(n = -2.7 \pm 0.1\).

5. Преодоление дифракционного предела: Суперлинза Пендри и гиперлинзы

Обычная оптика ограничена дифракционным пределом Аббе (\(\Delta x = \lambda / (2\,\text{NA})\)). Для субволновых пространственных частот (\(k_\parallel > k_0 = 2\pi/\lambda\)) продольное волновое число становится чисто мнимым:

Эти эванесцентные (затухающие) волны экспоненциально спадают как \(\exp(-\kappa_z z)\). В 2000 году Пендри (Physical Review Letters) доказал, что слой с \(\varepsilon_r \to -1, \mu_r \to -1\) за счёт резонансного возбуждения поверхностных плазмон-поляритонов (SPP) экспоненциально усиливает эванесцентные волны:

В 2005 году группа Сяна Чжана (Xiang Zhang) в UC Berkeley продемонстрировала серебряную суперлинзу толщиной \(35\,\text{нм}\) на длине волны \(\lambda = 365\,\text{нм}\) с разрешением \(60\,\text{нм}\) (\(\sim \lambda/6\)), а затем цилиндрические гиперлинзы с гиперболической дисперсией (\(\frac{k_r^2}{\varepsilon_\theta} - \frac{k_\theta^2}{|\varepsilon_r|} = \frac{\omega^2}{c^2}\)), проецирующие субволновые изображения в дальнюю зону.

6. Трансформационная оптика и математика плаща-невидимки

В мае 2006 года Пендри, Шуриг и Смит, а также Ульф Леонхардт опубликовали в Science основы трансформационной оптики (TO), опирающейся на форм-инвариантность уравнений Максвелла относительно гладких преобразований координат \(x^i \mapsto x^{i'}\) (матрица Якоби \(\Lambda^{i'}_i = \frac{\partial x^{i'}}{\partial x^i}\)):

  flowchart TD
    S1["1. Однородное виртуальное пространство (r, θ, φ)<br/>Плоская евклидова сетка, прямолинейные лучи<br/>Начало координат r = 0 — безразмерная точка"] --> S2["2. Радиальное сжатие координат<br/>r' = R_1 + r (R_2 - R_1) / R_2<br/>Раздувание точки r = 0 в сферическую полость радиуса R_1"]
    S2 --> S3["3. Отображение тензора Якоби<br/>ε' = Λ ε Λ^T / det(Λ), μ' = Λ μ Λ^T / det(Λ)<br/>Расчёт анизотропных радиальных и угловых профилей среды"]
    S3 --> S4["4. Физический метаматериальный плащ (R_1 ≤ r' ≤ R_2)<br/>Внутренняя полость (r' < R_1): Полностью изолированная скрытая область<br/>Внешняя область (r' > R_2): Волновые фронты выходят без тени и отражений"]

Сжимая шар \(0 \le r \le R_2\) в сферический слой \(R_1 \le r' \le R_2\) преобразованием \(r' = R_1 + r\frac{R_2 - R_1}{R_2}\), получаем параметры сферического плаща Пендри:

На внутренней границе \(r' = R_1\) радиальные компоненты обращаются в нуль (\(\varepsilon_r = \mu_r \to 0\)), плавно огибая светом скрытую полость \(r' < R_1\). В октябре 2006 года в Университете Дьюка был создан первый микроволновый плащ-невидимка на частоте \(8.5\,\text{ГГц}\). В дальнейшем были разработаны ковровые плащи (Carpet Cloaks), плазмонная компенсация рассеяния (Алу и Энгета) и геометрические параксиальные плащи (Рочестерский плащ).

7. Почему «плащ-невидимка Гарри Поттера» запрещён законами физики

Всенаправленный пассивный плащ-невидимка без потерь для макроскопического 3D-объекта в широком видимом спектре (\(400\text{–}700\,\text{нм}\)) невозможен из-за трёх фундаментальных запретов:

1. Релятивистская причинность и предел полосы пропускания Крамерса — Кронига: Чтобы обогнуть полость по дуге \(\pi R_1 > 2R_1\) в фазе с внешней волной, фазовая скорость внутри плаща должна быть сверхсветовой (\(v_p = c/n(\omega) > c\)). Однако групповая скорость переноса энергии ограничена условием \(v_g \le c\): \( v_g(\omega) = \frac{c}{n(\omega) + \omega \frac{dn(\omega)}{d\omega}} \le c \quad \Longrightarrow \quad \omega \frac{dn(\omega)}{d\omega} \ge 1 - n(\omega) > 0 \) В сочетании с соотношениями Крамерса — Кронига: \( \text{Re}\,\varepsilon_r(\omega) - 1 = \frac{2}{\pi} \mathcal{P} \int_0^\infty \frac{\omega' \, \text{Im}\,\varepsilon_r(\omega')}{\omega'^2 - \omega^2} d\omega' \) относительная полоса частот маскировки ограничена величиной \(\Delta\omega / \omega_0 \lesssim \lambda_0 / R_1\) (\(\sim 5\times 10^{-7}\) для \(R_1 = 1\,\text{м}\) в видимом свете; правило сумм Монтиконе и Алу).
2. Плазмонные омические потери и насыщение кинетической индуктивности: На оптических частотах (\(400\text{–}750\,\text{ТГц}\)) кинетическая индуктивность электронов (\(L_{\text{kin}} \propto m_e / (n_e e^2)\)) насыщает магнитный резонанс, а поглощение в металлах снижает добротность до \(\text{FOM} = |\text{Re}(n)|/\text{Im}(n) \sim 1\text{–}3\).
3. Технологический барьер 3D-нанофабрикации: Послойное создание тысяч выровненных 3D-нанослоёв с разрешением \(10\text{–}50\,\text{нм}\) на метровых площадях недоступно современной литографии.

8. От 3D-метаматериалов к 2D-диэлектрическим метаповерхностям и металинзам

В 2010-х годах произошёл переход к двумерным метаповерхностям (Metasurfaces) на основе резонаторов Ми из высокопреломляющих диэлектриков без потерь (\(\text{Si}, \text{TiO}_2, \text{GaN}\)).

В 2011 году Федерико Капассо (Federico Capasso) и Наньфан Юй (Гарвард) сформулировали обобщённый закон Снеллиуса, где пространственный градиент фазы \(\frac{d\Phi(x)}{dx}\) управляет преломлением:

При выполнении первого условия Керкера (\(\alpha_e = \alpha_m\)) метаповерхности Гюйгенса подавляют обратное отражение, обеспечивая почти \(100\%\)-ное пропускание. Формируя гиперболический фазовый профиль идеальной линзы:

и управляя групповой задержкой \(\frac{\partial\Phi}{\partial\omega}\), плоские ахроматические металинзы толщиной менее \(1\,\text{мкм}\) внедряются в 3D-распознавание лиц смартфонов, твердотельные лидары (LiDAR), очки AR/VR и медицинские эндоскопы, а концепция метаматериалов успешно расширяется на акустические, сейсмические, тепловые и пространственно-временные метаматериалы.