Explore the vast range of titles available.

A cosmological model treating dark matter as a coherent quantum wave agrees well with conventional dark-matter theory on an astronomical scale. But on smaller scales, the quantum nature of wave-like dark matter can explain dark-matter cores that are observed in dwarf galaxies, which standard theory cannot. The conventional cold-particle interpretation of dark matter (known as ‘cold dark matter’, or 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 1 , 2 , 3 , 4 , 5 . In contrast, galaxies predicted by CDM extend to much lower masses, with steeper, singular profiles 6 , 7 , 8 , 9 . 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 10 , 11 , 12 . Here we achieve cosmological simulations of this quantum state at unprecedentedly high resolution capable of resolving dwarf galaxies, with only one free parameter, m B , the boson mass. We demonstrate the large-scale structure is indistinguishable from CDM, as desired, but differs radically inside galaxies where quantum interference forms solitonic cores surrounded by extended haloes of fluctuating density granules. These results allow us to determine 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. The onset of galaxy formation is substantially delayed relative to CDM, appearing at redshift z ≲ 13 in our simulations.

In the exascale computing era, handling and analyzing massive data sets have become extremely challenging. In situ analysis, which processes data during simulation runtime and bypasses costly intermediate disk input and output steps, offers a promising solution. We present libyt (https://github.com/yt-project/libyt), an open-source C library that enables astrophysical simulations to analyze and visualize data in parallel computation with yt or other Python packages. libyt can invoke Python routines automatically or provide interactive entry points via a Python prompt or a Jupyter Notebook. It requires minimal intervention in researchers’ workflows, allowing users to reuse job submission scripts and Python routines. We describe libyt’s architecture for parallel computing in high-performance computing environments, including its bidirectional connection between simulation codes and Python, and its integration into the Jupyter ecosystem. We detail its methods for reading patch-based adaptive mesh refinement simulations and handling in-memory data with minimal overhead, and procedures for yielding data when requested by Python. We describe how libyt maps simulation data to yt front ends, allowing postprocessing scripts to be converted into in situ analysis with just two lines of change. We document libyt’s application programming interface (API) and demonstrate its integration into two astrophysical simulation codes, GAMER and Enzo, using examples including core-collapse supernovae, isolated dwarf galaxies, fuzzy dark matter, the Sod shock tube test, Kelvin–Helmholtz instability, and the AGORA galaxy simulation. Finally, we discuss libyt’s performance, limitations related to data redistribution, extensibility, architecture, and comparisons with traditional postprocessing approaches.

Fuzzy dark matter (FDM), consisting of ultralight bosons, is an intriguing alternative to cold dark matter. Numerical simulations solving the Schrödinger–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 density fluctuations of order unity 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 × 109 M ⊙ FDM halo with boson mass m b = 8 × 10−23 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 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, nonconcentric isodensity contours. These telltale characteristics of dynamical heating of dwarf galaxies in FDM halos can potentially be used to constrain the boson mass.

This paper introduces a hybrid numerical scheme for the fuzzy dark matter (FDM) model: it combines a wave-based approach to solve the Schrödinger 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 FDM. 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–1 box, and a hybrid cosmological simulation in a 5.6 Mpc h–1 box. The corresponding code is published as part of the GAMER project on https://github.com/gamer-project/gamer.

The Fermi Gamma-Ray Space Telescope reveals two large bubbles in the Galaxy, extending nearly symmetrically ∼50° 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 emanating from the GC at various angles ≤45° 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 multiwavelength features suggest that forming the Galactic bubbles by oblique AGN failed jets is a plausible scenario.

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 deg2. 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 two to four 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.

The cold dark matter (CDM) paradigm successfully explains the cosmic structure over an enormous span of redshifts. However, it fails when probing the innermost regions of dark matter halos and the properties of the Milky Way’s dwarf galaxy satellites. Moreover, the lack of experimental detection of Weakly Interacting Massive Particle (WIMP) favors alternative candidates such as light axionic dark matter that naturally arise in string theory. Cosmological N-body simulations have shown that axionic dark matter forms a solitonic core of size of ≃150 pc in the innermost region of the galactic halos. The oscillating scalar field associated to the axionic dark matter halo produces an oscillating gravitational potential that induces a time dilation of the pulse arrival time of ≃400 ns/(m B /10 − 22 eV) for pulsar within such a solitonic core. Over the whole galaxy, the averaged predicted signal may be detectable with current and forthcoming pulsar timing array telescopes.

Unveiling the true nature of dark matter, which manifests itself only through gravity, is one of the principal quests in physics. Leading candidates for dark matter are weakly interacting massive particles 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 dark matter weakly interacting massive particles behave like discrete particles (ϱDM), quantum interference between dark matter axions is manifested as waves (ψ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 ϱDM or ψDM halo. Whereas ϱDM lens models leave well documented anomalies between the predicted and observed brightnesses and positions of multiply lensed images, ψDM lens models correctly predict the level of anomalies remaining with ϱ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, ψDM is able to reproduce all aspects of this system whereas ϱDM often fails. The ability of ψDM to resolve lensing anomalies even in demanding cases such as HS 0810+2554, together with its success in reproducing other astrophysical observations, tilt the balance toward new physics invoking axions.Modelling of the gravitationally lensed system HS 0810+2554 with wavelike dark matter resolves brightness and position anomalies remaining after the standard massive-particle dark matter treatment.

The Fanaroff–Riley Class II radio galaxy Cygnus A hosts jets that 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 used relativistic hydrodynamic simulations of a collimated jet interacting with an inclined interface between lobe and cluster plasma to model the basic processes that may lead to such a feature. We found 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 produced 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.

The origin of the cool-core (CC) and non-cool-core (NCC) dichotomy of galaxy clusters remains uncertain. Previous simulations have found that cluster mergers are effective in destroying CCs but fail to prevent overcooling in cluster cores when radiative cooling is included. Feedback from active galactic nuclei (AGN) is a promising mechanism for balancing cooling in CCs; however, the role of AGN feedback in CC/NCC transitions remains elusive. In this work, we perform 3D binary cluster merger simulations incorporating AGN feedback and radiative cooling, aiming to investigate the heating effects from mergers and AGN feedback on CC destruction. We vary the mass ratio and impact parameter to examine the entropy evolution of different merger scenarios. We find that AGN feedback is essential in regulating the merging clusters, and that CC destruction depends on the merger parameters. Our results suggest three scenarios regarding CC/NCC transitions: (1) CCs are preserved in minor mergers or mergers that do not trigger sufficient heating, in which cases AGN feedback is crucial for preventing the cooling catastrophe; (2) CCs are transformed into NCCs by major mergers during the first core passage, and AGN feedback is subdominant; (3) in major mergers with a large impact parameter, mergers and AGN feedback operate in concert to destroy the CCs.