Научная статья на тему 'Компьютерное моделирование фазовых переходов в двумерных структурах, описываемых трехвершинной антиферромагнитной моделью Поттса на треугольной решетке'

Компьютерное моделирование фазовых переходов в двумерных структурах, описываемых трехвершинной антиферромагнитной моделью Поттса на треугольной решетке Текст научной статьи по специальности «Физика»

CC BY
75
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ФАЗОВЫЙ ПЕРЕХОД / PHASE TRANSITION / КРИТИЧЕСКИЕ ЯВЛЕНИЯ / CRITICAL PHENOMENA / ФРУСТРАЦИИ / FRUSTRATIONS

Аннотация научной статьи по физике, автор научной работы — Бабаев А. Б., Дибирова Г. А., Саидахмедова М. Б.

Методом Монте-Карло проведено компьютерное моделирование фазовых переходов в двумерных структурах, описываемых трехвершинной моделью Поттса на треугольной решетке с учетом антиферромагнитных обменных взаимодействий между ближайшими J 1 и вторыми соседями J 2. Изучены фазовые переходы в этой модели для отношения обменных взаимодействий вторых и ближайших соседей r=J 2/J 1 в диапазоне значений r=0÷2.0. На основе анализа энтропии, теплоемкости, и кумулянтов Биндера четвертого порядка показано, что в модели Поттса с величинами взаимодействий J 1<0 и J 2<0 в интервалах изменений величины 0≤r<0.2 и 1.0

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Бабаев А. Б., Дибирова Г. А., Саидахмедова М. Б.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

COMPUTER SIMULATION OF PHASE TRANSITIONS IN THE TWO-DIMENSIONAL STRUCTURES DESCRIBED THREE-VERTEX ANTIFERROMAGNETIC POTTS MODE ON A TRIANGULAR LATTICE

Using Monte-Carlo simulations, we investigated phase transitions and frustrations in the three-state Potts model on a triangular lattice with allowance for antiferromagnetic exchange interactions between nearestneighbors J 1 and nextnearestneighbors J 2. The ratio of the next-nearestneighbor and nearestneighbor exchange constants r=J 2/J 1 is chosen within the 0÷2 range. Based on the analysis of the entropy, specific heat, system state density function, and fourth order Binder cumulants, the phase transitions in the Potts model with interactions J 1<0 and J 2<0 are shown to be found in value ranges of 0≤r<0.2 and 1.0≤r≤2.0. In an intermediate range 0.2≤r≤1.0 the phase transition fails and the frustrations are revealed.

Текст научной работы на тему «Компьютерное моделирование фазовых переходов в двумерных структурах, описываемых трехвершинной антиферромагнитной моделью Поттса на треугольной решетке»

Для определения магнитуды землетрясения по величине I и а из выражений (1) и (8) имеем:

По известной магнитуде может быть определена интенсивность землетрясения в баллах по шкале МСК по выражению [2]:

где: В и у региональные коэффициенты, для условий Дагестана соответственно равны 1,5 и 3,6.

И - глубина залегания очага.

Вывод. Подключение сейсмографов к ЭВМ, позволит определить в реальном масштабе времени энергетические характеристики землетрясения по математическим зависимостям (8), (9), (10).

Библиографический список:

1. Асланов Г.К., Гаджиев М.М., Исмаилов Т.А., Магомедов Х.Д. О землетрясениях. (Прошлое, современность, прогноз). - Махачкала, ИПЦ ДГТУ, 2001. -98 стр.

2. Быстрицкая Ю.В. Соотношение и сопоставление макросейсмических и инструментальных данных (дагестанские землетрясения). В сборнике «Сейсмичность и гидрогеогазохимия территории Дагестана. Вып. 2 (17).- Махачкала, 1978.

УДК 536.7

Бабаев А.Б., Дибирова Г.А., Саидахмедова М.Б.

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ФАЗОВЫХ ПЕРЕХОДОВ В ДВУМЕРНЫХ СТРУКТУРАХ, ОПИСЫВАЕМЫХ ТРЕХВЕРШИННОЙ АНТИФЕРРОМАГНИТНОЙ МОДЕЛЬЮ ПОТТСА НА ТРЕУГОЛЬНОЙ РЕШЕТКЕ

Babaev A.B., Dibirova G.A., Saidakhmedova M.B.

COMPUTER SIMULATION OF PHASE TRANSITIONS IN THE TWO-DIMENSIONAL STRUCTURES DESCRIBED THREE-VERTEX ANTIFERROMAGNETIC POTTS MODE ON A TRIANGULAR LATTICE

M

l,935(lg a + 1,7341g l )-1,06 1,8

(9)

(10)

Методом Монте-Карло проведено компьютерное моделирование фазовых переходов в двумерных структурах, описываемых трехвершинной моделью Поттса на треугольной решетке с учетом антиферромагнитных обменных взаимодействий между ближайшими J1 и вторыми соседями J2. Изучены фазовые переходы в этой модели для отношения обменных взаимодействий вторых и ближайших соседей r=J2/J1 в диапазоне значений r=0+2.0. На основе анализа энтропии, теплоемкости, и кумулянтов Биндера четвертого порядка показано, что в модели Поттса с величинами взаимодействий J1<0 и J2<0 в интервалах изменений величины 0<r<0.2 и 1.0<r<2.0 наблюдаются фазовые переходы первого рода. В то время как в интервале 0.2<r<1.0 фазовый переход отсутствует, и наблюдаются фрустрации.

Ключевые слова: фазовый переход, критические явления, фрустрации.

Using Monte-Carlo simulations, we investigated phase transitions and frustrations in the three-state Potts model on a triangular lattice with allowance for anti-ferromagnetic exchange interactions between nearest- neighbors J1 and next- nearest- neighbors J2. The ratio of the next-nearest- neighbor and nearest- neighbor exchange constants r=J2/J1 is chosen within the 0^2 range. Based on the analysis of the entropy, specific heat, system state density function, and fourth order Binder cumulants, the phase transitions in the Potts model with interactions J1<0 and J2<0 are shown to be found in value ranges of 0<r<0.2 and 1.0<r<2.0. In an intermediate range 0.2 <r <1.0 the phase transition fails and the frustrations are revealed.

Key words: phase transition, critical phenomena, frustrations.

Исследование фазовых переходов (ФП) и критических явлений (КЯ) в спиновых системах с конкурирующим обменным взаимодействием является одной из центральных в физике конденсированного состояния [1]. Это обусловлено тем, что конкуренция обменного взаимодействия может привести к фрустрации, т.е. такому пространственному расположению магнитных атомов в кристалле, при котором невозможно одновременное антиферромагнитное упорядочение всех взаимодействующих спинов [2].

В последнее двадцатилетие исследованию спиновых систем описываемых ферро- и антиферромагнитными (АФ) моделями Поттса на треугольной решетке с взаимодействиями ближайших соседей было посвящено значительное число работ [3-5] в которых были получены ответы на многие вопросы.

Учет взаимодействия вторых ближайших соседей может привести к возникновению фрустраций, что существенно усложняет решение этой модели. В данной работе нами методом Монте-Карло проведено компьютерное моделирование ФП в АФ трехвершинной модели Поттса и предпринята попытка определить поведение ряда термодинамических параметров этой модели при учете взаимодействия вторых ближайших соседей.

Рисунок 1 - Двумерная трехвершинная АФ модель Поттса с учетом вторых ближайших соседей (при г<0.2) на треугольной решетке

Антиферромагнитная модель Поттса с числом состояний спина q=3 на треугольной решетке с учетом взаимодействия вторых ближайших соседей (см. рис.1) описывается следующим микроскопическим гамильтонианом [6]:

H = - , (1)

2 2

где первый член в формуле (1) учитывает обменное взаимодействие ближайших соседей с величиной Jl<0, а второй - вторых ближайших соседей с величиной J2<0. 61. , 6^ - углы между взаимодействующими спинами & - Б. и Б, - Бк соответственно.

Расчеты проведены для систем с периодическими граничными условиями и с линейными размерами ЬхЬ=Ы, Ь=20^ 144, где Ь - безразмерная величина в единицах межатомного расстояния. При этом соотношение обменного взаимодействия вторых и ближайших соседей r=J2/Jl=0+-2.0.

Для анализа характера ФП и особенностей поведения термодинамических характеристик вблизи критической точки, в такого рода исследованиях, наиболее эффективным методом зарекомендовал себя метод кумулянтов Бин-дера четвертого порядка [7]:

Е4

глт)=1 (2)

3 Е2

к

Ь

m

3 m Vr

(3)

где Е- энергия и т- параметр порядка системы с линейными размерами

Ь.

Методика определения критической температуры нами описана в работах [8-10]. Применение кумулянтов Биндера позволяет также хорошо определить род ФП в системе. ФП первого рода характеризуются следующими отличительными особенностями [11]: усредненная величина Уь,тт(Т) при Т=ТК

тг*

стремится к некоторому нетривиальному значению V согласно выражению

VL,mn(T) - V* + bL

-d

(4)

при Ь ^ да и Т = Тм (Ь), где V* отлично от 2/3, что и продемонстрировано на рис.2 для АФ модели Поттса с д=3 при г=0.1 и ТК=1.261.

Рисунок 2 - Температурная зависимость кумулянтов Биндера VL(T) двумерной трехвершинной АФ модели Поттса при r=0.1

Аналогичное поведение наблюдалось и для величин г=0, 0.167, 1.25, 1.30, 1.50 и 2.00.

Для наблюдения за температурным ходом поведения теплоемкости и восприимчивости использовались флуктуационные соотношения [12]:

С = (КК%и2) - и2), (5)

Ж = (КК )((т2) - (т)2), (6)

где К=| Л /квТ, К=Ь2-число магнитных узлов, [/-внутренняя энергия, т-параметр порядка системы с линейными размерами Ь, угловые скобки означают термодинамическое усреднение.

Рисунок 3 - Температурные зависимости теплоемкости в интервале изменений г, 1.0<г<2.0

Рисунок 4 - Температурные зависимости теплоемкости в интервале изменений г, 0<г<0.167

На рис.3 и 4 представлены характерные температурные зависимости теплоемкости для двумерной АФ модели Поттса при различных отношениях г

обменных параметров, г = —. Как видно из рис.3, уменьшение значения г в

/1

интервале 1.0<г<2.0 приводит к сдвигу максимумов в сторону более низких температур. При этом наблюдается и уменьшение абсолютных значений максимумов теплоемкости.

■"С

♦ ♦ ♦

ад

0.40.3 0.20.1

0.30357

0

0.5

1.5

Рисунок 5 - Зависимость точки фазового перехода в 3-вершинной модели Поттса на треугольной решетке от параметра г/2/1

Рисунок 6 - Зависимость энтропии (при Т=0) в 3-вершинной модели Поттса на треугольной решетке от параметра г=/2//1.

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

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

Кроме того, из рисунка 3 видно, что при r=1.0 очевидной расходимости в критической области не наблюдается, что характерно для фрустрированных систем. И, наоборот, для систем с r=1.25, 1.30, 1.50 и 2.0 в критической области наблюдается явная расходимость. При r=1.25 и 1.3 наблюдается расщепление теплоемкости. Расщепление теплоемкости обычно характерно вблизи точки фрустрации. На рис. 4 такое поведение теплоемкости было обнаружено для значения r=0.167.

Основываясь на результатах данной работы и предыдущих исследований [13, 14], мы построили зависимость точки фазового перехода (рис.5) в 3-вершинной модели Поттса на треугольной решетке от параметра r=J2/J1.

Зависимость нуль-температурной энтропии от величины r приведена на рисунке 6. Как видно из рисунка 5 выделяются две фазы с соответствующими интервалами изменений величины r, 0<r<0.2 и 1.0<r<2.0 для которых имеется дальний порядок по ориентации спинов в каждой из подрешетках. Энтропия этих упорядоченных фаз как видно из рис. 6 равна нулю. В тоже время в интервале изменений величины 0.2<r<1.0 энтропия отлична от нуля, что является характерным признаком отсутствия упорядочения в рассматриваемой магнитной системе.

Вывод. При рассмотренных нами соотношениях между ближайшими и вторыми соседями (Ji<0 и J2<0) в двумерной антиферромагнитной модели Поттса с q=3 на треугольной решетке фрустрации наблюдаются не в одной точке величины r, а на линии - на отрезке 0.2<r<1 (рис.5). В этом интервале изменения r в системе возникают фрустрации. В то время как при других значениях наблюдается фазовый переход, а фрустрация не существует. Работа поддержана грантом РФФИ №13-02-00220.

Библиографический список:

1. Diep H.T. Frustrated spin systems. Singapore: World Scientific Publishing Co. Pte. Ltd., 2004.

2. Гехт Р.С. Успехи физических наук. 1989. Т. 159. С. 261.

3. Murtazaev A.K., Babaev A.B. Phase Transitions in the Two-Dimensional Ferro- and Antiferromagnetic Potts Models on a Triangular Lattice // Journal of Experimental and Theoretical Physics. 2012. V. 142. P. 1189.

4. Муртазаев А.К., Бабаев А.Б., Атаева Г.Я. Фазовые переходы в двумерной ферромагнитной модели Поттса при q=3 на треугольной решетке // Физика низких температур. 2013. Т. 39. С. 194.

5. Муртазаев А.К., Бабаев А.Б. Фазовые переходы в двумерной антиферромагнитной модели Поттса на треугольной решетке // Известия Ран. Серия физическая. 2013. Т. 77. С. 1476.

6. Wu F.Y. The Potts model // Rev. Mod. Phys. 1982. Т. 54. С. 235.

7. Eichhorn K., Binder K.J. Monte Carlo investigation of the three-dimensional random-field three-state Potts model // Phys.: Condens. Matter. 1996. V. 8. P. 5209.

8. Murtazaev A.K., Babaev A.B. Critical properties of the three-dimensional Ising model with quenched disorder // Journal of Magnetism and Magnetic Materials. 2009. V. 321. P. 2630.

9. Murtazaev A.K., Babaev A.B. Phase transitions and critical phenomena in a three-dimensional site-diluted Potts model // Journal of Magnetism and Magnetic Materials. 2012. V. 324. P. 3870.

10. Муртазаев А.К., Бабаев А.Б. Трикритическая точка трехмерной модели Поттса (#=4) c вмороженным немагнитным беспорядком // Письма ЖЭТФ. 2014. Т. 99. С. 618.

11. Loison D., Schotte K.D. First and second order transition in frustrated XY systems // The European Physical Journal B. 1998. V. 5. P. 735.

12. Peczac P., Ferrenberg A.M., Landau D.P. High-accuracy Monte Carlo study of the three-dimensional classical Heisenberg ferromagnet // Phys.Rev. B. 1991. V. 43. P. 6087.

13. Муртазаев А.К., Бабаев А.Б., Магомедов М.А., Касан-Оглу Ф.А., Прош-кин А.И. Фрустрации и фазовые переходы в трехвершинной модели Поттса на треугольной решетке с взаимодействиями вторых ближайших соседей // Письма ЖЭТФ. 2014. Т. 100. С. 267.

14. Murtazaev A.K., Babaev A. B., Kassan-Ogly F.A. The investigation of phase transitions in two-dimensional 3-state antiferromagnetic Potts model on a triangular lattice with interaction of next nearest neighbors // Solid State Phenomena. 2014. V. 215. P. 52.

УДК 621.22 Курбанова З.А.

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ И ЧИСЛЕННЫЙ ЭКСПЕРИМЕНТ ПО ИССЛЕДОВАНИЮ ВЛИЯНИЯ СОСТАВА СООРУЖЕНИЙ ВОДОСЛИВНОГО ФРОНТА ПЛОТИНЫ НА ЗНАЧЕНИЕ УДЕЛЬНЫХ РАСХОДОВ

Kurbanova Z.A.

COMPUTER MODELING AND NUMERICAL EXPERIMENT TO STUDY THE INFLUENCE OF COMPOSITION STRUCTURES SPILLWAY FRONT OF A DAM ON THE VALUE OF THE SPECIFIC EXPENSES

Из-за отсутствия в нормативных документах методик и рекомендаций для определения удельных экономичных расходов водосливных плотин, с уче-

i Надоели баннеры? Вы всегда можете отключить рекламу.