Search

Elucidating the nature of Dark Matter (DM), which does not interact with light and which interacts with known matter primarily or 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 growing success of $\psi$DM in reproducing astrophysical observations tilt the balance toward new physics invoking axions., Comment: This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-023-01943-9 or https://www.nature.com/articles/s41550-023-01943-9. Please check out the final version on Nature Astronomy's website for the improved main text

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 axion is a hypothetical, well-motivated dark-matter particle whose existence would explain the lack of charge-parity violation in the strong interaction. In addition to this original motivation, an `axiverse' of ultra-light axions (ULAs) with masses $10^{-33}\,{\rm eV}\lesssim m_{\rm a}\lesssim 10^{-10}\,{\rm eV}$ also emerges from string theory. Depending on the mass, such a ULA contributes to the dark-matter density, or alternatively, behaves like dark energy. At these masses, ULAs' classical wave-like properties are astronomically manifested, potentially mitigating observational tensions within the $\Lambda$CDM paradigm on local-group scales. ULAs also provide signatures on small scales such as suppression of structure, interference patterns and solitons to distinguish them from heavier dark matter candidates. Through their gravitational imprint, ULAs in the presently allowed parameter space furnish a host of observational tests to target in the next decade, altering standard predictions for microwave background anisotropies, galaxy clustering, Lyman-$\alpha$ absorption by neutral hydrogen along quasar sightlines, pulsar timing, and the black-hole mass spectrum., Comment: 5 pages, 1 figure, Astro2020 Decadal Survey science white paper

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., Preprint submitted to J. Comp. Phys

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

Wave dark matter ($��$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 $��$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 $��$DM cosmological simulations. Moreover, the composite (star+$��$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., 15 pages, 11 figures

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

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., 14 pages, 12 figures, 2 tables. Correcting minor typos. Accepted for publication in ApJ

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., 6 pages, 4 figures, accepted for publication in PRL