A Stable Finite-Volume Method for Scalar-Field Dark Matter
Abstract
We describe and test a family of new numerical methods to solve the Schrodinger equation in self-gravitating systems, e.g. Bose-Einstein condensates or 'fuzzy'/ultra-light scalar field dark matter. The methods are finite-volume Godunov schemes with stable, higher-order accurate gradient estimation, based on a generalization of recent mesh-free finite-mass Godunov methods. They couple easily to particle-based N-body gravity solvers (with or without other fluids, e.g. baryons), are numerically stable, and computationally efficient. Different sub-methods allow for manifest conservation of mass, momentum, and energy. We consider a variety of test problems and demonstrate that these can accurately recover solutions and remain stable even in noisy, poorly-resolved systems, with dramatically reduced noise compared to some other proposed implementations (though certain types of discontinuities remain challenging). This is non-trivial because the "quantum pressure" is neither isotropic nor positive-definite and depends on higher-order gradients of the density field. We implement and test the method in the code GIZMO.
Keywords
Cite
@article{arxiv.1811.05583,
title = {A Stable Finite-Volume Method for Scalar-Field Dark Matter},
author = {Philip F. Hopkins},
journal= {arXiv preprint arXiv:1811.05583},
year = {2020}
}
Comments
9 pages, 3 figures. Updated to published (MNRAS) version. Expanded from first draft to present a family of related numerical schemes, with additional discussion of nodes