Про ланцюжок Фермі–Пасти–Улама–Цингу
Цей симулятор моделює ланцюжок із N точкових мас, з'єднаних однаковими нелінійними пружинами, обидва кінці якого закріплені нерухомо. Рух кожної маси i описується рівнянням mẍᵢ = V′(xᵢ₊₁−xᵢ) − V′(xᵢ−xᵢ₋₁), яке інтегрується симплектичною схемою швидкісного Верле, завдяки чому повна енергія добре зберігається навіть за тривалого прогону. Доступні три моделі пружини: α-FPUT (V = ½x² + ⅓αx³ — невелика кубічна поправка до закону Гука), β-FPUT (V = ½x² + ¼βx⁴ — квартична поправка) і ланцюжок Тода (V = (e^(αx) − αx − 1)/α² — експоненційна пружина, яка інтегрується точно). Якщо поставити повзунок нелінійності на 0, ланцюжок стає суто гармонійним: енергія, вкладена в одну нормальну моду, залишається в ній назавжди — нормальні моди лінійного ланцюжка ніколи не обмінюються енергією.
Під час запуску (і при натисканні «Скинути») ланцюжок збуджується в одній нормальній моді k за формулою xᵢ(0) = A·sin(kπi/(N+1)); для цього слугують повзунки моди збудження, початкової амплітуди й кількості мас N. Верхній холст малює миттєві зміщення кожної маси у вигляді ланцюжка «пружина–кулька», а колір кульки кодує її кінетичну енергію (синій — холодна, червоний — гаряча). Нижній холст — це спектр енергії мод: на кожному кадрі поточний стан проєктується на N дискретних синусних нормальних мод лінійного ланцюжка (дискретне синус-перетворення), і енергія кожної моди Eₖ = ½pₖ² + ½ωₖ²qₖ² накопичується з часом, показуючи, чи «витікає» енергія із початково збудженої моди k до сусідніх мод. Панель статистики показує час симуляції, повну (збережену) енергію, енергію, що лишилась у моді k, і максимальне зміщення в ланцюжку.
Поширені запитання
Що таке проблема FPUT?
У 1953 році Фермі, Паста, Улам і (визнана значно пізніше) Цингу провели один із перших комп'ютерних експериментів: змоделювали ланцюжок із 64 мас із невеликою кубічною (α) нелінійністю в пружинах. Вони очікували, що енергія, вкладена в одну моду, швидко розподілиться рівномірно між усіма модами (термалізація). Натомість вони спостерігали майже-повернення: енергія ненадовго ділилася із сусідніми модами, а потім майже повністю поверталася до вихідної моди — знову і знову, без жодних ознак термалізації. Цей симулятор відтворює саме цей експеримент: збудіть одну моду, дивіться на спектр енергій мод і перевірте, чи повертається енергія назад.
Чим відрізняються моделі α-FPUT, β-FPUT і Тода?
У всіх трьох ідеальна лінійна пружина V=½x² замінюється трохи нелінійною. α-FPUT додає кубічний член (асиметричний — м'якший в один бік і жорсткіший в інший); β-FPUT додає квартичний член (симетричний, жорсткіший в обидва боки); ланцюжок Тода використовує експоненційний потенціал, який є точно розв'язним (інтегрованим) методом оберненого розсіювання й підтримує точні солітонні розв'язки. Збільшення повзунка нелінійності (α або β) підсилює взаємодію між модами; модель Тода майже ідеально повертається до початкового стану за будь-якої амплітуди, тоді як α- і β-FPUT врешті переходять у справжній хаос і термалізацію при достатньо сильній нелінійності.
Чому енергія повертається до початкової моди, а не розсіюється?
Розгадку знайшли Забуський і Краскал у 1965 році: у неперервній границі ланцюжок FPUT поводиться як рівняння Кортевега–де Фріза (KdV), розв'язки якого визначаються солітонами — самоузгодженими хвилями, що проходять одна крізь одну, не втрачаючи форми. Майже-повернення — це дискретний відголосок солітонної динаміки: система достатньо близька до інтегрованої (Тода, KdV), тому обмеження теореми КАМ на термалізацію все ще діють. Лише за сильної нелінійності, коли достатньо сусідніх мод перекриваються й резонують, ланцюжок втрачає цю майже-інтегровану структуру і справді термалізується в хаотичний, рівнорозподілений рух.
Що саме показує панель спектра енергій мод?
На кожному кадрі анімації симулятор проєктує поточні координати й швидкості на N нормальних мод лінійного ланцюжка за допомогою дискретного синус-перетворення (qₖ = Σᵢ xᵢ·sin(kπi/(N+1)), і аналогічно для імпульсів), а потім обчислює гармонійну енергію кожної моди Eₖ = ½pₖ² + ½ωₖ²qₖ² із дисперсійним співвідношенням ωₖ = 2sin(kπ/2(N+1)). Ці енергії мод накопичуються для перших 24 мод і показуються стовпчиками, причому початково збуджена мода k виділена кольором — це саме той діагностичний інструмент, який використали Фермі, Паста, Улам і Цингу у своєму звіті 1955 року.
Чому використовується саме інтегрування Верле, а не простіший метод?
Швидкісний Верле — симплектичний інтегратор: він наближено зберігає гамільтонову структуру рівнянь руху, тому повна енергія коливається навколо сталого значення, а не «дрейфує» на тривалих проміжках часу. Для хаотичної чи майже-інтегрованої системи, як ланцюжок FPUT, несимплектичний метод (наприклад, звичайний Ейлер) вносив би штучну втрату чи набуття енергії й міг би породити хибну термалізацію або хибне повернення. Показник «Повна енергія» у цьому симуляторі має лишатися практично незмінним — це пряма перевірка якості інтегратора.
Що станеться, якщо поставити повзунок нелінійності на нуль?
При α=β=0 усі три моделі зводяться до суто гармонійного ланцюжка, точний розв'язок якого — сума незалежних нормальних мод, що ніколи не обмінюються енергією. Збудження моди k залишить весь спектр енергій зосередженим у стовпчику k, тоді як усі інші стовпчики будуть нульовими — жодних биттів, жодного повернення і жодного хаосу, бо лінійні нормальні моди просто не взаємодіють. Це базовий випадок, з яким порівнювали свій результат автори експерименту FPUT.
Яка роль теореми КАМ у цій моделі?
Теорема Колмогорова–Арнольда–Мозера (КАМ) описує, як за малого нелінійного збурення інтегрованої гамільтонової системи (як-от лінійний ланцюжок чи інтегрований ланцюжок Тода) більшість квазіперіодичних інваріантних торів зберігається, а не одразу розпадається в хаос. Саме тому ланцюжки α- і β-FPUT демонструють довготривале майже-повернення при слабкій і помірній нелінійності замість миттєвої термалізації: система залишається близькою до інтегрованої. Лише коли повзунок нелінійності підняти достатньо високо, перекриття резонансів руйнує тори, що вціліли, і виникає справжній хаотичний обмін енергією між модами.
Як довжина ланцюжка N і мода збудження k впливають на динаміку?
Більше N дає тонше просторове й спектральне розрізнення та дозволяє існувати вищим модам, але й збільшує кількість майже-резонансних трійок/четвірок мод, які можуть зв'язувати нелінійні члени, зазвичай полегшуючи термалізацію за тієї самої нелінійності. Збудження низької моди (мале k, довга хвиля) зазвичай дає найчистіше повернення FPUT, оскільки її найближчі сусіди в спектрі мод добре розділені за частотою; збудження високої моди (k близьке до N) наближає систему до краю зони Бріллюена, де дисперсійне співвідношення сплощується, а сусідні моди легше входять у резонанс.
Чому кульки на верхньому холсті змінюють колір?
Колір кожної кульки кодує її миттєву кінетичну енергію ½mvᵢ², масштабовану відносно найгарячішої кульки в ланцюжку в цю мить (синій → низька кінетична енергія, червоний → висока). Спостереження за кольоровим патерном — швидкий візуальний індикатор того самого питання про розподіл енергії, на яке відповідає спектр мод: у режимі майже-повернення «гаряча точка» досить впорядковано перекочується туди-сюди відповідно до довжини хвилі збудженої моди, а в хаотичному/термалізованому режимі гарячі й холодні кульки стають просторово безладними й приблизно однорідними в середньому за часом.
Чи це та сама система, яку моделювали Фермі, Паста, Улам і Цингу?
Так, структурно: одновимірний ланцюжок мас на однакових нелінійних пружинах із закріпленими кінцями, запущений з одного синусоїдального збудження нормальної моди й проінтегрований у часі на комп'ютері — саме така постановка задачі була у звіті Лос-Аламоської лабораторії 1955 року (спочатку з N=32 або 64 масами й невеликим членом α чи β). Цей симулятор додатково дозволяє підставити точно інтегрований потенціал Тода, який запропонували пізніше (1967) саме для пояснення — через теорію солітонів — чому початковий ланцюжок FPUT відмовлявся термалізуватися.