English

A Stable Finite-Volume Method for Scalar-Field Dark Matter

Cosmology and Nongalactic Astrophysics 2020-02-11 v2

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

R2 v1 2026-06-23T05:14:43.050Z