Формирование пространственных неоднородностей в распределении изоморфных компонентов при кристаллизации твердых растворов

Штукенберг, Александр Григорьевич

Бесплатный автореферат и диссертация по наукам о земле на тему
Формирование пространственных неоднородностей в распределении изоморфных компонентов при кристаллизации твердых растворов
ВАК РФ 25.00.05, Минералогия, кристаллография

Автореферат диссертации по теме "Формирование пространственных неоднородностей в распределении изоморфных компонентов при кристаллизации твердых растворов"

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

003467296

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

Штукенберг Александр Григорьевич

ДЕ-

ФОРМИРОВАНИЕ ПРОСТРАНСТВЕННЫХ НЕОДНОРОДНОСТЕЙ В РАСПРЕДЕЛЕНИИ ИЗОМОРФНЫХ КОМПОНЕНТОВ ПРИ КРИСТАЛЛИЗАЦИИ ТВЕРДЫХ РАСТВОРОВ

Специальность 25.00.05 - минералогия, кристаллография

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

2 3АПР22-9

Санкт-Петербург - 2009

003467296

Работа выполнена на кафедре кристаллографии геологического факультета Санкт-Петербургского государственного университета

Научный консультант:

доктор геолого-минералогических наук, профессор Юрий Олегович Пунин

Официальные оппоненты: чл.-корр. РАН доктор геолого-

минералогических наук, профессор Асхаб Магомедович Асхабов

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

доктор геолого-минерапогических наук Григорий Юрьевич Иванюк

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

Московский государственный университет имени М.В. Ломоносова

Защита состоится 16 апреля 2009 г. в 15 час. на заседании совета Д.212.232.25 по защите докторских и кандидатских диссертаций при Санкт-Петербургском государственном университете по адресу: 199034 г. Санкт-Петербург, Университетская наб., 7/9, геологический факультет, ауд. № 52.

E-mail: m-char@yandex.ru

С диссертацией можно ознакомиться в Научной библиотеке им. A.M. Горького Санкт-Петербургского государственного университета.

Автореферат разослан « Мв.рг^ 2009 г.»

Ученый секретарь диссертационного совета, доктор геол.-мин. наук

М.В. Чарыкова

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Научная новизна.

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

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

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

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

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

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

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

Практическая значимость. Созданная количественная модель взаимодействий в системах «твердый раствор - водный раствор» дает возможность управления составом кристалла и его пространственной однородностью. Результаты анализа конкретных типов пространственных неоднородностей распределения компонентов по кристаллу могут быть использованы для реконструкции условий минералообразования. Полученный материал используется в курсах лекций «Аномальная кристаллооптика», «Деформационные дефекты кристаллов», «Кристаллогенезис», читаемых на геологическом факультете СПбГУ.

Защищаемые положения.

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

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

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

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

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

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

Апробация работы и публикации. Материалы работы опубликованы в 3 монографиях, 33 статьях в реферируемых журналах, 58 тезисах докладов. Они были представлены на конференциях разного уровня, в том числе: Международное совещание по рентгенографии и кристаллохимии минералов (Рентгенография минерального сырья) (1992, 1999, 2003, 2007); Международная конференция «Спектроскопия, рентгенография и кристаллохимия минералов» (1997); International Conference of Crystal Growth (1998); European Crystallographic Meeting (1998); Congress and General Assembly of the International Union of Crystallography (1993,1996, 1999); Уральское кристаллографическое совещание (1998); Национальная конференция по росту кристаллов (2000, 2002, 2004, 2006); Международная конференция "Кристаллогенезис и минералогия" (2001,2007); Совещание по экспериментальной минералогии (1995,2001); International Symposium on Experimental Mineralogy, Petrology and Geochemistry (2002); General Meeting of the International Mineralogical Association (1994, 2002); Международный симпозиум «Минералогические музеи» (2002, 2005); Национальная кристаллохимическая конференция (2003, 2006); Международная конференция по физике кристаллов «Кристаллофизика 21-го века» (2006).

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

Благодарности. Прежде всего, выражаю глубокую благодарность моему научному консультанту Ю.О. Пунину. Хочу поблагодарить П.Я. Азимова, О.В. Франк-Каменецкую, И.В. Рождественскую, Г. Клаппера, А. Кирфеля, X. Ойлера, Д.Ю. Попова, М.В. Морозова за плодотворное сотрудничество; И.И. Баннову, О.И. Артамонову, E.H. Котельникову, А. Путниса, Х.М. Астиллероса, А.Э. Волошина за помощь при выполнении отдельных этапов работы; П.Б. Соколова, A.A. Золотарева, А.Н. Зайцева за предоставление материалов для исследования, а также студентов и аспирантов кафедр кристаллографии и минералогии геологического факультета СПбГУ за помощь в получении экспериментальных данных.

СОДЕРЖАНИЕ РАБОТЫ

Глава 1. Особенности кристаллизации в системах «твердый раствор - водный раствор»

Прежде чем рассматривать формирование различного рода пространственных неоднородностей в распределении изоморфных компонентов, необходимо исследовать особенности кристаллизации твердых растворов. В данной работе мы, главным образом, ограничиваемся кристаллами твердых растворов, образующимися из водных растворов, но принципиально результаты применимы к кристаллизации из любых растворов. Хотя теория кристаллизации простых солей из водных растворов разработана достаточно хорошо, процессы кристаллизации твердых растворов поняты гораздо хуже, и надежные данные в этой области ограничиваются возможностью моделировать фазовые равновесия (Lippmann, 1980; Glynn, 2000). Весьма небольшое количество работ посвящено кинетике кристаллизации, взаимодействиям твердого и неравновесного с ним водного раствора

(Гликин, 2004) и влиянию напряжений несоответствия на рост кристаллов (Во1к1кт1уапоу, 1982; Мошкин и др., 2000).

Нами рассматривается простая система (В,С)А - Н20, в которой катионы В,г+ и С"т и анион А"" присутствуют в водном растворе и образуют непрерывный ряд твердых растворов. Мы не учитываем процессы зародышеобразования и массопереноса и ограничиваем каше рассмотрение процессами на фронте роста кристалла. Рост кристалла представляет собой переход вещества водного раствора, (обозначаемого индексом "ч) в твердый раствор (обозначаемый индексом'"), и часто описывается как гетерогенная химическая реакция:

хвлВ"++хслС"++А"-=В1|мС,1лА, (1)

где правая часть уравнения относится к твердому раствору, а левая - к водному; а хИЛ и*са представляют собой мольные доли компонентов В А и СА в твердом растворе соответственно. С точки зрения термодинамики такая запись реакции не корректна (Пригожин, Дефей, 1954), и более правильно рассматривать две самостоятельные реакции, связанные только условием баланса массы х!1Л + хСл =

В^'+Л"" =ВА" (2)

С"" + А"1' = СА" (3)

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

ап е^м, = = /ВЛ11!Лх,!Л (4)

ас«,аич = = /сл к-лх гл' <5)

где ¿,А - произведения растворимости конечных членов изоморфного ряда (/ = В, С),/1А - коэффициент активности /-го компонента в твердом

растворе, для регулярных твердых растворов j- _ еХр

-xj

- энтальпийный параметр смешения в твердом растворе.

1П = ал (йг„ + ас) - ионное произведение активности, у, - мольная доля i-

го компонента в солевой части раствора, а а, - активность /-го иона в водном растворе. Уравнения (4) и (5) позволяют рассчитывать фазовую диаграмму в общем виде (Lippmann, 1980; Glynn, Reardon, 1990; Glynn et

al, 1990; Glynn, 2000). Из них также сразу следует выражение для равновесного коэффициента распределения

К - Ус _ ^СЛ . /сл . (6)

4 Г ' V I f

СА .У В ЬВЛ -/ ВЛ

Нами получена удобная модификация этих уравнений для хорошо растворимых соединений и произведено распространение данного метода расчета фазовых диаграмм на общий случай сложного электролита (Bi,B2,B3...)¿A/Zm с несколькими типами ионов, где B¿ -изоморфные компоненты. Как показало сопоставление с экспериментальными данными по примерно 50 водно-солевым системам, хорошее соответствие с расчетами достигается при использовании дополнительного подгоночного параметра у/. Введение в расчеты коэффициентов активности в водном растворе позволяет достичь соответствия между экспериментом и теорией, не прибегая к подгоночному параметру (//(рис. 1а). Вместе с тем, созданный полуэмпирический метод расчета фазовых диаграмм оказывается удобным для описания имеющихся экспериментальных данных.

Если в целом ситуация с равновесием достаточно понятна, то существенно большие трудности и споры (Astilleros et al., 2003; Гликин, 2004; Prieto et al., 2007) возникают при описании взаимодействий между водным и твердым растворами произвольного состава. Рассмотрим равновесие твердого раствора фиксированного составах^. = 1 - хСЛх с серией водных растворов по каждому из компонентов в отдельности, согласно реакциям (2) и (3), и исследуем поведение частных (для каждого компонента) произведений активности как функций состава заводного раствора ЕП(С? = Llj^f¡s,.xiKJyi (i = В, С). На фазовой диаграмме

Липпмана эти функции представлены двумя гиперболами (рис. 16). Точка их пересечения (точка О на рис. 16) соответствует истинному равновесию = = ЕП (твердый раствор находится в

равновесии с обоими компонентами). Выше (ниже) соответствующих линий химические потенциалы компонентов /А = ВА, СА в водном растворе у, больше (меньше), чем химические потенциалы этих компонентов в данном твердом растворе хсл,. Соответственно, полученные кривые делят фазовое пространство на четыре области,

обозначенные на рис. 16 римскими цифрами. Относительно заданного твердого раствора водные растворы в области I пересыщены в отношении обоих компонентов ВА и СА (область роста). В области Ш водные растворы недосыщены в отношении обоих компонентов (область растворения). Две оставшиеся области соответствуют водным растворам, пересыщенным и отношении одного компонента (ВА в области II и СА в области IV) и недосыщенным в отношении другого (области обменных реакций). Данный подход позволяет теоретически описать взаимодействия между водным и твердым растворами, изучение которых ведется пока только на эмпирическом уровне (Гликин, 2004).

Рис. 1. а. Фазовая диаграммы систем с непрерывным рядом твёрдых растворов при 25°С. Система (NH4)2(Cu,Ni)(S04)2'6H20-H20; Точки -экспериментальные данные (Hill, Taylor, 1938), кривые - расчет: штрих-пунктирные линии рассчитаны без учета коэффициентов активности в водном растворе и без использования подгоночных параметров, пунктирные - без учета коэффициентов активности в водном растворе, но с использованием подгоночного параметра у/, сплошные - с учетом коэффициентов активности в водном растворе и без использования подгоночных параметров, б. Частные равновесия в системе «твердый раствор - водный раствор». Идеальный водный и твердый растворы, £вд = 4, LCA = 10. Сплошные линии показывают равновесия между твердым раствором jCAi = 0.4 и непрерывными сериями водных растворов по компонентам ВА и СА. Равновесная фазовая диаграмма обозначена штриховыми линиями.

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

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

Кристаллизация вдали от равновесия уже не может быть корректно описана в рамках термодинамического подхода, и необходимо использовать уравнения химической кинетики. Мы рассматриваем независимые реакции присоединения/отрыва компонентов ¿А (/' = В, С) в соответствии с уравнениями (2) - (3). Потоки компонентов в кристалл могут быть записаны как

•Л-Л = Кл ( - ¿сл/сл*сл ) = Кл (УсШ - £сл/сл-*сл ) ' (8)

где - кинетические коэффициенты. В случае стационарного процесса роста выражения (7) - (8) должны быть дополнены уравнением

Л,л _ £вд.. (9)

Результаты модельных расчетов показывают (рис. 2а, нижняя линия), что при увеличении пересыщения (скорости роста) коэффициент распределения отклоняется от равновесного значения, асимптотически стремясь к постоянной величине Кх - /¿С+Л . Важным выводом из

данного рассмотрения является возможность перехода коэффициента распределения через единицу при увеличении скорости роста кристалла (такое поведение наблюдалось в эксперименте и до сих пор не имело удовлетворительного объяснения). Количественная проверка полученных соотношений выполнена на примере системы К2(№,гпХ504)2-6Н20 - Н20, исследованной в (Жмурова, Хаимов-Мальков, 1970). Наблюдается неплохое соответствие между теорией и экспериментом (рис. 26).

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

Рис. 2. Зависимость коэффициента распределения К от пересыщения, а. Система, показанная на рис. 16. Около кривых подписана величина несоответствия решеток затравки и наросшего слоя ~bCA)/bt •> упругие

постоянные E¡{\- к) = Юи Н/м2, молярный объем w = 40 см3/моль (остальные

численные параметры указаны в подписи к рис. 16). Кинетические коэффициенты к*л = 10, к^л = 16 усл. ед., состав водного раствораус = 0.4. Штриховая линия соответствует величине коэффициента распределения при бесконечно большой скорости роста Кг, = к'рл ¡k¿.A =0.625- б. Система

K2(Ni,Zn)(S04)2-6H20 - Н20. Треугольники и кружки - экспериментальные данные (Жмурова, Хаимов-Мальков, 1970) для секторов роста {001} и {110}, соответственно, и состава солевой части водного раствора yN¡ = 0.01 - 0.137. Оттуда же взяты кинетические данные. Сплошные линии - расчет для >!Ni = 0.01, 0.05,0.1. Штриховая линия соответствует величине коэффициента распределения при бесконечно большой скорости роста = кщ jк-/п - 0.389 •

отношению коэффициентов диффузии или комбинации коэффициентов диффузии и кинетических коэффициентов соответственно.

Здесь также уместно обратиться к понятию пересыщения для системы «твердый раствор-водный раствор», поскольку по этому важному вопросу в литературе нет единого мнения (Kasatkin et al., 1995; Prieto et al., 1993,2007; Astilleros et al., 2003). Мы полагаем, что не существует универсальной величины, которая бы характеризовала пересыщение (или недосыщение) в общем случае, так как процесс управляется двумя независимыми разностями химических потенциалах компонентов А//ил и ÁfiCA. Обращаем внимание, что пересыщение не может быть определено только из состава водного раствора, как для однокомпонентных кристаллов, поскольку необходимо учитывать также

состав образующегося кристалла, который может быть определен из решения системы уравнений (7) - (9).

До сих пор мы предполагали, что после отложения нескольких новых слоев затравочный кристалл больше не влияет на процесс роста. Эксперименты, однако, показывают, что такое влияние может быть существенным за счет возникновения напряжений несоответствия (гетерометрии), хорошо известных для гетероэпитаксиальных пленок и практически не изученных для кристаллов, растущих из растворов. Смещение фазовых равновесий и изменение состава кристалла под действием напряжений несоответствия рассматривалось в ряде работ (Дорфман, 1979; Во1кЬоукуапоу, 1981, 1982; Болховитянов, 1990; Мошкин и др., 2000). В нашу задачу входило получить явное выражение для коэффициента распределения между водным раствором и напряженным кристаллом, исследовать зависимость коэффициента распределения от скорости роста кристалла и учесть влияние процессов релаксации напряжений при росте кристалла.

Величина упругой энергии ¿/тонкого слоя, отложенного на толстую подложку, может быть вычислена исходя из параметров элементарной ячейки подложки Ьг и наросшего слоя Ъ (Ольсен, Этгенберг, 1978). Эта энергия характеризует твердый раствор как целое и подобно другим средним мольным величинам может быть представлена как сумма двух парциальных величин, относящихся к компонентам ВА и СА: и = хикиа„ +х,.л(У„л. Вводя эти величины в

лЛ НА (-Л С А

выражения для химических потенциалов, получим новые выражения для

произведений растворимости конечных членов ряда ¿л;/ = ¿л ехр

(лг)'

где индекс и означает наличие напряжений несоответствия. Эти величины, подставленные в уравнения (7) - (8), определяют процессы сокристаллизации компонентов с учетом напряжений несоответствия. Равновесное значение коэффициента распределения теперь определяется

как г -к РУП__\_JlLL I, соответственно меняется и зависимость

Р[ ЯТдх^)

коэффициента распределения от скорости роста кристалла (рис. 2а), но

остается его асимптотическое стремление к К^ при бесконечно больших

скоростях роста. Напряжения несоответствия увеличивают свободную энергию кристалла, то есть они увеличивают растворимости твердых растворов. В связи с этим возникает дополнительное переохлаждение АТ(-„ наблюдаемое в эксперименте (Гликин, 2004; УЫовЫп е£ а1., 2004) и необходимое для кристаллизации из насыщенного раствора на затравку из неравновесного с раствором кристалла твердого раствора. При значительном несоответствии решеток срастание слоев может быть не вполне когерентным или сопровождаться релаксацией возникающих напряжений, что ведет к уменьшению АТс;. Наличие этого эффекта было продемонстрировано нами путем обсчета данных для системы бифталат калия-рубидия, полученных в работе (УоЬбЫп с\ а1., 2004).

Глава 2. Неоднородности распределения изоморфных компонентов на микроуровне (ростовая диссимметризация)

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

Хотя само по себе явление ростовой диссимметризации известно более 40 лет (Шубников, 1961; Цинобер, Самойлович, 1975; Ви1ка е! а!., 1980), его систематическое описание и анализ отсутствуют, поэтому нами проводилось детальное изучение основных аспектов

Табл. 1. Исследованные в работе твердые растворы, обладающие кинетическим упорядочением изоморфных атомов___

Соединение Сектор N Пр. гр.

роста Идеал. Набл.

Гранаты (гроссуляр-андрадит), Саз(А1,Ре)2(8Ю4)3 {110} 3 Ia3d Fddd, Л

Гранаты (гроссуляр-уваровит-гольдманит), Са3(А1,Сг,У)2(5Ю4)3 {110} 1 ¡a3d Л

Иа(С1,Вг)Оз {100} 2 />2,3 р\

Квасцы, (К,МН4)(А1,Сг)(804)2-12Н20 {111} 5 РаЗ R3

Квасцы, (НН4,ЯЬ)А1(804)2- 12НгО {111} 1 РаЗ РаЗ

(РЬ,8г)(Ш3)2 {111} 1 РаЗ R3

{100} 2 РаЗ РаЗ

(РЬ,Ва)(М03)2 {111} 1 РаЗ R3

Турмалины (эльбаит-лиддикоатит) (Са,№)(А1,и)зА16(816018)(В0з)з(0Н)з(0Н,Р) {0221} 3 Ют R\

Ы— число уточненных кристаллических структур

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

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

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

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

При ростовой диссимметризации центрированность ячейки не меняется, сверхструктуры не образуются, а метрика решетки искажается незначительно. Показателен пример фанатов уграндитовой серии Саз(А1,Ре,Сг)2(5Ю4)з, кубическая метрика решетки которых нарушается только незначительным отклонением угла [1 от 90° (рис. За).

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

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

0.2 0.4 xv,-xYr ат. ед.

б.

Рис. 3. а. Зависимость величины угла моноклинности от различий в заселенности позиций для гранатов гроссуляр-андрадитового (квадратики; по данным (Takeuchi et al., 1982; Allen, Buseck, 1988; Shtukenberg et al., 2005; Frank-Kamenetskaya et al., 2007)) и гроссуляр-уваровитового (кружки; по данным (Wildner, Andrut, 2001)) изоморфных рядов, сектор роста {110}; б. Характеристика диссимметризованных кристаллических структур фанатов Ca2Y2(Si04)3 (Y = AI,Fe,Cr) в координатах степень степень упорядоченности

Г - Г max(|x„-X. \,\x2i -*2J) ..

ij = i 2 ~ степень триклинности _ м ч\ i L>\/} ij =

pc. -Xy

1,2,..4, х\ и х2 - заселенности позиций У1 и У2. соответственно (уточнение в ромбической пр.гр. /чШ); хи ихъ~заселенности позиций УН, УН, соответственно (уточнение в триклинной пр.гр. /1). Вертикальные линии условно разграничивают области (слева направо) ромбических, псевдоромбических и триклинных гранатов

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

Численные оценки для гранатов гроссуляр-андрадитового изоморфного ряда показывают, что для разумных величин энергии активации отрыва трехвалентных катионов из поверхностного слоя (70100 кДж/моль) и нормальной скорости роста (10~8 - Ю'10 м/с) при температурах роста 150-500°С могут наблюдаться все возможные кристаллические структуры фанатов - от ромбических до полностью упорядоченных триклинных, что свидетельствует в пользу данной гипотезы. Другим ее подтверждением служит тот факт, что кристалл, обнаруживающий ромбическую симметрию, не проявляет субсекториальной неоднородности в распределении аномального

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

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

Табл. 2. Измеренные энергии активации перехода из двупреломляющего состояния в изотропное (снятие кинетического упорядочения атомов)___

Изоморфный ряд Энергия активации, кДж/моль

Низкие Т Высокие Т Лит-ра

Квасцы, (К,КН4)А1(804)2-12Н20 116(2) - 70-250

Квасцы, К(А1,Сг)(504)г12Н20 196(13) - 70-250

(РЬ^гХШз), 111(5) 359(17) 70-250

Гранаты, Са.;(А1,Ре)2(5Ю4)з 249(12) «676 155-290

Переход от низкой энергии активации к высокой в квасцах не наблюдается вплоть до температуры плавления; в нитратах он происходит при Т = 370-400°С, в гранатах - при Г» 1090°С. В графе «литература» для сравнения приведены энергии активации самодиффузии ионов в схожих соединениях

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

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

Ли

0.015- -0.015

/ • \

0.010- •/ • ° • \ -0.010 а •

.V V

0.005- tt М» -0.005

• • • •

□ _с • V»

*• п а *

О.ООО-у * " ?з uf 9 . |» ,_¡У ■«, Щ- 0.000

0 20 40 60 80 100 в Гроссуляр х, мол. % Андрадит

Рис. 4. Зависимость силы аномального двупреломлеиия (пропорциональной разнице в заселенностях позиций) от состава кристаллов: а. Квасцы (K1.i(NH4)J.)Al(S04)2-12H20, сектора роста октаэдра {111} (кружки), ромбододекаэдра {110} (квадратики) и куба {100} (треугольники). Линия -парабола; б. Хлорат-бромат натрия, сектора роста куба {100}. Сплошная линия -парабола; пунктирная линия - наилучшее описание экспериментальных данных с учетом линейного изменения межпозиционного коэффициента распределения; в. Гранаты гроссуляр-андрадитового и гроссуляр-уваровитового изоморфных рядов (х - средняя доли железа в октаэдрических позициях). Зачерненные значки - литературные данные разных авторов, полые значки - наши данные, треугольники - данные, полученные на одном зональном кристалле. Линия -парабола

смешения между позициями А2//?Г= 0.4). Кроме того, для случая гранатов гроссуляр-андрадитового изоморфного ряда (рис. 4в) вместо одной линии наблюдается поле точек, ограниченное сверху примерно параболической кривой. Подобное более сложное поведение нами связывается с тем фактом, что на упорядоченное состояние сильное влияние оказывает не только состав кристалла, но и условия его роста, прежде всего температура и скорость роста.

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

Аи, Ш5

Рис. 5. Зависимость силы аномального двупреломления от температуры роста для кристаллов квасцов (К0.5М).4о(^Н4)о.41)-о.(1о)А1(504)2-12Н20 в трех секторах роста

Дя, 10 t

2.0-1

1.5

1.0

0.5-

An, 10"4

•i

V. hm/c

r—^—i-

0.0

—i— 0.1

V, мм/сутки

0.2

Рис. 6. Зависимость силы аномального двупреломления кристаллов от скорости роста, а. Квасцы К(Л1].1Сгд)(804)2-12Н20, треугольники соответствуют х = 0.03 - 0.05 (Т = 24 - 37°С); квадратики -х = 0.04 - 0.05 (Г= 34 - 36 °С). Рост в режиме свободной (открытые символы) и вынужденной (зачерненные символы) конвекции; б. Нитраты (Pbo.sSro^XNChh. Рост в режиме вынужденной конвекции при температуре 37-390С

Глава 3. Секториальная неоднородность кристаллов

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

Первая группа не требует селективной адсорбции компонентов при встраивании в кристалл. Наиболее известна диффузионная гипотеза (Burton et al., 1953). Грани, растущие с разной скоростью, формируют различные градиенты концентрации компонентов в диффузионном пограничном слое, обеспечивая различные коэффициенты распределения компонентов для этих граней. Такой механизм действительно работает в некоторых случаях при выращивании

кристаллов из расплава, однако, как показывают оценки, он не может быть эффективным при росте кристаллов из растворов и должен быть исключен из рассмотрения для большей части минералов. Другой вариант (кинетическая гипотеза) рассматривался только в одной работе (Таусон, 2005). Так как коэффициент распределения зависит от скорости роста, то вследствие различия кинетических коэффициентов возможен разный захват компонентов в разные грани. Хотя примеры действия этого механизма в чистом виде нам неизвестны, численные оценки показывают, что он в принципе может работать при росте из растворов кристаллов квасцов, солей Туттона и т.д. Кроме того, этот эффект дает вклад в формирование секториальной неоднородности в тех кристаллах, где существенную роль играет адсорбционный механизм (см. ниже).

Наиболее популярна вторая группа механизмов, когда секториальность возникает за счет избирательной адсорбции компонентов на разных гранях кристалла (Леммлейн, 1948; Оолту, 1976). По своей сути адсорбционный механизм близок к механизму ростовой диссимметризации, только различия возникают не между позициями одной правильной системы точек на одной грани кристалла, а между позициями на разных гранях кристалла. Поэтому, так же как и для ростовой диссимметризации, межсекториальная разница в концентрациях компонентов должна примерно параболически увеличиваться к середине ряда. Действительно, такое поведение наблюдается в эксперименте (рис. 6а). В некоторых случаях максимум может быть смещен от центра ряда (рис. 66), что можно связать с неидеальностью смешения в поверхностных твердых растворах (для воспроизведения наблюдаемой кривой (рис. 66) достаточно взять различные энергии смешения в твердых растворах на разных гранях £>!/ЯТ= 0.1, 0^1 КГ= 0.3 ; при этом анализ фазовой диаграммы дает среднюю величину 0)КГ= 0.2). Возможно также гораздо более сложное и необычное поведение, что наблюдается, например, для ряда (РЬ,8г)(ЫОз)2. Здесь имеет место инверсия знака межсекториальной разницы в концентрациях компонентов Ах (рис. 7). При этом область инверсии распространяется примерно на половину ряда для роста при малых пересыщениях (« 0.01) и практически отсутствует для роста при

больших пересыщениях (0.1-0.33). Полученные зависимости Ах(х) могут быть смоделированы в рамках развиваемого феноменологического подхода (уравнения (7) - (9))-см. рис. 7. При низких пересыщениях зависимость Ах(х) определяется только вкладом адсорбционной

Рис. 6. Зависимость межсекториальных различий концентрации изоморфного компонента от его средней концентрации для изоморфных рядов солей Туттона (температура роста 36.5°С, пересыщение < 0.2, рост в режиме естественной конвекции), а. Ряд (Г*Щ4)2(№,СоХ804)2-6Н20, штриховая линия -парабола; б. Ряд (КН4)2(2п,Н1)(804)2'бН20 (квадратики). Для сопоставления треугольниками показаны результаты для схожего изоморфного ряда К.2(2п,№)(504)2-6Н20 (пересчет данных (Жмурова, Хаимов-Мальков, 1970)). Сплошная линия - расчет в предположении близости к равновесию по модели адсорбционной секториальной селективности. Пунктирная линия — расчет по модели кинетической секториальной селективности

Рис. 7. Зависимость межсекториальной разницы в концентрациях компонентов твердого раствора (8ггРЬ,^)(К'03)2 от состава кристалла, а. Низкие пересыщения я 0.01; б. Высокие пересыщения 0.1-0.33.Треугольники - рост в режиме естественной конвекции при Т = 31.6-32°С, квадратики - рост при принудительном перемешивании раствора при 7, = 32-40°С, звездочка - данные работы (\Уи е! а!., 2000) для роста путем испарения воды в режиме естественной конвекции при Т = 40-45°С. Линии - результаты моделирования

воспроизводится, если принять различные параметры смешения компонентов твердого раствора на поверхности каждой грани: (<2\ -<2г)/ЯТ= 0.23. При высоких пересыщениях начинает расти вклад кинетической секториальной селективности.

Еще одним важным аспектом адсорбционной секториальности является ее связь с атомной структурой растущей грани. Собственно, именно этой связи и посвящено подавляющее большинство работ о секториальности кристаллов (Но1Ш1ег, Оапсагг, 1971; Оои^у, 1976; СЬакЬпюигасНап е£ а1., 2007). Подобный анализ был выполнен нами для 12 групп твердых растворов, и было установлено, что имеющиеся простые оценки предпочтительности встраивания в данную грань того или иного атома, основанные на правиле полярности Гольдшмидта, дают корректные результаты только для существенно ионных кристаллов и не работают для ионно-ковалентных соединений. В частности, для нескольких изоморфных рядов солей Туттона (МН^Ме^О^-бН^О (Ме = N4, Со, Хп, Г^) было показано, что более крупные ионы преимущественно входят в более свободные позиции на грани {110} по сравнению с более тесными позициями на грани {001}, причем максимальная межсекториальная разница в концентрациях компонентов коррелирует с различиями в ионных радиусах соответствующих ионов.

Глава 4. Зональная неоднородность кристаллов

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

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

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

Данная модель применялась нами для анализа механизмов возникновения зональности фанатов (Fe,Al,Mn,Ca)3Al2(Si04)3 из метаморфических пород, практически всегда обладающих непрерывной трендовой зональностью в распределении концентрации компонентов (Tracy, 1982). Наиболее распространена (особенно в низкотемпературных <600 °С породах) прямая зональность, когда концентрация марганца падает, а магния растет от центра к краю зерна. В литературе предполагается, что такая зональность возникла либо за счет фракционирования компонентов в закрытой системе, либо за счет изменения температуры при профессивном метаморфизме. В более высокотемпературных породах обычна обратная зональность, которую связывают с рефессивным метаморфизмом. Наконец, иногда наблюдаются сложные (например, с выполаживанием на краях зерна) или даже немонотонные профили зональности, возникновение которых приписывают полиметаморфизму.

Расчёты профилей зональности были выполнены для гранатов, входящих в изоморфный ряд "альмандин Fe3Al2Si3012 - пироп Mg3Al2Si30I2". Растворимости миналоь фаната в водном флюиде были вычислены нами с помощью профамм SUPCRT92 (Johnson et al., 1992), GBFLOW (Гричук, 1996) и FLUID (Skvirsky, 1995) по термодинамическим данным Бермана и Арановича (Berman, Aranovich, 1996) и Хелгесона (Shock, Helgeson, 1988; Tanger, Helgeson, 1988).

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

большие скорости нагрева (2000°С/Ма) дают незначительные различия в концентрациях (Ах < 1 мол. %). Обратная зональность может быть воспроизведена при моделировании кристаллизации в щелочных или в высокотемпературных (>600°С) нейтральных растворах, когда имеет место инверсия коэффициента распределения. Сложные (в том числе немонотонные) профили зональности также могут быть получены за счет фракционирования компонентов в проточной системе. Таким образом, зональность кристаллов метаморфических гранатов определяется не столько вариациями температуры и давления, сколько фракционированием компонентов и составом поступающего флюида. Такое объяснение сложной зональности фанатов выглядит гораздо более правдоподобным, чем распространенные в литературе предположения о сложных перемещениях значительных участков земной коры для обеспечения необходимых трендов изменения температуры и давления - см., например, (Schumacher et al., 1999; Garciä-Casco et al., 2002).

Если способы образования трендовой зональности понятны, то механизмы возникновения осцилляционной зональности совершенно не ясны. Считается, что такая зональность может возникнуть либо за счет внешних колебаний состава среды кристаллизации или условий роста, либо за счет автоколебательных процессов на фронте роста кристалла (Shore, Fowler, 1996). Однако детальный анализ литературных данных показал, что практически нет ни одного примера, где был бы выявлен и подтвержден конкретный механизм возникновения колебаний состава кристалла. В задачу данной работы входит поиск подобных механизмов.

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

периодичность изменения концентраций компонентов. Величина скейлинговых параметров в среднем несколько ниже 0.5, что указывает на слабое антиперсистентное поведение и достаточно высокую роль случайных процессов при формировании зональности. Размерность странного аттрактора колеблется в интервале от 2 до 4, то есть процесс формирования осцилляционной зональности может быть описан 3-5 обыкновенными дифференциальными уравнениями первого порядка (Мун, 1990). Все эти данные свидетельствуют о том, что возникновение

Табл. 3. Статистические характеристики осцилляционной зональности кристаллов_______

Образец N L, цш Л, цш щ Hi. Мтах

Квасцы (Ko.77R.bo r,)Al(S04)r 12Н20 980 2238 113 0.52 0.18 2.5/3.1

Квасцы (К„ ,(NH4)n 0Al(SO4)212H2O 1482 2964 87 0.55 0.14 2.9 / 3.7

(Sr,Pb)(N03)2, образец N39.1 904 3470 154 «1 0.11 2.5 / 2.7

(Sr,Pb)(N03)2, образец N26.2tw 847 3710 59 0.55 0.13 2.6

Основные плагиоклазы

Толбачик, образец №1 (Фрагмент 1) 1800 1931 47 0.40 0.53 3.6/4.0

Толбачик, образец №1 (Фрагмент 2) 2006 2152 145 0.39 0.63 3.2/3.8

Толбачик, образец №2 2500 1700 36 0.18 0.28 3.8/4.3

Ключевская сопка, образец №5 2600 4000 60 0.07 0.27 4/4.7

Ключевская сопка, образец №6 559 1000 41 0.25 0.14 3.5/3.9

Коростеньский плутон, образец №7 260 5000 - 0.36 0.9 2.6

Пирохлор

Керимаси, образец LL1 500 42 - 0.55 0.55 2.8

Керимаси, образец LL2 341 49 2.9 0.34 0.4 2.7

Керимаси, образец KZ2 285 45 6.5 0.33 0.39 3.2

Ковдор, образец RU 366 47 6.7 0.42 0.32 2.7/3.0

Ковдор, образец LD 519 98 11.4 0.52 0.34 2.5/3.8

Обсчет данных (Hogarth et al., 2000) 257 275 17.4 0.29 0.68 2.9

Гранаты, Cai(Al,Fe)2(Si04)3. Дальнегорское месторождение

Образец Grtl 1746 1000 8.6 0.68 0.14 4.7

Образец Grt4 1627 1000 8.5 0.56 0.08 5.1

Блеклые руды, обсчет данных (Плотинская и др., 2005)

Образец El lb 563 12.4 0.55 0.41 0.40 4.0

Образец Oz2e 535 158 - >0.1 0.37 5.5

Образец Ozla 401 136 - 0.38 0.48 3.3

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

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

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

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

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

Другим более общим способом возникновения осцилляционной зональности является механизм микроморфологической неустойчивости грани под действием сильно адсорбирующихся примесей, тормозящих рост. В результате периодического исчерпания и накопления пересыщения на фронте роста с образованием и распадом макроступеней и отравлением и «очисткой» поверхности возникают колебания скоростей роста. Зависимость коэффициента распределения от скорости роста (см. главу 1) приводит к формированию осцилляционной зональности. Таким механизмом, например, объясняется зональность кристаллов группы дигидрофосфата калия, что подтверждается наличием у этих кристаллов гистерезиса скоростей роста (Пунин, Артамонова, 1986). Мы предполагаем действие этого механизма и при формировании осцилляционной зональности кристаллов (РЬ,5г)^03)2, а также считаем, что с его помощью можно объяснить возникновение осцилляционной зональности и во многих минералах, в частности, в плагиоклазах (отравляющие примеси - Бе, К).

Глава 5. Неоднородности состава и внутренние напряжения в кристаллах

Выше были рассмотрены различные типы конституционной неоднородности кристалла на разных масштабных уровнях. Хорошо известно, что неоднородности состава кристалла приводят к сопряжению в нем областей с различными параметрами элементарной ячейки (явление гетерометрии (Штернберг, 1962)). Вследствие этого в кристалле развиваются внутренние напряжения. В главе 1 мы показали роль упругих напряжений в смещении равновесия «твердый раствор-водный раствор», изменении коэффициента распределения и скорости роста кристалла. Однако влияние внутренних напряжений на рост кристалла этим не ограничивается. Задача данной главы - проследить

пути возникновения, эволюции и релаксации напряжений конституционной гетерометрии.

Точный расчет поля напряжений гетерометрии - трудная задача, так как распределение по кристаллу участков с разным составом обычно сложное, а изменение параметров элементарной ячейки на границах участков скачкообразное. Нами были созданы приближенные модели, позволяющие рассчитать напряжения гетерометрии для наиболее важных случаев зональной неоднородности кристаллов (плоская зональность в тонкой пластинке и высокой балке; концентрическая зональность в диске, цилиндре и шаре). С помощью этих моделей нами рассчитывались величины и распределение напряжений гетерометрии для ряда природных (скарновые фанаты, полихромные турмалины, бериллы) и синтетических (дигидрофосфат калия-аммония, квасцы) соединений. Полученные результаты сопоставлялись с напряжениями, непосредственно определенным по оптическим аномалиям с помощью поляризационного-оптического метода (Сиротин, Шаскольская, 1979). Было обнаружено, что расчетное распределение напряжений гетерометрии достаточно хорошо соответствует экспериментальным данным (см., например, рис. 8), однако по абсолютной величине измеренные напряжения обычно ниже расчетных. Таким образом, в кристаллах значительная часть напряжений гетерометрии либо релаксирует во время и после роста, либо вообще не возникает.

ат. % 8 6 4-1 2

А1

ттах,ГПа Ап, 10

■ . - V _

.Да/"

Ре

______'

0.4

0.2

• /

,__ЛА Л

5 6 /-, ММ

3 4 5 6 7 8 а. Лмм б.

Рис. 8. Аномальное двупреломление, вызываемое межзональными напряжениями конституционной гетерометрии в кристалле граната, а. Изменение концентраций А1 и Ре (зачерненные значки) и величина максимальных скалывающих напряжений (полые значки) в зависимости от расстояния от центра кристалла; б. Рассчитанная (полые значки) и экспериментально измеренная (зачерненные значки) силадвупреломления

Действительно, величины возможных напряжений гетерометрии достигают 10 ГПа и более, что должно давать существенную добавку к величине свободной энергии кристалла. Поэтому кристалл при росте стремится минимизировать уровень внутренних напряжений. При значительных несоответствиях решеток рост происходит островками, в промежутках между которыми образуются различные дефекты. Поэтому слой нового состава сопрягается с предыдущим не по всей площади, что ведет к снижению уровня напряжений. Такой островковый рост описан в литературе для ряда объектов. Нами он наблюдался in situ методом ACM при гетероэпитаксиальном росте хашемита BaCrCU на барите BaS04 (несоответствие параметров решетки до 2.6 %).

Таким образом, уже исходный уровень напряжений гетерометрии может быть существенно ниже теоретически возможного. Далее напряжения могут релаксировать в процессе роста с образованием пластических дефектов. Как правило, дефекты, снижающие уровень исходных внутренних напряжений, порождают новые дальнодействующие напряжения, которые, в свою очередь, релаксируют с образованием новых дефектов. Запускается автокаталитический процесс автодеформации, с возникновением «ростовых» дислокаций, блокования, ростового изгиба-кручения, расщепления, двойникования, ростовых мартенситных превращений (Пунин, Штукенберг, 2008).

Важно отметить, что степень релаксации напряжений тем больше, чем выше исходный уровень напряжений. Соответственно остаточные напряжения будут выше в тех кристаллах, где исходные напряжения были относительно небольшими. Этот парадоксальный вывод следует из того, что уровень напряжений, необходимый для запуска процессов релаксации, превосходит уровень напряжений, необходимый для развития и эволюции дефектной структуры. Подобное поведение наблюдалось нами для кристаллов целого ряда соединений. Показательным примером являются бизональные кристаллы квасцов. При нарастании квасцов состава К(А 10ИСг(,07)(8О4)2-12Н20 на алюмо калиевые квасцы KA1(S04)2-12H20 (несоответствие решеток А а!а = 0.0005) теоретические напряжения равны 6 МПа. Измерения

аномального дву прело мления дают распределение напряжений по кристаллу близкое к расчетному, но значения напряжений в несколько раз меньше расчетных (максимально 1.6 МПа). Поскольку дефекты на границе срастания не обнаруживаются, можно предположить существенную некогерентность срастания слоев. С другой стороны, при нарастании алюмоаммонийиых квасцов ЫН4Л1(504)2- 12Н20 на злюмскалиевыс КА1(504)2-12Н20 (несоответствие решеток Да/а = 0.007) теоретические напряжения гораздо больше и равны 90 МПа. Однако поляризационно-оптический метод показывает отсутствие упругих напряжений в кристалле. При этом на границе срастания возникает большое количество ростовых дислокаций.

Напряжения гетерометрии могут релаксировать не только пластически, но и хрупко, причем хрупкие автодеформации (образование трещин) изучены очень плохо. Они происходят при достижении предела прочности кристалла на растяжение сг^,, равного по имеющимся оценкам, 10'3£, где Е - модуль Юнга. Для зональных кристаллов фанатов гроссуляр-андрадитового ряда нами показано (рис. 9а), что коэффициент трещиноватости, характеризующий интенсивность хрупкой автодеформации, прямо пропорционален величине нормальных растягивающих напряжений гетерометрии. При этом образуются системы параллельных трещин, ориентированных нормально или почти нормально к границам зон разного состава. Распределение трещин в таком ансамбле подчиняется закону Вейбулла (распределение слабого звена) - рис. 96. Другим характерным морфологическим типом хрупкой автодеформации является возникновение трещин в местах концентрации напряжений, особенно на изломах зональных и секториальных границ. В целом, в отличие от пластических, хрупкие автодеформации обычно не порождают новых дальнодействующих напряжений, то есть автокаталитическое развитие дефектности и образование сложных структур в этом случае отсутствуют.

К , 1/мм

тр'

О < О

кр.

ст > а,

кр.

Р(х) 0.30.20.1-1

0.0-=^

т.

0.3 0.6

х, мм

т-1-1-г

-2-1 0 1

а ог,ГПа б.

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

Выводы

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

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

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

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

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

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

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Монографии и обзорные статьи

Штукенберг А.Г., Пунин Ю.О. Оптические аномалии в кристаллах. СПб. Наука. 2004. 264 с.

Shtukenberg A.G., Punin Yu.O. Optically anomalous crystals (Ed. B. Kahr) Springer. Dordrecht. 2007. 253 p.

Пунин Ю.О., Штукенберг А.Г. Автодеформационные дефекты кристаллов. Изд-во СПбГУ. 2008. 318 с.

Штукенберг А.Г., Пунин Ю.О. Оптические аномалии в кристаллах // Записки Всероссийского Минералогического Общества. 1996. Ч. 125. №4. С. 104-120.

Штукенберг А.Г., Пунин Ю.О., Франк-Каменецкая О.В. Кинетическое упорядочение и ростовая диссимметризация кристаллических твердых растворов // Успехи химии. 2006. Т. 75. С. 1212-1236.

Shtukenberg A.G., Punin Yu.O., Azimov P. Crystallization kinetics in binary solid solution - aqueous solution systems // American Journal of Science. 2006. Vol. 306. P. 553-574.

Статьи в реферируемых журналах

Штукенберг А.Г., Пунин Ю.О., Котельникова E.H. Железистые мусковиты "барботов глаз" из Ильменских гор (о природе изгиба) // Записки Всероссийского Минералогического Общества. 1993. Ч. 122. № 5. С. 53-63.

Штукенберг А.Г., Пунин Ю.О., Котельникова Е.Н., Сухаржевский С.М. Аномальная оптика изоморфно-смешанных кристаллов // Вестник СПбГУ. 1994. Сер. 4. Вып. 2. С. 109-110.

Штукенберг А.Г., Пунин Ю.О., Котельникова Е.Н., Сухаржевский С.М. Оптическая диссимметризация кристаллов в связи с неоднородностью распределения изоморфных компонент // Журнал структурной x'iMiivi. 1S94. Т. 35. С. 60-69.

Штукенберг А.Г., Соловьев В.Н., Пунин Ю.О. Аномальное двупреломление в кристаллах изоморфных серий квасцов // Физика кристаллизации. Тверь. 1994. С. 76-81.

Штукенберг А.Г., Пунин Ю.О., Ковалев О.Г.Температурное поведение оптических аномалий в кристаллах квасцов // Кристаллография. 1998. Т. 43. С. 505-508.

Ivanova Т.1., Shtukenberg A.G., Punin Yu.O., Frank-Kamenetskaya O.V., Sokolov P.B. On the complex zonality in grandite garnets and implications H Mineralogical Magazine. 1998. Vol. 62. P. 857-868.

Shtukenberg A.G., Punin Yu.O., Soloviev V.N. Effect of growth conditions on the birefringence of mixed crystals revealed in alum solid solutions// Mineralogical Magazine. 2000. Vol. 64. P. 837-845.

Азимов П.Я., Штукенберг А.Г. Расчёт фазовых диаграмм водно-солевых систем с твёрдыми растворами // Журнал Неорганической Химии. 2000. Т. 45. С. 1424-1432.

Штукенберг А.Г., Франк-Каменецкая О.В., Баннова И.И., Рождественская И.В., Кулакова С.В., Иванова Т.И., Пунин Ю.О. Диссимметризация кристаллических структур аномально-двупреломляющих квасцов // Кристаллография 2000. Т. 45. С. 999-1005.

Рождественская И.В., Франк-Каменецкая О.В., Штукенберг А.Г., Баннова И.И. Триклинная структура двупреломляющего кристалла квасцов К(А10<j5Cro.o5)(S04)2- 12Н20 //Журнал структурной химии. 2001. Т. 42. С. 753-765.

Shtukenberg A.G., Punin Yu.O., Frank-Kamenetskaya O.V., Kovalev O.G. and Sokolov P.B. On the origin of anomalous birefringence in grandite garnets // Mineralogical Magazine. 2001. Vol. 65. P. 445-459.

Штукенберг А.Г., Азимов П.Я. Учет неидеальности твердого раствора при расчете фазовых диаграмм водно-солевых систем с

непрерывными твердыми растворами // Журнал Неорганической Химии. 2001. Т. 46. С. 1755-1758.

Азимов П.Я., Штукенберг А.Г. Термодинамический анализ факторов, определяющих ростовую зональность в метаморфических гранатах // Доклады РАН. 2001. Т. 376. С. 235-237.

Shtukenberg A.G., Punin Yu.O., Haegele E., Klapper H. On the origin of inhomogeneity of anomalous birefringence in mixed crystals: an example of alums // Physics and Chemistry of Minerals. 2001. Vol. 28. P. 665-674.

Shtukenberg A.G., Popov D.Yu., Punin Yu.O. An application of the point-dipole model to the problem of optical anomalies in grandite garnets // Mineralogical Magazine. 2002. Vol. 66. P. 275-286.

Shtukenberg A.G., Punin Yu.O. Trend of anomalous birefringence in alum crystals: a new effect of the growth ordering phenomenon // Neues Jahrbuch fuer Mineralogie. Monatshaefte. 2003. P. 241-253.

Azimov P., Shtukenberg A. Numerical modelling of growth zoning at nonstationary crystallization of solid solutions - Metamorphic garnets // Mathematical Geology. 2003. Vol. 35. P. 405-430.

Shtukenberg A.G., Rozhdestvenskaya I.V., Popov D.Yu., Punin Yu.O. Kinetic ordering of atoms in sodium chlorate-bromate solid solutions // Journal of Solid State Chemistry. 2004. Vol. 177. P. 4732-4742.

Shtukenberg A.G., Astilleros J.M., Putnis A. Nanoscale observations of epitaxial growth of hashemite on barite (001) // Surface Science. 2005. Vol. 590. P. 212-223.

Shtukenberg A.G. Metastability of atomic ordering in lead-strontium nitrate solid solutions // Journal of Solid State Chemistry. 2005. Vol. 178. P. 2608-2612.

Shtukenberg A.G., Popov D.Yu., Punin Yu.O. Growth ordering and anomalous birefringence in ugrandite garnets // Mineralogical Magazine. 2005. Vol. 69. P. 537-550.

Денисов A.B., Пунин Ю.О., Габриелян B.T., Грунский О.С., Штукенберг А.Г. Физические и физико-химические процессы, сопровождающие синтез шихты, выращивание и отжиг кристаллов РЬМо04 в различных средах. Часть 4. Аномальная двуосность и

термопластические напряжения в кристаллах РЬМо04 в их взаимосвязи с условиями выращивания // Кристаллография. 2006. Т. 51. С. 136-142.

Денисов A.B., Лунин Ю.О., Штукенберг А.Г., Грунский О.С. Исследование оптических аномалий и расчет внутренних напряжений в кристаллах молибдата свинца, выращенных методом Чохральского // Вестник СПбГУ. 2006. Сер. 7. Вып. 1. С. 106-110.

Shtukenberg A.G., Calci H., Kiríei A., Popov D.Yu. Symmetry reduction and cation ordering in solid solutions of strontium-lead and barium-lead nitrates// Zeitschrift fur Kristallographie. 2006. Vol. 221. P. 681-688.

Shtukenberg A.G., Euler H., Kirfel A. Symmetry reduction and cation ordering in alum solid solutions // Zeitschrift fiir Kristallographie. 2007. Vol. 222. P.73-82.

Shtukenberg A., Rozhdestvenskaya I., Frank-Kamenetskaya O., Bronzova J., Euler H., Kirfel A., Bannova I., Zolotarev A. Symmetry and crystal structure of biaxial elbaite-liddicoatite tourmaline from the Transbaikalia region, Russia // American Mineralogist. 2007. Vol. 92. P. 675686.

Frank-Kamenetskaya O.V., Rozhdestvenskaya I.V., Shtukenberg A.G., Bannova I.I., Skalkina Yu.A. Dissymmetrization of crystal structures of grossular-andradite garnets Ca3(Al,Fe)2(Si04)3 // Structural Chemistry. 2007. Vol. 18. P. 493-503.

Shtukenberg A.G., Punin Yu.O., Azimov P.Ya. Reply to comment: Supersaturation in binary solid solution — aqueous solution systems by M. Prieto, J.M. Astilleros, C.M. Pina, L. Fernández-Díaz and A. Putnis // American Journal of Science. 2007. Vol. 307. P. 1046-1050.

Штукенберг А.Г., Пунин Ю.О., Азимов ПЛ. Феноменологический анализ процессов взаимодействия кристалла и среды в системах «твердый раствор - водный раствор» // Записки Российского минералогического общества. 2007. Ч. 136. № 7. С. 36-52.

Uher Р., Kovácik M., Kubis M., Shtukenberg A., Ozdín D. Metamorphic vanadian-chromian silicate mineralization in carbon-rich amphibole schists from the Male Karpaty Mountains, Western Carpathians, Slovakia // American Mineralogist. 2008. Vol. 93. P. 63-73.

Подписано в печать 12.01.2009 г. Формат 60x84/16 Бумага офсетная. Печать офсетная. Усл. печ. л. 2,3. Тираж 150 экз. Заказ №

Типография Издательства СПбГУ. 199061, С,- Петербург, Средний пр., 41.

Содержание диссертации, доктора геолого-минералогических наук, Штукенберг, Александр Григорьевич

Введение.

1. Особенности кристаллизации в системах «твердый раствор-водный раствор»

1.1. Введение.

1.2. Описание модельной системы.

1.3. Равновесие, простые системы (В,С)А - Н2О.

1.4. Равновесие, общий случай.

1.5. Квазиравновесия.'.

1.6. Рост кристалла: термодинамическое рассмотрение.

1.7. Рост кристалла: кинетическое рассмотрение, кинетический режим кристаллизации.

1.8. Рост кристалла: кинетическое рассмотрение, модификации для смешанного режима кристаллизации.

1.9. Взаимодействие между раствором и кристаллом.

1.10. Пересыщение.

1.11. Влияние напряжений несоответствия.

1.12. Заключительные замечания.

2. Неоднородности распределения изоморфных компонентов на микроуровне ростовая диссимметризация).V.

2.1. Введение.

2.2. Методы выявления ростовой диссимметризации кристаллов.

2.2.1. Кристаллооптика.

2.2.2. Электронный парамагнитный резонанс.

2.2.3. Инфракрасная и рамановская спектроскопия.

2.2.4. Электронная микроскопия.

2.2.5. Дифракционные методы.

2.3. Симметрийный аспект диссимметризации.

2.4. Оптическая и структурная неоднородность кристаллов.

2.4.1. Секториальность.

2.4.2. Субсекториальность.

2.4.3. Взаимоотношения между секториальностью и субсекториальностью.

2.4.4. Зональность.

2.4.5. Тренд аномального двупреломления.

2.5. Изменения кристаллической структуры при диссимметризации.

2.6. Метастабильность кинетического упорядочения.

2.7. Кинетическое упорядочение и атомная структура растущей грани.

2.8. Диффузионная релаксация.

2.9. Факторы, управляющие кинетическим упорядочением.

2.9.1. Состав твердого раствора.

2.9.2. Состав среды кристаллизации.

2.9.3. Температура и давление.

2.9.4. Скорость роста.

2.9.5. Гидродинамика и массоперенос в растворе.

2.10. Заключительные замечания.

3. Секториальная неоднородность кристаллов.

3.1. Сущность явления и классификация.

3.2. Общая схема возникновения секториальности.

3.3. Механизм 1: Диффузионная секториальная селективность.

3.4. Механизм 2а: Кинетическая секториальная селективность.

3.5. Механизм 26: Адсорбционная секториальная селективность.

3.6. Механизм За: Диффузионная релаксация на поверхности грани.

3.7. Механизм 36: Диффузионная релаксация в приповерхностных слоях кристалла

3.8. Закономерности формирования секториального кристалла.

3.8.1. Вклад различных механизмов.V.

3.8.2. Влияние состава кристалла и условий кристаллизации.

3.9. Кристаллохимический аспект адсорбционной секториальности.

4. Зональная неоднородность кристаллов.

4.1. Введение.

4.2. Трендовая зональность.

4.2.1. Математическая модель образования трендовой зональности.

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

4.2.3. Моделирование зональности кристаллов метаморфических гранатов.

4.2.4. Испарение растворителя.

4.3. Осцилляционная зональность.

4.3.1. Введение.

4.3.2. Механизмы возникновения осцилляционной зональности.

4.3.3. Микроморфологическая неустойчивость растущей грани.

4.3.4. Ингибирование роста кристалла микропримесью.

4.3.5. Описание осцилляционной зональности.

4.3.6. Зональность в плагиоклазах.

4.3.7. Зональность в пирохлоре.

5. Неоднородности состава и внутренние напряжения в кристаллах.

5.1. Введение.

5.2. Расчет напряжений гетерометрии.

5.2.1. Основные положения.

5.2.2. Плоская зональность в балке или пластинке.

5.2.3. Плоская зональность в толстых пластинах.

5.2.4. Концентрическая зональность в цилиндре.

5.2.5. Концентрическая зональность в шаре.

5.2.6. Сопоставление результатов расчетов по разным моделям.

5.2.7. Общий случай.

5.2.8. Секториальность.

5.3. Ограничения на величину напряжений гетерометрии.

5.3.1. Кристаллохимические ограничения. S J \

5.3.2. Влияние условий кристаллизации.

5.3.3. Смещение фазовых равновесий.

5.3.4. Формирование некогерентной границы.".Г:.>.

5.3.5. Способы релаксации внутренних нарпряжений.

5.4. Внутренние напряжения и искажения оптической индикатрисы.

5.4.1. Пьезооптический эффект.

5.4.2. Дигидрофосфат калия-аммония.

5.4.3. Гранаты гроссуляр-андратитового изоморфного ряда.

5.4.4. Бериллы.

5.4.5. Секториальность.

5.4.6. Общие особенности оптических аномалий, связанных с напряжениями гетерометрии.

5.5. Хрупкая релаксация напряжений.

5.5.1. Пределы прочности, связь с внутренними напряжениями.

5.5.2. Морфология хрупкой автодеформации.

5.5.3. Геометрия расположения трещин в зональных кристаллах.

5.5.4. Распределение расстояний между трещинами в ансамбле.

5.5.5. Влияние различных факторов на возникновение трещин.

5.6. Общая схема возникновения и релаксации напряжений гетерометрии.

Выводы.

Введение Диссертация по наукам о земле, на тему "Формирование пространственных неоднородностей в распределении изоморфных компонентов при кристаллизации твердых растворов"

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

Объекты исследования. Экспериментальные исследования проводились на кристаллах твердых растворов, выращиваемых в контролируемых условиях из низкотемпературных водных растворов (квасцы, соли Туттона, дигидрофосфат калия-аммония, хлорат-бромат натрия, нитраты двухвалентных металлов и др.), а также на ряде природных минералов различного генезиса (гранаты, берилл, турмалины, плагиоклаз, пирохлор и др.). Для проверки разрабатываемых теоретических моделей использовались как оригинальные, так и литературные данные по различным соединениям.

Научная новизна.

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

Защищаемые положения.

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

Апробация работы и публикации. Материалы работы опубликованы в 3 монографиях, 33 статьях в реферируемых журналах, 58 тезисах докладов. Они были представлены на конференциях разного уровня, в том числе: Международное совещание по рентгенографии и кристаллохимии минералов (Рентгенография минерального сырья) (1992, 1999, 2003, 2007); Международная конференция «Спектроскопия, рентгенография и кристаллохимия минералов» (1997); International Conference of Crystal Growth (1998); European Crystallographic Meeting (1998); Congress and General Assembly of the International Union of Crystallography (1993,1996,1999); Уральское кристаллографическое совещание (1998); Национальная конференция по росту кристаллов (2000, 2002, 2004, 2006); Международная конференция "Кристаллогенезис и минералогия" (2001,2007); Совещание по экспериментальной минералогии (1995, 2001); International Symposium on Experimental Mineralogy, Petrology and Geochemistry (2002); General Meeting of the International Mineralogical Association (1994, 2002); Международный симпозиум «Минералогические музеи» (2002, 2005); Национальная кристаллохимическая конференция (2003, 2006); Международная конференция по физике кристаллов «Кристаллофизика 21-го века» (2006).

Объем и структура работы. Диссертация содержит 414 страниц текста, 131 рисунок, 27 таблиц, список литературы, насчитывающий 684 наименования, и состоит из 5 глав. Первая глава посвящена процессам образования кристаллов твердых растворов из водных растворов. В ней анализируется поведение коэффициентов распределения компонентов. Во второй - четвертой главах рассматриваются основные типы пространственных неоднородностей в распределении изоморфных компонентов по кристаллу — кинетическое упорядочение атомов, секториальность и зональность, 9 соответственно. Наконец, пятая глава анализирует процессы возникновения и релаксации внутренних напряжений, возникающих при росте неоднородных по составу кристаллов.

Благодарности. Прежде всего, выражаю глубокую благодарность моему научному консультанту Ю.О. Пунину. Хочу поблагодарить П.Я. Азимова, О.В. Франк-Каменецкую, И.В. Рождественскую, Г. Клаппера, А. Кирфеля, X. Ойлера, Д.Ю. Попова, М.В. Морозова за плодотворное сотрудничество; И.И. Баннову, О.И. Артамонову, Е.Н. Котельникову, А. Путниса, Х.М. Астиллероса, А.Э. Волошина за помощь при выполнении отдельных этапов работы; П.Б. Соколова, А. А. Золотарева, А.Н. Зайцева за предоставление материалов для исследования, а также студентов и аспирантов кафедр кристаллографии и минералогии геологического факультета СПбГУ за помощь в получении экспериментальных данных.

Заключение Диссертация по теме "Минералогия, кристаллография", Штукенберг, Александр Григорьевич

Выводы

1. Разработка феноменологической модели кристаллизации в системах «твердый раствор - водный раствор»

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

1.2. При малых отклонениях от равновесия стационарному росту отвечает равенство разностей химических потенциалов компонентов; при больших отклонениях от равновесия это равенство нарушается (тем сильнее, чем выше степень неравновесности). Коэффициент распределения вблизи равновесия определяется фазовой диаграммой системы, при увеличении отклонения от равновесия он меняется, асимптотически стремясь к постоянному значению. Это значение определяется соотношением кинетических коэффициентов в кинетическом режиме роста кристалла, комбинацией кинетических коэффициентов и коэффициентов диффузии в смешанном режиме и соотношением коэффициентов диффузии — в диффузионном. При удалении от равновесия величина коэффициента распределения может перейти через единицу.

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

1.4. Напряжения гетерометрии (несоответствия), возникающие при срастании слоев с разным соотношением компонентов, способны существенно изменить растворимость твердого раствора, уменьшить скорость его роста и привести к изменению коэффициента распределения.

2. Экспериментальное и теоретическое исследование явления ростовой диссимметризации

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

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

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

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

2.5. Из феноменологического рассмотрения процессов встраивания-отрыва индивидуальных компонентов на фронте роста кристалла определены основные факторы, влияющие на кинетическое упорядочение атомов. Степень упорядочения: (1) зависит от состава кристалла и увеличивается к середине изоморфного ряда; (2) падает при увеличении температуры кристаллизации и (3) в целом падает при увеличении скорости роста кристалла. Эксперимент целиком подтверждает теоретические выводы.

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

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

3. Исследование механизмов образовании секгориальной неоднородности кристаллов

3.1. Систематизированы и проанализированы возможные причины возникновения и стирания секториальности кристаллов, включающие в себя: (1) различия в скорости диффузионного переноса компонентов к граням кристалла; (2) различия в скорости поверхностной реакции на разных гранях кристалла; (3) избирательную адсорбцию компонентов на разных гранях при встраивании в кристалл; (4) диффузионную релаксацию в поверхностном слое и в объеме кристалла. Подтверждена доминирующая роль механизма избирательной адсорбции компонентов в образовании секториальной неоднородности кристаллов.

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

361 феноменологического рассмотрения процессов встраивания-отрыва индивидуальных компонентов в разные грани кристалла определены основные факторы, влияющие на секториальную неоднородность кристаллов. Межсекториальные различия в концентрациях компонентов: (1) зависят от состава кристалла и увеличиваются к середине изоморфного ряда; (2) уменьшаются при увеличении температуры кристаллизации и (3) в целом падают при увеличении скорости роста кристалла. Экспериментальные данные согласуются с теоретическими результатами.

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

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

4. Анализ путей возникновения трендовой и осцилляционной зональности в кристаллах

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

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

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

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

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

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

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

5.2. Согласно расчетам, внутренние напряжения в кристаллах могут достигать очень больших величин (до Ю10 н/м2). Однако по экспериментальным данным они гораздо меньше, так как кристалл сопротивляется развитию значительных напряжений на всех этапах своей эволюции. Теоретическое рассмотрение и эксперименты дают следующие способы саморегуляции и снижения ростовых внутренних напряжений: (1) сдвиг фазового равновесия, за счет чего участки кристалла с высоким уровнем напряжений либо не возникают, либо растут медленней и имеют меньшие размеры (что облегчает релаксацию напряжений в них); (2) некогерентное сопряжение участков разного состава (зон, секторов, субсекторов); (3) образование дефектов, препятствующих формированию дальнодействующих полей напряжений (дислокации несоответствия, трещины, включения); (4) образование дефектов, уменьшающих напряжения несоответствия, но при этом создающих новые дальнодействующие напряжения (мартенситные превращения, блокование, расщепление, ростовое изгиб-кручение, ростовое двойникование).

Библиография Диссертация по наукам о земле, доктора геолого-минералогических наук, Штукенберг, Александр Григорьевич, Санкт-Петербург

1. Авченко О.В. Петрогенетическая информативность гранатов метаморфических пород. М.: Наука. 1982.

2. Азимов П.Я. Формирование ростовой зональности в метаморфических гранатах (роль Р-Т условий и флюидного режима) // Автореф. дисс. к.г.-м.н. Санкт-Петербург. 2000.

3. Азимов П.Я., Штукенберг А.Г. Расчёт фазовых диаграмм водно-солевых систем с твёрдыми растворами // Журн. неорг. химии. 2000. Т. 45. С. 1424-1432.

4. Азимов П.Я., Штукенберг А.Г. Термодинамический анализ факторов, определяющих ростовую зональность в метаморфических гранатах // Доклады РАН. 2001. Т. 376. № 2. С. 235-237.

5. Азимов П.Я., Бушмин С.А. Растворимость минералов метаморфических и метасоматических пород в гидротермальных растворах разной кислотности: Термодинамическое моделирование при 400-800°С и 1-5 кбар // Геохимия. 2007а. № 12. С. 1305-1330.

6. Ананьев С.А., Коноваленко С.И. Влияние внутренних напряжений на оптические свойства еремеевита // АН СССР, Сибирское отделение. Геология и геофизика. 1984. № 9. С. 97-103.

7. Аполлонов В.Н. Механизм и условия роста ритмичнопостроенных кристаллов // Доклады РАН. 1999. Т. 364. № 1. С. 94-96.

8. Баларев X. Коэффициент активности в насыщенных тройных водно-солевых растворах, подчиняющихся правилу аддитивности // Журн. неорг. химии. 1987. Т. 32. С. 3068-3074.

9. Барабанов В.Ф. Минералогия вольфрамитовых месторождений Забайкалья. Т. 2. Л.: Изд-во ЛГУ. 1975. 360 с. '12.15

Информация о работе

Штукенберг, Александр Григорьевич
доктора геолого-минералогических наук
Санкт-Петербург, 2009
ВАК 25.00.05

Диссертация

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

Автореферат

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

Похожие работы