УДК 532.783
А. Н. Захлевных, Д. В. Макаров
ПЕРЕХОД ФРЕДЕРИКСА ПЕРВОГО РОДА В ФЕРРОНЕМАТИКАХ FIRST ORDER FREEDERICKSZ TRANSITION IN FERRONEMATICS
Пермский государственный университет,
614990 Пермь, ул. Букирева, д. 15
В работе исследуется влияние сегрегационных эффектов на индуцированные магнитным полем ориентационные фазовые переходы в ферронематике. Рассмотрена ориентация магнитного поля, ортогональная плоскости слоя. Найдено распределение директора, намагниченности и концентрации магнитных частиц при различных значениях напряженности магнитного поля, энергии сцепления магнитных частиц с нематиком, упругих констант Франка и сегрегационного параметра. Обнаружено, что переход Фредерикса в ферронематиках может быть фазовым переходом как первого, так и второго рода.
Ключевые слова: ферронематик, жидкий кристалл, фазовые переходы, три-критическая точка.
The influence of segregation effects on the magnetic field induced orientational phase transitions in a ferronematic is investigated. The magnetic field orientation orthogonal to layer surface are examined. Distributions of the director, the magnetization and the concentration of magnetic particles are found at various values of magnetic field strength, energy of cohesion of magnetic particles with a nematic, Frank elastic constants and segregation parameter. We found out, that Freedericksz transition in ferronematics could be the phase transition both of the first and the second order.
Key words: ferronematic, liquid crystal, phase transitions, tricritical point.
Введение
Как известно, жидкие кристаллы обладают анизотропией магнитных свойств, что позволяет менять ориентацию их оптической оси (директора) под влиянием внешнего магнитного поля. Этот эффект в ограниченных геометриях называется переходом Фредерикса. Требуемые для такого перехода магнитные поля оказывается возможным существенно уменьшить, если внедрить анизометричные магнитные частицы в жидкокристаллическую матрицу [1]. Такие суспензии магнитных частиц, приготовленные на основе нематического жидкого кристалла, получили название ферронематиков и интенсивно исследуются в последние годы как теоретически, так и экспериментально [2 - 28]. Их отличительной особенностью является наличие ориентационной связи (называемой сцеплением) между магнитными частицами и директором. Воздействие магнитного поля на ферронематик существенно зависит от типа сцепления, который характеризуется величиной энергии сцепления и взаимной ориентацией вытянутых осей
© Захлевных А. Н., Макаров Д. В., 2010
магнитных частиц и молекул жидкого кристалла. Ниже мы рассматриваем ферронематик с так называемым гомеотропным сцеплением между магнитными частицами и жидкокристаллической матрицей, при котором в отсутствие магнитного поля намагниченность ферронематика перпендикулярна его директору [2, 5, 14, 17, 20, 24, 28]. В этом случае магнитное поле, приложенное к ферронематику с положительной анизотропией диамагнитной восприимчивости жидкокристаллической матрицы, оказывает конкурирующее действие на его ориентацию: оно стремится ориентировать магнитные частицы и директор в одном направлении, чему препятствуют гомеотропные условия сцепления.
В данной работе исследуется индуцированный внешним магнитным полем переход Фредерикса в плоском слое ферронематика, директор которого абсолютно жестко закреплен в плоскости границ. Магнитное поле направлено поперек слоя. В работе показано, что при наличии достаточно сильных сегрегационных эффектов, заключающихся в индуцированном магнитным полем перераспределении магнитной примеси по толщине слоя, меняется характер перехода Фредерикса от второго рода (который этот переход имеет в «чистом» нематическом жидком кристалле и в ферронематике со слабой сегрегацией), к первому.
Уравнения равновесия
Рис. 1. Ориентация слоя ферронематика в магнитном поле Н
Рассмотрим слой ферронематика толщиной D, заключенный между двумя параллельными пластинами (рис. 1). Введем прямоугольную систему координат, ось x направим вдоль плоскости пластин, ось z - перпендикулярно пластинам; начало координат выберем в середине слоя. Сцепление директора n на границах слоя будем считать жестким и планарным (ось легкого ориентирования направлена вдоль x), а анизотропию диамагнитной восприимчивости %а положительной. Если к ферронематику приложено однородное магнитное поле H = H(cos рн, 0, sin рн), то директор n и единичный вектор намагниченности m можно искать в виде
n = (cosp(z), 0, sin p(z)), m = (- sin /(z), 0, cos/(z)). (1)
Свободная энергия ферронематика имеет вид Ф = { FdV, где [1, 8]
F = Fl + F2 + Fз + F4 + F5,
(2)
Кі п)2 + К 2 (п • rot п)2 + К з (п х rot п)2
1 2
F2 =-Msf тН, Fз =--Ха (пН )2,
2
F4 = кВТ/1п /, F5 = ^/(пт )2.
V d
Здесь К]_, К2, К3- модули ориентационной упругости нематического жидкого кристалла (константы Франка), п - директор ферронематика (единичный вектор, характеризующий направление преимущественной ориентации длинных осей молекул нематика), М <$ - намагниченность насыщения материала магнитных частиц, f - объемная доля магнитных частиц в суспензии, т - единичный вектор намагниченности суспензии, ха - анизотропия диамагнитной восприимчивости нематика (далее всюду предполагается, что ха > 0), V - объем феррочастицы, kв - постоянная Больцмана, Т - температура, н - поверхностная плотность энергии сцепления молекул нематического жидкого кристалла с поверхностью магнитных частиц, d - диаметр феррочастицы. Значение н выбирается положительным, так что в отсутствие магнитного поля минимуму энергии р соответствуют гомеотропные условия сцепления на частицах (т ± п).
Путем минимизации термодинамического потенциала по ф(г), щ(г) и f (г) при условии постоянства числа N магнитных частиц в суспензии | f dV = ^ приходим к следующей системе уравнений равновесия [27]:
1 ф0
2 = ±1 RFN (^ V)йф 2 о
1 фо ЛП
2 = ±1 gRFN(Ф,v)dФ, 2о
bhсоб(^ -фн ) = ° ^п 2(ф - V),
/ (г) Г Ь^ а . 2 Л
% = ^г = QехР1-----8тО-Фн)-------8т (ф^к
/ I к к і
где функция RF7 (ф, V) определена соотношением
=
h
2 2 СОБ ф + к БІП ф
22 Бгп (ф - Фн ) - 8гп (ф0 -Фн ) + 2к(%0 - %)
(3)
(4)
(5)
(6)
(7)
и введено обозначение % о = % (Фо, Vо) для концентрации магнитных частиц в середине слоя. Здесь / = 7 / V - средняя объемная доля магнитных частиц в суспензии.
Система уравнений (3) - (6) определяет зависимости углов ориентации директора фо = ф(о) и намагниченности Vо = V(0) в середине слоя ферронематика, а также
функцию распределения магнитных частиц / от напряженности внешнего магнитного поля h = ха / К , угла ориентации поля Ри, анизотропии ориентационной упругости k = Кз / К1, энергии сцепления магнитных частиц с жидкокристаллической матрицей с = н/ D2 /(К^) , параметра Ь = МD / */ К\Ха и различных значений сегре-
гационного параметра к = £ВТ / D2 /(К^).
х* ю-3
Рис. 2. Трикритическое значение сегрегационного параметра к* как функция энергии сцепления с при Ь = 0.1 (квадрупольный режим) для различной анизотропии упругости £
Рис. 3. Трикритическое значение сегрегационного параметра к* как функция энергии сцепления с при Ь = 10 (дипольный режим) для различной анизотропии упругости £
Поясним смысл безразмерного параметра Ь. Выше определена безразмерная напряженность магнитного поля И, где в качестве единицы ее измерения выбрана величина
Ич = Б_17К1/Ха . Она определена из условия, что при И « Иq энергия упругих деформаций р и диамагнитный (квадрупольный) р вклад в свободную энергию ферронематика (2) оказываются одного порядка. Аналогичное сопоставление упругого р и ферромагнитного (дипольного) р2 вкладов дает величину напряженности магнитного поля И d = К11(М3/02). Параметр Ь представляет собой отношение указанных выше полей Ь = Иq|Иd и потому характеризует механизм влияния магнитного поля на ФН [7]. При Ь >> 1 (Иd << И q ) искажение ориентационной структуры ФН в слабых полях осуществляется дипольным механизмом, а при Ь << 1 (Иq << Иd) - квад-
рупольным механизмом. Смена режима воздействия от дипольного к квад-рупольному (и наоборот) происходит в полях, для которых вклады р и Рз в свободную энергию становятся одного порядка, т. е. при И « И о = М / ха .
Заметим, что уравнение (6) описывает так называемый эффект сегрегации [1], заключающийся в росте концентрации магнитных частиц в тех местах образца, где минимальна сумма их магнитной энергии в поле Н и ориентационной энергии частиц в ЖК-матрице.
При ориентации магнитного поля под углом рн = ж/ 2 (рис. 1) и жестких планарных условиях сцепления директора с границами слоя система уравнений (3) - (6) допускает однородное решение, отвечающее ферронематику с гомеотропным (п ± т || Н ) сцеплением магнитных частиц с директором (р = щ = 0). При любых других ориентациях магнитного поля в рассматриваемой геометрии уравнения равновесия допускают лишь решение, для которого угол между п и т отличен от ж/2 (т. н. угловая фаза - см. [19, 24]). Далее мы рассматриваем случай рн = ж/2 .
Неоднородные решения для полей директора и намагниченности возникают
ЕМ
лишь в полях, выше порогового (поле Фредерикса hc ), которое в случае слабых сегрегационных эффектов определяется уравнением [6, 27]
'2 ' (8)
(а™ )2 = ж2 + ЬА™,
где введено обозначение
2а
2а + bh
ЕМ
(9)
В работе [27] показано, что вблизи поля Фредерикса hc равновесия имеет вид
ЕМ
Р0 =±1
у-
h - И
ЕМ
решение уравнений
(10)
к-к
где
8к
у=
4ем )
3Ь
5(2 - 5)ИЕМ f + 4ж2 шЕМ
> 0,
(11)
*
к=
(а™ )2
- ж
3Ь5 2(2 - 5)2 И™ + 4ж2 k
(12)
Формула (10) показывает, что переход из однородного состояния ферронематика в неоднородное (переход Фредерикса) обнаруживает трикритическое поведение: при
*
к > к (случай слабой сегрегации магнитной фазы) этот переход является переходом
*
второго рода, а при к < к (сильная сегрегация) - переходом первого рода. Значение
*
к = к , при котором меняется характер перехода, отвечает трикритической точке.
Если сегрегационные эффекты достаточно сильные и переход Фредерикса между однородным и неоднородным состояниями ферронематика осуществляется по типу
ЕМ
перехода первого рода, пороговое поле Ас теряет смысл поля перехода. В этом слу-
ЕМ
чае критическое поле А , при котором происходит равновесный переход Фредерикса,
*
2
находится из условия равенства полных свободных энергий Ф = | FdV однородного и
ЕМ
неоднородного состояний. При А = Ас углы ориентации директора и намагниченности испытывают конечные скачки, характерные для перехода первого рода.
Зависимость трикритического значения к* сегрегационного параметра от материальных параметров ферронематика показана на рис. 2 и 3. Как видно из рис. 2 и 3, характерные значения сегрегационного параметра к* в дипольном режиме (Ь = 10) в 103 раз больше, чем в квадрупольном (Ь = 0.1). Учет анизотропии ориентационной упругости показывает, что увеличение модуля продольного изгиба (увеличение k) приводит к уменьшению к* для обоих режимов (рис. 2 и 3, штриховые линии).
Полагая для ферронематика, приготовленного на основе жидкого кристалла 5СВ
[26, 29], при Т = 298 К: К^ = 0.64-10 -6дин, К 3 = 1.0-10 -6дин, %а = 1.7-10 -7,
М5 = 500 Гс, / = 2.0-10-7, d = 0.75-10-5см, = 10-3 ^ 10-1эрг/см2, у= 0,88-10-15 см3
и считая D = 2.5 -10-2 см, находим k « 1.6, а « 2.6 - (10-2 -И), Ь « 7.6 и к « 9.1-10-3 . В этом случае трикритическое значение сегрегационного параметра (12) меняется в до-
* -5 -1 -3
вольно широких пределах к «10 ^ 10 . Указанное выше значение к «9.1-10
лежит внутри этого интервала, поэтому в реальных ферронематиках можно ожидать переход Фредерикса по типу фазового перехода как первого, так и второго рода. Заметим, что задача в подобной геометрии была рассмотрена раньше в некоторых экспериментальных [9, 10] и теоретических работах [6], однако трикритическое поведение не было обнаружено. Это связано с тем, что в работах [9, 10] материальные параметры ис*
следуемых ферронематиков соответствовали к > к , а в работе [6] анализировалось только пороговое поле (8).
2
Заметим, что сегрегационный параметр может быть записан в виде к = ^ / Л) ,
где А = (к^/kвT/)1/2 - так называемая сегрегационная длина [1]. С учетом этого формула (10) принимает вид
Р0 =±_
о А - аЕМ
у Л2 2 с 2 , (13)
D2 - D2
где D* = Л4к представляет собой трикритическую толщину слоя. Для D < D* (тонкие слои) сегрегационные процессы существенны и переход Фредерикса осуществляется по типу перехода первого рода. Для указанных выше значений параметров ферронематика
находим D* « 2.5 -10-2см. Таким образом, трикритическое поведение ферронематика при переходе Фредерикса может быть достигнуто экспериментально в достаточно тонких слоях.
На рис. 4 и 5 показаны результаты численного решения уравнений ориентационного состояния ферронематика, которые определяют зависимости р(z), щ(z) и /(z)
*
при различных значениях напряженности поля. В случае слабой сегрегации (к > к ) в поле А « Ас распределение концентрации магнитных частиц по толщине слоя почти однородно (рис. 4). Увеличение напряженности поля вызывает миграцию магнитных частиц от середины слоя к его границам, т. е. эффект сегрегации - магнитные частицы накапливаются в тех областях слоя, где минимальна сумма их энергии в поле и ориен-
тационной энергии в нематической матрице. В случае сильной магнитной сегрегации
*
(к < к ) при И « А имеется существенная концентрационная стратификация магнитной примеси, так что подавляющая доля магнитных частиц локализована вблизи границ слоя (рис. 5).
Рис. 4. Ориентационное (слева) и концентрационное (справа) распределения в слое ферронематика для а = 10, & = 1.56, Ь = 10, к = 5 и различных значений поля И.
Здесь к « 2.07, Ис « 4.91
Рис. 5. Ориентационное (слева) и концентрационное (справа) распределения в слое ферронематика для а = 10, & = 1.56, Ь = 10, к = 0.1 и различных значений поля И.
Здесь к « 2.07, Ис « 4.91, А « 3.72
Заключение
В работе проанализировано влияние сегрегационных эффектов на индуцированный магнитным полем переход Фредерикса в ферронематике. Рассмотрен плоский слой ферронематика, ортогонально которому приложено однородное магнитное поле. На границах слоя были заданы условия жесткого планарного сцепления, на поверхности магнитных частиц сцепление предполагалось мягким и гомеотропным. Путем минимизации термодинамического потенциала ферронематика получены решения для плоских полей директора и единичного вектора намагниченности с учетом эффекта магнитной сегрегации. Произведен численный расчет углов поворота директора и намагниченности для различных значений напряженности приложенного магнитного поля и сегрегационного параметра.
Установлено, что переход Фредерикса в ферронематике может быть переходом как первого, так и второго рода в зависимости от степени сегрегации. При слабой магнитной сегрегации этот переход является переходом второго рода, как и переход Фредерикса для чистых жидких кристаллов. При сильной сегрегации переход Фредерикса в ферронематике становится переходом первого рода, обуславливая бистабильное поведение директора. Найдено аналитическое выражение для трикритического значения сегрегационного параметра, при котором происходит смена характера фазового перехода.
Работа выполнена при частичной поддержке гранта РФФИ № 10-02-96030.
Список литературы
1. BrochardF., de Gennes P. G. // J. Phys. (France). 1970. Vol. 31. P. 691 - 708.
2. Chen S. -H, Amer N. M. // Phys. Rev. Lett. 1983. Vol. 51. P. 2298 - 2301.
3. Figueiredo Neto A. M., SabaM. M. F. // Phys. Rev. A. 1986. Vol. 36. P. 3483 - 3485.
4. Райхер Ю. Л., Бурылов С. В., Захлевных А. Н. // Журн. эксперим. и теорет. физики. 1986. Т. 91. С. 542 - 551.
5. LiangB. J., Chen S.-H. // Phys. Rev. A. 1989. Vol. 39. P. 1441 - 1446.
6. Burylov S. V., Raikher Yu. L. // J. Magn. Magn. Mater. 1993. Vol. 122. P. 62 - 65.
7. Zakhlevnykh A. N., Sosnin P. A. // J. Magn. Magn. Mater. 1995. Vol. 146. P. 103 - 110.
8. Burylov S. V., Raikher Y. L. // Mol. Cryst. Liq. Cryst. 1995. Vol. 258. P. 107 - 122.
9. Koneracka M., Kellnerova V., Kopcansky P., Kuczynski T. // J. Magn. Magn. Mater. 1995. Vol. 140/144. P. 1455 - 1456.
10. Koneracka M., Zavisova V., Kopcansky P. et al. // J. Magn. Magn. Mater. 1996. Vol. 157/158. P. 589 - 590.
11. Berejnov V., Bacri J.-C., Cabuil V. et al. // Europhys. Lett. 1998. Vol. 41. P. 507 - 512.
12. Berejnov V., Cabuil V., Perzynski R., Raikher Yu. // J. Phys. Chem. B. 1998. Vol. 102. P. 7132 - 7138.
13. Potocova I., Koneracka M., Kopcansky P. et al. // J. Magn. Magn. Mater. 1999. Vol. 196/197. P. 578 - 580.
14. Kopcansky P., KonerackaM., Potocova I. et al. // Czech. J. Phys. 2001. Vol. 51. P. 59 - 63.
15. Burylov S. V., ZadorozhniiV. I., Pinkevich I. P. et al. // J. Magn. Magn. Mater. 2002. Vol. 252. P. 153 - 155.
16. Buluy O., Ouskova E., Reznikov Yu. et al. // J. Magn. Magn. Mater. 2002. Vol. 252. P. 159 - 161.
17. Kopcansky P., Potocova I., Koneracka M. et al. // Phys. Stat. Sol. (B). 2003. Vol. 236. P. 450 - 453.
18. Buluy O., Ouskova E., Reznikov Yu., Litvin P. // Ukr. J. Phys. 2004. Vol. 49. P. A48 -A50.
19. Zakhlevnykh A. N. // J. Magn. Magn. Mater. 2004. Vol. 269. P. 238 - 244.
20. Kopcansky P., Potocova I., Koneracka M. et al. // J. Magn. Magn. Mater. 2005. Vol. 289. P. 101 - 104.
21. Zadorozhnii V. I., Vasilev A. N., Reshetnyak V. Yu. et al. // Europhys. Lett. 2006. Vol. 73. P. 408 - 414.
22. Zakhlevnykh A. N., Makarov D. V. // Mol. Cryst. Liq. Cryst. 2007. Vol. 475. P. 233 -245.
23. Zadorozhnii V. I., Reshetnyak V. Yu., Kleshchonok A. V. et al. // Mol. Cryst. Liq. Cryst.
2007. Vol. 475. P. 221 - 231.
24. Tyszkiewicz C., Pustelny T., Nowinowski-Kruszelnicki E. // Eur. Phys. J. Special Topics.
2008. Vol. 154. P. 221 - 224.
25. Makarov D. V., Zakhlevnykh A. N. // J. Magn. Magn. Mater. 2008. Vol. 320. P. 1312 -1321.
26. Zadorozhnii V. I., Sluckin T. J., Reshetnyak V. Yu., Thomas K. S. // SIAM J. Appl. Math. 2008. Vol. 68. P. 1688 - 1716.
27. Захлевных А. Н., Макаров Д. В. // Вестник Пермского университета. Физика. 2009. Вып. 1(27). С. 62 - 68.
28. Tomasovicova N., Koneracka M., Kopcansky P. et al. // Acta Physica Polonica A. 2009. Vol. 115. P. 336 - 338.
29. Blinov L. M., Chigrinov V. G. Electrooptic Effects in Liquid Crystal Materials. Springer-Verlag, New York, 1994.
Поступила в редакцию 17.02.2010 г.