Многомасштабное моделирование процессов и структур в нанотехнологиях

Подождите немного. Документ загружается.

ной сетке с равномерной дискретизацией по осям x и y (с шагом 5нм). Для

ускорения расчет магнитостатического взаимодействия осуществлялся в

пространстве Фурье, для перехода в него и обратно использовалась тех-

ника быстрого преобразования Фурье. Расчёт проводился в режиме релак-

сации при заданном постоянном поле, создаваемым током, текущим пер-

пендикулярно плоскости объекта и пронизывающим

структуру. Возмож-

ности пакета также позволяли учитывать температурные флуктуации в

материалах.

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

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

ный момент, рассчитанный относительно центра объектов. Показано, что

для объектов, находящихся в вихревых состояниях, тороидный момент

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

объ-

ектов из одного вихревого состояния в другое характеризуется изменени-

ем знака тороидного момента.

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

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

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

протекании токов без спиновой поляризации и токов, поляризованных по

спину. В частности, наличие спиновой поляризации

существенно умень-

шает величину критического тока, при котором происходит переключение

объектов. Также исследовано влияние различных направлений спиновой

поляризации на процесс перемагничивания.

Список литературы

1. Yang T., Hara M., Hirohata A., Kimura T., Otani Y. Appl. Phys. Lett. 90, 022504 (2007).

2. Sheka D.D., Gaididei Y., Mertens F.G. Appl. Phys. Lett. 91, 082509 (2007).

3. Slonczewski J.C. J. Magn. Magn. Mater. 159, L1 (1996).

Е.Г. ГРИГОРЬЕВ

Московский инженерно-физический институт (государственный университет)

МОДЕЛИРОВАНИЕ МАКРОСКОПИЧЕСКИХ ПРОЦЕССОВ

В ПОРОШКОВЫХ СРЕДАХ ПРИ МОЩНОМ

ЭЛЕКТРОИМПУЛЬСНОМ ВОЗДЕЙСТВИИ

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

клика порошковых проводящих сред на электроимпульсное воздействие исследо-

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

при электроимпульсном спекании материалов. Выявлен волновой характер про-

цесса консолидации порошкового материала. Установлены соотношения между

характерными временными масштабами процессов электроимпульсного спекания

порошковых материалов.

Актуальной проблемой современного материаловедения является раз-

работка экономически выгодных технологических процессов изготовле-

ния консолидированных нано материалов с регулируемыми и заданными

параметрами состава, структуры и свойств [1]. В настоящее время активно

исследуются методы консолидации, использующие воздействие давления

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

К ним относятся электроразрядное и электроимпульсное

спекание [2],

электроимпульсное прессование [3], развиваемые в странах СНГ, а за ру-

бежом – это Field Assisted Sintering Technique (FAST), Plasma Assisted

Sintering (PAS), Spark Plasma Sintering (SPS) и другие [4]. Главное досто-

инство этих методов заключается в кратковременности процессов консо-

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

Наиболее принципиальным отличием электроимпульсных способов кон-

солидации порошковых материалов является концентрация высокой

плотности выделяемой энергии на контактах между

частицами порошка.

Высокая концентрация энергии в контактных областях вызывает интен-

сивные локальные деформации и фазовые превращения. Физические про-

цессы в контактных зонах характеризуются пространственной неодно-

родностью и нестационарностью во времени. Их изучение и моделирова-

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

частиц позволяет определить требования к оптимальным параметрам про-

цесса электроимпульсного спекания. Наряду с процессами, протекающи-

ми в масштабах отдельных частиц, необходимо установить закономерно-

сти макроскопического поведения материала в масштабе всего спекаемого

изделия. Факторами, определяющими макроскопические процессы, явля-

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

териал, пространственное и временное распределение вводимой энергии в

течение процесса спекания, характер и

величина давления, воздействую-

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

ложена математическая модель физических процессов протекающих при

электроимпульсном спекании в порошковом материале как в масштабе

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

порошка (в зонах контактов частиц). В основе модели лежат законы

со-

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

намики для компактных и порошковых проводящих сред. Макроскопиче-

ское уравнение состояния порошковой среды определяется с учетом про-

цессов в зонах контактов частиц порошка. Изучение процессов консоли-

дации порошковых материалов на основе разработанной модели позволи-

ло выявить основные закономерности электроимпульсного спекания и

установить безразмерные параметры, определяющие характер протекаю-

щих процессов. Результаты расчетов, выполненные на основе предложен-

ной модели, в совокупности с экспериментальными данными выявили

волновой характер процесса уплотнения материала.

Важным результатом данной работы является определение соотноше-

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

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

рошковые среды.

Список литературы

1. . Андриевский Р.А. Наноматериалы: концепция и современные проблемы // Россий-

ский химический журнал, т. 46, №5, с.50 (2002).

2. Райченко А.И. Основы процесса спекания порошков пропусканием электрического

тока// М.: Металлургия, 1987. с. 129.

3. Баланкин С.А., Горбачев Л.П., Григорьев Е.Г. и др. Тепловые процессы при электро-

импульсном

прессовании порошков// Физика и химия обработки материалов, N2, с.124.

(1984).

4. Zhang J., Zavaliangos A., Groza J. FIELD ACTIVATED SINTERING TECHNIQUES: A

COMPARISON AND CONTRAST// P/M Science & Technology Briefs., Vol. 5, No. 4, p.5 (2003)

Е.В. ПУСТОВАЛОВ

, Б.Н. ГРУДИН, В.С. ПЛОТНИКОВ,

Е.Б. МОДИН, О.В. ВОЙТЕНКО

Дальневосточный государственный университет, Владивосток

МОДЕЛИРОВАНИЕ СТРУКТУРЫ

И ЭЛЕКТРОННАЯ ТОМОГРАФИЯ НАНОКЛАСТЕРОВ

В СТРУКТУРЕ БЫСТРОЗАКАЛЕННЫХ АМОРФНЫХ

СПЛАВОВ

В работе исследованы методы визуализации структуры атомных нанокласте-

ров методом электронной микроскопии и томографии. Использование методов

обработки изображений Фурье, косинус и вейвлет фильтрации позволяет визуали-

зировать на модельных высокоразрешающих электронно-микроскопических изо-

бражениях атомные нанокластеры различной структуры и формы. Высокоразре-

шающая электронная томография расчетных электронно-микроскопических изо-

бражений позволяет исследовать геометрию

упорядоченных нанокластеров в

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

рицы из 2500 атомов Fe

и упорядоченные кластеры Fe

B и Fe

B с количест-

вом атомов от 20 до 500. Показана возможность использования косинус и вейвлет-

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

нокластеров.

Атомная структура аморфных твердых тел представляет значительный

научный интерес, поскольку определяет уникальность свойств, сущест-

венно отличающихся от кристаллических материалов. Особое внимание

уделяется не идеальной аморфной структуре, а прошедшей частичную

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

нанокластеров с упорядоченной структурой [1]. Актуальной задачей явля-

ется расшифровка структуры нанокластеров с использованием методов

высокоразрешающей электронной микроскопии (ВРЭМ) и электронной

томографии (ЭТ).

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

стер из 2500 атомов состава Fe

, полученный методом Ишикава [2] с

внедрением атомов металлоида и последующей релаксацией. Модели на-

нокластеров представляют собой искаженные нанокристаллы Fe

B и Fe

состоящие из тригональных призм в центре атом B и по узлам атомы Fe

[3]. Смешанные модели представляли собой ячейки аморфной матрицы

3*3*3 нм с внедренными нанокластерами различной структуры. Расчет

электронно-микроскопических изображений осуществлялся слоевым ме-

тодом с заданием характеристик конкретных электронных микроскопов.

Для моделирования электронной томографии высокого разрешения осу-

ществлялся поворот ячейки вокруг

произвольной оси на углы ±70º с ша-

гом 1º и расчет серии электронно-микроскопических изображений. Для

создания томограммы использовался пакет IMOD [4]. Для обработки и

парметризаци ВРЭМ изображений использовались ортогональные преоб-

разования: Фурье, Хартли, косинус и непрерывное вейвлет-преобразо-

вание [5,6].

В результате моделирования получено, что обнаружение и расшиф-

ровка структуры атомных кластеров с упорядоченной структурой

аморфной матрице возможно при соотношении толщин аморф-

ный/упорядоченный 1:1 или при большей толщине упорядоченной струк-

туры.

Показано, что для фильтрации и улучшения на ее основе психофизи-

ческого восприятия изображений предпочтительнее использовать коси-

нусное преобразование. Рассмотрены возможности реализации согласо-

ванной фильтрации на основе косинусного преобразования. Предложена и

исследована процедура декомпозиции

ВРЭМ изображений с использова-

нием различных ортогональных преобразований и фильтров, построенных

на основе фурье-образов вейвлетов Гаусса и Морле. Показана перспек-

тивность анализа структуры на ВРЭМ изображениях по скалограммам

ортогональных преобразований изображений [7].

Список литературы

1. Grudin B., Dolzhikov S., Plotnokov V., Pustovalov E. High resolution electron microscopy

of amorphous metallic alloys nanostructure // Proc.of Asian Symp. Adv.Mater., Vladivostok, 1-4

Oct. 2007, p.O-14.

2. Ichikawa T. The assembly of hard spheres as a structure model of amorphous iron. //

Phys.Stat.Sol., v.29a, N1, p.293 (1975).

3. Gaskell P.H. Medium range order and random networks. // J.of Non-Cryst.Sol., v.293-295,

p.146 (2001).

4. imod

5. Грудин Б.Н., Кисленок Е.Г., Плотников В.С., Фищенко В.К. Анализ, фильтрация и де-

композиция микроскопических изображений на основе ортогональных преобразований. //

Автометрия, т.43, с. 24 (2007).

6. Грудин Б.Н., Должиков С.В., Кисленок Е.Г., Плотников В.

С. Методы анализа струк-

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

// Известия РАН. Серия Физическая, т.71, №10, с.1486 (2007).

7. Grudin B., Dolzhikov S., Plotnokov V., Pustovalov E. Electron microscopy investigation of

amorphous alloys // Abstr. of Electron microscopy and multiscale modeling conf. Moscow, 3-

7Sept. p.38.

В.А. ДАНИЛОВ

, Т.А. ПРОКОПЬЕВА, С.С. КАНТОРОВИЧ

Уральский государственный университет им. А.М. Горького, Екатеринбург

МИКРОСТРУКТУРА НАНОСЛОЯ МАГНИТНОЙ

ЖИДКОСТИ В ОСНОВНОМ СОСТОЯНИИ

В работе исследуется микроструктура нанослоя магнитной жидкости при 0 К.

Для анализа топологии образующихся агрегатов используется компьютерный

эксперимент и теоретическое моделирование. Исследованы наиболее вероятные

агрегаты из наночастиц: кольца, цепочки и их комбинации. По результатам моде-

лирования идеальное кольцо из феррочастиц является наиболее вероятной тополо-

гической структурой агрегата для феррожидкости в нанослое в

основном состоя-

нии.

Синтезированные в середине 60-х годов 20-ого столетия феррожидко-

сти (магнитные жидкости, ферроколлоиды) представляют собой взвеси

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

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

вышают 20-30 нм, и каждая частица, покрытая защитной оболочкой, об-

ладает собственным магнитным моментом. Не только теоретически, но и

экспериментально микроструктура магнитных

жидкостей остается до

конца не изученной. Из-за того, что феррочастицы обладают столь малы-

ми размерами, непосредственное наблюдение микроструктур - технически

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

костей оптически непрозрачны. В 2003 году в лаборатории Вант Гоффа

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

рожидкостей – нанослоев (толщина слоя около 40 нм) [1], анализ структу-

ры которых может проводиться электронной микроскопией. Для эффек-

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

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

ется анализ структур в основном состоянии.

В данной работе

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

тиц методом Монте Карло (МК). Последовательно понижая температуру,

изучается процесс рекомбинации частиц. На данном этапе было выявлено

несколько наиболее часто встречающихся топологических структур агре-

гатов: идеальное кольцо, два соосных кольца, цепочка и два кольца, рас-

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

с работой [2].

Теоретический расчет полной энергии различных наност-

руктур показал, что наиболее энергетически выгодным является одиноч-

ное идеальное кольцо из феррочастиц (рис. 1). Также на основе компью-

терного эксперимента и теоретического анализа было показано, что близ-

кой энергией обладают два соосных кольца, если их диаметры макси-

мально близки (см. рис. 2). Оказалось, что образование двух

соосных ко-

лец является «необратимым», так как процесс рекомбинации наталкивает-

ся на высокий энергетический барьер, который возникает, когда количе-

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

ческой модели и компьютерного эксперимента совпадают с высокой точ-

ностью.

Список литературы

1. Klokkenburg M, Dullens RPA, Kegel WK, Erné BH, Philipse A. Quantitative Real-Space

Analysis of Self-Assembled Structures of Magnetic Dipolar colloids // Phys. Rev. Lett., V. 96, P.

037203, (2006).

Рис. 1. Идеальное кольцо. Снэпшот

из МК симуляций

Рис. 2. Два соосных кольца.

Снэпшот из МК симуляций

2. Morimoto H, Maekawa T, Matsumoto Y. Statistical analysis of 2-dimensional cluster

structures composed of ferromagnetic particles based on a flexible chain model // Phys. Rev. E, V.

68, P. 061505, (2003).

А.Л. ДАНИЛЮК

, А.В. АНДРЕЕНКО, Е.Л. ПРУДНИКОВА,

Б.Г. ШУЛИЦКИЙ, К.И. ЯНУШКЕВИЧ

Белорусский государственный университет информатики и радиоэлектроники,

Минск

Объединенный институт физики твердого тела и полупроводников НАН

Беларуси, Минск

МОДЕЛИРОВАНИЕ ВИХРЕВОГО СОСТОЯНИЯ И

МАГНИТНЫХ ПЕРЕХОДОВ НАНОДИСКОВ ЦЕМЕНТИТА

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

нанодисков цементита с помощью пакета OOMMF. Проведены расчеты распреде-

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

области температур 300-470 K. Выявлены неустойчивые области магнитных со-

стояний. Определены размеры вихря в зависимости от радиуса, толщины наноди-

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

разра-

ботки элементов магнитной памяти.

Наночастица в вихревом магнитном состоянии является перспектив-

ным кандидатом для бита информации в магнитных накопителях повы-

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

датчиков [1]. Вихревое магнитное состояние наночастиц является предме-

том многочисленных теоретических и экспериментальных исследований,

но остаются еще нерешенными вопросы, связанные с моделированием

вихревой намагниченности, особенно при повышенных температурах.

Целью

данной работы является моделирование магнитных переходов и

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

пературном диапазоне 300-470 K.

Для преодоления суперпарамагнитного порога в магнитной записи по-

требуются стабильные к тепловым флуктуациям ферромагнитные нано-

размерные материалы. Для этих целей необходимы высококоэрцитивные

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

свойств в условиях быстрого нагрева и

охлаждения. Наночастицы цемен-

тита, капсулированные в углеродных нанотрубках (УНТ), проявляют вы-

сокую стабильность физико-химических свойств. Их коэрцитивная сила

~0.4 кЭ при температуре 300 K (точка Кюри 480-483 К). При нагреве до

температуры блокировки 420-450 K коэрцитивная сила спадает практиче-

ски до нуля [2]. Проведенные нами экспериментальные исследования по-

казали, что нагрев и охлаждение УНТ, содержащих цементит, позволяет

сохранить обратимость магнитных свойств цементита в области темпера-

тур от 78 до 720 K. Установлено, что рост температуры ведет к уменьше-

нию намагниченности согласно закону Блоха в области температур от 78

до 400 K. В области 300-475 K температурная зависимость намагниченно-

сти описывается моделью Стонера.

Микромагнитное моделирование нанодисков цементита проводили

помощью пакета OOMMF (Object Oriented Micromagnetic Framework),

разработанного национальным институтом стандартов и технологий

(США) [3]. Код программы реализован на основе решения системы диф-

ференциальных уравнений Ландау-Лифшица-Гильберта для спиновой

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

ления магнитостатического (размагничивающего) поля.

Проведены расчеты намагниченности и магнитных переходов в нано-

дисках цементита с характерными размерами

от 5 до 100 нм. Использова-

лись следующие параметры для цементита: намагниченность насыщения

=(1.2-1.5)⋅10

А/м (300 K), обменная константа A=(0.42-0.49)⋅10

-11

Дж/м.

Определены характерные размеры нанодиска, при которых его маг-

нитное состояние является устойчивым однородным

z- состоянием, пла-

нарным и вихревым состоянием. Рассчитана диаграмма основных состоя-

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

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

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

состояний. Показано, что вихревое начинается при диаметре нанодиска 20

нм и его толщине 10 нм. Показано, что с ростом

толщины нанодиска на-

блюдается линейный рост диаметра вихря, переходящей в закономерность

1/3.4

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

переходов между вихревым и

z-состоянием в области температур 300-470

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

из вихревого в перпендикулярное

z-состояние из-за роста критического

радиуса нанодиска. Определена область размеров нанодисков цементита,

в которой такой переход происходит в диапазоне 300-470 К.

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

грамма магнитных переходов между вихревым магнитным и однородным

z-состоянием в нанодисках цементита в области температур 300-470 K.

Список литературы

1. Cowburn R.P. // J. Magn. Magn. Mater. 242-245, 505 (2002).

2. Колмогорцев С.В., Исхаков Р.С., Балаев А.Д., и др. Нанонити Fe в углеродных нанот-

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

ФТТ, т. 49, с. 700 (2007).

3. The object oriented micromagnetic framework. http://math.nist.gov/oommf/

А.В. ДЗЯБЧЕНКО

Научно-исследовательский физико-химический институт им. Л.Я.Карпова

ПРЕДСКАЗАНИЕ СТРУКТУРЫ И СТАБИЛЬНОСТИ

ОРГАНИЧЕСКИХ КРИСТАЛЛОВ И НАНОРАЗМЕРНЫХ

СИСТЕМ

Цель данной работы показать, что решение задачи предсказания структуры

кристаллов является важным и абсолютно необходимым этапом в разработке ме-

тодов и программ решения аналогичной задачи для наносистем. В докладе рас-

смотрены основные этапы решения задачи конструирования кристаллов из орга-

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

ция которых существенно зависит

от упаковки.

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

использовать известные структурные и термодинамические свойства в

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

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

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

лов путем проведения соответствующих расчетов и сравнения моделей с

экспериментом

. Уточненные потенциалы с полным правом можно ис-

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

которых наблюдать экспериментальными методами затруднительно или

вообще невозможно. Существующие методы квантовой химии в непо-

средственном виде вряд ли непригодны для глобальной минимизации

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

машинных ресурсов, требуемых на проведение вычислений. Однако

они

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

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

систем. Ландшафт поверхности потенциальной энергии (ППЭ) в глобаль-

ном параметрическом пространстве ответственен не только за потенци-

ально возможные структуры вещества при данных условиях, но и за набор

свойств и явлений, происходящих в твердом теле как

реакция на измене-

ние внешних условий (прежде всего давления и температуры). Для пони-

мания этих свойств важно знание не только минимумов ППЭ, но также и

седловых точек (СТ), разделяющих соседние минимумы, поскольку высо-

та СТ характеризует энергетический барьер фазового перехода или спо-

собность системы к появлению в ней элементов частичного

беспорядка.

Наряду с фундаментальной теоретической значимостью проблемы в

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

тентоспособных полиморфов лекарственных веществ и рентгенострук-

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

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

граммных систем для конструирования молекулярных материалов, кото-

рые включают методы теоретической симуляции структуры, молекуляр-

но-графические программы, базы структурных данных и программы рас-

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

том. В течение последних лет Кембриджский центр кристаллографиче-

ских структурных данных

регулярно организует "слепые тесты" предска-

заний структур кристаллов органических молекул, имеющие цель объек-

тивной оценки программ независимых разработчиков [1].

В докладе сообщаются основные характеристики комплекса программ,

применяемых нами на разных этапах решения задачи структурного про-

гноза. Так, программа PMC [2] осуществляет глобальный поиск опти-

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

онных параметров.

Программа CRYCOM [3] служит для сравнения кри-

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

тификации группы симметрии данной структуры, а также для определе-

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

грамма FitMEP [4] служит для определения характеристик точечных

мультиполей до квадруполя включительно путем аппроксимации непре-

рывных распределений МЭП, рассчитанных ab initio (Gaussian 98 [5]).

Доклад иллюстрируется

примерами успешных предсказаний структур,

в числе которых новая полиморфная форма бензола, полимеры С

и С

кембриджские тестовые структуры, комплекс [Co(NH

)

(NO

)]Cl

а также

примеры определений органических структур методом рентгенографии на

порошке с помощью моделей, найденных энергетическим расчетом.

Работа выполнена в рамках проекта РФФИ 05-03-32808

Список литературы

1. Day, G.M., et al. A third blind test of crystal structure prediction. //Acta crystallogr. v.

B61, 511 (2005).

2. Dzyabchenko, A. V. PMC, Version 2005. // Karpov Institute of Physical Chemistry, Mos-

cow: (2006).

3. Dzyabchenko A.V. Method of crystal-structure similarity searching //Acta Crystallogr. v.

B50, 414 (1994).

4. Дзябченко А.В. Мультипольная аппроксимация электростатического потенциала

можекул // Журнал физической химии, т. 82, вып. 5 (2008).

5. Frisch M. J., Trucks G. W., Schlegel H. B., et al.// Gaussian 98, revision A.11.3, Gaussian,

Inc., Pittsburgh PA, 2002