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

1 Применение суперкомпьютерных технологий и технологии CUDA для PIMC моделирования в физике конденсированного состояния вещества на примере металлического водорода Александр Новоселов (ИТПМ МГУ) Олег Павловский (ИТПМ МГУ)

2 Что такое PIMC: Path Integral теоретически в численном моделировании

3 Что такое PIMC: Path Integral среднее значение наблюдаемой

4 Что такое PIMC: Monte Carlo (по всему бесконечному числу конфигураций) - большое, но конечное число случайных конфигураций с вероятностью

5 Процессор 1... Процессоры...Процессор N proc N conf /N proc x P(x) Генерация N conf /N proc случайных конфигураций x с вероятностью P(x) … A Вычисление и усреднение наблюдаемой A по этим конфигурациям … N proc Усреднение данных по процессорам (среднее арифметическое N proc чисел) PIMC: распараллеливание 1

6 PIMC: генерация конфигураций x 0 Любая начальная конфигурация x 0 Марковский процесс Условие детального баланса Предельное распределение «нетермализованные» конфигурации «термализованные» конфигурации

7 PIMC: распараллеливание 2 Процессор 1... Процессоры...Процессор N proc Начальная конфигурация Генерация нетермализованных конфигурации до достижения предельного распределения «первая термализованная» конфигурация N conf /N proc x P(x) Генерация N conf /N proc случайных конфигураций x с вероятностью P(x) … A Вычисление и усреднение наблюдаемой A по этим конфигурациям … N proc Усреднение данных по процессорам (среднее арифметическое N proc чисел)

8 Алгоритм: проблема автокорреляций

9 Алгоритм: многоуровневый

10 Алгоритм: многоуровневый выбор действия уровня

11 Алгоритм: многоуровневый CUDA-распараллеливание Блок 1 Блок 2 Блок N blocks Координаты разных слоёв не связаны У блоков нет общего кэша Координаты разных частиц одного слоя связаны У нитей в блоке есть общий кэш Нить 1 Нить 2 Нить N threads

12 Металлический водород: на Земле Ячейки с алмазными наковальнями – Diamond Anvil Cells Eremets M.I., Troyan I.A. Conductive dense hydrogen // Nature Materials

13 Металлический водород: в космосе

14 Металлический водород: модель Вигнеровский кристалл протонов

15 Металлический водород: постановка задачи Энергия –Потенциальная –Кинетическая –Полная внутренняя Давление Отношение Линдеманна Конфигурации Уравнения состояния Фазовый переход

16 Результаты: потенциальная энергия r s =200 (ρ=2,1*10^6 кг/м^3)

17 Результаты: кинетическая энергия r s =200 (ρ=2,1*10^6 кг/м^3)

18 Результаты: давление r s =200 (ρ=2,1*10^6 кг/м^3)

19 Результаты: отношение Линдеманна

20 Результаты: конфигурации r s =200 (ρ=2,1*10^6 кг/м^3) T=14*10^3 К

21 Результаты: конфигурации r s =200 (ρ=2,1*10^6 кг/м^3) T=13*10^3 К Объёмно-центрированная кубическая (bcc) решётка

22 Фазовый переход

23 Фазовый переход: проблема термализации «метастабильное состояние» Один процессор успел дотермализоваться

24 Фазовый переход: распараллеливание термализации Процессор 1 Процессор 2 Процессор 3 Полная термализация Процессор 1 Процессор 2 Процессор 3 Полная термализация Раздача термализованной конфигурации время Эффективная набранная статистика

25 Выводы: физика Получены уравнения состояния и – полное термодинамическое описание bccОбнаружен и исследован фазовый переход жидкость – bcc кристалл

26 Выводы: вычислительные технологии Реализованы высокопроизводительные алгоритмы, применимые и в других моделях mpiИспользование суперкомпьютера – почти идеальное mpi-распараллеливание CUDA GPUИспользование технологии CUDA, распараллеливание на GPU – ускорение в 2*30 раз (на процессор) Полный объём использованных вычислительных ресурсов процессоро-часов