Quadratic Scaling Bosonic Path Integral Molecular Dynamics

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Feldman, Yotam M. Y., Hirshberg, Barak
Natura: Preprint
Pubblicazione: 2023
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866913691572633600
author Feldman, Yotam M. Y.
Hirshberg, Barak
author_facet Feldman, Yotam M. Y.
Hirshberg, Barak
contents Bosonic exchange symmetry leads to fascinating quantum phenomena, from exciton condensation in quantum materials to the superfluidity of liquid Helium-4. Unfortunately, path integral molecular dynamics (PIMD) simulations of bosons are computationally prohibitive beyond $\mathord{\sim} 100$ particles, due to a cubic scaling with the system size. We present an algorithm that reduces the complexity from cubic to quadratic, allowing the first simulations of thousands of bosons using PIMD. Our method is orders of magnitude faster, with a speedup that scales linearly with the number of particles and the number of imaginary time slices (beads). Simulations that would have otherwise taken decades can now be done in days. In practice, the new algorithm eliminates most of the added computational cost of including bosonic exchange effects, making them almost as accessible as PIMD simulations of distinguishable particles.
format Preprint
id arxiv_https___arxiv_org_abs_2305_18025
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Quadratic Scaling Bosonic Path Integral Molecular Dynamics
Feldman, Yotam M. Y.
Hirshberg, Barak
Computational Physics
Other Condensed Matter
Quantum Gases
Chemical Physics
Bosonic exchange symmetry leads to fascinating quantum phenomena, from exciton condensation in quantum materials to the superfluidity of liquid Helium-4. Unfortunately, path integral molecular dynamics (PIMD) simulations of bosons are computationally prohibitive beyond $\mathord{\sim} 100$ particles, due to a cubic scaling with the system size. We present an algorithm that reduces the complexity from cubic to quadratic, allowing the first simulations of thousands of bosons using PIMD. Our method is orders of magnitude faster, with a speedup that scales linearly with the number of particles and the number of imaginary time slices (beads). Simulations that would have otherwise taken decades can now be done in days. In practice, the new algorithm eliminates most of the added computational cost of including bosonic exchange effects, making them almost as accessible as PIMD simulations of distinguishable particles.
title Quadratic Scaling Bosonic Path Integral Molecular Dynamics
topic Computational Physics
Other Condensed Matter
Quantum Gases
Chemical Physics
url https://arxiv.org/abs/2305.18025