Глава 9. Расчётная программа и численная верификация

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

Такой подход принципиально усиливает монографию. Сильная научная теория отличается не только глубиной идеи, но и способностью точно указать, каким именно расчётом её утверждения должны быть подтверждены или опровергнуты. Поэтому задача настоящей главы не сводится к краткому перечислению будущих работ. Её цель состоит в построении строгой дорожной карты, в которой Monte Carlo, Helmholtz, Maxwell, Schrödinger и карты чувствительности выступают как последовательные уровни проверки, а не как случайный набор разрозненных симуляций.

Важнейшая особенность псевдопараболоидов состоит в том, что для них не существует одной-единственной метрики, автоматически закрывающей вопрос об аттракторном механизме. Удержание, выход через апертуру, направленность, модовая селекция, добротность, время жизни и робастность описывают разные стороны работы резонатора. Глава 9 служит ещё и инструментом научной прозрачности. Она заранее фиксирует, какой именно вычислительный пакет должен сопровождать будущие результаты, какие карты нужно строить, какие данные сохранять и в каких случаях режим может считаться действительно закрытым.

9.1. Место главы 9 в общей логике монографии и принцип поэтапного закрытия критериев

Необходимо показать, каким образом каждый вычислительный блок связан с критериями C4_pp, C5_pp, C7_pp и C8_pp. Только при таком подходе список будущих задач превращается в строгую лестницу доказательства, а не в набор пожеланий. Сначала должен быть закрыт вопрос о совместимости удержания и вывода в открытом режиме. Затем — вопрос о направленности и модовой селекции. После этого — вопрос о межфизической переносимости между Helmholtz, Maxwell и Schrödinger. И лишь затем — вопрос о запасе устойчивости и практической жизнеспособности режима.

Эта последовательность имеет физический смысл. Невозможно обсуждать межфизическую универсальность до тех пор, пока хотя бы одна волновая постановка не дала полноценных карт η_out и θ_div. Точно так же невозможно содержательно вводить робастность, если не определено, какие именно режимы считаются рабочими и какие наблюдаемые должны сохраняться под возмущением. Особое значение имеет различение между численной демонстрацией и численной верификацией. Демонстрацией может быть отдельная удачная симуляция в одном режиме. Верификация же требует сканирования по безразмерным параметрам, контроля сеточной и внешней сходимости, проверки чувствительности к малым возмущениям и явной фиксации критериев принятия или отклонения режима. Эта оговорка нужна потому, что именно в теории псевдопараболоидов соблазн очень велик: красивая геометрия и выразительные поля легко создают иллюзию уже доказанной универсальности. Глава 9 сознательно снимает эту иллюзию и заменяет её формальным вычислительным стандартом.

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

9.2. Open-Monte Carlo как первый этап открытой верификации: что именно он должен измерять и чего не может доказать

Первым этапом после замкнутой Monte Carlo-калибровки должен стать переход к open-Monte Carlo, в котором часть отражающей границы переводится в выходную. Этот шаг не даёт полного волнового решения, но он необходим как промежуточный геометрический мост между главой 4 и полноволновыми расчётами раздела II. Его задача состоит в том, чтобы впервые на лучевом уровне проверить компромисс между удержанием и выводом: насколько сильно апертура разрушает аттракторный режим, как меняется вероятность достижения активной зоны и какие угловые направления начинают доминировать у траекторий, покидающих полость.

Open-Monte Carlo не должен восприниматься как быстрый способ закрыть критерии C4_pp и C5_pp. По своей природе он не видит дифракции на кромке щели, интерференционного распределения поля по апертуре, поляризационной селективности и векторной структуры моды. Поэтому его физическая функция иная: он должен выделить перспективные области параметров и отсеять заведомо бесперспективные геометрические режимы до запуска полноволновых расчётов.

Чтобы этот этап был действительно полезным, он должен извлекать не одну, а несколько согласованных величин. Для каждой траектории следует регистрировать достижение аттракторной зоны, число последовательных отражений в активном режиме, факт выхода через апертуру, координату точки выхода и угловое направление покидающей траектории. Из этих данных уже можно построить лучевые аналоги η_center, η_out и D_axis, а также карту плотности выхода по апертуре. Последняя особенно важна, поскольку она сразу показывает различие между локальными полярными окнами вертикальной топологии и пространственно распределённой кольцевой щелью горизонтальной формы.

Open-Monte Carlo должен запускаться не в одной точке параметров, а в сетке по K и χ, как минимум для нескольких сценариев возбуждения: центрального, объёмного и специально ориентированного на апертуру. Только в этом случае он даёт первые карты Φ_open, ray(K, χ), пригодные для выбора диапазонов, в которых запуск Helmholtz и Maxwell действительно имеет смысл. Иными словами, первый этап нужен не для окончательного ответа, а для рациональной фильтрации пространства параметров.

η_center = E_attractor / E_in

η_out = P_out / P_in

D_axis(θ_max) = P(θ ≤ θ_max) / P_out

S = P_side_lobes / P_main_lobe

Хотя окончательное значение этих величин определяется уже в полноволновой постановке, именно они должны стать сквозным языком всей расчётной программы. Даже промежуточные методы обязаны работать не в произвольных показателях, а в тех наблюдаемых, которые сохранятся до финального уровня Helmholtz-, Maxwell- и Schrödinger-верификации.

9.3. Полноволновой блок Helmholtz как первый обязательный эталон для карт открытый режим

Вторая ступень расчётной программы — систематическое сканирование открытой Helmholtz-задачи. Именно этот блок должен стать первым полноволновым эталоном теории, поскольку он уже включает интерференцию, дифракцию и квазимодовую структуру, но ещё не перегружен векторными осложнениями Maxwell-постановки. Если после open-Monte Carlo следующий шаг не ведёт к Helmholtz-картам, теория остаётся в промежуточном состоянии между геометрооптикой и настоящей волновой физикой.

Глава требует, чтобы Helmholtz-блок выполнялся не как набор отдельных примеров, а как сканирование по трём безразмерным параметрам: K, χ и ka. Для каждой точки сетки должны извлекаться распределение интенсивности в активной зоне, поток через апертуру, дальнеполевая диаграмма, ширина главного лепестка, уровень боковых лепестков и признаки существования устойчивой квазимоды. Только после этого становится возможным строить карты η_center(K, χ, ka), η_out(K, χ, ka), θ_div(K, χ, ka) и S_dB(K, χ, ka).

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

С вычислительной точки зрения Helmholtz-блок обязан включать и собственную, и вынужденную постановки. Собственная задача выявляет спектральные окна и квазимоды. Вынужденная — показывает реальный поток через щель и внешнюю диаграмму при конкретном возбуждении. Только их совместное использование позволяет избежать ошибки, когда собственное поле ошибочно принимается за доказательство полезного открытого режима. Для псевдопараболоидов это особенно важно, потому что сильная внутренняя локализация ещё не гарантирует инженерно полезного вывода наружу.

9.4. Maxwell-блок как решающий тест направленности, модовой селекции и 3D-правдоподобия теории

Третья ступень программы — полный электродинамический расчёт. В логике всей монографии именно он является решающим, поскольку проверяет самые сильные и самые уязвимые утверждения теории: возможность узкого осевого вывода из кольцевой щели, реальность модового условия m = 0, влияние поляризации и устойчивость режима к векторным и трёхмерным эффектам. Если Helmholtz ещё можно рассматривать как скалярный эталон, то Maxwell показывает, выживает ли геометрически предложенный механизм в максимально строгой практической волновой постановке.

Недостаточно просто повторить сканирование с учётом поляризации. Требуется полноценная 3D полновекторная постановка, разложение поля по азимутальным гармоникам, вычисление доли компоненты m = 0 на апертуре, анализ боковых лепестков и обязательная проверка параметрической устойчивости PML. Только такой пакет расчётов способен превратить апертурную оценку главы 6 в подтверждённый или опровергнутый результат.

Особенно критично это для горизонтальной кольцевой щели. Теория уже показала, что геометрически большое кольцо может давать узкий осевой максимум только при сохранении квазиосесимметричной фазы. Однако само по себе большое значение kρ_s ещё не гарантирует, что именно m = 0 будет доминировать в реальной 3D-моде. Это должно быть извлечено из расчёта как самостоятельная наблюдаемая величина. Следовательно, Maxwell-блок здесь выполняет методологическую функцию: он отделяет геометрически соблазнительные режимы от действительно физически реализуемых.

Вычислительный протокол Maxwell должен включать сканирование по K, χ и ka, варьирование типа возбуждения, контроль сеточной сходимости, тесты по параметрам PML и проверку устойчивости диаграммы при малом нарушении симметрии. Режим можно считать вычислительно закрытым только тогда, когда одновременно стабилизированы Q или комплексная модовая добротность, η_out, θ_div, S_dB и доля m = 0. Эта позиция усиливает монографию, поскольку заранее исключает выдачу частично подтверждённого результата за уже доказанную универсальность.

9.5. Schrödinger-блок и условие межфизической переносимости: от отдельных расчётов к пересечению U_AC, U_EM и U_Q

Четвёртая ступень программы — квантовая постановка Schrödinger. В контексте всей монографии она играет не роль экзотического дополнения, а роль самого строгого теста межфизической переносимости. Если для одной и той же безразмерной геометрии удаётся обнаружить совместимое поведение в Helmholtz, Maxwell и Schrödinger-постановках, тогда утверждение о геометрически организованном аттракторном механизме становится значительно сильнее. Однако именно поэтому квантовый блок нельзя проводить формально или эпизодически: он должен быть встроен в ту же систему параметров и наблюдаемых, что и две предыдущие физики.

Базовым вариантом для межфизического сравнения должна оставаться геометрическая ловушка с жёсткими границами и с открытым каналом, максимально близким по смыслу к апертуре классических волн. Только при этом условии можно честно сравнивать долю локализованной вероятности, поток через канал и время жизни квазистационарного состояния с аналогами η_center, η_out и Q в двух других физиках.

Главная цель Schrödinger-блока состоит не в построении квантового примера, а в локализации множества U_Q в пространстве параметров. Это множество должно определяться теми же условиями, что и для остальных физик: наличие локализации в активной зоне, ненулевая, но контролируемая утечка через канал, приемлемая структура потока и положительный запас устойчивости к малым геометрическим возмущениям. Только после этого становится содержательным главный вопрос всей монографии: существует ли непустое пересечение U_AC ∩ U_EM ∩ U_Q.

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

9.6. Карты чувствительности, вычислимый объект ε* и критерии завершённости программы

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

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

В вычислительном смысле ε* удобно трактовать как максимальную норму допустимого совокупного возмущения, при которой базовые наблюдаемые остаются выше заранее заданных порогов: η_center не падает ниже минимума, η_out остаётся ненулевой и контролируемой, θ_div не превышает допустимого значения, а S_dB не выходит за пределы установленной инженерной нормы. Такая постановка делает ε* естественным завершающим объектом всей монографии, связанным и с открытым режимом, и с направленностью, и с межфизической переносимостью.

После построения карт чувствительности можно сформулировать строгие критерии завершённости всей численной программы. Режим считается вычислительно закрытым только тогда, когда выполнены одновременно пять условий: сканирование по K, χ и ka завершено с контролем сходимости; построены карты η_center, η_out, θ_div и S_dB; для Maxwell извлечена доля m = 0 и проверена устойчивость к PML и малой асимметрии возбуждения; для Schrödinger построена область U_Q; и, наконец, вычислен положительный запас ε* в согласованной норме возмущений. Лишь такое сочетание делает выводы монографии действительно зрелыми.

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

9.7. Выводы по главе

Глава 0 показывает, что окончательная научная сила теории псевдопараболоидов зависит не от одной удачной симуляции, а от последовательной программы расчётов, идущей от open-Monte Carlo к Helmholtz, затем к 3D Maxwell, далее к Schrödinger и, наконец, к картам чувствительности и вычислению ε*.

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

Следующий логический шаг после настоящей главы состоит в демонстрации того, как эта программа работает на конкретных проектных параметрах.