Теоретико-полевые и численные исследования критического поведения сложных однородных и структурно неупорядоченных систем, описываемых многовершинными моделями тема автореферата и диссертации по физике, 01.04.07 ВАК РФ

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

485023»

ПРУДНИКОВ ПАВЕЛ ВЛАДИМИРОВИЧ

ТЕОРЕТИКО-ПОЛЕВЫЕ И ЧИСЛЕННЫЕ ИССЛЕДОВАНИЯ КРИТИЧЕСКОГО ПОВЕДЕНИЯ СЛОЖНЫХ ОДНОРОДНЫХ И СТРУКТУРНО НЕУПОРЯДОЧЕННЫХ СИСТЕМ, ОПИСЫВАЕМЫХ МНОГОВЕРШИННЫМИ МОДЕЛЯМИ

01.04.07 - физика конденсированного состояния

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

1 6 июн 2011

Махачкала - 2011

4850239

Работа выполнена на кафедре теоретической физики Омского государственного университета им. Ф.М. Достоевского

Официальные оппоненты: доктор физико-математических наук,

член-корреспондент РАН, профессор Муртазаев Акай Курбанович

Ведущая организация:

Казанский физико-технический институт им. Е.К. Завойского КазНЦ РАН

Защита состоится " 28" июля 2011 г. в 16-00 на заседании диссертационного совета Д 002.095.01 при Институте физики ДагНЦ РАН но адресу: 367003, г. Махачкала, пр. Шамиля, 39 А

Отзывы на автореферат просьба направлять по адресу:

367003, г. Махачкала, ул. М. Ярагского, д. 94,

Институт физики ДагНЦ РАН, секретарю диссертационного совета

С диссертацией можно ознакомиться в библиотеке Института физики ДагНЦ РАН

Автореферат разослан мая 2011 г.

Ученый секретарь диссертационного совета,

доктор физико-математических наук, профессор Суслов Игорь Михайлович

доктор физико-математических наук, профессор Аплеснин Сергей Степанович

доктор физико-математических наук

А.Б. Батдалов

Общая характеристика работы

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

Выявленная общность свойств фазовых переходов второго рода в различных веществах позволила сформулировать принцип универсальности критических явлений и предложить модель, в основе которой лежала гипотеза масштабного подобия флуктуаций (А.З.Паташинский, В.Л.Покровский, 1964,1966; Л.Кадаиов, 1966). Идеи использования метода ренормализационной группы и последующая их иллюстрация с помощью метода £ - разложения (К.Вильсон, М.Фишер, 1971, 1972) позволили сделать еще несколько шагов в качественном понимании фазовых переходов и в их количественном описании. Дальнейшее развитие этих идей привело к появлению теоретико-полевого подхода к описанию критических явлений (С.Ди Кастро, Г.Иона-Лазинио, Г.Паризи, С.Л.Гинзбург, 1975), дающему более точные количественные результаты для непосредственно трехмерных систем при применении методов суммирования асимптотических рядов (Г.Бейкер, В.Никел и др., 1976; Е.Брезин и др., 1977; И.М.Суслов, 2001, 2008).

В критической точке наряду с особенностями равновесных характеристик сингулярное поведение демонстрируют кинетические коэффициенты и динамические функции отклика, что обусловлено аномально большими временами релаксации сильно флуктуирующих величин. Критическая динамика исследовалась ренормгрупповыми методами, совмещенными с £-разложением, в работах ряда авторов (Б.Гальперин, П.Хоэнберг, Ш.Ма, 1972, 1974; Ж.Мазенко, К.Кавасаки, 1975 и др.). Однако исследование динамических свойств критических флуктуаций сталкивается с трудностями более сложными, чем при описании равновесных свойств. Это вызвано как необходимостью учета взаимодействия флуктуаций параметра порядка с другими долгоживу-щими возбуждениями (Б.Гальперин, П.Хоэнберг, 1967; А.М.Поляков, К.Кавасаки, 1969), так и более слабой сходимостью получаемых рядов по £ и большим числом существенных диаграмм уже в низших порядках теории возмущений. В динамике фазовых переходов и понятие универсальности приобретает более широкий смысл - становятся существенными законы сохранения для локальной плотности долгоживущих переменных. Для описания критической динамики в однородных и неупорядоченных системах был развит теоретико-полевой вариант ренормгруп-пы, позволяющий исследовать динамику трехмерных и двумерных систем без использования г-разложения (В.В.Прудников, 1992, 1997, 1998; К.Ёрдинг, 1995).

Одной из нерешенных задач теории критических явлений остается описание неравновесного критического поведения макроскопических систем, далеких от состояния равновесия. Это, прежде всего, относится к явлениям критической релаксации однородных и структурно неупорядоченных систем при фазовых переходах второго рода и фазовых переходах первого рода близких ко второму. Критическое замедление времени релаксации и аномально большие времена корреляции различных состояний для данных систем приводят к реализации динамического скейлингового поведения даже когда системы находятся в состояниях, далеких от состояния равновесия. Значительный интерес к исследованию подобных систем обусловлен предсказываемыми и наблюдаемыми в них свойствами старения при медленной эволюции систем из неравновесного начального состояния и нарушениями флуктуационно-диссипативной теоремы (Г.Янсен и др.,1989; Б.Ченг и др.,1999; П.Калабрезе,2005).

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

Наибольших успехов исследователи достигли при изучении влияния некоррелированных дефектов с эффектами случайной локальной температуры на критическое поведение неупорядоченных систем. Ренормгрупповой анализ с использованием е - разложения (Д.Е.Хмельницкий, Т.Лубенски, 1975), а затем в рамках более точного теоретико-полевого подхода (А.И.Соколов, Б.Н.Шалаев, 1981,1989; И.О.Майер, 1989; Д.В.Пахнин, А.И.Соколов, 1999) выявил, что поведение неупорядоченных изингоподобных систем характеризуется новым набором критических индексов. Эксперимент (Р.Биржено и др., 1983,1988) подтвердил численное отличие статических критических индексов для неупорядоченных систем от их значений для однородных систем и показал хорошее согласие с теоретическими результатами. Однако влияние дефектов структуры на критическую динамику значительно менее исследовано (У.Крей, Г.Гринстейн и др., 1977, Прудников В.В., 1983,1992), несмотря на ожидаемое более сильное их проявление в неравновесных характеристиках. По критической динамике разбавленных систем до сих пор существует скудная экспериментальная информация (Д.Беланджер и др.,1988). При этом достигнутая точность результатов низка для достоверной проверки результатов теоретических расчетов. Нет и

достаточно обоснованных теоретических оценок динамического индекса г.

Для описания критического поведения неупорядоченных систем используются многовершинные модели, для которых предсказания, сделанные на основе применения метода е- разложения, не являются надежными. Это объясняется конкуренцией различных типов критического поведения в многопараметрическом пространстве модели, что делает протяжку € —> 1,2 невозможной без пересечения областей стабильности различных фиксированных точек. Для получения достоверных результатов требуется разработка более надежных методов описания.

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

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

Цель работы

1. Развитие методики и осуществление теоретико-полевого описания неравновесного критического поведения однородных и структурно неупорядоченных систем с замороженными дефектами структуры в многопетлевом приближении с применением методов суммирования. Исследование влияния неравновесных начальных состояний и создаваемых ими нарушений трансляционной симметрии во времени на медленную эволюцию систем с сильно коррелированными состояниями в критической точке. Теоретико-полевое и численное исследование статического и динамического критического поведения неупорядоченных трехмерных систем с дальнодейству-ющей пространственной корреляцией дефектов структуры.

2. Численное исследование критического поведения неупорядоченной ферромагнитной модели Изинга численными методами Монте-Карло. Компьютерное моделирование неравновесных процессов критической релаксации трехмерных систем методом коротковременной динамики при изменении концентрации дефектов структуры в широком интервале от уровня слабого разбавления до концентраций, близких к порогу перколяции. Проведение сопоставления результатов теоретико-полевого расчета с результатами компьютерного моделирования критической динамики однородных и слабо неупорядоченных систем. Выяснение вопроса об универсальности неравновесного критического поведения неупорядоченных систем.

3. Ренормгрупповое исследование критического и мультикритического поведения сложных однородных и структурно неупорядоченных спиновых систем. Разработка методики теоретического описания влияния дефектов структуры и эффектов их корреляции на аномальное поведение характеристик распространения ультразвука в твердых телах при температуре фазо-

вого перехода второго рода. Исследование влияния поверхности и ее ориентации относительно направления анизотропии кристалла на мультикритическое поведение системы в окрестности обобщенной точки Лифшица т-го порядка.

Научим новизна результатов

1. Впервые проведен расчет значений динамического критического индекса г, определяющего релаксационные свойства системы, в высоких порядках теории для однородных и структурно неупорядоченных трехмерных и двумерных изингоподобных систем с последовательным применением различных методов суммирования асимптотических рядов.

2. Осуществлено теоретико-полевое описание влияния неравновесных начальных значений параметра порядка на его эволюцию в критической точке для однородных систем. В рамках диссипативной модели впервые в рекордном трехпетлевом приближении проведен расчет критического индекса коротковременной эволюции в1 с использованием метода е-разложения.

3. Разработано и осуществлено исследование одновременного влияния неравновесных начальных состояний, дефектов структуры и эффектов их корреляции на эволюцию структурно неупорядоченных систем в критической точке. Впервые реализовано теоретико-полевое описание неравновесного критического поведения непосредственно трехмерных систем и проведен расчет критического индекса коротковременной эволюции в' в двухпетлевом приближении без использования г-разложения. Значения динамических критических индексов 0' и г, полученные с применением методов суммирования, находятся в лучшем соответствии с результатами компьютерного моделирования, чем результаты применения метода е-разложения.

4. Осуществлены численные Монте-Карло исследования равновесного критического поведения неупорядоченной ферромагнитной модели Изинга в широком интервале изменения концентрации немагнитных атомов примеси. Для термодинамических и корреляционных характеристик методом конечноразмерного скейлинга впервые определены универсальные скейлин-говые функции, вычислены критические температуры и статические критические индексы с учетом ведущих поправок к скейлингу. На основе выявленных концентрационных изменений скейлинговых функций и значений критических индексов сделан вывод о существовании двух универсальных классов критического поведения разбавленной модели Изинга с различными характеристиками для слабо и сильно неупорядоченных систем.

5. Осуществлено численное исследование методом коротковременной динамики неравновесного критического поведения для анизотропных ферромагнетиков, описываемых моделью Изинга, в широкой области изменения концентрации точечных дефектов структуры для различных начальных неравновесных состояний системы. На основе выявленной концентрационной зависимости динамических и статических характеристик впервые показано, что неравновесное поведение слабо и сильно неупорядоченных систем также принадлежит к различным классам универсальности критического поведения с несовпадающими в пределах статистических погрешностей значениями динамических показателей в' и 2.

С. Впервые в рамках теоретико-полевого подхода в двухпетлевом приближении с применением методов суммирования осуществлено описание критического поведения трехмерных структурно неупорядоченных систем с нарушенной репличной симметрий (НРС) без использования £-разложения. Для систем с произвольной размерностью от 3 до 4 проведен ренормгрупповой

анализ эффективного гамильтониана модели с НРС. Выявлена устойчивость критического поведения трехмерных систем относительно эффектов НРС.

7. Впервые осуществлено теоретико-полевое описание статического и неравновесного критического поведения непосредственно трехмерных систем с дальнодействующей корреляцией дефектов для различных значений числа компонент параметра порядка и показателя корреляции а с применением методов суммирования. Полученная картина областей устойчивого критического поведения данных систем и значения критических индексов существенно отличаются от предсказанных ранее в рамках двухпараметрического е,в - разложения. Показано, что эффекты дальнодействующей корреляции дефектов меняют критическое поведение всех систем, для описания критического поведения которых используются базовые теоретические модели Изинга, ХУ и Гейзенберга с многокомпонентным параметром порядка, в то время как некоррелированные дефекты модифицируют критическое поведение только изингоподобных систем.

8. Проведено численное исследование методом коротковременной динамики критического поведения различных модельных трехмерных систем с линейными протяженными дес{>ектами в широкой области изменения концентрации дефектов. Выявлены особенности влияния дальнодействующей корреляции дефектов на неравновесную критическую динамику. Впервые получено подтверждение предсказаний ренормгрупповой теории о существенности влияния дальнодействующей корреляции дефектов на поведение систем, описываемых моделями Изинга, ХУ и Гейзенберга. Показано, что поведение слабо и сильно неупорядоченных систем с протяженными дефектами принадлежит к различным классам универсальности.

9. Разработано и впервые осуществлено ренормгрупяовое описание особенностей влияния дефектов структуры и эффектов их дальнодействующей корреляции на характеристики распространения ультразвука в твердых телах при фазовых переходах второго рода с учетом как флуктуационного, так и релаксационного механизмов рассеяния при фиксированной размерности системы (I = 3 в двухпетлевом приближении с применением методов суммирования асимптотических рядов. Выделено асимптотическое поведение коэффициента поглощения и дисперсии скорости ультразвука от частоты и температуры в гидродинамической, переходной и критической областях для однородной и структурно неупорядоченных систем.

10. Впервые проведено теоретико-полевое описание фазовых превращений в сложных неупорядоченных системах с двумя взаимодействующими параметрами порядка без применения метода ^-разложения. В двухпетлевом приближении с применением техники суммирования рядов проведен ренормгрупповой анализ многовершинной модели, выделены фиксированные точки, соответствующие устойчивому мультикритическому поведению и проведено исследование влияния дефектов структуры на характер фазовых диаграмм и свойства систем в окрестности мультикритических точек. Выявлено существенное изменение областей устойчивости различных типов мультикритического поведения по сравнению с результатами ^-разложения.

11. Впервые осуществлено описание мультикритического поведения ¿-мерной системы с многокомпонентным параметром порядка с учетом влияния поверхности системы, перпендикулярной оси анизотропии кристалла, в окрестности обобщенной точки Лифшица пг-го порядка. Проведен расчет независимого поверхностного критического индекса для трехмерных систем, полученные значения находятся в хорошем согласии с результатами проведенных ранее числен-

ных исследований модели Изинга для одноосной точки Лифшица (тп = 1).

Научно-практическая значимость работы

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

Исследования, составившие основу диссертации, были поддержаны грантами РФФИ № 00-02-16455-мас, 02-02-06181-мас, 05-02-16188,10-02-00787 (2000-2012 гг.), грантами Президента РФ № МК-8738.2006.2, МК-3815.2010.2 (2006-2011 гг.), грантом совместной программы "Михаил Ломоносов" Минобрнауки РФ и Германской службы академических обменов (2004-2005 гг.) и грантами Минобрнауки РФ № А03-2.9-73 (2003-2004 гг.), № 2.1.1/930 и 2.1.1/13956 программы "Развитие научного потенциала высшей школы" (2009-2012 гг.), № 02.740.11.0541 программы "Научные и научно-педагогаческие кадры инновационной России" (2010-2012 гг.)

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

Основные положения, выносимые на защиту

1. Методика теоретико-полевого описания критического поведения структурно неупорядоченных систем, характеризуемых многовершинными моделями, с фиксированной пространственной размерностью с последовательным применением к рядам теории методов суммирования асимптотических рядов.

2. Результаты теоретического исследования влияния начальных состояний на неравновесную эволюцию однородных и структурно неупорядоченных систем в критической точке. Теоретико-полевое описание влияния на критическое поведение двумерных и трехмерных неупорядоченных систем спинстекольных эффектов, вызванных НРС.

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

4. Методика и результаты теоретического описания влияния дефектов структуры и эффектов их пространственной корреляции на характеристики распространения ультразвука в твердых

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

5. Методика расчета и результаты численного исследования равновесных критических характеристик неупорядоченных систем с применением процедуры конечноразмерного скейлинга и учета ведущих поправок к асимптотической зависимости измеряемых величин.

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

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

8. Методика и результаты теоретико-полевого описания мультикритического поведения сложных однородных и структурно неупорядоченных систем с двумя взаимодействующими параметрами порядка и исследования влияния поверхности на мультикритическое поведение в окрестности точки Лифшица.

Апробация работы

Основные результаты диссертационной работы докладывались на конференциях: Second International Workshop on Cooperative Phenomena In Condenced Matter "Quantum Phases and Phase Transitions" (Pamporovo, Bulgaria, 2001); Международной зимней школе по теоретической физике "Коуровка" (Екатеринбург, 2002, 2004, 200G, 2008, 2010); Seminar der Stipendiaten des "Michail Lomonosov" Programms (Bonn, Germany, 2005); III международной конференции "Фундаментальные проблемы физики" (Казань, 2005); IX конференции молодых ученых "Проблемы физики твердого теда-и высоких давлений" (Туапсе, 2006); 3rd International Workshop on Simulational Physics (Hangzhou, China, 2006); Молодежном семинаре по проблемам физики конденсированного состояния вещества (Екатеринбург, 2006, 2007, 2010); Международной конференции "Фазовые переходы, критические и нелинейные явления в конденсированных средах" (Махачкала, 2007, 2009, 2010); International conference "Renormalization Group & Related Topics" (Dubna, 2008); семинаре "Методы численного моделирования актуальных задач" (Москва, 2009).

Публикации По теме диссертации опубликовано 94 печатных работ, из которых 34 - статьи в ведущих международных и российских журналах, монография и 2 обзора в научных сборниках.

Структура и объем диссертации Диссертация состоит из введения, семи глав и заключения. Объем диссертации - 336 страниц машинописного текста, в том числе 85 рисунков, 48 таблиц и список цитируемой литературы из 332 наименований.

Краткое содержание работы

Во введении обоснована актуальность выбранной темы диссертационной работы и сформулированы основные цели исследований.

В первой главе, носящей обзорный характер, в краткой форме излагается содержание ря-

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

Во второй главе осуществлено теоретико-полевое описание неравновесного критического поведения однородных и неупорядоченных систем с некоррелированными дефектами структуры.

В первом параграфе проведен расчет динамического критического индекса г, определяющего критическое замедление времени релаксации системы т,*^ ~ ~ \Т — Тс\~*" вблизи температуры Тс фазового перехода второго рода (£ - корреляционная длина, и - индекс корреляционной длины) в высоких порядках теории с применением методов суммирования.

Индекс г для однородной системы определяется на основе ренормгрупповых (РГ) функций

Д(з)И7лЫ: 2 = 2 + 7л(д'), (1)

где неподвижная точка д' находится из уравнения /Цд") = 0. РГ функции ¡3(д) и у\(д) могут быть вычислены в виде рядов по д. Для размерности пространства д., близкой к четырем, координата неподвижной точки д* функции /%) принимает малые значения. В этом случае применимы методы теории возмущений по константе связи д ~ 4 — с1 и можно провести расчет критических индексов. Для реальных систем с <1 = 3,2 ряды по д являются асимптотическими и для их суммирования нужно применять специальные методы, не основанные на представлениях теории возмущений. Явный вид Р{д) для трехмерных систем в шестипетлевом приближении и для двумерных систем в четырехпетлевом приближении был получен в работах (Байкер Г.А и др., 1976,1978) ', (Ле Гуйилоу, Зин-Жустен, 1977,1978)2. В работе (Прудников В.В. и др.,1997)3 был проведен расчет динамических скейлинговых функций ■)>, для двумерной и трехмерной модели Изинга в четырехпетлевом приближении.

Асимптотическое поведение коэффициентов рядов теории подчиняется уравнению:

оо

/%) = 2>-"<Л ^ с(—а)ппьп\ (1 + 0(\/п)\ (п -> оо). (2)

п=0

Сшивка асимптотики (2) со значениями первых коэффициентов дает информацию обо всех членах ряда и позволяет приближенно восстановить функцию ¡){д) при применении процедуры суммирования асимптотических рядов. Разработаны специальные методы суммирования подобных рядов, из которых наиболее эффективными показали себя методы Паде-Бореля, Паде-Бореля-Лероя и конформного отображения.

Данные методы суммирования были применены в диссертации к расчету динамического критического индекса г для однородной модели Изинга, описывающей критическое поведение систем с размерностями г1 = 2,3. Процедура усреднения результатов с учетом аппроксимакт порядка N > 4 позволила получить следующие значения заряда д" в неподвижной точке и индекса г при применении методов Паде-Бореля (ПБ) д* = 1.7987(285), г = 2.0847(27) для

'G.A. Bilker, B.G. Nickel, et.al. // Phys. Rev. Lett. 1976. - V.36. - P. 1351; Phys. Rev. B. 1978. - V.1T. P.1365. 3J.C. Le Guillou, J, Z inn-Justin // Phys. Rev. Lett. 1977. - V.39. - P.95; Phys. Rev. B. - 1980. - V. 21. - P. 3976.

3B.B. Прудников, A.B. Иванов, A.A. Федоренко // Письмя в ЖЭТФ. 1997. - Т. 66. - С. 793.

d = 2, g' = 1.4270(79), л = 2.0171(2) для d = 3; Паде-Бореля-Лероя (ПБЛ) д' = 1.7012(163), г = 2.0757(15) для d = 2, g* = 1.4125(49), г = 2.0168(1) для d = 3; и конформного Паде-Бореля (КБ) = 1.8193, г = 2.0922 для d = 2, g' = 1.4231(21), г = 2.0372(1) для d = 3. Полученные значения заряда д' в неподвижной точке для d = 3 хорошо согласуются с результатом работы Ле Гийу, Зин-Жустина: д' = 1.416(5). Усреднение приведенных выше значений, полученных разными методами, дает итоговые величины: д* = 1.4205(38), г = 2.0213(18) для d = 3; g' = 1.7638(278), г = 2.0826(33) для d = 2.

При описании критического поведения неупорядоченных систем с замороженными некоррелированными дефектами структуры возникает дополнительная вершина взаимодействия v, определяющая эффекты взаимодействия флуктуаций параметра порядка через поле дефектов. В этом случае РГ функции ¡}д и /3„ представляются в виде двухпараметрических асимптотических рядов. Явный вид РГ функций /З3, /3„ и 7а был получен в работах (Пелисетто А., Вика-ри Е., 2000)4 и (Прудников В.В. и др., 1998)5, соответственно, в шестипетлевом и трехпетлевом приближениях. В диссертации в модифицированных на двухпараметрический случай методах Паде-Бореля и Паде-Бореля-Лероя при суммировании ряда для функции 2 +7а использовалась наилучшая для них аппроксиманта третьего порядка [2/1], а в методе конформного Паде-Бореля - аппроксиманта [1/2]. Процедура усреднения координат неподвижной точки и индекса z при N > 3 дала следующие значения:

ПБ : v' = -0.7019(119), g* = 2.2569(52), г = 2.1788(30);

ПБЛ : v' = -0.7059(36), д' = 2.2563(29), z = 2.1802(9);

КБ: и* = -0.7070(108), д" = 2.2410(245), 2 = 2.1786(29).

Полученные различными методами близкие значения координат неподвижной точки и индекса г, практически совпадающие в пределах погрешностей, говорят о достоверности полученных результатов. Усреднение данных значений, полученных разными методами, дало итоговые величины: d* = -0.7050 ± 0.0052, g* = 2.2514 ± 0.0081, 2 = 2.1792 ± 0.0013.

Во втором и третьем параграфах проведено исследование влияния неравновесных начальных состояний на эволюцию намагниченности m(t) однородных и структурно неупорядоченных ферромагнитных систем. Представлены методика и результаты исследования процессов критической релаксации из начального неравновесного состояния системы, созданного при температурах много больших критической и характеризуемого поэтому малой корреляционной длиной, в сильно коррелированное состояние при критической температуре или близких к ней.

Как показано в работе (Янсен Г. и др., 1989)6, критическая эволюция системы из начального неравновесного состояния с малой намагниченностью то = т(0) <£ 1 приводит к универсальному скейлинговому поведению для m(t) на коротковременном этапе ее критической эволюции и характеризуется аномальным увеличением намагниченности системы со временем, описываемым степенным законом m(t) ~ i"'. Вычисление показателя № было осуществлено в работе (Янсен Г. и др.,1989) методами ренормгруппы и г-разложения в двухпетлевом приближении.

Последующее исследование неравновесного поведения трехмерной модели Изинга в коротко-

"Pelissctto A., Vicari Е. // Phys. Rev. В. - 2000. - V. 62. - Р. 6393.

5Прудников В. В., Белим С. В., Иванов А. В., и др.Ц ЖЭТФ. - 1998. - Т. 114. - С. 972.

6Jansseu Н.К., Scliaub В., Schraittœaim В. // Z. Phys. В. - 1989. - V. 73. - P. 539.

Рис. 1. График неравновесной эволюции намагниченности m(t) из начального состояния с тщ = 0.03 для трехмерной модели Изинга при критической температуре Т,: = 3.49948 для р = 0.8 (Прудников П.В. и др., Phys.Rev.E, 2010).

временном режиме методами компьютерного моделирования (Ястер А. и др., 1999)7 подтвердило предсказание теории о степенном характере эволюции намагниченности, однако определенное значение показателя 9' = 0.108(2) находится в неудовлетворительном соответствии с теоретическим значением в' — 0.130, получаемым при непосредственной подстановке е = 1, или значением Н' = 0.138, получаемым при использовании метода Пзде-Бореля для суммирования имеющегося очень короткого ряда по е. В диссертации впервые был осуществлен расчет показателя корот-ковременной эволюции 0' в следующем трехпетлевом приближении РГ теории в рамках метода f-разложения и получено хорошее согласие с результатами компьютерного моделирования при применении метода Паде-Вореля для суммирования трехчленного ряда.

В соответствии с теорией скейлинга сингулярная часть потенциала Гиббса $Sm6(t, т, h, m0), определяющая состояние системы в критической области, характеризуется обобщенной однородностью относительно основных термодинамических переменных

$SmS(í,r,/i,mo) = ЬФ^г(Ь"Ч,Ьагт,Ьа'-Н,Ьатт0), времени i, приведенной температуры т, поля h и начальной намагниченности т0, здесь Ь - фактор подобия, n¿ - показатели подобия. Как следствие этого, в критической точке (г = 0, h. = 0) намагниченность m = —¿Ф/Sh характеризуется следующей временной зависимостью

m(i, то0) = r<°"+1)/"'Fm(m0rt'"'/'"). (3)

Разложение правой части в (3) по малой величине moi-""'/0' приводит к степенной зависимости m(t) ~ t-(«k+«m+i)/<tt ^ ¿0 Все «i, за исключением ат, можно связать с известными критическими индексами, описывающими поведение системы без эффектов влияния неравновесных начальных состояний. Поэтому в работе (Янсен Г. и др., 1989) были введены новые динамические критические индексы ti и в', характеризующие в коротковременном режиме неравновесное поведение функции отклика G(t,t') ~ (t/t')e и параметра порядка m(t) ~ t0' и связанные скейлинговым соотношением в' = № + (2 — г — rj)/z. Выло показано, что при временах í > tcr ~ начальный режим, характеризующийся увеличением намагниченности

m(í), переходит в традиционный режим критической релаксации системы к равновесному состоянию. Этот этап критической релаксации характеризуется степенной временной зависимостью намагниченности m ~ í-"/"* (рис. 1).

'Jaster A., MaiiivUle J., Schulke L., Zlieng В. // J. Pliys. A. - 1999. - V. 32. - P. 1395.

Для описания критического поведения однородных систем в состоянии равновесия используется модельный гамильтониан Ландау-Гинзбурга-Вильсона

ад = /'^{¿¿Ï[(^.М)2 + т4(х)] + !(£J,

(4)

где s(x) - поле n-компонентного параметра порядка, г - приведенная температура фазового перехода второго рода, g - амплитуда взаимодействия флуктуаций параметра порядка.

Пусть реализация в системе любой конфигурации параметра порядка в момент времени t определяется условием, что в начальный момент t = 0 для системы с начальной намагниченностью та распределение для поля параметра порядка s(x, 0) = so(x) характеризуется функцией распределения PJso] ~ ехр(-Яо[зо)), где

Я„Ы = Jd''i | jso(x) - m0(x)j2. (5)

Будем рассматривать случай чисто релаксационной динамики параметра порядка (модель А)8, задаваемой уравнением Ланжевена

Ôtsa(x, t) = И + С(х, t), (6)

где Л - кинетический коэффициент, Ç(x,f) - гауссовая случайная сила, моделирующая коротко-живущие возбуждения, с моментами (С»(г, ()) = 0; «Q(i, t7)) = 2ASali S{x - x')8(t - t").

В рамках теоретико-полевого описания динамики критических явлений (Баущ Р. и др., 1976)' вводится вспомогательное поле 5(х) и производящий функционал lV[/i, Л] для динамических корреляционных функций C(xi,ti,xî, t2) и функций отклика G(xt, tux2, t2) в виде

М=о Sh(xutl)Sh^,t2) hMo'

C(rb ti ,x2,t2)=----—-1

&h(xi,ti)&h(x2,t2)

W[h, h] = ln | j P(s, is) exp (-£[s, s] - tf„[.v0]) exp ^Jddx jdt Çhasa + /i„s„)j |, (7)

щонал действия С системы характеризуется выражением: = /dt /** + Х(Т - + у«, - As-a J.

Рассмотрение гауссовой составляющей функционала (8) при </ = 0 позволяет при граничном условии Дирихле (т0 = оо) получить выражения для затравочной функции отклика С?0(р, 4 -и затравочной корреляционной функции :

Со(р, г - О = ехр[—А(р2 + т)|4 - <']], С'°>(Р. О = - (') + С<"(р,4 + О, (9)

С(«>(р, (_ ,) = «РН^ + тЖ-^ ^ г + ^ = _ехР[-А(^ + г)(4 + П) (1())

При РГ анализе модели для устранения возникающих в пределе т 0 при учете взаимодействия критических флуктуаций параметра порядка расходимостей в динамических корреляционных функциях и функциях отклика была применена процедура размерной регуляризации и схема

"Hohenberg P.C. and Halpetül B.I. // Rev.Mwl.Phys. - 1977. - V. 49. - P. 435. 0Bausch R., Janssen H.K., Wagner H. // Z. Phys. B. - 1976. - V. 24,- P. 113.

Рис. 2. Диаграммы, определяющие вклад в вершинную функцию Г^. Линиям соответствует равновесный коррелятор (10), линиям со стрелкой - затравочные функции отклика. "Поверхность" t = 0 обозначена вертикальной чертой.

минимальных вычитаний с последующим переопределением параметров гамильтониана и мультипликативной перенормировкой полей в (7):

»-»г,"2», 80-» (г^0)1/2 ¿о, ,

А —> (г,/г})Ч2х, т -> д -> ( '

где е = 4 — й, /1 - размерный параметр. За счет введения в теорию начальных условий вида (5), возникает необходимость в перенормировке функции отклика (¡(р, ¿).?0(—р, 0)), задающей влияние начальных состояний системы. Поправочные слагаемые в собственно-энергетической части функции отклика, возникающие за счет эффектов взаимодействия флуктуаций параметра порядка, характеризуются приводимыми динамическими диаграммами Фейнмана, поскольку их вычисление осуществляется с использованием коррелятора (9), не обладающим свойством трансляционной инвариантности во времени. В работе (Г.Янсен и др.,1989) было введено следующее представление для данной функции отклика:

I

с10,(Р,0 = («СМ)«.(-Р,0)> = 1<И'С1Л(р,1,1')гЧ1{р,1')ы. (12)

о

Одночастичная вершинная функция Г''£,(р, {)[;„( с одной вставкой поля ё0 в трехпетлевом приближении описывается диаграммами, характеризуемыми требованием, чтобы они содержали хотя бы один коррелятор Множитель <5и(р, М') определяется равновесной составляющей коррелятора С'¡р в (9). Отметим, что он отличен от равновесной функции отклика по причине интегрирования в (12) по времени от начального момента с t = 0 вместо ( = -оо. Между ними можно установить функциональную связь, вводя вместо функционала (5) функционал На\зо\ (4) с новой вершиной взаимодействия в функционале действия (8)

^ Лй /л (^«о,,) (®од)2 ■ (13)

За счет усреднения по начальным полям возникает дополнительная вершинная функция локализованная на "поверхности " Ь — 0. Подобно (12), имеет место следующее выражение:

Решив интегральное уравнение ¿(4—4') = Л" К {у, I", О!"!/]' ('/> '•)[»(<")! в каждом порядке теории было найдено его ядро 4", ('), флуктуационные поправки к которому для данной модели возникают только в третьем порядке теории (рис.2). Используя выражения (12) и (14), а также осуществляя перенормировку полей в соответствии с (11), можно установить следующее выражение для определения перенормировочной константы Zg:

Таблица 1. Результаты расчета индекса в1 для моделей Изинга и ХУ

2-петлевое приближение 3-петлевое приближение Подст. е = 1 Метод ПБ Подст. с — 1 Метод ПВ Комп. моделирование

Модель

Модель Изинга XY-модель

0.130 0.138 0.0791 0.1078(22)

0.154 0.170 0.0983 0.1289(23)

0.108(2)п 0.144(10p

7-1/2

ОО

У

dte'^ltifi(q = 0, i)[äol,

(15)

из требования устранения полюсов по е в данном выражении в каждом порядке теории, т.е. конечности выражения (15) в пределе е 0. Последовательная реализация изложенной процедуры и расчет диаграмм с использованием метода £-разложения позволили вычислить константу перенормировки 2о в трехпетлевом приближении:

2 "+2 + 216

(п + 5)(п + 6) +

а

1,

(9In2 — 19)(п + 2) + -(61п2 - 7)(п + 8)

18679(п + 3.13882)£!

(f)3- 0«)

Для коротковременого режима неравновесной критической релаксации принципиально новой является лишь РГ функция 70, которая в трехпетлевом приближении, как показали наши расчеты, принимает следующее выражение

(17)

п + 2 7о ---—3

|l+ ^ln 2 - 0 д- 0.0988989 (п+ 3.13882) g2j +0(г/).

При определении значения индекса в' были использованы результаты работы (Кляйнерт Г. и др., 1991)10 по определению координаты устойчивой неподвижной точки </ и полученное нами выражение для динамического критического индекса г:

г = 2+-

(п + 2)е ; 4(п + 8)

1 +

6е

(п + 8)2

("Иг?

>ит]

^п + 3 + (п + 8) In

1 + е

/6(Зп + 14)_0438481Л

V (« + 8)2 )

2 \ 3 7(л + 8)2 В итоге выражение для критического индекса в' принимает следующий вид:

(18)

э-

7.2985ег (п + 8)4

(п3 + 17.312п2 + 153.267п + 383.552)

(19)

К ряду (19) для индекса & был применен метод суммирования Паде-Бореля (ПБ). Полученные результаты представлены в табл. 1.

Проведенное сопоставление значений & с результатами компьютерного моделирования моделей Изинга (Ястер А. и др., 1999)11 и ХУ (Прудников П.В. и др., 2007) 12 наглядно демонстрирует, что результаты расчета в трехпетлевом приближении находятся в гораздо лучшем согласии с результатами компьютерного моделирования, чем результаты двухпетлевого приближения.

Третий параграф главы посвящен исследованию одновременного влияния эффектов нарушения пространственной трансляционной симметрии системы, создаваемых присутствием дефектов структуры, и эффектов нарушения временной трансляционной симметрии, обусловлен-

l0Kleinert Н., Neu J„ Schulte-ftohlinde U., et. al. // Phys.Lett. В. - 1991. - V. 272. - P. 39.

"Jaster A., Mainvilte J., Schulke L„ Zheng B. // J.Phys.A: Math.Geii. - 1999. - Vol. 32. - P. 1395.

"Прудников П.В., Гергертд E.A. и др. // Вестник ОмГУ. - 2007. - Вып. 4. - С. 32.

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

Для описания критического поведения структурно неупорядоченных изинговских систем в состоянии равновесия используется модельный гамильтониан Ландау-Гинзбурга-Вильсона

HclM = [(Vs(x))2 + t(x).s2(x)] + |/(x)j, (20)

где s(x) - поле параметра порядка, т(х) - приведенная случайная локальная температура фазового перехода второго рода, g - амплитуда взаимодействия флуктуаций параметра порядка. Случайную температуру можно представить в виде т(х) = г + V(x), где т - приведенная температура фазового перехода для однородной системы, a V(x) - потенциал случайного поля дефектов. Пространственное распределение системы замороженных некоррелированных дефектов характеризуется гауссовым распределением P\V\ и полностью определяется значениями первого и второго моментов для случайных величин V(x): ^

«V(x)»=0, <{K(x)V(y)» = wr(x-y). (21)

где t> - положительная константа, пропорциональная концентрации дефектов и квадрату величины их потенциала.

В соответствии с теоретико-полевым подходом вводится производящий функционал W[h, h]

оо

W[/i,h] — \n(J P(s,is)P[V]exp(—£v[s, s, V] — H0[su\) exp^J J dt'Jii + |, (22)

о

в котором функционал действия s> V] неупорядоченной системы характеризуется выражением: у ,

_

£v[s, S, V] = Jdt Jddx i

о

(23)

В (22) можно осуществить усреднение но случайным полям дефектов У(х)

Ур[К]ехр(-£и[5,-5,^])=ехр(-£[5,.;]) (24)

и получить функционал действия независящий от случайных полей У(х) и являющийся

трансляционно инвариантным:

ТО öo

£[5,5)= fétjd*x~s | + A(r-V2)s+^S3-AI -v^lfàtJdAx

(25)

Вычисление всех констант перенормировкй при фиксированной размерности системы <1— 3, кроме проведено в работе (Прудников В.В. и др., 1998)13. В диссертации представлен расчет 2о Для неупорядоченных систем в двухпетлевом приближении теории при <1 = 3.

За счет введения в теорию начальных условий вида (5), возникает необходимость в перенормировке функции отклика 4)50(—р,0)), задающей влияние начальных"состояний системы. По аналогии с однородными системами для описания поправочных слагаемых в собственно-энергетической части функции отклика за счет эффектов взаимодействия флуктуаций параметра порядка вводится в соответствии с (12) одночастичная вершинная функция Г'дДр, 0[»„]

"Прудников В.В., Белим C.B., Иванов A.B. и др. // ЖЭТФ. - 1998. - Т. 114. - С. 972.

lalQQ Le- кЮ-ЦЭ.

Рис. 3. Диаграммы, определяющие вклад в вершинные функции Г^д. Линиям соответствуют корреляторы линиям со стрелкой - функции отклика £?о, вершинам д - жирная точка, v -волнистая линия.

с одной вставкой поля so- Она в двухпетлевом приближении описывается диаграммами, представленными на рис. 3 и характеризуемыми требованием, чтобы они содержали хотя бы один коррелятор Cq1'. Множитель G\t\(p,t,t') в (18) определяется равновесной составляющей коррелятора Ср"' в (9). Отметим, что он отличен от равновесной функции отклика G*¡*j\p,t — t'} по причине интегрирования в (12) по времени от начального момента с ¿ = 0 вместо t = -ca. Между ними можно установить функциональную связь, если вместо функционала (5) ввести функционал Hc¿[so] (20) с новыми вершинами взаимодействия в функционале действия (25)

оо 2

Jdt Jd'x - Vj У dt já'x (SoS„)j . (26)

За счет усреднения по начальным полям возникает дополнительная вершинная функция Г'^', локализованная на "поверхности" i = 0. От первого слагаемого в (26), как показано в (Прудников П.В. и др., 2008)14, флуктуационные поправки в Г^' возникают только начиная с трехпет-левого приближения, в то время как за счет второго слагаемого в (2G), обусловленного влиянием структурных дефектов, флуктуационные поправки в возникают уже начиная с двухпетле-вого приближения (рис. 4). Для определения константы Zo зададим следующее нормировочное соотношение :

V2lto(í> = 0, ш/2А = ¡i2) = 1, (27)

где rf0(p, ш) j Фурье-образ перекормироеанной одночастичной вершинной функции Г^о(p,t), рассчитываемой в удобной для нормировки точке с импульсом р = 0 и частотой iu/2\ — /А Последовательная реализация изложенной процедуры и расчет диаграмм при d — 3 позволили вычислить константу перенормировки Za и РГ функцию 7о в двухпетлевом приближении:

Z0 = l + \gR + 0.012682gl - 0.608932дякд, То = ~\qr + 0.457127g2n - 0.614995дК"я, (28)

где уя и vR - перенормированные константы связи. Неподвижная тючка РГ преобразований (g",v') определяется из уравнений: ß9(g',v') = 0, А,(д*,к*) = 0. В результате, для неупорядоченной модели Изинга были получены следующие выражения для динамических критических индексов: г = 2- 0.25v' + 0.00840(g*)2 + 0.030862g\>* + 0.053240(и')2, (29)

0 = lg' - 0.1142817(g*)2 + 0.174582g V,

ff' = ig- + 0.125u* - 0.123968(g*)2 + 0.14680608g*v* - 0.0156245(u*)2.

6

"Прудников П.В., Прудников B.B., Калашников H.A., Циркин C.C. // ЖЭТФ. - 2008. - Т. 133. - С. 1251.

Для дальнейших вычислений нами были использованы значения констант связи в неподвижной точке с д' = 2.2514(42), v' = 0.7049(13), определенные в работе (Прудников П.В. и др.,2006) и представленные во второй главе диссертации. Для получения физически разумных значений критических индексов к рядам (29) были применены методы суммирования Паде-Бореля, Паде-Бореля-Лероя и конформного отображения. В результате были вычислены следующие значения г и в': z = 2.2009867, в' = 0.091898 (ПБ); г = 2.198340, в' = 0.120284 (ПБЛ при значении параметра Ь = 2.221426 (Прудников П.В. и др., 2006)); г = 2.205156, (f = 0.104441 (КБ).

Итоговые средние значения критических индексов:

г = 2.202(2); в' = 0.106(8). Отметим, что данные значения динамического критического индекса г, вычисленные в двух-петлевом приближении, превышают его значения (со средним г = 2.1792(13)), рассчитанными в (Прудников П.В. и др., 2006; Прудников В.В. и др., 1998) в трехпетлевом приближении теории с применением различных методов суммирования. Это служит отражением наблюдаемых в рядах теории критических явлений осциллирующих в каждом последующем порядке теории отклонений суммы ряда от его асимптотического значения, уменьшающихся с ростом порядка теории за счет применения обоснованных методов суммирования. Полученные значения критического индекса z находятся в прекрасном согласии со значением z = 2.18(10), выявленным экспериментально в (Розов Н. и др., 1992)15 методом меесбауэровской спектроскопии в результате высоко прецезионного измерения динамического уширения мессбауэровских линий в слабо разбавленном образце изинговского магнетика Feo.9Zno.1F2.

Третья глава посвящена компьютерному моделированию критического поведения трехмерных спиновых систем с некоррелированными дефектами структуры.

В первом параграфе главы осуществлено компьютерное моделирование неупорядоченной трехмерной модели Изинга, определены функциональная форма скейлинговых функций и значения критических температур и критических индексов.

Рассматривается модель неупорядоченной спиновой системы в виде кубической решетки с линейным размером L и наложенными граничными условиями. Микроскопический гамильтониан неупорядоченной модели Изинга записывается в виде

H=-JY,P.vASj, (30)

<У>

где J > 0 - короткодействующее обменное взаимодействие между закрепленными в узлах решетки спинами Si, принимающими значения ±1. Немагнитные атомы примеси образуют пустые узлы. Числа заполнения р, принимают значения 0 или 1 и описываются функцией распределения

P(Pi) = (l-p)S(p,)+pS(l-pt) (31)

с р = 1 — с, где с - концентрация атомов примеси. Примесь равномерно распределяется по всей системе, и при моделировании ее положение фиксируется для отдельной примесной конфигу-

15Rosov N., Hoheuemser С., Eibschutz М. // Pliys. Rev. В. - 1992. - V. 46. - P. 3452.

рации. В процессе исследования были рассмотрены неупорядоченные системы со значениями спиновых концентраций р = 0.95,0.80,0.60,0.50.

Для снижения влияния эффектов критического замедления и корреляции различных спиновых конфигураций в работе был применен наиболее эффективный в этом смысле однокла-стерный алгоритм Вольфа. За один шаг Монте-Карло на спин (МСБ/я) принималось от 10 до 20 переворотов кластера Вольфа в зависимости от линейного размера моделируемой решетки, спиновой концентрации системы и близости температуры к критической точке. Процедуре установления термодинамического равновесия в' системе отводилось 104 МСЭ/з и 105 МСЗ/я отводилось на статистическое усреднение моделируемых величин при заданной примесной конфигурации. Для определения средних значений термодинамических и корреляционных функций наряду со статистическим усреднением применялось усреднение по различным примесным конфигурациям: для систем с р = 0.95 усреднение проводилось по 3000 образцов, для р = 0.80 по 5000 образцов, для р — 0.60,0.50 по 10000 образцов.

В процессе моделирования различных спиновых систем на решетках с линейным размером Ь осуществлялся расчет корреляционной длины и восприимчивости XI в соответствии со следующими соотношениями: _

г'-го узла решетки, (...) означает статистическое усреднение по шагам Монте-Карло, а [...] -усреднение по примесным конфигурациям. Была определена температурная зависимость £.l{T) и Xl(T) В интервале т = (Т - Тс)/Тс ~ 5 • 10~4 - 10~2 для образцов с р = 0.95 и линейными размерами в интервале L ~ 20 — 400, для образцов с остальными спиновыми концентрациями температуры и линейные размеры решеток выбирались в интервалах т ~ III'1 — К)"2 и L = 20 - 300. Зн£1чение £гпах для каждой температуры ограничивалось решеткой, для которой корреляционная длина и восприимчивость выходили на свои асимптотические значения.

Т.к. аномалии характеристик фазового перехода второго рода проявляются лишь в термодинамическом пределе, то для определения асимптотических значений термодинамических величин А(Т) по их значениям А[.(Т), определяемым на конечных решетках, применяются различные методы конечноразмерного скейлинга. В диссертации был применен метод, предложенный в (Ким Ж.К. и др., 1996)ш для анализа результатов моделирования критического поведения двумерной и трехмерной однородных моделей Изинга.

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

Это приводит к соотношению, позволяющему определять асимптотическое значение А{Т) на основе измеряемых значений Al и скейлинговой функции Qa от хс(Т) — {¿(Т)/Х/. Скейлинго-вая функция Qa(xl) в критической области должна удовлетворять асимптотическому условию lim Qa(x) 1. С учетом этого скейлинговые функций для восприимчивости и корреляционной

"Kim J.K., de Souza A.J., Landau D.P. // Pliys. Rev. E. - 1996. - V. 54. - P. 2291.

Al{T) = A(T)Qa(Xl(T)), Xl(T) = ÎL(T)/L.

(33)

OtWu

p«0.95 - pal

p«0.BO • pci

P>0.a0 * pd

P'0-50 * pal

"HllBj

Рис. 5. Усредненные скейлинговые функции для корреляционной длины (а) и восприимчивости (6), полученные с использованием полиномиальной от х (символы) и от ехр(—Х/х) (сплошные линии) аппроксимаций

длины выбирались в виде полиномиальной зависимости от х: Qa (х) = 1 + f^a'", и полиномиальной функции от экспоненты е-1/1: Qa (х) = 1 + <-п('~"'х, с подбираемыми по методу наименьших квадратов коэффициентами с,, для каждой температуры Т. На рис. 5 представлены усредненные по различным температурам скейлинговые функции для корреляционной длины (рис. 5(a)) и восприимчивости (рис. 5(6)) для различных полиномиальных аппроксимаций и спиновых концентраций р. Данные скейлинговые функции наглядно демонстрируют существование двух классов универсального критического поведения для разбавленной модели Изинга с различным характером их поведения для слабо (]> = 0.95,0.80) и сильно (р = 0.60,0.50) неупорядоченных систем. На основе усредненных скейлинговых функций был проведен расчет асимптотических значений ((Т) и х(Г), отвечающих Т = 4.265 ч- 4.335 для р = 0.95, Т = 3.51 н- 3.57 для р = 0.80, Т = 2.430 -f 2.460 для р = 0.60 и Г = 1.851 1.874 для р = 0.50.

Асимптотический критический индекс S термодинамической величины А(т) описывается выражением

А(т) = А±\т\-', (34)

МО 1п|т]

где А+ и - критические амплитуды выше и ниже критической точки, соответственно. Степенной закон типа (34) является точным лишь в пределе т -> 0. Для расчета критических индексов в промежуточном неасимптотическом режиме необходимо вводить дополнительные поправочные слагаемые к степенному закону (34). В соответствии с разложением Вегнера

А(т) = (Ао + AiT"" + Агт2и" + ...) т-' (т > 0), (35)

где Ai - неуниверсальные амплитуды, ш - критический индекс поправки к скейлингу. В диссертации для расчета характеристик учитывалась только первая ведущая поправка к асимптотическому поведению для корреляционной длины и восприимчивости:

«т) = г-(4 + Л1г"'), Х(т) = т->(А% + А1 г"), а = ик (36)

Был проведен расчет значений критических индексов к, 7 и о, « также критических температур, используя метод наименьших квадратов для наилучшей аппроксимации асимптотических значений f (Г) и %{Т) выражениями (36). В табл. 2 представлены усредненные по применяемым аппроксимациям значения критических характеристик для различных спиновых концентраций р. Видно, что критические индексы образуют две группы, близкие по значениям в пределах погрешностей вычисления: одна группа с р = 0.95,0.80, т. е. для слабо неупорядоченных систем со спиновыми концентрациями р, большими порога примесной перколяции Pirap (для кубиче-

Таблица 2. Усредненные значения критических характеристик

р V 7 Е7 Тс

р= =0.95 0.6909(33) 1.3385(54) 0.137(56) 4.26267(4)

р= =0.80 0.6956(29) 1.3447(40) 0.178(87) 3.49948(18)

р= =0.60 0.7253(36) 1.4154(107) 0.199(103) 2.42413(9)

р= =0.50 0.7370(33) 1.4283(33) 0.207(100) 1.84509(6)

ских систем ршф = 0.69), другая с р = 0.60,0.50, т. е. для сильно неупорядоченных систем с Pc < Р < Pimp, где Pc - порог спиновой перколяции (для кубических систем рс = 0.31), когда в системе существуют два взаимопроникающих протекающих кластера - спиновый и примесный. Фрактальные эффекты этих двух пронизывающих друг друга кластеров могут явиться причиной изменения характера критического поведения для сильно неупорядоченных систем.

В качестве итоговых можно рассматривать усредненные значения критических индексов и = 0.693(5), 7 = 1.342(7), та = 0.157(92) для слабо неупорядоченных систем и и = 0.731(11), 7 = 1.422(12), га = 0.203(106) для сильно неупорядоченных систем. Отметим, что значения индексов для слабо неупорядоченных систем хорошо соотносятся со значениями f = 0.678(10), 7 = 1.330(17), £¡7 = 0.170(71) (и = 0.25(10)), рассчитанными в (Пелисегго А., Викари Э., 2000) РГ методами в шестипетлевом приближении.

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

Во втором параграфе главы представлено численное исследование влияния неравновесных начальных состояний на эволюцию намагниченности m(t) ферромагнитной структурно неупорядоченной системы в критической точке.

В соответствии с работами (Янсен Г., 1989; Ястер А., 1999) при описании неравновесного критического поведения спиновой системы с начальным значением намагниченности гпц выражение для к-го момента намагниченности системы уже после микроскопически малого щюмежутка времени хаотизации trnLC может быть описано следующим скейлинговым выражением:

m(*> (i, г, L, rn0) = b-W'mM (Ь~Ч, Ь^т, b~lL, b™m0) , (37)

где t - время, т — (Т ~ Тс)/Тс, Ь - произвольный масштабный фактор, L - линейный размер решетки, ß, v, z - известные критические индексы, хо - новый независимый критический индекс, задающий масштабную размерность начального значения намагниченности то- На ранней стадии эволюции системы корреляционная длина еще достаточна мала и конечность размера моделируемой системы оказывается не существенной. Полагая в (37) Ь = г1/*, дня первого момента намагниченности (к = 3) и малой величины mot1/' получаем следующее выражение:

m(i, т, ma) ~ т0/F(tl'"T, f^'nu,) = m</'( 1 + at^'r) + 0(r2, mj),

где d — (хо — ß/v)/z. Для т -» 0 и достаточно малых t получаем асимптотическое поведение m(t) ~ te . Временной интервал увеличения намагниченности t„ ~ заметно растет с

17Фольк Р., Головач Ю„ Яворский Т. // УФН. - 2003. - Т. 173. - С. 175.

уменьшением ту. С течением времени коротковременная динамика увеличения параметра порядка сменяется на привычную долговременную динамику уменьшения параметра порядка со временем по степенному закону т(4) ~

Для численного определения показателя в' была рассмотрена неупорядоченная модель Изин-га (30), задаваемая на кубической решетке с линейным размером Ь = 128 и наложенными граничными условиями. Для моделирования спиновых конфигураций в системе был применен алгоритм Метрополиса. Алгоритм Метрололиса, реализующий динамику односпиновых переворотов, наилучшим образом соответствует релаксационной модели в классификации Гальперина-Хоэнберга и позволяет провести сравнение результатов моделированиия и проведенного во второй главе ренормгруппового описания.

В данном параграфе представлены результаты численного исследования неравновесного поведения слабо неупорядоченной системы со спиновой концентрацией р = 0.80 и сильно неупорядоченной с р = 0.60. Осуществлялось численное определение временной зависимости к-го момента намагниченности в виде

N.

(38)

А{1) =

о)

(39)

где ^ = р1? - число спинов в решетке. Проводилось усреднение вычисляемых величин по 4000 различным примесным конфигурациям с 25 прогонками для каждой примесной конфигурации. Для независимого вычисления динамических критических индексов в' и г, а также отношения статических критических индексов р/и на каждом этапе эволюции системы наряду с намагниченностью (к — 1) определялись автокорреляционная функция

N.

и второй момент намагниченности т'2'(<). Их скейлинговый анализ показывает (Ястер А. и др., 1999), что при то = 0 и Т = Тс данные величины характеризуются степенной зависимостью от

времени Л(е)~г-, (40)

где са — <1/г — в', {'2 = {(I — 2/3/и)/г, <1 - размерность системы.

Моделирование осуществлялось при критических температурах Тс = 3.49948(18) и 2.42413(9) для систем с р = 0.80 и 0.60, соответственно, определенными при Монте-Карло исследованиях неупорядоченной модели Изинга в равновесном состоянии (Прудников П.В. и др., 2007). Временное поведение намагниченности с начальными значениями то = 0.01, 0.02 и 0.03 исследовалось на временах до 1000 МСв/й. На рис. 6а представлены в двойном логарифмическом масштабе временные зависимости для намагниченности системы т(1) для различных спиновых концентраций при моделировании из начального состояния с тц = 0.01 для случая р = 0.8 и с т0 = 0.005 для р = 0.6. Они позволяют определять показатели в'(т0) и их асимптотическое значение в' (тд —V 0) на основе линейной аппроксимации значений в (от0) при т0 —> 0.

На рис. 6 для данных систем, стартующих из неравновесного начального состояния с близким к нулю значением то = 0.0001, представлены временные зависимости для второго момента намагниченности (Рис. 66) и автокорреляционной функции (Рис. 6в), изображенные в двойном логарифмическом масштабе. Анализ данных зависимостей позволяет определять значения

Рис. 6. Временные зависимости для намагниченности m(t) при эволюции из начального состояния

то (а), второго момента m(2'(t) (6) и автокорреляционной функции A(t) (в), показателей с„ и c¡ в соответствии с (40). На временных зависимостях данных величин могут быть выделены по два линейных участка: для временных интервалов в среднем от 10 до 80 MCS/s и от 300 до 800 MCS/s. Это связано с наблюдаемым уже при моделировании структурно неупорядоченных систем с линейными дефектами (Прудников П.В. и др., 2007) 13 явлением кроссовера, т.е. перехода от критического поведения, характерного для однородных систем, к поведению, характеризуемого влиянием дефектов структуры. Был проведен расчет показателей для каждого линейного участка исследуемых величин.

При анализе был также осуществлен учет поправок к асимптотической зависимости измеряемых величин за счет влияния конечности моделируемых систем и неточности в определении их критической температуры. Учет данных поправок к скейлингу позволяет получать корректные значения критических индексов. Для этого были применены следующие выражения для временной зависимости наблюдаемых величин:

X(t)~t\l + Axt-"'*), (41)

где Ах - неуниверсальные амплитуды, ы является хорошо известным критическим индексом поправки к скейлингу, а показатель S = в' в случае X = m(f), á = —с„ в случае X s A(t) и ó = с2 в случае X = m'2'(i). Для расчета значений критических индексов Ö , с„, с2 и w/z на временном интервале, соответствующем влиянию структурного беспорядка, был использован метод наименьших квадратов для наилучшей аппроксимации значений m(í), A{t) и m(2,(í) выражением (41). Процедура заключалась а следующем: 1) временной интервал проявления влияния дефектов структуры разбивался на всевозможные участки At, начиная от участков с At = 50 до участков с At — 550; 2) на каждом из участков At осуществлялось определение значения показателя ó при фиксированном значении oj/z; 3) найденные значения <? усреднялись по выбранным участкам с определением среднего значения (á) и погрешности аппроксимации Д<5; 4) показатель u/z определялся из условия минимальности значений относительных погрешностей проведенных аппроксимаций.

Наряду с аппроксимационной погрешностью AS для показателей S определялась их статистическая погрешность. Для этого общее количество используемых для усреднения примесных конфигураций делилось на 4 группы. Для каждой из групп вычислялись показатели в, с„ и с'2, а затем определялись отклонения от показателей, найденных ири использовании усредненных по общему числу примесных конфигураций значений т, А и т'2>. Полученные значения критических индексов представлены в табл. 3.

В случае исследования влияния неравновесного состояния с начальным значением намагни-

I8Prudnikov P.V., Prudmkov V.V., Zheng В. еЫ. // Progr. Tlieor. Pliys. - 2007. - V. 117. - Р. 973.

"Ч,

(а)

"'■'¡а-»-

[г-Ю.95

мм

"да

//Orí» 1«

UWM

1

■í!' (б)

dr]nM(t,T) р-0.95

>>р-о-зо

/

А

а"

t (MCS/s)""

Рис. 7. Временные зависимости для намагниченности m(t) (а), кумулянта ¡7гШ (6) и логарифмической

производной намагниченности 91nm(t)|T_0 (в) для различных спиновых концентраций р. ченности тпо = 1 не возникает зависимости характеристик неравновесного критического поведения от нового динамического показателя В . Исследования показали, что динамический процесс, начинающийся с полностью упорядоченного состояния {то = 1), более предпочтителен из-за меньшего влияния флуктуаций на результаты. В критической области т 1 для временного поведения намагниченности справедлива следующая зависимость

m(¿,т) = Г^тЦ^'т) = t'V" (1 + ««'/"г + 0(т2)), (42)

где в пределе г -> 0, оно приобретает вид:

m(¿) ~ Г*'". (43)

Представляя (42) в виде ln m(t, т) = {-fi/vz) ln t + ln m(l, t^"zr) и дифференцируя по т можно получить соотношение для логарифмической производной намагниченности

дгЫт&т)],^-?'". (44)

Другой характеристической величиной, определяющей неравновесное критическое поведение является кумулянт Биндера, характеризуемый выражением

U2(t) = m'2' /(m)2 - 1. (45)

Размерный анализ показывает, что в критической точке поведение кумулянта Биндера описывается степенным законом u2(t) ~ td/z. Численное определение намагниченности, ее логарифмической производной и кумулянта Биндера позволяет рассчитать динамический индекс z и статические индексы fi и и.

Осуществлялось моделирование кубических решеток с размерами L = 128 при критических температурах Тс = 4.26267,3.49948, 2.42413 и 1.84509 для систем со спиновыми концентрациями р — 0.95, 0.80, 0.60 и 0.50, соответственно. Временное поведение намагниченности и кумулянта Биндера исследовалось для слабо неупорядоченных систем на временах до 1200 MCS, а для сильно неупорядоченных систем на временах до 3000 MCS. Для систем с р = 0.95 проводилось усреднение вычисляемых величин по 6000 различным примесным конфигурациям, с р = 0.80 - по 50000, р = 0.60 - по 10000 и р — 0.50 - по 7000 примесным конфигурациям. Для вычисления логарифмической производной дт ln m осуществлялся расчет намагниченности для двух температур, смещенных относительно Тс на интервал AT = ±0.005.

На рис. 7 представлены в двойном логарифмическом масштабе усредненные временные зависимости для намагниченности m(t) (рис. 7а), кумулянта U2{t) (рис. 7Ь) и логарифмической производной намагниченности д\пm(t) для спиновых концентраций р = 0.95, 0.80, 0.60, 0.50.

Таблица 3. Значения критических индексов и сравнение их с результатами компьютерного модслиро

вания (МС), теоретико-полевого описания (FTM) и эксперимента (схр)

z в' p/ir v 0 и

р = 0.95, то = 1 2.185(25) 0.533(13) 0.668(14) 0.356(6) 0.369(96)

р = 0.80, m0 = 1 2.208(32) 0.508(17) 0.685(21) 0.348(11) 0.404(110)

р = 0.80, то 1 2.191(21) 0.127(16)0.504(14) 0.256(55)

р = 0.60, то ~ 1 2.525(15) 0.496(15) 0.682(13) 0.339(12) 0.286(10)

р -- 0.50, то = 1 2.664(37) Ç.443(7) 0.760(30) 0.337(19) 0.242(4)

р = 0.60, то « 1 2.589(52) 0.167(18)0.485(18) 0.461(55)

Pelissetto, Vicati, 2000; (FTM) 0.515(15) 0.G78(10) 0.349(5) 0.25(10)

Prudnikov, et .al., 2006; (FTM) 2.1792(13)

Rosov, et.al., 1988,1992 FepZni-pF2 p = 0.9; (exp) 2.18(10) 0.350(9)

Slanic, et.al., 1999, FepZm_ pF2 p = 0.93 (exp) 0.70(2)

Prudnikov, Vakilov, 1992, = 0.95, 2.19(7)

= 0.80, 2.20(8)

V = 0.60, 2.58(9)

P = 0.40; (MC) 2.65(12)

Wiseman, Domany, 1998, p = 0.80, 0.505(2) 0.682(2)

P = 0.60; (MC) 0.437(21)0.717(6)

Ballesteros, et.al., 1998, p = 0.90 -r 0.40; (MC) 0.519(8) 0.684(5) 0.355(3) 0.370(63)

Parisi, et.al., 1999, p = 0.90 -r 0.40 (MC) 2.02(7) 0.50(13)

Murtazaev, et.al., 2004, p = 0.95, 0.646(2) 0.306(3)

P = 0.9, 0.664(3) 0.308(3)

P = O.S, 0.683(4) 0.310(3)

P = 0.6 (MC) 0.725(6) 0.349(4)

Schehr, Paul, 2006; (MC) 0.10(2)

Hasenbusch, et.al., 2007, p = 0.8; (MC) 2.35(2)

Prudnikov, et.al., 2007, p = 0.95 -b 0.80, 0.532(12)0.693(5) 0.26(13)

P = 0.60 4- 0.50; (MC) 0.524(13)0.731(11) 0.28(15)

В слабо неупорядоченных системах с р — 0.95; 0.80 в отличие от поведения однородных си-

стем может быть выявлено два универсальных динамических режима со степенным временным изменением m(t), U2(t) и ôlnm(i), а именно: на раннем временном интервале t. = [20,200] реализуется поведение, соответствующее поведению однородной системы, определяемое индексом 2 = 2.03(1), а лишь затем, проходя через режим кроссоверного поведения, реализуется режим поведения неупорядоченной системы. Для сильно неупорядоченных систем динамического режима с характеристиками однородной системы не наблюдается. При определении показателей был проведен учет ведущих поправок к асимптотической зависимости измеряемых величин, согласно (41). В табл. 3 приведены полученные итоговые значения критических показателей.

На основе проведенных в данном параграфе численных исследований (табл. 3) можно сделать следующие выводы: слабо и сильно неупорядоченные системы принадлежат к различным классам универсальности критического поведения с несовпадающими в пределах статистических погрешностей значениями критических индексов в' и г; полученные значения индексов для слабо неупорядоченных систем находятся в хорошем согласии с результатами теоретико-полевого описания и компьютерного моделирования другими методами, а также согласуются с результатами экспериментальных исследований изинговских неупорядоченных магнетиков.

В четвертой главе осуществлено ренормгрупповое описание модели слабо неупорядоченной системы с введенным потенциалом взаимодействия четвертого порядка по флуктуациям

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

/f I п т 1 п т . п т Ч

л {2 Е Е+г Е Е +J Е Е ?<4«(*)i2K(*)i2 . (46)

I ¿=1 ¡=1 «=t .j=la,4=l J

где индекс i определяет число компонент параметра порядка ф{х), индекс а нумерует реплики (образы) однородной составляющей в исходном гамильтониане неупорядоченной системы, матрица (;„/, задает эффективное взаимодействие флуктуаций (т х га) - компонентного параметра порядка через поле дефектов. Данная статистическая модель термодинамически эквивалентна исходной неупорядоченной модели в пределе т —> 0. Наличие дефектов структуры приводит к реализации в системе большого числа локальных минимумов энергии и матрица даь уже не является реплично-симметричной, а имеет структуру НРС Паризи (Доценко Вик. и др., 1995)19. Так, в пределе т О, матрица /)„4 характеризуется диагональными элементами д и недиагональными элементами, задаваемыми функцией д[х), которая определена на интервале 0 < х < 1: дл (д,д(я))- Реплично-симметричной ситуации соответствует д(х) = const, не зависящая от х. В работе (Доценко Вик. и др., 1995) была выявлена ступенчатая структура функции д{х). Мы ограничились рассмотрением функции д(х) одноступенчатого вида', д(х) — да для 0 < х < хо, и д(х) = д1 для .(•(> < х < 1, где координата ступеньки 0 < х0 < 1 остается произвольным параметром, который не меняется при масштабных преобразованиях. В результате РГ преобразования реп личного гамильтониана с НРС задаются тремя параметрами g,go,Vi-

РГ описание модели, задаваемой репличным гамильтонианом (46), нами было осуществлено в рамках теоретико-полевого подхода в двухпетлевом приближении как непосредственно для т])ехмерного случая, так и для 3 < d < 4. Возможные типы критического поведения и их устойчивость определяются коэффициентами fio,Si) РГ уравнений для вершинных частей неприводимых функций Грина. Для их определения был применен метод диаграмм Фейнмана и процедура перенормировки. Искомые /J-функции были получены в виде степенных рядов по параметрам g,go,gi. Для их пересуммирования был применен метод Паде-Бореля.

Природа критического поведения определяется существованием устойчивой фиксированной точки (ФТ), удовлетворяющей системе уравнений ft(a*\So>!/!) = 0, {i = 1,2,3). В результате решения данной системы для значений числа компонент параметра порядка п = 1,2,3 было выделено три типа нетривиальных ФТ в представляющей физический интерес области значений -параметров д',ц'„д] > 0. Так, фиксированная точка первого типа ФТ1 с д' ф 0,</5 = </? = 0 соответствует критическому поведению однородной системы, ФТ2 с д' ф 0, д'й ~ д\ ф 0 -критическому поведению неупорядоченной системы с репличной симметрией, а ФТЗ с д' ф 0, Но — д\ ф0 - критическому поведению неупорядоченной системы с НРС. При этом значения параметров в ФТЗ зависят от координаты ступеньки х0.

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

"Doteraiko Vik.S., Felilmau D.E. /,/ J. Pliys. A. - 1095. - V. 28. - P. 5183.

сводится к условию, чтобы действительные части собственных значений А,- матрицы устойчивости Bij — d/3,(g',ffg,gl) /dgj были положительны. Для каждого типа ФТ был проведен численный расчет собственных значений матрицы устойчивости при изменении размерности системы от 3 до 4. Пороговые размерности системы, разделяющие области различных типов устойчивого критического поведения, определялись но изменению знака хотя бы одного из Уточнение пороговой размерности системы осуществлено с точностью до 0.001.

Анализ устойчивости полученных ФТ позволил сделать следующие выводы:

1) для модели Изинга (п ~ 1) при размерности системы d ниже пороговой d,: -- 3.98G устойчива ФТ2. А т.к. во всем интервале изменения размерности системы 3 < d < 4 остальные фиксированные точки остаются неустойчивыми (n = 1), следовательно, при Л > 3.98G в системе за счет эффектов НРС вообще не реализуется устойчивое критическое поведение.

2) для трехмерной XY-модели получаемые малые положительные значения А, указывают на слабую устойчивость реплично-симметричной ФТ2 (п — 2). Однако уже при d, = 3.1 устойчивой становится ФТЗ с эффектами НРС. При этом критическое поведение, определяемое ФТЗ, является пеуниверсальньш и зависящим от величины параметра х0, а, следовательно, от концентрации примесей. Выявлено, что ФТЗ устойчива лишь для интервала 0 < хц < x,.(d), где л\. - пороговое значение параметра, зависящее от размерности системы. Так, для d = 3.1 х,. = 0.1, а для d — 3.999 хс = 0.3. В интервале x,.(d) < х» < 1 ни одна из ФТ не является устойчивой;

3) для изотропной трехмерной модели Гейзенберга (п = 3) устойчивой становится ФТ1, в то время как в других ФТ константы y„,gl принимают нефизические отрицательные значения. Лишь при d > 4.0 значения констант g(*,</f для ФТЗ принимают физические значения и одновременно ФТЗ становится устойчивой для интервала 0 < х0 < 0.4. В интервале 0.4 < хп < 1 ни одна из ФТ не является устойчивой.

Таким образом, проведенные исследования показали, что критическое поведение двумерных и трехмерных систем устойчиво относительно влияния эффектов НРС. В системах с п = 1 реализуется критическое поведение, определяемое структурным беспорядком с реплично-симметричной фиксированной точкой. Наличие слабого беспорядка не влияет на критическое поведение систем с п > 2.

В пятой главе осуществлено теоретико-полевое описание и численное исследование влияния эффектов корреляции дефектов на критическое поведение трехмерных систем.

В первом параграфе главы проведено теоретико-полевое описание критического поведения систем с изотропной дальнодействующей корреляцией пространственного распределения дефектов на основе модели Вейнриба-Гальперина (Вейнриб А. и Гальперин В.И., 1983)20. Модель характеризуется эффективным гамильтонианом

= + + - Е/<l"xd"y .у(х-у)^(х)^(у) (47)

где фJ = IC/fci^)2 и Ф1 (х) " ('пхп)-компонентный параметр порядка, и(х-у) - корреляционная функция потенциала поля дефектов, которая, для модели с дальнодействующей корреляцией дефектов задается в виде д(х — у) ~ |х — у|~". При этом, фурье-образ вершины у(х — у) имеет

aoWeinrib A., Halperin B.I. // PUys. Rev. В. - 1983. - V. 27. - P. 413.

n

2

III

(2,84; 3) \ 1

\O.Pj

I 11

(2.93; 1) i -— ov^lrt-S-O a

1

H—'-1-1—I-r—1—'—l' I '—I—'—I-1—I—I—I—•—1

2,0 2,2 2,4 2.6 2,8 3.0 3.2 3,4 3,6 3.8 4,0

Рис. 8. Области устойчивого критического поведения системы с дальнодействующей корреляцией дефектов определенные в настоящей работе в рамках теоретико-полевого подхода дня d = 3.

при малых волновых векторах к. вид д{к) = vo + woka~''. Из положительной определенности д{к) следует, что для а > d вклад слагаемого с iu0 Ф 0 несущественен, поэтому v<¡ > 0, и гамильтониан (47) описывает примесную модель с некоррелированными дефектами. Для а < d определяющим становится вклад второго слагаемого, пропорционального w0, поэтому п>0 > 0. Таким образом, ренормгрупповые преобразования гамильтониана (47) задаются тремя параметрами ио.^о,™»-Система с дальнодействующей корреляцией дефектов структуры является типичным примером многовершинной модели, для которой предсказания, сделанные на основе применения метода e-разложения (Вейнриб А. и Гальперин В.И., 1983), не являются надежными. Это объясняете конкуренцией различных типов фиксированных точек в многопараметрическом пространстве модели,'что делает протяжку с 1 невозможной без пересечения областей стабильности различных фиксированных точек.

Теоретико-полевое описание модели в двухпетлевом приближении и расчет соответствующих диаграмм Фейнмана, дающих вклад в двух- и четырехточечные вершинные функции, непосредственно при <1—3 и численных значениях параметра корреляции 2 < а < 3, проведенный с шагом Дн = 0,01, позволили определить ß- и 7 - функции - коэффициенты РГ уравнений, в виде [>ядов по неренормированным константам связи u, v, w. С применением обобщенного на трех-параметрический мучай метода Паде-Бореля были определены фиксированные точки модели и условия их устойчивости. Было выделено три типа фиксированных точек в представляющей физический интерес области значений параметров и, v, ta: I - соответствует критическому поведению однородной системы (и* ф 0,v*,u.<* = 0), II - системы с некоррелированными дефектами (и", у* ф 0, i/t* = 0) и III - системы с дальнодействующей корреляцией дефектов (u',v", w' ф 0). Анализ устойчивости фиксированных точек позволил определить области существования различных типов критического поведения на плоскости (а, п), где п - число компонент параметра порядка (рис.8). Картина областей устойчивого критического поведения существенно отличается от полученной ранее в (Вейнриб А. и Гальперин Б.И., 1983) в рамках с, 5 - разложения. В диссертации подсказаны возможные изменения областей критического поведения в более высоких порядках теории.

Для различных п и 2 < а < 3 были вычислены значения критических индексов (табл. 4), характеризующих как равновесное, так и неравновесное критическое поведение (Прудников П.В. и др., 2000, 2011)21. Полученные значения индексов существенно отличаются от предсказывае-

-'PiudniW V.V., Pniduibv P.v.. Feíloreuku A.A. П Pliys. Rev. В. - 200U. - V. 62. - P. 8777; Прудников П.В., Куликов Д.Н.// Письма в ЖЭТФ -'2011 - Т. 83. - Вып. 2. - С. 105

Таблица 4. Критические индексы для системы с дальнодействующей корреляцией д<ч)мктов

n = 1 71 = 2 и = 3

а ц f г в' т/ u г в' г/ V z 0'

3.0 0.021(15)0.869(02)2.168(02)0.329(10) 2.8 0.018(13)0.688(04)2.252(04)0.129(11) 2.6 0.014(10)0.701(05)2.317(06)0.131(12) 2.4 0.009(11)0.711(05)2.377(11)0.133(13) 2.2 0.001(12)0.725(06)2.436(15)0.135(14) 2.0-0.018(14)0.735(06)2.505(19)0.138(14) 0.019(13)0.655(09)2.000(01)0.172(11) 0.014(10)0.717(04)2.150(01)0.171(14) 0.000(06)0.736(06)2.202(00)0.175(10) -0.004(07)0.763(10)2.268(08)0.178(17) -0.014(11)0.775(12)2.324(12)0.182(18) -0.032(23) 0.825(37) 2.376(20)0.185( 19) 0.018(13)0.696(07)2.000(01)0.204(13) 0.013(09)0.734(11)2.109(06)0.206(17) 0.007(06)0.768(11)2.143(14)0.212(19) -0.010(08)0.789(10)2.205(24)0.217(20) -0<Ш(16)0.79«(10)2.253(29)<).222(22/ -0.020(15)0.802(14)2.3119(48)0.226(23)

мых в рамках метода e,S- разложения. Дополнительно к результатам ранней работы (Прудников П.В. и др., 2000) в работе (Прудников П.В. и др., 2010)" для пересуммирования рядов теории наряду с методом Паде-Бореля были применены методы Пад^Бореля-Лероя и автомодельного приближения (Юкалов В.И., 1997)23. В табл. 4 представлены средние значения критических индексов и их погрешности, полученные с применением данных методов.

Вычисленные в диссертации значения динамических индексов г и ff демонстрируют, что с увеличением пространственной корреляции дефектов (уменьшением параметра а) происходит значительное замедление процессов критической релаксации в системе и усиление влияния неравновесных начальных состояний по сравнению с однородными системами и системами с некоррелированными дефектами.

Во втором параграфе главы осуществлено численное исследование трехмерных систем с линейными дефектами, изотропно ориентированными в образце.

Была рассмотрена неупорядоченная модель с гамиьтонианом (30), где S, - это n-мерный единичный вектор в узле. Общая спиновая концентрация в системе была выбрана равной р-=0.80. Полагается, что дальнодействующие эффекты корреляции между точечными дефектами реализуются в виде случайно ориентированных линий с корреляционными характеристиками, спадающими по степенному закону с показателем а = 2. Для этого был использован следующий способ создания примесных конфигураций: из заполненной спинами трехмерной решетки случайным образом удалялись линии, параллельные осям координат, до достижения заданной концентрации примесей. Для обеспечения изотропности распределения дефектов в кристалле число удаляемых линий в каждом из трех направлений поддерживалось одинаковым.

Было осуществлено определение критической температуры системы методом кумулянтов Биндера 4-го порядка. Данный кумулянт задается следующим выражением

£/4(Г, ¿) = i (3 - [(MW)] /([(М<2»)])'2) , (48)

U4(L, Т) имеет важную для описания поведения конечных систем скейлинговую форму Ut(L, Т) = — Те)\, которая позволяет определить критическую температуру TC(L = ос) для бесконечной системы через координату точки пересечения кривых, задающих температурную зависимость Ui(L, Т) для различных L. Для уменьшения эффектов критического замедления использовался кластерный алгоритм моделирования Вольфа. За один шаг Монте-Карло на спин (MCS/s) принималось 5 переворотов кластера Вольфа. Процедуре установления термодинами-

^Прудников П.В., Яковлев М.И., Бакланов A.B. и др. // Весгн. Омск. Унив. - 2010. - Вып. 2. - С. 62.

MGluimaii S., Yukalov V.l. // Pliys. Rev. Lett. - 1997. - V. 79. - P. 4.

3,922 3,926 3,930 Г 1,8605 1,8620 1,8635 1,8650 1,190 1,200 1,210 Т

Рис. 9. Температурные зависимости кумулянта (Т, L) для трехмерных моделей Изинга (a), XY (6) и

Гейзенберга (в) с линейными дефектами для различных линейных размеров решетки L ческого равновесия в системе отводилось 104 MCS/s и 105 MCS/s отводилось на статистическое усреднение вычисляемых характеристик системы при заданной примесной конфигурации. Для определения средних значений термодинамических и корреляционных функций наряду со статистическим усреднением применялось усреднение по различным примесным конфигурациям: для решеток с L = 16 и 32 - по 15000 образцам, для решеток с £ = 64 и 128 - по 10000 образцам.

В результате анализа температурного поведения кумулянтов (рис. 9) и отношения Ç/L (32) для кубических решеток с размерами L от 16 до 128 были определены критические температуры: для трехмерной модели Изинга с линейными дефектами - Тс = 3.9281(1), ХУ-модели - Т, — 1.8626(5), модели Гейзенберга - Т, = 1.197(2). При данных критических температурах для каждой модельной системы было осуществлено численное моделирование релаксации из полностью упорядоченного начального состояния с m0 = 1, а также эволюции из начальных состояний с m(, « 1. Временная зависимость логарифмической производной намагниченности была получена на основе квадратичной интерполяции по трем кривым m(t) при температурах Т — 3.9250, 3.9281, 3.9310. Результирующие кривые были получены усреднением по 3000 образцам с различными конфигурациями линейных дефектов.

Анализ зависимости кумулянта U2{t) показал, что во временном интервале [50,150] MCS/s, степенному характеру зависимости U2(t) соответствует значение динамического индекса г 2.02, описывающее критическое поведение однородной модели Изинга (Ченг Б., 1998)24, а влияние линейных дефектов начинает проявляться лишь на временах t > 400 MCS/s. Выявленные динамические кроссоверные явления были учтены и при анализе временных зависимостей намагниченности и ее логарифмической производной. При определении показателей, определяющих режим критического поведения неупорядоченной системы, был проведен учет ведущих поправок к асимптотической зависимости измеряемых величин, согласно (41).

При исследовании критической эволюции систем из начальных состояний с тщ <g. 1 осуществлялось численное определение временной зависимости следующих величин: намагниченности m(t) (38), второго момента намагниченности m'^(i) и автокорреляционной функции A(t) (40). На основе показателей, характеризующих степенное поведение данных величин, были определены критические индексы (i/v и

В слабо неупорядоченных системгис с дальнодействующей корреляцией дефектов в отличие от поведения однородных систем было выявлено два универсальных динамических критических режима: на раннем временном интервале t. = [10,100) реализуется неравновесное критическое

24Zlieng В. // Int.. .Т. Mod. Phys. В. - 1Э98. - V. 12. - Р. 1410.

Таблица 5. Критические индексы для системы с дальнодействующей корреляцией дефектов

модель г в' и

Изинга та = 1 2.489(21) 0.719(22) 0.375(45) 0.8(1)

ХУ та = 1 2.364(7) 0.778(26) 0.37(30) 1.05(4)

Гейзенберга то = 1 2.245(60) 0.757(26) 0.388(15) 0.780(45)

Изинга то <К I 2.517(32) 0.149(11) 0.492(28)

ХУ то < 1 2.342(57) 0.374(11) 0.534(35)

Гейзенберга то <8 1 2.330(72) 0.365(71) 0.457(00) 0.275(14)

Изинга ( < 100 МСЭ/в то « 1 2.049(27) 0.101(10) 0.501(27)

Изинга, однородная 22 ТЩ « 1 2.041(18) 0.108(2) 0.510(14)

поведение, соответствующее однородной системе (с критическими индексами, например, для модели Изинга в' = 0.101(10), г = 2.049(27), ¡}/и = 0.501(27)), а затем, проходя на временах { = [100,200] МСЗ/б через режим кроссоверного поведения, реализуется динамический режим критического поведения неупорядоченной системы (табл. 5).

Итоговые значения критических индексов для рассмотренных трехмерных моделей Изинга, ХУ и Гейзенберга с линейными дефектами при моделировании как из полностью упорядоченного начального состояния то = 1. так и из начального состояния с та -С 1 приведены в табл. 5. Сопоставление полученных значений с результатами теоретико-полевого подхода для случая а — 2 (табл. 4) показывает их хорошее согласие, особенно для модели Изинга и ХУ-модели, в пределах статистических погрешностей моделирования и аппроксимнций.

С целью проверки применимости метода коротковременной динамики для численного описания структурно неупорядоченных систем и полученных с его помощью результатов в диссертации было реализовано также компьютерное моделирование трехмерной модели Изинга с линейными дефектами традиционными методами Монте-Карло в состоянии равновесия. Для моделирования использовался однокластерный алгоритм Вольфа. В состоянии равновесия при критической температуре Тс = 3.9281 был проведен расчет различных термодинамических и корреляционных функций, таких как намагниченность, восприимчивость, корреляционная длина, теплоемкость и кумулянт Биндера £/4, для решеток с линейными размерами Ь = 16+ 128 и спиновой концентрацией р — 0.80. Применение к рассматриваемым функциям процедуры учета конечноразмерных скейлинговых поправок С(Ь) ~ Ь"!"( 1 + аЬ~"), М(Ь) ~ Ь'!,/"( 1 + Иг"), ХЩ ~ + сЬ~и), ~ £''"(1 + ЛЬ"") позволил определить критические индексы а,

/3, 7, V и ш с помощью статистической обработки данных моделирования. Для этого проводилась линейная аппроксимация зависимостей (А'/Гд) от где X - исследуемая функция, Д - соответствующий критический показатель. Затем исследовались графики зависимостей погрешности аппроксимации а функций XЬ~Л(Ь~") при варьировании значений параметров Д и ш. По минимуму и определялись значения показателей. Минимум погрешности аппроксимации наблюдался при и = 0.90 для а/и, и> = 0.65 для р/и, ш = 0.70 для 7/у, и = 0.80 для и. Выло взято среднее значение ш = 0.76(5), при котором и были рассчитаны окончательные значения показателей: а - -0.078(30), = 0.362(20), 7 = 1.441(15), и = 0.710(10). Полученные значения находятся, в хорошем согласии в пределах погрешностей измерений с результатами применения метода коротковременной динамики и теоретико-полевого описания.

В третьем параграфе главы представлены результаты компьютерного моделирования трехмерных моделей Изинга и ХУ с линейными дефектами в области их сильной неупорядоченности с р = 0.60. Ставилось целью проверить, проявляется ли в критическом поведении систем с дальней пространственной корреляцией дефектов зависимость от их концентрации.

Для получения значений критических температур для сильно неупорядоченных моделей Изинга и ХУ было осуществлено вычисление равновесных значений кумулянта Биндера и^ (48). Из пересечения температурных зависимостей иА(Т) для решеток с размерами Ь = 64,96,128 были получены значения Тс = 3.1956(34) для модели Изинга, и Тс = 1.4455(5) для ХУ-модели. В этом случае для моделирования применялся однокластерный алгоритм Вольфа (25 переворотов кластера за МСЭ/.ч). Было использовано 5000 МСБ/з для термализации и 50000 МСЭ/а для статистического усреднения по спиновым конфигурациям. Окончательные результаты были получены усреднением по 500 образцам с различными конфигурациями линейных дефектов.

Далее методом коротковременной динамики была исследована релаксация данных моделей из начального состояния с т<> = 1. Анализировались временные зависимости намагниченности т(<) (43), ее. логарифмической производной по температуре Э1п тп(Ь) (44) и кумулянта {/2(г) (45). Результирующие кривые были получены усреднением по 80 образцам для модели Изинга и по 200 образцам для ХУ-модели с различными конфигурациями распределения линейных дефектов в решетке и по 25 прогонкам для каждого образца. Зависимость логарифмической производной намагниченности при критической температуре была получена на основе разностной схемы по двум кривым >п(4) при температурах Т = 3.1956 и Т — 3.2180.

Анализ данных зависимостей с применением процедуры учета поправок к скейлингу позволил рассчитать значения критических показателей для модели Изинга: /5/у = 0.396(52), г = 2.707(34), ¡) = 0.420(97), и = 1.061(155) и для ХУ-модели: ¡З/и = 0.51(1), г = 2.589(52), /9 = 0.511(72), »/= 1.001(143).

При исследовании эволюции системы из начальных состояний осуществлялся ана-

лиз временных зависимостей следующих величин: намагниченности т(() (38), второго момента намагниченности ш<2>(<) и автокорреляционной функции АЦ) (40). При расчете критического индекса в' было реализовано компьютерное моделирование при значениях начальной намагниченности щ, - 0.002 и т-ц = 0.01 для модели Изинга и при значениях т0 = 0.004, т0 — 0.003 и = 0.002 для ХУ-модели с последующей линейной аппроксимацией результатов к то = 0. Данные усреднялись по 30 для модели Изинга и по 1300 для ХУ-модели различным примесным конфигурациям при 25 "прогонках" для каждой примесной конфигурации. Анализ временных зависимостей позволил определить значения показателей в -- 0.210(19), /З/к = 0.409(52), г = 2.699(90) для модели Изинга и в' = 0.393(45), = 0.564(48), г = 2.511(95) для ХУ-модели.

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

Таблица 6. Значения критических показателей для систем с линейными дефектами при р - 0.00

/3 V 0

модель Изинга ХУ-модель 2.703(48) 2.550(54) 0.403(38) 0.537(25) 0.420(97) 0.511(72) 1.001(155) 1.001(143) 0.210(19) 0.393(45)

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

Шестая глава диссертации посвящена исследованию влияния дефектов структуры и :*)> фектов их пространственной корреляции на аномальные особенности распространения ультразвука в твердых телах при фазовых переходах второго рода.

В первом параграфе главы представлено теоретико-полевое описание процесса распространения ультразвука в твердых телах вблизи температуры фазового перехода второго рода для однородных и неупорядоченных систем с некоррелированными де^ктами структуры. Рассматривается гамильтониан трехмерной неупорядоченной сжимаемой модели Изинга в

ЕИДе Я = Яы + Яор + Яы + Д,„р. (49)

Вклад деформационных степеней свободы определяется следующим образом

= (с{\ £г£а + 2С»2£■«„„«,, + 4С"4 £, (50)

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

н^ = ] а"* [¿т„5г + \ (УЗ)2 + , (51)

где 5(х) - спиновый параметр порядка; щ - положительная константа взаимодействия, Та = (Т — Тос)/Тос- Составляющая задает спин-фононное взаимодействие

= IА |й»£«а

где да параметр квадратичной стрикции. описывает влияние дефектов структуры на критическое поведение системы

Я„лр = У Л*х [Дт(х).?2] + I <1'>х

где посредством случайных и гауссовски распределенных переменных Дт(.т) и Ь,(х) задаются флуктуации лохальной температуры фазового перехода и случайные поля напряжений. Функция Дт(х) характеризует потенциал случайного поля дефектов в точке со средним значением по распределению дефектов равным нулю ((Дт(.г))) = 0 и вторым моментом

((Дт(.т) Дг(у))) = у(х — у), (54)

характеризующим корреляционные свойства пространственного распределения дефектов. В случае рассмотрения однородных систем у(х - у) = 0, а в случае некоррелированных точечных

(52)

(53)

дефектов <j(x — у) = vS'l{x — у), где v - вершина, характеризующая взаимодействие флуктуаций параметра порядка через поле дефектов.

После перехода к фурье-образам переменных и проведения интегрирования по несущественным для критического поведения системы в упруго-изотропной среде недиагональным компонентам тензора деформации искомый эффективный гамильтониан принимает вид:

Й = \ J^'l (Ъ + Ч2) d'V Ar^S,,^ J ddqqhvQq,x + aQ J ddq<fQq,xQ-q,x-

- К J d'1, (S, (S, 5_„) - <j0 j <1"V q S„ 5,_„ + Jd6<i S„ S„a . (55)

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

Релаксационная критическая динамика сжимаемых систем описывается динамическими уравнениями типа обобщенных уравнений Ланжевена

= + = + + (56)

где Го и Д, - затравочные кинетические коэффициенты; i) и rj^(x,t) - гауссовски распределенные величины, имеющие характер случайной силы; ft® и - поля, термодинамически сопряженные соответственно спиновым и деформационным переменным.

В результате итерационного решения системы нелинейных уравнений (56) с гамильтонианом (55) могут быть выделены функция отклика для упругих переменных D{q, и) и функция отклика для спиновых переменных G(q,u)

Щд, ш) = = ((<?„.,xQ-„_„,,)], G(4, и) = М = KS,uS-„-„>], (57)

где (...) обозначает статистическое усреднение по случайным ланжевеновским силам, [...] - усреднение по флуктуациям случайных полей Дт_,; и /i,„ определяемых дефектами структуры. Диа-. грамыное представление собственно-энергетической части £(<?, ш) функции отклика D(q, w) для системы с некоррелированными структурным беспорядком в двухпетлевом приближении имеет следующий вид

ЕМ = 4'Г - 96<?ч

+ 16<fv + 16 ,fv (^У + 16<j2v

^ lGfv ^^ г 1Су21) 4 1С,/» ^^

где сплошной линии соответствует Git(q, ui) = 1 ¡{iui /Г0 + (т0 + q2)), линии с крестом - Co(q, ш) = 2Г,7' [(ai/Го)2 + ((/2 + Гц)2]-1, вершине с прерывистой линией - v = [(Дг)2]. Слагаемые пропорциональные ЛЯ описывают релаксационный вклад, обусловленный тем, что при Т <ТС спиновая плотность = MSqji + ip,t наряду с флуктуационной частью ¡рч содержит намагниченность системы М — В\Т — Tcf, где В - феноменологический релаксационный параметр.

1ц у Ц SÍ

Рис. 10. Скейлинговые функции Ф(у) (а) и /(у) (6) для однородной (1) и неупорядоченной (2) систем при Т > Тс, (Г) и (2') соответственно при Т < Тс.

Таблица 7. Асимптотическое поведение коэффициента поглощения и дисперсии скорости звука Коэффициент поглощения сг(ю, г)

Режим Однородная Неупорядоченная

Т < Г,. Т> Г, Т<1). Т>Т,

Критический у = 101 -г 10J Предкритический у = 1 -г 102 Гидродинамический у = 10~3 -г Ю-1 , .0.98——0.08 W т Ш2Г-1.38 ш1.(15г-0.17 W2T-1.38 ш1.22т-0.25 Ш1.37Т-0М8

Дисперсия скорости звука <?(ы, т) — c'fO, т)

Режим Однородная Неупорядоченная

Т<Т,. Т>ТГ Т<ТГ т>тг

Критический у = 1Ü1 -г 10-1 Предкритический у = 1 -г 103 Гидродинамический у — 10~3 -г 10_| U(U1T-(I.2 5 LJ i ш2 т-гж, ыомт-нм W(Í.«JT-I>S4 ш].«1г-1.45

Собственноэнергетическая часть непосредственно связана с характеристиками рас-

пространения ультразвука: коэффициентом поглощения и дисперсией скорости звука

а(и,т) ~ш1т£(и>), 1?(и>,т)-с*(1Кт)~ЯеС£(и>)-£(0)). (58)

Из теории рассеяния ультразвука в твердых телах вблизи температуры фазового перехода известно, что в асимптотическом пределе (т -> 0, и 0) могут быть введены следующие выражения для мнимой и действительной составляющих £(ы, т)

Ие(Е(0, ш) — £(0,0)) = т~" (/(;(/) — /(О)), (59)

и*

со скейлинговыми функциями ф(у) и ](у), зависящими от единственной обобщенной переменной у = шт~г" /Го - В диссертации непосредственно для трехмерных систем в двухпетлевом приближении был осуществлен расчет динамических скейлинговых функций с применением метода суммирования Паде-Бореля. Их поведение представлено на рис. 10.

Результаты проведенных расчетов асимптотических зависимостей коэффициент« поглощения и дисперсии скорости звука для критической (у ~ 3> 1) и гидродинамической (у ~ Ыр ~ (</£)* -С 1) областей представлены в табл. 7. Показатели их частотной и температурной зависимостей для гидродинамического режима определялись для интервала 10~3 < у < Ю-1, а критического режима для 10 < у < 103. Реальной температурной области с 10_:! <т< Ю-1 для ультразвуковых исследований фазовых переходов соответствует интервал 1 < у < 102, т.е. он захватывает кроссоверную область и начало критической области (предкритический режим).

Рис. 11. Температурное поведение коэффициента поглощения а(и;,т) для систем: однородной (а1), с некоррелированными дефектами (а2), с дальнодействукмцей корреляцией дефектов (6) Из представленных в табл. 7 результатов следует, что наличие дефектов структуры приводит к более сильному, чем для однородной системы, увеличению коэффициента поглощения и дисперсии скорости звука по мере приближения к критической температуре уже в гидродинамической области. В критической области для структурно неупорядоченных систем должна наблюдаться как более сильная частотная, так и температурная зависимость акустических характеристик по сравнению с их однородными аналогами. Особенно важным представляется предсказываемое проявление динамических эффектов влияния дефектов структуры в аномальном поглощении и дисперсии скорости звука в более широком температурном интервале относительно критической температуры, чем в других экспериментальных методах, в которых для выявления данных эффектов необходимо проводить исследования в узком температурном интервале вплоть до г ~ Ю-4. Данные выводы находят подтверждение в модельном представлении результатов расчета температурного критического поведения коэффициента поглощения и дисперсии скорости звука (рис. 11) для однородной и неупорядоченной систем, проведенных при значениях параметров В = 0.3 и ш/Ги = 0.0015.

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

Выла применена модель изотропной дальнодействующей корреляции дефектов, задаваемая в (54) корреляционной функцией у(х — у) ~ \x-y\~" с а - показателем корреляции. При а < 3 эффекты корреляции являются существенными и определяют новый тип критического поведения. При а >3 реализуется критическое поведение, характерное для некоррелированных дефектов. Случай с. а = 2 соответствует поведению систем со случайно ориентированными линейными дефектами, с 2 < а < 3 - системам с примесными фрактальными структурами.

Диаграммное представление дополнительных по сравнению с рассмотренным выше случаем системы с некоррелированными дефектами вкладов в собственно-энергетическую часть £(?,ш) динамической функции отклика £>(г/,о>), возникающих в двухпетлевом приближении, может быть представлено в виде:

11>ц2т + 1С'/«' I ^^ I 1йд2т

где волнистая линия в диаграммах характеризует влияние эффектов корреляции дефектов через вершину взаимодействия те с передачей импульса qл~,i.

Ш*,)

«лад

04

-си

4.4 а . Т>Те

4.1

■1.6

•13 Т< Тс > \\Х «-2.4

Л Л а-3.7

•и о-З.О.

•и

40 \ и-2.7-

•гл 4 •■з.о-

...

о-2.<

: б

а-2.4 /

4*2.7 ¡Г Ру ^ а'3.0

а'3.0 [/

Т < Тс / / Т>Тс

|ду

|ду

Рис. 12. Скейлинговые функции ф(у) (а) и ¡(у) (б) дли различных значений а

В диссертации был осуществлен расчет в двухпетлевом приближении динамических скейлин-говых функций для различных значений параметра корреляции 2 < а < 3. Поведение скейлин-говых функций ф(у) и ¡(у) для системы с дальнодейсгвующей корреляцией дефектов, рассчитанное с применением метода суммирования Паде-Бореля, представлено на рис. 12 в двойном логарифмическом масштабе. Для различных интервалов изменения переменной у в поведении Ф{у) и ¡(у) были выделены асимптотические области гидродинамического, предкритического и критического поведения системы, а затем проведен расчет коэффициента поглощения а(и,т) и дисперсии скорости звука г?(ы, т) — <?(0, т) для соответствующих областей и выделены значения показателей их частотной и температурной зависимостей.

Анализ значений данных показателей позволил сделать следующие выводы: для систем с дальнодействующей корреляцией дефектов предсказывается более сильное, чем для систем с некоррелированными дефектами, увеличение коэффициента поглощения по мере усиления эффектов корреляции (характеризуется уменьшением параметра а в интервале 2 < а < 3) при приближении к температуре фазового перехода уже в гидродинамической области; в критической области эффекты корреляции дефектов проявляются в виде более сильной как частотной, так и температурной зависимостей коэффициента поглощения по сравнению со случаем некоррелированного структурного беспорядка; аналогичные выводы справедливы для частотной и температурной зависимостей дисперсии скорости звука в критической области.

Сказанное находит подтверждение в модельном представлении 1>езультатов расчета температурного критического поведения коэффициента поглощения для неупорядоченных систем с различными значениями параметра корреляции, проведенных при значениях параметров В = 0.3 и ш/Г0 = 0.0015 (рис.116).

Седьмая глава посвящена изложению методики и результатов описания мультикритичес.ко-го поведения однородных и неупорядоченных систем с двумя параметрами порядка и выявлению влияния поверхности на мультикритическое поведение в окрестности точки Лифшица.

В первом параграфе проведено в двухпетлевом приближении описание критического поведения непосредственно трехмерных систем с двумя взаимодействующими параметрами порядка, которое выявило существенное изменение областей с различным типом устойчивого мультикри-тического поведения на плоскости (п — тге) - числа компонент для параметров порядка, по сравнению с полученными ранее результатами (рис. 13) при использовании с-разложения. Показано, что устойчивое мультикритическое поведение, соответствующее изотропной ФТ1 с флуктуаци-

онно индуцированной асимптотической симметрией системы 50(п + т), возможно только для п + т < 2.91, т.е. наивысшей асимптотической симметрией системы является 50(2). Значительное изменение претерпели области стабильности и остальных двух типов фиксированных точек: ФТ2 (с симметрией 30(п) х 50(т)) и ФТЗ (с симметрией 50(п) © 50(т)), соответствующих устойчивому тетракритическому поведению. Показано, что изменение областей стабильности фиксированных точек приводит к заметному изменению типов фазовых диаграмм систем во флуктуационной области. Устойчивое бикритическое поведение предсказывается в диссертации только для взаимодействующих однокомпонентных параметров порядка (п,т = 1) с критическими индексами, соответствующими ХУ-модели. Тетракритическое же поведение должно иметь более широкую реализацию среди систем с многокомпонентными параметрами порядка. Обсуждены эффекты флуктуационной неустойчивости мультикритического поведения.

Во втором параграфе проведено исследование влияния точечных замороженных примесей, создающих в системах с двумя параметрами порядка флуктуации случайной локальной температуры, на характер фиксированных точек и их устойчивость. Показано, что присутствие примесей приводит к флуктуационному расцеплению связи параметров порядка и осуществлению единственного типа устойчивого мультикритического поведения - тетракритического с общей симметрией системы 80(п) © 80(т). В случае однокомпонентных параметров порядка (п = т = 1) наличие примесей существенно и приводит к критическому поведению с индексами, соответствующими индексам неупорядоченной модели Изинга. Для систем с многокомпонентными параметрами порядка присутствие примесей не сказывается и мультикритическое поведение носит тетракритический характер однородной системы. Наличие примесей приводит к существенному сокращению по сравнению с однородными системами возможных типов фазовых диаграмм. Принципиальный момент изменения связан с невозможностью реализации в неупорядоченных системах фазовой диаграммы, содержащей бикритическую точку.

В третьем параграфе главы проведено теоретико-полевое описание мультикритического поведения (¿-мерной системы с многокомпонентным параметром порядка с учетом влияния поверхности системы, перпендикулярной оси анизотропии кристалла, в окрестности обобщенной точки Лифшица 171-го порядка.

Хорошо известно, что в критической точке параметр порядка характеризуется изотротизй масштабной инвариантностью ф{х) = 1А* ф(Ь~1х) с = {<1-2 + г?)/2. В случае критического поведения вблизи точки Лифшица т-го порядка проявляется анизотропная масштабная

инвариантность, и координатные оси X могут быть разделены на две группы х = (а:„.а;()) с различным пространственным масштабом: ф(хи,хр) = 1А* ф(Ь~^х0,Ь~1хц), где <р - индекс анизотропии, индексом а отмечены оси модуляции, /У - смодулированные оси. с» принимает значения от 1 до т., а ¡3 от т + 1 до

Критическое поведение системы чувствительно к влиянию границы образца, т.к. дапьнодей-ствующий характер взаимодействия флуктуации параметра порядка модифицируется вблизи поверхности. В окрестности точки Лифшица возможно два сценария влияния поверхности на критическое поведение: параллельный случай, когда все оси Модуляции и параллельны поверхности (Диль Г.В. и др., 20 03)25 и, более сложный случай, когда оси модуляции перпендикулярны поверхности. Для описания критического поведения полубесконечной системы был применен гамильтониан Н = / Съ(х) дГх^хц + / С^х) , (СО) с объемной 1

У„ф)1 + и> + 1\ф\\ (01)

Ц=т+1

и поверхностной плотностью

^ = ^ 4? + ^ Х><М2 + ЪФЬФ + ^ (д„ф)\ (62)

Из действия определяемого гамильтонианом (СО) могут быть получены граничные условия

тп-1 ...-1

+ Ф - р]дп + с± - [Л|, + (/-а)д„] £ а2\ф = о, { - <7 + лха„+ь - (/+¿о £ о2 }</> = о. (ез)

Для описания критического поведения локальной намагниченности в точке Лифшица вблизи поверхности (г = 0) может быть введен критический индекс

т^ = 0) ~, г-И), (04)

причем значения данного индекса существенно отличаются от значения индекса для бесконечной системы. Для устранения ультрафиолетовых расходимостей в корреляционных функциях была применена процедура размерной регуляризации с последующим переопределением параметров гамильтониана (60) и мультипликативной перенормировкой:

ф = г1'2ф"", а = 2„о, иа~'"1л = Е„^цс2пи, (65)

Т - тьр = [г + Ат ,,2], (¡, - р,.г)а-1/2 = ,,2„,,, (66)

где е = <Г - /1, <Г = 4 + т/2 - верхняя критическая размерность системы, ^ = [Г(1 + с/2) Г2(1 - с/2) Г(т/4)]/[(47г)'8+,"~2е'/'1 Г(2 - е) Г(т/2)].

За счет введения в теорию двух граничных условий вида (63), возникает необходимость в перенормировке величины д2ф, задающей поверхностное поведение параметра порядка ф\

д2пф = [2г 1/2 (а2 ф)"" = г\'2с. (67)

Рассчитанное в однопетлевом приближении выражение для константы перенормировки 2-1 име-

25КеЫ НДУ., СагошИ А., ШикеухсЬ Б. // РЬув. В. - 2(103. - V. 08. - Р. 224428.

Взятие производной -«<?„ от Zi позволяет определить РГ функцию ife. Подстановка значения и в фиксированной точке и* — 6е/(п + 8) + 0(с2) дает выражение

4Î = SbH ~ и' + О [(и-)2] = да(ш) п±| £ + 0(£Г2) . (69)

Используя известное ¿--представление для индекса и может быть получено итоговое выражение для независимого поверхностного критического индекса ¡З^ в аиде ряда по е:

АШ = ("/2)[4-е + »й] + 0(г2) = 1 + -

« - 4 + —— gm(m)

4 (п + 8)

Для случая одноосной модели Изинга m = n = 1 выражение упрощается

11

1 + ^£ + 0(е2)

+ 0(е2). (70)

(71)

Просуммированное с Паде-апроксимантом |0/1) в точке е = 3/2 для с2 = 3 значение индекса = 0.59 хорошо согласуется с результатами компьютерного моделирования 0.62(1) (Плейм-линг М, 2002)21'.

В заключении сформулированы основные результаты и выводы диссертации.

Основные результаты и выводы

1. В рамках единого подхода, основанного на формализме динамического производящего функционала для корреляционных функций и функций отклика, осуществлено теоретико-полевое описание неравновесного критического поведения с учетом влияния начальных состояний системы, нарушающих трансляционную инвариантность во времени, и квазиравновесной критической динамики, когда влияние неравновесных начальных состояний самоусреднилось, для однородных и неупорядоченных систем с некоррелированными дефектами и дефектами с даль-нодейетвующей пространственной корреляцией.

2. Впервые показано, что в неравновесной критической релаксации однородных систем только начиная с трехпетлевого приближения в динамической функции отклика возникают флукту-ационные поправки за счет влияния начальных неравновесных состояний. С учетом данных флуктуационных эффектов проведен расчет независимого динамического критического индекса в', задающей) эволюцию 71-компонентного параметра порядка системы в коротковременном |>ежиме, и при применении метода ^-разложения были получены значения (У, согласующиеся с результатами компьютерного моделирования.

3. Впервые реализовано тео^тико-полевое описание неравновесного критического поведения структурно неупорядоченных систем с некоррелированными дефектами с учетом влияния начальных состояний при фиксированной размерности системы <1 = 3 и в двухпетлевом приближении проведен расчет динамичес.ких'критических индексов & и г с последовательным применением к рядам теории различных методов суммирования. Численные значения критических индексов находятся в лучшем соответствии с результатами компьютерного моделирования, чем результаты применения метода г-разложения.

4. Впервые осуществлено последовательное применение методов суммирования Паде-Бореля, Паде-Бореля-Лероя и конформного Паде-Бореля к рядам теории, полученным в многопетлевом

'•"'Plrimlins M. // РЬу.ч. Rev. В. - 2002. - V. 65. - P. 184406.

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

5. Разработана методика теоретико-полевого описания критического поведения структурно неупорядоченных спиновых систем с эффектами нарушения репличной симметрии (НРС) при их фиксированной размерности. На примере данных систем, описываемых многовершинными моделями, показано, что предсказания, основанные на применении метода £- разложения, не являются надежными. Конкуренция различных типов фиксированных точек в многопараметрическом пространстве вершин взаимодействия не позволяет осуществлять протяжку с -+ 1,2 без пересечения областей стабильности различных фиксированных точек.

6. Для двумерных, трехмерных систем и систем произвольной размерности от трех до четырех в двухпетлевом приближении проведен ренормгрупповой анализ эффективного репличного гамильтониана модели с потенциалом взаимодействия, не являющимся реплично-симмктричным и определяемым тремя вершинами взаимодействия флуктуации. Для случая однсктупенчатсь го НРС с применением техники суммирования Паде-Бореля были выделены возможные типы критического поведения и осуществлен анализ возможности их реализации. Было установлено, что критическое поведение двумерных и трехмерных систем устойчиво относительно влияния эффектов НРС. Показано, что наличие слабого беспорядка не влияет на критическое поведение систем с многокомпонентным параметром порядка, а в системах с однокомпонентным параметром порядка реализуется критическое поведение, определяемое структурным б«,иорядком с реплично-симметричной фиксированной точкой.

7. Выявлено: эффекты НРС проявляются лишь при размерностях неупорядоченной системы больших трех, при этом пороговые значения <1,- зависят от числа компонент параметра порядка п и величины параметра х0, определяющею координату ступеньки для потенциала взаимодействия с НРС. Определены значения пороговых размерностей с', (п), отделяющих область критического поведения с эффектами НРС с/,..(п) < 4 < 4 от области, в которой данные эффекты несущественны. Эти пороговые разх№рности задают одновременно и нижнюю границу области применимости результатов ^-разложения к описанию модели слабо неупорядоченных систем с эффектами НРС. Показано, что возможные типы устойчивого критического поведения и критические индексы, полученные для неупорядоченных систем в рамках стандартного метода реплик для размерностей систем ниже пороговых г1,.(п), являются достоверными и реализуется прежний сценарий влияния задтрожеяного беспорядка на критическое поведение.

8. Осуществлено теоретико-полевое описание статических свойств критического поведения и неравновесной критической динамики неупорядоченных трехмерных систем с: дальнодействую-щей корреляцией дефектов в двухпетлевом приближении с использованием методов пересуммирования рядов теории. Для различных значений числа компонент параметра порядка и показателя корреляции а определены типы устойчивого критического поведения и значения критических индексов. Полученная картина областей устойчивого критического поведения и значения критических индексов существенно отличаются от предсказанных ранее в рамках двухпара-

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

9. Осуществлены численные Монте-Карло исследования равновесного критического поведения трехмерной спиновой модели Изинга в широком интервале изменения концентрации немагнитных атомов примеси. Для образцов с различными спиновыми концентрациями и различными линейными размерами была определена температурная зависимость корреляционной длины и магнитной восприимчивости. С помощью метода конечноразмерного скейлинга для данных величин были определены скейлинговые функции, демонстрирующие универсальное поведение в критической области, а также вычислены критические температуры и статические критические индексы с использованием ведущих поправок к асимптотической зависимости данных величин. На основе выявленных концентрационных изменений скейлинговых функций и значений критических индексов сделан вывод о существовании двух универсальных классов критического поведения разбавленной модели Изинга с различными характеристиками для слабо и сильно неупорядоченных систем.

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

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

11. На основе проведенных численных исследований методом коротковременной динамики критического поведения трехмерной модели Изинга со спиновыми концентрациями р = 0.95, 0.80, О.вО и 0.50 можно сделать следующие выводы: поведение слабо и сильно неупорядоченных систем принадлежит к различным универсальным классам критического поведения с несовпадающими в пределах статистических погрешностей значениями динамических критических индексов V и ;; в неравновесном поведении слабо неупорядоченных систем с р = 0.95 и 0.80 выявлено два режима со степенным временным изменением измеряемых величин: на раннем вре-

менном интервале реализуется критическое поведение с характеристикам» однородной системы, и затем, проходя через интервал кроссоверного поведения реализуется режим поведения неупорядоченной системы; в неравновесном поведении сильно неупорядоченных систем с р — 0.60 и 0.50 не выявлено режима с характеристиками однородной системы: для слабо неупорядоченных систем значения критических индексов, полученные с учетом ведущих поправок к'скейлингу, находятся в хорошем согласии с результатами теоретико-нолевого описания, результатами моделирования другими методами, а также согласуются с результатами экспериментальных исследований слабо неупорядоченных изинговских магнетиков; для сильно неупорядоченных систем значения динамического критического индекса полученные при анализе эволюции из разных начальных неравновесных состояний, находятся в хорошем согласии друг с другом, согласуются со значениями, определенными другими численными методами, а значение индекса ко)хл ковре-менной эволюции намагниченности в' = 0.186(39) получено впервые.

12. С целью исследования влияния дальнодействующей корреляции дефектов структуры на критическое поведение впервые осуществлено численное исследование методом коротковремен-ной динамики поведение трехмерных модельных спиновых систем (Изинга, ХУ и Гейзенберга) с линейными дефектами при общей спиновой концентрации р = 0.80. В неравновесном критическом поведении данных систем выявлено два динамических режима со степенным временным изменением измеряемых величин: на раннем временном интервале реализуется поведение с характеристиками однородной системы, а после интервала кроссоверного поведения реализуется динамический режим неупорядоченной системы. Определены значения совокупности динамических и статических критических индексов при применении методики учета ведущих поправок к скейлингу, соответствующих режиму неупорядоченной системы. Полученные значения критических индексов находятся в хорошем согласии с результатами теоретико-полевого описания на основе модели Вейнриба-Гальперина с показателем корреляции и, = 2 и результатами проведенного моделирования равновесного критического поведения. На примере исследования ХУ-модели и модели Гейзенберга впервые получено численное подтверждение существенности влияния дальней пространственной корреляции дефектов на критическое поведение систем с многокомпонентным параметром порядка.

13. Впервые осуществлено численное исследование неравновесного критического поведения трехмерных сильно неупорядоченных при р = 0.60 спиновых моделей Изинга и ХУ с даль-нодействующей корреляцией дефектов, моделируемой изотропно распределенными линейными дефектами. Сопоставление полученных значений критических индексов со значениями для слабо неупорядоченных с р = 0.80 моделей Изинга и ХУ позволяют сделать вывод о существовании различных универсальных классов критического поведения для рассматриваемых систем, соответствующих областям слабой и сильной структурной неупорядоченности.

14. Осуществлено теоретико-полевое описание аномальных свойств распространения ультразвука в однородных и структурно неупорядоченных твердых телах вблизи температуры фазового перехода второго рода с учетом как флуктуационного, так и релаксационного механизмов рассеяния. Разработана методика и проведен расчет скейлинговых функций для коэ<М>ицпен-тов поглощения и дисперсии скорости звука, а также самих коэффициентов при фиксированной размерности системы в. = 3 в двухпетлевом приближении с применением методов суммирова-

ния асимптотических рядов. Представлены численные значения показателей, характеризующих частотную и температурную зависимости скейлинговых функций, коэффициентов поглощения и дисперсии ультразвука в гидродинамической (Ю-3 <у = ит~с"/Т0 < Ю-1), нредкритической (1 < .'/ 5 Ю2) и критической (10 < у < Ю3) областях.

Рассчитанное для однородных систем поведение коэффициента поглощения а(ш,т) демонстрирует аномально сильное поглощение ультразвука в критической области и хорошо согласуется с результатами экспериментальных исследований в образцах ^е/^.

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

15. Впервые исследовано влияние эффектов дальнодействующей корреляции дефектов на аномальное рассеяние ультразвука в твердых телах при фазовых переходах второго рода. Показано, что усиление корреляции дефектов приводит к увеличению поглощения ультразвука и усилению дисперсии скорости звука в критической области. Предсказывается более сильное, чем для однородных систем и систем с некоррелированными дефектами, температурное изменение характеристик ультразвука при приближении к критической температуре уже в гидродинамической области. В результате, экспериментальное исследование критической динамики ультразвуковыми методами позволяет выявить влияние дефектов и эффектов их корреляции в более широком температурном интервале относительно критической температуры (10~3 < т < Ю-1), чем в других экспериментальных методах, в которых для выявления данных эффектов необходимо проводить исследования в узком температурном интервале вплоть до г ~ 10~4.

1С. Теоретико-полевое описание мультикритического поведения однородных систем с двумя параметрами порядка в рамках метода с фиксированной 4 = 3 размерностью выявило в двух-нетлевом приближении существенное изменение областей различного типа устойчивого мультикритического поведения на плоскости (и — т) - числа компонент данных параметров порядка по сравнению с результатами применения е - разложения. Это приводит к изменению характеристик мультикритического поведения и возможных типов фазовых диаграмм системы во флуктуационной области.

Показано, что присутствие дефектов структуры в системе приводит к флуктуационному расцеплению связи параметров порядка и осуществлению единственного типа устойчивого мультикритического поведения - тетракритического. В случае однокомпонентных параметров порядка (п,т = 1) наличие дефектов существенно и приводит к мультикритическому поведению с индексами, соответствующими индексам неупорядоченной модели Изинга. Для систем с п,тп > 2 присутствие дефектов не сказывается на мультикритическом поведении и оно носит тетракри-тический характер однородной системы. Влияние дефектов проявляется в сокращении по сравнению с: однородными системами возможных типов фазовых диаграмм.

17. Впервые исследовано влияние поверхности, перпендикулярной оси нннзотрогши кристалла, на мультикритическое поведение (i-мерной системы вблизи точки Лифшица m-го порядка.

Полученное значение нового независимого поверхностного критического индекса ¡I1 = 0.50 подтверждает результаты численного исследования модели Изингл для случая точки Лифшица с

т, 71 = 1.

Основные результаты диссертации опубликованы в работах

1. Prudnikov P.V., Pruiteikov V.V., Fedorenko А.А. Static and dynamic critical properties of 3D systems with long-range correlated quenched defects. // .1. Pliys. A: Mat.li.Gcn. - 1999. - V. 32. - N. 49. -P. 8587-8600.

2. Прудников П.В., Прудников В.В., Федоренко А.А. Мультикритическое поведение слабо неупорядоченных систем с двумя параметрами порядка. // ЖЭТФ. - 1999. - Т. 11G. - Вып. 2. - С. 611-019.

3. Prudnikov P.V., Prudnikov V.V., Fedorenko А.А. Field-theory approach to critical behavior of systems with long-range correlated defects. // Phys. Rev. B. - 2000.'- V. 62. - N. 13. - P. 8777-8780.

4. Прудников П.В., Прудников B.B., Федоренко А.А. Мультикритическое поведение неупорядоченных систем с двумя параметрами порядка. //' Физика твердого тела. - 2000. - Т. 42. -■ Вып. 1. -С. 158-162.

5. Prudnikov P.v., Prudnikov V.V., Fedorenko А.А. Stability of critical behaviour of weakly disordered systems to introduction of potentials with replica symmetry breaking. // Л. Phys. A: Math.Gen. ~ 2001.

- V. 34. - N. 12. - P. L145-L152.

0. Prudnikov P.V., Prudnikov V.V., Fedorenko A.A. Stability of critical lx:haviour of weakly disordered systems with respect to the replica symmetry breaking. // Phys. Rev. B. - 2001. - V. 03. -- N. 18. -P. 184201-184206.

7. Прудников П.В., Прудников В.В., Федоренко А.А. Устойчивость критического поведения слабо неупорядоченных систем к нарушению репличной симметрии. // Письма в ЖЭТФ. - 2001. - Т. 73.

- Вып. 3. - С. 153-158.

8. Прудников П.В., Прудников В.В., Федоренко А.А. Устойчивость критического поведения слабо неупорядоченных систем к введению потенциала взаимодействия с нарушенной репличной симметрией. // Физика твердого тела. - 2001. - Т. 43. - Вып. 9. - С. 1088-1G92.

9. Prudnikov P.V., Prudnikov V.V. Critical behaviour of weakly disordered systems with replica symmetry breaking potentials. // J. Phys. Stud. - 2001. - V. 5. - N. ¡3/4. - P. 280-292.

10. Прудников П.В., Прудников В.В. Критическое поведение неупорядоченных систем с эффектами нарушения репличной симметрии. // ЖЭТФ. - 2002. - Т. 122. - Вьш. 3. - С. G30-040.

И. Prudnikov P.V., Prudnikov V.V., Dorofeev S.V., Kolesnikov V.Yu. Monte Carlo studies of critical behaviour of systems with long-range correlated disorder. // Condensed Matter Physic s. - 2005. -V. 8. - N. 1. - P. 213-224.

12. Prudnikov P.V., Prudnikov V.V. The influence of disorder on the critical sound attenuation in solids. // J. Phys.: Condens. Matter. - 2005. - V. 17. - P. L485-L492.

13. Прудников П.В., Прудников В.В., Криницын А .С. Расчет динамического критического индекса методом сз'ммирования асимптотических рядов. // Теоретическая и математическая физика. -2006. - Т. 147. - Вып. 1. - С. 138-155.

14. Diehl H.W., Shpot М.А., Prudnikov P.V. Boundary critical behaviour at m-axial Lifshitz points of semiinfinite systems with a surface plane perpendicular to a modulation axis. // .1. Phys. A: Math.Gen. -2006. - V. 39. - P. 7927-7942.

15. Prudnikov P.V., Prudnikov V.V. Critical sound attenuation of three-dimensional Ising systems. // Condensed Matter Physics. - 200S. - V. 9. - N. 2. - P. 403-410.

16. Prudnikov P.v., Prudnikov V.V., Zheng В., Dorofeev S.V., Kolesnikov V.Yu. Short-time critical dynamics of the three-dimensional systems with long-range correlated disorder. // Progress of Theoretical Physics.

- 2007. - V. 117. - N. 0. - P. 973-991.

17. Прудников B.B., Прудников П.В., Вакилов A.H., Криницын А.С. Компьютерное моделирование критического поведения трехмерной неупорядоченной модели Изинга. // ЖЭТФ. - 2007. -- Т. 132.

- Вып. 2. - С. 417-425.

18. Прудников П.В., Прудников В.В., Носихин Е.А. Аномальное поглощение ультразвука в твердых телах при фазовых переходах второго рода с учетом эффектов релаксации. // Физика металлов и металловедение. - 2007. - Т. 104. - Вып. 3. - С. 235-240.

19. Прудников П.В., Прудников В.В. Влияние структурного беспорядка па аномальное критическое поглощение ультразвука в твердых телах. // Известия вузов. Физика. - 2007. - Вып. 5. - С. 43-49.

20. Прудников П.В., Прудников В.В., Носихин Е.А. Влияние дефектов структуры па аномальны» особенности распространения ультразвука в твердых телах при фазовых переходах второго рода. // ЖЭТФ. - 2008. - Т. 133. - Вып. 5. - С. 1027-1035.

21. Прудников П.В., Прудников В.В., Калашников И.А., Циркин С.С. Ренорм-гругшовое описание процессов неравновесной критической релаксации в коротко-временном режиме: трехпетлевое приближение. // ЖЭТФ. - 2008. - Т. 133. - Вып. 6. - С. 1251-1257.

22. Prudnikov P. V., Prudnikov V.V. Influence of long-range correlated defects on critical ultrasound propagation ill solids. // Pliys. Rev. B. - 2009. - V. 80. - P. 024115-1 - 024115-11.

23. Прудников В.В., Прудников П.В., Вакилов А.Н., Криницын A.C., Рычков М.В. Компьютерное моделирование неравновесной критической динамики структурно неупорядоченных ферромагнетиков. - Научное издание. Труды Семинара по вычислительным технологиям в естественных науках. Вып. 1. Вычислительная физика / Под ред. Р. Р. Назирова.- М. : Изд-во Книжный Дом Университет, 2009. - С. 240-263.

24. Прудников П.В., Прудников В.В., Колесников В.Ю., Медведева М.А., Желтышев П.А. Численное исследование влияния протяженных дефектов структуры на критическое поведение трехмерных систем методом коротковременной динамики. - Научное издание. Труды Семинара по вычислительным технологиям в естественных науках. Вып. 1. Вычислительная физика / Под ред. Р. Р. Назирова,- М. : Изд-во Книжный Дом Университет, 2009. - С. 264-278.

25. Прудников В.В., Вакилов А.Н., Прудников П.В. Фазовые переходы и методы их компьютерного моделирования. Москва: ФИЗМАТЛИТ, 2009. - 224 с.

26. Прудников П.В., Медведева М.А., Желтышев П.А. Численное исследование неравновесного критического поведения трехмерной модели Гейзенберга с линейными дефектами. // Вестник Омского университета. - 2009. - Вып. 4. - С. 108-113.

27. Прудников П.В.. Прудников В.В., Калашников И.А., Рычков М.В. Неравновесная критическая релаксация структурно неупорядоченных систем в коротко-временном режиме: ренорм-групповое описание и компьютерное'моделирование. // ЖЭТФ. - 2010. - Т. 137. - Вып. 2 - С. 287-300.

28. Prudnikov P.V., Prudnikov V.V., Krimtsyn A.S., Vakilov A.N., Pospelov E.A., Rychkov M.V. Short-time dynamics and critical behavior of the tliree-dimensional site-diluted Ising model, // Physical Review E. - 2010. - V. 81. - P. 011130-1 - 011130-11.

29. Прудников П.В., Анкилов H.H., Анкилова Г.А. Теоретико-полевое описание мультикритическо-го поведения n-компонеитных сжимаемых систем. // Вестник Омского университета. - 2010. -Вып. 2. - С. 57-61.

30. Прудников П.В., Яковлев М.И., Бакланов A.B., Воронина А.О., Горохова О.В. Теоретический расчет критических характеристик неупорядоченной системы с дальнодействующей корреляцией дефектов. // Вестник Омского университета. - 2010. - Вып. 2. - С. 62-66.

31. Прудников П.В., Кормилов В.К. Суммирование многопараметрических асимптотических рядов в теории критических явлений методом конформного отображения. // Вестник Омского университета. - 2010. - Вып. 2. - С. 67-70.

32. Прудников П.В., Рычков М.В., Кузнецова Ю.С. Численное исследование неравновесного критического поведения неупорядоченной антиферромагнитной модели Изинга с эффектами случайных полей. /У Вестник Омского университета. - 2010. - Вып. 2. - С. 77-80.

■33. Прудников П.В., Медведева М.А., Желтышев П.А. Численное исследование влияния эффектов корреляции дефектов структуры на критическую динамику модели Гейзенберга. // Вестник Омского университета. - '2010. - Вып. 4. - С. 64-69.

34. Прудников П.В., Куликов Д.Н. Неравновесная критическая динамика структурно неупорядоченных систем с дальнодействующей корреляцией дефектов. // Письма в ЖЭТФ. - 2011. - Т. 93. -Вып. 2-С'. 106-111.

Введение

1 Фазовые переходы второго рода и критические явления

Введение.

1.1 Теория Ландау-Гинзбурга-Вильсона.

1.2 Критические индексы. Гипотеза подобия.

1.3 Метод ренорыгруппы и е - разложения.

1.4 Динамические критические явления.

1.5 Влияние дефектов структуры па критическое поведение.

1.6 Теоретико-полевой подход к описанию критического поведения.

1.6.1 Теоретико-полевой вариант ренормгруппы.

1.6.2 Производящий функционал для функций Грина и вершинных функций

1.6.3 Уравнение ренормгруппы. Асимптотическое поведение функций Грина.

1.7 Суммирование асимптотических рядов.

1.7.1 Суммирование многопараметрических рядов

1.8 Компьютерное моделирование критического поведения спиновых систем

1.8.1 Моделирование методом Монте-Карло канонического ансамбля.

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

1.8.3 Кластерные методы моделирования.

1.8.4 Модификация метода Монте-Карло для неупорядоченных систем.

1.8.5 Динамическая интерпретация процесса моделирования

1.9 Метод реплик и нарушение репличной симметрии.

1.10 Распространение ультразвука вблизи критической температуры.

1.11 Выводы и задачи исследования

2 Теоретико-полевое описание неравновесного критического поведения однородных и неупорядоченных систем с некоррелированными дефектами структуры

Введение.

2.1 Динамическое критическое поведение однородных и неупорядоченных систем

2.2 Теоретико-полевое описание неравновесной критической релаксации однородной системы в трехпетлевом приближении.

2.2.1 Введение.

2.2.2 Модель.

2.2.3 Ренорм-групповой анализ модели.

2.2.4 Анализ результатов.

2.3 Теоретико-полевое описание неравновесной критической релаксации структурно неупорядоченных систем в коротко-временном режиме.

2.3.1 Введение.

2.3.2 Ренормгрупповое описание неравновесного критического поведения структурно неупорядоченных систем.

2.4 Выводы главы.

3 Численные исследования критического поведения неупорядоченных систем с некоррелированными дефектами структуры

Введение.

3.1 Компьютерное моделирование равновесного критического поведения трехмерной неупорядоченной модели Изинга.

3.1.1 Метод конечноразмерного скейлинга.

3.1.2 Расчет критических характеристик.

3.1.3 Анализ результатов и выводы.

3.2 Компьютерное моделирование неравновесного критического поведения неупорядоченной трехмерной модели Изинга с некоррелированными дефектами структуры.

3.2.1 Исследование влияния начального неравновесного состояния с то

3.2.2 Исследование влияния начального неравновесного состояния с mQ = 1 143 3.3 Основные результаты и выводы главы.

4 Эффекты влияния нарушений репличной симметрии на критическое поведение структурно неупорядоченных систем с замороженными дефектами структуры.

Введение.

4.1 Определение модели. Методика расчетов.

4.2 Уравнение Каллана-Симанзика и скейлинговые функции системы.

4.3 Фиксированные точки и различные типы критического поведения.

4.4 Критическое поведение неупорядоченной двумерной модели Изинга с НРС

4.5 Критическое поведение систем с произвольной размерностью d от 3 до

5.2.1 Исследование неравновесной критической динамики модели Изинга с дальнодействующей корреляцией дефектов .212

5.2.2 Компьютерное моделирование равновесного критического поведения неупорядоченной модели Изинга с дальнодействующей корреляцией дефектов.222

5.2.3 Численное исследование неравновесной критической динамики ХУ модели с линейными дефектами.224

5.2.4 Численное исследование неравновесной критической динамики модели Гейзенберга с линейными дефектами .229

5.3 Компьютерное моделирование критического поведения сильно неупорядоченных систем с дальнодействующей корреляцией дефектов .237

5.4 Основные результаты и выводы главы.240

6 Теоретико-полевое описание влияния дефектов структуры и эффектов их корреляции на характеристики распространения ультразвука в твердых телах вблизи температуры фазового перехода второго рода. 243

Введение.243

6.1 Теоретическое описание влияния точечных дефектов структуры на характеристики распространения ультразвука в твердых телах.247

6.1.1 Описание модели.247

6.1.2 Расчет характеристик распространения ультразвука и скейлинговых функций с учетом влияния некоррелированных дефектов структуры 252

6.1.3 Анализ результатов и выводы.254

6.2 Исследование влияния дальнодействующей корреляции дефектов на характеристики распространения ультразвука в твердых телах.257

6.3 Основные результаты и выводы главы.266

7 Теоретико-полевое описание мультикритического поведения однородных и неупорядоченных систем 269

Введение.269

7.1 Теоретико - полевое описание мультикритического поведения однородных систем с двумя параметрами порядка.271

7.2 Исследование влияния неупорядоченности па мультикритическое поведение систем с двумя параметрами порядка.286

7.3 Влияние поверхности на поведение систем в окрестности мультикритической точки Лифшица.293

7.4 Основные результаты и выводы главы.297

Заключение 299

Литература 307

Введение

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

Экспериментальные исследования выявили общность свойств фазовых переходов второго рода в различных веществах. Это позволило сформулировать принцип универсальности критических явлений [10,27,35,45,56,180] и предложить модель, в основе которой лежала гипотеза масштабного подобия флуктуаций [52-54,218]. Идеи использования метода ренормализационной группы и последующая их иллюстрация с помощью метода разложения по отклонению размерности системы от четырех (d = 4 — е) [13,317,319] позволили сделать еще несколько шагов в качественном понимании фазовых переходов и в их количественном описании. Дальнейшее развитие этих идей привело к появлению более надежного теоретико-полевого подхода к описанию критических явлений [16,122,151,257], позволяющему исследовать критическое поведение непосредственно трехмерных систем и дающему более точные количественные результаты при применении методов суммирования асимптотических рядов [122,232].

Большой интерес вызывает проблема исследования влиянии дефектов структуры на критическое поведение. Рассеяние флуктуаций на дефектах структуры, вызывающих нарушение трансляционной инвариантности системы, обусловливает дополнительное взаимодействие флуктуаций параметра порядка через поле дефектов, характеризующееся специфическими законами сохранения. Особенно интересно влияние замороженных дефектов структуры, чье присутствие проявляется или как случайное возмущение локальной температуры (как это происходит, например, в ферромагнитных и антиферромагнитных системах в отсутствие внешнего магнитного поля) или как случайные поля, сопряженные параметру порядка (например, в антиферромагнитных системах в однородном магнитном поле). Наибольших успехов в качественном понимании и количественном описании исследователи достигли при изучении влияния точечных некоррелированных (¿-коррелированных) дефектов с эффектами типа случайной локальной температуры на критическое поведение неупорядоченных систем. Так, в работе Харриса [190] был сформулирован эвристический критерий существенности некоррелированных точечных дефектов, согласно которому присутствие замороженных точечных дефектов изменяет свойства системы вблизи критической точки, если теплоемкость соответствующей однородной системы характеризуется критическим индексом теплоемкости а'о > 0. В противном случае присутствие дефектов не сказывается на значении критических индексов.

Согласно последним исследованиям критических явлений [94.171.216,236,237]. данному критерию удовлетворяют только неупорядоченные системы, эффективный гамильтониан которых вблизи критической точки изоморфен гамильтониану модели Изинга. Ре-нормгрупповой анализ с использованием ег-разложения [105,182,234,242] выявил, что критическое поведение неупорядоченных изингоподобных систем действительно характеризуется новым набором критических индексов. Однако, асимптотическая сходимость рядов с-разложения для систем с дефектами еще более слабая, чем для однородных. Поэтому для их исследования был применен теоретико-полевой подход [94,171,216,236,237.250], в рамках которого были получены более точные значения критических индексов и доказано, что маргинальная размерность параметра порядка ?>,., для которого существенны точечные дефекты, действительно меньше двух (пс < 2) [161,216].

При ренорм-групповом описании критического поведения неупорядоченных систем с замороженным беспорядком для восстановления трансляционной симметрии эффективного гамильтониана, описывающего взаимодействие флуктуации, используется метод реплик [162,163,182]. Однако в ряде работ [21,158,159] были высказаны идеи о возможности нарушения репличной симметрии (НРС) в системах с замороженным беспорядком. Для систем с числом компонент параметра порядка п, меньшем четырех, в рамках метода е

- разложения в низшем порядке теории, было выявлено определяющее влияние эффектов НРС па критическое поведение. Несмотря на столь интересные выводы данных работ результаты проведенных ранее исследований по теоретико-полевому описанию однородных и неупорядоченных систем в двухпетлевом и более высоких порядках приближения с применением методов суммирования асимптотических рядов показали [58,88,89,267], что анализ устойчивости различных типов критического поведения в первом порядке е -разложения можно рассматривать лишь в качестве грубой оценки, особенно для многовершинных статистических моделей [11,288,295,296]. Поэтому результаты исследований эффектов НРС, полученные в работах [21,158,159]. требуют детальной переоценки с позиций применения более точного подхода.

Смещение исследований в современной физике твердого тела в область микромасштабов вызвало необходимость понимания тонких явлений связанных с наличием протяженных дефектов структуры типа дислокаций, границ зерен, примесных комплексов и т.д. Все эти особенности представляют собой проявление пространственно скоррелирован-ных неоднородностей. Вопрос влияния эффектов дальнодействующей пространственной корреляции замороженных дефектов структуры на статическое и, в особенности, на динамическое критическое поведение несмотря на его важность при описании критического поведения реальных неупорядоченных систем мало исследован, а полученные результаты носят оценочный характер [18,20,223,225,314]. Поэтому существует потребность в развитии более надежных методов описания критического поведения систем с дальнодействующей корреляцией дефектов структуры и получении более точных сведений об условиях устойчивости различных типов критического поведения таких систем и характеристиках этого поведения.

Важным представляется исследование влияния неупорядоченности, создаваемой присутствием примесей, на характер фазовых диаграмм систем в окрестности мульти-критических точек. Теоретические модели, соответствующие данным системам, являются многовершинными. Исследование возможных типов устойчивого мультикритического поведения, описываемых данными моделями, в рамках метода £ - разложения [37,42,212] в однопетлевом приближении нельзя считать достоверными. Уже при исследовании мультикритического поведения однородной системы в [58] было наглядно показано слабое соответствие предсказаний однопетлевого приближения реальному мультикритическому поведению. В случае неупорядоченных систем можно ожидать еще более существенных отличий. Поэтому необходимо развитие теоретико-полевого описания мультикритическо-го поведения неупорядоченных систем в более высоких порядках приближения теории с применением эффективных математических методов для суммирования асимптотических рядов, получаемых в реальном пространстве в различных многопетлевых приближениях.

В связи с этим целыо настоящей диссертации является:

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

2. Численное исследование критического поведения неупорядоченной ферромагнитной модели Изинга с дефектами типа случайная температура численными методами Монте-Карло. Ставится задача провести:

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

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

- численное определение с использованием методов Монте-Карло равновесных критических индексов и критических температур для трехмерной модели Изинга в широком диапазоне изменения спиновых концентраций (0.95 < р < 0.50). применяя процедуру конеч-норазмерного скейлинга с учетом асимптотических поправок к скейлингу;

- выяснение вопроса об универсальности неравновесного критического поведения неупорядоченных систем.

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

4. Теоретико-полевое и численное исследование статического и динамического критического поведения неупорядоченных трехмерных систем с дальнодействующей пространственной корреляцией замороженных дефектов структуры. Предполагается осуществить:

- разработку теоретико-полевого описания трехмерных систем с изотропной дальнодействующей пространственной корреляцией дефектов структуры в двухпетлевом приближении с применением методов суммирования асимптотических рядов;

- определение фиксированных точек ренормгрупповых преобразований модели и условий устойчивости различных типов критического поведения;

- вычисление статических критических индексов и динамического критического индекса;

- исследование влияния дальнодействующей корреляции дефектов на критическое поведение систем с различным числом компонент параметра порядка посредством численного изучения методами Монте-Карло ферромагнитных трехмерных моделей Изинга и ХУ;

- численное исследование неравновесного критического поведения трехмерных моделей Изинга и ХУ с линейными дефектами при спиновых концентрациях р = 0.80 и р = 0.60 методом коротковременной динамики при рассмотрении эволюции систем из разных начальных неравновесных состояний.

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

6. Ренормгрупповое исследование мультикритического поведения сложных однородных и структурно неупорядоченных систем с двумя взаимодействующими параметрами порядка и выявление влияния поверхности на мультикритическое поведение в окрестности точки Лифшица. Ставится задача провести:

- развитие методики и осуществление теоретико-полевого описания мультикритического поведения трехмерных однородных систем с двумя параметрами порядка;

- исследование влияния структурной неупорядоченности на характер фазовых диаграмм и свойства мультикритического поведения трехмерных систем с двумя взаимодействующими параметрами порядка;

- исследование влияния поверхности и ее ориентации относительно направления анизотропии кристалла на мультикритическое поведение системы в окрестности обобщенной точки Лифшица ?тт,-го порядка.

Результаты диссертации были опубликованы в работах [58-85,266-279]

7.4 Основные результаты и выводы главы

В заключение можно выделить следующие основные результаты данной главы.

1. Теоретико-по левое описание мультикритического поведения непосредственно трехмерных систем с двумя параметрами порядка выявило в двухпетлевом приближении существенное изменение областей различного типа устойчивого мультикритического поведения на плоскости (п — т) - числа компонент данных параметров порядка по сравнению с полученными ранее результатами. Показано, что устойчивое мультикритическое поведение, соответствующее изотропной фиксированной точке с флуктуационно индуцированной асимптотической симметрией системы Б О (п + т), возможно только для п + т < 2.9088, т.е. наивысшей асимптотической симметрией системы является 80(2). Значительное изменение претерпели области стабильности и остальных двух типов фиксированных точек, соответствующих устойчивому тетракритическому поведению.

2. Показано, что изменение областей стабильности фиксированных точек приводит к заметному изменению типов фазовых диаграмм систем во флуктуационной области. Устойчивое бикритическое поведение предсказывается в диссертации только для взаимодействующих однокомпонентных параметров порядка (п = т = 1) с критическими индексами, соответствующими изотропной XV-модели. Тетракритическое же поведение должно иметь более широкую реализацию среди систем с многокомпонентными параметрами порядка. Обсуждены эффекты флуктуационной неустойчивости мультикритического поведения.

3. Исследование влияния точечных замороженных примесей, создающих в системах с двумя параметрами порядка флуктуации случайной локальной температуры, на характер фиксированных точек и их устойчивость показало, присутствие примесей в системе приводит к флуктуационному расцеплению связи параметров порядка и осуществлению лишь единственного типа устойчивого мультикритического поведения - тетракритическо-го с общей симметрией системы БО (п) ® Б О (т).

В случае однокомпонентных параметров порядка (п = т = 1) наличие примесей существенно и приводит к критическому поведению с индексами, соответствующими индексам неупорядоченной модели Изинга. Когда параметры порядка системы характеризуются числом компонент большим или равным двум, присутствие примесей не сказывается на характеристиках их критического поведения, а мультикритическое поведение носит тетракритический характер однородной системы.

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

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

5. Впервые исследовано влияние поверхности, перпендикулярной оси анизотропии кристалла, на мультикритическое поведение ¿-мерной системы вблизи точки Лифшица 771-го порядка. Полученное значение нового независимого поверхностного критического индекса /3~ = 0.59 подтверждает результаты численного исследования модели Изинга для случая точки Лифшица с т, п = 1.

Заключение

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

1. В рамках единого подхода, основанного на формализме динамического производящего функционала для корреляционных функций и функций отклика, осуществлено теоретико-полевое описание неравновесного критического поведения с учетом влияния начальных состояний системы, нарушающих трансляционную инвариантность во времени, и квазиравновесной критической динамики, когда влияние неравновесных начальных состояний самоусреднилось, для однородных и неупорядоченных систем с некоррелированными дефектами и дефектами с дальнодействующей пространственной корреляцией.

2. Впервые показано, что в неравновесной критической релаксации однородных систем только начиная с трехпетлевого приближения в динамической функции отклика возникают флуктуационные поправки за счет влияния начальных неравновесных состояний. С учетом данных флуктуационных эффектов проведен расчет независимого динамического критического индекса в', задающего эволюцию п-компонентного параметра порядка системы в коротковременном режиме, и при применении метода е-разложения были получены значения в', согласующиеся с результатами компьютерного моделирования.

3. Впервые реализовано теоретико-полевое описание неравновесного критического поведения структурно неупорядоченных систем с некоррелированными дефектами с учетом влияния начальных состояний при фиксированной размерности системы с1 = 3 и в двухпетлевом приближении проведен расчет динамических критических индексов в' и 2 с последовательным применением к рядам теории различных методов суммирования. Численные значения критических индексов находятся в лучшем соответствии с результатами компьютерного моделирования, чем результаты применения метода е-разложения.

4. Впервые осуществлено последовательное применение методов суммирования Паде-Бореля, Паде-Бореля-Лероя и конформного Паде-Бореля к рядам теории, полученным в многопетлевом приближении с целью определения значений динамического критического индекса 2 для однородных трехмерных и двумерных систем, описываемых моделью Изинга, а также для неупорядоченных изингоподобных систем. Вычисленные значения индекса z находятся в хорошем согласии с результатами экспериментальных исследований критической динамики изинговских антиферромагнетиков и результатами компьютерного моделирования критической динамики.

5. Разработана методика теоретико-полевого описания критического поведения структурно неупорядоченных спиновых систем с эффектами нарушения репличной симметрии (НРС) при их фиксированной размерности. На примере данных систем, описываемых многовершинными моделями, показано, что предсказания, основанные на применении метода е- разложения, не являются надежными. Конкуренция различных типов фиксированных точек в многопараметрическом пространстве вершин взаимодействия не позволяет осуществлять протяжку £-—>1,2 без пересечения областей стабильности различных фиксированных точек.

6. Для двумерных, трехмерных систем и систем произвольной размерности от трех до четырех в двухпетлевом приближении проведен ренормгрупповой анализ эффективного репличного гамильтониана модели с потенциалом взаимодействия, не являющимся реплично-симметричным и определяемым тремя вершинами взаимодействия флуктуаций. Для случая одноступенчатого НРС с применением техники суммирования Паде-Бореля были выделены возможные типы критического поведения и осуществлен анализ возможности их реализации. Было установлено, что критическое поведение двумерных и трехмерных систем устойчиво относительно влияния эффектов НРС. Показано, что наличие слабого беспорядка не влияет на критическое поведение систем с многокомпонентным параметром порядка, а в системах с однокомпонентным параметром порядка реализуется критическое поведение, определяемое структурным беспорядком с реплично-симметричной фиксированной точкой.

7. Выявлено: эффекты НРС проявляются лишь при размерностях пеупорядоченной системы больших трех, при этом пороговые значения dc зависят от числа компонент параметра порядка п и величины параметра xq, определяющего координату ступеньки для потенциала взаимодействия с НРС. Определены значения пороговых размерностей dc(n), отделяющих область критического поведения с эффектами НРС dc(n) < d < 4 от области, в которой данные эффекты несущественны. Эти пороговые размерности задают одновременно и нижнюю границу области применимости результатов е-разложения к описанию модели слабо неупорядоченных систем с эффектами НРС. Показано, что возможные типы устойчивого критического поведения и критические индексы, полученные для неупорядоченных систем в рамках стандартного метода реплик для размерностей систем ниже пороговых dc(n), являются достоверными и реализуется прежний сценарий влияния замороженного беспорядка на критическое поведение.

8. Осуществлено теоретико-полевое описание статических свойств критического поведения и неравновесной критической динамики неупорядоченных трехмерных систем с дальнодействующей корреляцией дефектов в двухпетлевом приближении с использованием методов пересуммирования рядов теории. Для различных значений числа компонент параметра порядка и показателя корреляции а определены типы устойчивого критического поведения и значения критических индексов. Полученная картина областей устойчивого критического поведения и значения критических индексов существенно отличаются от предсказанных ранее в рамках двухпараметрического е,5 - разложения. Показано, что корреляция дефектов приводит к проявлению влияния неупорядоченности в поведении более широкого круга систем, вызывая существенное изменение статических и динамических характеристик критического поведения. В частности, с усилением корреляции дефектов происходит замедление процессов критической релаксации в системе и усиление влияния неравновесных начальных состояний по сравнению с однородными системами и системами с некоррелированными дефектами.

9. Осуществлены численные Монте-Карло исследования равновесного критического поведения трехмерной спиновой модели Изинга в широком интервале изменения концентрации немагнитных атомов примеси. Для образцов с различными спиновыми концентрациями и различными линейными размерами была определена температурная зависимость корреляционной длины и магнитной восприимчивости. С помощью метода конечноразмер-ного скейлинга для данных величин были определены скейлинговые функции, демонстрирующие универсальное поведение в критической области, а также вычислены критические температуры и статические критические индексы с использованием ведущих поправок к асимптотической зависимости данных величин. На основе выявленных концентрационных изменений скейлинговых функций и значений критических индексов сделан вывод о существовании двух универсальных классов критического поведения разбавленной модели Изинга с различными характеристиками для слабо и сильно неупорядоченных систем.

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

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

11. На основе проведенных численных исследований методом коротковременной динамики критического поведения трехмерной модели Изинга со спиновыми концентрациями р = 0.95, 0.80, 0.60 и 0.50 можно сделать следующие выводы: поведение слабо и сильно неупорядоченных систем принадлежит к различным универсальным классам критического поведения с несовпадающими в пределах статистических погрешностей значениями динамических критических индексов в' и г; в неравновесном поведении слабо неупорядоченных систем с р = 0.95 и 0.80 выявлено два режима со степенным временным изменением измеряемых величин: на раннем временном интервале реализуется критическое поведение с характеристиками однородной системы, и затем, проходя через интервал крос-соверного поведения реализуется режим поведения неупорядоченной системы; в неравновесном поведении сильно неупорядоченных систем с р = 0.60 и 0.50 не выявлено режима с характеристиками однородной системы; для слабо неупорядоченных систем значения критических индексов, полученные с учетом ведущих поправок к скейлингу, находятся в хорошем согласии с результатами теоретико-полевого описания, результатами моделирования другими методами, а также согласуются с результатами экспериментальных исследований слабо неупорядоченных изинговских магнетиков; для сильно неупорядоченных систем значения динамического критического индекса г, полученные при анализе эволюции из разных начальных неравновесных состояний, находятся в хорошем согласии друг с другом, согласуются со значениями, определенными другими численными методами, а значение индекса коротковременной эволюции намагниченности в' = 0.186(39) получено впервые.

12. С целью исследования влияния дальнодействующей корреляции дефектов структуры на критическое поведение впервые осуществлено численное исследование методом коротковременной динамики поведение трехмерных модельных спиновых систем (Изинга, ХУ и Гейзенберга) с линейными дефектами при общей спиновой концентрации р = 0.80. В неравновесном критическом поведении данных систем выявлено два динамических режима со степенным временным изменением измеряемых величин: на раннем временном интервале реализуется поведение с характеристиками однородной системы, а после интервала кроссоверного поведения реализуется динамический режим неупорядоченной системы. Определены значения совокупности динамических и статических критических индексов при применении методики учета ведущих поправок к скейлингу, соответствующих режиму неупорядоченной системы. Полученные значения критических индексов находятся в хорошем согласии с результатами теоретико-полевого описания на основе модели Вейнриба-Гальперина с показателем корреляции а = 2 и результатами проведенного моделирования равновесного критического поведения. На примере исследования ХУ-модели и модели Гейзенберга впервые получено численное подтверждение существенности влияния дальней пространственной корреляции дефектов на критическое поведение систем с многокомпонентным параметром порядка.

13. Впервые осуществлено численное исследование неравновесного критического поведения трехмерных сильно неупорядоченных при р = 0.60 спиновых моделей Изинга и ХУ с дальнодействующей корреляцией дефектов, моделируемой изотропно распределенными линейными дефектами. Сопоставление полученных значений критических индексов со значениями для слабо неупорядоченных с р = 0.80 моделей Изинга и ХУ позволяют сделать вывод о существовании различных универсальных классов критического поведения для рассматриваемых систем, соответствующих областям слабой и сильной структурной неупорядоченности.

14. Осуществлено теоретико-полевое описание аномальных свойств распространения ультразвука в однородных и структурно неупорядоченных твердых телах вблизи температуры фазового перехода второго рода с учетом как флуктуационного, так и релаксационного механизмов рассеяния. Разработана методика и проведен расчет скейлинговых функций для коэффициентов поглощения и дисперсии скорости звука, а также самих коэффициентов при фиксированной размерности системы й — 3 в двухпетлевом приближении с применением методов суммирования асимптотических рядов. Представлены численные значения показателей, характеризующих частотную и температурную зависимости скейлинговых функций, коэффициентов поглощения и дисперсии ультразвука в гидродинамической (10~3 < у = шг-2"/Го < Ю-1), предкритической (1 < у < 102) и критической (10 < у < Ю3) областях.

Рассчитанное для однородных систем поведение коэффициента поглощения а(ш, г) демонстрирует аномально сильное поглощение ультразвука в критической области и хорошо согласуется с результатами экспериментальных исследований в образцах ГеГ2.

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

15. Впервые исследовано влияние эффектов дальнодействующей корреляции дефектов на аномальное рассеяние ультразвука в твердых телах при фазовых переходах второго рода. Показано, что усиление корреляции дефектов приводит к увеличению поглощения ультразвука и усилению дисперсии скорости звука в критической области. Предсказывается более сильное, чем для однородных систем и систем с некоррелированными дефектами, температурное изменение характеристик ультразвука при приближении к критической температуре уже в гидродинамической области. В результате, экспериментальное исследование критической динамики ультразвуковыми методами позволяет выявить влияние дефектов и эффектов их корреляции в более широком температурном интервале относительно критической температуры (Ю-3 < т < Ю-1), чем в других экспериментальных методах, в которых для выявления данных эффектов необходимо проводить исследования в узком температурном интервале вплоть до т ~ Ю-4.

16. Теоретико-полевое описание мультикритического поведения однородных систем с двумя параметрами порядка в рамках метода с фиксированной й = 3 размерностью выявило в двухпетлевом приближении существенное изменение областей различного типа устойчивого мультикритического поведения на плоскости (п — т) - числа компонент данных параметров порядка по сравнению с результатами применения е - разложения. Это приводит к изменению характеристик мультикритического поведения и возможных типов фазовых диаграмм системы во флуктуационной области.

Показано, -что присутствие дефектов структуры в системе приводит к флуктуацион-ному расцеплению связи параметров порядка и осуществлению единственного типа устойчивого мультикритического поведения - тетракритического. В случае однокомпонентных параметров порядка (?г, т = 1) наличие дефектов существенно и приводит к мультикри-тическому поведению с индексами, соответствующими индексам неупорядоченной модели Изинга. Для систем с п, т > 2 присутствие дефектов не сказывается на мультикрити-ческом поведении и оно носит тетракритический характер однородной системы. Влияние дефектов проявляется в сокращении по сравнению с однородными системами возможных типов фазовых диаграмм.

17. Впервые исследовано влияние поверхности, перпендикулярной оси анизотропии кристалла, на мультикритическое поведение ¿-мерной системы вблизи точки Лифшица 771-го порядка. Полученное значение нового независимого поверхностного критического индекса ¡3х = 0.59 подтверждает результаты численного исследования модели Изинга для случая точки Лифшица с га, п = 1.

Разработанные в диссертации методы и полученные результаты вносят существенный вклад в обоснование и развитие представлений теории критических явлений в неупорядоченных спиновых системах.

1. Александров К.С., Анистратов А.Т., Безносиков Б.В., Федосеева Н.В. Фазовые переходы в кристаллах галоидных соединений АВХ3. - Новосибирск: Наука. - 1981.

2. Алиев Х.К., Камилов И.Х., Омаров A.M. Критическая динамика гадолиния. -ЖЭТФ. 1989. - Т. 95. - С. 1896-1907.

3. Анисимов М.А., Городецкий Е.Е., Запрудский В.М. Фазовые переходы с взаимодействующими параметрами порядка. УФН. - 1981. Т. 133. - Вып. 1. - С. 103-137.

4. Антонов Н.В., Васильев А.Н. Критическая динамика как теория поля. ТМФ. - 1984.- Т. 60. С. 59-71.

5. Аплеснин С.С., Москвин А.И. Влияние сильных электронных корреляций и взаимодействия электронов с решеткой на орбитальное упорядочение электронов. Письма в ЖЭТФ - Т. 92. - Вып. 4. - Р. 254-259.

6. Аплеснин С.С. Исследование магнитных свойств слабовзаимодействующих антиферромагнитных цепочек с альтернированным обменным взаимодействием со спином s = 1/2 при помощи квантового метода Монте-Карло. ЖЭТФ. - 2000. - Т. 117.- Вып. 1. С. 218-226.

7. Бейкер Г. Аппроксимация Паде. М.: Мир. - 1986. - 336С.

8. Боголюбов H.H., Ширков Д.В. Введение в теорию квантованных полей. М.: Наука.- 1984. 540С.

9. Вакилов А.Н., Прудников В.В. Компьютерное моделирование критической динамики разбавленных магнетиков. Письма в ЖЭТФ. - 1992. - Т. 55. - Вып. 12. - С. 709-712.

10. Вакс В.Г., Ларкин А.И. О фазовых переходах второго рода. ЖЭТФ. - 1965. - Т. 49.- Вып. 3. С. 975-989.

11. Варнашев К.Б., Соколов А.И. Критическая термодинамика кубических и тетрагональных кристаллов с многокомпонентными параметрами порядка. ФТТ. - 1996. -Т. 38. - С. 3665.

12. Васильев А.Н. Квантовополевая ренорм-группа в теории критического поведения и стохастической динамике. Санкт-Петербург: ПИЯФ. - 1998.

13. Вильсон К., Когут Дж. Ренормализационная группа и ^-разложение. М.: Мир. -1975. - 256С.; - УФН. - 1985. - Т. 146. - Вып. 3. - С. 459-491.

14. Владимиров A.A., Казаков Д.И., Тарасов О.В. О вычислении критических индексов методами квантовой теории поля. ЖЭТФ. - 1979. - Т. 77. - Вып. 3. - С. 1035-1045.

15. Гинзбург В.Л. Несколько замечаний о фазовых переходах второго рода и микроскопической теории сегнетоэлектриков. ФТТ. - 1960. - Т. 2. - Вып. 9. - С. 2034-2043.

16. Гинзбург С.Л. Определение фиксированной точки и критических индексов. ЖЭТФ.- 1975. Т. 68. - Вып. 1. - С.273-286.

17. Дейген М.Ф., Глинчук М.Д. Параэлектрический резонанс нецентральных ионов. -УФН. 1974. - Т. 114. - Вып. 2. - С.185-211.

18. Дороговцев С.Н. Фазовый переход в системе с протяженными дефектами. ФТТ. -1980. - Т. 22. - Вып. 2. - С. 321-327.

19. Дороговцев С.Н. Критические свойства систем с протяженными дефектами. Анизотропия критических индексов. ФТТ. - 1980. - Т. 22. - Вып. 12. - С. 3658-3664.

20. Дороговцев С.Н. Критические свойства магнетиков с дислокациями и точечными примесями. ЖЭТФ. - 1981. - Т. 80. - Вып. 5. - С. 2053-2067.

21. Доценко B.C. Критические явления в спиновых системах с беспорядком. УФН. -1995. - Т. 165. - Вып. 5. - С. 481-528.

22. Доценко B.C. Физика спин-стекольного состояния. УФН. - 1993. - Т. 163. - Вып. 6.- С.1-37.

23. Займан Дж. Модели беспорядка. Теоретическая физика однородно неупорядоченных систем. М.: Мир. - 1982. - 591С.

24. Иванченко Ю.М., Лисянский A.A., Филиппов А.Э. Флуктуационные эффекты в системах с конкурирующими взаимодействиями. Киев: Наука думка. - 1989. - 280С.

25. Изюмов Ю.А., Сыромятников В.Н. Фазовые переходы и симметрия кристаллов. М.: Наука. - 1984. - 248С.

26. Кавасаки К. Динамическая теория флуктуаций вблизи критических точек . Квантовая теория поля и физика фазовых переходов. М.: Мир. - 1975. - С. 101-148.

27. Каданов Л.П. Критические явления, гипотеза универсальности, скейлинг и капельная модель. Квантовая теория поля и физика фазовых переходов. М.: Мир. - 1975. -С. 7-32.

28. Казаков Д.И., Тарасов О.В., Ширков Д.В. Аналитическое продолжение результатов теории возмущений модели gip4 в область д > 1. ТМФ. - 1979. - Т. 38. - С. 15-25.

29. Казаков Д.И., Попов B.C. О суммировании расходящихся рядов теории возмущений в квантовой механике и теории поля. ЖЭТФ. - 2002. - Т. 122. - С. 675.

30. Камилов И.К., Алиев Х.К. Исследование критической динамики магнитоупорядочен-ных кристаллов ультразвуковыми методами. УФН. - 1998. - Т. 168. - С. 953.

31. Камилов И.К., Муртазаев А.К., Алиев Х.К. Исследование фазовых переходов и критических явлений методами Монте-Карло. УФН. - 1999. - Т. 169. - С. 773-795.

32. Корженевский А.Л. Регулярные крупномасштабные сверхструктуры вблизи фазовых переходов в кристаллах. ФТТ. - 1984. - Т. 26. - Вып. 4. - С. 1223-1225.

33. Ландау Л.Д. К теории фазовых переходов. ЖЭТФ. - 1937. - Т. 7. - Вып. 1. - С.19.

34. Ландау Л.Д., Халатников И.М. Об аномальном поглощении звука вблизи точек перехода второго рода. ДАН СССР. - 1954. - Т. 96. - С. 469.

35. Ландау Л.Д., Лифшиц Е.М. Статистическая физика. М.: Наука. - 1976. - 584С.

36. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. М.: Наука, 1988. - 736С.

37. Лаптев В.М., Скрябин Ю.Н. Фазовые диаграммы разупорядоченных систем со связанными параметрами порядка. ФТТ. - 1980. - Т. 22. - Вып. 10. - С. 2949-2955.

38. Ларкин А.И., Пикин С.А. О фазовых.переходах первого рода, близких ко второму. -ЖЭТФ. 1969. - Т. 56. - С. 1664-1674.

39. Леванюк А.П. К теории рассеяния света вблизи точек фазового перехода второго рода. ЖЭТФ. - 1959. - Т. 36. - Вып. 3. - С. 810-818.

40. Леванюк А.П., Собянин A.A. О фазовых переходах второго рода без расходимостей во вторых производных термодинамического потенциала. Письма в ЖЭТФ. - 1970.- Т. И. Вып. И. - С. 540-543.

41. Липатов Л.Н. Расходимость ряда теории возмущений и квазикласика. ЖЭТФ. -1977. - Т.72. - Вып. 2. - С. 411-427.

42. Лисянский A.A., Филиппов А.Э. Критическая термодинамика примесных систем со связанными флуктуирующими полями. УФЖ. - 1987. - Т. 32. - Вып. 4. - С. 626-634.

43. Люксютов И.Ф., Покровский В.Л., Хмельницкий Д.Е. Пересечение линий переходов второго рода. ЖЭТФ. - 1975. - Т. 69. - Вып. 5. - С. 1817-1824.

44. Ма Ш. Современная теория критических явлений. М.: Мир. - 1980. - 298С.

45. Мигдал A.A. Диаграммная техника вблизи точки Кюри и фазовый переход в бозе-жидкости. ЖЭТФ. - 1968. - Т. 55. - Вып. 5. - С. 1964-1979.

46. Муртазаев А.К., Камилов И.К., Магомедов М.А. Кластерные алгоритмы метода Монте-Карло, конечно-размерный скейлинг и критические индексы сложных решеточных моделей. ЖЭТФ. - 2001. - Т. 120. - Вып. 6. - С. 1535

47. Муртазаев А.К., Камилов И.К., Бабаев А.Б. Критическое поведение трехмерной модели Изинга с вмороженным беспорядком на кубической решетке. ЖЭТФ. - 2004.- Т. 126. С. 1377-1383.

48. Муртазаев А.К., Мутайламов В.А. Исследование динамического критического поведения моделей ферромагнитного гадолиния. ЖЭТФ. - 2005. - Т. 128. - Вып. 2. -С. 344.

49. Муртазаев А.К., Рамазанов М.К., Бадиев М.К. Статическое критическое поведение трехмерной фрустрированной модели гейзенберга на слоистой треугольной решетке с переменным межслойным обменным взаимодействием. ЖЭТФ. - 2007. - Т. 132. -Вып. 5. - С. 1152.

50. Муртазаев А.К., Бабаев А.Б., Азнаурова Г.Я. Особенности фазовых переходов в трехмерных разбавленных структурах, описываемых моделью Поттса. ЖЭТФ. - 2009.- Т. 136. Вып. 3. - С. 516.

51. Найш В.Е., Скрябин Ю.Н., Сыромятников В.II. Фазовые переходы с взаимодействующими параметрами порядка в соединениях NiAs типа. ФММ. - 1981. - Т. 52. -Вып. 6. - С. 1147-1155.

52. Паташинский А.З. Гипотеза подобия в теории фазовых переходов второго рода. -ЖЭТФ. 1967. - Т. 53. - Вып. 6. - С. 1987-1996.

53. Паташинский А.З., Покровский B.JI. Фазовый переход второго рода в бозе-жидкости.- ЖЭТФ. 1964. - Т. 46. - Вып. 3. - С. 994-1016.

54. Паташинский А.З., Покровский B.JI. О поведении упорядочивающихся систем вблизи точки фазового перехода. ЖЭТФ. - 1966. - Т. 50. - Вып. 2. - С. 439-447.

55. Паташинский А.З., Покровский В.Л. Флуктуационная теория фазовых переходов. -М.: Наука. 1982. - 383С.

56. Поляков A.M. Микроскопическое описание критических явлений. ЖЭТФ. - 1968. -Т. 55. - Вып. 3. - С. 1026-1038.

57. Поляков A.M. Свойства далеких и близких корреляций в критической области. -ЖЭТФ. 1969. - Т. 57. - Вып. 1 - С. 271-284.

58. Прудников П.В., Прудников В.В., Федоренко A.A. Мультикритическое поведение слабо неупорядоченных систем с двумя параметрами порядка. ЖЭТФ. - 1999. - Т. 116.- Вып. 2. С. 611-619.

59. Прудников П.В., Прудников В.В., Федоренко A.A. Мультикритическое поведение неупорядоченных систем с двумя параметрами порядка. Физика твердого тела. -2000. - Т. 42. - Вып. 1. - С. 158-162.

60. Прудников П.В., Прудников В.В., Федоренко A.A. Устойчивость критического поведения слабо неупорядоченных систем к нарушению репличной симметрии. Письма в ЖЭТФ. - 2001. - Т. 73. - Вып. 3. - С. 153-158.

61. Прудников П.В., Прудников В.В., Федоренко A.A. Устойчивость критического поведения слабо неупорядоченных систем к введению потенциала взаимодействия с нарушенной репличной симметрией. Физика твердого тела. - 2001. - Т. 43. - Вып. 9. -С. 1688-1692.

62. Прудников П.В., Прудников В.В. Критическое поведение неупорядоченных систем с НРС. Вестник Омского университета. - 2001. - Вып. 3. - С. 26-28.

63. Прудников П.В., Прудников В.В. Критическое поведение неупорядоченных систем с эффектами нарушения репличной симметрии. ЖЭТФ. - 2002. - Т. 122. - Вып. 3. -С. 636-646.

64. Прудников П.В, Носихин Е.А. Особенности описания аномальных свойств распространения ультразвука при фазовых переходах. Вестник Омского Университета. -2005. - Вып. 4. - С. 42-45.

65. Прудников П.В., Прудников В.В., Криницын A.C. Расчет динамического критического индекса методом суммирования асимптотических рядов. Теоретическая и математическая физика. - 2006. - Т. 147. - Вып. 1. - С. 138-155.

66. Прудников В.В., Прудников П.В., Вакилов А.Н., Криницын A.C. Компьютерное моделирование критического поведения трехмерной неупорядоченной модели Изинга. -ЖЭТФ. 2007. - Т. 132. - Вып. 2. - С. 417-425.

67. Прудников П.В., Прудников В.В., Носихин Е.А. Аномальное поглощение ультразвука в твердых телах при фазовых переходах второго рода с учетом эффектов релаксации. Физика металлов и металловедение. - 2007. - Т. 104. - Вып. 3. - С. 235-240.

68. Прудников П.В., Прудников В.В. Влияние структурного беспорядка на аномальное критическое поглощение ультразвука в твердых телах. Известия вузов. Физика. -2007. - Вып. 5. - С. 43-49.

69. Прудников П.В., Прудников В.В., Вакилов А.Н., Криницын A.C. Исследование неравновесной критической релаксации в трехмерной неупорядоченной модели Изинга. -Вестник Омского университета. 2007. - Вып. 3. - С. 15-19.

70. Прудников П.В., Гергертд Е.А., Колесников В.Ю., Прудников В.В. Вестник Омского университета. - 2007 - Вып. 4. - С. 32.

71. Прудников П.В., Прудников В.В., Носихин Е.А. Влияние дефектов структуры на аномальные особенности распространения ультразвука в твердых телах при фазовых переходах второго рода. ЖЭТФ. - 2008. - Т. 133. - Вып. 5. - С. 1027-1035.

72. Прудников П.В., Прудников В.В., Калашников И.А., Циркин С.С. Ренорм-групповое описание процессов неравновесной критической релаксации в коротко-временном режиме: трехпетлевое приближение. ЖЭТФ. - 2008. - Т. 133. - Вып. 6. - С. 1251-1257.

73. Прудников П.В, Прудников В.В., Носихин Е.А. Влияние дальнодействуюгцей корреляции дефектов на аномальное критическое поглощение ультразвука в твердых телах. Вестник Омского университета. - 2008. - Вып. 2. - С. 37-41.

74. Прудников П.В., Прудников В.В., Поспелов Е.А., Рычков М.В., Шляхтин А.О. Исследование неравновесной критической эволюции структурно неупорядоченной трехмерной модели Изинга. Вестник Омского университета. - 2008. - Вып. 4. - С. 35-39.

75. Прудников В.В., Вакилов А.Н., Прудников П.В. Фазовые переходы и методы их компьютерного моделирования. Москва: ФИЗМАТЛИТ, 2009. 224 с.

76. Прудников П.В., Медведева М.А., Желтышев П.А. Численное исследование неравновесного критического поведения трехмерной модели Гейзенберга с линейными дефектами. Вестник Омского университета. - 2009. - Вып. 4. - С. 108-113.

77. Прудников П.В., Анкилов H.H., Анкилова Г.А. Теоретико-полевое описание муль-тикритического поведения тг-компонентных сжимаемых систем. Вестник Омского университета. - 2010. - Вып. 2. - С. 57-61.

78. Прудников П.В., Яковлев М.И., Бакланов A.B., Воронина А.О., Горохова О.В. Теоретический расчет критических характеристик неупорядоченной системы с дальнодей-ствующей корреляцией дефектов. Вестник Омского университета. - 2010. - Вып. 2. - С. 62-66.

79. Прудников П.В., Кормилов В.К. Суммирование многопараметрических асимптотических рядов в теории критических явлений методом конформного отображения. -Вестник Омского университета. 2010. - Вып. 2. - С. 67-70.

80. Прудников П.В., Рычков М.В., Кузнецова Ю.С. Численное исследование неравновесного критического поведения неупорядоченной антиферромагнитной модели Изинга с эффектами случайных полей. Вестник Омского университета. - 2010. - Вып. 2. -С. 77-80.

81. Прудников П.В., Медведева М.А., Желтышев П.А. Численное исследование влияния эффектов корреляции дефектов структуры на критическую динамику модели Гейзенберга. Вестник Омского университета. - 2010. - Вып. 4. - С. 64-69.

82. Прудников П.В., Куликов Д.Н. Неравновесная критическая динамика структурно неупорядоченных систем с дальнодействующей корреляцией дефектов. Письма в ЖЭТФ. - 2011. - Т. 93. - Вып. 2. - С. 106-111.

83. Прудников В.В., Вакилов А.Н. Критическая динамика разбавленных магнетиков. -ЖЭТФ. 1992. - Т. 101. - Вып. 6. - С. 1853-1861.

84. Прудников В.В., Вакилов А.Н. Компьютерное моделирование критической динамики разбавленных магнетиков. ЖЭТФ. - 1993. - Т. 103. - Вып. 3. - С. 962-969.

85. Прудников В.В., Иванов A.B., Федоренко A.A. Критическая динамика спиновых систем в четырехпетлевом приближении. Письма в ЖЭТФ. - 1997. - Т. 66. - Вып. 12.- С. 793-798.

86. Прудников В.В., Белим C.B., Иванов A.B., Осинцев Е.В., Федоренко A.A. Критическая динамика слабо неупорядоченных спиновых систем. ЖЭТФ. - 1998. - Т. 114.- Вып. 3. С. 972-984.

87. Прудников В.В., Белим C.B., Осинцев Е.В., Федоренко A.A. Критическая динамика неупорядоченных магнетиков в трехпетлевом приближении. ФТТ. - 1998. -.Т. 40. -Вып. 8. - С. 1526-1531.

88. Прудников В.В., Белим C.B. Трикритическое поведение сжимаемых систем с замороженными дефектами структуры. ФТТ. - 2001. - Т. 43. - С. 1299.

89. Райдер JI. Квантовая теория поля. М. :Мир. - 1987. - 512С.

90. Сиперли Д., Кейлос М. Квантовые многочастичные задачи. Методы Монте-Карло в статистической физике. - Москва: Мир, 1982. - С. 162-219.

91. Соколов А.И., Шалаев Б.Н. О критическом поведение модели Изинга с примесями. -ФТТ. 1981. - Т. 23. - Вып. 7. - С. 2058-2063.

92. Стенли Г. Фазовые переходы и критические явления. М.: Мир. - 1973.-342С.

93. Стоунхэм A.M. Теория дефектов в твердых телах. М.: Мир. - 1978. - Т. 1. - 569С.

94. Суслов И.М. Суммирование расходящихся рядов теории возмущений в пределе сильной связи. Функция Гелл-Манна-JIoy теории с/Л ЖЭТФ. - 2001. - Т. 120. - С. 5-30.

95. Суслов И.М. Расходящиеся ряды теории возмущений. ЖЭТФ. - 2005. - Т. 127. -С. 1350-1402.

96. Суслов И.М., Погорелов A.A. Ренормгрупповые функции для двумерных фазовых переходов: к вопросу о сингулярных вкладах. ЖЭТФ. - 2007. - Т. 132. - С. 406-416.

97. Суслов И.М., Погорелов A.A. Оценка критических индексов из теоретико-полевой ренормгруппы: математический смысл «стандартных значений». ЖЭТФ. - 2008. -Т. 133. - С. 1277-1289.

98. Суслов И.М. Ренормгрупповые функции теории с/?4 в пределе сильной связи: аналитические результаты. ЖЭТФ. - 2008. - Т. 134. - С. 490-508.

99. Суслов И.М. Ренормгрупповые функции теории с/?4 из высокотемпературных разложений. ЖЭТФ. - 2011. - Т. 139. - С. 319-333.

100. Фольк Р., Головач Ю., Яворский Т. Критические показатели трехмерной слабо разбавленной замороженной модели Изинга. УФН. - 2003. - Т. 173. - С. 175-200.

101. Харди Г. Расходящиеся ряды. М.: Иностр. литература. 1951. - 240С.

102. Хмельницкий Д.Е. Фазовый переход второго рода в неоднородных телах. ЖЭТФ.- 1975. Т. 68. - Вып. 5. - С. 1960-1968.

103. Хоепберг П.С. Динамические явления в окрестности критической точки: жидкий гелий и антиферромагнетики. Квантовая теория поля и физика фазовых переходов.- М.: Мир. 1975. - С. 149-218.

104. Шалаев Б.Н. ФТТ. - 2010. - Т. 52. - С. 83.

105. Эллиот Р., Крамхансл Дж., Лис П. Теория и свойства неупорядоченных материалов.- М.: Мир. 1977. - 300С.

106. Юкалов В.И. Автомодельные приближения для проблемы собственных значений. -Дубна: ОИЯИ. 1991. - 15С.

107. Юхновский И.Р. Фазовые переходы второго рода. -Киев: Наук, думка. 1985. - 224С.

108. Aeppli G., Guggenheim H., Uemura Y.J. Spin dynamics near the magnetic percolation threshold. Phys. Rev. Lett. - 1984. - V. 52. - N. 11. - P. 942-945.

109. Aharony A. Critical phenomena in disordered systems. J. Magn. Magn. Mater. - 1978. -V. 7. -N. 1. -P.198-206.

110. Albano E.V., Bab M.A., Baglietto G., Borzi R.A., Grigera T.S., Loscar E.S., Rodriguez D.E., Rubio Puzzo M.L., Saracco G.P. Study of phase transitions from short-time non-equilibrium behaviour. Rep. Prog. Phys. - 2011. - V. 74. - P. 026501.

111. Alexandrowicz Z. Critically branched chains and percolation clusters. Phys. Lett. A. -1980. - V. 80. - N. 4. - P. 284-286.

112. Amit D. Field theory the renormalization group and critical phenomena. New York: Acad.press: McGraw-Hill. - 1978. - 333P.

113. Antonenko S.A., Sokolov A.I. Phase transitions in anisotropic superconducting and magnetic systems with vector order parameters: Three-loop renormalization-group analysis.- Phys.Rev. B. 1994. - V. 49. - N. 22. - P. 15901-15912.

114. Antonenko S.A., Sokolov A.I. Critical exponents for three-dimensional 0(n)-symmetric model with n > 3. Phys.Rev. B. - 1995. - V. 51. - N. 3. - P. 1894-1898.

115. Aplesnin S.S., Moskvin A.I Magnetic structures upon ordering of eg orbitals in a square lattice. Journal of Physics: Condensed Matter. - 2008. - V. 20. - P. 325202-325209.

116. Aplesnin S.S. Static and dynamic magnetic properties of coupled spin-1/2 antiferromagnetic chains. Journal of Phys.: Cond. Matt. - 2000. - V. 12. - N. 37. -P. 8191-207.

117. Baker G.A., Nickel B.G., Green M.S., Meiron D.I. Ising-model critical indices in three dimensions from the Callan-Symanzik equation. Phys. Rev. Lett. - 1976. - V. 36. - N. 23.- P.1351-1354.

118. Baker G.A., Nickel B.G., Meiron D.I. Critical indices from perturbation analysis of the Callan-Symanzik equation. Phys. Rev. B. - 1978. - V. 17. - N. 3. - P.1365-1374.

119. Ballesteros H.G., Fernández L.A., Martín-Mayor V., Sudupe M. New universality class in three dimensions?: the antiferromagnetic Mp2 model. Phys. Lett. B. - 1996. - V. 378. - P. 207-212.

120. Ballesteros H.G., Fernández L.A., Martín-Mayor V., Sudupe M. Finite size effects on measures of critical exponents in d = 3 O(N) models. Phys. Lett. B. - 1996. - V. 387. -P. 125-131;

121. Ballesteros H.G., Fernández L.A., Martín-Mayor V., Sudupe M. Critical properties of the antiferromagnetic Rp2 model in three dimensions. Nucl. Phys. B. - 1997. - V. 483. -P. 707-736.

122. Ballesteros H.G., Fernández L.A., Martín-Mayor V., Sudupe M. Critical exponents of the three-dimensional diluted Ising model. Phys. Rev. B. - 1998. - V. 58. - N 5. -P. 2740-2747.

123. Ballesteros H.G., Parisi G. Site-diluted three-dimensional Ising model with long-range correlated disorder. Phys. Rev. B. - 1999. - V. 60. - P. 12912-12917.

124. Bausch R. Dohm V., Janssen H.K., Zia R.K. Critical dynamics of an interface in 1+e dimensions. Phys. Rev. Lett. - 1981. - V. 47. - N. 25. - P. 1837-1840.

125. Bausch R., Janssen H.K., Wagner H. Renormalized field theory of critical dynamics. Z. Phys. B. - 1976. - V. 24. - P. 113-127.

126. Belanger D.P., King A.R., Jaccarino V. Crossover from random-exchange to random-field critical behavior in FexZnixF2 Phys. Rev. B. - 1986. - V. 34. - P. 452.

127. Belanger D.P., Birgeneau R.I., Shirane G., Yoshizawa H., King A.R., Jaccarino V. Critical dynamics of site-diluted three dimensional Ising magnet. J. de Physique Collque C8. -1988. - V. 49. - N. 7. - P. 1229-1238.

128. Belanger D.P., Young A.P. The random field Ising model. J. Magn. Magn. Mater. -1991. - V. 100. - N. 1-3. - P. 272-291.

129. Belanger D.P., Slanic Z., Fernandez-Baca J.A. Random-field critical scattering at high-magnetic concentration in the Ising antiferromagnet Fe0 93Zn0.o7F2. J. Magn. Magn. Mater. - 1998. - V. 177-181. - P. 171-172.

130. Benber C.M., Wu T.T. Anharmonic oscillator. Phys. Rev. - 1969. - V. 184. - P. 12311260.

131. Binder K., Reger J.D. Theory of orientational glasses. Models, concepts, simulations. -Adv. Phys. 1992. - V. 41. - P. 547-627.

132. Birgeneau R.I., Cowley R.A., Shirane G., Yoshizawa H., Belanger D.P., King A.R., Jaccarino V. Critical behaviour of site-diluted three dimensional Ising magnet. Phys. Rev. B. - 1983. - V. 27. - N. 12. - P. 6747-6757.

133. Blavats'ka V., Ferber C., Holovatch Yu. Entropy-induced separation of star polymers in porous media. Phys.Rev. E. - 2001. - V. 64. - P. 041102.

134. Boyanovsky D., Cardy J.L. Critical behavior of m-component magnets with correlated impurities. Phys. Rev. B. - 1982. - V. 26. - N. 1. - P.154-170.

135. Bresin E., Le Guillou J.C., Zinn-Justin J. Field theoretical approach to critical phenomena. Phase transition and critical phenomena. - ed. Domb C. and Lebowitz J.L.- New York: Acad, press. 1976. - V. 6. - P.127-249.

136. Calabrese P., Gambassi A. Aging in ferromagnetic systems at criticality near four dimensions. Phys Rev. E. -2002. - 65. - P. 066120.

137. Calabrese P., Martin-Mayor V., Pelissetto A., Vicari E. Three-dimensional randomly dilute Ising model: Monte Carlo results. Phys. Rev. E. - 2003. - V. 68. - P. 036136.

138. Calabrese P., Gambassi A. Ageing properties of critical systems. - J. Phys. A. - 2005.- V. 38. P. R133.

139. Cizek J., Vrscay E.R. Large order perturbation theory in the context of atomic and molecular physics-interdisciplinary aspects. Int. J. Quantum Chem. - 1982. - V. 21. -P. 27-68.

140. Chatelain C., Berche B. Universality and multifractal behaviour of spin-spin correlation functions in disordered Potts models. Nucl. Phys. B. - 2000. - V. 572. - N. 3. - P. 626-650.

141. Crisanti A., Ritort F. Violation of the fluctuation-dissipation theorem in glassy systems: basic notions and the numerical evidence. J. Phys. A. - 2003. - V. 36. - P. R181.

142. Chen K., Ferrenberg A.M., Landau D.P. Static critical behavior of three-dimensional classical Heisenberg models: A high-resolution Monte Carlo study. Phys. Rev. B. - 1993.- V. 48. P. 3249-3256.

143. Cooper F., Freedman B., Preston D. Solving <j)\ 2 field theory with Monte Carlo. Nucl. Phys. B. - 1989. - V. 210. - P. 210-228.

144. De Dominicis C., Brezin E., Zinn-Justin J. Field-theoretic techniques and critical. I. Ginzburg-Landau stochastic models without energy conservation. Phys. Rev. B. - 1975.- V. 12. N. 11. - P. 4945-4952.

145. De Dominicis C. Techniques de renormalisation de la théorie des champs et dynamique des phénomènes critiques. J. Physique (France). - 1976. - V. 37. - Suppl. 1. - P.l-247.

146. De Dominicis C. Peliti L. Field-theory renormalization and critical dynamics above Tc: Helium, antiferromagnets, and liquid-gas systems. Phys. Rev. B. - 1978. - V. 18. -P. 353-376.

147. Di Castro C. The multiplicative renormalization group and the critical behavior in d = 4 — e dimensions. Lett, nuovo cim. - 1972. - V. 5. - N. 1. - P. 69-74.

148. Di Castro C., Jona-Lasinio G. Renormalization group approach to critical phenomena. -Phase transition and critical phenomena. ed. Domb C. and Lebowitz J.L.- New York: Acad, press. - 1976. - V. 6. - P. 508-558.

149. Diehl H.W. in Phase Transitions and Critical Phenomena, edited by Domb C., Lebowitz J.L. London:Academic, 1986. - V. 10. - P. 75-267.

150. Diehl H.W. The theory of boundary critical phenomena. Intern. J. Modern Phys. B. -1997. - V. 11. - N 30. - P.3503-3523.

151. Diehl H.W. Critical behavior at m-axial Lifshitz points. Acta physica slovaca. - 2002.- V. 52. P. 271-283; E-print arXiv.org:cond-mat/0205284.

152. Diehl H.W., Gerwinski A., Rutkevich S. Boundary critical behavior at m-axial Lifshitz points for a boundary plane parallel to the modulation axes. Phys. Rev. B. - 2003. -V. 68. - P. 224428-1-224428-30.

153. Dorogovtsev S.N. The critical behaviour of systems with correlated defects. J. Phys. A.- 1984. V. 17. - P. L677-L679.

154. Dotsenko Vik.S., Harris A.B., Sherrington D., Stinchcombe R.B. Replica-symmetry breaking in the critical behaviour of the random ferromagnet. J. Phys. A. - 1995. -V .28. - P. 3093.

155. Dotsenko Vik.S., Feldman D.E. Replica symmetry breaking and the renormalization group theory of the weakly disordered ferromagnet. J, Phys. A. - 1995. - V. 28. -P. 5183.

156. Dotsenko V.S., Dotsenko V.S. Critical behaviour of the 2D-Ising model with impurity bonds. J. Phys. C. - 1982. - V. 15. - N. 3. - P. 495-507.

157. Dudka M., Holovatch Yu., Yavorskii T. A marginal dimension of a weakly diluted quenched m-vector model. J. Phys. Stud. - 2001. - V. 5. - N. 3. - P. 233-239.

158. Edwards S.F., Anderson P.W. Theory of spin glasses. J. Phys. F. - 1975. - V. 5. -P. 965.

159. Emery V.J. Critical properties of many-component systems. Phys. Rev. B. - 1975. -v.ll. - p.239

160. Ferrell R.A., Mirhashem B., Bhattacharjee J.K. Sound propagation in liquid helium near the A point. II. Ultrasonic attenuation. Phys. Rev. B. - 1987. - V. 25. - P. 4662-4668.

161. Feldman D.E., Izyumov A.V., Dotsenko Vik.S. Stability of the Renormalization Group in the 2D Random Ising and Baxter Models with respect to the Replica Symmetry Breaking.- e-print cond-mat/9512158. 1995.

162. Fisher M.E. The theory of equilibrium critical phenomena. Rep. Progr. Phys. - 1967. -V. 30. - P. 615-730.

163. Fisher M.E. Renormalization of critical exponent by hidden variables. Phys. Rev. -1968. - V. 176. - N. 1. - P. 257-272.

164. Fisher M.E. The renormalization group and the theory of critical behavior. Rev. Mod. Phys. - 1974. - V. 46. - N. 4. - P. 597-616.

165. Fisher M.E., Nelson D.R. Spin flop, supersolids, and bicritical and tetracritical points. -Phys. Rev. Lett. 1974. - V. 32. - N. 24. - P. 1350-1353.

166. Fishman S., Aharony A. Random field effects in disordered anisotropic antiferromagnets.- J. Phys. C. 1979. - V. 12. - N. 8. - P. L729-733.

167. Folk R., Holovatch Yu., Yavors'kii T. The correction-to-scaling exponent in dilute systems.- Pis'ma v ZETF. 1999. - V. 69. - N. 10. - P. 698-702.

168. Folk R., Iro H., Schwabl F. Critical dynamics and stactics of uniaxial dipolar magnets. -Z. Phys. B. V. 27. - P. 169.

169. Lubensky T.C., Rubin M.H. Critical phenomena in semi-infinite systems. II. Mean-field theory. Phys. Rev. B. - 1975. - V. 12. - P. 3885-3901.

170. Murtazaev A.K., Ramazanov M.K. Critical properties of the three-dimensional frustrated Ileisenberg model on a layered-triangular lattice with variable interplane exchange interaction. Phys. Rev. B. - 2007. - V. 76. - P. 174421-1-174421-6.

171. Freedman R., Mazenko G.F. Critical dynamics of antiferromagnets. Phys. Rev. B. -1976. - V. 13. - N. 12. - P. 4967-4983.

172. Ganton J.D., Kawasaki K. Renormalization group equations in critical dynamics. Progr. Theor. Phys. - 1976. - V. 56. - N. 1. - P. 61-76.

173. Gell-Mann M., Low F.E. Quantum electrodynamics of small distances. Phys. Rev. -1954. - V. 95. - N. 5. - P. 1300-1312.

174. Guida R., Zinn-Justin J. Critical exponents of the N-vector model. J. Phys. A: Math. Gen. - 1998. - V. 31. - P. 8103-8122.

175. Grassberger P. Damage spreading and critical exponents for "model A" Ising dynamics.- Physica A. 1995. - V. 214. - P. 547-559; Erratum. - Physica A. - 1995. - V. 217. -P. 227.

176. Griffiths R.B. Termodynamic function for fluids and ferromagnets near the critical point.- Phys. Rev. 1967. - V. 158. - N. 1. - P. 176-189.

177. Grinstein G., Fernandez J.F. Equilibration of random-field Ising systems. Phys. Rev.

178. B. 1984. - V. 29. - N. 12. - P.6389-6398.

179. Grinstein G., Luther A. Application of the renormalization group to phase transition in disordered systems. Phys. Rev. B. - 1976. - V. 13. - N. 3. - P. 1329-1343.

180. Grinstein G., Ma S.K., Mazenko G.F. Dynamics of spin interacting with quenched random impurities. Phys. Rev. B. - 1977. - V. 15. - N. 1. - P. 258-272.

181. Gropengiesser U. Damage spreading and critical exponent's for "model A" Isingtlynamics.- Physica A. 1995. - V. 215. - P. 308-310.

182. Halperin B.I., Hohenberg P.C. Calculation of dynamic critical properties. Phys. Rev. Lett. - 1967. - V. 19. - N. 2. - P. 700-703.

183. Halperin B.I., Hohenberg P.C., Ma S. Calculation of dynamic critical properties using Wilson's expansion methods. Phys. Rev. Lett. - 1972. - V. 29. - N. 23. - P. 1548-1551.

184. Halperin B.I., Hohenberg P.C., Ma S. Renormalization-group methods for critical dynamics. Phys. Rev. B. - 1974. - V. 10. - N. 1. - P. 139-153.

185. Halperin B.I., Hohenberg P.C., Siggia E.D., Ma S. Renormalization-group treatment of the critical dynamics of the binary-fluid and gas-liquid transition. Phys. Rev. B. - 1976.- V .13. N. 5. - P.2110-2123.

186. Hammersley J.M., Handscomb D.C. Monte Carlo Methods. London: Methuen & Co. -1964. 107 p.

187. Harris A.B. Effect of random defects on the critical behaviour of Ising models. J. Phys.

188. C. 1974. - V. 7. - N. 6. - P. 1671-1692.

189. Harris'A.B., Lubensky T.C. Renormalization-Group Approach to the Critical Behavior of Random-Spin Models. Phys. Rev. Lett. - 1974. - V. 33. - P. 1540.

190. Harris C.K., Stinchcombe R.B. Critical dynamics of diluted Ising systems. Phys. Rev. Lett. - 1986. - V. 56. - N. 8. - P. 869-872.

191. Hasenbusch M., Pelissetto A., Vicari E. J. Stat. Mech.: Theory Exp. - 2007. - P. P11009.

192. Heermann D.W., Burkitt A.N. Parallel Algorithms in Computational Science. Berlin: Springer-Verlag, 1991.

193. Hennecke M., Heyken U. Critical dynamics of cluster algorithms in the dilute Ising model.- J. Stat. Phys. 1993. - V. 72. - P. 829.

194. Heuer H.-O. Monte Carlo simulation of strongly disordered Ising ferromagnets. -Europhys. Lett. 1990. - V. 12. - N. 6. - P. 551-556.

195. Heuer H.-O. Monte Carlo simulation of strongly disordered Ising ferromagnets. Phys. Rev. B. - 1990. - V. 42. - N. 10. - P. 6476-6484.

196. Heuer H.-O. Monte Carlo simulation of disordered 2-dimensional Ising" systems. -Europhys. Lett. 1991. - V. 16. - N. 5. - P. 503-508.

197. Heuer H-O. Critical slowing down in local dynamics simulations. -vJ.Phys. A. 1992. -V. 25. - N. 9. - P.L567-L573.

198. Heuer H.-O. Critical crossover phenomena in disordered Ising systems. J. Phys. A. -1993. - V. 26. - N. 6. - P. L333-L339.

199. Heuer H.-O. Dynamic scaling of disordered Ising systems. J. Phys. A. - 1993. - V. 26.- N. 6. P. L341-L346.

200. Hohenberg P.C., Halperin B.I. Theory of dynamic critical phenomena. Rev. Mod. Phys.- 1977. V. 49. - P. 435-479.

201. Honkonen J., Komarova M.V., Nalimov M.Yu. Nucl. Phys. B. - 2005. - V. 714. - P. 292. e-print arXiv:hep-th/0406168.

202. Hornreich R.M. The Lifshitz point: Phase diagrams and critical behavior. J. Magn. Magn. Mater. - 1980. - V. 15-18. - P. 387-392.

203. Huse D.A. Remanent magnetization decay at the spin-glass critical point: A new dynamic critical exponent for nonequilibrium autocorrelations. Phys. Rev. B. - 1989. - V. 40. -P. 304-308.

204. Ikushima A., Feigelson R. Acoustic study of the critical phenomena in FeF2 near the Neel temperature. J. Phys. Chem. Solids. - 1971. - V. 32. - P. 417-425.

205. Imry Y. Tricritical points in compressible magnetic systems. Phys. Rev. Lett. - 1974.- V. 33. P. 1304.

206. Iro H., Schwabl F. Damping and dispersion of sound at structural transition. Solid State Communications. - 1983. - V. 46. - N. 2. - P. 205.

207. Ito N. Non-equilibrium relaxation and interface energy of the Ising model. Physica A.- 1993. V. 196. - P. 591-600.

208. Ivaneyko D., Ilnytskyi J., Berche B., Holovatch Yu. Criticality of the random-site Ising model: Metropolis, Swendsen-Wang and Wolff Monte Carlo algorithms. Condens. Matter Phys. - 2005. - V. 8. - P. 149.

209. Ivaneyko D., Ilnytskyi J., Berche B., Holovatch Yu. Local and cluster critical dynamics of the 3d random-site Ising model. Physica A. - 2006. - V. 370. - P. 163-178. e-print cond-mat/0603521.

210. Izyumov Y.A., Skryabin Y.N., Laptev V.M. Critical behaviour near the intersection of second-order phase transition lines in a random system. Phys. stat. solidi (b). - 1978. -V. 87. - N. 2. - P. 441-445.

211. Janssen H.K., Oerding K., Sengespeick E. On the crossover to universal criticality in dilute Ising systems. J. Phys. A. - 1995. - V. 28. - N. 21. - P. 6073-6085.

212. Jaster A., Mainville J., Schulke L., Zheng B. Short-time critical dynamics of the three-dimensional Ising model J. Phys. A: Math.Gen. -1999. - V. 32. - P. 1395.

213. Jayaprakash C., Katz H.J. Higher-order corrections to the varepsilon- expansions of the critical behaviour of the random Ising system. Phys. Rev. B. - 1977. - V. 16. - N. 9. -P. 3987-3990.

214. Jug G. Critical behaviour of disordered spin systems in two and three dimensions. Phys. Rev. B. - 1983. - V. 27. - N. 1. - P. 607-612.

215. Jug G. Critical singularities of the random two-dimensional Ising model. Phys. Rev. B.- 1983. V. 27. - N. 7. - P. 4518-4521.

216. Kadanoff L.P. Scaling laws for Izing models near Tc. Physics. - 1966. - V. 2. - N. 6. -P. 263-273.

217. Kalle C. Vectorised dynamics Monte Carlo renormalisation group for the Ising model -J. Phys. A. 1984. - V. 17. - P. L801.

218. Kawasaki K. Dynamical theory of fluctuations near critical points. Proceedings of the International school of physics Enrico Fermi course LI. ed. M.S.Green (Academic Press, New York and London. - 1971. - P. 342-379.

219. Kawasaki K., Ikushirna I. Velocity of sound in MnF2 near the Neel temperature. Physical Review B. - 1970. - V. 1. - N. 7. - P. 3143-3151.

220. Kleinert H., Neu J., Schulte-Frohlinde U., Chetyrkin K.G., Larin S.A. Five-loop renormalization group functions of 0(n)-symmetric </?4-theory and £-expansions of critical exponents up to e5. Phys. Lett. B. - 1991. -V. 272. - P. 39.

221. Korucheva E.R., De La Rubia F.J. Dynamical properties of the Landau-Ginzburg model with longe-range correlated quenched impurities. Phys. Rev. B. - 1998. - V. 58. - N. 9.- P. 5153-5156.

222. Korucheva E.R., Uzunov D.I. On the longe-range random critical behaviour. Phys. status solidi (b). - 1984. - V. 126. - P. K19-K22.

223. Korzhenevskii A.L., Luzhkov A.A., Heuer H.-O. Critical behaviour of systems with longerange correlated quenched defects. Europhys. Lett. - 1995. - V. 32. - P. 19-24.

224. Korzhenevskii A.L., Luzhkov A.A., Schirmacher W. Critical behavior of crystals with long-range correlations caused by point defects with degenerate internal degrees of freedom. Phys. Rev. B. - 1994. - V. 50. - N. 6. - P. 3661-3666.

225. Kosterlitz J.M., Nelson D.R., Fisher M.E. Bicritical and tetracritical points in anisotropic antiferromagnetic systems. Phys. Rev. B. - 1976. - V. 13. - N. 1. - P 412-433.

226. Kim J.K.,de Souza A.J., Landau D.P. Numerical Computation of Finite Size Scaling Functions: An Alternative Approach to Finite Size Scaling. Phys. Rev. E. - 1996. -V. 54. - P. 2291.

227. Kissner J.G. Nonequilibrium critical relaxation in the presence of random impurities. -Phys. Rev. B. -1992. V. 46. - P. 2676-2685.

228. Lawrie I.D., Prudnikov V.V. Static and dynamic properties of systems with extended defects: two-loop approximation. J. Phys. C. - 1984. - V. 17. - P. 1655-1668.

229. Le Guillou J.C., Zinn-Justin J. Critical exponents for the n-vector model in three dimensions from field theory. Phys. Rev. Lett. - 1977. - V. 39. - N. 2. - P. 95-98.

230. Le Guillou J.C., Zinn-Justin J. Critical exponents from field theory. Phys.Rev. B. -1980. - V. 21. - N. 7. - P. 3976-3998.

231. Leath P.L. Cluster size and boundary distribution near percolation threshold. Phys. Rev. B. - 1976. - V. 14. - P. 5046-5055.

232. Lubensky T.C. Critical properties of random-spin models from of the e expansion. -Phys.Rev. B. 1975. - V. 11. - N. 9. - P. 3573-3580.

233. Ma S-k., Mazenko G.F. Critical dynamics of ferromagnets in 6 —£ dimension. Phys.Rev. B. - 1975. - V. 11. - N. 11. - P. 4077-4100.

234. Mayer I.O. Critical exponents of the dilute Ising model from four-loop expansion. J. Phys. A. - 1989. - V. 22. - P. 2815-2823.

235. Mayer I.O. Sokolov A.I., Shalaev B.N. Critical exponents for cubic and impure uniaxial crystals: most accurate theoretical values. Ferroelectries. - 1989. - V. 95. - N. 1. -P. 93-96.

236. Metropolis N., Rosenbluth A.W., Rosenbluth M.N., Teller A.H., Teller E. Equation of State Calculations by Fast Computing Machines. Journal of Chemical Physics. - 1953.- V. 21. P. 1087-1092.

237. Mezard M., Parisi G., Virasoro M. Spin-Glass Theory and Beyond. Singapore: World Scientific. - 1987.

238. Mitchell P.W., Cowely R.A., Yoshizawa H., Boni P., Uemura Y.J. Critical behavior of the three-dimensional site-random Ising magnet: MnxZn!xF2. Phys. Rev. B. - 1986. -V. 34. - P. 4719.

239. Mukamel D. Tetracritical points in antiferromagnetic systems. Phys. Rev. B. - 1976. -V. 14. - N. 3. - P. 1303-1306.

240. Mukamel D., Grinstein G. Critical behavior of random systems. Phys. Rev. B. - 1981.- V. 25. N. 1. - P. 381-388.

241. Moran T.J., Luthi B. High-Frequency Sound Propagation near Magnetic Phase Transitions. Phys. Rev. B. - 1971. - V. 66. - P. 122-132.

242. Mori M., Tsuda Y. Vectorized Monte Carlo simulation of large Ising models near the critical point. Phys. Rev. B. - 1988. - V. 37. - P. 5444.

243. Nightingle M.P. Scaling theory and finite systems. Physica A. - 1975. - V. 83. - N. 3.- P. 561-572.

244. Nightingale M.P., Blote H.W. Monte Carlo computation of correlation times of independent relaxation modes at criticality. Phys. Rev. B. - 2000. - V. 62. - P. 1089.

245. Nelson D.R., Fisher M.E. Renormalization-group of metamagnetic tricritical behaviour.- Phys. Rev. B. 1975. - v. 11. - N. 3. - p. 1030-1039.

246. Nelson D.R., Kosterlitz J.M., Fisher M.E. Renormalization-group analysis of bicritical and tetracritical points. Phys. Rev. Lett. - 1974. - V. 33. - N 14. - P. 813-816.

247. Nelson D.R. Two-point correlations near four dimensions. Phys. Rev. B. - 1976. - V. 14.- P. 1123.

248. Newman K.E., Riedel E.K. Cubic N-vector model and randomly dilute Ising model in general dimensions. Phys. Rev. B. - 1982. - V. 25. - N. 1. - P. 264-280.

249. Oerding K. The dynamic critical exponent of dilute and pure Ising systems. J. Phys. A. - 1995. - V. 28. - P. L639-L643.

250. Ohta T., Kawasaki K. Mode coupling theory of dynamic critical phenomena for classical liquids. Progr. Theor. Phys. - 1976. - V. 55. - N. 5. - P. 1384-1395.

251. Onsager L. A two-dimensional model with an order-disorder transition. Phys. Rev. -1944. - V. 65. - N. 1. - P. 117-149.

252. Parisi G. The order parameter for spin glasses: a function on the interval 0-1. J. Phys. A. - 1980. - V. 13. - P. 1101.

253. Parisi G. A sequence of approximated solutions to the S-K model for spin glasses. J. Phys. A. - 1980. - V. 13. - P. L115.

254. Parisi G. Magnetic properties of spin glasses in a new mean field theory. J. Phys. A. -1980. - V. 13. - P. 1887.

255. Parisi G. Field-theoretic approach to second-order phase transitions in two- and three-dimensional systems. J. Stat. Phys. - 1980. - V. 23. - P. 49-82.

256. Parisi G., Ricci-Tersenghi F., Ruiz-Lorenzo J.J. Universality in the off-equilibrium critical dynamics of the three-dimensional diluted Ising model. Phys. Rev. E. - 1999. - V. 60. -P. 5198.

257. Pawlak A. Sound propagation in ammonium halides near the tricritical point. Phys. Rev. B. - 1991. - V. 44. - P. 5296.

258. Pawlak A., Fechner B. Sound attenuation and dispersion in a diluted Ising model. Phys. Rev. B. - 1989. - V. 40. - N. 13. - P. 9324.

259. Pearson R.B., Richardson J.L., Toussaint D. Dynamic correlations in the three-dimensional Ising model. Phys. Rev. B. - 1985. - V. 31. - N. 7. - P. 4472-4475.

260. Pelissetto A., Vicari E. Randomly dilute spin models: A six-loop field-theoretic study. -Phys. Rev. B. 2000. - V. 62 - P. 6393.

261. Pleimling M. Surface critical exponents at a uniaxial Lifshitz point. Phys. Rev. B. -2002. - V. 65. - P. 184406-1-184406-9.

262. Poole P.H., Jan. N. Dynamical properties of the two- and three-dimensional Ising models by "damage spreading" J. Phys. A. - 1990. - V. 23. - P. L453.

263. Preis T., Virnau P., Paul W., Schneider J.J. GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model. Journal of Computational Physics. - 2009. - V. 228 - P. 4468-4477.

264. Prudnikov P.V., Prudnikov V.V., Fedorenko A.A. Static and dynamic critical properties of 3D systems with long-range correlated quenched defects. J. Phys. A: Math.Gen. -1999. - V. 32. - N. 49. - P. 8587-8600.

265. Prudnikov P.V., Prudnikov V.V., Fedorenko A.A. Field-theory approach to critical behavior of systems with long-range correlated defects. Phys. Rev. B. - 2000. - V. 62. -N. 13. - P. 8777-8786.

266. Prudnikov P.V., Prudnikov V.V., Fedorenko A.A. Stability of critical behaviour of weakly disordered systems with respect to the replica symmetry breaking. Phys. Rev. B. - 2001. - V. 63. - N. 18. - P. 184201-184206.

267. Prudnikov P.V., Prudnikov V.V., Fedorenko A.A. Stability of critical behaviour of weakly disordered systems to introduction of potentials with replica symmetry breaking. J. Phys. A: Math.Gen. - 2001. - V. 34. - N. 12. - P. L145-L152.

268. Prudnikov P.V., Prudnikov V.V. Critical behaviour of weakly disordered systems with replica symmetry breaking potentials. J. Phys. Stud. - 2001. - V. 5. - N. 3/4. - P. 285-292.

269. Prudnikov P.V., Prudnikov V.V. The influence of disorder on the critical sound attenuation in solids. J. Phys.: Condens. Matter. - 2005. - V. 17. - P. L485-L492.

270. Prudnikov P.V., Prudnikov V.V., Dorofeev S.V., Kolesnikov V.Yu. Monte Carlo studies of critical behaviour of systems with long-range correlated disorder. Condensed Matter Physics. - 2005. - V. 8. - N. 1. - P. 213-224.

271. Prudnikov P.V. Free propagator for system at Lifshitz points: ordinary transition for perpendicular surface orientation. Materialien zum wissenschaftlichen Seminar der Stipendiaten des "Michail Lomonosov" . Bonn, 2005. - P.92-95.

272. Prudnikov P.V., Diehl H.W., Shpot M.A. Boundary critical behaviour at m-axial Lifshitz points of semi-infinite systems with a surface plane perpendicular to a modulation axis. -J. Phys. A: Math.Gen. 2006. - V. 39. - P. 7927-7942.

273. Prudnikov P.V., Prudnikov V.V. Critical sound attenuation of three-dimensional Ising systems. Condensed Matter Physics. - 2006. - V. 9. - N. 2. - P. 403-410.

274. Prudnikov P.V., Prudnikov V.V., Zheng B., Dorofeev S.V., Kolesnikov V.Yu. Short-time critical dynamics of the three-dimensional systems with long-range correlated disorder. -Progress of Theoretical Physics. 2007. - V. 117. - N. 6. - P. 973-991.

275. Prudnikov P.V., Prudnikov V.V. Influence of long-range correlated defects on critical ultrasound propagation in solids. Phys. Rev. B. - 2009. - V. 80. - P. 024115-1 - 02411511.

276. Prudnikov P.V., Prudnikov V.V., Krinitsyn A.S., Vakilov A.N., Pospelov E.A., Rychkov M.V. Short-time dynamics and critical behavior of the three-dimensional site-diluted Ising model. Physical Review E. - 2010. - V. 81. - P. 011130-1 - 011130-11.

277. Prudnikov V.V. On the critical dynamics of disordered spin systems with extended defects. J. Phys. C. - 1983. - V. 16. - N. 19. - P. 3685-3691.

278. Prudnikov V.V., Markov O.N. Monte Carlo renormalization group of dilute 2D Ising dynamics. Evrophys. Lett. - 1995. - V. 29. - P. 245.

279. Racz Z., Collins M.F. Linear and nonlinear critical slowing down in the kinetic Ising model: high-tempurature series. Phys. Rev. B. - 1976. - V. 13. - N. 11. - P. 3074-3077.

280. Raedt D.H., Lagendijk A. Monte Carlosimulations of quantum statistical lattice models.- Phys. Reports. 1985. - V. 127. - N. 4. - P. 233-307.

281. Rosov N., Kleinhammes A., Lidbjork P., Hohenemser C., Eibschutz M. Single-crystal Mossbauer measurement of the critical exponent /3 in the random-exchange Ising system Fe0 gZn0 iF2. Phys. Rev. B. - 1988. - V. 37. - P. 3265.

282. Rosov N., Hohenemser C., Eibschutz M. Dynamic critical behavior of the random-exchange Ising system Feo.9Zno.1F2 determined via Mossbauer spectroscopy. Phys. Rev. B. - 1992. - V. 46. - P. 3452.

283. Schorgg A.M., Schwabl F. Theory of ultrasonic attenuation at incommensurate phase transitions. Phys. Rev. B. - 1994. - V. 49. - P. 11682.

284. Shalaev B.N. Critical behavior of the two-dimensional Ising model with random bonds.- Phys. Rep. 1994. - V. 237. - N. 3. - P. 129-188.

285. Shalaev B.N., Antonenko S.A., Sokolov A.I. Five-loop y^-expansions for random Ising model and marginal spin dimensionality for cubic systems. Phys. Lett. A. - 1997. -V. 230. - P. 105.

286. Shapira Y. Experimental studies of bicritical points in 3D antiferromagnets. Multicritical phenomena. London-New York: Plenum press. - 1984. - P. 35-50.

287. Schehr G., Paul R. Universal aging properties at a disordered critical point. Phys. Rev. E. - 2005. - V. 72. - P. 016105-1-016105-7.

288. Schehr G., Paul R. Non-equilibrium critical dynamics in disordered ferromagnets. J. Phys: Conf. Series. - 2006. - V. 40. - P. 27. e-print arXiv:cond-mat/0511571.

289. Slanic Z., Belanger D.P. and Fernandez-В аса J.A. Equilibrium random-field Ising critical scattering in the antiferromagnet Fe(0.93)Zn(0.07)F2. Phys. Rev. Lett. - 1999. - V. 82.- P. 426.

290. Sokolov A.I., Varnashev K.B., Mudrov A.I. Critical exponents for the model with unique stable fixed point from three-loop RG expansions. Int. J. Mod. Phys. B. - 1998. - V. 12.- N. 12-13. P. 1365-1377.

291. Sokolov A.I., Varnashev K.B. Critical behavior of three-dimensional magnets with complicated ordering from three-loop renormalization -group expansions. Phys. Rev. B. - 1999. - V. 59. - P. 8363.

292. Stauffer D. Scaling theory of percolation clasters. Physics Reports. - 1979; - V. 54. -N. 1. - P. 1-78.

293. Stauffer D. Introduction to percolation theory. Taylor & Fransis. - 1985. - 294P.

294. Stinchcombe R.B. Dilute magnetism. Phase transitions and critical phenomena. ed. Domb C. and Lebowitz J.L. - New York: Acad, press. - 1983. - V. 7. - P. 151-191.

295. Stueckelberg E.C.G., Peterman A. La normalization des constantes dans la theorie des quanta. Helv. Phys. Acta. - 1951. - V. 25. - N. 5. - P. 499-520.

296. Suzuki M., Komatsubara T. Ultrasonic attenuation study on the critical dynamics of MnP near the Curie temperature. J. Phys. C. - 1982. - V. 15. - P. 4559.

297. Suzuki M. Quantum statistical monte carlo methods and applications to spin systems. J. Stat. Phys. 1986. - V. 43. - N. 5-6. - P. 883-909.

298. Swendsen R.H., Wang J.-S. Nonuniversal critical dynamics in Monte Carlo simulations.- Phys. Rev. Lett. 1987. — V. 58. — P. 86-88.

299. Talapov A.L., Shchur L.N. The critical region of the random-bond Ising model. J. Phys.: CM. - 1994. - V. 6. - P. 8295-8308.

300. Thurston T.R., Peter C.J., Birgeneau R.J., Horn P.M. Critical behaviour of site-diluted three dimensional Ising magnet. Phys. Rev. B. - 1988. - V. 37. - P. 9559-9563.

301. Tsypin M.M.-Effective potential for a scalar field in three dimensions: Ising model in the ferromagnetic phase. Phys. Rev. B. - 1997. - V. 55. - N. 14. - P. 8911.

302. Vasquez C.R., Paredes R.V., Hasmy A., Jullien R. New universality class for the three-dimensional XY model with correlated impurities: Application to 4He in aerogels.- Phys. Rev. Lett. 2003. - V. 90. - P. 170602.

303. Vincent E., Hammann J., Ocio M., Bouchaud J.P., Cugliandolo L.F. Slow dynamics and aging in spin glasses. Lect. Notes Phys. - 1997. - V.492. - P. 184-219.

304. Wang J.S., Selke W., Dotsenko VI.S., Andreichenko V.B. The two-dimensional random bond Ising model at criticality a Monte Carlo study. Europhys. Lett. - 1990. - V. 11. -N. 4. - P. 301-305.

305. Wang J.S., Selke W., Dotsenko Vl.S., Andreichenko V.B. The critical behaviour of the two-dimensional dilute Ising magnet. Physica A., 1990, v.164, p.221-239.

306. Wang J.S., Swendsen R.H. Monte Carlo and high-temperature-expansion calculations of a spin-glass effective hamiltonian. Phys. Rev. B. — 1988. — V. 38. - N 13. — P. 9086-9092.

307. Wang F., Hatane N., Suzuki M. Study on dynamical critical exponents of the Ising model using the damage spreading method. J. Phys. A. - 1995. - V. 28. - P. 4543-4552.

308. Wansleben S., Landau D.P. Monte Carlo investigation of critical dynamics in three-dimensional Ising model. Phys. Rev. B. - 1991. - V. 43. - P. 6006-6014.

309. Weinrib A., Halperin B.I. Critical phenomena in systems with long-range-correlated quenched disorder. Phys. Rev. B. - 1983. - V. 27. - P. 413-427.

310. Wegner F.J. Corrections to Scaling Laws. Phys. Rev. B. - 1972. - V. 5. - P. 4529.

311. Widom B. Equation of state in the neighbourhood of the critical point. J. Chem. Phys.- 1965. V. 43. - N. 11. - P. 3898-3916.

312. Wilson K.G. Feynmann-graph expansion for critical exponents. Phys. Rev. Lett. - 1972.- V. 28. N. 9. - P. 548-551.

313. Williams J.K. Monte Carlo estimate of the dynamical critical exponent of the 2D kinetic Ising model. J. Phys. A. - 1985. - V. 18. - P. 49.

314. Wilson K.G., Ficher M.E. Critical exponent in 3.99 dimensions. Phys. Rev. Lett. - 1972.- V. 28. N. 4. - P. 240-241.

315. Wiseman S., Domany E. Self-averaging, distribution of pseudo-critical temperatures and finite size scaling in critical disordered systems. Phys. Rev. Lett. - 1998. - V. 81. - P. 22; Phys. Rev. E. - 1998. - V. 58. - P. 2938.

316. Wolf U. Collective Monte Carlo updating for spin systems. Phys. Rev. Lett. - 1989. -V. 62. - P. 361-364.

317. Wolf U. Collective Monte Carlo updating in a high precision study of the x — y model. -Phys. Rev. Lett. 1989. - V. 62. - P. 361-364.

318. Yin J.Q., Zheng B., Prudnikov V.V., Trimper S. Short-time dynamics and critical behavior of three-dimensional bond-diluted Potts model. Eur. Phys. J. B. 2006. - V. 49.- P. 195.

319. Yoshizawa H., Belanger D.P. Mean-field simulation of field-induced domains and hysteretic behavior in dilute Ising antiferromagnets. Phys. Rev. B. - 1984. - V. 30.- N. 11. P. 5220-5228.

320. Yukalov V.I., Yukalova E.P. Spherical anharmonic oscillator in self-similar approximation.- J.Phys. A: Math.Gen. 1993. - V. 26. - P. 2011-2019.

321. Yukalov V.I., Gluzman S. Critical indices as limits of control functions. -Phys.Rev.Letters. 1997. - V. 79. - N. 3. - P. 4.

322. Yukalov V.I., Gluzman S. Algebraic self-similar renormalization in the theory of critical phenomena. Phys. Rev. E. - 1997. - V. 55. - N. 4. - P. 17.

323. Yukalov V.I., Yukalova E.P., Gluzman S. Extrapolation and interpolation of asymptotic series by self-similar approximants. arXiv.math-ph:1004.1041. - 2010. - P. 26.

324. Yukalov V.I., Yukalova E.P. Calculation of critical exponents by self-similar factor approximants. Eur. Phys. J. B. - 2007. - V. 55. -P. 93-99; e-print arXiv.cond-mat.stat-mech: 0704.2125. - 2007. - P. 17.

325. Zainullina R.I., Bebenin N.G., Burkhanov A.M., Ustinov V.V. Longitudinal sound velocity and internal friction in ferromagnetic La!xSrxMn03 single-crystal manganites model. Physical Review B. - 2002. - V. 66. - P. 64421.

326. Zheng B. Monte Carlo simulations of short-time critical dynamics. Int. J. Mod. Phys. B. - 1998. - V. 12. - P. 1419-1484.

327. Zinn-Justin J. Quantum field theory and critical phenomena. -- 1996. 1008 p.1. Oxford: Clarendon Press.