Issue 
A&A
Volume 628, August 2019



Article Number  A100  
Number of page(s)  9  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201936184  
Published online  13 August 2019 
An ALMA+ACA measurement of the shock in the Bullet Cluster
^{1}
MaxPlanckInstitut für Astrophysik (MPA), KarlSchwarzschildStrasse 1, Garching 85741, Germany
email: lucadim@mpagarching.mpg.de
^{2}
European Southern Observatory (ESO), KarlSchwarzschildStrasse 2, Garching 85748, Germany
^{3}
Space Research Institute, Profsoyuznaya 84/32, Moscow 117997, Russia
^{4}
Astrophysics Science Division, NASA Goddard Space Flight Center, Greenbelt, MD, 20771, USA
^{5}
Argelander Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
^{6}
Naval Research Laboratory, 4555 Overlook Avenue SW, Code 7213, Washington, DC 20375, USA
^{7}
Department of Physics and Astronomy, University of Pennsylvania, 209 South 33rd Street, Philadelphia, PA 19104, USA
^{8}
National Radio Astronomy Observatory (NRAO), 520 Edgemont Road, Charlottesville, VA 22903, USA
^{9}
HarvardSmithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
^{10}
Moorpark College, 7075 Campus Rd., Moorpark, CA 93021, USA
^{11}
Department of Physics and Astronomy, The University of Utah, Salt Lake City, Utah, USA
Received:
26
June
2019
Accepted:
19
July
2019
Context. The thermal SunyaevZeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through Xray observations. Widely regarded as the “textbook example” of a cluster merger bow shock, the western, mostprominent shock front in the Bullet Cluster (1E065756) represents the ideal test case for such an SZ study.
Aims. We aim to characterize the shock properties using deep, highresolution interferometric SZ effect observations in combination with priors from an independent Xray analysis.
Methods. Our analysis technique relies on the reconstruction of a parametric model for the SZ signal by directly and jointly fitting data from the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA) in Fourier space.
Results. The ALMA+ACA data are primarily sensitive to the electron pressure difference across the shock front. To estimate the shock Mach number ℳ, this difference can be combined with the value for the upstream electron pressure derived from an independent Chandra Xray analysis. In the case of instantaneous electronion temperature equilibration, we find ℳ = 2.08_{−0.12}^{+0.12}, in ≈ 2.4σ tension with the independent constraint from Chandra, M_{X} = 2.74 ± 0.25. The assumption of purely adiabatic electron temperature change across the shock leads to ℳ = 2.53_{−0.25}^{+0.33}, in better agreement with the Xray estimate ℳ_{X} = 2.57 ± 0.23 derived for the same heating scenario.
Conclusion. We have demonstrated that interferometric observations of the thermal SZ effect provide constraints on the properties of the shock in the Bullet Cluster that are highly complementary to Xray observations. The combination of Xray and SZ data yields a powerful probe of the shock properties, capable of measuring ℳ and addressing the question of electronion equilibration in cluster shocks. Our analysis is however limited by systematics related to the overall cluster geometry and the complexity of the postshock gas distribution. To overcome these limitations, a simultaneous, jointlikelihood analysis of SZ and Xray data is needed.
Key words: galaxies: clusters: individual: 1E 065756 / galaxies: clusters: intracluster medium / cosmic background radiation
© L. Di Mascolo et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Open Access funding provided by Max Planck Society.
1. Introduction
Mergers play a crucial role in the formation of galaxy clusters, which are situated at intersections of the Cosmic Web. These spectacular events can have a profound impact on the intracluster medium (ICM) and the galaxies within these environments (see, e.g., Kravtsov & Borgani 2012). Mergers provide largescale astrophysical laboratories for plasmas where the mean free path can be substantial (see, e.g., Markevitch & Vikhlinin 2007) and for measuring the selfinteraction crosssection of dark matter (Markevitch et al. 2004; Randall et al. 2008; Wittman et al. 2018; Tulin & Yu 2018). Furthermore, the very existence of dark matter was conclusively demonstrated through the merging cluster 1E0657–56, or “Bullet Cluster”, which exhibits spatial offsets between its baryonic and total mass peaks in the Xray and gravitational lensing maps (Clowe et al. 2006; Bradač et al. 2006; Paraficz et al. 2016).
Key to identifying merging clusters is the detection of shocks in the ICM. A “textbook example of a bow shock” is observed in the Xray image of the Bullet Cluster (Markevitch et al. 2002). Using 500 ks of Chandra Xray data, Markevitch (2006) reported a Mach number for the western, mostprominent shock in the Bullet Cluster, an estimate largely determined by the density jump conditions. We also note that Shimwell et al. (2015) revealed a second shock, on the eastern (opposite) side of the cluster, which we do not consider here.
Here we present deep, continuum Atacama Large Millimeter/submillimeter Array (ALMA) observations, sensitive to the thermal SunyaevZeldovich (tSZ; Sunyaev & Zeldovich 1972) effect, of the main shock in the Bullet Cluster. These observations include data from both the 12 m array (hereafter “ALMA”) and 7 m Atacama Compact (Morita) Array (ACA). As the tSZ effect is linearly sensitive to the lineofsight integral of the electron thermal pressure (see, e.g., Carlstrom et al. 2002; Mroczkowski et al. 2019a, for reviews), these observations complement the Xray constraints on plasma density and, less accurately, electron temperature, yielding a groundbased, mmwave view of the shock properties.
All the results presented in this work have been derived assuming a flat ΛCDM cosmology with Ω_{m} = 0.30, Ω_{Λ} = 0.70, and H_{0} = 70.0 km s^{−1} Mpc^{−1}. At the redshift of the Bullet Cluster (z = 0.296), 1″ corresponds to a physical scale of 4.41 kpc. Unless stated differently, all reported bestfitting parameters and their respective uncertainties are obtained from the 50th, 16th, and 84th percentiles of the marginalized posterior distributions, corresponding to the 68% credibility interval of the distribution^{1}.
2. Data and analysis overview
2.1. ALMA/ACA observations
As part of ALMA Cycle 2 operations ALMA and the ACA observed the Bullet Cluster for a total of 3.1 and 5.9 h integration time respectively in Band 3 (project code 2013.1.00760.S). These wideband observations span the frequency range 84−100 GHz in four 2 GHzwide spectral windows, centered at approximately 85, 87, 97, and 99 GHz. Our strategy employed a single, deep observation centered approximately on the nose of the shock front, as inferred from the Xray observations (Clowe et al. 2006; Markevitch 2006). The ALMA and ACA observations were respectively performed in 4 and 11 separate executions spanning 2014, obtaining rootmeansquare noise levels of approximately 5 μJy and 45 μJy respectively (as measured in naturallyweighted imaging), and a synthesized beam with a main lobe of 4.01″ × 3.07″ FWHM (PA 81°).
We rereduced the data using the ALMA pipeline (Shinnaga et al. 2015; Humphreys et al. 2016) in CASA 4.7 (McMullin et al. 2007), producing results consistent with the previous calibration using the script provided on data delivery. Our rereduction provided a crosscheck of the earlier reduction, and was necessary due to backwardscompatibility issues and bug fixes in subsequent CASA releases. The data were calibrated using the default calibration strategy of the ALMA observatory, which has nominal uncertainty ≤5%. However, since the flux calibrators, which included quasars, differed for each execution, we performed a manual crosscheck of the values for the designated flux calibrators as well as the phase calibrators, finding they were consistent for the dates spanned by the observations.
An interferometer behaves as a spatial filter, sampling the sky only in Fourier modes corresponding to projected baselines in the array (see, e.g., Condon & Ransom 2016; Di Mascolo et al. 2019). This provides clean imaging free from atmospheric structure, but also leads to two major complexities: incomplete sampling in Fourier space even for the modes accessible to the array, and the lack of recovery of angular scales larger than those corresponding to the shortest projected distances between array elements (i.e., the “missing flux issue”). Based on the uvspace coverage of the ALMA and ACA data presented here, the largest recoverable scales are respectively ∼40″ and ∼55″. As detailed in the next section, we choose to forwardmodel the observed SZ signal using Xraymotivated priors to address such issues. Additionally, to avoid known deconvolution biases intrinsic to the CLEAN algorithm (Högbom 1974; Thompson et al. 2015), we perform our analysis directly in visibility (uv) space. We extend the interferometric SZ analysis techniques presented in the appendix of Di Mascolo et al. (2019) to allow fitting pressure discontinuities due to shocks. Our approach builds upon the work of Basu et al. (2016), but incorporates several advances in the parameterspace sampling technique as well as more sophisticated and flexible models allowed by the deeper Xray and SZ observations. In brief, we build an imagespace model of the SZ signal by integrating numerically the threedimensional pressure distribution model, and applying the proper SZ frequency scaling. The dependence of the SZ signal on the electron temperature is taken into account when modelling the SZ spectrum (Itoh & Nozawa 2004; Chluba et al. 2012, and end of Sect. 2.2.2). The pixel scale is chosen to fulfill the Nyquist sampling criterion for the smallest scales probed. The SZ model image is then Fourier transformed and sampled to the position of the sparse interferometric data. The resulting synthetic visibilities are then employed in combination with the observed ones to evaluate the likelihood at each step of the Bayesian inference procedure.
However, we choose not to model the raw postcalibration data, instead binning the data in each spectral window following the optimal averaging scheme described in Hobson et al. (1995). This is crucial for gaining a significant reduction in data volume and hence computational time.
2.2. SunyaevZeldovich model
A summary of the model priors introduced in this section can be found in Table 1. We test for biases in the parameter reconstruction arising from our specific choice for the distribution of priors by performing a prioronly run, which is done by setting the likelihood to a constant value regardless of model fit (see, e.g., Di Mascolo et al. 2019). As expected, the result of this test simply returns the input distribution of priors.
Priors on the model parameters employed in our analysis.
2.2.1. Shock front
The common approach employed in the study of Xray observations of shock fronts consists in describing them as spherical sectors within a specific region of the cluster image. This takes advantage of the imagespace nature of the Xray data to select a spatial region narrow enough to allow one to locally approximate the shock front as spherical. However, among the complexities of studying interferometric data is the difficulty of applying any spatial masking. This would entail convolution of the visibilities, inducing a nontrivial correlation between them. To avoid this, a complete twodimensional model of the observed field is then required.
In order to allow more freedom in the description of the shock front than in the case of a spherical model, we describe the shock front as an axiallysymmetric hyperbolic surface (see the dashed line in Fig. 1), with central axis coincident with the direction of the merger and lying in the plane of the sky. Since the interferometric data alone cannot constrain the lineofsight distribution of pressure, we consider the curvature of the front to be symmetric with respect to the lineofsight and the planeofsky direction. Although this is likely a reasonable assumption, any deviations from cylindrical symmetry may introduce nonnegligible systematic errors into our results. In particular, the derived downstream pressure P_{e, ds} will be related to the true pressure as , where ℓ_{LoS} and are respectively the assumed and true lineofsight extents of the shock front. We present specific estimates for the ratio when discussing our results in the next sections. It is worth noting that, as already discussed by Wang et al. (2018), this is expected to scale approximately with the square root of the curvature radius of the shock surface. For lineofsight and planeofsky curvatures radii r_{LoS} and r_{PoS}, we would then get . We refer to Appendix A for a discussion.
A preliminary attempt to determine the shape of the shock front using SZ data alone shows that the parameters defining its geometry are heavily degenerate, and the small extent of the ALMA+ACA field of view does not allow for meaningful constraints. We would like to note this is a consequence of the sole parametrization of the shock geometry, since ALMA+ACA has proven to be able to identify edge positions with a beamscale precision (Basu et al. 2016). We therefore derive a description of the morphology of the bow shock by finding the hyperbola that best describes the discontinuity observed in the Chandra surface brightness map. Analogous to Ueda et al. (2017), we find the bestmatching shock geometry by minimizing the variance of the Xray image within a defined region. The model is assumed to be simply given by a step function in which the discontinuity has a hyperbolic shape. The values of the function inside and outside the front itself are set equal to the mean photon counts in the respective regions of the Xray image. In order to gain better leverage on the azimuthal geometry of the shock front, we further split the region in several angular sectors (see Fig. 1). The resulting maximumaposterior model for the hyperbolic surface is then employed for describing the profile of the shock front in all the following analyses. To account for a possible mismatch in the shock coordinates (e.g. due to astrometry errors) from the Chandra modelling with respect to ALMA+ACA, we allow for some additional freedom in the nose coordinates (RA, Dec)_{shock} and axis orientation θ. Specifically, we assign each one priors based on the respective marginalized posteriors derived in the Xraymatching step described above. Any asymmetry in the recovered parameter uncertainties is modelled by means of splitnormal distributions (Wallis 2014).
Fig. 1.
Cutout of the 0.5−2.0 keV Chandra Xray surface brightness map of the Bullet Cluster. The solid wedges represent the sectors employed to derive the hyperbolic shape bestmatching the shock front geometry. The maximumaposteriori model is shown as a dashed line, while the dotted contours indicate the corresponding 95% credible interval. Further, the dashed and solid yellow circles denote respectively the ACA and ALMA fullwidthathalfmaximum fields of view. For reference, the upstream and downstream gas lie respectively west (right) and east (left) of the shock front. 
Apart from its morphology, the main parameter defining the shock is the ratio of downstream to upstream pressure at the jump itself. In practice, for a fixed lineofsight geometry, the available ALMA+ACA data are mainly sensitive to the absolute difference of the downstream and upstream electron pressures near the nose of the shock, i.e., ΔP_{e} = P_{e, ds} − P_{e, us}; the SZ signal associated with the largescale distribution of the gas is effectively filtered out (see, e.g., Basu et al. 2016, and discussion in Sect. 3.2). Thus, the modelling of the ALMA signal remains only weakly sensitive to the assumed largescale model. The immediate downside is the relative pressure jump at the shock, x_{P} = P_{e,ds}/P_{e,us} = 1 + ΔP_{e}/P_{e,us}, which serves as a proxy for ℳ, is poorly constrained by the interferometric SZ data alone. In fact, due to the lack of information on the pressure normalization, the marginalized posterior distribution of the Mach number inferred when performing an SZonly analysis are found to entirely span the corresponding prior interval. To get a meaningful measure of the pressure jump from the ALMA+ACA data, we therefore employ an Xrayinformed analysis of the ALMA+ACA SZ observations as in Sect. 3.3 of Basu et al. (2016), and set the upstream electron pressure P_{e, us} to the value derived by modelling the Chandra data in a narrow sector centered on the shock nose (Markevitch 2006).
We do not include any model components describing the “bullet” itself (i.e., the contact discontinuity or “cold front”) or the subtle additional cold front between the bullet and the main shock reported by Markevitch (2006) and Markevitch & Vikhlinin (2007). The former lies outside the ALMA field of view, precluding any interesting constraints on the pressure difference (or lack thereof) across the cold front, while the latter is intrinsically faint, and is expected to be in thermal pressure equilibrium. For simplicity, we thus assume these features have a negligible effect on the measurements of the shock itself, and therefore ignore them. Further, we assume a single powerlaw profile for the downstream electron pressure (see below). A future analysis, joint with Xrays, will allow more model freedom for trying to build description of such features.
2.2.2. Bulk pressure distribution
To model the pressure distribution in the downstream region, we employ a power law radial profile with slope α, centered along the merger axis at a distance from the shock nose equal to the front curvature radius r_{c},
On the other hand, we consider the preshock pressure distribution to be relaxed, thus to be described by a spherical generalized NavarroFrenkWhite (gNFW) profile (Nagai et al. 2007):
where r_{500} and P_{500} are functions of M_{500}, the total mass contained within an average overdensity 500× the critical density of the Universe at that redshift (Arnaud et al. 2010). We constrain M_{500} so that the gNFW model always reproduces the Xray value for the upstream pressure, P_{e, us}. For the main results reported here, the pressure normalization P_{0}, concentration parameter c_{500}, massdependence index a_{P}, and slopes (a, b, c) are fixed to the values reported in Arnaud et al. (2010) for the universal pressure profile. However, we show below that our results for ℳ are insensitive to the choice of gNFW parametrization and position of the gNFW model centroid. We therefore simply fix the gNFW centroid’s coordinates (RA, Dec)_{gNFW} to the position of main lensing κmap peak inferred by Clowe et al. (2006), which we note does not coincide with the center of the postshock profile.
Though pressure perturbations driven by the primary merger are confined to be within the shock front, it is possible for the passage of its associated dark matter component to affect, through infall, the ICM ahead of the shock (Springel & Farrar 2007). This may undermine our choice of the universal gNFW profile, reliable in the case of relaxed clusters, when describing the bulk pressure distribution. We tested against possible systematics introduced by this assumption. We found no significant deviations in the reconstructed parameters after changing either the slopes of the profiles or the position of the assumed centroids. The same applies to the structure of the gas on the downstream side as long as it is smooth, even though it may differ from expectations for a solid body moving through homogeneous fluid (see, e.g., Zhang et al. 2019). The results of the above tests are summarized in Fig. 2. The net result is that, as a consequence of the interferometric filtering, our shock model is largely sensitive to the pressure conditions right at the front, and not to the properties of the bulk pressure distribution (see Sect. 3.2). Thus, we consider wide uninformative priors on both the Mach number and postshock slope, and marginalize over the latter.
Fig. 2.
Bivariate posterior density function (left panel) for the inferred Mach number ℳ and slope of the downstream pressure distribution α for a set of upstream pressure profiles and respective centroids. To facilitate comparison, the bestfitting ℳ estimates are also plotted (right panel). We consider three different cases for the underlying gNFW profile by setting the slopes equal to the values reported in Arnaud et al. (2010) for the coolcore (cc; red), morphologicallydisturbed (md; peach), or ensemble (up; blue) cluster samples. For each of them, the distribution centroid is then fixed to a number of different positions: far downstream and far upstream (right panel, left and midleft points), respectively to arbitrary distances of 10′ east and 3′ west of the Xrayderived shock nose coordinates; APEXSZ centroid (midright point; Halverson et al. 2009); peak of the κmap (right point; Clowe et al. 2006). In all the above cases, we assumed instantaneous shock heating of the electrons. The gray line in the right panel denotes the corresponding bestfitting ℳ reported in Sect. 3.1, while the darker and lighter bands the respective 68% and 95% credibility intervals. 
In order to account for the high temperatures measured in the system, relativistic corrections to the SZ spectrum (Itoh & Nozawa 2004) are included in our modelling. In fact, variations in the measured SZ signal of the order of 5% up to 15% are expected for an electron gas with temperature ranging from 9 keV to 30 keV as measured from the Xray data. As for P_{e, us}, we employ an Xraymotivated prior on the upstream temperature T_{e, us}. In addition, we incorporate the 5% uncertainties on the ACA and ALMA flux calibration by introducing normalization hyperparameters κ_{ACA} and κ_{ALMA} (Di Mascolo et al. 2019). For all the modelling runs presented in the following sections, κ_{ACA} and κ_{ALMA} have been found to not deviate significantly from unity.
Given the planeofsky geometry of the merger involving the Bullet Cluster, any contribution from the kinetic SZ effect (Sunyaev & Zeldovich 1980) to the observed signal due to the motion of the single subclusters should be subdominant with respect to the thermal SZ effect. In fact, if we assume the velocity v ≈ 3000 inferred from the shock Mach number (Springel & Farrar 2007) to be measured with respect to the CMB rest frame and the merger direction to be oriented by around 8° with respect to the plane of sky (Markevitch et al. 2004), we find that the contribution of the kinetic SZ effect to the total SZ signal from the postshock region would be of the order of 3% of the corresponding thermal component. This would induce a systematic error on the estimate of the Mach number ℳ lower than 2%. Given the small effect as well as the lack of robust constraints on the merger proper velocities and orientation, we then decide to not include the kinetic SZ effect in our model, keeping the merger axis aligned with plane of sky (see Sect. 2.2.1).
2.3. Implementation details
As noted earlier, the modelling algorithm and its specific implementation are detailed in Di Mascolo et al. (2019, see also references therein). However, we improved the posterior sampling algorithm by adopting the dynamic nested sampling by Higson et al. (2018). In particular, we employ the purePython implementation provided by dynesty (Speagle 2019).
3. Results
3.1. Instantaneous electronion temperature equilibration
We first consider a model for the shock front under the standard assumption of an instantaneous electronion temperature equilibration, i.e., T_{e} = T_{i}, at the shock front. This is consistent with the Chandra analysis by Markevitch (2006), who derived ℳ_{X} = 3.0 ± 0.4 from the density jump, and measured an electron temperature jump as expected if electrons reach the average postshock temperature nearinstantly. A more recent analysis of the same Chandra data by Markevitch (in prep.) provides an electron density jump x_{n} = 2.86 ± 0.16, assuming instantaneous electron heating (T_{e, ds} ≈ 30 keV) for converting the Xray surface brightness to the density jump by accounting for the different emissivity in the post and preshock regions (see, e.g., Ettori 2000 for a discussion about Xray brightness modelling in the presence of temperature gradients). This corresponds to ℳ_{X} = 2.74 ± 0.25. The slight difference from the older estimate is due to a bettercentered shock model and the inclusion of the azimuthal decline in x_{n} for angles away from the shock nose when projecting along the l.o.s. (as in Wang et al. 2018), based on the amplitude of the density jumps measured from the Xray data in different sectors of the shock (Markevitch, in prep.).
The assumption of instantaneous heating implies that the nominal RankineHugoniot condition can be used to relate ℳ to the measured amplitude of the electron pressure jump x_{P} relative to the upstream value P_{e, us} as
Here, γ is the polytropic exponent, which we assume to be γ = 5/3, appropriate for nonrelativistic fullyionized gas. We further allow for the azimuthal variation of ℳ along the shock front, whose scaling with the azimuthal angle is derived using the same density jump decline discussed above. The omission of such azimuthal dependence would cause the Mach number to be averaged down with respect to its maximum value due to the effect of the wings with lower ℳ. For the results provided in this and the following sections, we estimate that the inclusion of the Xraybased model for the azimuthal variation of the shock pressure jump increases the value of the inferred Mach number by only 5–7%. However, more severe effects should be expected for observations with larger field of views, which would include values from farther in the wings.
We obtain (Fig. 3). While the model relies on the Xray priors on the preshock pressure and temperature, the derived Mach number is inconsistent at a 2.4σ level with the Xray estimate ℳ_{X} = 2.74 ± 0.25. Projection effects may play a nonnegligible role in biasing the SZbased measurement of the Mach number. However, a strong ellipticity of the shock front shape (see Sect. 2.2.1) would be required to bridge the gap between SZ and Xray estimates. In reality, an even larger ellipticity would be necessary, given that the Xray estimates would also be affected by geometry, albeit with a different dependence. Another potential source of bias is the Xraymotivated prior on the upstream pressure, which comes from deprojected density and temperature estimates, used to compute the relative pressure jump. While the definition of a centroid for Xray deprojection remains ambiguous, we found only extreme choices would alter our results significantly. A jointlikelihood Xray+SZ analysis may be required to find a consistent geometry that fully reconciles such discrepancies.
3.2. Collisional electronion temperature equilibration
Here we consider the possibility that the electron and ion temperatures do not equilibrate instantaneously in cluster shocks (i.e., T_{e} ≠ T_{i} immediately inside the shock front; see, e.g., Fox & Loeb 1997; Markevitch 2006; Russell et al. 2012; Wang et al. 2018). Ions carry the majority of the gas bulk kinetic energy in collisionless shocks, and are heated dissipatively on scales comparable to their gyroradii, while electrons might remain much colder (Vink et al. 2015), unless there is some process that equilibrates the ion and electron temperatures. The upper limit on the equilibration time scales is set by Coulomb collisions (Zeldovich & Raizer 1966), which for the downstream density and temperature in the Bullet Cluster is long (∼few 10^{8} yr), occurring over a distance comparable to the offset between the shock and the cold front.
Under the assumption of conservation of the enthalpy flux, electrons equilibrate with ions to the RankineHugoniot downstream temperature at a rate driven by Coulomb collisions (Fox & Loeb 1997)
where t_{eq} is the Coulomb collisional timescale (Spitzer 1962), n_{e} and n_{i} are respectively the electron and ion densities, and x_{T} is the temperature ratio across the shock front. To build our SZ model, we convert the above equation in terms of the distance from the shock front by means of the downstream gas velocity u_{ds} = (ℳ/x_{n}) c_{us}, with c_{us} given by the upstream sound speed. Furthermore, we assume that electrons are first heated adiabatically (Vink et al. 2015), so that the electron temperature immediately inside the shock front equals . As required by the conservation of charge neutrality across the shock front, the density jump x_{n} is also set to follow the RankineHugoniot condition. In this case, the pressure jump x_{P} cannot be directly related to ℳ as in Eq. (3), and instead must be derived as the product of the density and temperature ratios at each threedimensional model coordinate.
In the case of collisional equilibration, we find (Fig. 3). This is consistent with ℳ_{X} = 2.57 ± 0.23, coming from the Chandra Xray brightness fit if one uses the adiabaticcompression postshock temperature (T_{e, ds} ≈ 20 keV) to convert to the density jump (Markevitch, in prep.). Unfortunately, due to the severe filtering of large spatial scales as well as the limited field of view, we are not able to put any significant constraint on the specific equipartition timescale when treating t_{eq} as a free parameter. Instead, we find that assuming the electronion equilibration to be driven by Coulomb collisions is practically equivalent to setting t_{eq} = ∞.
Fig. 3.
Marginalized posterior distributions for the shock Mach number ℳ derived under the assumptions of instantaneous (red) and collisional (blue) electronion equilibration. The dashed and dotted lines indicates the median of the posterior distributions and the 68% credible intervals, respectively. 
For illustrative purposes, we present in Fig. 4 the dirty^{2} images of the raw ALMA+ACA data employed in our analysis, the interferometric model corresponding to the collisional electronion equilibration scenario, and the respective modelsubtracted data. As shown in the right panel, it is not possible to identify residuals that differ at a significant level from noiselike features.
Fig. 4.
Dirty images of the raw (left), model (middle), and residual (right) ALMA+ACA interferometric data. They are generated by jointly gridding the ACA and ALMA data using a multifrequency natural weighting scheme. For reference, this provides a synthesized beam of 4.07″ × 3.01″ FWHM (PA 81°; bottomleft corner of right panel). We cut off the fields at the 0.2 gain level of the joint ALMA+ACA antenna pattern. To better highlight the largescale shock features, we avoid correcting for the primary beam attenuation and apply an additional 30 kλ taper. We note that the model subtraction is performed directly in Fourier space. The dashed line in the center and right panels indicate the inferred position of the shock front. Xray contours are overlaid on the left panel from Fig. 1. We also note the positive signal at the shock front is not due to an increment of the SZ signal; rather it is an artifact of the highpass filtering effects of ALMA+ACA (see also Fig. 5). 
The analogous image for the instantaneous case is visually identical to Fig. 4, and therefore is not shown. This is evident in Fig. 5, where it is not possible to identify any significant difference between the filtered SZ models for the instantaneous and collisionalequilibration scenarios. This confirms that the ALMA+ACA data are only sensitive to the properties of the pressure distributions near the shock edge, thus providing a direct estimate of the pressure difference across the shock front, i.e., ΔP_{e} rather than x_{P}. Moreover, our result is found to be practically independent of the specific assumption about the underlying gNFW profile, which is entirely filtered by the interferometric response (Fig. 5).
Fig. 5.
Comparison of SZ signal profiles (upper panel) and the corresponding electron pressure profiles (lower panel) across the shock nose. The vertical gray line denotes the position of the shock front. We note that the upper panel contains both unfiltered input model fits for the SZ signal and the corresponding filtered (observed) profiles (see legend for details). As discussed in Sect. 3.2, the instantaneous (dashdotted yellow) and collisional (solid red) shock models are indistinguishable after the spatial filtering of the interferometric ALMA+ACA observation, reflecting the fundamental limitation of ALMA+ACA to constrain any largescale (≳1′) component of the SZ signal. Again, we note that while we measure a decrement due to the SZ effect, the filtered (observed) profiles can exhibit both positive and negative excursions (analogous to the Gibbs phenomenon; see, e.g., Bracewell 1978). For comparison, we also report the input unfiltered model (solid yellow) for the instantaneous equilibration, and both the raw and filtered underlying gNFW profile (dashed and dotted blue lines, as noted in the legend). We also present the Xray expectation for both the filtered and raw SZ signal profiles (dashed and solid green lines) and corresponding pressure profiles that we would expect for the case of instantaneous equilibration, using the value ℳ_{X} = 2.74 ± 0.25 derived from fits to the Xray data. 
The fact ALMA+ACA is only sensitive in practice to the electron pressure difference ΔP_{e} across the shock is also reflected in the lack of any significant difference between the Bayesian evidences of the instant equilibration and adiabatic heating models (). In fact, given that we cannot observe any largescale feature in the SZ signal induced by the slow increase of the postshock electron temperature in the case of collisional equilibration, the different heating scenarios practically differ only in the way we convert the pressure difference to an estimate of the shock ℳ.
We recall here a subtle cold front is observed in the Chandra Xray image between the shock and the bullet boundary, ∼15″ east from the shock. Since the total pressure across a cold front is expected to be approximately continuous, this was not included in our modeling of the downstream pressure profile. In fact, no apparent signature of such feature can be distinguished in the ALMA+ACA SZ observation. However, because the flow of the postshock gas would not cross the cold front, it is unlikely that any electronion temperature nonequilibrium would extend past it. If indeed there is an electronproton temperature difference in the postshock region, we should expect the electron temperature (and, hence, pressure) to reach its equilibrium value. This would in turn result in a discontinuity in the SZ signal in the direction of the front itself. Future, more sensitive ALMA observations may search for such a feature.
4. Conclusions
We further demonstrate the ability of using deep, highresolution ALMA+ACA observations of the SZ effect to characterize shocks in merging clusters (see, for comparison, Basu et al. 2016). For this purpose, we studied the SZ effect across the shock in the Bullet Cluster, chosen as it is widely regarded as the “textbook example” of a cluster merger bow shock.
The application of our interferometric modelling technique – using Xraymotivated priors – has allowed us to place constraints on the electron pressure discontinuity across the shock. Assuming a RankineHugoniot shock adiabat, our pressure jump implies a Mach number , which is significantly lower than the one derived from Chandra data using the same geometric assumptions (ℳ = 2.74 ± 0.25). An interesting physical possibility to reconcile the two measurements is to allow that the electron and ion temperatures do not equilibrate instantly after the shock passage has heated the electrons adiabatically. For a given Mach number, this would lower the postshock electron temperature and thus the observed electron pressure jump. Our Mach number would then become , in agreement with the Xray estimate that assumes the adiabatic temperature jump for conversion between the Xray brightness and density. We note that Chandra Xray data constrain the gas density (from which the Mach number is derived) and electron temperature across the shock separately, and its postshock temperature prefers instant equilibration over adiabatic heating of the electrons (at ∼2σ confidence; Markevitch 2006). However, while Chandra is free from ALMA+ACA interferometric limitations and can probe the upstream and downstream gas directly, the Bullet postshock temperature is above the range where Chandra can measure electron temperatures reliably, and hence suffer significant systematic uncertainties.
To summarize, ALMA+ACA has proven to provide a clean measurement of the differential jump in pressure due to the shock, and, in combination with data that can access larger scales, can provide compelling constraints on shock properties such as the Mach number. In particular:

Interferometric observations cleanly measure the projected pressure jump due to the shock. However, due to the inherent spatial filtering of ALMA and the ACA, which recover scales ∼0.5−1.1′ in Band 3 (compared to θ_{500}∼10′, corresponding to r_{500} for the Bullet Cluster), Xray priors on both the model geometry and upstream pressure are necessary in order to infer ℳ. On the other hand, the combination of SZ observations covering a broader range of spatial scales (i.e., from 0.1″ − 30′) could provide an SZonly view of the shock properties.

Once the geometry is fixed, the key quantities which drive the analysis and the interpretations of the results are the pressure difference across the shock front, the normalization of the preshock pressure, and the independent Xray estimates of ℳ. We show their statistical uncertainties are small enough to allow us to differentiate between the instantaneous and adiabatic heating scenarios. Nevertheless, neither model is unambiguously preferred. Although the two scenarios result in SZbased estimates for ℳ that deviate one from the other by ∼2σ, we find that the difference of the respective Bayesian logevidence is not significant enough to completely rule out one versus another.

We extensively tested our modelling choices – varying the geometry, pre and postshock pressure slopes, and underlying pressure distribution – and find our results to be robust for a broad range of possible assumptions motivated by the Xray analyses. However, our model does not fully describe the complex morphology observed in the Xray surface brightness. Together with the uncertainties on the threedimensional morphology of the cluster, this may limit our ability to elucidate the nature of electron heating across the shock front.
Together, these illustrate the fundamental complementarity of Xray and SZ effect observations in the study of the physics of galaxy clusters. It is then clear that a simultaneous, jointlikelihood analysis of the SZ and Xray data on the Bullet Cluster, extending the approach of incorporating Xray information in the form of priors (see our discussion in Sect. 2.2.1, and Sect. 3.3 of Basu et al. 2016), would benefit our understanding of the morphology of the galaxy cluster, as well as provide further insights into the physical mechanisms for shock heating of the intracluster medium. A forthcoming paper will present the results of a full jointlikelihood analysis of interferometric SZ and Xray observations, as well as singledish SZ measurements, building on the methodology discussed in Di Mascolo et al. (2019). Meanwhile, upcoming results from NuSTAR (Wik et al. in prep.) will better access the high photon energies corresponding to the high temperatures inferred from Chandra. Further, the number of observations of shocks with unambiguous geometry and sufficiently high Mach number that allow the detection of deviations from instantaneous electron heating is limited. Thus, along with improved modelling, observations of a larger sample of cluster shocks will be needed to improve our understanding. And finally, both improved spatial and spectral resolution, larger instantaneous field of view, and the ability to recover zerospacing information will vastly improve future SZonly studies. However, in order to provide sufficient overlap with the interferometric data in Fourier space, while also probing higher frequencies and spatial scales > 10′, a new widefield (> 1°) singledish facility, such as the Atacama Large Aperture Submillimeter/millimeter Telescope (AtLAST; see, e.g., Klaassen et al. 2019; Mroczkowski et al. 2019b) is required.
Acknowledgments
We thank the anonymous referee for the insightful comments. EC and RS acknowledges partial support from the Russian Science Foundation grant 191200369. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00760.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. Basic research in radio astronomy at the US Naval Research Laboratory is funded by 6.1 Base funding, supporting T. Clarke. The scientific results reported in this article are obtained using the following softwares: AstroPy, a communitydeveloped core Python package for Astronomy (Astropy Collaboration et al. 2018); dynesty (Speagle 2019); galario (Tazzari et al. 2018).
References
 Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Astropy Collaboration (PriceWhelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
 Basu, K., Sommer, M., Erler, J., et al. 2016, ApJ, 829, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Bracewell, R. 1978, The Fourier Transform and its Applications, 2nd edn. (Tokyo: McGrawHill Kogakusha, Ltd.) [Google Scholar]
 Bradač, M., Clowe, D., Gonzalez, A. H., et al. 2006, ApJ, 652, 937 [NASA ADS] [CrossRef] [Google Scholar]
 Carlstrom, J. E., Holder, G. P., & Reese, E. D. 2002, ARA&A, 40, 643 [NASA ADS] [CrossRef] [Google Scholar]
 Chluba, J., Nagai, D., Sazonov, S., & Nelson, K. 2012, MNRAS, 426, 510 [NASA ADS] [CrossRef] [Google Scholar]
 Clowe, D., Bradač, M., Gonzalez, A. H., et al. 2006, ApJ, 648, L109 [Google Scholar]
 Condon, J. J., & Ransom, S. M. 2016, Essential Radio Astronomy (Princeton University Press) [Google Scholar]
 Di Mascolo, L., Churazov, E., & Mroczkowski, T. 2019, MNRAS, 487, 4037 [CrossRef] [Google Scholar]
 Ettori, S. 2000, MNRAS, 311, 313 [NASA ADS] [CrossRef] [Google Scholar]
 Fox, D. C., & Loeb, A. 1997, ApJ, 491, 459 [NASA ADS] [CrossRef] [Google Scholar]
 Halverson, N. W., Lanting, T., Ade, P. A. R., et al. 2009, ApJ, 701, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Higson, E., Handley, W., Hobson, M., & Lasenby, A. 2018, Stat. Comput., https://doi.org/10.1007/s1122201898440 [Google Scholar]
 Hobson, M. P., Lasenby, A. N., & Jones, M. 1995, MNRAS, 275, 863 [NASA ADS] [CrossRef] [Google Scholar]
 Högbom, J. A. 1974, A&AS, 15, 417 [Google Scholar]
 Humphreys, E., Miura, R., Brogan, C. L., et al. 2016, Proc. 2016 ALMA Conf., 1 [Google Scholar]
 Itoh, N., & Nozawa, S. 2004, A&A, 417, 827 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Klaassen, P., Mroczkowski, T., Bryan, S., et al. 2019, ArXiv eprints [arXiv:1907.04756] [Google Scholar]
 Kravtsov, A. V., & Borgani, S. 2012, ARA&A, 50, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Markevitch, M. 2006, in The Xray Universe 2005, ed. A. Wilson, ESA Spec. Publ., 604, 723 [NASA ADS] [Google Scholar]
 Markevitch, M., & Vikhlinin, A. 2007, Phys. Rep., 443, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Markevitch, M., Gonzalez, A. H., David, L., et al. 2002, ApJ, 567, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Markevitch, M., Gonzalez, A. H., Clowe, D., et al. 2004, ApJ, 606, 819 [NASA ADS] [CrossRef] [Google Scholar]
 McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in Astronomical Data Analysis Software and Systems XVI, eds. R. A. Shaw, F. Hill, & D. J. Bell, ASP Conf. Ser., 376, 127 [Google Scholar]
 Mroczkowski, T., Nagai, D., Basu, K., et al. 2019a, Space Sci. Rev., 215, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Mroczkowski, T., Nagai, D., Andreani, P., et al. 2019b, BAAS, 51, 124 [Google Scholar]
 Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007, ApJ, 668, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Paraficz, D., Kneib, J.P., Richard, J., et al. 2016, A&A, 594, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Randall, S. W., Markevitch, M., Clowe, D., Gonzalez, A. H., & Bradač, M. 2008, ApJ, 679, 1173 [NASA ADS] [CrossRef] [Google Scholar]
 Russell, H. R., McNamara, B. R., Sanders, J. S., et al. 2012, MNRAS, 423, 236 [NASA ADS] [CrossRef] [Google Scholar]
 Shimwell, T. W., Markevitch, M., Brown, S., et al. 2015, MNRAS, 449, 1486 [NASA ADS] [CrossRef] [Google Scholar]
 Shinnaga, H., Humphreys, E., Indebetouw, R., et al. 2015, in Revolution in Astronomy with ALMA: The Third Year, eds. D. Iono, K. Tatematsu, A. Wootten, L. Testi, et al., ASP Conf. Ser., 499, 355 [NASA ADS] [Google Scholar]
 Speagle, J. S. 2019, MNRAS, submitted [arXiv:1904.02180] [Google Scholar]
 Spitzer, L. 1962, Physics of Fully Ionized Gases (New York: Interscience ) [Google Scholar]
 Springel, V., & Farrar, G. R. 2007, MNRAS, 380, 911 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, I. B. 1980, MNRAS, 190, 413 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comm. Astrophys. Space Phys., 4, 173 [Google Scholar]
 Tazzari, M., Beaujean, F., & Testi, L. 2018, MNRAS, 476, 4527 [Google Scholar]
 Thompson, R., Davé, R., & Nagamine, K. 2015, MNRAS, 452, 3030 [NASA ADS] [CrossRef] [Google Scholar]
 Tulin, S., & Yu, H.B. 2018, Phys. Rep., 730, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Ueda, S., Kitayama, T., & Dotani, T. 2017, ApJ, 837, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Vink, J., Broersen, S., Bykov, A., & Gabici, S. 2015, A&A, 579, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wallis, K. F. 2014, Statist. Sci., 29, 106 [Google Scholar]
 Wang, Q. H. S., Giacintucci, S., & Markevitch, M. 2018, ApJ, 856, 162 [NASA ADS] [CrossRef] [Google Scholar]
 Wittman, D., Golovich, N., & Dawson, W. A. 2018, ApJ, 869, 104 [NASA ADS] [CrossRef] [Google Scholar]
 Zeldovich, Y. B., & Raizer, Y. P. 1966, Elements of Gasdynamics and the Classical Theory of Shock Waves. [Google Scholar]
 Zhang, C., Churazov, E., Forman, W. R., & Jones, C. 2019, MNRAS, 482, 20 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Interferometric view of the SZ effect
The amplitude of the thermal SZ effect in the direction x of a galaxy cluster characterized by an electron pressure distribution P_{e} is proportional to Compton y,
where σ_{T}, m_{e}, and c are respectively the Thomson crosssection, the electron mass, and the speed of light, while ℓ is the coordinate along the line of sight.
As shown in Sect. 3.2 and Fig. 5, the available ALMA+ACA data (i) only probe a small region near the tip of the shock and (ii) any extended structures in this region are effectively filtered out. Therefore, a model that includes a jump of pressure across the shock front and is smooth otherwise should capture the most of the information contained in the ALMA+ACA data. To this end, we represent the electron pressure distribution P_{e}(x, ℓ) as a combination of two spatially smooth components, P_{e, 0}(x, ℓ) and P_{e, 1}(x, ℓ), where their sum is
Here f_{ds} ≡ f_{ds}(x, ℓ) is equal to 1 in the downstream region and 0 in the upstream region (i.e., is a Heaviside step function). We note both that here and below, for simplicity, we omit the explicit coordinate dependence of P_{e} or y. Rearranging the terms and integrating along the line of sight we obtain
Since the first term in the above expression corresponds to a smooth, largescale pressure distribution its contribution to y is filtered out from the ALMA+ACA data. Moreover, the function (P_{e, 1} − P_{e, 0}) in the second term is also spatially smooth and would be filtered too without the step function f_{ds}. Therefore, the signal measured by ALMA+ACA in the vicinity of the shock tip is effectively defined by the second term in Eq. (A.3), which is set by the pressure jump at the shock front and the lengthscale of the downstream region. Thus,
where ℓ_{LoS} is the lineofsight extent of the probed postshock region, and P_{e, us} and P_{e, ds} are the electron pressures measured just outside and inside the shock front, respectively (Sect. 2.2). Therefore, ALMA+ACA data effectively constrain a product of the electron pressure difference at the shock ΔP_{e} and the physical size of the region ℓ_{LoS}. The latter quantity can be easily determined if the merger is in the plane of the sky and the shock front possesses rotational symmetry. If the shape of the front can be approximated by a sphere with a radius R, then along the symmetry axis , where h ≪ r is the distance from the tip of the shock. While the calculations in the paper were done without these simplifying assumptions, the Eq. (A.4) is useful to estimate the uncertainty introduced by the (unknown) geometry of the shock along the line of sight. In particular, if the curvatures in the sky plane r_{PoS} and along the line of sight r_{LoS} differ, the estimate of the pressure difference ΔP_{e}, which assumes r_{LoS} = r_{Pos}, will be biased by a factor (the same argument is discussed in Wang et al. 2018).
All Tables
All Figures
Fig. 1.
Cutout of the 0.5−2.0 keV Chandra Xray surface brightness map of the Bullet Cluster. The solid wedges represent the sectors employed to derive the hyperbolic shape bestmatching the shock front geometry. The maximumaposteriori model is shown as a dashed line, while the dotted contours indicate the corresponding 95% credible interval. Further, the dashed and solid yellow circles denote respectively the ACA and ALMA fullwidthathalfmaximum fields of view. For reference, the upstream and downstream gas lie respectively west (right) and east (left) of the shock front. 

In the text 
Fig. 2.
Bivariate posterior density function (left panel) for the inferred Mach number ℳ and slope of the downstream pressure distribution α for a set of upstream pressure profiles and respective centroids. To facilitate comparison, the bestfitting ℳ estimates are also plotted (right panel). We consider three different cases for the underlying gNFW profile by setting the slopes equal to the values reported in Arnaud et al. (2010) for the coolcore (cc; red), morphologicallydisturbed (md; peach), or ensemble (up; blue) cluster samples. For each of them, the distribution centroid is then fixed to a number of different positions: far downstream and far upstream (right panel, left and midleft points), respectively to arbitrary distances of 10′ east and 3′ west of the Xrayderived shock nose coordinates; APEXSZ centroid (midright point; Halverson et al. 2009); peak of the κmap (right point; Clowe et al. 2006). In all the above cases, we assumed instantaneous shock heating of the electrons. The gray line in the right panel denotes the corresponding bestfitting ℳ reported in Sect. 3.1, while the darker and lighter bands the respective 68% and 95% credibility intervals. 

In the text 
Fig. 3.
Marginalized posterior distributions for the shock Mach number ℳ derived under the assumptions of instantaneous (red) and collisional (blue) electronion equilibration. The dashed and dotted lines indicates the median of the posterior distributions and the 68% credible intervals, respectively. 

In the text 
Fig. 4.
Dirty images of the raw (left), model (middle), and residual (right) ALMA+ACA interferometric data. They are generated by jointly gridding the ACA and ALMA data using a multifrequency natural weighting scheme. For reference, this provides a synthesized beam of 4.07″ × 3.01″ FWHM (PA 81°; bottomleft corner of right panel). We cut off the fields at the 0.2 gain level of the joint ALMA+ACA antenna pattern. To better highlight the largescale shock features, we avoid correcting for the primary beam attenuation and apply an additional 30 kλ taper. We note that the model subtraction is performed directly in Fourier space. The dashed line in the center and right panels indicate the inferred position of the shock front. Xray contours are overlaid on the left panel from Fig. 1. We also note the positive signal at the shock front is not due to an increment of the SZ signal; rather it is an artifact of the highpass filtering effects of ALMA+ACA (see also Fig. 5). 

In the text 
Fig. 5.
Comparison of SZ signal profiles (upper panel) and the corresponding electron pressure profiles (lower panel) across the shock nose. The vertical gray line denotes the position of the shock front. We note that the upper panel contains both unfiltered input model fits for the SZ signal and the corresponding filtered (observed) profiles (see legend for details). As discussed in Sect. 3.2, the instantaneous (dashdotted yellow) and collisional (solid red) shock models are indistinguishable after the spatial filtering of the interferometric ALMA+ACA observation, reflecting the fundamental limitation of ALMA+ACA to constrain any largescale (≳1′) component of the SZ signal. Again, we note that while we measure a decrement due to the SZ effect, the filtered (observed) profiles can exhibit both positive and negative excursions (analogous to the Gibbs phenomenon; see, e.g., Bracewell 1978). For comparison, we also report the input unfiltered model (solid yellow) for the instantaneous equilibration, and both the raw and filtered underlying gNFW profile (dashed and dotted blue lines, as noted in the legend). We also present the Xray expectation for both the filtered and raw SZ signal profiles (dashed and solid green lines) and corresponding pressure profiles that we would expect for the case of instantaneous equilibration, using the value ℳ_{X} = 2.74 ± 0.25 derived from fits to the Xray data. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.