НАЦІОНАЛЬНА АКАДЕМІЯ НАУК УКРАЇНИ

ІНСТИТУТ ФІЗИКИ КОНДЕНСОВАНИХ СИСТЕМ

На правах рукопису

ПРИТУЛА Орест Олегович

УДК 536.7, 530.145

ВПЛИВ ЗОВНІШНЬОГО МАГНІТНОГО

ПОЛЯ НА КРИТИЧНУ ПОВЕДІНКУ

ТРИВИМІРНОГО ІЗИНГОПОДІБНОГО МАГНЕТИКА

01.04.02 – теоретична фізика

А В Т О Р Е Ф Е Р А Т

дисертації на здобуття наукового ступеня кандидата

фізико-математичних наук

Львів – 2006

Дисертацією є рукопис.

Робота виконана в Інституті фізики конденсованих систем Національної академії наук України.

Науковий керівник – | | доктор фізико-математичних наук, професор Козловський Михайло Павлович, Інститут фізики конденсованих системи НАН України, завідувач відділу статистичної теорії конденсованого стану.

Офіційні опоненти – | | член-кореспондент НАН України, доктор фізико-математичних наук, професор Стасюк Ігор Васильович, Інститут фізики конденсованих системи НАН України, завідувач відділу квантової статистики. –

| доктор фізико-математичних наук, професор Козицький Юрій Васильович, Інститут математики Університету Марії Кюрі-Склодовської (м.Люблін, Польща), завідувач кафедрою інформатики.

Провідна організація – | | Київський національний університет імені Тараса Шевченка, фізичний факультет, м. Київ

Захист відбудеться ?4 cічня 2007 року о 15 год. 30 хв. на засіданні спеціалізованої вченої ради Д 35.156.01 при Інституті фізики конденсованих систем Національної академії наук України за адресою: 79011 Львів, вул. Свєнціцького, 1.

З дисертацією можна ознайомитись у науковій бібліотеці Інституту фізики конденсованих систем НАН України за адресою: 79026 Львів, вул. Козельницька, 4.

Автореферат розіслано 1 грудня 2006 року.

Вчений секретар

спеціалізованої вченої ради Д .156.01,

кандидат фіз.-мат. наук??? ??Т.Є. Крохмальський

ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

Актуальність теми. Аналіз критичних явищ до цього часу викликає значний інтерес дослідників, що працюють в галузях теоретичної і експериментальної фізики. Актуальність даної тематики спричинена ріноманітністю застосування теорії фазових переходів до опису ряду систем як у фізиці, так і у інших галузях науки [PelissettoVicariPhys. Reports, 2002, 368, 549.]. Критична поведінка тієї чи іншої реальної системи виникає внаслідок колективної взаємодії величезної кількості частинок. Врахувати ці взаємодії за допомогою отримання розв’язків рівнянь руху для кожної окремої частинки є надто складним завданням. Тому в теоретичній фізиці для вирішення таких завдань використовують різноманітні спрощені моделі і наближення, що дозволяють описати основні риси корпоративної поведінки. В загальному випадку модель визначається гамільтоніаном системи, що дає можливість розрахувати статистичну суму

(1)

і, як наслідок, решти фізичних характеристик.

Однією з найбільш простих і наочних моделей, що дозволяє описати поведінку цілого ряду реальних фізичних систем, є добре відома модель Ізинга [Ising E., Zeitschrift fur Physik, 1925, 31, 253]. Одновісні магнетики, прості рідини, бінарні сплави, міцелярні системи це далеко неповний перелік об'єктів, яким притаманна ізингоподібна поведінка. Критична поведінка класу універсальності тривимірної моделі Ізинга виявлена в системах із сильними і електрослабкими взаємодіями, що спостерігаються у фізиці високих енергій [Rummukainen K., Tsypin M., Kajantie K., Laine M., Shaposhnikov M., Nucl. Phys. B, 1998, 532, 283]. Вдалим є також застосування моделі Ізинга для опису явищ в суспільстві, зокрема, в економіці та соціології [W.-X. Zhou, D. Sornette, Preprint www.arxiv.org/physics/0503230]. Ще одним фактором, що спричинив популярність цієї моделі, є наявність точних розв'язків у двовимірному випадку, що дозволило апробувати ряд наближених методів для опису більш складних моделей. Не зважаючи на простоту моделі Ізинга, відсутність точного розв'язку у тривимірному випадку надалі залишає задачу актуальною.

Вирішальний вплив на критичну поведінку системи має зовнішнє поле. Для спінової моделі Ізинга воно є магнітним, але в інших статистичних системах роль поля може відігравати хімічний потенціал, як це має місце, зокрема, в системі рідина-газ або температура - при описі електрослабких взаємодій в стандартній моделі. Присутність зовнішнього поля приводить до розмивання фазового переходу ІІ роду (ФП2). Крім того, в таких системах існує фазовий перехід І роду (ФП1) за полем. Отже у цьому випадку поведінка фізичного об'єкта є більш цікавою. Проте, розрахунок критичної поведінки таких систем на мікроскопічному рівні суттєво ускладнюється і вимагає математичних удосконалень відомих на сьогодні теоретичних методів.

В даній роботі для опису тривимірного ізингоподібного магнетика в зовнішньому полі використовується метод колективних змінних, який був розвинутий для дослідження поведінки ряду рівноважних статистичних систем І.Р. Юхновським [Yukhnovskii I. R., Phase Transitions of the Second Order. Collective Variables Method, World Scientific, Singapore, 1987]. Метод дозволяє отримати на мікроскопічному рівні повний набір фізичних характеристик, зокрема, критичних показників, вільної енергії, рівняння стану системи в околі критичної точки, виходячи з мікроскопічного представлення гамільтоніану.

Зв'язок роботи з науковими програмами, планами і темами. Дисертаційна робота виконана в ІФКС НАН України згідно з планами робіт в рамках держбюджетних тем 0102U000218 "Розвиток кількісної теорії фазових переходів у конденсованих системах" та 0105U002081 "Особливості критичної поведінки конденсованих систем під впливом зовнішнього поля, структурного безладу, фрустрацій та анізотропії".

Метою даної дисертації є мікроскопічний опис критичної поведінки ізингоподібного магнетика, що перебуває в зовнішньому магнітному полі, зокрема, дослідження універсальних і неуніверсальних критичних характеристик системи, що включає розрахунок вільної енергії, рівняння стану та сприйнятливості як функцій поля і температури. Об’єктом дослідження у цій роботі є тривимірна ізингоподібна система, що перебуває в зовнішньому магнітному полі в околі точки фазового переходу. Предметом дослідження є явища, зумовлені впливом зовнішнього поля на таку систему в околі точки фазового переходу другого роду. Методом дослідження є мікроскопічний підхід до опису ізингоподібних систем в околі точки фазового переходу [Юхновський І.Р., КозловськийП., Пилюк І.В. Мікроскопічна теорія фазових переходів у тривимірних системах, Євросвіт, Львів, 2001, 592 c.], який ґрунтується на узагальненому на випадок спінових систем методі колективних змінних (КЗ).

Наукова новизна одержаних результатів.

· На мікроскопічному рівні здійснено опис критичної поведінки тривимірної ізингоподібної системи з врахуванням зовнішнього поля. В рамках моделі отримано ефективні гамільтоніани системи як із врахуванням кубічного доданку, так і при його відсутності, і показано, що вплив поля може бути зведений до лінійного доданку.

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

· Запропоновано метод опису ізингоподібної системи поблизу точки фазового переходу в присутності поля, яке залежить від температури. Значення такого поля характеризується параметром . В залежності від знаку цього параметра поведінка такої системи визначається температурною або польовою змінною. Для різних значень цього параметра отримано явні аналітичні вирази для вільної енергії, параметра порядку і сприйнятливості при температурах вищих і нижчих за температуру фазового переходу. Знайдено значення відповідних критичних показників.

· Вперше отримано аналітичні вирази для вільної енергії, параметра порядку і сприйнятливості ізингоподібної системи в постійному зовнішньому полі як функції поля і температури для високо- і низькотемпературної області. Досліджено кросовер між температурозалежною критичною поведінкою системи і поведінкою, що визначається зовнішнім полем.

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

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

Запропонований підхід застосовується на мікроскопічному рівні і дозволяє аналізувати як колективні ефекти, так і специфіку мікроскопічних взаємодій між частинками. Такі задачі є актуальним для дослідження фазових переходів типу рідина-газ або переходи типу розшарування в рідинах.

Дослідження максимума сприйнятливості системи на температурній шкалі при наявності зовнішнього поля становлять інтерес для експериментаторів, оскільки дають можливість обчислювати одночасно два критичні показники і, як наслідок, всі інші, без надточних вимірювань критичної температури системи.

Особистий внесок здобувача. В спільних публікаціях автору дисертації належить:

· обчислення вільної енергії, параметра порядку і сприйнятливості системи при наявності поля, залежного від температури [1]. Побудова рівняння стану у випадку критичного значення температури [2];

· обчислення вкладів до вільної енергії системи від коротко- і довгохвильових флуктуацій, а також побудова рівняння стану для високотемпературної [3] та низькотемпературної [4] області у випадку, коли поле і температура є незалежними змінними;

· дослідження кросоверу між температурозалежною поведінкою системи і поведінкою, що визначається польовою змінною [3,4,5];

· здійснення поетапного інтегрування статистичної суми з ефективним гамільтоніаном взаємодії, в якому поле представлене лише лінійним доданком; розрахунок ренормгрупових співвідношень і визначення точки виходу з критичного режиму як функції поля і температури. Обчислення вільної енергії системи, а також, побудова скейлінгових функцій для рівняння стану і сприйнятливості [5];

· дослідження максимуму сприйнятливості системи і його поведінку при зміні зовнішнього поля [4,5];

Автор брав безпосередню участь у проведенні числових розрахунків, аналізі та інтерпретації усіх результатів, отриманих в даній роботі.

Апробація роботи здійснена під час доповідей і обговорень основних результатів дисертації на семінарах Інституту фізики конденсованих систем НАН України. Ці результати також доповідалися, дискутувалися на таких конференціях: VI Українсько-польська та II Східноєвропейська нарада з фізики сегнетоелектриків UPEMFP’ 2002 (Ужгород-Синяк, 6-10 вересня 2002р.), II Міжнародна конференція "Фізика рідкої речовини: сучасні проблеми" (Київ, 12-15 вересня 2003р.), конференція імені Боголюбова "Сучасні проблеми математичної і теоретичної фізики" (Київ, 2004р.), Робоча нарада НАТО “Просторові ефекти та нелінійність у фероїках” (Львів, 2004р.), III міжнародна конференція "Фізика рідкої речовини: сучасні проблеми" (Київ, 27-31 травня 2005р.), конференція "Статистична фізика 2005: актуальні проблеми та новітні застосування" (Львів, 28-30 серпня 2005р.) і є опубліковані у їхніх матеріалах.

Результати, що викладено в дисертації, опубліковано в п’яти статтях, виданих у реферованих журналах, зазначених у переліках ВАК України, чотирьох препринтах, а також, в матеріалах та тезах шести міжнародних конференцій.

Структура і об’єм дисертації. Дисертація складається зі вступу, чотирьох розділів, висновків і списку використаних джерел. Робота викладена на 105 сторінках (разом з літературою – 118 сторінок), включає бібліографічний список, що містить 134 найменування у вітчизняних та закордонних виданнях.

ОСНОВНИЙ ЗМІСТ РОБОТИ

У вступі описано обґрунтування актуальності виконаних досліджень, визначено мету роботи і метод дослідження, відзначено новизну отриманих результатів і особистий внесок здобувача у виконання даних досліджень. Подано також зв’язок роботи з науковими програмами і темами, в розробці яких автор приймав безпосередню участь.

У першому розділі "Огляд літератури" на основі літературних джерел проаналізовано основні результати досліджень, виконаних з допомогою інших підходів. Зокрема, проаналізовано роботи, що виконано в рамках теорії збурень та з допомогою комп’ютерного моделювання. Досліджено напрямок незбуреної ренормалізаційної групи (РГ), розглянуто результати, що були отримані з допомогою розв’язку точного рівняння РГ з використанням різного роду наближень для ефективного гамільтоніану системи. Здійснено огляд також попередніх робот, присвячених дослідженню впливу поля на модель Ізинга та особливостям поетапного інтегрування статистичної суми з допомогою методу КЗ.

Другий розділ під назвою "Обчислення статистичної суми в скейлінговій області для різних значень поля і температури" присвячено опису процедури поетапного інтегрування статистичної суми в рамках моделі 4 та вивченню однокомпонентного магнетика у магнітному зовнішньому полі, що залежить від температури. Зокрема, визначено модель, гамільтоніан якої представлений у вигляді

(2)

Тут - змінна -компоненти спіна, що приймає значення . Взаємодія між частинками на вузлах i простої кубічної гратки визначається короткосяжним експонентноспадним потенціалом , де - відстань між частинками. Його наближений фур’є-образ має вигляд

(3)

і містить в собі залежність від сталої гратки і ефективного радіусу взаємодії , що є мікроскопічними параметрами системи. Другий доданок у гамільтоніані (2) описує вплив зовнішнього поля . Визначено множину колективних змінних і представлено статистичну суму системи у початковому вигляді

(4)

що відповідає наближенню моделі . Тут , вимірність простору, s0 =2 є параметром, що визначає параболічну апроксимацію (3) Фур’є-образу потенціалу взаємодії, і символ Кронекера. В процесі розрахунків використовується так зване наближення LPA (local potential approximation), яке не враховує на кожному етапі розрахунку залежності Фур’є-образу потенціалу взаємодії від хвильового вектора. Дане наближення не дозволяє отримати значення критичного показника парної кореляційної функції, проте, дозволяє розрахувати основні риси поведінки даної моделі поблизу точки фазового переходу.

Значення коефіцієнтів у виразі (4) визначаються з допомогою рівності

(5)

Тут , i . Кумулянти , , , де , а є оберненою температурою.

Доданки, що пропорційні непарним степеням КЗ у підінтегральному виразі (4), виникають внаслідок присутності зовнішнього поля. В цьому розділі представлено результати РГ перетворень, що отримано внаслідок поетапного інтегрування статистичної суми (4). Процедура поетапного інтегрування передбачає поетапне відінтегровування змінних, що пов’язані з великими значеннями хвильового вектора. Після здійснення кроків інтегрування коефіцієнти , статистичної суми (4) перенормовуються і визначаються як , на -ій блочній гратці з допомогою рекурентних співвідношень, розв’язки яких в лінійному наближенні відхилень від фіксованої точки мають наступний вигляд

(6)

Тут , , , , . Величини , , , є власними значеннями матриці РГ перетворення, зведена температура (), , , , - координати фіксованої точки, а решта коефіцієнтів виразу (6) є незалежними від поля величинами.

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

(7)

при домінуванні температурної змінної , і

(8)

у випадку, коли визначальною у формуванні критичної поведінки є польова змінна . Ці вирази отримано з другого і першого рівняння рекурентних співвідношень (6) відповідно.

На підставі виразів (7) i (8) визначено області слабких і сильних полів по відношенню до величини відхилень температури Т від її критичного значення Тс. Графічне зображення таких областей показано на рис.1. Умова відповідає випадку слабких полів (області І та IV), де температурна змінна є домінуючою у формуванні критичної поведінки системи. У випадку маємо критичну поведінку, що контролюється польовою змінною, і система перебуває в області сильних полів (області ІI та III). Рівність визначає так звану псевдокритичну лінію,

(9)

що позначена на рисунку числами i для i відповідно. Тут , де є критичним показником кореляційної довжини. Поблизу псевдокритичної лінії використання виразів (7) i (8) є некоректним, оскільки обидві змінні є однаково важливими у формуванні критичної поведінки системи.

В розділі розглянуто випадок, коли поле змінюється при наближенні температури до її критичного значення . Така залежність визначається співвідношенням

(10)

де введено показник , що характеризує траєкторію руху системи до критичної точки. Така задача є актуальною, зокрема, для дослідження псевдоспін-електронної моделі, фазового переходу "рідина-газ" в околі критичної точки. У випадку слабких полів параметр є додатнім (), а для області сильних полів цей параметр є від’ємним (). Представлення (10) дозволяє визначити траєкторію руху системи до точки фазового переходу, що не перетинає псевдокритичної лінії при наближенні до температури і, таким чином, для опису поведінки системи можна користуватися однією точкою виходу ((7) або (8)) зі скейлінгової області. Приймаючи до уваги залежність (10), у розділі представлено розрахунки неуніверсальних характеристик системи, проведені з врахуванням зовнішнього поля, залежного від температури. Зокрема, розраховано явні аналітичні вирази для вільної енергії, параметра порядку і сприйнятливості системи.

Третій розділ має назву "Неуніверсальні характеристики системи як функції незалежних змінних поля і температури". Його присвячено випадку поведінки тривимірної ізингоподібної системи, коли зовнішнє поле і температура є незалежними змінними. Після виконання ітерацій розмір ефективної блочної гратки стає співмірним з кореляційною довжиною системи. Цей факт дозволяє в подальших розрахунках користуватися перенормованим гаусовим розподілом флуктуацій параметра порядку системи. В дисертаційній роботі продемонстровано особливості інтегрування статистичної суми після виходу зі скейлінгової області при наявності поля. Умовою використання гаусового розподілу для здійснення заключного кроку інтегрування є додатній знак і домінуюча роль коефіцієнта біля квадратичного доданку колективної змінної у порівнянні з рештою коефіцієнтів біля вищих степенів КЗ. Показано, що у випадку слабких полів і температур, вищих за критичну, така умова виконується після здійснення додаткового кроку інтегрування. Коли поля є сильними або система перебуває нижче критичної точки при будь-яких значення поля, для виконання вищезгадуваної умови достатньо провести зміщення змінних на величину, пропорційну до параметра порядку.

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

(11)

У області температур, нижчих за критичну, для вільної енергії отримано схожий за структурою вираз

(12)

Коефіцієнти , , i є сталими величинами. Проте, коефіцієнти і у (11) і , у (12) містять додаткову залежність від поля внаслідок присутності лінійного доданку у експоненті статистичної суми. Ця залежність дозволяє розширити межі використання виразу для опису критичної поведінки у порівнянні з розкладами за скейлінговою змінною, але у вузькому околі псевдокритичної лінії застосування цих виразів є неможливим. У випадку сильних полів для вільної енергії системи отримано вираз

(13)

де знаки відповідають випадкам і .

Виходячи з виразів (12) і (11), для вільної енергії знайдено величини параметра порядку і сприйнятливості системи. Проте, як показали результати досліджень величини вкладів до вільної енергії системи від змінних, пов’язаних з різними значенням хвильового вектора, є різними. Основну роль у формуванні фізичних характеристик системи (параметра порядку, сприйнятливості) відіграє вклад, пов’язаний зі змінною . Такий результат є фізично обгрунтованим, оскільки з параметром порядку пов’язане середнє значення саме цієї змінної.

На рис. 2 побудовано графічні залежності для параметра порядку від поля при фіксованому додатньому значенні зведеної температури. Криві 1 і 1’ відповідають повному параметру порядку системи, що співпадає з середнім значенням змінної . Вклади, пов’язані з ненульовими значеннями хвильового вектора компенсують один одного.

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

На рис. 3 побудовано залежність величини як параметра порядку системи від температури для різних значень зовнішнього поля.

У четвертому розділі дисертації під назвою "Узагальнений підхід до опису критичної поведінки ізингоподібної системи" побудовано узагальнений опис критичної поведінки системи, що є дійсним у всій площині (поле-температура), включаючи окіл псевдокритичної лінії. Внаслідок виконання заміни змінних, що призводить до занулення коефіцієнта біля кубічного доданку у виразі для експоненти в статистичній сумі системи, остання має простіший вигляд

Тут

Тут . Отже, в цьому випадку вплив поля є представлений лише лінійним доданком 0 з коефіцієнтом , що пропорційний до поля. Слід зазначити, що при використанні моделей вищих наближень (, ...) процедура "знищення" кубічного доданку також призводить до зникнення решти непарних доданків КЗ окрім лінійного, що дозволяє точніше врахувати вплив зовнішнього поля. Внаслідок здійснення процедури поетапного інтегрування статистичної суми було отримано рекурентні співвідношеня для коефіцієнтів у виразі для експоненти (див. (14)) і їхні розв’язки у лінійному наближенні за відхиленнями від фіксованої точки (). У порівнянні з результатом, отриманим в попередніх розділах (див. (6)), відхилення від фіксованої точки пов’язані з польовою змінною присутні лише у першому рівнянні для величини.

Наступною ключовою відмінністю результатів цього розділу у порівнянні з попередніми є визначення точки виходу зі скейлінової області. Оскільки при скінчених значеннях поля і температури відхилення від фіксованої точки формуються, в основному, в перших двох рівняннях РС (16), то точку виходу слід визначати як функцію обох змінних. Для цього запропоновано рівняння

(17)

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

З огляду на попередні результати, для опису поведінки неуніверсальних характеристик в околі критичної точки враховано вклад від інтегрування за колективною змінною . Для випадку температур вищих від критичної, вільна енергія системи представлена у вигляді так званого мікроскопічного аналога вільної енергії Ландау

(18)

Відповідне рівняння стану можна представити у вигляді

(19)

Варто відзначити, що коефіцієнт містить також залежність від поля, як корінь кубічного рівняння з лінійним доданком, пропорційним до зовнішнього поля.

На рис. 4 представлена залежність параметра порядку як функція поля, побудована на основі виразу (19) (суцільна лінія), і як похідна від повної енергії системи (штрихована лінія). Результат порівняний з даними, отриманими в рамках Монте-Карло моделювання.

Незвичною у високотемпе-ратурній області є поведінка сприйнятливості системи. При скінчених значеннях поля і тем-ператури спостерігається макси-мум сприйнятливості на тем-пературній шкалі. На рис. 5 побудована дана залежність як похідна від величини за польовою змінною при різних значеннях поля. При збільшені величини поля максимум зміщується в ви-сокотемпературну область і роз-мивається. Дослідження макси-муму сприйнятливості є актуаль-ними для експе-риментаторів, оскільки з їх допомогою можна обчислювати значення двох критичних показників без надточних вимірювань критичної температури.

Для впорядкованої фази вільна енергія системи має вигляд

(20)

з рівнянням стану

(21)

Для порівняння отриманих результатів доцільно використати представлення неуніверсальних характеристик у скейлінговій формі, оскільки такі представлення є універсальними. Зокрема, для параметра порядку матимемо

(22)

як функції скейлінгової змінної . Тут , , і . Тут , .

Скейлінгова функція для сприйнятливості має вигляд

(23)

де .

Нар рис. 6, 7 подано графічні представлення скейлінгових функцій (22) i (23), відповідно. Результати порівнянні з даними, отриманими в рамках Монте-Карло моделювання і теоретико-польового підходу з використанням теорії збурень. Як видно з рисунків, на якісному рівні отримані в дисертації результати досить добре узгоджуються з найбільш точними на даний час даними. Розбіжності на кількісному рівні пояснюються насамперед наближенням LPA, що впливає на значення критичних показників через відсутність показника . Крім того, як було показано в попередніх роботах, виконаних в рамках методу колективних змінних при відсутності поля, з допомогою моделі можна покращити результат в кількісному відношенні. Проте, на відміну від інших результатів, результати цієї роботи отримано, виходячи з гамільтоніану системи, що містить мікроскопічні характеристики системи, не містять феноменологічних припущень і підгоночних параметрів.

Основні результати та висновки.

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

2. В рамках мікроскопічного підходу знайдено вирази для вільної енергії, параметра порядку та сприйнятливості моделі Ізинга в зовнішньому полі як функції температури для фіксованого значення поля. Знайдені значення відповідних критичних показників добре узгоджуються з результатами, отриманими в рамках незбуреної теорії критичних явищ при розв’язанні точного ренормгрупового рівняння в наближенні локального потенціалу (LPA).

3. Досліджено величини вкладів до вільної енергії системи поблизу точки фазового переходу від колективних змінних, пов’язаних з різними значеннями хвильового вектора. Встановлено, що вклад пов’язаний з нульовим значенням хвильового вектора є визначальним у формуванні величини параметра порядку.

4. Досліджено випадок зовнішнього магнітного поля, залежного від температури. Введено показник , що дозволяє розглядати різні типи залежності і використовувати відповідні визначення точки виходу. Показано, що для додатніх значень цього показника поведінка системи визначається температурною змінною і при від’ємних критичну поведінку визначає польова змінна. Отримано явні аналітичні вирази для вільної енергії параметра порядку і сприйнятливості системи як функції поля, що залежить від температури.

5. Побудовано скейлінгові фукції рівняння стану і сприйнятливості системи. Досліджено поведінку сприйнятливості системи як функції температури при різних значеннях зовнішнього поля. Показана наявність максимума цієї функції і її скейліногової форми у високотемпературній області. Його положення на шкалі скейлінгової змінної добре узгоджується з результатами, отриманими в рамках теоретико-польового підходу і МК моделювання. Показано, що зі збільшенням значень зовнішнього поля, максимум сприйнятливості розмивається і зсувається у високотемпературну область.

Результати дисертації опубліковано в таких роботах:

1. Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Behaviour of a three-dimensional uni-axial magnet near the critical point in an external field.// Condens. Matter Phys. — 2004. — Vol. . — p. –382.

2. Kozlovskii M. P., Pylyuk I. V., Prytula O. O. Microscopic description of the critical behavior of three-dimensional Ising-like systems in an external field.// Phys. Rev. B — 2006. — Vol . p. 1-13.

3. Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Critical behaviour of a three-dimensional one-component magnet in strong and weak external fields at T>Tc.// Physica A — 2006. — Vol. . — p. .

4. Kozlovskii M. P., Pylyuk I. V. , Prytula O. O. Behaviour of the order parameter of the simple magnet in an external field.// Condens. Matter Phys. — 2005. — Vol. . — p. .

5. Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Free energy and equation of state of Ising-like magnet near the critical point.// Nucl. Phys. B 2006. — Vol . — p. 242-251.

6. Козловський М.П., Пилюк І.В., Притула О.О. Поведінка тривимірного магнетика поблизу критичної точки за наявності зовнішнього поля: Препринт ICMP-03-21U. — Львів: ІФКС НАН України, 2003. – 42 с.

7. Козловський М.П., Пилюк І.В., Притула О.О. Термодинамічні функції однокомпонентної моделі магнетика при Т>Тс за наявності зовнішнього поля. Препринт ICMP-04-03U. — Львів: ІФКС НАН України, 2004. – 33 с.

8. Козловський М.П., Пилюк І.В., Притула О.О. Термодинамічні функції однокомпонентної моделі магнетика при Т<Тс за наявності зовнішнього поля. Препринт ICMP-04-06U. — Львів: ІФКС НАН України, 2004. – 33 с.

9. Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Free energy and equation of state of Ising-like magnet near the critical point. Preprynt www.arxiv.org/hep-th/0609105, 2006 – 13p.

10.

10 Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Microscopic analog of the Landau Free Energy of Three-Dimensional Ising-Like Systems.// VI Ukrainian - Polish and II East-European Meeting on Ferroelectrics Physics, Uzhgorod-Synjak, Ukraine, September 6-10, 2002. — p. .

11. Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Influence of an external field on the critical behaviour of a three-dimensional magnet.// 2nd International Conference "Physics of Liquid Matter: Modern Problems", Kyiv, Ukraine, September 12-15, 2003. — - p. 58.

12. Kozlovskii M.P., Prytula O.O. Behaviour of 3D statistical systems near the phase transition point in an external field.// Bogolyubov Kyiv conference "Modern Problems of Mathematics and Theoretical Physics", Kyiv, Ukraine, 2004. — p. 

13. Pylyuk I.V., Kozlovskii M.P. and Prytula O.O. Analytic Method of Calculating Thermodynamic Functions for Ising-Like System in an External Field.// Conference "Dimesionality Effects Non-linearity in Ferroics" Ukraine, Lviv, October 19-22, 2004. — p. .

14. Kozlovskii M.P., Pylyuk I.V., Prytula O.O. Influence of an external field on the critical behaviour of a three-dimensional magnet.// 3nd International Conference "Physics of Liquid Matter: Modern Problems", Kyiv, Ukraine, May 27-31, 2005. — p. 101.

15. Kozlovskii M.P, Pylyuk I.V., Prytula O.O Behaviour of the order parameter of the simple magnet in an external field.// Conference "Statistical physics 2005", Ukraine, Lviv, August 28-30, 2005. — p. .

Притула О.О. Вплив зовнішнього магнітного поля на критичну поведінку тривимірного ізингоподібного магнетика. – Рукопис.

Дисертація на здобуття наукового ступеня кандидата фізико-математичних наук за спеціальністю 01.04.02 – теоретична фізика, Інститут фізики конденсованих систем Національної академії наук України, Львів, 2006.

В дисертаційній роботі представлено мікроскопічний опис критичної поведінки тривимірного ізингоподібного магнетика в присутності зовнішнього магнітного поля. Підхід базується на використанні методу колективних змінних, що включає ренормгрупові перетворення коефіцієнтів ефективного гамільтоніану системи як результат поетапного інтегрування статистичної суми. Обчислення виконано з використанням четвірного розподілу флуктуацій параметра порядку. Використовуючи різні способи визначення необхідної кількості кроків інтегрування (точки виходу зі скейлінгової області), отримано відповідні вирази для вільної енергії, параметра порядку і сприйнятливості системи як функції поля і температури. Розглянуто також критичну поведінку системи у випадку зовнішнього поля, що залежить від температури. Досліджено поведінку системи в околі псевдокритичної лінії і встановлено, що для отримання якісних результатів необхідно враховувати одночасне обмеження кореляційної довжини польовою і температурною змінними. Тому точка виходу зі скейлінгової області повинна бути функцією поля і температури. Побудовано скейлінгові функції для параметра порядку і сприйнятливості як універсальні характеристики системи. Результати якісно узгоджуються з даними, отриманими в рамках Монте-Карло моделювання і теоретико-польового підходу в теорії збурень, проте їх отримано без використання феноменологічних припущень і підгоночних параметрів.

Ключові слова: критична точка, параметр порядку, модель Ізинга.

Притула О.О. Влияние внешнего магнитного поля на критическое поведение трёхмерного изингоподобного магнетика. – Рукопись.

Диссертация на соискание ученой степени кандидата физико-математических наук по специальности 01.04.02 – теоретическая физика, Институт физики конденсированных систем Национальной академии наук Украины, Львов, 2006.

В дисертационной работе представлено микроскопическое описание критического поведения трёхмерного изингоподобного магнетика в присутствии внешнего магнитного поля. Подход базируется на использовании метода коллективных переменных, который включает ренормгрупповые преобразования коэффициентов эффективного гамильтониана системы как результат поэтапного интегрирования статистической суммы. Расчеты выполнены с использованием четверного распределения флуктуаций параметра порядка. Используя разные способы определения необходимого числа шагов интегрирования (точки выхода из области скейлинга), получено соответсвующие выражения для свободной энергии, параметра порядка и восприимчивости системы как функции поля и температуры. Рассмотрено также критическое поведение системы в случае поля, зависящего от температуры. Исследовано поведение системы в непосредственной близости к псевдокритической линии и установлено, что для получения качественных результатов необходимо принимать во внимание одновременное ограничение корреляционной длины полевой и температурной переменными. Следственно, точка выхода должна быть функцией поля и температуры. Построены также скейлинговые функции для параметра порядка и восприимчивости как универсальные характеристики системы. Результаты качественно согласуются с данными, полученными в рамках Монте-Карло моделирования и теоретико-полевого подхода в теории возмущений, но получены без использования феноменологических предположений и подгоночных параметров.

Ключевые слова: критическая точка, параметр порядка, модель Изинга.

PrytulaInfluence of an external magnetic field on the critical behaviour of the tree-dimensional Ising-like magnet. – Manuscript.

Thesis on search of the scientific degree of candidate of physical and mathematical sciences, speciality 01.04.02 – theoretical physics, Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, Lviv, 2006.

We study the behavior of a three-dimensional Ising-like system in an external field near the critical point by using the non-Gaussian spin-density fluctuations, namely, the quartic measure density with the odd part of the effective Hamiltonian. The basic idea of the new analytic method for deriving complete expressions of the thermodynamic characteristics is described. Using the renormalization group transformations, the scaling region size is defined as a function of temperature and field. The critical exponents, which are comparable with ones obtained within the framework of LPA, are calculated. The expressions for the total free energy, order parameter, susceptibility, entropy, and specific heat of the system are obtained as functions of the temperature and field. The contributions to the order parameter, which are related to the integration with respect to the variables with different values of the wave vector, are investigated. It is established that the main contribution is formed by the average value of the variable .

The calculations are also performed in the case, when the field and temperature are dependent and related by some expression (the system tends to the critical point along some trajectory), which contains the parameter . The critical behaviour of the system is considered for different values of this parameter .

The crossover between temperature-dependent and field-dependent critical behaviour is investigated. It is established, for obtaining valid result one should consider the point of exit, which determines the scaling region size, as function of two variables, namely, field and temperature. The equation for finding such a function is proposed and solved numerically. From this point of view, the free energy, order parameter and susceptibility are obtained as functions of the temperature and field.

The behaviour of the susceptibility as function of the temperature is investigated for different values of the external field. The characteristic exhibit the maximum, which is shifted in the high-temperature region and smeared with increasing the external field.

The results agree qualitatively with ones obtained within the framework of the parametric representation of the equation of state and Monte-Carlo simulations. The calculations do not involve any parametrization, phenomenological assumptions and adjustable parameters. The approach can be extended to models with a multicomponent order parameter.

Keywords: critical point, order parameter, Ising model.