Головна Нейронауки та Біофізика Нейральні Осцилятори

🧠 Нейральні Осцилятори

Нейральні осцилятори з моделлю Курамото, ритми мозку, гамма-осциляції та баланс збудження-гальмування.

Нейронауки та Біофізика2DСередній60 FPS
neural-oscillators ↗ Відкрити окремо
DRAG · SCROLL · CLICK — керуйте прямо у вікні симуляції.

Про Нейральні Осцилятори

Модель Курамото

Модель Курамото описує N зв'язаних осциляторів, кожен із природною частотою ωᵢ, взятою з розподілу g(ω). Зв'язок K зближує фази: dθᵢ/dt = ωᵢ + K/N·Σsin(θⱼ−θᵢ). Параметр порядку r = |1/N·Σe^(iθⱼ)| вимірює синхронність: r≈0 — некогерентність, r≈1 — повна синхронізація. Фазовий перехід відбувається при критичному зв'язку K_c = 2/(πg(0)) — аналогічно магнітному фазовому переходу.

Ритми мозку

Мозок генерує характерні осциляторні патерни, які вимірюються за допомогою ЕЕГ. Дельта (1-4 Гц) домінує під час глибокого сну; тета (4-8 Гц) з'являється під час кодування пам'яті та REM-сну; альфа (8-13 Гц) відображає розслаблене неспання; бета (13-30 Гц) супроводжує активне мислення; гамма (30-80 Гц) корелює зі сприйняттям і свідомістю. Ці ритми виникають із синхронізованих популяцій нейронів — подібно до моделі Курамото.

Гамма та баланс збудження-гальмування

Гамма-осциляції (30-80 Гц) виникають унаслідок балансу між збуджувальними (AMPA, NMDA) і гальмівними (GABA) синаптичними струмами. Швидкоспайкові інтернейрони забезпечують точно синхронізоване гальмування, яке задає ритм збуджувальним клітинам. Цей баланс збудження-гальмування (E-I) критично важливий: надлишок збудження спричиняє епілептичні напади; надлишок гальмування пригнічує будь-яку активність. Порушення гамма-ритму пов'язують із шизофренією та хворобою Альцгеймера.

Про цю симуляцію

Ця симуляція показує, як популяція осциляторів Курамото — кожен із власною природною частотою — спонтанно об'єднується у спільний ритм зі зростанням сили зв'язку між ними. Це той самий механізм синхронізації, який, як вважають, лежить в основі ритмів мозку, таких як дельта, тета, альфа, бета й гамма, які тут вимірюються наживо за допомогою параметра порядку та спектра потужності середнього поля популяції.

🔬 Що показано

N фазових осциляторів розташовані на кільці, кожен еволюціонує за формулою dθᵢ/dt = ωᵢ + K/N·Σⱼsin(θⱼ−θᵢ) плюс опціональний шум, де ωᵢ взято з розкидом σ навколо середньої частоти f0. Параметр порядку r = |1/N·Σe^(iθⱼ)| відстежує синхронність від 0 (некогерентність) до 1 (повна синхронізація), а живе ДПФ середнього поля показує, у якому діапазоні ритмів мозку — δ, θ, α, β чи γ — коливається популяція.

🎮 Як користуватись

Регулюйте кількість осциляторів N (5–80), Зв'язок K (0–10), Розкид частот σ (0.1–5) і Рівень шуму (0–2), або перейдіть одразу до п'яти пресетів: Десинхронізовано, Тета (6 Гц), Альфа (10 Гц), Гамма (40 Гц) і Глибокий сон (δ). Спостерігайте, як точки на фазовій діаграмі збираються разом, а результуючий вектор (стрілка параметра порядку) росте зі зростанням K вище критичного значення; панелі Порядок r(t) і Спектр потужності оновлюються наживо внизу.

💡 Чи знали ви?

Модель Курамото має точний фазовий перехід при критичному зв'язку K_c = 2/(π·g(0)), де g(0) — щільність розподілу природних частот у середній точці — нижче K_c осцилятори дрейфують незалежно, вище нього частина з них спонтанно синхронізується, подібно до зв'язаних маятникових годинників або блискання світлячків. Вважається, що гамма-синхронність (30–80 Гц) у реальній корі залежить від подібного балансу, але між збуджувальними та гальмівними (E-I) нейронами, а не від єдиної константи зв'язку.

Часті питання

Що таке модель Курамото?

Це математична модель N зв'язаних фазових осциляторів, кожен із власною природною частотою ωᵢ, взятою з розподілу. Фаза кожного осцилятора θᵢ еволюціонує за формулою dθᵢ/dt = ωᵢ + K/N·Σⱼsin(θⱼ−θᵢ), тож кожен осцилятор тягнеться до середньої фази всіх інших із силою, пропорційною константі зв'язку K. Це одна з найпростіших розв'язних моделей спонтанної синхронізації, яку використовують у фізиці, біології та нейронауці.

Що насправді вимірює параметр порядку r?

r = |1/N·Σe^(iθⱼ)| трактує фазу кожного осцилятора як одиничний вектор на комплексній площині й усереднює їх. Якщо фази розкидані випадково, вектори скасовують одне одного, і r лишається близьким до 0; якщо осцилятори синхронізуються в одній фазі, вектори додаються, і r наближається до 1. Симуляція малює цей результуючий вектор безпосередньо на кільці фазової діаграми, тож можна спостерігати, як синхронність наростає в реальному часі зі зростанням K.

Чому осцилятори раптово синхронізуються, коли K перевищує поріг?

Модель Курамото має аналітично виведений критичний зв'язок K_c = 2/(π·g(0)), де g(0) — значення розподілу природних частот у його центрі. Нижче K_c розкид природних частот σ домінує, і осцилятори розходяться, тримаючи r близьким до 0. Вище K_c дедалі більша частка осциляторів синхронізується на спільній частоті, і r різко зростає — фазовий перехід, математично подібний до появи намагніченості у феромагнетику під час зниження температури.

Як пресети співвідносяться з реальними ритмами мозку?

Кожен пресет задає середню частоту f0 і комбінацію зв'язку/розкиду, налаштовану на класичний діапазон ЕЕГ: Тета цілиться в 6 Гц (кодування пам'яті, REM-сон), Альфа — у 10 Гц (розслаблене неспання), Гамма — у 40 Гц (сприйняття, високий зв'язок K=8 для щільної синхронії), а Глибокий сон — у 2 Гц із сильним зв'язком і низьким розкидом, імітуючи високосинхронізовану повільнохвильову активність глибокого не-REM сну. Пресет Десинхронізовано використовує низький K і високий σ, тож популяція ніколи не синхронізується, ілюструючи некогерентний базовий стан.

Як визначається домінантна частота й діапазон із симуляції?

Симуляція веде плинну історію середнього поля популяції (дійсна частина вектора параметра порядку) і виконує дискретне перетворення Фур'є за останніми 512 вибірками. Частотний бін із найбільшою величиною стає показником «домінантної частоти», який потім класифікується за діапазоном — δ нижче 4 Гц, θ нижче 8 Гц, α нижче 13 Гц, β нижче 30 Гц, γ вище цього — відображаючи те, як нейронауковці класифікують спектри потужності ЕЕГ за канонічними діапазонами ритмів.

Схожі симуляції