No Arabic abstract
We examine an ensemble of 48 simulated clusters to determine the effects of small-scale density fluctuations and large-scale substructure on X-ray measurements of the intracluster medium (ICM) mass. We measure RMS density fluctuations in the ICM which can be characterized by a mean mass-weighted clumping factor C = <rho^2>/<rho>^2 between 1.3 and 1.4 within a density contrast of 500 times the critical density. These fluctuations arise from the cluster history of accretion shocks and major mergers, and their presence enhances the clusters luminosity relative to the smooth case. We expect, therefore, that ICM mass measurements utilizing models which assume uniform density at a given radius carry a bias of order sqrt(C) = 1.16. We verify this result by performing ICM mass measurements on X-ray images of the simulations and finding the expected level of bias. The varied cluster morphologies in our ensemble also allow us to investigate the effects of departures from spherical symmetry on our measurements. We find that the presence of large-scale substructure does not further bias the resulting gas mass unless it is pronounced enough to produce a second peak in the image of at least 1% the maximum surface brightness. We analyze the subset of images with no secondary peaks and find a bias of 9% and a Gaussian random error of 4% in the derived mass.
We describe a method to extract resonant orbits from N-body simulations exploiting the fact that they close in a frame rotating with a constant pattern speed. Our method is applied to the N-body simulation of the Milky Way by Shen et al. (2010). This simulation hosts a massive bar, which drives strong resonances and persistent angular momentum exchange. Resonant orbits are found throughout the disc, both close to the bar itself and out to the very edges of the disc. Using Fourier spectrograms, we demonstrate that the bar is driving kinematic substructure even in the very outer parts of the disc. We identify two major orbit families in the outskirts of the disc that make significant contributions to the kinematic landscape, namely the m:l = 3:-2 and 1:-1 families resonating with the pattern speed of the bar. A mechanism is described that produces bimodal distributions of Galactocentric radial velocities at selected azimuths in the outer disc. It occurs as a result of the temporal coherence of particles on the 3:-2 resonant orbits, which causes them to arrive simultaneously at pericentre or apocentre. This resonant clumping, due to the in-phase motion of the particles through their epicycle, leads to both inward and outward moving groups which belong to the same orbital family and consequently produce bimodal radial velocity distributions. This is a possible explanation of the bimodal velocity distributions observed towards the Galactic anti-Centre by Liu et al. (2012). Another consequence is that transient overdensities appear and dissipate (in a symmetric fashion) on timescales equal to the their epicyclic period resulting in a periodic pulsing of the discs surface density.
(Abridged)The X-ray measurements of the ICM metallicity are becoming more frequent due to the availability of powerful X-ray telescope with excellent spatial and spectral resolutions. The information which can be extracted from the measurements of the alpha-elements, like Oxygen, Magnesium and Silicon with respect to the Iron abundance is extremely important to better understand the stellar formation and its evolutionary history. In this paper we investigate possible source of bias connected to the plasma physics when recovering metal abundances from X-ray spectra. To do this we analyze 6 simulated galaxy clusters processed through the new version of our X-ray MAp Simulator, which allows to create mock XMM-Newton EPIC MOS1 and MOS2 observations. By comparing the spectroscopic results to the input values we find that: i) Fe is recovered with high accuracy for both hot (T>3 keV) and cold (T<2 keV) systems; at intermediate temperatures, however, we find a systematic overestimate which depends on the number counts; ii) O is well recovered in cold clusters, while in hot systems its measure may overestimate by a factor up to 2-3; iii) Being a weak line, the measurement of Mg is always difficult; despite of this, for cold systems (T<2 keV) we do not find any systematic behavior, while for very hot systems (T>5 keV) the spectroscopic measurement may be strongly overestimated up to a factor of 4; iv) Si is well recovered for all the clusters in our sample. We investigate in detail the nature of the systematic effects and biases found. We conclude that they are mainly connected with the multi-temperature nature of the projected observed spectra and to the intrinsic limitation of the XMM-Newton EPIC spectral resolution that does not always allow to disentangle among the emission lines produced by different elements.
We have compared Monte Carlo photoionization models of H II regions with a uniform density distribution with models with the same central stars and chemical compositions but with 3-D hierarchical clumps. We compare the abundances of He, N, O, Ne, and S obtained from emission line strengths and [O III] and [N II] temperatures to those in our models. We consider stellar temperatures in the range 37.5 -- 45kK and ionizing luminosities from 10^{48} to 10^{51} photons/s. Clumped models have different ionic abundances than uniform. For hot stars, He^0/He^+ is 2 -- 3%, much larger than with uniform models. This amount of He I is independent of metallicity and so impacts the determination of the primordial abundance of He. The total abundances of O, Ne, and S obtained by the usual methods of analysis, using T([OIII) for high stages of ionization and T([NII]) for low, are about as accurate for clumped models as for uniform and within about 20% of the true values. If T([OIII]) is used for analyzing all ions, the derived (O/H) is 40 to 60% too large for cool stars but is good for hot stars. Uniform models have similar errors, so the clumping does not change the accuracy of abundance analysis. The physical causes of the ionic abundance errors are present in real nebulae. In clumped models, helium ionizing radiation from zones of high ionization (low He^0 and low UV opacity) can penetrate nearby regions near the edge of the ionized zone. This effect allows He^0 to absorb more stellar photons than in uniform or radially symmetrical geometries. In turn, these absorptions compete with O+, etc., for those energetic stellar photons.
Gravitational encounters between small-scale dark matter substructure and cold stellar streams in the Milky Way halo lead to density perturbations in the latter, making streams an effective probe for detecting dark matter substructure. The Pal 5 stream is one such system for which we have some of the best data. However, Pal 5 orbits close to the center of the Milky Way and has passed through the Galactic disk many times, where its structure can be perturbed by baryonic structures such as the Galactic bar and giant molecular clouds (GMCs). In order to understand how these baryonic structures affect Pal 5s density, we present a detailed study of the effects of the Galactic bar, spiral structure, GMCs, and globular clusters on the Pal 5 stream. We estimate the effect of each perturber on the stream density by computing its power spectrum and comparing it to the power induced by a CDM-like population of dark matter subhalos. We find that the bar and GMCs can each individually create power that is comparable to the observed power on large scales, leaving little room for dark matter substructure, while spirals are subdominant on all scales. On degree scales, the power induced by the bar is small, but GMCs create small-scale density variations that are similar in amplitude to the dark-matter induced variations but otherwise indistinguishable from it. These results demonstrate that Pal 5 is a poor system for constraining the dark matter substructure fraction and that observing streams further out in the halo will be necessary to confidently detect dark matter subhalos.
(abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. The effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps, are examined. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming TE chemistry, such as HCN in O-rich outflows, H$_2$O in C-rich outflows, and NH$_3$ in both. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O- and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H$_2$O in C-rich outflows and of NH$_3$ in O- and C-rich outflows are however underpredicted.