EN KR JP CN RU IN
Physical Review Research

Динамическая остановка в турбулентности активных нематиков

Authors Ido Lavi, Ricard Alert, Jean-François Joanny, Jaume Casademunt
Original Paper Published 2026-03-17
ISOM Posted 2026-03-19 20:10 UTC
Read Time 7M
Open PDF Open MICCAI page

Предыстория и академическая родословная

Происхождение и академическая родословная

Исследование активных жидкостей — материалов, приводимых в движение внутренними компонентами, такими как молекулярные моторы, клетки или синтетические частицы, — является динамичной областью исследований на протяжении десятилетий. Эта внутренняя активность порождает спонтанные, часто хаотические потоки, коллективно называемые активной турбулентностью. Это явление наблюдалось в широком спектре систем, включая бактериальные суспензии, сперматозоиды, цитоскелетные смеси, клеточные монослои и искусственные самодвижущиеся частицы [1, 2-31]. Исторически эти потоки, несмотря на то, что они происходят при низких числах Рейнольдса, где инерция пренебрежимо мала, поразительно напоминают классическую инерционную турбулентность.

Недавние достижения в области активной нематической турбулентности — специфического типа активных жидкостей, где удлиненные частицы имеют тенденцию к выравниванию — выявили универсальные законы масштабирования для определенных особенностей, таких как спектр мощности скорости [22,32,33]. Например, было обнаружено, что масштабирование спектра мощности скорости $q^{-1}$ при малых волновых числах $q$ является удивительно устойчивым, независимым от специфических свойств материала, таких как вязкость или активность, при условии сохранения низкого числа Рейнольдса [34]. Это универсальное поведение наблюдалось как в нематиках с дефектами, так и в сильно упорядоченных нематиках без дефектов [32,33].

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

Интуитивные термины предметной области

Чтобы помочь полному новичку интуитивно понять специализированные термины, используемые в этой статье, здесь приведены некоторые повседневные аналогии:

  • Активные жидкости: Представьте себе миску супа, где каждый крошечный ингредиент (например, кусочек макарон или овощей) имеет свой собственный миниатюрный мотор и постоянно движется, перемешивая суп изнутри. В отличие от обычного супа, который нужно перемешивать ложкой, этот суп полностью самоперемешивающийся.
  • Нематический порядок: Представьте поле высокой травы. В нормальном, неупорядоченном состоянии стебли могут указывать в разные стороны. Но при нематическом порядке большинство стеблей травы имеют тенденцию выравниваться примерно в одном направлении, как будто они все вместе колышутся на легком бризе. Это предпочтительное коллективное направление ученые называют "директором".
  • Активная нематическая турбулентность: Теперь объедините самоперемешивающийся суп с выравнивающейся травой. Если внутренние моторы очень мощные и хаотичные, они не просто мягко выравниваются; они создают вихревые, непредсказуемые и высокодинамичные течения, очень похожие на бурные речные пороги или хаотичный шлейф дыма, но все это генерируется самими ингредиентами. Это и есть активная нематическая турбулентность.
  • Топологические дефекты: В нашем поле выровненной травы топологический дефект подобен небольшому локализованному "узлу" или "водовороту", где стебли травы внезапно теряют свое коллективное выравнивание и закручиваются вокруг центральной точки или указывают в противоречивых направлениях. Это точка интенсивного беспорядка в иначе упорядоченной системе.
  • Динамическая остановка: Представьте себе оживленное шоссе с постоянно движущимися и перестраивающимися машинами, создающими плавный, постоянно меняющийся дорожный узор. Динамическая остановка похожа на то, как этот транспортный поток внезапно замирает в фиксированной, древовидной сети "дорог" и "потоков". Машины все еще движутся по этим путям, но сама общая структура дорожной сети становится стабильной и перестает драматически эволюционировать. Хаотичная эволюция системы "замораживается" в устойчивый, организованный паттерн.

Таблица обозначений

Обозначение Тип Описание
$\mathbf{v}$ Переменная Поле скорости потока
$\mathbf{n}$ Переменная Поле нематического директора (локальная средняя ориентация стержнеобразных частиц)
$\psi$ Переменная Функция тока (описывает несжимаемый поток)
$s$ Переменная Скалярный параметр порядка (квантифицирует локальную силу нематического выравнивания)
$Q_{\alpha\beta}$ Переменная Тензор нематического параметра порядка (описывает нематический порядок и выравнивание)
$K$ Параметр Упругая постоянная (измеряет сопротивление искажениям поля директора)
$\eta$ Параметр Сдвиговая вязкость (сопротивление жидкости сдвиговому потоку)
$\zeta$ Параметр Параметр активного напряжения (сила внутреннего активного привода)
$v$ Параметр Параметр выравнивания потока (тенденция нематиков переориентироваться под действием сдвига)
$\gamma$ Параметр Вращательная вязкость (сопротивление вращению поля директора)
$A$ Параметр Число активности (безразмерная мера активной силы относительно размера системы)
$R$ Параметр Отношение вязкостей (отношение вращательной к сдвиговой вязкости, $\gamma/\eta$)
$S$ Параметр Знак активного напряжения ($\pm 1$ для растягивающих/сжимающих систем)
$\epsilon$ Параметр Размер ядра дефекта (характерный размер ядер топологических дефектов)

Определение проблемы и ограничения

Формулировка основной проблемы и дилемма

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

Отправная точка (входное состояние/текущее состояние) для данной работы — активная нематическая турбулентность, в частности, в режиме "без дефектов". Это состояние без дефектов достигается либо путем явного построения модели (модель на основе директора, где поле директора имеет фиксированный модуль, исключающий образование дефектов), либо в условиях, когда энергетические затраты на ядра дефектов непомерно высоки (в модели Q-тензора с малым размером ядра дефекта). Система характеризуется такими параметрами, как параметр активного напряжения $\zeta$ и параметр выравнивания потока $v$, который описывает тенденцию жидкокристаллического материала переориентироваться под действием сдвига.

Желаемый конечный результат (выходное состояние/целевое состояние) — всестороннее понимание того, как параметр выравнивания потока $v$ и подавление топологических дефектов влияют на пространственно-временную структуру и динамику активных нематических потоков. Авторы стремятся выявить и охарактеризовать новые "остановленные паттерны", которые возникают в активной нематической турбулентности без дефектов, отличающиеся от хаотических потоков, обычно наблюдаемых при наличии дефектов. В конечном итоге цель состоит в том, чтобы прояснить лежащий в основе механизм этой "динамической остановки", в частности, как она возникает из эмерджентной эффективной топологии нематических доменных стенок и их правил связности.

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

Болезненный компромисс или дилемма, которая удерживала предыдущих исследователей, пытавшихся решить эту конкретную проблему, заключается в самой природе активной турбулентности. Активная турбулентность часто характеризуется устойчивыми, хаотическими, доминирующими вихрями потоками, при этом топологические дефекты часто выступают в качестве ключевых драйверов и организаторов этого хаоса. Дилемма заключается в том, что, хотя дефекты являются источниками беспорядка, они также обеспечивают устойчивую турбулентность в определенных режимах. Работа показывает, что в растягивающих стержнеобразных нематиках отсутствие топологических дефектов, вместо того чтобы приводить к более простой форме турбулентности, приводит к "динамической остановке". Это подразумевает контринтуитивный компромисс: устранение источника хаоса (дефектов) не обязательно приводит к более упорядоченному турбулентному состоянию, а скорее к совершенно другому, остановленному и менее динамичному режиму. Это бросает вызов общепринятому мнению о том, что дефекты являются исключительно разрушительными, и подчеркивает их роль в поддержании турбулентной динамики.

Ограничения и режимы отказа

Проблема понимания динамической остановки в активной нематической турбулентности становится невероятно сложной из-за нескольких суровых, реалистичных ограничений:

  • Физические/материальные ограничения:

    • Требование отсутствия дефектов: Основное внимание уделяется активной нематической турбулентности "без дефектов" — это значительное ограничение. Топологические дефекты повсеместны в большинстве активных нематических систем и часто являются основными драйверами хаотических потоков. Экспериментальная реализация и поддержание истинно бездефектной активной нематической системы чрезвычайно сложна. Работа решает эту проблему, используя модель на основе директора, которая запрещает дефекты по конструкции (фиксированный модуль $|n|=1$), и модель Q-тензора, где дефекты не образуются из-за высокой энергетической стоимости, связанной с их ядрами.
    • Чувствительность к параметрам: Возникновение динамической остановки сильно зависит от взаимодействия свойств материала, в частности, параметра выравнивания потока $v$ и параметра активного напряжения $\zeta$. Остановленное состояние наблюдается только в режиме растягивающего выравнивания ($\zeta v < 0$ для стержнеобразных нематиков с $v \le 0$). Это означает, что явление не является универсальным для всех активных нематических систем, а ограничено определенным, узким пространством параметров.
    • Низкое число Рейнольдса: Теоретическая основа и симуляции строго работают при "исчезающем числе Рейнольдса", что подразумевает пренебрежимо малые инерционные эффекты. Хотя это упрощает гидродинамику, это ограничивает прямую применимость результатов к системам, где инерция играет значительную роль, потенциально изменяя паттерны потока и стабильность.
    • Пороговый размер ядра дефекта: В более общей модели Q-тензора отношение размера ядра дефекта $\epsilon$ к активной длине $l_a$ является критическим управляющим параметром. Динамическая остановка сохраняется только при $\epsilon \ll l_a$, что означает энергетическое подавление нуклеации дефектов. Если $\epsilon$ превышает критический порог (приблизительно $0.25 l_a$ для фиксированного $Sv = -1$), происходит нуклеация дефектов, что приводит к разрушению остановленного состояния и переходу к турбулентности с дефектами. Это накладывает строгое физическое условие для наблюдения явления.

  • Вычислительные ограничения:

    • Крупномасштабные симуляции: Наблюдение эмерджентных паттернов и долгосрочной динамики динамической остановки требует "крупномасштабных численных симуляций". Это означает существенные вычислительные затраты и время, особенно при исследовании различных режимов параметров и размеров систем.
    • Численная жесткость модели Q-тензора: Симуляция полной модели Q-тензора, особенно в пределе малого размера ядра дефекта ($\epsilon \to 0$), становится "все более жесткой и вычислительно дорогой". Эта жесткость возникает из-за необходимости разрешения резких градиентов и малых пространственных масштабов, связанных с потенциальными ядрами дефектов, что делает более простую модель на основе директора более практичной альтернативой в этом пределе.
    • Сложная гибридная численная схема: Симуляции используют сложную гибридную численную схему, сочетающую псевдоспектральный метод для баланса импульса и метод конечных элементов для эволюции поля угла. Эта сложная настройка, включая такие методы, как характеристики-Галеркин для конвективных членов, неявное рассмотрение лапласианов и Адамса-Башфорта для временной эволюции, необходима для обеспечения численной стабильности и точности в разнообразных физических явлениях.
    • Ограничения размера системы: Хотя авторы провели симуляции при "еще больших размерах систем" для проверки устойчивости остановки, они "требовали значительно больших вычислительных усилий". Это указывает на то, что вычислительные ресурсы накладывают практическое ограничение на масштаб, в котором эти явления могут быть изучены, потенциально влияя на наблюдение очень дальнодействующих корреляций или стабильности остановленных состояний в истинно макроскопических системах.

  • Ограничения, основанные на данных/экспериментальные ограничения:

    • Отсутствие экспериментального подтверждения: Основным ограничением является то, что "экспериментальная проверка наших предсказаний для активных нематиков без дефектов представляет собой вызов", поскольку "динамическая остановка — наше основное предсказание — еще не была реализована" экспериментально. Хотя некоторые элементы, такие как псевдодефекты и лабиринтообразные паттерны, были наблюдаемы, устойчивое, охватывающее всю систему остановленное состояние остается неуловимым в экспериментах. Это требует поиска новых активных материалов, где стоимость нуклеации дефектов значительно выше, чем в в настоящее время изучаемых системах.
    • Сложность контроля дефектов: Способность точно контролировать или подавлять топологические дефекты в экспериментальных системах активных нематиков является значительным препятствием. Большинство экспериментальных установок по своей природе производят дефекты, что затрудняет изоляцию и изучение режима без дефектов, предсказанного моделью.

Почему такой подход

Неизбежность выбора

Выбор авторами модели на основе директора без дефектов был не просто предпочтением, а неизбежностью, продиктованной самой проблемой, которую они поставили перед собой: понимание "динамической остановки в активной нематической турбулентности" в отсутствие топологических дефектов. Традиционные "SOTA" (State-Of-The-Art, передовые) методы, такие как стандартные модели Q-тензора, разработаны для учета нуклеации и динамики топологических дефектов, которые повсеместны в большинстве активных нематических систем. Однако основная новизна данной работы заключается в исследовании режима, где эти дефекты либо энергетически подавляются, либо полностью запрещены по конструкции.

В тот момент, когда авторы осознали, что традиционные методы недостаточны, они решили исследовать состояние без дефектов. Их выбранная модель, обобщение предыдущей минимальной модели [33], явно накладывает фиксированный модуль $|n|=1$ для поля нематического директора. Это ограничение, как указано в Разделе II, "исключает генерацию топологических дефектов [33,38]". Этот преднамеренный выбор конструкции был необходим, поскольку явление динамической остановки, характеризующееся стабильными, древовидными паттернами доменных стенок, принципиально изменяется или даже предотвращается присутствием подвижных топологических дефектов. Без модели, которая гарантирует среду без дефектов, изоляция и характеристика этого нового остановленного состояния были бы невозможны. Модель Q-тензора, хотя и более общая, ввела бы дефекты, затеняя специфическую динамику, которую авторы желали выявить.

Сравнительное превосходство

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

  1. Раскрытие динамической остановки: Наиболее значительным качественным преимуществом является ее способность выявлять само состояние "динамической остановки". В растягивающих нематиках с выравниванием ($S\nu < -1$) модель показывает, что нематические доменные стенки "сильно стабилизируются потоком, который они увлекают", приводя к "заполняющему пространство древовидному паттерну" [Рис. 1(d)-1(f)]. Это остановленное состояние является новым открытием, отличным от хаотической, дефектной турбулентности, обычно наблюдаемой и изучаемой другими методами.
  2. Обработка высокоразмерного шума и долгосрочной стабильности: Остановленный режим демонстрирует "гораздо более слабые" крупномасштабные хаотические потоки и "гораздо более длительные времена корреляции, особенно поля нематического тензора в больших масштабах [Рис. 1(i), синий]". Это указывает на качественное превосходство в захвате стабильных, долгоживущих структур и динамики, которые часто подавляются высокоразмерным шумом и короткими временами корреляции в системах с дефектами. Модель эффективно снижает "шум" от динамики дефектов, позволяя четко проявиться механизму формирования паттерна.
  3. Вычислительная эффективность для целевого режима: Хотя это и не является основной целью, статья отмечает практическое преимущество. Обсуждая модель Q-тензора в Разделе V.A, авторы заявляют: "Следует отметить, что симуляции Q-тензора становятся все более жесткими и вычислительно дорогими в этом пределе [малый размер ядра дефекта], что делает модель $\theta$ более эффективной и практичной альтернативой". Для режима без дефектов (который соответствует очень малому размеру ядра дефекта) модель на основе директора имеет подавляющее преимущество по вычислительным затратам и стабильности, позволяя проводить более обширные симуляции остановленного состояния.

Соответствие ограничениям

Выбранная модель на основе директора идеально соответствует неявным ограничениям изучения активной нематической турбулентности без сложностей, вносимых топологическими дефектами. "Слияние" суровых требований проблемы и уникальных свойств решения очевидно в нескольких аспектах:

  1. Принудительное состояние без дефектов: Важнейшим ограничением является отсутствие топологических дефектов. Формулировка модели, фиксируя модуль поля директора на $|n|=1$, по своей сути предотвращает нуклеацию дефектов [33,38]. Это прямое и идеальное соответствие, гарантирующее, что любые наблюдаемые явления обусловлены исключительно взаимодействием активности, выравнивания потока и нематической упругости, а не динамикой дефектов.
  2. Фокус на выравнивании потока: Цель работы — исследовать, как "параметр выравнивания потока $\nu$" влияет на активную нематическую турбулентность. Модель жидкокристаллического материала Эриксена-Лесли, которая лежит в основе их подхода, естественным образом включает этот параметр, позволяя напрямую изучать его влияние на пространственно-временную структуру потоков и возникновение динамической остановки.
  3. Минималистский подход для ясности: Устраняя сложности динамики дефектов, модель придерживается минималистского подхода, позволяя более четко идентифицировать лежащие в основе механизмы формирования паттернов и остановки. Это соответствует цели интуитивно объяснить сложные явления читателю с нулевой базой знаний. Модель упрощает систему, чтобы выделить принципиально иные поведения, которые могут возникать при отсутствии дефектов.

Отклонение альтернатив

Работа неявно и явно отклоняет альтернативные подходы, в частности, полную модель Q-тензора, как основной инструмент для их первоначального исследования активных нематиков без дефектов. Причины этого отклонения ясны и многогранны:

  1. Вычислительные затраты и жесткость: Как упомянуто в Разделе V.A, когда размер ядра дефекта $\epsilon$ очень мал (приближаясь к пределу отсутствия дефектов), "симуляции Q-тензора становятся все более жесткими и вычислительно дорогими". Для режима, где дефекты энергетически подавляются или отсутствуют, модель на основе директора $\theta$ является "более эффективной и практичной альтернативой". Это веская практическая причина не использовать модель Q-тензора в качестве основного подхода для данной конкретной проблемы.
  2. Смешивающие эффекты дефектов: Основная задача заключалась в изучении явлений в отсутствие топологических дефектов. Модели Q-тензора по своей природе допускают нуклеацию и динамику дефектов. Использование такой модели ввело бы дополнительный уровень сложности, затрудняя изоляцию специфических механизмов, ведущих к динамической остановке в среде без дефектов. Авторы явно заявляют, что "большинство экспериментов и симуляций в активных нематиках проводились в режимах с дефектами" (Раздел VII), подразумевая, что эти предыдущие подходы не подходили для их нового исследования без дефектов.
  3. Валидация, а не замена: Модель Q-тензора вводится позже в Разделе V, не как замена, а для того, чтобы "очертить режим, в котором наши выводы для систем без дефектов остаются применимыми, и (2) исследовать переход к более привычному поведению с дефектами, которое возникает, когда эти условия больше не выполняются". Это демонстрирует, что модель Q-тензора служит ценным инструментом для контекстуализации их выводов для систем без дефектов и исследования переходов, но не является идеальным инструментом для первоначального, сфокусированного исследования самого остановленного состояния без дефектов. Работа показывает, что модель $\theta$ сходится к модели Q-тензора в пределе отсутствия дефектов, подтверждая справедливость более простой модели для данного конкретного режима. Способность отделить эффекты дефектов от других параметров имеет решающее значение для выводов работы.
FIG. 1. Strong and arrested regimes of active nematic turbulence. Snapshots from simulations of defect-free active nematic turbulence in contractile [panels (a)–(c)] and extensile [panels (d)–(f)] flow-aligning systems. Parameter values were set to R = 1, ν = −1.1, and A = 3.2 × 105. Top panels (a) and (d) show the flow field; black curves are streamlines, and the color indicates the speed (see Movies S1 and S2 of the Supplemental Material [35]). Middle panels (b) and (e) show the Frank free energy density ∼|∇θ|2, with high-intensity lines corresponding to nematic domain walls (see Movies S1 and S2 of the Supplemental Material [35]). Bottom panels (c) and (f) are zooms highlighting the type of nematic distortion as well as the interplay between nematic walls and flows. The gray-scale background is the line integral convolution representation of the director field n. Magenta and cyan intensities, respectively, represent splay (∇· n)2 and bend |∇× n|2 contributions to the Frank energy density. The black arrows represent the flow field v, which localizes along the nematic walls in the arrested regime. Black circles indicate stagnation points of the flow. White scale bar represents the selected wavelength λi. (g)–(i) Spectra characterizing fully developed active nematic turbulence (see details in Appendix B). The lines in panels (h) and (i) represent a smoothed (Gaussian) interpolation of the computed data points. We compare the flow-aligning contractile [red, as in panels (a)–(c) and Movie S1 of the Supplemental Material [35]] and extensile [blue, as in panels (d)–(f) and Movie S2 of the Supplemental Material [35]] cases with the ν = 0 case (black, as in Fig. S1 and Movie S3 of the Supplemental Material [35]), for which contractile and extensile stresses are equivalent up to a rotation [38,42,43]. (g) Velocity power spectrum on a log-log scale, showing (1) the universal low-q scaling law and (2) the distinct organization of flows across scales in the different cases. The wider scaling regime in the contractile case captures the strong large-scale jets [see panel (a)]. The peak in the extensile case is representative of wall streams [see panel (d)]. (h) Frank energy spectrum, showing that (1) the selected wavelength (peak position) depends on ν but not on the sign of active stress and (2) the peak width depends on the sign of active stress when ν ̸= 0. (i) Spectrum of correlation times associated with the flow v (light colored points and lines) and the nematic tensor ˆqαβ (darker points and lines). This log-log plot reveals strong differences in decay times between the regimes, as well as the differences between the flow and nematic tensor within a regime. Correlation times are extracted from exponential fits to the corresponding space-time autocorrelation functions in Fourier space (see Appendix B)

Математический и логический механизм

Основное уравнение

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

Безразмеренное уравнение баланса импульса (полученное из ротора уравнения Навье-Стокса при исчезающем числе Рейнольдса) имеет вид:
$$ \nabla^4\psi = -S\partial_\alpha\partial_\beta\tilde{g}_{\alpha\beta} + R\nu\partial_\alpha\partial_\beta(\tilde{g}_{\alpha\beta}h_\parallel - \tilde{g}_{\alpha\beta}h_\perp) \quad \text{(6)} $$

Безразмеренное уравнение динамики директора (описывающее, как ориентация нематика изменяется со временем) имеет вид:
$$ \partial_t\theta - \epsilon_{\alpha\beta}\partial_\alpha\psi \partial_\beta\theta + \frac{1}{A}\nabla^2\theta = h_\parallel + \nu\tilde{g}_{\alpha\beta}\partial_\alpha\partial_\beta\psi \quad \text{(8)} $$
Эти два уравнения, (6) и (8), составляют математический двигатель модели активной нематической турбулентности без дефектов. Они связаны, поскольку течение жидкости (представленное $\psi$) влияет на нематический директор ($\theta$), а конфигурация нематического директора влияет на течение жидкости.

Потермный анализ

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

Уравнение (6): Баланс импульса

$$ \nabla^4\psi = -S\partial_\alpha\partial_\beta\tilde{g}_{\alpha\beta} + R\nu\partial_\alpha\partial_\beta(\tilde{g}_{\alpha\beta}h_\parallel - \tilde{g}_{\alpha\beta}h_\perp) $$

  • $\nabla^4\psi$:

    • Математическое определение: Это бигармонический оператор, действующий на функцию тока $\psi$. Он эквивалентен $\nabla^2(\nabla^2\psi)$. Поскольку $\omega = -\nabla^2\psi$ (вихрь), этот член фактически является $-\nabla^2\omega$.
    • Физическая/логическая роль: В гидродинамике этот член представляет собой вязкую диссипацию вихря. Он действует как механизм затухания, сопротивляясь изменениям в потоке жидкости и стремясь сгладить градиенты скорости. Использование $\nabla^4$ (вместо просто $\nabla^2$) возникает из-за двойного взятия ротора уравнения баланса импульса для исключения давления и выражения всего через функцию тока.
    • Почему $\nabla^4$: Этот оператор естественным образом возникает при взятии ротора уравнения Стокса (баланс импульса при низком числе Рейнольдса), когда он выражен через функцию тока. Он обеспечивает несжимаемость потока жидкости и учитывает вязкие силы.

  • $S$:

    • Математическое определение: Безразмерный параметр, $S = \text{sign}(\zeta)$, где $\zeta$ — параметр активного напряжения.
    • Физическая/логическая роль: Этот параметр определяет природу активного напряжения. $S = +1$ для растягивающих напряжений (например, стержнеобразные компоненты, отталкивающиеся наружу) и $S = -1$ для сжимающих напряжений (например, моторные белки, тянущие внутрь). Это ключевой управляющий параметр поведения системы, влияющий на то, проявляет ли система сильную турбулентность или динамическую остановку.
    • Почему знак: Знак отражает фундаментальное различие в том, как активные компоненты генерируют напряжение, что имеет глубокие последствия для результирующих паттернов потока.

  • $\partial_\alpha\partial_\beta\tilde{g}_{\alpha\beta}$:

    • Математическое определение: Это дивергенция тензора ориентации нематика $\tilde{g}_{\alpha\beta}$. Тензор $\tilde{g}_{\alpha\beta}$ определяется как $\tilde{g}_{\alpha\beta} = n_\alpha n_\beta - \frac{1}{2}\delta_{\alpha\beta}$, где $n_\alpha$ — компонент нематического директора $\mathbf{n} = (\cos\theta, \sin\theta)$, а $\delta_{\alpha\beta}$ — символ Кронекера.
    • Физическая/логическая роль: Этот член представляет собой активное напряжение, генерируемое нематическими компонентами. Активные нематики преобразуют химическую энергию в механическую работу, генерируя напряжения, которые приводят в движение поток жидкости. Дивергенция $\tilde{g}_{\alpha\beta}$ отражает распределение этих напряжений и то, как они приводят в движение жидкость.
    • Почему дивергенция: Дивергенция тензора напряжений (например, $\sigma^A$) представляет собой плотность силы, действующей со стороны этого напряжения на жидкость. Это стандартный способ связи напряжения с импульсом жидкости.

  • $R$:

    • Математическое определение: Безразмерное отношение вязкостей, $R = \gamma/\eta$.
    • Физическая/логическая роль: Этот параметр сравнивает вращательную вязкость ($\gamma$) нематического директора со сдвиговой вязкостью ($\eta$) жидкости. Он влияет на то, насколько легко нематический директор переориентируется в ответ на поток и упругие моменты.
    • Почему отношение: Это отношение, поскольку оно сравнивает два различных типа вязкостного сопротивления в системе, что имеет решающее значение для понимания баланса между вращением директора и течением жидкости.

  • $\nu$:

    • Математическое определение: Безразмерный параметр выравнивания потока.
    • Физическая/логическая роль: Этот параметр характеризует тенденцию директора жидкокристаллического материала переориентироваться под действием сдвигового потока. Если $\nu < 0$, нематик имеет тенденцию выравниваться по потоку (режим выравнивания потока). Если $\nu > 0$, он имеет тенденцию к вращению (режим вращения). Этот параметр является центральным для выводов работы о динамической остановке.
    • Почему коэффициент: Это коэффициент, который масштабирует силу эффекта выравнивания потока, являющегося ключевым свойством жидкокристаллических материалов.

  • $\partial_\alpha\partial_\beta(\tilde{g}_{\alpha\beta}h_\parallel - \tilde{g}_{\alpha\beta}h_\perp)$:

    • Математическое определение: Этот член представляет собой дивергенцию напряжений, возникающих из-за выравнивания потока и напряжений Эриксена. $h_\parallel$ и $h_\perp$ — компоненты поля ориентации $h_\alpha$, определенного в Уравнении (7) как $h_\parallel = \frac{1}{A}\nabla^2\theta$ и $h_\perp = \nu(\sin 2\theta d_1\psi + \cos 2\theta d_2\psi)$, где $d_1 = \frac{1}{2}(\partial_x^2 - \partial_y^2)$ и $d_2 = \partial_x\partial_y$. Полное $h_\alpha$ — это упругий момент, действующий на директора.
    • Физическая/логическая роль: Этот член захватывает упругие напряжения (напряжение Эриксена) и напряжения, связанные с выравниванием потока. Упругие напряжения возникают из-за искажений поля нематического директора, в то время как выравнивание потока описывает, как нематический директор переориентируется в ответ на поле деформации жидкости (связанное с $\partial_\alpha\partial_\beta\psi$). Эти напряжения возвращаются в баланс импульса жидкости, влияя на поле потока.
    • Почему разность: Конкретная комбинация $\tilde{g}_{\alpha\beta}h_\parallel - \tilde{g}_{\alpha\beta}h_\perp$ возникает из симметричной части девиаторного тензора напряжений (Уравнение 4), который включает вклады от упругих моментов и выравнивания потока. Дивергенция этого тензора представляет собой плотность силы, действующей на жидкость.

Уравнение (8): Динамика директора

$$ \partial_t\theta - \epsilon_{\alpha\beta}\partial_\alpha\psi \partial_\beta\theta + \frac{1}{A}\nabla^2\theta = h_\parallel + \nu\tilde{g}_{\alpha\beta}\partial_\alpha\partial_\beta\psi $$

  • $\partial_t\theta$:

    • Математическое определение: Частная производная нематического угла $\theta$ по времени $t$.
    • Физическая/логическая роль: Это член временной эволюции, представляющий скорость изменения ориентации нематического директора в фиксированной точке пространства. Это ядро динамики.
    • Почему производная: Это фундаментальный компонент любого динамического уравнения, описывающий скорость изменения переменной состояния.

  • $\epsilon_{\alpha\beta}\partial_\alpha\psi \partial_\beta\theta$:

    • Математическое определение: Это конвективный член, представляющий адвекцию нематического угла $\theta$ течением жидкости. $\epsilon_{\alpha\beta}$ — тензор Леви-Чивиты (полностью антисимметричный), а $\partial_\alpha\psi$ представляет компоненты скорости жидкости $\mathbf{v}$.
    • Физическая/логическая роль: Этот член описывает, как нематический директор переносится течением жидкости. Если элемент жидкости движется, он уносит с собой свою нематическую ориентацию. Это критический нелинейный член, который связывает динамику директора с скоростью жидкости.
    • Почему $\epsilon_{\alpha\beta}$ и умножение: Тензор Леви-Чивиты используется для формирования ротора функции тока, который дает компоненты скорости. Член $\mathbf{v} \cdot \nabla\theta$ представляет материальную производную, отражающую изменения, обусловленные как локальной временной эволюцией, так и переносом потоком.

  • $\frac{1}{A}\nabla^2\theta$:

    • Математическое определение: Лапласиан нематического угла $\theta$, масштабированный обратным числом активности $A$.
    • Физическая/логическая роль: Этот член представляет собой упругий момент, сопротивляющийся пространственным вариациям (искажениям) поля нематического директора. Он стремится сгладить резкие изменения ориентации, действуя как восстанавливающая сила. Число активности $A = L^2/l_a^2$ сравнивает размер системы $L$ с активной длиной $l_a = \sqrt{K/(|\zeta|R)}$, которая является балансом между активными и упругими нематическими напряжениями. Большее $A$ означает, что упругие силы относительно слабее по сравнению с активными силами.
    • Почему $\nabla^2$: Лапласиан — стандартный математический оператор для диффузии или упругих сил, стремящихся минимизировать пространственные градиенты.

  • $h_\parallel$:

    • Математическое определение: Параллельная компонента поля ориентации $h_\alpha$, которая связана с свободной энергией Франка. В частности, $h_\parallel = \frac{1}{A}\nabla^2\theta$ (из Уравнения (7) и контекста Уравнения (8)).
    • Физическая/логическая роль: Этот член представляет собой упругий момент, действующий на директора, направляя его к состоянию минимального упругого искажения. Это та же самая упругая восстанавливающая сила, что и член $\frac{1}{A}\nabla^2\theta$ в левой части, но здесь она явно показана как член источника. Формулировка статьи в (8) эффективно объединяет члены упругого момента.
    • Почему сложение: Члены в правой части являются источниками или стоками вращения директора, внося вклад в общее изменение $\theta$.

  • $\nu\tilde{g}_{\alpha\beta}\partial_\alpha\partial_\beta\psi$:

    • Математическое определение: Член выравнивания потока, масштабированный параметром выравнивания потока $\nu$ и включающий тензор ориентации нематика $\tilde{g}_{\alpha\beta}$ и вторые производные функции тока $\psi$.
    • Физическая/логическая роль: Этот член описывает переориентацию нематического директора из-за растягивающего потока (сдвига). Он отражает, как поле деформации жидкости (связанное с $\partial_\alpha\partial_\beta\psi$) вызывает выравнивание или вращение стержней нематика, в зависимости от знака $\nu$. Это критический связующий член.
    • Почему умножение: Это произведение параметра выравнивания потока, ориентации нематика и тензора скорости деформации (полученного из $\partial_\alpha\partial_\beta\psi$), отражающее, как эти факторы объединяются для создания переориентирующего момента.

Ключевые безразмерные параметры:

  • $L$: Размер системы (пространственный масштаб).
  • $\tau_a = \eta/|\zeta|$: Активное временное масштабирование.
  • $K$: Упругая постоянная Франка.
  • $\eta$: Сдвиговая вязкость.
  • $\gamma$: Вращательная вязкость.
  • $\zeta$: Параметр активного напряжения.
  • $A = L^2/l_a^2 = RL^2|\zeta|/K$: Число активности. Сравнивает размер системы с активной длиной $l_a$. Большое $A$ указывает на сильное активное воздействие по сравнению с упругим сопротивлением.
  • $R = \gamma/\eta$: Отношение вязкостей.
  • $\nu$: Параметр выравнивания потока.

Пошаговый поток

Представьте себе одну абстрактную точку данных, представляющую состояние активной нематической системы в определенном месте и времени. Это состояние в основном определяется локальным углом нематического директора $\theta$ и функцией тока жидкости $\psi$. Система итеративно эволюционирует посредством численного моделирования, подобно механической сборочной линии, обрабатывающей эти точки данных.

  1. Входные данные начального состояния: В данный момент времени $n$ у нас есть текущее поле угла нематика $\theta^n$ и поле функции тока $\psi^n$ по всей пространственной сетке.

  2. Расчет баланса импульса (Уравнение 6 / A1):

    • Система сначала фокусируется на обновлении потока жидкости. Она берет текущий угол нематика $\theta^n$ и использует его для расчета активных напряжений (член $S(d_1 \sin 2\theta^n + d_2 \cos 2\theta^n)$) и напряжений выравнивания потока/Эриксена (член $-\frac{R\nu}{A}(d_1 \cos 2\theta^n \nabla^2\theta^n - d_2 \sin 2\theta^n \nabla^2\theta^n)$, который является частью правой части (A1)). Эти члены действуют как "функции возбуждения", которые приводят в движение жидкость.
    • Эти функции возбуждения затем подаются в уравнение баланса импульса, которое является бигармоническим уравнением для $\psi$. Это уравнение решается с использованием псевдоспектрального метода. Это означает, что пространственные производные вычисляются очень эффективно в Фурье-пространстве.
    • Критически важно, что член $G(\theta^n, \psi^n)$ в левой части (A1) (который является частью напряжения выравнивания потока) зависит от самого $\psi$. Для обработки этого используется итерация неподвижной точки: начальное предположение для $\psi$ (часто из предыдущего временного шага) используется для расчета $G$, затем уравнение решается для нового $\psi$. Этот процесс повторяется до тех пор, пока рассчитанное $\psi$ не сойдется к стабильному значению для текущего временного шага, обеспечивая самосогласованность между потоком и напряжениями, которые он генерирует.
    • Как только $\psi^n$ определено, вычисляются вспомогательные поля потока, такие как компоненты скорости ($v_x, v_y$), вихрь ($\omega$) и вращения выравнивания потока ($C^n$) из $\psi^n$.

  3. Расчет динамики директора (Уравнение 8 / A4):

    • Затем система переходит к обновлению нематического директора. Она берет недавно вычисленные поля потока ($v_x, v_y$) и текущий угол нематика $\theta^n$.
    • Уравнение динамики директора (8) описывает, как $\theta$ изменяется. Конвективный член ($\epsilon_{\alpha\beta}\partial_\alpha\psi \partial_\beta\theta$) определяет, как директор адвектируется течением жидкости. Упругий член ($\frac{1}{A}\nabla^2\theta$) сглаживает искажения директора. Члены в правой части ($h_\parallel + \nu\tilde{g}_{\alpha\beta}\partial_\alpha\partial_\beta\psi$) представляют упругие моменты и переориентацию из-за растягивающего потока.
    • Это уравнение решается с использованием метода конечных элементов. Конвективный член обрабатывается отдельно с использованием метода Галеркина на характеристиках для стабильности, фактически отслеживая путь элементов жидкости. Упругий (лапласианский) член обрабатывается неявно, то есть он решается одновременно с новым $\theta^{n+1}$ для обеспечения численной стабильности, особенно для жестких задач.
    • Результатом является обновленное поле угла нематика $\theta^{n+1}$.

  4. Переход к следующему временному шагу: Затем система увеличивает время с $n$ до $n+1$, и весь процесс повторяется, используя $\theta^{n+1}$ и $\psi^{n+1}$ в качестве новых входных данных.

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

Оптимизационная динамика

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

  1. Численная сходимость на каждом временном шаге:

    • Итерация неподвижной точки для функции тока: В пределах каждого временного шага уравнение баланса импульса (A1) для функции тока $\psi$ решается итеративно методом неподвижной точки. Это означает, что начальное предположение для $\psi$ используется для расчета нелинейных членов, затем вычисляется новое $\psi$. Это новое $\psi$ затем становится входными данными для следующей итерации, и процесс повторяется до тех пор, пока разница между последовательными решениями $\psi$ не станет меньше малого вычислительного допуска (например, $10^{-8}$). Это гарантирует, что поле потока жидкости самосогласовано с напряжениями, генерируемыми нематическим директором в данный момент. "Ландшафт потерь" здесь — это остаточная ошибка уравнения, а итерации направлены на поиск его "минимума" (нулевого остатка).
    • Неявное рассмотрение для динамики директора: Для уравнения динамики директора (A4) лапласианский (упругий) член обрабатывается неявно. Это означает, что значение $\theta$ на следующем временном шаге, $\theta^{n+1}$, решается напрямую, а не явно вычисляется из $\theta^n$. Этот подход повышает численную стабильность, особенно при работе с жесткими членами (такими как упругие силы, которые могут вызывать быстрые изменения, если их обрабатывать неосторожно). Он эффективно "сходится" к стабильному $\theta^{n+1}$, которое удовлетворяет уравнению при текущем потоке.

  2. Физическая эволюция системы и "обучение":

    • Итеративное обновление состояния: Вся симуляция прогрессирует путем итеративного обновления полей $\theta$ и $\psi$ от одного временного шага к другому. Этот непрерывный цикл обратной связи между потоком и динамикой директора движет эволюцией системы.
    • Ландшафт потерь (неявный): Хотя явно не определена как функция потерь, динамика системы может рассматриваться как навигация по сложному "энергетическому ландшафту" (например, свободная энергия Франка, объемная свободная энергия). Система естественным образом эволюционирует к состояниям, минимизирующим эти энергии, или, в активных системах, достигает динамически устойчивых неравновесных состояний.
    • Градиенты и движущие силы: Различные члены в уравнениях действуют как "градиенты" или движущие силы. Например, член активного напряжения ($S\partial_\alpha\partial_\beta\tilde{g}_{\alpha\beta}$) приводит в движение жидкость, в то время как упругие члены ($\nabla^2\psi$ и $\nabla^2\theta$) действуют как восстанавливающие силы, стремясь сгладить искажения. Член выравнивания потока ($\nu\tilde{g}_{\alpha\beta}\partial_\alpha\partial_\beta\psi$) определяет, как директор переориентируется в ответ на поток, формируя реакцию системы.
    • Сходимость к динамической остановке: Центральный вывод работы заключается в том, что для определенных режимов параметров (например, растягивающих нематиков с выравниванием и $S\nu < -1$) система "сходится" к динамически остановленному состоянию. Это не оптимизация в смысле поиска глобального минимума, а скорее система, оседающая в стабильном, неравновесном паттерне (древовидные сети доменных стенок), где хаотические потоки подавляются, а динамика становится очень медленной, напоминая явления старения в стеклообразных системах. Система "обучается" формировать эти стабильные паттерны посредством непрерывного взаимодействия активного возбуждения, вязкой диссипации, упругих моментов и выравнивания потока. "Градиенты" эффективно направляют систему в эти остановленные конфигурации.
    • Роль шума: Небольшая амплитуда гауссова белого шума добавляется к динамике директора (Приложение A.4) для ускорения эволюции к полностью развитой активной турбулентности. Этот шум действует как возмущение, помогая системе исследовать "ландшафт" и избегать потенциальных локальных минимумов или метастабильных состояний, обеспечивая достижение характерных турбулентных или остановленных режимов.

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

FIG. 4. Motifs of arrested wall networks and experimental snapshots of microtubule-kinesin active nematics. In all diagrams, lines, nodes, and colors are as defined in Fig. 3. (a)–(d) Basic network motifs. The anchoring motif [panel (a)] is made of an endpoint and branchpoint that meet head-on, with the endpoint trapped between the two outgoing walls of the branchpoint. In the motif depicted in panel (b), the endpoint meets the branchpoint from one of its sides, i.e., between the incoming wall and an outgoing wall. The dashed line traces a weak distortion, indicating that the wall associated with the endpoint tends to align its direction with the outgoing wall on the opposite side of the branchpoint. The motifs shown in panels (c) and (d) involve a single pseudodefect interacting with a bare wall. These, along with panel (b), do not follow the tendency to have strictly antiparallel walls. (e) Composite motif schematic (left) and one formed spontaneously in a simulation (right). The stream plot on the right represents the flow, with black indicating maximal |v| and full transparency indicating |v| = 0. The gray background is the line-integral-convolution representation of the nematic director n. Parameter values and color legend for splay and bend distortions are as in Fig. 1(f). (f), (g) Raw fluorescence images from experiments (left panels) and overlaid schematic drawings (right panels) depicting domain walls, pseudodefects, and actual ±1/2 defects in white. (f) Taken from a movie in Ref. [50] (courtesy of Guillaume Duclos), which shows the evolution of the microtubule-based nematic following the bending instability of the aligned state. (g) Taken from a movie in Ref. [51] (courtesy of Pau Guillamat), which shows a turbulent transient with all types of pseudodefects and actual nematic defects. Note how walls may also originate from true +1/2 defects and be absorbed by true −1/2 defects FIG. 5. Ageing of the arrested wall network. Skeleton of the domain walls (black) with startpoints, branchpoints, and endpoints (green, blue, and red triangular nodes) at an early time [panel (a)] and a late time [panel (b)]. The detection of the network skeleton and its nodes is described in Appendix E (Fig. 13). (c) Evolution of the number of startpoints (green), branchpoints (blue), and endpoints (red). In the initial transient, sequential “zigzag” instabilities result in the proliferation of both branchpoints and endpoints. Once the wall pattern establishes a wavelength, the system ages slowly as some endpoints retract and annihilate with their connected branchpoint, while others extend (Movie S5 of the Supplemental Material [35]). Throughout the simulation, there are frequent transitions between branchpoints and startpoints, though the number of startpoints re- mains low. Additionally, the detection algorithm is not perfect, occasionally misidentifying endpoints or branchpoints as startpoints and vice versa. Parameter values were set to R = 1, S = 1, ν = −0.9, and A = 3.2 × 105

Результаты, ограничения и заключение

Экспериментальный дизайн и базовые уровни

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

Во-первых, использовалась модель на основе директора без дефектов, которая обобщала предыдущие работы путем включения выравнивания потока и напряжения Эриксена. Эта модель, управляемая уравнениями (6)–(8), по своей сути исключает образование топологических дефектов, обеспечивая непрерывное поле директора с фиксированным модулем. Этот архитектурный выбор был критически важен для изоляции динамики доменных стенок без искажающего влияния дефектов. Ключевые безразмерные параметры в этой модели включали число активности $A = L^2/l_a^2$, отношение вязкостей $R = \gamma/\eta$ и параметр выравнивания потока $\nu$. Для большинства симуляций были установлены $R=1$ и $A=3.2 \times 10^5$ — режим, известный проявлением крупномасштабной активной турбулентности при $\nu=0$. Эксперименты были специально разработаны для сравнения двух различных режимов:
1. Сжимающий режим с выравниванием (Sv > 1), например, $S=-1$ (сжимающее напряжение) и $\nu=-1.1$ (стержнеобразные элементы с выравниванием по потоку). Это служило базовым уровнем для "сильной крупномасштабной турбулентности".
2. Растягивающий режим с выравниванием (Sv < -1), с $S=+1$ (растягивающее напряжение) и $\nu=-1.1$. Предполагалось, что этот режим будет проявлять "остановленную турбулентность".
В качестве нейтрального базового уровня для сравнения также был включен эталонный случай с $\nu=0$ (без выравнивания), где сжимающие и растягивающие напряжения эквивалентны с точностью до поворота.

Во-вторых, использовалась неограниченная модель Q-тензора для расширения контекста их выводов и исследования перехода к динамике с дефектами. Эта модель, основанная на рамках Ландау-де Женна (модель Бериса-Эдвардса), подробно описанная в Приложении D, явно допускала нуклеацию истинных топологических дефектов $\pm 1/2$. Критическим управляющим параметром здесь было отношение размера ядра дефекта $\epsilon$ к активной длине $l_a$. Симуляции с этой моделью использовали параметры $R=1$, $S=1$, $A=10000$ и $\nu=-1$, начиная с покоящегося нематического состояния с небольшой угловой пертурбацией.

Численное интегрирование основывалось на гибридной схеме, сочетающей псевдоспектральный метод для баланса импульса (Уравнение (6)) с методом конечных элементов (FEM) для эволюции поля угла (Уравнение (8)), как описано в Приложении A. Для ускорения эволюции к полностью развитой активной турбулентности вводился небольшой амплитудный гауссов белый шум. "Жертвами" (базовыми моделями), от которых авторы стремились отличиться или победить, были в основном случай с $\nu=0$ без выравнивания и сильная крупномасштабная турбулентность, наблюдаемая в сжимающем режиме с выравниванием. Модель Q-тензора также служила строгим тестом на надежность предсказаний модели без дефектов.

Что доказывают доказательства

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

Основной механизм динамической остановки был безжалостно доказан путем сравнения растягивающего режима с выравниванием (Sv < -1) со сжимающим режимом с выравниванием (Sv > 1) и эталонным случаем $\nu=0$. В сжимающем режиме с выравниванием симуляции (Рис. 1a-1c, Видео S1) показали высокохаотические потоки с сильными, переходными крупномасштабными струями и фрагментированными, динамически реорганизующимися нематическими доменными стенками. Это представляло собой базовый уровень "сильной крупномасштабной турбулентности". В явном контрасте, растягивающий режим с выравниванием (Рис. 1d-1f, Видео S2) продемонстрировал значительно более слабые крупномасштабные потоки, с доминирующими потоками, локализованными вдоль изгибающихся доменных стенок. Критически важно, что эти доменные стенки были сильно стабилизированы, росли и ветвились в устойчивый, заполняющий пространство древовидный паттерн, который стал "заблокированным" — состояние, названное динамической остановкой.

Окончательные, неоспоримые доказательства этой остановки были получены из спектрального анализа (Рис. 1g-1i):
* Спектр мощности скорости (Рис. 1g): Хотя оба режима демонстрировали универсальное масштабирование $q^{-1}$ при низких волновых числах, сжимающий случай показал широкий диапазон масштабов, связанных с крупномасштабными струями. Растягивающий остановленный случай, однако, имел более узкий диапазон $q^{-1}$ и отчетливый пик на более коротких пространственных масштабах, подтверждая доминирование локализованных потоков стенок над крупномасштабными хаотическими потоками.
* Спектр упругой энергии (Рис. 1h): Растягивающий (остановленный) случай показал гораздо более узкий и резкий пик, указывающий на высокоорганизованный и стабильный паттерн стенок, в отличие от более широкого пика в хаотическом сжимающем режиме.
* Спектры времени корреляции (Рис. 1i): Это, возможно, было самое убедительное доказательство. Остановленный режим продемонстрировал значительно более длительные времена корреляции как для поля потока, так и для поля нематического тензора, особенно в больших масштабах. Это напрямую количественно определило снижение динамики и "запирание" паттерна, доказав, что система действительно претерпела динамическую остановку.

Дополнительные доказательства подчеркнули эмерджентную топологию этих остановленных состояний. Доменные стенки образовывали "псевдодефекты" (точки начала, точки ветвления, конечные точки), которые, хотя и не являлись истинными топологическими дефектами, несли сохраняющийся псевдозаряд (Рис. 3). Симуляции проиллюстрировали динамическую эволюцию от начальных полосатых паттернов через зигзагообразные неустойчивости и ветвление стенок к конечному заблокированному древовидному паттерну (Рис. 2). "Мотив якорения" (Рис. 4a), специфическая конфигурация псевдодефектов, был идентифицирован как особенно стабильная структура, действующая как "ловушка" для потока. Медленная релаксация числа псевдодефектов за длительные промежутки времени (Рис. 5c), напоминающая старение в стеклообразных системах, еще больше подчеркнула остановленный характер состояния. Более того, было показано, что остановленные паттерны заключают в себе "уникурсальные лабиринты" (Рис. 6b), редкий паттерн в других системах, демонстрируя уникальное явление формирования паттернов.

Надежность динамической остановки затем была подтверждена с использованием модели Q-тензора. Эта модель подтвердила, что остановка не является артефактом ограничения без дефектов, а является общим свойством, когда нуклеация дефектов энергетически подавлена (т.е. малый $\epsilon/l_a$). По мере увеличения размера ядра дефекта $\epsilon$ сверх критического порога ($\approx 0.25 l_a$) начинали нуклеироваться истинные дефекты $\pm 1/2$ (Рис. 7b, 7e, 7f), нарушая остановленные паттерны стенок и переводя систему в неупорядоченное, доминируемое вихрями турбулентное состояние. Этот переход сопровождался резким снижением общей упругой свободной энергии (Рис. 7b, 12c), поскольку дефекты эффективно растворяли высокоэнергетические доменные стенки. Симуляции также выявили, что псевдодефекты ветвления действовали как "горячие точки" для нуклеации дефектов (Видео S8), устанавливая четкую связь между режимами без дефектов и с дефектами. Количественное согласие между моделью Q-тензора и ограниченной моделью $\theta$ при малом $\epsilon/l_a$ (Рис. 11, Видео S6) еще больше укрепило справедливость более простой модели в этом пределе.

Ограничения и будущие направления

Хотя данная работа дает глубокое понимание активной нематической турбулентности, она также выявляет несколько ограничений и открывает многочисленные направления для будущих исследований.

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

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

Нелинейный механизм выбора длины волны остается открытым вопросом. Работа признает, что этот механизм, балансирующий укрупнение и изгиб/складывание стенок, является нелинейным и двумерным, и его точная зависимость от различных параметров еще не полностью понята. Необходимы более глубокие теоретические и вычислительные исследования для раскрытия этого сложного взаимодействия.

Кроме того, наблюдаемая немонотонная зависимость плотности дефектов от размера ядра (уменьшение с увеличением размера ядра в некоторых режимах) интригует, но не полностью объяснена и отложена для будущей работы. Это предполагает более сложную связь между энергетикой дефектов и формированием паттернов, чем в настоящее время понимается. Вычислительные затраты симуляций Q-тензора при очень малых размерах ядер дефектов также представляют собой практическое ограничение, делая более простую модель $\theta$ необходимой альтернативой в этом режиме.

Заглядывая вперед, возникает несколько тем для обсуждения:

  • Проектирование активных нематических систем без дефектов: Как мы можем использовать понимание выравнивания потока и размера ядра дефекта для разработки активных нематических систем, которые изначально подавляют образование дефектов, тем самым способствуя динамической остановке? Это может включать настройку свойств материала, уровней активности или геометрии ограничения.
  • Взаимодействие псевдодефектов и истинных дефектов: Переходный режим, где сосуществуют псевдодефекты и истинные топологические дефекты, богат и заслуживает дальнейшего изучения. Как истинные дефекты, которые могут действовать как точки начала или окончания стенок, взаимодействуют с сетью псевдодефектов? Можно ли контролировать это взаимодействие для переключения между остановленным и турбулентным состояниями?
  • Природа дальнодействующего топологического порядка: Уникурсальные лабиринты, наблюдаемые в остановленном состоянии, предполагают форму дальнодействующего топологического порядка. Каковы математические и физические свойства этого порядка? Как он соотносится с другими формами топологического порядка в системах конденсированных сред, и каковы его последствия для свойств материалов?
  • Аналогии со стеклообразными системами и старением: Медленная динамика релаксации остановленных сетей стенок имеет сходство с явлениями старения в стеклообразных системах. Можно ли использовать идеи из физики стекол для лучшего понимания стабильности, релаксации и реакции остановленных активных нематиков на возмущения?
  • Влияние истощения составляющих: В работе упоминается, что нематические составляющие наблюдаются истощенными у дефектов в экспериментах, что может снизить стоимость образования дефектов. Включение этого эффекта в будущие модели могло бы дать более точную картину нуклеации дефектов и ее влияния на стабильность сети стенок.
  • Обобщение на три измерения: Хотя данное исследование сосредоточено на 2D системах, расширение анализа на три измерения было бы естественным следующим шагом. Возникли бы подобные явления динамической остановки или псевдотопологические структуры в 3D активных нематиках, и чем отличались бы их свойства?

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

FIG. 8. Verification of the numerical scheme. (a) Consecutive snapshots from a simulation initiated with a small bending modulation about the uniformly aligned state: θ(r,t = 0) = ϵ cos(2πx), with ϵ = 0.01. Parameter values were set to R = 1, ν = −1.1, A = 1000, and S = 1 (extensile stress). Numerical parameters were set to N = 64, dt = 0.01, and D = 0 (zero noise). Color indicates the Frank energy, |∇θ|2, and white lines trace the director n. (b) We plot the numerical data corresponding to a horizontal slice, θ(x, y = 0), for increasing times, starting from t = 0 (lightest gray) up to t = 40τa (black) in intervals of 2.5τa. It is shown that the angle profile relaxes on the predicted 1D steady state (dashed blue), obtained for the same parameter values as explained in our recent work [38]. (c) The maximum angle is plotted as a function of time in our simulation (black). First, it is shown that the modulation grows exponentially in time with the growth rate ? matching the linear stability dispersion relation [38] (cyan). The gray inset is a log-scale plot comparing the simulation results with this explicit prediction for small perturbations. Additionally, it is demonstrated that the maximal angle in the simulation relaxes at long times to the value predicted by the 1D steady state [38] (dashed blue) FIG. 11. Agreement with the constrained model at small defect core size. For all simulations shown in Fig. 7, which share the same initial conditions and parameters except for ϵ, we plot the space- averaged deviation of the nematic tensor Q from that corresponding to the reference θ model simulation. Shaded regions indicating the standard deviation. As ϵ decreases, quantitative agreement (low deviation) persists systematically for longer integration times. A deviation of ∼0.5 corresponds to complete statistical decorrelation. See also Movie S6 of the Supplemental Material [35]

Изоморфизмы с другими областями

Структурный каркас

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

Дальние родственники

  1. Целевая область: Нейронаука (Мозговые сети)
    Связь: Основная логика данной работы, исследующая, как локальные взаимодействия между активными элементами (нематическими частицами) приводят к формированию сложных, древовидных сетей доменных стенок и псевдодефектов, которые в конечном итоге оседают в "остановленном" состоянии, является зеркальным отражением давних проблем в нейронауке, касающихся динамики мозговых сетей. В мозге отдельные нейроны и синапсы (аналогичные активным элементам) взаимодействуют локально, порождая эмерджентные функциональные цепи и сетевые паттерны. Эти сети могут проявлять состояния "динамической остановки", такие как синхронизированная, заблокированная активность, наблюдаемая при эпилептических приступах, или стабильные, устойчивые паттерны активности, связанные с консолидацией памяти. "Псевдодефекты" в нематической системе можно рассматривать как критические узлы или центры в нейронных сетях, которые определяют поток информации и общую стабильность сети, в то время как "доменные стенки" представляют собой границы между различными функциональными областями мозга или путями.

  2. Целевая область: Градостроительство и управление дорожным движением
    Связь: Исследование данной работы того, как активные жидкости посредством внутреннего привода и выравнивания потока переходят от хаотической турбулентности к "остановленному" состоянию, характеризующемуся стабильными, древовидными паттернами доменных стенок, которые направляют поток и подавляют хаос, напрямую отражает проблемы в градостроительстве. В городских транспортных системах отдельные транспортные средства (активные элементы) взаимодействуют локально, порождая сложные и часто турбулентные транспортные потоки. При определенных условиях эти потоки самоорганизуются в "пробки" или "заторы" (динамическая остановка), где формируются "доменные стенки" (очереди, заблокированные перекрестки), которые направляют или препятствуют дальнейшему движению. "Уникурсальные лабиринты", наблюдаемые в данной работе, могут быть аналогичны запутанным дорожным сетям или односторонним системам, которые, хотя и кажутся организованными, могут приводить к общесистемному застою или неэффективному потоку. "Псевдодефекты" подобны критическим перекресткам или точкам слияния, которые из-за своих локальных правил могут вызывать глобальную остановку движения.

Сценарий "Что если"

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

Прорывом стала бы возможность предсказывать и предотвращать заторы с беспрецедентной точностью. Манипулируя параметрами, аналогичными параметру выравнивания потока $v$ и параметру активного напряжения $\zeta$ (например, посредством динамического назначения полос, адаптивных систем светофоров или даже ценообразования за пробки, которое изменяет "активное напряжение"), градостроители могли бы проектировать транспортные системы, которые изначально избегают условий, ведущих к "динамической остановке". Они могли бы идентифицировать топологии дорожных сетей, которые изначально устойчивы к формированию "древовидных" или "уникурсальных лабиринтных" паттернов транспортных заторов, тем самым поддерживая состояние "активной турбулентности" (плавного, текущего трафика), а не впадая в заблокированное состояние. Это революционизировало бы управление дорожным движением, перейдя от реактивных решений к проактивному, основанному на физике городскому дизайну, который обеспечивает оптимальный поток и предотвращает общесистемную остановку до того, как она произойдет.

Универсальная библиотека структур

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