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 и медицинские эндоскопы, а концепция метаматериалов успешно расширяется на акустические, сейсмические, тепловые и пространственно-временные метаматериалы.