Source author record

Hsi-Yu Schive

Hsi-Yu Schive appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2026arXiv

Universal numerical convergence criteria for subhalo tidal evolution

Dark matter subhalos and satellite galaxies in state-of-the-art cosmological simulations still suffer from the ``overmerging'' problem, where inadequate force and/or mass resolution cause artificially enhanced tidal mass loss and premature disruption. Previous idealized simulations addressing this issue have been restricted to a small subset of the subhalo orbital parameter space, and all assumed subhalos to be isotropic. Here, we present the first extensive simulation suite that quantifies numerical convergence in the tidal evolution of anisotropic subhalos under varying numerical resolutions and orbits. We report a universal force resolution criterion: the subhalo's instantaneous tidal radius must always be resolved by at least 20 cells in adaptive mesh refinement (AMR)-based simulations, or by 20 softening lengths (Plummer equivalent) in tree-based simulations, regardless of refinement details or subhalo physical properties such as concentration or velocity anisotropy. We also report a universal expression for the discreteness-noise-driven scatter in the bound-mass fraction of subhalos that depends only on the subhalo mass resolution at infall and the instantaneous bound mass fraction, agnostic of any further subhalo properties. Such stochastic discreteness noise causes both premature disruption and, notably, spurious survival of poorly mass-resolved subhalos. We demonstrate that as many as 50 percent of all subhalos in state-of-the-art cosmological simulations are likely to be either force and/or mass unresolved. Our findings advocate for adaptive softening or grid refinement based on the instantaneous tidal radius of the subhalo.

preprint2020arXiv

Dynamical Evidence of a Solitonic Core of $10^{9}M_\odot$ in the Milky Way

A wavelike solution for the non-relativistic universal dark matter (wave-DM) is rapidly gaining interest, following distinctive predictions of pioneering simulations of cosmic structure as an interference pattern of coherently oscillating bosons. A prominent solitonic standing wave is predicted at the center of every galaxy, representing the ground state, that has been identified with the wide, kpc scale dark cores of common dwarf-spheroidal galaxies, providing a boson mass of, $\simeq 10^{-22}$ eV. A denser soliton is predicted for Milky Way sized galaxies where momentum is higher, so the de Broglie scale of the soliton is smaller, $\simeq 100$ pc, of mass $\simeq 10^9 M_\odot$. Here we show the central motion of bulge stars in the Milky Way implies the presence of such a dark core, where the velocity dispersion rises inversely with radius to a maximum of $\simeq 130$ km/s, corresponding to an excess central mass of $\simeq 1.5\times 10^9 M_\odot$ within $\simeq 100$ pc, favouring a boson mass of $\simeq 10^{-22}$ eV. This quantitative agreement with such a unique and distinctive prediction is therefore strong evidence for a light bosonic solution to the long standing Dark Matter puzzle, such as the axions generic in String Theory.

preprint2020arXiv

Multiple Images and Flux Ratio Anomaly of Fuzzy Gravitational Lenses

Extremely light bosonic wave dark matter ($ψ$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 $ψ$DM halos. $ψ$DM halos are distinctly filled with large-amplitude, small-scale density fluctuations with $δρ/ρ_{\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.

preprint2020arXiv

Soliton Random Walk and the Cluster-Stripping Problem in Ultralight Dark Matter

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.

preprint2020arXiv

Testing the Prediction of Fuzzy Dark Matter Theory in the Milky Way Center

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.

preprint2016arXiv

Contrasting Galaxy Formation from Quantum Wave Dark Matter, $ψ$DM, with $Λ$CDM, using Planck and Hubble Data

The newly established luminosity functions of high-z galaxies at $4 \lesssim z \lesssim 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_ψ \ge 1.2 \times 10^{-22}{\,\rm eV}$ ($2σ$). In this case, the halo mass function is suppressed below $\sim 10^{10}{\,M_\odot}$ 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 $M_{\rm UV} \gtrsim -16$ at $z \gtrsim 4$. We also show that for galaxies magnified $\mathord{>}10\times$ in the Hubble Frontier Fields, $ψ$DM predicts an order of magnitude fewer detections than cold dark matter at $z \gtrsim 10$ down to $M_{\rm UV} \sim -15$, allowing us to distinguish between these very different interpretations for the observed coldness of dark matter.

preprint2014arXiv

Cosmic Structure as the Quantum Interference of a Coherent Dark Wave

The conventional cold, particle interpretation of dark matter (CDM) still lacks laboratory support and struggles with the basic properties of common dwarf galaxies, which have surprisingly uniform central masses and shallow density profiles. In contrast, galaxies predicted by CDM extend to much lower masses, with steeper, singular profiles. This tension motivates cold, wavelike dark matter ($ψ$DM) composed of a non-relativistic Bose-Einstein condensate, so the uncertainty principle counters gravity below a Jeans scale. Here we achieve the first cosmological simulations of this quantum state at unprecedentedly high resolution capable of resolving dwarf galaxies, with only one free parameter, $\bf{m_B}$, the boson mass. We demonstrate the large scale structure of this $ψ$DM simulation is indistinguishable from CDM, as desired, but differs radically inside galaxies. Connected filaments and collapsed haloes form a large interference network, with gravitationally self-bound solitonic cores inside every galaxy surrounded by extended haloes of fluctuating density granules. These results allow us to determine $\bf{m_B=(8.1^{+1.6}_{-1.7})\times 10^{-23}~eV}$ using stellar phase-space distributions in dwarf spheroidal galaxies. Denser, more massive solitons are predicted for Milky Way sized galaxies, providing a substantial seed to help explain early spheroid formation. Suppression of small structures means the onset of galaxy formation for $ψ$DM is substantially delayed relative to CDM, appearing at $\bf{z\lesssim 13}$ in our simulations.

preprint2014arXiv

Understanding the Core-Halo Relation of Quantum Wave Dark Matter, $ψ$DM, from 3D Simulations

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter ($ψ$DM), described by the Schrödinger-Poisson (SP) equation with a particle mass $\sim 10^{-22} {\rm eV}$. A distinct gravitationally self-bound solitonic core is found at the center of every halo, with a profile quite different from cores modeled in the warm or self-interacting dark matter scenarios. Furthermore, we show that each solitonic core is surrounded by an extended halo composed of large fluctuating dark matter granules which modulate the halo density on a scale comparable to the diameter of the solitonic core. The scaling symmetry of the SP equation and the uncertainty principle tightly relate the core mass to the halo specific energy, which, in the context of cosmological structure formation, leads to a simple scaling between core mass ($M_c$) and halo mass ($M_h$), $M_c \propto a^{-1/2} M_h^{1/3}$, where $a$ is the cosmic scale factor. We verify this scaling relation by (i) examining the internal structure of a statistical sample of virialized halos that form in our 3D cosmological simulations, and by (ii) merging multiple solitons to create individual virialized objects. Sufficient simulation resolution is achieved by adaptive mesh refinement and graphic processing units acceleration. From this scaling relation, present dwarf satellite galaxies are predicted to have kpc sized cores and a minimum mass of $\sim 10^8 {M_\odot}$, capable of solving the small-scale controversies in the cold dark matter model. Moreover, galaxies of $2\times10^{12} {M_\odot}$ at $z=8$ should have massive solitonic cores of $\sim 2\times10^9 {M_\odot}$ within $\sim 60 {\rm pc}$. Such cores can provide a favorable local environment for funneling the gas that leads to the prompt formation of early stellar spheroids and quasars.

preprint2011arXiv

Directionally Unsplit Hydrodynamic Schemes with Hybrid MPI/OpenMP/GPU Parallelization in AMR

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.

preprint2011arXiv

Vortex Turbulence in Linear Schroedinger Wave Mechanics

Quantum turbulence that exhibits vortex creation, annihilation and interactions is demonstrated as an exact solution of the time-dependent, free-particle Schrödinger 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.

preprint2009arXiv

GAMER: a GPU-Accelerated Adaptive Mesh Refinement Code for Astrophysics

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.