= ФИЗИКА КОНДЕНСИРОВАННОГО СОСТОЯНИЯ ВЕЩЕСТВА =
Электрофизические свойства монокристаллов метиламмониевых квасцов: роль сетки водородных связей в формировании температурных аномалий
Н.Д. Гаврилова,1, * И.А. Малышкина,1, * О.Д. Новик2, *
1 Московский государственный университет имени М.В.Ломоносова, физический факультет, кафедра физики полимеров и кристаллов Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2 2Московский государственный университет имени М.В.Ломоносова, химический факультет Россия, 119991, Москва, Ленинские горы, д. 1, стр. 3 (Поступила в редакцию 19.09.2022; после доработки 05.02.2023; принята к публикации 25.04.2023)
Пироэлектрические, диэлектрические свойства и токи термодеполяризации монокристаллов кристаллогидрата метиламмониевых квасцов были исследованы в интервале температур 120-200К. Температурные зависимости указанных величин имели максимумы при ^143, ^160 и 176.4К (точка Кюри). Поведение, общее для многих кристаллов с сеткой водородных связей с расстоянием О-Н...О ~2.7±0.1А, рассмотрено с точки зрения квантово-химической модели сетки водородной связи, содержащей, помимо отдельных молекул Н2О, кластеры воды. Перемещение мостиковых протонов в молекулярных кластерах в электрических полях является когерентным и скоррелированным, в него вовлечены все мостиковые протоны системы. Эти перемещения связаны с колебаниями кислородного каркаса типа сжатие-растяжение. Дискретность температурных максимумов электрофизических свойств связывается с дискретностью уровней энергии двухминимумного потенциала, вызванной расщеплением возбужденных уровней.
РЛСЯ: 77.22.Ej УДК: 537.226
Ключевые слова: дипольное упорядочение воды, квантово-химические модели кластеров водородных связей, пиротоки, токи термодеполяризации, дискретные уровни энергии, спутанные фотонные состояния.
БОТ: 10.55959/М8Ш579-9392.78.2340503
ВВЕДЕНИЕ
Алюмометиламмониевые сульфатные квасцы Ж3СН3А1(804)2 • 12Н20 (МАББ) являются классическим сегнетоэлектриком с фазовым переходом первого рода вблизи температуры Кюри Тс = 177К. Структура кристалла при комнатной температуре определена рентгеновским методом [1, 2]. Вокруг атома алюминия находится группировка Н2 О, а водородную связь образуют атомы кислорода этой группировки с кислоро-дами сульфатных групп длиной ¿оо ~ 2.64 А. Два других иона кислорода из образующегося тетраэдра участвуют в образовании слабых водородных связей кристаллизационной воды длиной ¿оо ~ 2.78 А. Структура состоит из положительных ([А1(Н20)6]3+, (СН3№Н3)+) и отрицательных ([В04]2-) ионов и шести молекул Н2О. Комплексный ион [А1(Н20)б]3+ является почти правильным октаэдром, симметрия группы [Б04]2- — С3щ, ось третьего порядка проходит через атом серы и один из атомов кислорода. Шесть ионов (СН3МН3)+ распределены вокруг пространственных диагоналей ячейки и каждый из них окружен шестью
* Е-таП: gavril@polly.phys.msu.ru ^ Е-таП: malysh@polly.phys.msu.ru ^ Е-таП: onovik99@gmail.com
молекулами Н20 на расстоянии ¿оо ~ 2.8 А, а на расстоянии 3.3 А от них находятся еще атомы О- групп [Б04]2-. Вокруг атома А1 находится также группировка Н2О, причем О—Н...О связи образуют связь с кислородами других молекул воды длиной ¿оо ~ 2.64 А. Два других иона кислорода тетраэдра сульфатной группы участвуют в слабых водородных связях с кислородами кристаллизационной воды длиной выше 2.78 А.
Рентгеновские исследования Флетчера [3, 4] показали, что наибольшие смещения при изменении симметрии ниже точки Кюри испытывают атомы метиламмониевых групп ^0.36 А, окружение которых также существенно меняется. Если в пара-электрической фазе метиламмониевый ион окружали 6 молекул воды почти октаэдрически со средним расстоянием ^3.55 А [1], то ниже 177К из них можно выделить две группы с расстояниями ^2.98 и 4.21 А соответственно до иона аммония и 4.19 и 3.24 А до метильной группы. При переходе изменяется и координация сульфатной группы относительно (СНзМНз)+ так, что каждая из двух групп (СНзМНз)+ примерно равно отстоит от двух сульфат-ионов. Сами сульфатные группы не меняют своей симметрии при фазовом переходе.
Для уточнения механизма фазового перехода Кришнаном с сотрудниками [5], О'Рейли [6], а также физиками школы Александрова [7] были проведены спектроскопические исследова-
ния методами ИК- и КР-спектроскопии и ЯМР на ядрах 'Ни
Изучение КР-спектров показало изменение симметрии иона (СН3МН3) + с Сз^ в парафазе до С' в сегнетофазе, а также появление в сегнетофазе новой полосы 3450 см-1, связанной с колебаниями растяжения группировок Н2О. ЯМР исследования на ядрах 'Н и 27А1 также показали, что наиболее сильное изменение при фазовом переходе присходит в группе СНзМНз. Изучение изоморфного замещения 804 на 8е04 [7] показало, что величина спонтанной поляризации в МАБеБ на ^40 К выше, чем в МАББ, а также существенно увеличены времена релаксации е. Это указывает на связь реориентации метиламмониевой группы с Б04 решеточными колебаниями.
Таким образом, главное изменение при фазовом переходе МАББ претерпевают Б04 и СН3МН3 группы, связанные Н2 0, а основной вклад в формирование спонтанной поляризации принадлежит ме-тиламмониевой группе. При этом возможно, что водородные связи не испытывают существенных изменений.
Исследования дейтерированных аналогов МАББ [8-10] показали, что температура Кюри и величина спонтанной поляризации изменяются незначительно, что будет для нас важно при обсуждении общей идеи работы.
Целью настоящей работы было связать дискретность наблюдаемых температурных аномалий электрофизических свойств (диэлектрической проницаемости е, пирокоэффициента 7 и токов термодеполяризации ]) с переходом системы на соответствующий возбужденный энергетический уровень двух-минимумного потенциала водородной связи. Для этого была привлечена теоретическая модель молекулярной динамики с использованием квантово-химических представлений о поведении мостико-вых протонов сетки водородной связи О-Н.. .О и их роль в поведении электрических параметров е, 7 и ] в зависимости от температуры.
1. ОБРАЗЦЫ И МЕТОДИКА ИЗМЕРЕНИИ
После ориентации по естественной огранке или рентгеновским снимкам кристаллы распиливались на плоскопараллельные пластины, которые затем полировались на увлажненном полотне. На рабочие поверхности образцов наносились электроды из сусального серебра или пасты Дегуса. Площадь образцов составляла 0.4-1 см2, толщина — 0.3-1.2 мм.
Пироэлектрический коэффициент 7 и токи термодеполяризации ] измеряли статическим методом. Изменение пироэлектрического заряда на поверхности образцов при изменении температуры наблюдали при помощи электрометра УА-^51. Ток разряда, появляющийся при изменении температуры, был зарегистрирован самопишущим потенциометром. Значения 7 определяли как изменение плотности поверхностного заряда с изменени-
ем температуры по сравнению с началом и концом ступенчатого нагрева. Величину заряда определяли через площадь под кривой ¿(t) на потенциограм-ме с учетом скорости движения ленты. Из-за того, что ток термодеполяризации j (t) на порядок медленнее пиротока ¿(t), мы успешно можем разделить эти токи, используя метод постепенного ступенчатого повышения температуры. Преимуществом данной техники является ступенчатый нагрев образца и его выдерживание в каждой точке в течении 10-20 мин [11], что позволило получить значения 7 и j, близкие к равновесным.
Диэлектрическую проницаемость е и диэлектрические потери tan S измеряли при помощи моста переменного тока на частоте 1 кГц в поле В/см.
2. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
Для объяснения механизма фазового перехода порядок-беспорядок в монокристаллах МАББ при 176.4 К мы использовали модель сетки водородных связей в воде, исходящей из идеи, основанной на присутствии в сетке О-Н... О двух видов формирований: среди множества свободных молекул Н2О находятся кластеры в виде цепей и колец (Н20)п (где п = 2, 3, 4, 6). Квантовая химия при изменении температуры учитывает, кроме колебания ядра водорода, взаимные колебания пар О-Н.. .О (рис. 1).
Рис. 1. Квантовая модель кластера О-Н.. .О
По литературным данным [12, 13], диэлектрическая проницаемость е МАББ может испытывать при Тс скачок как вверх (рис. 2), так и вниз. Кроме того, в работе Монтаньера [14] было отмечено наличие аномалий е в районе 150 К, а также влияние на них механического состояния образцов, постоянного электрического измерительного поля.
Температурная зависимость е(Т) в парафазе измерялась неоднократно [12, 13], она всегда удовле-
1 -1 t
165 175 185 Т, К
Рис. 2. Температурные зависимости диэлектрической проницаемости е пластин кристаллов МАЯБ различной толщины: 1.0 (1), 0.75 (2), 0.5 (3) и 0.35 мм (по данным [13])
творительно воспроизводится и демонстрирует стабильные значения. Закон Кюри-Вейсса выполняется в интервале 176.4^180.0 K с константой Кюри C -430 K и в интервале 180^188 K с C -540 K. При первом охлаждении е обычно скачком изменяется в Tc от 55 до 12 единиц и незначительно изменяется в сегнетофазе (рис. 3, кривая 2). При последующем нагревании наблюдается вновь скачкообразное изменение е(Т) до значений -50 с небольшим сдвигом Tc (—0.8 K). Если образец тут же нагреть выше Tc, то удается получить скачок вверх до е — 800 и при этом е быстро уменьшается при понижении температуры (рис. 3, кривая 1), причем от образца к образцу путь изменения е(Т) неоднозначен. В районе Torder — 145.3 K при этом получается воспроизводимая аномалия е без скачков и более устойчивая слева, чем справа от Torder. Если еще раз отогреть образец выше Tc, значение emax устанавливается медленнее, при стабилизированной температуре между 176.4 и 169 K время достижения равновесного значения е возрастает втрое и на зависимости 1/е при 169 K имеется излом. Ниже 145.3 K при всех изменениях мы наблюдаем четкое выполнение закона Кюри-Вейсса с константой Кюри C — 210 K и To —150 K. Отметим, что отношение констант Кюри выше 177 и ниже 145 K —2.05.
Традиционные измерения e(T) в постоянных электрических полях, параллельных (Ец ) и перпендикулярных (E^)Ps, подтвердили гипотезу о существенной доле ориентации Ps относительно ориентировки образца и переменного измерительного поля (EL). Результаты измерений приведены на рис. 4. Поле Ец незначительно повышает е в па-рафазе, сдвигает температуру максимума е вверх, причем AT/ДЕ — 4.0 х 10-5 град-см/В.
Согласно работе [15], при фазовом переходе первого рода скачок поляризации PS в точке Кюри связан со скачком энтропии Дй уравнением Клапейрона-Клазиуса dT/dE = PS /Д£.
Если Ps — 0.5 мкКл/см2, по нашим измерениям, а — 0.81 кал/моль-град [16], то дТ/дЕ = 4.2 х 10-5 град-см/В, что хорошо совпадает с полученным экспериментальным значением.
В случае наложения электрического поля перпендикулярно измерительным граням е также падает до значений 10^12.
Интересно влияние различных электродов на температурное поведение и величину е. Все измерения на рис. 4 выполнены на образцах с электродами из сусального серебра. Из литературы известно о выпрямляющем действии электродов некоторых видов на свойства сегнетоэлектриков [17]. Мы исследовали влияние на диэлектрические свойства электродов из эвтектической смеси индия с га-лием пасты на основе серебра и пасты Дегусса (рис. 5). Затвердевающие пасты, по-видимому, сжимают образец, что приводит к образованию 90-градусных доменов и резкому уменьшению величины е при охлаждении за точку Кюри. В парафазе паста Дегусса практически не влияет на величину е и незначительно сдвигает Тс вверх по температуре на —0.4К. В сегнетофазе наблюдается явное размытие фазового перехода. Серебряная паста снижает значение е как в пара-, так и в сегнетофазе, явно «размывая» переход. После снятия пасты на кристалле наблюдаются ямки травки. Эвтектическая смесь 1п-Са влияет на образец аналогично внешнему электрическому полю Ец —5 кВ/см. Таким образом, электроды из сусального серебра являются наиболее подходящими для точного исследования диэлектрических свойств, т.к. меньше других материалов изменяют механическое и электрическое состояние образца.
Особенно интересно проследить воздействие радиационных дефектов на диэлектрические свойства МАББ. Как и у всех водородосвязанных кристаллов, даже малые дозы радиации способны разрушить связи О-Н.. .О и М-Н.. .0, особенно если эти связи слабые, как это имеет место в МАББ. Из экспериментов по измерению е(Т) в 7-облученных кристаллах МАББ удалось установить, что наибольшим изменения подверглись области температур 145 К (рис. 6).
Влияние 7-облучения на е(Т) без поляризующего поля Е различно в различных температурных интервалах в области 140-180 К:
1. В парафазе Т > 176.4 К не обнаружено какой-либо закономерности для доз 1-10 кРн.
2. В области 176.4< Т <145.3 К с ростом дозы облучения увеличивается значение етах и общий уровень е в сегнетофазе.
3. В области Т <145.3 К при увеличении дозы от 0 до 30 кРн етах подавляется, размытия кривой не отмечено. Ниже Тог^ег сохраняется линейная зависимость 1/е(Т). Константы Кюри для доз 0, 15 и 30 кРн составили соответственно 224, 210 и 167 К.
Рис. 3. Температурная зависимость диэлектрической проницаемости е (а) и ее обратной величины (б) для монокристаллов МАЯБ
Рис. 4. Температурная зависимость диэлектрической проницаемости е монокристаллов МАЯБ при наложении внешнего электрического поля, параллельного спонтанной поляризации, величиной 0 (1), 5 (2), 10 (3) и 18 кВ/см (4)
Рис. 5. Температурная зависимость диэлектрической проницаемости е монокристаллов МАЯБ с электродами и различных материалов: эвтектической смеси индий-галлий (1), пасты «Дегусса» (2) и серебряной пасты (3)
В ряде случаев при облучении наблюдалась тенденция к увеличению температуры Кюри Тс при увеличении дозы, что в сегнетоэлектричестве связывают обычно с радиационным отжигом дефектов и снятием механического напряжения. Если аппроксимировать зависимость Тс от дозы облучения В прямой, то АТ/АБ = 4.7 • 10-3 град/кРн.
Другим чувствительным к ближнему порядку параметром сегнетоэлектриков является пироэлектрический коэффициент. Как известно, измерения пироэлектрических свойств сегнетоэлектрических монокристаллов возможны лишь при стабилизации монодоменного состояния каким-либо внешним воз-
действием (электрическим поляризующим полем, механическим полем для сегнетоэластиков, созданием 7-радиационных или других дефектов в поле поляризации или внутреннем поле кристалла, введением поляризующих дефектов на стадии роста и тд). В случае МАББ мы изучали пироэлектрические свойства на образцах без примесей при наложении электрических поляризующих полей или при 7-облучении кристаллов.
Типичные температурные зависимости пироко-эффициента 7 (Т) для различных поляризующих полей от 1 до 14 кВ/см представлены на рис. 7, а. На всех кривых имеются две области аномального поведения 7: в области точки сегнетоэлектри-
Рис. 6. Температурная зависимость диэлектрической проницаемости е монокристаллов МАЯБ вблизи 145 К, облученных 7-лучами дозами 1, 15 и 30 кРн
ческого фазового перехода ^177.4 K и в районе Torder ^143.5 K. По мере увеличения поля величина Y возрастает как в Tc, так и в Torder. В нашем случае при поле 8.2кВ/см ymax в Torder достигла своей максимальной величины. При дальнейшем увеличении E ymax в точке Кюри продолжает увеличиваться, а в Torder пик «размывается» и спадает.
Согласно данным [18], коэрцитивное поле в Tc кВ/см, что, возможно, определяет «порог» 8.2 кВ/см в нашем случае. Для наглядности методом графического интегрирования по изложенным данным были построены зависимости спонтанной поляризации Ps от температуры (рис. 7, б). В точке Кюри величина скачка составила 0.48 мкКл/см2, что в пределах ошибки измерений совпало с аналогичной величиной, полученной нами на тех же образцах по петлям диэлектрического гистерезиса. Изменение ступени APs при 143 K составило, по нашим данным, 0.05 мкКл/см2.
Сравнение данных Ps (T) в малых и больших (~14 кВ/см) полях показывает, что большое поле размывает ступеньку при 143 K и способствует появлению «хвоста» выше Tc, что типично для сегне-тоэлектрических переходов первого рода. Кривая Y(T) между 143 и 177 K ведет себя нестабильно, можно отметить еще две повторяющиеся области в районе ~158 и 165 K, где мы склонны относить аномалии y к перестройке доменной структуры.
Имеет свои особенности поведение пирокоэффи-циента y-облученных кристаллов MASD. На тех же образцах были измерены температурные зависимости пирокоэффициента и по ним методом графического интегрирования построена зависимость PS(T). При исследовании y-облученных образцов без наложения электрического поляризующего поля для малых доз облучения удавалось монодоме-низировать образец и измерить dPS/дТ. Облучение привело к значительному разбросу значений
Y и размытию структурного перехода и аномалии при 145K. На рис. 8 приведена также зависимость Ps(T) необлученного образца. Влияние облучения проявляется в уменьшении скачка Ps в Tc от 0.48 до 0.31 мкКл/см2 и появлении небольшого «хвоста» PS в парафазе. В раИоне 145 K наблюдается размытие «ступеньки» на кривой Ps(T) облученного кристалла.
Облучение кристаллов дозами 15, 30, 50 кРн при одновременном наложении поля 5 кВ/см при температуре кипения азота привело к следующим результатам. В парафазе уже не выполняется закон Кюри-Вейсса: зависимость 1/e(T) не является линейной. Этот факт указывает на изменение влияния Y-облучения на упорядочения движений ме-тиламинной группы, ответственной за фазовый переход. При дозах 15 и 30 кРн наблюдается скачок диэлектрической проницаемости всегда вниз, а при 50 кРн наблюдается скачок всегда вверх до ^800 кРн в необлученном образце. По-видимому, свыше 50 кРн величина внутреннего поля дефектов так велика, что закрепляется монодоменное состояние.
Аналогично облучению без поля в заполяризо-ванных образцах сдвиг Tc в сторону высоких температур составляет AT/AD ~ 3.2 • 10-3 град/кРн (рис. 8, а).
В области 140-150 K наблюдается увеличение Ymax с увеличением дозы до ~104 кРн и «размывание» аномалий при дальнейшем увеличении дозы (рис. 8, б).
Поскольку радиационные дефекты возникают при Y-облучении преимущественно за счет вторичных электронов, то можно предположить повреждение слабых связей в решетке, что типично для водородсодержащих сегнетоэлектриков [19]. Эта радиация должна в первую очередь обрывать самые слабые водородные связи, что обычно сказывается на таких характеристиках, как тангенс угла диэлектрических потерь tan 5, удельное сопротивление р.
На рис. 9 представлены температурные зависимости tan 5. Видно, что сопротивление образца очень мало в районе 270 и 150 K. Область высоких значений tan5 смещается вверх по температуре при облучении с увеличением дозы до кРн. На чистых образцах проводилось измерение удельного сопротивления от температуры p(T). Линейный ход p(T) выше 170 K и ниже 140 K хорошо согласуется с результатами для кристаллов диглициннитрата группы KDP.
Остановимся на экспериментальных данных по термодеполяризационному анализу. Все результаты по измерению токов термодеполяризации (ТДП) j(T) в температурном интервале 90-200 K представлены на рис. 10. Температурная зависимость токов ТДП в объектах с дискретным спектром электрически активных дефектов имеет вид ряда частично перекрывающихся пиков [20]. Каждому энергетическому уровню соответствует свой пик тока ТДП, который может дать величины энергии активации дефекта и его время релаксации.
Рис. 8. Температурная зависимость пирокоэффициента 7 и спонтанной поляризации Рв чистого (1) и облученного 7-лучами дозой 8 кРн (2) образца кристалла МАЯБ (а) и монокристаллов МАЯБ при различных дозах облучения в постоянных поляризующих полях (б)
Для МАББ зависимости токов ТДП ^(Т) снимались при различных условиях эксперимента (напряжение поля поляризации, температуры поляризации, скорости нагрева, дозе 7-облучения Б). Варьирование условий измерений позволяет при интерпретации экспериментальных данных разделить ди-польный и объемный зарядовый механизм релаксационного процесса, ответственного за каждый тип.
Как видно на рис. 10, а, для необлученного образца в полях 2.3 кВ/см2 можно выделить две области с максимумами ](Т): 175 и 145К; при увеличении поляризующего поля до 4.6 кВ/см2 можно устойчиво получить пик при —160 К. При облучении кристаллов дозами —2 и 6 кРн все три пика значительно усиливаются в малых полях — 0.5 кВ/см2 и уплощаются в полях выше 2-3 кВ/см2. При этом пик _7'(Т), отвечающий точке Кюри, четко прослеживается на всех без исключения кривых, в то время как максимум при 145 К сильно дефор-
мируется от внешних воздействий. По существу, вместо узкой области вблизи 145 К мы имеем широкую область выше —120 К, где происходит главное нарастание ^(Т).
Увеличение радиационных дефектов при разрыве водородных связей особенно сильно отразилось на температурном изменении токов деполяризации при малых дозах облучения. Значительных сдвигов областей аномалий (смещение Тс на десяток градусов), как это мы наблюдали в кристаллах триглицинсульфата (ТГС) и его изоморфов (триг-лицинтеллурата (ТГТел)), в МАББ не зафиксировано. Возможно, это связано с превалированием в МАББ сильных связей О-Н.. .О по сравнению с ТГС и ТГТел (см. статью Н.Д. Гавриловой и др. в этом журнале № 3. 46. (2020)).
Наличие особых температурных точек в кристаллах и полимерах с водородными связями можно объяснить при помощи моделей молекулярной ди-
У/-
г, к
Рис. 9. Температурная зависимость тангенса угла диэлектрических потерь tan S монокристаллов MASD при различных дозах Y-облучения: D = 0 (1), 1940 (2), 6040 Рн (3)
100 120 140 160 180
7, К
Рис. 10. Температурные зависимости токов термодеполяризации ] кристаллов МЛЯБ при различных напря-женностях поля Е и дозах 7-облучения Б
намики и квантово-химических моделей, заселением дискретных энергетических уровней системы при нагревании образца медленным ступенчатым образом. В настоящей работе принята модель воды, учитывающая когерентные коррелированные состояния системы (в том числе учитывающая наличие мостиковых протонов), что, вероятно, можно распространить на системы с сетками О—Н.. .О с ¿оо такой же, как у воды.
При наложении поляризующего поля перемещение протонов в водородных связях (мостиковых
протонов) в молекулярных кластерах воды является согласованным движением, в котором при наличии цепей и колец с упорядоченными водородными связями вовлекаются все мостиковые протоны. Колебания мостиковых протонов в молекулярных цепях и кольцах (в том числе входящих в состав объемных клеток) являются движениями больших амплитуд и связаны с колебаниями сжатие—растяжение кислородного каркаса молекул воды и зависит от расстояния между атомами кислорода в водородной связи ¿оо.
При несимметричном окружении того молекулярного кольца, в котором осциллируют мостико-вые протоны, конфигурация, отвечающая вершине потенциального барьера на пути изучаемого процесса, включает заряженные фрагменты И3О+5, И5О+5 и ОИ-5, число и взаимное расположение которых определены конфигурацией ближайшего окружения кольца.
Возможно, что в каждом температурном интервале превалирует определенный вид кластеров (цепи, димеры, оксимеры и т.д.). В области «пороговых» точек за счёт скачка энергии системы и протонов О—Н.. .О сетки при переходе на более высокий уровень двухминимумного потенциала меняются е, 7 и ], что проявляется в виде максимумов данных физических величин.
Перемещение мостиковых протонов в молекулярных кластерах сетки водородных связей является когерентным и скоррелированным, и в него вовлечены все мостиковые протоны системы. Колебания мостиковых протонов в молекулярных кластерах, входящих в объемные клетки, как и индивидуальных протонов, являются движениями большой амплитуды и связаны с колебаниями (сжатие—растяжение) кислородного каркаса структуры. Базовые частоты движения при значительном возбуждении, приводящем к существенным смещениям мостико-вых протонов, в 20—40 раз ниже частот их гармонических колебаний. При менее симметричном окружении кластера меньше заряженных фрагментов формируется в окружении. При менее симметричном окружении кластера—кольца меньше заряженных фрагментов Н3О5+ и ИбО^+ОН5- формируется в окружении, в котором перемещаются мостико-вые протоны на большие расстояния между ними и ниже потенциальный барьер двухямного потенциала. В небольших изученных кластерах, согласно анализу модельных одномерных квантовых задач, при энергиях возбуждения меньше 500 см-1 ниже вершины барьера формируются также состояния, для которых времена жизни составляют около пи-косекунды и меньше. И для самых мелких кластеров процесс оказывается возможным уже при энергии возбуждения 17—18 ккал/моль. Когда система приближается к верхним уровням энергии, сокращается время пребывания в этом состоянии.
В этом случае мы можем температурные точки аномалий 7 и ] в МЛББ (-176.4, 160, 140 и 120 К) связать с «пороговыми точками» сетки водородных связей О—Н.. .О с ¿оо -2.8 Л. Эти температуры
также подтверждены спектроскопическими данными (методами ЯМР, ИК и дс) [17, 19].
Малые дозы 7-радиации (особенно в районе 146 К и в точке Кюри) в МАББ, как и в других кристаллогидратах, способны разрушить водородные связи О-Н.. .О, О-Н.. N Измерения е(Т), 7(Т) и ^(Т) позволили установить, что облучение сильнее всего сказалось на области температур —145 К, а в пара-фазе вблизи 176.4 К не обнаружены заметные изменения Тс при дозах от 1 до 10 кРн.
При радиационном облучении образцов происходит разрыв слабых водородных связей, вследствие чего возникает много свободных протонов, оксонов и других кислородно-водородных фрагментов, что повышает прыжковую проводимость носителей тока и поэтому возрастает диэлектрическая проницаемость, пироток и ток термодеполяризации в пороговых температурных точках (для МАББ: 120, 143, 160, 177К).
В таких точках эти соединения и другие водо-родосодержащие диэлектрики испытывают аномалии электрофизических свойств: например, у фто-робериллата аммония ТС=177К, у изоморфа КБР С8Н2Ав04 — Тс=147 К и др.
По-видимому, для кристаллов с сеткой водородных связей длиной ¿оо —2.70 А между кислоро-дами характерны уровни энергии в двухминимум-ном потенциале, отвечающие этим «пороговым» температурам. Эти уровни образовываются вследствие расщепления возбужденного уровня основного состояния системы в поле дипольной матрицы других анионов.
Замечательной особенностью перемещения протонов в воде является их необычайно высокая подвижность, в 2 раза превосходящая подвижность ионов щелочных металлов, что связывают с согласованным (когерентным) перемещением мости-ковых протонов по цепочке водородных связей, при котором каждый из протонов перемещается на небольшие расстояния по сравнению с возможными перемещениями зарядов. При низких температурах (—100-150 К) такие явления связывают с увеличением доли туннелирования, без учета которого мала вероятность достижения высоких значений токов проводимости и деполяризации.
Ярким проявлением когерентных скоррелирован-ных колебаний пар молекул воды в цепочках или кольцах О-Н.. .О с углом —106° (симметрия гексагональных звездочек) проявляется в снежинках или призмах-столбиках, которые образуются в облаках в атмосфере в четких температурных интервалах: разнообразные звездочки от -40 °С, плоские шестигранники от -12 до -8° С, палочки-столбики от -8 до -5°С. При более низких температурах с неба сыплется снежная крупа или образуется колючий морозный туман.
Перечисленные выше температуры соответствуют нижним возбужденным энергетическим уровням двухминимумного потенциала или расщеплениям этих уровней во внешних полях. В наших опытах другие внутренние поля, которые соответству-
ют другим уровням энергии и температуры перестройки системы.
На основе развитой квантово-химической модели показано, что температурное поведение моды может быть описано откликом ансамбля классических ротаторов, взаимодействующим электроди-польным образом [21, 22]. Дипольное упорядочение при компьютерном моделировании, предсказанное в [23], подтвердилось при сегнетоэлектрическом переходе в кордиерите типа порядок-беспорядок при гелиевых температурах (—3К). Аналогичные макроскопические температурные зависимости пирото-ка и диэлектрической проницаемости мы наблюдали при 11 К в кристаллогидрате лития-таллий тарт-рате [24] в подрешетке О-Н...О водородных свзя-ей кристаллизационной воды, что подтверждалось отсутствием петель гистерезиса и отсутствием признаков структурного фазового перехода, что говорит о том, что природа аномалии е, 7 и ] связана только с переходом порядок-беспорядок в сетке О-Н. . .О-связей.
Возможно, что подавление структурного перехода в пороговых точках при температурах ниже 160 К в МАББ связано с присутствием поляризующих электрических полей. В работах [22, 23] выяснено, что с построением двухминимумного модельного потенциала возможно оценить вероятности и времена туннельного перемещения протонов для энергетического спектра состояний, а также что в случае высоколежащих (вблизи вершины потенциального барьера) дискретных состояний возможны проявления структур, сходных с переходными конфигурациями с временем существования —10-12 с. Таким образом, реализуется решетка взаимодействующих электрических диполей.
Особенности согласованных движений ядер кислорода и водорода в кластерах водородной сетки в целом коррелируют с известными коллективными свойствами систем, стабилизированных водородными связями.
Изучение свойств систем с включениями воды в каналы, мембраны, клетки или между слоями матрицы, а также в мембранную оболочку гидрати-рованных белков имеет широкое распространение, что важно, в том числе для понимания функционирования живых организмов. Феномен диполь-ного электрического упорядочения воды, помимо фундаментальных аспектов, может найти применение в исследовании биосовместимых устройств. Когда в системе велика роль дипольных взаимодействий (например, в сегнетоэлектриках), все частицы, «размазанные по пространству» — формулировка Фейнмана — вместе рассматриваются как область сильных корреляций в квантовой системе.
ЗАКЛЮЧЕНИЕ
В работе представлены результаты изучения температурных зависимостей диэлектрической проницаемости е, пирокоэффициента 7 и токов термо-
деполяризации у в широком диапазоне температур от 90 до 180 К для монокристаллов сегнето-электрика МЛББ. Экспериментальные зависимости демонстрируют ряд максимумов вблизи дискретных температур: -176.4 (ТС), 160, 143, 120К. Особенностью этих аномалий является совпадение с «пороговыми» точками воды, асимметричность А-зависимости и влияния замены водорода на дейтерий для ТС, что часто является признаком квантовых взаимодействий в системе, имеющей сетку О—Н.. .О связей. В квантовой химии это связывают с наличием дискретных уровней энергии системы и свободных радикалов р+, ОН-, И3О+ в двухмини-мумном потенциале, а также состояние структуры
матрица-сетка О—Н.. .О-связей и структур кластеров свободных молекул воды.
Полученные в работе результаты способствуют углублению знаний о свойствах нанолокализован-ных молекул воды, что важно для понимания свойств биологических объектов (анализе процессов в мембранной оболочке белков и общем метаболизме растительных организмов) и развития представлений об окружающей среде. Кроме того, в работе продемонстрировано, что взаимодействующие нанолокализованные полярные молекулы воды проявляют квантовые свойства, что представляет интерес с точки зрения создания перепутанных состояний в сегнетоэлектриках [21, 22, 25, 26].
[1] Lipson H.S. // Philosophical Magazine Series. 1, N 19. 887. (1935).
[2] Okaya Y., Ahmed M.S., Pepinsky R., Vand V. // Zeitschrift für Kristallographie - Crystalline Materials. 109, N 1-6. 367. (1957).
[3] Fletcher R.O.W., Steeple H. // Acta Crystallographica. 17, N 3. 290. (1964).
[4] Fletcher R.O.W., Steeple H. // Acta Crystallographica. 14, N 8. 891. (1961).
[5] Krishnan R.S., Narayanan R.S., Venkatesh G.M., Shmidt P. // J. Phys. Soc. Jap. 28. Supplement. 233. (1970).
[6] O'Reilly D.E., Tsang T. // Phys. Rev. 157, N 2. 417. (1967).
[7] Александрова И.П., Серюков И.В., Жеребцова М.И. // ФТТ. 27, N 11. 3438. (1985).
[8] Зайцева М.П., Жеребцова Л.И., Кожин В.И. et al. // Изв. АН СССР, сер. физ. 31, N 7. 1168. (1967).
[9] Зайцева М.П., Анистратов А.Т., Круминь И.М. // Кристаллография. 19, N 5. 1020. (1974).
[10] Зайцева М.П., Жеребцова Л.И., Виноградова И.С. // Изв. АН СССР, сер.физ. 29, N 6. 914. (1965).
[11] Gavrilova N.D., Novik V.K., Malyshkina I.A. // J. Non-Cryst. Solids. 483. 60. (2018).
[12] Mansingh A., Krupanidis S.B., Prasad E. // J. Chem. Phys. 69. 3039. (1978).
[13] Магатаев В.Н., Глушков В.Ф., Шувалов А.А. et al. // Кристаллография. 28, N 6. 1214. (1983).
[14] Le Montagner S., Rousselot M.-M. // J. Phys.
Radium. 21, N 10. 756. (1960).
[15] Струков Б.А., Сонин А.С. Введение в сегнетоэлек-тричество. М., 1970.
[16] Makita Y. //J. Phys. Soc. Jap. 19. 576. (1964).
[17] Lines M., Glass A. Principles and Applications of Ferroelectrics and Related Materials. Oxford, 2001.
[18] Иона Ф., Ширане Д. Сегнетоэлектрические кристаллы. М., 1965.
[19] Пешиков Е.В. Действие радиации на сегнетоэлек-трики. Ташкент, 1972.
[20] Верхопрахов В.Н., Сальман Е.Г. Термостимулиро-ванные токи в неорганических веществах. М., 1979.
[21] Бедняков А.С., НоваковскаяЮ.В. // Журнал Физ. Химии. 90, N 9. 1358. (2016). (BednyakovA.S., NovakovskayaY.V. // Russ. J. Phys. Chem. A. 90, N 9. 1813. (2016).)
[22] Бедняков А.С., Степанов Н.Ф., Новаковская Ю.В. // Журнал Физ. Химии. 88, N 2. 297. (2014). (Bednyakov A.S., Stepanov N.F. Novakovskaya Y.V. // Russ. J. Phys. Chem. A. 88, N 2. 287. (2014).)
[23] Belyanchikov M.A., Savinov M., Bedran Z.V. et al. // Nature Comm. 11. 3937. (2020).
[24] Gavrilova N.D., Malyshkina I.A., Novik O.D. // Ferroelectrics. 582, N 1. 1. (2021).
[25] Gavrilova N.D., Malyshkina I.A., Novik O.D. // Ferroelectrics. 585. N 1. 40. (2021).
[26] Novik O.D., Malyshkina I.A., Gavrilova N.D. // Ferroelectrics. 600, N 1. 124. (2022).
Electrophysical properties of methylammonium alum single crystals: the role of the hydrogen bond network in the formation of temperature anomalies
N.D. Gavrilova10, I.A. Malyshkina1'6, O.D. Novik2'0
1 Department of polymers and crystals physics, Faculty of Physics, Lomonosov Moscow State University. Moscow
119991, Russia
2 Faculty of Chemistry, Lomonosov Moscow State University. Moscow 119991, Russia E-mail: agavril@polly.phys.msu.ru, bmalysh@polly.phys.msu.ru, conovik99@gmail.com
Pyroelectric, dielectric properties, and thermodepolarization currents of methylammonium alum crystallohydrate single crystals were studied in temperature range 120-200 K. Temperature dependences of these quantities demonstrated peaks at —143, —160, and 176.4 K (Curie point). The behavior, which is typical to many crystals with a network of hydrogen bonds with an O-H...O distance of —2.7 ± 0.1 A, is
considered from the point of view of a quantum chemical model of a hydrogen bond network containing water clusters in addition to individual H2O molecules. The movement of bridging protons in molecular clusters in electric fields is coherent and correlated, and all bridging protons of the system are involved in it. These displacements are associated with oscillations of the oxygen frame of the compression-tension type. The discreteness of the temperature peaks of the electrophysical properties is associated with the discreteness of the energy levels of the double-well potential, caused by the splitting of the excited levels.
PACS: 77.22.Ej.
Keywords: dipole ordering of water, quantum chemical models of hydrogen bond clusters, pyroelectric currents, thermodepolarization currents, discrete energy levels. Received 19 .September 2022.
English version: Moscow University Physics Bulletin. 2023. 78, No. . Pp. 531-540. Сведения об авторах
1. Гаврилова Надежда Дмитриевна — доктор физ.-мат. наук, профессор, инженер; тел.: (495) 932-22-07, e-mail: gavril@polly.phys.msu.ru.
2. Малышкина Инна Александровна — канд. физ.-мат. наук, доцент, ст. науч. сотрудник; тел.: (495) 939-01-75, e-mail: malysh@polly.phys.msu.ru.
3. Новик Ольга Дмитриевна — студент; e-mail: onovik99@gmail.com.