Лекция 6 Quantum-Espresso.. Экзамен Удовлетворительно 1. Тест; 2. Защита реферата Хорошо 1. Тест; 2. GAMESS (или Espresso) Отлично 1. Тест; 2. Espresso;

Презентация:



Advertisements
Похожие презентации
Внутренняя структура веществ: Кристаллические решетки, Решетки Браве Соколов Алексей Гр
Advertisements

Исследование систем уравнений, содержащих параметр.
ОБОЗНАЧЕНИЕ НАПРАВЛЕНИЙ В КРИСТАЛЛЕ В кристаллографии возникает необходимость в определении направления отдельных атомных рядов; или атомных плоскостей.
n=1 n=2 n=3 n=4 n=5 E r -- Решением данного уравнения является функция Блоха. Состояние с энергией может быть описано не только функцией, но и. В одномерном.
ОБОЗНАЧЕНИЕ НАПРАВЛЕНИЙ В КРИСТАЛЛЕ В кристаллографии возникает необходимость в определении направления отдельных атомных рядов; или атомных плоскостей.
Влияние прижимающего потенциала на устойчивость электронного кристалла над поверхностью жидкого гелия. В.В.Славин, A.A.Кривчиков.
ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ ПО КУРСУ ФИЗИЧЕСКИЕ ОСНОВЫ ЭКСПЕРИМЕНТАЛЬНЫХ МЕТОДОВ ИССЛЕДОВАНИЯ РЕАЛЬНОЙ СТРУКТУРЫ КРИСТАЛЛОВ Д.ф.-м.н., проф. Э.В.Суворов.
Модели – уравнения квантовой механики. Модели – уравнения квантовой механики. Методы численного исследования: метод функционала плотности, метод Хартри-Фока.
Структурный фактор Повторение материала предыдущей лекции.
Основные свойства синхротронного излучения Синхротронное излучение (СИ) это магнитотормозное излучение релятивистских электронов с энергией где Е – энергия.
БЕСПОРЯДОК В ТВЕРДЫХ ТЕЛАХ Выполнил: Митруков Максим, 553 гр.
Лекция 3 Кинематический анализ рычажных механизмов Задачей кинематического анализа рычажных механизмов является определение кинематических параметров и.
Дифракция медленных электронов Энергии – эВ. Образцы – монокристаллы Глубина снятия информации – один моноатомный слой.
Триангуляция Делоне Выполнил: Е.И. Мишкин Научный руководитель: Пузанкова А.Б.
Одноэлектронное приближение. Атом водорода Почему электрон не падает на ядро? Почему спектры поглощения и излучения атомов и молекул имеют полосчатый.
Атом водорода по Бору. Квантовые постулаты. © В.Е. Фрадкин, А.М.Иконников, 2004.
Тема 2 СТРОЕНИЕ АТОМА. ПЕРИОДИЧЕСКИЙ ЗАКОН И ПЕРИОДИЧЕСКАЯ СИСТЕМА Д.И. МЕНДЕЛЕЕВА (в лекциях использованы материалы преподавателей химического факультета.
4.3.Б. Метод валентных связей Молекула водорода Первый - кинетическая энергия электронов Волновая функция объединенной системы Второй – кулоновское взаимодействие.
Технологическое моделирование (TCAD) Лабораторная работа 1.
Структурная биоинформатика Лекция 1, первая часть Электронная плотность Кристаллы С.А.Спирин 1 октября 2013.
Транксрипт:

Лекция 6 Quantum-Espresso.

Экзамен Удовлетворительно 1. Тест; 2. Защита реферата Хорошо 1. Тест; 2. GAMESS (или Espresso) Отлично 1. Тест; 2. Espresso; 3. GAMESS

Алгоритм Scf (самосогласованный расчет) Оптимизация геометрии Scf Расчет характеристик

How is a crystal structure defined in QE ?

Электронная структура кристалла MgO calculation = 'scf', prefix = mgo, SYSTEM ibrav = 2, celldm(1) = 7.96, nat = 2, ntyp = 2, ecutwfc = 30, ecutrho = 240

celldm(1) – значение равно параметру решетки элементарной ячейки кристалла, выраженному в борах (атомная единица длины). Параметр ячейки для кристалла MgO равен 4.21 ангстрем (Å). 1 бор = радиусу первой боровской орбиты в атоме водорода ( Å). Таким образом, параметр элементарной ячейки MgO равен 7.96 бор. Массив celldm может содержать до 6 величин, которые нужны для задания параметров решетки в самом общем случае низкосимметричного кристалла. Для решетки Браве MgO (ibrav=2) требуется только одно значение. В QE существует и другой способ задать параметры элементарной ячейки, который соответствует традиционному кристаллографическому описанию с использованием параметров элементарной ячейки a, b, c и углов между ними. В этом случае значения celldm(i) задаются равными, соответственно A, B, C, cosAB, cosAC, cosBC. Здесь A, B, C – длины ребер элементарной ячейки в ангстремах, за которыми следуют косинусы указанных углов.

Рис. 1. Сходимость полной энергии E tot кристалла MgO в зависимости от энергии обрезки плоских волн E cut для двух наборов псевдопотенциалов: NC – 1-й набор, US – 2-й набор.

Рис. 2. Сходимость полной энергии E tot кристалла MgO в зависимости от числа k-точек. Расчеты с ультрамягким псевдопотенциалом для атома О и энергией обрезки 30 Ry.

Расчет зонной структуры Рис. 3. Зонный спектр MgO.

Расчет и визуализация распределения электронной плотности Рис. 7. Распределение электронной плотности в MgO.

XCrySDen Charge density difference for the horizontal (left) and tilted (right) N 2 O moieties adsorbed on Pd(110) surface. The charge flows from blue to red regions. Ref: A. Kokalj et al., J. Phys. Chem. B 107, (2003).

Содержание Введение3 1. Возможности пакета Quantum ESPRESSO4 2. Установка Quantum ESPRESSO5 3. Графический интерфейс PWgui6 4. Электронная структура кристалла MgO Задание параметров самосогласованного расчета Контроль сходимости результатов расчета полной энергии Расчет зонной структуры Расчет электронной плотности состояний Расчет и визуализация распределения электронной плотности27 5. Расчет оптических функций30 6. Генерация псевдопотенциалов37 Приложение45 A. Основные команды Linux45 B. Виртуальная машина49

Литература E-Tutorial/