Search

This paper introduces a hybrid numerical scheme for the fuzzy dark matter model: It combines a wave-based approach to solve the Schr\"odinger equation using Fourier continuations with Gram polynomials and a fluid-based approach to solve the Hamilton-Jacobi-Madelung equations. This hybrid scheme facilitates zoom-in simulations for cosmological volumes beyond the capabilities of wave-based solvers alone and accurately simulates the full nonlinear dynamics of fuzzy dark matter. We detail the implementation of a Hamilton-Jacobi-Madelung solver, the methodology for phase matching at fluid-wave boundaries, the development of a local pseudospectral wave solver based on Fourier continuations, new grid refinement criteria for both fluid and wave solvers, an interpolation algorithm based on Fourier continuations, and the integration of these building blocks into the adaptive mesh refinement code GAMER. The superiority of the scheme is demonstrated through various performance and accuracy tests, tracking the linear power spectrum evolution in a 10 Mpc/h box, and a hybrid cosmological simulation in a 5.6 Mpc/h box. The corresponding code is published as part of the GAMER project on https://github.com/gamer-project/gamer., Comment: 25 pages, 18 figures, submitted to ApJS

It is common practice, both in dynamical modelling and in idealised numerical simulations, to assume that galaxies and/or dark matter haloes are spherical and have isotropic velocity distributions, such that their distribution functions are ergodic. However, there is no good reason to assume that this assumption is accurate. In this paper we use idealised $N$-body simulations to study the tidal evolution of subhaloes that are anisotropic at infall. We show that the detailed velocity anisotropy has a large impact on the subhalo's mass loss rate. In particular, subhaloes that are radially anisotropic experience much more mass loss than their tangentially anisotropic counterparts. In fact, in the former case, the stripping of highly radial orbits can cause a rapid cusp-to-core transformation, without having to resort to any baryonic feedback processes. Once the tidal radius becomes comparable to the radius of the core thus formed, the subhalo is tidally disrupted. Subhaloes that at infall are tangentially anisotropic are far more resilient to tidal stripping, and are never disrupted when simulated with sufficient resolution. We show that the preferential stripping of more radial orbits, combined with re-virialisation post stripping, causes an isotropisation of the subhalo's velocity distributions. This implies that subhaloes that have experienced significant mass loss are expected to be close to isotropic, which may alleviate the mass-anisotropy degeneracies that hamper the dynamical modelling of Milky Way satellites., Comment: 16 pages, 16 figures

The Fermi Gamma-Ray Space Telescope reveals two large bubbles in the Galaxy, extending nearly symmetrically $\sim50^{\circ}$ above and below the Galactic center (GC). Previous simulations of bubble formation invoking active galactic nucleus (AGN) jets have assumed that the jets are vertical to the Galactic disk; however, in general, the jet orientation does not necessarily correlate with the rotational axis of the Galactic disk. Using three-dimensional special relativistic hydrodynamic simulations including cosmic rays (CRs) and thermal gas, we show that the dense clumpy gas within the Galactic disk disrupts jet collimation ("failed jets" hereafter), which causes the failed jets to form hot bubbles. Subsequent buoyancy in the stratified atmosphere renders them vertical to form the symmetric Fermi and eROSITA bubbles (collectively, Galactic bubbles). We find that (1) despite the relativistic jets emanated from the GC are at various angles $\le45^{\circ}$ with respect to the rotational axis of the Galaxy, the Galactic bubbles nonetheless appear aligned with the axis; (2) the edge of the eROSITA bubbles corresponds to a forward shock driven by the hot bubbles; (3) followed by the forward shock is a tangling contact discontinuity corresponding to the edge of the Fermi bubbles; (4) assuming a leptonic model we find that the observed gamma-ray bubbles and microwave haze can be reproduced with a best-fit CR power-law spectral index of 2.4; The agreements between the simulated and the observed multi-wavelength features suggest that forming the Galactic bubbles by oblique AGN failed jets is a plausible scenario., Comment: 29 pages, 18 figures

Fuzzy dark matter (FDM), an attractive dark matter candidate comprising ultralight bosons (axions) with a particle mass $m_a\sim10^{-22}$ eV, is motivated by the small-scale challenges of cold dark matter and features a kpc-size de Broglie wavelength. Quantum wave interference inside an FDM halo gives rise to stochastically fluctuating density granulation; the resulting gravitational perturbations could drive significant disc thickening, providing a natural explanation for galactic thick discs. Here we present the first self-consistent simulations of FDM haloes and stellar discs, exploring $m_a=0.2-1.2\times10^{-22}$ eV and halo masses $M_\text{h} = 0.7-2.8\times10^{11}$ M$_\odot$. Disc thickening is observed in all simulated systems. The disc heating rates are approximately constant in time and increase substantially with decreasing $m_a$, reaching $dh/dt \simeq 0.04$ ($0.4$) kpc Gyr$^{-1}$ and $d\sigma_z^2/dt \simeq4$ ($150$) km$^2$s$^{-2}$Gyr$^{-1}$ for $m_a=1.2$ ($0.2$) $\times10^{-22}$ eV and $M_\text{h} =7\times10^{10} \text{M}_\odot$, where $h$ is the disc scale height and $\sigma_z$ is the vertical velocity dispersion. These simulated heating rates agree within a factor of two with the theoretical estimates of Chiang et al., confirming that the rough estimate of Church et al. overpredicts the granulation-driven disc heating rate by two orders of magnitude. However, the simulation-inferred heating rates scale less steeply than the theoretically predicted relation $d\sigma^2_z/dt \propto m_a^{-3}$. Finally, we examine the applicability of the Fokker-Planck approximation in FDM granulation modelling and the robustness of the $m_a$ exclusion bound derived from the Galactic disc kinematics., Comment: 20 pages, 26 figures; Accepted for publication in MNRAS

Superclusters, encompassing environments across a wide range of overdensities, can be regarded as unique laboratories for studying galaxy evolution. Although numerous supercluster catalogs have been published, none of them goes beyond redshift $z=0.7$. In this work, we adopt a physically motivated supercluster definition, requiring that superclusters should eventually collapse even in the presence of dark energy. Applying a friends-of-friends (FoF) algorithm to the CAMIRA cluster sample constructed using the Subaru Hyper Suprime-Cam survey data, we have conducted the first systematic search for superclusters at $z=0.5-1.0$ and identified 673 supercluster candidates over an area of 1027 deg$^2$. The FoF algorithm is calibrated by evolving $N$-body simulations to the far future to ensure high purity. We found that these high-$z$ superclusters are mainly composed of $2-4$ clusters, suggesting the limit of gravitationally bound structures in the younger Universe. In addition, we studied the properties of the clusters and brightest cluster galaxies (BCGs) residing in different large-scale environments. We found that clusters associated with superclusters are typically richer, but no apparent dependence of the BCG properties on large-scale structures is found. We also compared the abundance of observed superclusters with mock superclusters extracted from halo light cones, finding that photometric redshift uncertainty is a limiting factor in the performance of superclusters detection., Comment: Accepted by ApJ. 36 pages, 26 figures, 7 tables

Unveiling the true nature of Dark Matter (DM), which manifests itself only through gravity, is one of the principal quests in physics. Leading candidates for DM are weakly interacting massive particles (WIMPs) or ultralight bosons (axions), at opposite extremes in mass scales, that have been postulated by competing theories to solve deficiencies in the Standard Model of particle physics. Whereas DM WIMPs behave like discrete particles ($\varrho$DM), quantum interference between DM axions is manifested as waves ($\psi$DM). Here, we show that gravitational lensing leaves signatures in multiply-lensed images of background galaxies that reveal whether the foreground lensing galaxy inhabits a $\varrho$DM or $\psi$DM halo. Whereas $\varrho$DM lens models leave well documented anomalies between the predicted and observed brightnesses and positions of multiply-lensed images, $\psi$DM lens models correctly predict the level of anomalies left over by $\varrho$DM lens models. More challengingly, when subjected to a battery of tests for reproducing the quadruply-lensed triplet images in the system HS 0810+2554, $\psi$DM is able to reproduce all aspects of this system whereas $\varrho$DM often fails. The ability of $\psi$DM to resolve lensing anomalies even in demanding cases like HS 0810+2554, together with its success in reproducing other astrophysical observations, tilt the balance toward new physics invoking axions., Comment: Updated with accepted version. The Version of Record of this article is published in Nature Astronomy, and is available online at https://www.nature.com/articles/s41550-023-01943-9. Free online version : https://rdcu.be/damHj

Fuzzy Dark Matter (FDM), consisting of ultralight bosons, is an intriguing alternative to Cold Dark Matter. Numerical simulations solving the Schr\"odinger-Poisson (SP) equation, which governs FDM dynamics, show that FDM halos consist of a central solitonic core (representing the ground state of the SP equation), surrounded by a large envelope of excited states. Wave interference gives rise to order unity density fluctuations throughout the envelope and causes the soliton to undergo density oscillations and execute a confined random walk in the central region of the halo. The resulting gravitational potential perturbations are an efficient source of dynamical heating. Using high-resolution numerical simulations of a $6.6 \times 10^{9} \rm M_{\odot}$ FDM halo with boson mass, $m_{\rm b}=8 \times 10^{-23} \ \rm eV$, we investigate the impact of this dynamical heating on the structure and kinematics of spheroidal dwarf galaxies of a fixed mass but different initial sizes and ellipticities. The galaxies are set up in equilibrium in the time-and-azimuthally averaged halo potential and evolved for $10 \ \rm Gyr$ in the live FDM halo. We find that they continuously increase their sizes and central velocity dispersions. In addition, their kinematic structures become strongly radially anisotropic, especially in the outskirts. Dynamical heating also causes initially ellipsoidal galaxies to become more spherical over time from the inside out and gives rise to distorted, non-concentric isodensity contours. These tell-tale characteristics of dynamical heating of dwarf galaxies in FDM halos can potentially be used to constrain the boson mass., Comment: 18 pages, 11 figures, accepted for publication in ApJ

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

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

Fuzzy Dark Matter (FDM), consisting of ultralight bosons ($m_{\rm b} \sim 10^{-22}\ \rm eV$), is an intriguing alternative to Cold Dark Matter. Numerical simulations that solve the Schr\"odinger-Poisson (SP) equation show that FDM halos consist of a central solitonic core, which is the ground state of the SP equation, surrounded by an envelope of interfering excited states. These excited states also interfere with the soliton, causing it to oscillate and execute a confined random walk with respect to the halo center of mass. Using high-resolution numerical simulations of a $6.6 \times 10^9 M_{\odot}$ FDM halo with $m_{\rm b} = 8 \times 10^{-23}\ \rm eV$ in isolation, we demonstrate that the wobbling, oscillating soliton gravitationally perturbs nuclear objects, such as supermassive black holes or dense star clusters, causing them to diffuse outwards. In particular, we show that, on average, objects with mass $\lesssim 0.3 \%$ of the soliton mass ($M_{\rm sol}$) are expelled from the soliton in $\sim 3\ \rm Gyr$, after which they continue their outward diffusion due to gravitational interactions with the soliton and the halo granules. More massive objects ($\gtrsim 1 \% M_{\rm sol}$), while executing a random walk, remain largely confined to the soliton due to dynamical friction. We also present an effective treatment of the diffusion, based on kinetic theory, that accurately reproduces the outward motion of low mass objects and briefly discuss how the observed displacements of star clusters and active galactic nuclei from the centers of their host galaxies can be used to constrain FDM., Comment: 20 pages, 9 figures, minor revisions, accepted for publication in ApJ

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

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 $10^6$ 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\times 10^{7}$ 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.

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

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., Comment: Accepted for publication in ApJ. 14 pages, 7 figures

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

Without an existing large scale coherent magnetic field in the early Universe, Population III (PopIII) stars would likely rotate at or near break-up speed. In this work, focusing on the accretion phase of PopIII stars, we investigate the possibility of generating a coherent magnetic field through large scale dynamo processes, as well as the corresponding field saturation level. Using results from hydrodynamic simulations, we demonstrate that primordial accretion disks are turbulent with a Shakura-Sunyaev disk parameter $\alpha_{ss} \gtrsim 10^{-3}$, and evidence helical turbulence with a dynamo number $\vert D_{\alpha \Omega} \vert \gg 10$. The presence of helical turbulence at these levels allows large scale dynamo modes to grow, and the saturation level is determined by the amount of net helicity remaining in the dynamo-active regions (a.k.a. the quenching problem). %We demonstrate that the magnetic field can reach approximate equipartition, with $B/B_{\rm eq} \sim 3$, indicating that the dynamo quenching problem could be alleviated through an accretion flow. We demonstrate that, if the accretion could successfully alleviate the quenching problem, the magnetic field can reach approximate equipartition with $B/B_{\rm eq} \sim 3$.

Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general preserves both the energy and the symplectic structure of a Hamiltonian system. One is thus normally forced to make a choice. Nevertheless, a symmetric integration formula, obtained by Lanczos-Dyche via two-point Taylor expansion (or Hermite interpolation), is shown here to preserve both energy as well as symplectic structure for linear systems. This formula shares similarities with the Euler-Maclaurin formula, but is superconvergent rather than asymptotically convergent. For partial differential equations, the resulting evolution methods are unconditionally stable, i.e, not subject to a Courant-Friedrichs-Lewy limit. Although generally implicit, these methods become explicit for linear systems., Comment: Preprint submitted to J. Comp. Phys

Using cosmological particle hydrodynamical simulations and uniform ultraviolet backgrounds, we compare Lyman-$\alpha$ forest flux spectra predicted by the conventional cold dark matter (CDM) model, the free-particle wave dark matter (FP$\psi$DM) model and extreme-axion wave dark matter (EA$\psi$DM) models of different initial axion field angles against the BOSS Lyman-$\alpha$ forest absorption spectra with a fixed boson mass $m_b\sim 10^{-22}$eV. We recover results reported previously (Ir\v{s}i\v{c} et al. 2017b; Armengaud et al. 2017) that the CDM model agrees better with the BOSS data than the FP$\psi$DM model by a large margin, and we find the difference of total $\chi^2$'s is $120$ for $420$ data bins. These previous results demand a larger boson mass by a factor $>10$ to be consistent with the date and are in tension with the favoured value determined from local satellite galaxies. We however find that such tension is removed as some EA$\psi$DM models predict Lyman-$\alpha$ flux spectra agreeing better with the BOSS data than the CDM model, and the difference of total $\chi^2$'s can be as large as $24$ for the same bin number. This finding arises with no surprise since EA$\psi$DM models have unique spectral shapes with spectral bumps in excess of the CDM power near the small-scale cutoff typical of $\psi$DM linear matter power spectra as well as more extended cutoffs than FP$\psi$DM (Zhang & Chiueh 2017a,b).

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., Comment: 5 pages, 2 figures, minor typos fixed. Accepted for publication on Physics of Dark Universe

Acting as powerful gravitational lenses, the strong lensing galaxy clusters of the deep Hubble Frontier Fields (HFF) program permit access to lower-luminosity galaxies lying at higher redshifts than hitherto possible. We analyzed the HFF to measure the volume density of Lyman-break galaxies at $z > 4.75$ by identifying a complete and reliable sample up to $z \simeq 10$. A marked deficit of such galaxies was uncovered in the highly magnified regions of the clusters relative to their outskirts, implying that the magnification of the sky area dominates over additional faint galaxies magnified above the flux limit. This negative magnification bias is consistent with a slow rollover at the faint end of the UV luminosity function, and indicates a preference for Bose-Einstein condensate dark matter with a light boson mass of $m_\mathrm{B} \simeq 10^{-22} \, \mathrm{eV}$ over standard cold dark matter. We emphasize that measuring the magnification bias requires no correction for multiply lensed images (with typically three or more images per source), whereas directly reconstructing the luminosity function will lead to an overestimate unless such images can be exhaustively matched up, especially at the faint end that is accessible only in the strongly lensed regions. In addition, we detected a distinctive downward transition in galaxy number density at $z \gtrsim 8$, which may be linked to the relatively late reionization reported by Planck. Our results suggests that JWST will likely peer into an "abyss" with essentially no galaxies detected in deep NIR imaging at $z > 10$., Comment: 46 pages, 23 figures, accepted for publication in ApJ

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

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

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)

Light Axionic Dark Matter, motivated by string theory, is increasingly favored for the "no-WIMP era". Galaxy formation is suppressed below a Jeans scale, of $\simeq 10^8 M_\odot$ by setting the axion mass to, $m_B \sim 10^{-22}$eV, and the large dark cores of dwarf galaxies are explained as solitons on the de-Broglie scale. This is persuasive, but detection of the inherent scalar field oscillation at the Compton frequency, $\omega_B= (2.5{\rm \, months})^{-1}(m_B/10^{-22}eV)$, would be definitive. By evolving the coupled Schr\"odinger-Poisson equation for a Bose-Einstein condensate, we predict the dark matter is fully modulated by de-Broglie interference, with a dense soliton core of size $\simeq 150pc$, at the Galactic center. The oscillating field pressure induces General Relativistic time dilation in proportion to the local dark matter density and pulsars within this dense core have detectably large timing residuals, of $\simeq 400nsec/(m_B/10^{-22}eV)$. This is encouraging as many new pulsars should be discovered near the Galactic center with planned radio surveys. More generally, over the whole Galaxy, differences in dark matter density between pairs of pulsars imprints a pairwise Galactocentric signature that can be distinguished from an isotropic gravitational wave background., Comment: 6 pages, 3 figures. Accepted for publication in Physics Review Letter

Radio observations suggest that 3C 75, located in the dumbbell shaped galaxy NGC 1128 at the center of Abell 400, hosts two colliding jets. Motivated by this source, we perform three-dimensional hydrodynamical simulations using a modified version of the GPU-accelerated Adaptive-MEsh-Refinement hydrodynamical parallel code ($\mathit{GAMER}$) to study colliding extragalactic jets. We find that colliding jets can be cast into two categories: 1) bouncing jets, in which case the jets bounce off each other keeping their identities, and 2) merging jets, when only one jet emerges from the collision. Under some conditions the interaction causes the jets to break up into oscillating filaments of opposite helicity, with consequences for their downstream stability. When one jet is significantly faster than the other and the impact parameter is small, the jets merge; the faster jet takes over the slower one. In the case of merging jets, the oscillations of the filaments, in projection, may show a feature which resembles a double helix, similar to the radio image of 3C 75. Thus we interpret the morphology of 3C 75 as a consequence of the collision of two jets with distinctly different speeds at a small impact parameter, with the faster jet breaking up into two oscillating filaments., Comment: 13 pages, 9 figures, accepted for publication in ApJ

We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome., Comment: 20 pages, 8 figures, accepted for publication in MNRAS. For more info, visit grackle.readthedocs.io

We present a time domain waveform model that describes the inspiral-merger-ringdown (IMR) of compact binary systems whose components are non-spinning, and which evolve on orbits with low to moderate eccentricity. The inspiral evolution is described using third order post-Newtonian equations both for the equations of motion of the binary, and its far-zone radiation field. This latter component also includes instantaneous, tails and tails-of-tails contributions, and a contribution due to non-linear memory. This framework reduces to the post-Newtonian approximant TaylorT4 at third post-Newtonian order in the zero eccentricity limit. To improve phase accuracy, we incorporate higher-order post-Newtonian corrections for the energy flux of quasi-circular binaries and gravitational self-force corrections to the binding energy of compact binaries. This enhanced inspiral evolution prescription is combined with an analytical prescription for the merger-ringdown evolution using a catalog of numerical relativity simulations. This IMR waveform model reproduces effective-one-body waveforms for systems with mass-ratios between 1 to 15 in the zero eccentricity limit. Using a set of eccentric numerical relativity simulations, not used during calibration, we show that our eccentric model accurately reproduces the features of eccentric compact binary coalescence throughout the merger. Using this model we show that the gravitational wave transients GW150914 and GW151226 can be effectively recovered with template banks of quasi-circular, spin-aligned waveforms if the eccentricity $e_0$ of these systems when they enter the aLIGO band at a gravitational wave frequency of 14 Hz satisfies $e_0^{\rm GW150914}\leq0.15$ and $e_0^{\rm GW151226}\leq0.1$., Comment: 31 pages, 19 figures, 3 appendices. Submitted to Phys Rev D. v2: direct comparison to eccentric numerical relativity simulations, not used for calibration of the model, included. References added. Accepted to Phys Rev D

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

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, $\psi$DM, appropriate for light bosons such as axions, does not suffer this problem, given a boson mass of $m_{\psi} \ge 1.2 \times 10^{-22}{\,\rm eV}$ ($2\sigma$). 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 $\psi$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, $\psi$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., Comment: 14 pages, 12 figures, 2 tables. Correcting minor typos. Accepted for publication in ApJ

Inspired by Barnard 68, a Bok globule, that undergoes stable oscillations, we perform multi-phase hydrodynamic simulations to analyze the stability of Bok globules. We show that a high-density soft molecular core, with an adiabatic index $\gamma$ = 0.7 embedded in a warm isothermal diffuse gas, must have a small density gradient to retain the stability. Despite being stable, the molecular core can still collapse spontaneously as it will relax to develop a sufficiently large density gradient after tens of oscillations, or a few $10^7$ years. However, during its relaxation, the core may abruptly collapse triggered by the impingement of small-amplitude, long-wavelength ($\sim$ 6 $-$ 36 pc) sound waves in the warm gas. This triggered collapse mechanism is similar to a sonoluminescence phenomenon, where underwater ultrasounds can drive air bubble coalescence. The collapse configuration is found to be different from both inside-out and outside-in models of low-mass star formation; nonetheless the mass flux is close to the prediction of the inside-out model. The condition and the efficiency for this core collapse mechanism are identified. Generally speaking, a broad-band resonance condition must be met, where the core oscillation frequency and the wave frequency should match each other within a factor of several. A consequence of our findings predicts the possibility of propagating low-mass star formation, for which collapse of cores, within a mass range short of one order of magnitude, takes place sequentially tracing the wave front across a region of few tens of pc over $10^7$ years., Comment: 9 pages, 8 figures

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter ($\psi$DM), described by the Schr\"{o}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., Comment: 6 pages, 4 figures, accepted for publication in PRL

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 ($\psi$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 $\psi$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 $\psi$DM is substantially delayed relative to CDM, appearing at $\bf{z\lesssim 13}$ in our simulations., Comment: 14 pages, 8 figures. Pre-submission version. Final version has appeared in Nature Physics (http://dx.doi.org/10.1038/nphys2996)

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 is a GPU-accelerated Adaptive-MEsh-Refinement code for astrophysical simulations. In this work, two further extensions of the code are reported. First, we have implemented the MUSCL-Hancock method with the Roe's Riemann solver for the hydrodynamic evolution, by which the accuracy, overall performance and the GPU versus CPU speed-up factor are improved. Second, we have implemented the out-of-core computation, which utilizes the large storage space of multiple hard disks as the additional run-time virtual memory and permits an extremely large problem to be solved in a relatively small-size GPU cluster. The communication overhead associated with the data transfer between the parallel hard disks and the main memory is carefully reduced by overlapping it with the CPU/GPU computations., Comment: 4 pages, 4 figures, conference proceedings of IAU Symposium 270 (eds. Alves, Elmegreen, Girart, Trimble)

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

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., Comment: 60 pages, 22 figures, 3 tables. More accuracy tests are included. Accepted for publication in ApJS

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