Магнітні Домени

2D модель Ізінга — домени Вейсса, доменні стінки та петля гістерезису B-H

Пресети Матеріалів

Гамільтоніан Ізінга:
H = −J ∑<ij> sisj − μH ∑i si
Метрополіс: приймати переворот, якщо ΔE ≤ 0 або rand < e−ΔE/kT

Управління Петлею B-H

Поточні Значення


0.00

0.00

300

–

50%

–
Про симулятор — Магнітні Домени та Модель Ізінга

Модель Ізінга

2D модель Ізінга — канонічна модель феромагнетизму. Кожен вузол i квадратної ґратки має спін si = ±1. Найближчі сусіди взаємодіють із константою зв'язку J. Зовнішнє магнітне поле H додає зеєманівську енергію. Гамільтоніан:

H = −J ∑ sisj − μH ∑ si

Точна 2D модель Ізінга розв'язана Онзагером (1944). Критична температура Кюрі kBTc = 2J/ln(1+√2) ≈ 2.269J. Нижче Tc з'являється спонтанна намагніченість; вище — система парамагнітна.

Алгоритм Метрополіса

Для симуляції рівноважних конфігурацій при температурі T ми використовуємо алгоритм Метрополіса-Гастінгса: вибираємо випадковий спін, обчислюємо зміну енергії ΔE при його переворті. Приймаємо переворт, якщо ΔE ≤ 0, або з імовірністю exp(−ΔE/kBT), якщо ΔE > 0. Це задовольняє детальному балансу та збігається до розподілу Больцмана.

Домени Вейсса та доменні стінки

Реальний феромагнетик розбивається на магнітні домени, де спіни скоординовані, розділені доменними стінками (стінки Блоха або Нееля). Ширина стінки δ ~ π√(A/K) (A — жорсткість обміну, K — константа анізотропії). Домени заліза — типово 10–300 мкм.

Петля гістерезису B-H

Ключові параметри петлі:

  • Насичення Ms: максимальна намагніченість (всі спіни паралельні)
  • Залишкова намагніченість Mr: M при H=0 після насичення
  • Коерцитивне поле Hc: H, необхідне для зведення M до нуля
  • Втрати на цикл: пропорційні площі петлі (закон Штейнметца)

🌐 English version

Про магнітні домени

Феромагнетик нижче температури Кюрі самовільно розбивається на магнітні домени — ділянки, де всі атомні спіни вирівняні, — розділені вузькими доменними стінками шириною 10–100 нм. Модель Ізінга відтворює цю фізику: кожен спін на вузлі гратки взаємодіє з сусідами через обмінний зв'язок J, що сприяє паралельному вирівнюванню, а теплові флуктуації kBT протидіють упорядкуванню. При критичній температурі відбувається фазовий перехід другого роду, і далекосяжний порядок зникає.

Симуляція використовує алгоритм Метрополіса Монте-Карло на 2D-гратці Ізінга. Можна регулювати температуру відносно точки Кюрі, прикладати зовнішнє поле, щоб спостерігати переміщення доменних стінок, і стежити за анімованою петлею B-H у процесі циклювання поля.

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

Чому феромагнетик розбивається на домени, а не намагнічується рівномірно?

Рівномірна намагніченість створювала б великі поля розсіювання поза матеріалом, що потребує значної магнітостатичної енергії. Матеріал зменшує повну енергію, утворюючи домени з протилежними напрямками намагніченості, що майже взаємно скасовуються зовні. Компромісом є енергія самих доменних стінок (~10⁻³ Дж/м²), тому рівноважний розмір доменів встановлюється з балансу між цими двома внесками.

Що таке модель Ізінга і як вона описує феромагнетизм?

Модель Ізінга ставить у відповідність кожному вузлу гратки спінову змінну si = ±1 і задає енергію E = −J Σ sisj − H Σ si, де J — обмінний інтеграл, H — зовнішнє поле. При J > 0 (феромагнітний зв'язок) і температурах нижче критичної Tc система самовільно намагнічується. Ларс Онзагер аналітично розв'язав двовимірну версію у 1944 р. — одне з видатних досягнень статистичної механіки.

Як алгоритм Метрополіса моделює теплові флуктуації?

На кожному кроці вибирається випадковий спін і обчислюється зміна енергії ΔE при його перекиданні. Якщо ΔE < 0 — перекидання приймається; якщо ΔE > 0 — приймається з імовірністю exp(−ΔE/kBT). Це задовольняє принцип детального балансу та забезпечує правильне больцманівське вибіркування. Після великої кількості кроків гратка приходить до теплової рівноваги з реалістичними доменними структурами.

Що відбувається при температурі Кюрі?

При T = Tc феромагнетик зазнає неперервного фазового переходу другого роду. Спонтанна намагніченість M спадає до нуля як M ∝ (Tc−T)^β з β ≈ 0,326 у 3D. Флуктуації розходяться за масштабом (довжина кореляції ξ → ∞) і часом (критичне сповільнення), породжуючи фрактальні візерунки спінів, видимі у симуляції поблизу Tc.

Що таке доменна стінка і яка її товщина?

Стінка Блоха — це поступовий поворот спінів від орієнтації одного домену до орієнтації сусіднього. Ширина стінки δ = π√(A/K), де A — жорсткість обміну, K — магнітокристалічна анізотропія. В залізі δ ≈ 40 нм; у жорстких магнетиках SmCo₅ з великою K стінки лише ~4 нм, що пояснює ефективне закріплення стінок у таких матеріалах.

Як спостерігають магнітні домени експериментально?

Прямі методи: декорація порошком Біттера (магнітний порошок збирається на доменних стінках), MOKE-мікроскопія (ефект Керра), магнітно-силова мікроскопія (MFM) з нанометровим розрізненням та рентгенівський магнітний круговий дихроїзм (XMCD) для елемент-специфічної візуалізації — необхідні інструменти для розробки носіїв магнітного запису та спінтронних пристроїв.

Чому прикладане поле переміщує стінки, а не перекидає всі спіни одразу?

Переміщення доменної стінки потребує подолання локальних енергетичних бар'єрів у місцях закріплення — набагато менше, ніж когерентний поворот усіх спінів одночасно. Тому при малих полях насамперед рухаються стінки; лише в дуже дрібних частинках (однодоменних зернах менше ~50 нм для Fe) когерентне перемикання є домінуючим, даючи значно вищу коерцитивну силу.

Що таке суперпарамагнетизм?

Коли феромагнітна частинка менша за критичний однодоменний розмір (~10–20 нм для Fe), теплова енергія достатня для випадкового перемикання магнітного моменту всієї частинки між легкими осями. Частинка поводиться як гігантський парамагнітний атом. Це явище є основою нанодисперсних МРТ-контрастних агентів і фундаментальним обмеженням густини запису на жорстких дисках.

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

Магнітострикція пов'язує магнітні та механічні ступені вільності: домени подовжуються або скорочуються вздовж M, а механічні напруження обертають орієнтацію доменів (ефект Віллярі). У нікелі магнітострикція від'ємна, у залізі — додатна. Цей зв'язок уможливлює магнітострикційні актуатори в гідроакустиці і пояснює гул трансформаторних осердь.

Що таке скірміон?

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

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

Це симуляція феромагнетика на основі 2D-моделі Ізінга, яка еволюціонує за алгоритмом Монте-Карло Метрополіса. Кожен вузол ґратки має спін +1 або −1, сусідні спіни взаємодіють через обмінну сталу J, а теплові флуктуації (температура T) протидіють цьому впорядкуванню. Коли прикладене поле H циклічно змінюється, симулятор у реальному часі малює петлю гістерезису B-H матеріалу — з насиченням, залишковою намагніченістю та коерцитивністю.

🔬 Що показано

Ліве полотно показує спінову ґратку кольоровою картою — сині ділянки для спінів угору, червоні для спінів униз, — а праве будує графік намагніченості M залежно від поля H. Спостерігаючи за обома одночасно, видно, як доменні стінки рухаються по ґратці при зміні знаку H і як саме це формує характерну петлю гістерезису.

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

Оберіть пресет матеріалу (залізо, нікель, кобальт, пермалой або м'який ферит) або встановіть Температуру, Прикладене поле H та Обмінну взаємодію J вручну повзунками. Кнопки «Випадкові спіни» та «Упорядкувати ↑» скидають ґратку; перемикайте авто-сканування поля (з регульованою швидкістю) або ручний режим, спостерігаючи, як миттєво оновлюються намагніченість, відсоток спінів угору та кількість доменних стінок.

💡 Чи знали ви?

Точний розв'язок 2D-моделі Ізінга, опублікований Ларсом Онзагером у 1944 році, став одним із перших строгих доведень того, що проста статистична модель здатна демонструвати справжній фазовий перехід — температуру Кюрі, вище якої спонтанна намагніченість зникає.

Поширені запитання

Що саме моделює тут модель Ізінга?

Кожна клітинка ґратки містить спін +1 або −1 — крихітний магнітний момент. Сусідні спіни взаємодіють через обмінну сталу J, що сприяє їхньому паралельному вирівнюванню, а зовнішнє поле H додає власний енергетичний внесок. Симулятор оновлює спіни алгоритмом Метрополіса, дозволяючи ґратці наближатися до теплової рівноваги при обраній температурі.

Чому на сітці з'являються кольорові домени?

Нижче температури Кюрі сусідні спіни об'єднуються в одному напрямку, щоб мінімізувати обмінну енергію, утворюючи видимі сині та червоні домени, розділені доменними стінками. Підвищення температури до точки Кюрі й вище поступово руйнує ці домени, доки вся ґратка не перетвориться на невпорядкований парамагнетик.

Що показує петля гістерезису B-H?

Коли прикладене поле H циклічно змінюється від додатного до від'ємного значення й назад, намагніченість M відстає від нього, а не точно повторює зміни поля, утворюючи замкнену петлю. Ширина петлі відображає коерцитивне поле, потрібне для розмагнічування, а намагніченість при H=0 — це залишкова намагніченість; обидва параметри видно прямо на графіку.

Що змінюють пресети матеріалів?

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

Чому підвищення температури звужує петлю гістерезису?

Вища температура збільшує ймовірність того, що алгоритм Метрополіса прийме енергетично невигідне перекидання спіна, тож теплове збудження протидіє впорядковувальному ефекту обмінної взаємодії J. Коли T наближається до температури Кюрі, спонтанна намагніченість спадає до нуля, і петля гістерезису звужується, аж поки повністю не зникне вище Tc.