Bullet Cluster Research Articles

A HYDRODYNAMICAL SOLUTION FOR THE “TWIN-TAILED” COLLIDING GALAXY CLUSTER “EL GORDO”

The distinctive cometary X-ray morphology of the recently discovered massive galaxy cluster "El Gordo" (ACT-CT J0102-4915; z=0.87) indicates that an unusually high-speed collision is ongoing between two massive galaxy clusters. A bright X-ray "bullet" leads a "twin-tailed" wake, with the SZ centroid at the end of the Northern tail. We show how the physical properties of this system can be determined using our FLASH-based, N-body/hydrodynamic model, constrained by detailed X-ray, Sunyaev-Zel'dovich (SZ), and Hubble lensing and dynamical data. The X-ray morphology and the location of the two Dark Matter components and the SZ peak are accurately described by a simple binary collision viewed about 480 million years after the first core passage. We derive an impact parameter of ~300 kpc, and a relative initial infall velocity of ~2250 km/sec when separated by the sum of the two virial radii assuming an initial total mass of 2.15x10^(15) Msun and a mass ratio of 1.9. Our model demonstrates that tidally stretched gas accounts for the Northern X-ray tail along the collision axis between the mass peaks, and that the Southern tail lies off axis, comprising compressed and shock heated gas generated as the massive component plunges through the main cluster. The challenge for LCDM will be to find out if this physically extreme event can be plausibly accommodated when combined with the similarly massive, high infall velocity case of the "Bullet cluster" and other such cases being uncovered in the new SZ based surveys.

The bullet cluster is not a cosmological anomaly

The Bullet Cluster (1E0657-56) merger is of exceptional interest for testing the standard cold-dark-matter plus cosmological constant cosmological model, and for investigating the possible existence of a long- or short-range ``fifth-force'' in the dark sector and possible need for modifications of general relativity or even of Newtonian gravity. The most recent previous simulations of the Bullet Cluster merger required an initial infall velocity far in excess of what would be expected within the standard cosmological model, at least in the absence of additional forces or modifications to gravity. We have recently carried out much more detailed simulations than previously had been done, making pixel-by-pixel fits to 2D data-maps of the mass distribution and X-ray emission, allowing for triaxial initial configurations and including MHD and cooling. Here, we compare the initial conditions of the Bullet Cluster merger as found in our new simulations to the initial conditions in similar-mass merging clusters in the Horizon cosmological simulation. We conclude that our initial infall velocity, 2900 km/s at a separation of 2.5 Mpc, is consistent with ΛCDM, given the inferred main cluster mass of 2 × 1015 M⊙. The initial concentration and shape found for the smaller (Bullet) cluster are typical for clusters of this mass range, but both quantities seem somewhat low for the larger (Main) cluster. Due to the paucity of examples of clusters with such a high mass in simulations, these features of the main cluster cannot presently be used to test ΛCDM.

Dark matter monopoles, vectors and photons

In a secluded dark sector which is coupled to the Standard Model via a Higgs portal interaction we arrange for the existence of 't Hooft-Polyakov magnetic monopoles and study their implications for cosmology. We point out that a dark sector which can accommodate stable monopoles will also contain massless dark photons gamma' as well as charged massive vector bosons W'. The dark matter in this scenario will be a combination of magnetically and electrically charged species under the unbroken U(1) subgroup of the dark sector. We estimate the cosmological production rate of monopoles and the rate of monopole-anti-monopole annihilation and conclude that monopoles with masses of few hundred TeV or greater, can produce sizeable contributions to the observed dark matter relic density. We scan over the parameter space and compute the relic density for monopoles and vector bosons. Turning to the dark photon radiation, we compute their contribution to the measured density of relativistic particles Neff and also apply observational constraints from the Bullet cluster and other large scale galaxies on long-range interactions for the self-interacting dark matter components made out of monopoles and out of dark vector bosons. At scales relevant for dwarf galaxies we identify regions on the parameter space where self-interacting monopole and vector dark mater components can aid solving the core-cusp and the too-big-to-fail problems.

NuSTAROBSERVATIONS OF THE BULLET CLUSTER: CONSTRAINTS ON INVERSE COMPTON EMISSION

The search for diffuse non-thermal inverse Compton (IC) emission from galaxy clusters at hard X-ray energies has been undertaken with many instruments, with most detections being either of low significance or controversial. Background and contamination uncertainties present in the data of non-focusing observatories result in lower sensitivity to IC emission and a greater chance of false detection. We present 266ks NuSTAR observations of the Bullet cluster, detected from 3-30 keV. NuSTAR's unprecedented hard X-ray focusing capability largely eliminates confusion between diffuse IC and point sources; however, at the highest energies the background still dominates and must be well understood. To this end, we have developed a complete background model constructed of physically inspired components constrained by extragalactic survey field observations, the specific parameters of which are derived locally from data in non-source regions of target observations. Applying the background model to the Bullet cluster data, we find that the spectrum is well - but not perfectly - described as an isothermal plasma with kT=14.2+/-0.2 keV. To slightly improve the fit, a second temperature component is added, which appears to account for lower temperature emission from the cool core, pushing the primary component to kT~15.3 keV. We see no convincing need to invoke an IC component to describe the spectrum of the Bullet cluster, and instead argue that it is dominated at all energies by emission from purely thermal gas. The conservatively derived 90% upper limit on the IC flux of 1.1e-12 erg/s/cm^2 (50-100 keV), implying a lower limit on B>0.2{\mu}G, is barely consistent with detected fluxes previously reported. In addition to discussing the possible origin of this discrepancy, we remark on the potential implications of this analysis for the prospects for detecting IC in galaxy clusters in the future.

Dark matter in cosmology

In this paper we review the main theoretical and experimental achievements in the field of dark matter from the cosmological and astrophysical point of view. We revisit it from the very first surveys of local astrophysical matter, up to the stringent constraints on matter properties, coming from the last release of data on cosmological scales. To bring closer and justify the idea of dark matter, we will go across methods and tools for measuring dark matter characteristics, and in some cases a combination of methods that provide one of the greatest direct proofs for dark matter, such as Bullet cluster.

On the interpretive role of theories of gravity and ‘ugly’ solutions to the total evidence for dark matter

Peter Kosso (2013) discusses the weak gravitational lensing observations of the Bullet Cluster and argues that dark matter can be detected in this system solely through the equivalence principle without the need to specify a full theory of gravity. This paper argues that Kosso gets some of the details wrong in his analysis of the implications of the Bullet Cluster observations for the Dark Matter Double Bind and the possibility of constructing robust tests of theories of gravity at galactic and greater scales. Even the Bullet Cluster evidence is not sufficiently detailed to allow precision tests of General Relativity that would distinguish it from its rivals at galactic and greater scales. Taking into account the total evidence available, we cannot rule out “ugly” solutions to the dynamical discrepancy in astrophysics that involve both a large quantity of dark matter and a theory of gravity whose predictions differ significantly from those of General Relativity for interactions taking place at galactic and greater scales.

Star formation in the massive cluster merger Abell 2744

We present a comprehensive study of star-forming (SF) galaxies in the HST Frontier Field recent cluster merger A2744 (z=0.308). Wide-field, ultraviolet-infrared (UV-IR) imaging enables a direct constraint of the total star formation rate (SFR) for 53 cluster galaxies, with SFR{UV+IR}=343+/-10 Msun/yr. Within the central 4 arcmin (1.1 Mpc) radius, the integrated SFR is complete, yielding a total SFR{UV+IR}=201+/-9 Msun/yr. Focussing on obscured star formation, this core region exhibits a total SFR{IR}=138+/-8 Msun/yr, a mass-normalised SFR{IR} of Sigma{SFR}=11.2+/-0.7 Msun/yr per 10^14 Msun and a fraction of IR-detected SF galaxies f{SF}=0.080(+0.010,-0.037). Overall, the cluster population at z~0.3 exhibits significant intrinsic scatter in IR properties (total SFR{IR}, Tdust distribution) apparently unrelated to the dynamical state: A2744 is noticeably different to the merging Bullet cluster, but similar to several relaxed clusters. However, in A2744 we identify a trail of SF sources including jellyfish galaxies with substantial unobscured SF due to extreme stripping (SFR{UV}/SFR{IR} up to 3.3). The orientation of the trail, and of material stripped from constituent galaxies, indicates that the passing shock front of the cluster merger was the trigger. Constraints on star formation from both IR and UV are crucial for understanding galaxy evolution within the densest environments.

CONSTRAINED SIMULATION OF THE BULLET CLUSTER

In this work, we report on a detailed simulation of the Bullet Cluster (1E0657-56) merger, including magnetohydrodynamics, plasma cooling, and adaptive mesh refinement. We constrain the simulation with data from gravitational lensing reconstructions and 0.5 - 2 keV Chandra X-ray flux map, then compare the resulting model to higher energy X-ray fluxes, the extracted plasma temperature map, Sunyaev-Zel'dovich effect measurements, and cluster halo radio emission. We constrain the initial conditions by minimizing the chi-squared figure of merit between the full 2D observational data sets and the simulation, rather than comparing only a few features such as the location of subcluster centroids, as in previous studies. A simple initial configuration of two triaxial clusters with NFW dark matter profiles and physically reasonable plasma profiles gives a good fit to the current observational morphology and X-ray emissions of the merging clusters. There is no need for unconventional physics or extreme infall velocities. The study gives insight into the astrophysical processes at play during a galaxy cluster merger, and constrains the strength and coherence length of the magnetic fields. The techniques developed here to create realistic, stable, triaxial clusters, and to utilize the totality of the 2D image data, will be applicable to future simulation studies of other merging clusters. This approach of constrained simulation, when applied to well-measured systems, should be a powerful complement to present tools for understanding X-ray clusters and their magnetic fields, and the processes governing their formation.

On the capture of dark matter by neutron stars

We calculate the number of dark matter particles that a neutron star accumulates over its lifetime as it rotates around the center of a galaxy, when the dark matter particle is a self-interacting boson but does not self-annihilate. We take into account dark matter interactions with baryonic matter and the time evolution of the dark matter sphere as it collapses within the neutron star. We show that dark matter self-interactions play an important role in the rapid accumulation of dark matter in the core of the neutron star. We consider the possibility of determining an exclusion region of the parameter space for dark matter mass and dark matter interaction cross section with the nucleons as well as dark matter self-interaction cross section, based on the observation of old neutron stars. We show that for a dark matter density of 103 GeV/cm3and dark matter mass mχ ≲ 10 GeV, there is a potential exclusion region for dark matter interactions with nucleons that is three orders of magnitude more stringent than without self-interactions. The potential exclusion region for dark matter self-interaction cross sections is many orders of magnitude stronger than the current Bullet Cluster limit. For example, for high dark matter density regions, we find that for mχ ∼ 10 GeV when the dark matter interaction cross section with the nucleons ranges from σχn ∼ 10−52 cm2 to σχn ∼ 10−57 cm2, the dark matter self-interaction cross section limit is σχχ ≲ 10−33 cm2, which is about ten orders of magnitude stronger than the Bullet Cluster limit.

Deep radio observations of the radio halo of the bullet cluster 1E 0657-55.8

We present deep 1.1-3.1 GHz Australia Telescope Compact Array observations of the radio halo of the bullet cluster, 1E 0657-55.8. In comparison to existing images of this radio halo the detection in our images is at higher significance. The radio halo is as extended as the X-ray emission in the direction of cluster merger but is significantly less extended than the X-ray emission in the perpendicular direction. At low significance we detect a faint second peak in the radio halo close to the X-ray centroid of the smaller sub-cluster (the bullet) suggesting that, similarly to the X-ray emission, the radio halo may consist of two components. Finally, we find that the distinctive shape of the western edge of the radio halo traces out the X-ray detected bow shock. The radio halo morphology and the lack of strong point-to-point correlations between radio, X-ray and weak-lensing properties suggests that the radio halo is still being formed. The colocation of the X-ray shock with a distinctive radio brightness edge illustrates that the shock is influencing the structure of the radio halo. These observations support the theory that shocks and turbulence influence the formation and evolution of radio halo synchrotron emission.

MEASURING THE STELLAR MASSES OF z ∼ 7 GALAXIES WITH THE SPITZER ULTRAFAINT SURVEY PROGRAM (SURFS UP)

We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

SPITZERULTRA FAINT SURVEY PROGRAM (SURFS UP). I. AN OVERVIEW

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey.

Direct detection portals for self-interacting dark matter

Dark matter self-interactions can affect the small scale structure of the Universe, reducing the central densities of dwarfs and low surface brightness galaxies in accord with observations. From a particle physics point of view, this points toward the existence of a 1--100 MeV particle in the dark sector that mediates self-interactions. Since mediator particles will generically couple to the Standard Model, direct detection experiments provide sensitive probes of self-interacting dark matter. We consider three minimal mechanisms for coupling the dark and visible sectors: photon kinetic mixing, $Z$ boson mass mixing, and the Higgs portal. Self-interacting dark matter motivates a new benchmark paradigm for direct detection via momentum-dependent interactions, and ton-scale experiments will cover astrophysically motivated parameter regimes that are unconstrained by current limits. Direct detection is a complementary avenue to constrain velocity-dependent self-interactions that evade astrophysical bounds from larger scales, such as those from the bullet cluster.

Dark matter, the Equivalence Principle and modified gravity

Dark matter (DM) is an essential ingredient of the present Standard Cosmological Model, according to which only 5% of the mass/energy content of our universe is made of ordinary matter. In recent times, it has been argued that certain cases of gravitational lensing represent a new type of evidence for the existence of DM. In a recent paper, Peter Kosso attempts to substantiate that claim. His argument is that, although in such cases DM is only detected by its gravitational effects, gravitational lensing is a direct consequence of Einstein's Equivalence Principle (EEP) and therefore the complete gravitational theory is not needed in order to derive such lensing effects. In this paper I critically examine Kosso's argument: I confront the notion of empirical evidence involved in the discussion and argue that EEP does not have enough power by itself to sustain the claim that gravitational lensing in the Bullet Cluster constitutes evidence for the DM Hypothesis. As a consequence of this, it is necessary to examine the details of alternative theories of gravity to decide whether certain empirical situations are indeed evidence for the existence of DM. It may well be correct that gravitational lensing does constitute evidence for the DM Hypothesis—at present it is controversial whether the proposed modifications of gravitation all need DM to account for the phenomenon of gravitational lensing and if so, of which kind—but this will not be a direct consequence of EEP.

THE IMPACT OF SUZAKU MEASUREMENTS ON ASTROPARTICLE PHYSICS

Results from the &lt;em&gt;Suzaku&lt;/em&gt; X-ray broad-band observations of clusters of galaxies are summarized. Aiming at understanding the physics of gas heating/particle acceleration and the cluster dynamical evolution, we search for non-thermal hard X-ray emission from merging clusters, particularly A2163 and the Bullet Cluster, based on the &lt;em&gt;Suzaku&lt;/em&gt; and &lt;em&gt;XMM-Newton/Chandra&lt;/em&gt; joint analyses. The observed hard X-ray emission is well represented by single- or multi-temperature thermal models, including super-hot (kT ~ 20 keV) gas. However, no significant non-thermal hard X-ray emission has been detected. Together with the presently available literature, the hard X-ray properties have been studied for about 10 clusters with &lt;em&gt;Suzaku&lt;/em&gt;. The present status on &lt;em&gt;Suzaku&lt;/em&gt; measurements of non-thermal X-ray emission and the cluster magnetic field are summarized and compared with those from the &lt;em&gt;RXTE, BeppoSAX&lt;/em&gt;, and &lt;em&gt;Swift&lt;/em&gt; satellites. The future prospects are briefly mentioned.

Statistics of extreme objects in the Juropa Hubble Volume simulation★

We present the first results from the JUropa huBbLE volumE (Jubilee) project, based a large N-body, dark matter-only cosmological simulation with a volume of $V=(6 h^{-1}\mathrm{Gpc})^3$, containing 6000$^3$ particles, performed within the concordance $\Lambda$CDM cosmological model. The simulation volume is sufficient to probe extremely large length scales in the universe, whilst at the same time the particle count is high enough so that dark matter haloes down to $1.5\times10^{12} h^{-1}\mathrm{M}_\odot$ can be resolved. At $z = 0$ we identify over 400 million haloes. The cluster mass function is derived using three different halofinders and compared to fitting functions in the literature. The distribution of clusters of maximal mass across redshifts agrees well with predicted masses of extreme objects, and we explicitly confirm that the Poisson distribution is very good at describing the distribution of rare clusters. The Poisson distribution also matches well the level to which cosmic variance can be expected to affect number counts of high mass clusters. We find that objects like the Bullet cluster exist in the far-tail of the distribution of mergers in terms of relative collisional speed. We also derive the number counts of voids in the simulation box for $z = 0$, $0.5$ and $1$.

THE PRE-MERGER IMPACT VELOCITY OF THE BINARY CLUSTER A1750 FROM X-RAY, LENSING, AND HYDRODYNAMICAL SIMULATIONS

Since the discovery of the "bullet cluster" several similar cases have been uncovered suggesting relative velocities well beyond the tail of high speed collisions predicted by the concordance LCDM model. However, quantifying such post-merger events with hydrodynamical models requires a wide coverage of possible initial conditions. Here we show that it is simpler to interpret pre-merger cases, such as A1750, where the gas between the colliding clusters is modestly affected, so that the initial conditions are clear. We analyze publicly available Chandra data confirming a significant increase in the projected X-ray temperature between the two cluster centers in A1750 consistent with our expectations for a merging cluster. We model this system with a self-consistent hydrodynamical simulation of dark matter and gas using the FLASH code. Our simulations reproduce well the X-ray data, and the measured redshift difference between the two clusters in the phase before the first core passage viewed at an intermediate projection angle. The deprojected initial relative velocity derived using our model is 1460 km/sec which is considerably higher than the predicted mean impact velocity for simulated massive haloes derived by recent LCDM cosmological simulations, but it is within the allowed range. Our simulations demonstrate that such systems can be identified using a multi-wavelength approach and numerical simulations, for which the statistical distribution of relative impact velocities may provide a definitive examination of a broad range of dark matter scenarios.

TANGENTIAL VELOCITY OF THE DARK MATTER IN THE BULLET CLUSTER FROM PRECISE LENSED IMAGE REDSHIFTS

We show that the fast moving component of the "bullet cluster" (1E0657-56) can induce potentially resolvable redshift differences between multiply-lensed images of background galaxies. The moving cluster effect can be expressed as the scalar product of the lensing deflection angle with the tangential velocity of the mass components, and it is maximal for clusters colliding in the plane of the sky with velocities boosted by their mutual gravity. The bullet cluster is likely to be the best candidate for the first measurement of this effect due to the large collision velocity and because the lensing deflection and the cluster fields can be calculated in advance. We derive the deflection field using multiply-lensed background galaxies detected with the Hubble Space Telescope. The velocity field is modeled using self-consistent N-body/hydrodynamical simulations constrained by the observed X-ray and gravitational lensing features of this system. We predict that the triply-lensed images of systems "G" and "H" straddling the critical curve of the bullet component will show the largest frequency shifts up to ~0.5 km/sec. This is within the range of the Atacama Large Millimeter/sub-millimeter Array (ALMA) for molecular emission, and is near the resolution limit of the new generation high-throughput optical-IR spectrographs. A detection of this effect measures the tangential motion of the subclusters directly, thereby clarifying the tension with LCDM, which is inferred from gas motion less directly. This method may be extended to smaller redshift differences using the Ly-alpha forest towards QSOs lensed by more typical clusters of galaxies. More generally, the tangential component of the peculiar velocities of clusters derived by our method complements the radial component determined by the kinematic SZ effect, providing a full 3-dimensional description of velocities.

Finslerian MOND versus the strong gravitational lensing of the early-type galaxies

The gravitational lensing of Bullet Clusters and early-type galaxies pose serious challenges on the validity of MOND. Recently, Finslerian MOND, a generalization of MOND in the framework of Finsler gravity, has been proposed to explain the mass discrepancy problem of Bullet Cluster 1E 0657\ 558. In this paper, we check the validity of the Finslerian MOND in describing the strong gravitational lensing of early-type galaxies. The investigation on ten strong lenses of the CASTLES samples shows that there is no strong evidence for the existence of dark matter.

THE DYNAMICS OF MERGING CLUSTERS: A MONTE CARLO SOLUTION APPLIED TO THE BULLET AND MUSKET BALL CLUSTERS

Merging galaxy clusters have become one of the most important probes of dark matter, providing evidence for dark matter over modified gravity and even constraints on the dark matter self-interaction cross-section. To properly constrain the dark matter cross-section it is necessary to understand the dynamics of the merger, as the inferred cross-section is a function of both the velocity of the collision and the observed time since collision. While the best understanding of merging system dynamics comes from N-body simulations, these are computationally intensive and often explore only a limited volume of the merger phase space allowed by observed parameter uncertainty. Simple analytic models exist but the assumptions of these methods invalidate their results near the collision time, plus error propagation of the highly correlated merger parameters is unfeasible. To address these weaknesses I develop a Monte Carlo method to discern the properties of dissociative mergers and propagate the uncertainty of the measured cluster parameters in an accurate and Bayesian manner. I introduce this method, verify it against an existing hydrodynamic N-body simulation, and apply it to two known dissociative mergers: 1ES 0657-558 (Bullet Cluster) and DLSCL J0916.2+2951 (Musket Ball Cluster). I find that this method surpasses existing analytic models - providing accurate (10% level) dynamic parameter and uncertainty estimates throughout the merger history. This coupled with minimal required a priori information (subcluster mass, redshift, and projected separation) and relatively fast computation (~6 CPU hours) makes this method ideal for large samples of dissociative merging clusters.