Search

Wave (fuzzy) dark matter ($\psi$DM) consists of ultralight bosons, featuring a solitonic core within a granular halo. Here we extend $\psi$DM to two components, with distinct particle masses $m$ and coupled only through gravity, and investigate the resulting soliton-halo structure via cosmological simulations. Specifically, we assume $\psi$DM contains $75$ per cent major component and $25$ per cent minor component, fix the major-component particle mass to $m_{\rm major}=1\times10^{-22}\,{\rm eV}$, and explore two different minor-component particle masses with $m_{\rm major}:m_{\rm minor}=3:1$ and $1:3$, respectively. For $m_{\rm major}:m_{\rm minor}=3:1$, we find that (i) the major- and minor-component solitons coexist, have comparable masses, and are roughly concentric. (ii) The soliton peak density is significantly lower than the single-component counterpart, leading to a smoother soliton-to-halo transition and rotation curve. (iii) The combined soliton mass of both components follows the same single-component core-halo mass relation. In dramatic contrast, for $m_{\rm major}:m_{\rm minor}=1:3$, a minor-component soliton cannot form with the presence of a stable major-component soliton; the total density profile, for both halo and soliton, is thus dominated by the major component and closely follows the single-component case. To support this finding, we propose a toy model illustrating that it is difficult to form a soliton in a hot environment associated with a deep gravitational potential. The work demonstrates the extra flexibility added to the multi-component $\psi$DM model can resolve observational tensions over the single-component model while retaining its key features., Comment: 19 pages, 24 figures, 1 table, accepted for publication in MNRAS

Ultralight axion-like particles $m_a \sim 10^{-22}$ eV, or Fuzzy Dark Matter (FDM), behave comparably to cold dark matter (CDM) on cosmological scales and exhibit a kpc-size de Broglie wavelength capable of alleviating established (sub-)galactic-scale problems of CDM. Substructures inside an FDM halo incur gravitational potential perturbations, resulting in stellar heating sufficient to account for the Galactic disc thickening over a Hubble time, as first demonstrated by Church et al. We present a more sophisticated treatment that incorporates the full baryon and dark matter distributions of the Milky Way and adopts stellar disc kinematics inferred from recent Gaia, APOGEE, and LAMOST surveys. Ubiquitous density granulation and subhalo passages respectively drive inner disc thickening and flaring of the outer disc, resulting in an observationally consistent `U-shaped' disc vertical velocity dispersion profile with the global minimum located near the solar radius. The observed age-velocity dispersion relation in the solar vicinity can be explained by the FDM-substructure-induced heating and places an exclusion bound $m_a \gtrsim 0.4\times10^{-22}$ eV. We assess non-trivial uncertainties in the empirical core-halo relation, FDM subhalo mass function and tidal stripping, and stellar heating estimate. The mass range $m_a\simeq 0.5-0.7\times10^{-22}$ eV favoured by the observed thick disc kinematics is in tension with several exclusion bounds inferred from dwarf density profiles, stellar streams, and Milky Way satellite populations, which could be significantly relaxed due to the aforesaid uncertainties. Additionally, strongly anisotropic heating could help explain the formation of ultra-thin disc galaxies., Comment: 19 pages, 9 figures; Accepted for publication in MNRAS. Matches published version

We present a new special relativistic hydrodynamics (SRHD) code capable of handling coexisting ultra-relativistically hot and non-relativistically cold gases. We achieve this by designing a new algorithm for conversion between primitive and conserved variables in the SRHD solver, which incorporates a realistic ideal-gas equation of state covering both the relativistic and non-relativistic regimes. The code can handle problems involving a Lorentz factor as high as 106 and optimally avoid the catastrophic cancellation. In addition, we have integrated this new SRHD solver into the code gamer (https://github.com/gamer-project/gamer) to support adaptive mesh refinement and hybrid OpenMP/MPI/GPU parallelization. It achieves a peak performance of 7 × 107 cell updates per second on a single Tesla P100 GPU and scales well to 2048 GPUs. We apply this code to two interesting astrophysical applications: (a) an asymmetric explosion source on the relativistic blast wave and (b) the flow acceleration and limb brightening of relativistic jets.

The Fanaroff-Riley class II radio galaxy Cygnus A hosts jets which produce radio emission, X-ray cavities, cocoon shocks, and X-ray hotspots where the jet interacts with the ICM. Surrounding one hotspot is a peculiar "hole" feature which appears as a deficit in X-ray emission. We use relativistic hydrodynamic simulations of a collimated jet interacting with an inclined interface between lobe and cluster plasma to model the basic processes which may lead to such a feature. We find that the jet reflects off of the interface into a broad, turbulent flow back out into the lobe, which is dominated by gas stripped from the interface at first and from the intracluster medium itself at later times. We produce simple models of X-ray emission from the ICM, the hotspot, and the reflected jet to show that a hole of emission surrounding the hotspot as seen in Cygnus A may be produced by Doppler de-boosting of the emission from the reflected jet as seen by an observer with a sight line nearly along the axis of the outgoing material., Comment: 16 pages, 11 figures, submitted to the journal Galaxies

Fuzzy dark matter (FDM) has been a promising alternative to standard cold dark matter. The model consists of ultralight bosons with mass $m_b \sim 10^{-22}$ eV and features a quantum-pressure-supported solitonic core that oscillates. In this work, we show that the soliton density oscillations persist even after significant tidal stripping of the outer halo. We report two intrinsic yet distinct timescales associated, respectively, with the ground-state soliton wavefunction $\tau_{00}$ and the soliton density oscillations $\tau_\text{soliton}$, obeying $\tau_\text{soliton} /\tau_{00} \simeq 2.3$. The central star cluster (SC) in Eridanus II has a characteristic timescale $\tau_\text{soliton} / \tau_\text{SC} \sim 2$ to $3$ that deviates substantially from unity. As a result, we demonstrate, both analytically and numerically with three-dimensional self-consistent FDM simulations, that the gravitational heating of the SC owing to soliton density oscillations is negligible irrespective of $m_b$. We also show that the subhalo mass function to form Eridanus II does not place a strong constraint on $m_b$. These results are contrary to the previous findings by Marsh & Niemeyer (2019)., Comment: 19 pages, 19 figures; Accepted for publication in PRD

Extremely light bosonic wave dark matter ($\psi$DM) is an emerging dark matter candidate contesting the conventional cold dark matter paradigm and a model subject to intense scrutiny of late. This work for the first time reports testable salient features pertinent to gravitational lenses of $\psi$DM halos. $\psi$DM halos are distinctly filled with large-amplitude, small-scale density fluctuations with $\delta\rho/\rho_{\rm halo}\sim 1$ in form of density granules. This halo yields ubiquitous flux ratio anomalies of a few tens of percent, as is typically found for lensed quasars, and may also produce rare hexad and octad images, for sources located in well-defined caustic zones. We have found new critical features appearing in the highly de-magnified lens center when the halo has sufficiently high surface density near a very compact massive core., Comment: 6 pages, 4 figures

Wave dark matter ($\psi$DM) predicts a compact soliton core and a granular halo in every galaxy. This work presents the first simulation study of an elliptical galaxy by including both stars and $\psi$DM, focusing on the systematic changes of the central soliton and halo granules. With the addition of stars in the inner halo, we find the soliton core consistently becomes more prominent by absorbing mass from the host halo than that without stars, and the halo granules become "non-isothermal", "hotter" in the inner halo and "cooler" in the outer halo, as opposed to the isothermal halo in pure $\psi$DM cosmological simulations. Moreover, the composite (star+$\psi$DM) mass density is found to follow a $r^{-2}$ isothermal profile near the half-light radius in most cases. Most striking is the velocity dispersion of halo stars that increases rapidly toward the galactic center by a factor of at least 2 inside the half-light radius caused by the deepened soliton gravitational potential, a result that compares favorably with observations of elliptical galaxies and bulges in spiral galaxies. However in some rare situations we find a phase segregation turning a compact distribution of stars into two distinct populations with high and very low velocity dispersions; while the high-velocity component mostly resides in the halo, the very low-velocity component is bound to the interior of the soliton core, resembling stars in faint dwarf spheroidal galaxies., Comment: 15 pages, 11 figures

Wave dark matter ($\psi \rm{DM}$) composed of extremely light bosons ($m_{\psi} \sim 10^{-22}\,\rm eV$), with quantum pressure suppressing structures below a kpc-scale de Broglie wavelength, has become a viable dark matter candidate. Compared to the conventional free-particle $\psi {\rm DM}$ (${\rm FP} \psi {\rm DM}$), the extreme-axion $\psi \rm{DM}$ model (${\rm EA} \psi {\rm DM}$) proposed by Zhang & Chiueh (2017) features a larger cut-off wavenumber and a broad spectral bump in the matter transfer function. Here we conduct cosmological simulations to compare the halo abundances and assembly histories at $z=4-11$ between three different scenarios: ${\rm FP} \psi {\rm DM}$, ${\rm EA} \psi {\rm DM}$, and cold dark matter (CDM). We show that ${\rm EA} \psi {\rm DM}$ produces significantly more abundant low-mass haloes than ${\rm FP} \psi {\rm DM}$ with the same $m_{\psi}$, and therefore could alleviate the tension in $m_{\psi}$ required by the Ly$\alpha$ forest data and by the kpc-scale dwarf galaxy cores. We also find that, compared to the CDM counterparts, massive ${\rm EA} \psi {\rm DM}$ haloes are on average $3-4$ times more massive at $z=10-11$ due to their earlier formation, undergo a slower mass accretion at $7 \lesssim z \lesssim 11$, and then show a rapidly rising major merger rate exceeding CDM by $\sim 50\%$ at $4 \lesssim z \lesssim 7$. This fact suggests that ${\rm EA} \psi {\rm DM}$ haloes may exhibit more prominent starbursts at $z \lesssim 7$., Comment: 5 pages, 4 figures, only minor updates (accepted for publication in MNRAS letters)

Although still under debate, observations generally suggest that dwarf spheroidal (dSph) galaxies exhibit large constant-density cores in the centers, which can hardly be explained by dissipationless cold dark matter simulations without baryonic feedback. Wave dark matter (${\psi {\rm DM}}$), characterized by a single parameter, the dark matter particle mass $m_{\psi}$, predicts a central soliton core in every galaxy arising from quantum pressure against gravity. Here we apply Jeans analysis assuming a soliton core profile to the kinematic data of eight classical dSphs so as to constrain $m_{\psi}$, and obtain $m_{\psi}=1.18_{-0.24}^{+0.28}\times10^{-22}{\,\rm eV}$ and $m_{\psi}=1.79_{-0.33}^{+0.35}\times10^{-22}{\,\rm eV}~(2\sigma)$ using the observational data sets of Walker et al. (2007) and Walker et al. (2009b), respectively. We show that the estimate of $m_{\psi}$ is sensitive to the dSphs kinematic data sets and is robust to various models of stellar density profile. We also consider multiple stellar subpopulations in dSphs and find consistent results. This mass range of $m_{\psi}$ is in good agreement with other independent estimates, such as the high-redshift luminosity functions, the reionization history, and the Thomson optical depth to the cosmic microwave background., Comment: 12 pages, 11 figures, 2 tables, add appendix B about the joint analysis, submitted to MNRAS

Simulations of ultralight, $\sim 10^{-22}\,\rm eV$, bosonic dark matter exhibit rich wave-like structure, including a soliton core within a surrounding halo that continuously self-interferes on the de Broglie scale. We show here that as an inherent consequence, the soliton undergoes a confined random walk at the base of the halo potential. This is significant for the fate of the ancient central star cluster in Eridanus II, as the agitated soliton gravitationally shakes the star cluster in and out of the soliton on a time scale of $\sim 100\,\rm Myr$, so complete tidal disruption of the star cluster can occur within $\sim 1\,\rm Gyr$. This destructive effect can be mitigated by tidal stripping of the halo of Eridanus II, thereby reducing the agitation, depending on its orbit around the Milky Way. Our simulations show the Milky Way tide affects the halo much more than the soliton, so the star cluster in Eridanus II can survive for over $5\,\rm Gyr$ within the soliton if it formed after significant halo stripping., Comment: Accepted for publication in Phys. Rev. Lett. 6 pages, 4 figures

The fuzzy dark matter model (FDM, also known as quantum wave dark matter model) argues that light bosons with a mass of $\sim10^{-22}{\;{\rm eV}}$ are a possible candidate for dark matter in the Universe. One of the most important predictions of FDM is the formation of a soliton core instead of a density cusp at the center of galaxies. If FDM is the correct theory of dark matter, then the predicted soliton core can help to form the Central Molecular Zone (CMZ) in the Milky Way. We present high-resolution hydrodynamical simulations of gas flow patterns to constrain the properties of the soliton core based on a realistic Milky Way potential. We find that a dense center is required to form a reasonable CMZ. The size and kinematics of the CMZ offer a relatively strong constraint on the inner enclosed mass profile of the Galaxy. If a soliton core is not considered, a compact nuclear bulge alone with a radially varying mass-to-light ratio can match the observed size and kinematics of the CMZ. A soliton core model with a mass of $\approx4.0\times10^8{\; {\rm M}_{\odot}}$ and a core radius of $\approx0.05{\;{\rm kpc}}$, together with a less massive nuclear bulge with a constant mass-to-light ratio, also agrees nicely with the current data. Such a FDM soliton core corresponds to a boson mass of $\sim2-7\times10^{-22}{\;{\rm eV}}$, which could be further constrained by the improved determination of the mass-to-light ratio in the Galactic center., Accepted for publication in ApJ. 14 pages, 7 figures

Wave dark matter ($\psi$DM), which satisfies the Schr\"odinger-Poisson equation, has recently attracted substantial attention as a possible dark matter candidate. Numerical simulations have in the past provided a powerful tool to explore this new territory of possibility. Despite their successes to reveal several key features of $\psi$DM, further progress in simulations is limited, in that cosmological simulations so far can only address formation of halos below $\sim 2\times 10^{11} M_\odot$ and substantially more massive halos have become computationally very challenging to obtain. For this reason, the present work adopts a different approach in assessing massive halos by constructing wave-halo solutions directly from the wave distribution function. This approach bears certain similarity with the analytical construction of particle-halo (cold dark matter model). Instead of many collisionless particles, one deals with one single wave that has many non-interacting eigenstates. The key ingredient in the wave-halo construction is the distribution function of the wave power, and we use several halos produced by structure formation simulations as templates to determine the wave distribution function. Among different models, we find the fermionic King model presents the best fits and we use it for our wave-halo construction. We have devised an iteration method for constructing the nonlinear halo, and demonstrate its stability by three-dimensional simulations. A Milky-Way-sized halo has also been constructed, and the inner halo is found flatter than the NFW profile. These wave-halos have small-scale interferences both in space and time producing time-dependent granules. While the spatial scale of granules varies little, the correlation time is found to increase with radius by one order of magnitude across the halo., Comment: 14 pages, 7 figures, accepted by PRD

We present GAMER-2, a GPU-accelerated adaptive mesh refinement (AMR) code for astrophysics. It provides a rich set of features, including adaptive time-stepping, several hydrodynamic schemes, magnetohydrodynamics, self-gravity, particles, star formation, chemistry and radiative processes with GRACKLE, data analysis with yt, and memory pool for efficient object allocation. GAMER-2 is fully bitwise reproducible. For the performance optimization, it adopts hybrid OpenMP/MPI/GPU parallelization and utilizes overlapping CPU computation, GPU computation, and CPU-GPU communication. Load balancing is achieved using a Hilbert space-filling curve on a level-by-level basis without the need to duplicate the entire AMR hierarchy on each MPI process. To provide convincing demonstrations of the accuracy and performance of GAMER-2, we directly compare with Enzo on isolated disk galaxy simulations and with FLASH on galaxy cluster merger simulations. We show that the physical results obtained by different codes are in very good agreement, and GAMER-2 outperforms Enzo and FLASH by nearly one and two orders of magnitude, respectively, on the Blue Waters supercomputers using $1-256$ nodes. More importantly, GAMER-2 exhibits similar or even better parallel scalability compared to the other two codes. We also demonstrate good weak and strong scaling using up to 4096 GPUs and 65,536 CPU cores, and achieve a uniform resolution as high as $10{,}240^3$ cells. Furthermore, GAMER-2 can be adopted as an AMR+GPUs framework and has been extensively used for the wave dark matter ($\psi$DM) simulations. GAMER-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome., Comment: 29 pages, 22 figures. MNRAS accepted. GAMER can be found at https://github.com/gamer-project/gamer

We present the implementation and performance of a class of directionally unsplit Riemann-solver-based hydrodynamic schemes on Graphic Processing Units (GPU). These schemes, including the MUSCL-Hancock method, a variant of the MUSCL-Hancock method, and the corner-transport-upwind method, are embedded into the adaptive-mesh-refinement (AMR) code GAMER. Furthermore, a hybrid MPI/OpenMP model is investigated, which enables the full exploitation of the computing power in a heterogeneous CPU/GPU cluster and significantly improves the overall performance. Performance benchmarks are conducted on the Dirac GPU cluster at NERSC/LBNL using up to 32 Tesla C2050 GPUs. A single GPU achieves speed-ups of 101(25) and 84(22) for uniform-mesh and AMR simulations, respectively, as compared with the performance using one(four) CPU core(s), and the excellent performance persists in multi-GPU tests. In addition, we make a direct comparison between GAMER and the widely-adopted CPU code Athena (Stone et al. 2008) in adiabatic hydrodynamic tests and demonstrate that, with the same accuracy, GAMER is able to achieve two orders of magnitude performance speed-up., Comment: 16 pages, 5 figures, accepted for publication in International Journal of High Performance Computing Applications (IJHPCA)

GAMER, a parallel Graphic-processing-unit-accelerated Adaptive-MEsh-Refinement hydrodynamic code, has been extended to support magnetohydrodynamics (MHD) with both the corner-transport-upwind (CTU) and MUSCL-Hancock schemes and the constraint transport (CT) technique. The divergent preserving operator for adaptive mesh refinement (AMR) has been applied to reinforce the divergence-free constraint on the magnetic field. GAMER-MHD has fully exploited the concurrent executions between the GPU MHD solver and other CPU computation pertinent to AMR. We perform various standard tests to demonstrate that GAMER-MHD is both second-order accurate and robust, producing results as accurate as those given by high-resolution uniform-grid runs. We also explore a new 3D MHD test, where the magnetic field assumes the Arnold-Beltrami-Childress (ABC) configuration, temporarily becomes turbulent with current sheets and finally settles to a lowest-energy equilibrium state. This 3D problem is adopted for the performance test of GAMER-MHD. The single-GPU performance reaches $1.2\times 10^8$ and $5.5\times 10^7$ cell-updates/sec for the single- and double-precision calculations, respectively, on Tesla P100. We also demonstrate a parallel efficiency of $\sim 70\%$ for both weak and strong scaling using $1,024$ XK nodes on the Blue Waters supercomputers., Comment: 29 pages, 22 figures, 1 table. GAMER-MHD will soon become available at https://github.com/gamer-project/gamer

We present the newly developed code, GAMER (GPU-accelerated Adaptive MEsh Refinement code), which has adopted a novel approach to improve the performance of adaptive mesh refinement (AMR) astrophysical simulations by a large factor with the use of the graphic processing unit (GPU). The AMR implementation is based on a hierarchy of grid patches with an oct-tree data structure. We adopt a three-dimensional relaxing TVD scheme for the hydrodynamic solver, and a multi-level relaxation scheme for the Poisson solver. Both solvers have been implemented in GPU, by which hundreds of patches can be advanced in parallel. The computational overhead associated with the data transfer between CPU and GPU is carefully reduced by utilizing the capability of asynchronous memory copies in GPU, and the computing time of the ghost-zone values for each patch is made to diminish by overlapping it with the GPU computations. We demonstrate the accuracy of the code by performing several standard test problems in astrophysics. GAMER is a parallel code that can be run in a multi-GPU cluster system. We measure the performance of the code by performing purely-baryonic cosmological simulations in different hardware implementations, in which detailed timing analyses provide comparison between the computations with and without GPU(s) acceleration. Maximum speed-up factors of 12.19 and 10.47 are demonstrated using 1 GPU with 4096^3 effective resolution and 16 GPUs with 8192^3 effective resolution, respectively.

The newly established luminosity functions of high-z galaxies at 4 . z . 10 can provide a stringent check on dark matter models that aim to explain the core properties of dwarf galaxies. The cores of dwarf spheroidal galaxies are understood to be too large to be accounted for by free streaming of warm dark matter without overly suppressing the formation of such galaxies. Here we demonstrate with cosmological simulations that wave dark matter, ψDM, appropriate for light bosons such as axions, does not suffer this problem, given a boson mass of mψ ≥ 1.2 × 10−22 eV (2σ). In this case, the halo mass function is suppressed below ∼ 10M at a level that is consistent with the high-z luminosity functions, while simultaneously generating the kpc-scale cores in dwarf galaxies arising from the solitonic ground state in ψDM. We demonstrate that the reionization history in this scenario is consistent with the Thomson optical depth recently reported by Planck, assuming a reasonable ionizing photon production rate. We predict that the luminosity function should turn over slowly around an intrinsic UV luminosity of MUV & −16 at z & 4. We also show that for galaxies magnified >10× in the Hubble Frontier Fields, ψDM predicts an order of magnitude fewer detections than cold dark matter at z & 10 down to MUV ∼ −15, allowing us to distinguish between these very different interpretations for the observed coldness of dark matter. Subject headings: cosmology: theory – dark ages, reionization, first stars – dark matter – galaxies: abundances – galaxies: evolution – galaxies: high-redshift

Graphical processing units (GPUs) have become widely used to accelerate a broad range of applications, including computational physics and astrophysics, image/video processing, engineering simulations, and quantum chemistry. This chapter presents results obtained from GPU clusters with previous generations of GPU accelerators, which have no or very limited double-precision support. It provides an astrophysical N-body application for star clusters and galactic nuclei, which is currently the well-tested and heavily used application. The chapter explains exemplary implementations of parallel codes using many GPUs as accelerators, so combining message passing parallelization with many-core parallelization. It discusses their benchmarks using up to 512 Fermi Tesla GPUs in parallel on the Mole-8.5 hardware of the Institute of Process Engineering of the Chinese Academy of Sciences (IPE/CAS) in Beijing, and the Laohu Tesla C1070 cluster of the National Astronomical Observatories of CAS in Beijing and smaller clusters in Germany and in the United States. galactic nuclei; graphics processing units; star clusters

In this paper, we describe the architecture and performance of the GraCCA system, a Graphic-Card Cluster for Astrophysics simulations. It consists of 16 nodes, with each node equipped with 2 modern graphic cards, the NVIDIA GeForce 8800 GTX. This computing cluster provides a theoretical performance of 16.2 TFLOPS. To demonstrate its performance in astrophysics computation, we have implemented a parallel direct N-body simulation program with shared time-step algorithm in this system. Our system achieves a measured performance of 7.1 TFLOPS and a parallel efficiency of 90% for simulating a globular cluster of 1024K particles. In comparing with the GRAPE-6A cluster at RIT (Rochester Institute of Technology), the GraCCA system achieves a more than twice higher measured speed and an even higher performance-per-dollar ratio. Moreover, our system can handle up to 320M particles and can serve as a general-purpose computing cluster for a wide range of astrophysics problems., Comment: Accepted for publication in New Astronomy

Quantum turbulence that exhibits vortex creation, annihilation and interactions is demonstrated as an exact solution of the time-dependent, free-particle Schr\"odinger equation evolved from a smooth random-phased initial condition. Relaxed quantum turbulence in 2D and 3D exhibits universal scaling in the steady-state energy spectrum as k-1 in small scales. Due to the lack of dissipation, no evidence of the Kolmogorov-type forward energy cascade in 3D or the inverse energy cascade in 2D is found, but the rotational and potential flow components do approach equi-partition in the scaling regime. In addition, the 3D vortex line-line correlation exhibits universal behaviour, scaled as \Deltar^-2, where \Deltar is the separation between any two vortex line elements, in fully developed turbulence. We also show that the quantum vortex is not frozen to the matter, nor is the vortex motion induced by other vortices via Biot-Savart's law. Thus, the quantum vortex is actually a nonlinear wave, propagating at a speed very different from a classical vortex., Comment: 9 pages, 14 figures