Search

We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase., 21 pages, 14 figures. Accepted for publication by MNRAS

To experimentally search for signals of the Gardner transition, an active quasi-thermal granular glass is constructed using a monolayer of air-fluidized star-shaped granular particles. The pressure of the system is controlled by adjusting the tension exerted on an enclosing boundary. Velocity distributions of the internal particles and the scaling of the pressure, density, effective-temperature, and relaxation time are examined, demonstrating that the system has important features of a thermal system. Using a pressure-based quenching protocol that brings the system into deeper glassy states, signals of the Gardner transition are detected via cage size and separation order parameters for both particle positions and orientations, offering experimental evidence of the Gardner transition for a system with in low spatial dimensions that is quasi-thermal., Comment: 5 pages, 4 figures

Modeling collisionless magnetic reconnection rate is an outstanding challenge in basic plasma physics research. While the seemingly universal rate of an order $\mathcal{O}(0.1)$ is often reported in the low-$\beta$ regime, it is not clear how reconnection rate scales with a higher plasma $\beta$. Due to the complexity of the pressure tensor, the available reconnection rate model is limited to the low plasma-$\beta$ regime, where the thermal pressure is arguably negligible. However, the thermal pressure effect becomes important when $\beta \gtrsim \mathcal{O}(1)$. Using first-principle kinetic simulations, we show that both the reconnection rate and outflow speed drop as $\beta$ gets larger. A simple analytical framework is derived to take account of the self-generated pressure anisotropy and pressure gradient in the force-balance around the diffusion region, explaining the varying trend of key quantities and reconnection rates in these simulations with different $\beta$. The predicted scaling of the normalized reconnection rate is $\simeq \mathcal{O}(0.1/\sqrt{\beta_{i0}})$ in the high $\beta$ limit, where $\beta_{i0}$ is the ion $\beta$ of the inflow plasma., Comment: 19 pages, 7 figures, accepted for publication in ApJ

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons., Comment: 11 pages, 7 figures, A&A in press. This is one of three companion papers (including Ghirardini et al. and Ettori et al.) on the X-COP sample. X-COP data products are available for download on https://www.astro.unige.ch/xcop

We explore the relation between the star formation rate surface density ($\Sigma$SFR) and the interstellar gas pressure for nearby compact starburst galaxies. The sample consists of 17 green peas and 19 Lyman break analogs. Green peas are nearby analogs of Ly$\alpha$ emitters at high redshift and Lyman break analogs are nearby analogs of Lyman break galaxies at high redshift. We measure the sizes for green peas using Hubble Space Telescope Cosmic Origins Spectrograph (COS) NUV images with a spatial resolution of $\sim$ 0.05$^{''}$. We estimate the gas thermal pressure in HII regions by $P = N_{total}Tk{_B} \simeq 2n_{e}Tk{_B}$. The electron density is derived using the [SII] doublet at 6716,6731 \AA, and the temperature is calculated from the [OIII] lines. The correlation is characterized by $\Sigma$ SFR = 2.40$\times$10$^{-3\,}$M$_{\odot\,}$yr$^{-1\,}$kpc$^{-2}$$\left(\frac{P/k_{B}}{10^{4}cm^{-3}K}\right)^{1.33}$. Green peas and Lyman break analogs have high $\Sigma$SFR up to 1.2 M$_{\odot\,}$yr$^{-1\,}$kpc$^{-2}$ and high thermal pressure in HII region up to P/k$_B$ $\sim$10$^{7.2}{\rm\, K\, cm}^{-3}$. These values are at the highest end of the range seen in nearby starburst galaxies. The high gas pressure and the correlation, are in agreement with those found in star-forming galaxies at z $\sim$ 2.5. These extreme pressures are shown to be responsible for driving galactic winds in nearby starbursts. These outflows may be a crucial in enabling Lyman-$\alpha$ and Lyman-continuum to escape., Comment: 16 pages, 6 figures, 1 table; Accepted for publication in ApJ

Clumping and turbulence are expected to affect the matter accreted onto the outskirts of galaxy clusters. To determine their impact on the thermodynamic properties of Abell 2142 we perform an analysis of the X-ray temperature data from XMM-Newton via our SuperModel, a state-of-the-art tool for investigating the astrophysics of the intracluster medium already tested on many individual clusters (since Cavaliere et al. 2009). Using the gas density profile corrected for clumpiness derived by Tchernin et al. (2016), we find evidence for the presence of a nonthermal pressure component required to sustain gravity in the cluster outskirts of Abell 2142, that amounts to about 30% of the total pressure at the virial radius. The presence of the nonthermal component implies the gas fraction to be consistent with the universal value at the virial radius and the electron thermal pressure profile to be in good agreement with that inferred from the SZ data. Our results indicate that the presence of gas clumping and of a nonthermal pressure component are both necessary to recover the observed physical properties in the cluster outskirts. Moreover, we stress that an alternative method often exploited in the literature (included Abell 2142) to determine the temperature profile k_BT = P_e/n_e basing on a combination of the Sunyaev-Zel'dovich (SZ) pressure P_e and of the X-ray electron density n_e does not allow to highlight the presence of nonthermal pressure support in the cluster outskirts., Comment: 8 pages, 4 figures, Accepted by MNRAS

Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks., Comment: Accepted for publication in ApJ on Dec. 20th, 2016; 17 pages; 7 figures and 1 table

Non-thermal pressure in the intracluster gas has been found ubiquitously in numerical simulations, and observed indirectly. In this paper we develop an analytical model for intracluster non-thermal pressure in the virial region of relaxed clusters. We write down and solve a first-order differential equation describing the evolution of non-thermal velocity dispersion. This equation is based on insights gained from observations, numerical simulations, and theory of turbulence. The non-thermal energy is sourced, in a self-similar fashion, by the mass growth of clusters via mergers and accretion, and dissipates with a time-scale determined by the turnover time of the largest turbulence eddies. Our model predicts a radial profile of non-thermal pressure for relaxed clusters. The non-thermal fraction increases with radius, redshift, and cluster mass, in agreement with numerical simulations. The radial dependence is due to a rapid increase of the dissipation time-scale with radii, and the mass and redshift dependence comes from the mass growth history. Combing our model for the non-thermal fraction with the Komatsu-Seljak model for the total pressure, we obtain thermal pressure profiles, and compute the hydrostatic mass bias. We find typically 10% bias for the hydrostatic mass enclosed within $r_{500}$., Comment: 12 pages, 9 figures, published in MNRAS. Discussions and references added. A factor of 2 corrected in t_dyn (Fig. 2), definition of t_d (Eq. 3) changed accordingly

UV absorption studies with FUSE have observed H2 molecular gas in translucent and diffuse clouds. Observations of the 158 micron [C II] fine structure line with Herschel also trace the same H2 molecular gas in emission. We present [C II] observations along 27 lines of sight (LOSs) towards target stars of which 25 have FUSE H2 UV absorption. We detect [C II] emission features in all but one target LOS. For three Target LOSs, which are close to the Galactic plane, we also present position-velocity maps of [C II] emission observed by HIFI in on-the-fly spectral line mapping. We use the velocity resolved [C II] spectra towards the target LOSs observed by FUSE to identify C II] velocity components associated with the H2 clouds. We analyze the observed velocity integrated [C II] spectral line intensities in terms of the densities and thermal pressures in the H2 gas using the H2 column densities and temperatures measured by the UV absorption data. We present the H2 gas densities and thermal pressures for 26 target LOSs and from the [C II] intensities derive a mean thermal pressure in the range 6100 to 7700 K cm^-3 in diffuse H2 clouds. We discuss the thermal pressures and densities towards 14 targets, comparing them to results obtained using the UV absorption data for two other tracers CI and CO., Comment: 16 pages; 13 figures

We present a numerical solver for plasma dynamics simulations in Hall magnetohydrodynamic (HMHD) approximation in one, two and three dimensions. We consider both isotropic and anisotropic thermal pressure cases, where a general gyrotropic approximation is used. Both explicit energy conservation equation and general polytropic state equations are considered. The numerical scheme incorporates second-order Runge-Kutta advancing in time and Kurganov-Tadmor scheme with van Leer flux limiter for the approximation of fluxes. A flux-interpolated constrained-transport approach is used to preserve solenoidal magnetic field in the simulations. The implemented code is validated using several test problems previously described in the literature. Additionally, we propose a new validation method for HMHD codes based on solitary waves that provides a possibility of quantitative rigorous testing in nonlinear (large amplitude) regime as an extension to standard tests using small-amplitude whistler waves. Quantitative tests of accuracy and performance of the implemented code show the fidelity of the proposed approach., Accepted for publication in Journal of Computational Physics, 20 pages, 15 figures

Galaxy cluster mass determinations achieved using X-ray and Sunyaev-Zel'dovich data combined with the assumption of hydrostatic equilibrium are generally biased. The bias exists for two main reasons: non-thermal pressure forces are expected to contribute to the overall pressure balance and deviations from spherical symmetry and hydrostatic equilibrium can be present. In this paper, we use a sample of zoom-in hydrodynamical simulations of galaxy clusters to measure the magnitude of hydrostatic bias and the non-thermal contribution to the total pressure. We propose a new empirical model for non-thermal pressure based on our simulations that can be applied to observations. We show that our model can be successfully applied to remove most of the bias related to neglection of non-thermal pressure, which is usually not included in hydrostatic cluster mass profile reconstructions. The use of this model may significantly improve the calibration of cluster scaling relations that are a key tool for cluster cosmology., Comment: 10 pages, 7 figures, 2 tables, submitted to MNRAS

Non-thermal pressure in galaxy clusters leads to underestimation of the mass of galaxy clusters based on hydrostatic equilibrium with thermal gas pressure. This occurs even for dynamically relaxed clusters that are used for calibrating the mass-observable scaling relations. We show that the analytical model for non-thermal pressure developed in Shi & Komatsu 2014 can correct for this so-called 'hydrostatic mass bias', if most of the non-thermal pressure comes from bulk and turbulent motions of gas in the intracluster medium. Our correction works for the sample average irrespective of the mass estimation method, or the dynamical state of the clusters. This makes it possible to correct for the bias in the hydrostatic mass estimates from X-ray surface brightness and the Sunyaev-Zel'dovich observations that will be available for clusters in a wide range of redshifts and dynamical states., 9 pages, 8 figures, published in MNRAS

Galaxy evolution depends strongly on the physics of the interstellar medium (ISM). Motivated by the need to incorporate the properties of the ISM in cosmological simulations we construct a simple method to include the contribution of non-thermal components in the calculation of pressure of interstellar gas. In our method we treat three non-thermal components - turbulence, magnetic fields and cosmic rays - and effectively parametrize their amplitude. We assume that the three components settle into a quasi-steady-state that is governed by the star formation rate, and calibrate their magnitude and density dependence by the observed Radio-FIR correlation, relating synchrotron radiation to star formation rates of galaxies. We implement our model in single cell numerical simulation of a parcel of gas with constant pressure boundary conditions and demonstrate its effect and potential. Then, the non-thermal pressure model is incorporated into RAMSES and hydrodynamic simulations of isolated galaxies with and without the non-thermal pressure model are presented and studied. Specifically, we demonstrate that the inclusion of realistic non-thermal pressure reduces the star formation rate by an order of magnitude and increases the gas depletion time by as much. We conclude that the non-thermal pressure can prolong the star formation epoch and achieve consistency with observations without invoking artificially strong stellar feedback., Comment: 18 pages, 14 figures, accepted to MNRAS. Updated to match final version

Turbulent gas motion inside galaxy clusters provides a non-negligible non-thermal pressure support to the intracluster gas. If not corrected, it leads to a systematic bias in the estimation of cluster masses from X-ray and Sunyaev-Zel'dovich (SZ) observations assuming hydrostatic equilibrium, and affects interpretation of measurements of the SZ power spectrum and observations of cluster outskirts from ongoing and upcoming large cluster surveys. Recently, Shi & Komatsu developed an analytical model for predicting the radius, mass, and redshift dependence of the non-thermal pressure contributed by the kinetic random motions of intracluster gas sourced by the cluster mass growth. In this paper, we compare the predictions of this analytical model to a state-of-the-art cosmological hydrodynamics simulation. As different mass growth histories result in different non-thermal pressure, we perform the comparison on 65 simulated galaxy clusters on a cluster-by-cluster basis. We find an excellent agreement between the modelled and simulated non-thermal pressure profiles. Our results open up the possibility of using the analytical model to correct the systematic bias in the mass estimation of galaxy clusters. We also discuss tests of the physical picture underlying the evolution of intracluster non-thermal gas motions, as well as a way to further improve the analytical modeling, which may help achieve a unified understanding of non-thermal phenomena in galaxy clusters., 10 pages, 13 figures, accepted to MNRAS

Cosmological constraints from X-ray and microwave observations of galaxy clusters are subjected to systematic uncertainties. Non-thermal pressure support due to internal gas motions in galaxy clusters is one of the major sources of astrophysical uncertainties. Using a mass-limited sample of galaxy clusters from a high-resolution hydrodynamical cosmological simulation, we characterize the non-thermal pressure fraction profile and study its dependence on redshift, mass, and mass accretion rate. We find that the non-thermal pressure fraction profile is universal across redshift when galaxy cluster radii are defined with respect to the mean matter density of the universe instead of the commonly used critical density. We also find that the non-thermal pressure is predominantly radial, and the gas velocity anisotropy profile exhibits strong universality when galaxy cluster radii are defined with respect to the mean matter density of the universe. However, we find that the non-thermal pressure fraction is strongly dependent on the mass accretion rate of the galaxy cluster. We provide fitting formulae for the universal non-thermal pressure fraction and velocity anisotropy profiles of gas in galaxy clusters, which should be useful in modeling astrophysical uncertainties pertinent to using galaxy clusters as cosmological probes., 9 pages, 6 figures; Accepted by ApJ

I present and discuss a unified scheme for jet launching that is based on stochastic dissipation of the accretion disk kinetic energy, mainly via shock waves. In this scheme, termed thermally-launched jet model, the kinetic energy of the accreted mass is transferred to internal energy, e.g., heat or magnetic energy. The internal energy accelerates a small fraction of the accreted mass to high speeds and form jets. For example, thermal energy forms a pressure gradient that accelerates the gas. A second acceleration stage is possible wherein the primary outflow stretches magnetic field lines. The field lines then reconnect and accelerate small amount of mass to very high speeds. This double-stage acceleration process might form highly relativistic jets from black holes and neutron stars. The model predicts that detail analysis of accreting brown dwarfs that launch jets will show the mass accretion rate to be larger than 10^{-9}-10^{-8} Mo/year, which is higher than present claims in the literature., To appear in the proceedings of Star-disk interaction in young stars, Grenoble 2007, ed. J. Bouvier

We found equilibrium conditions for a self-gravitating toroidal vortex by taking into account thermal pressure. These conditions are shown to significantly differ from those for a disk and a sphere. The evolution of a thin vortex turns it into a compact vortex that loses mechanical stability for low masses at a polytropic index $\gamma < 4/3$ but retains stability for sufficiently high masses and densities determined by the velocity circulation in the vortex., Comment: 11 pages, 4 figures

We have used the Deep Survey telescope of EUVE to investigate shadows in the diffuse EUV/Soft X-Ray background cast by clouds in the interstellar medium. We confirm the existence of a shadow previously reported, and provide evidence for two new shadows. We used IRAS data to identify the clouds producing these shadows and to determine their optical depth to EUV radiation. The EUV-absorbing clouds are optically thick in the EUV, and all EUV emission detected in the direction of these shadows must be produced from material in front of the clouds. We obtained new optical data to determine the distance to these clouds. We use a new differential cloud technique to obtain the pressure of the interstellar medium. These results do not depend on any zero level calibration of the data. Our results provide evidence that the pressure of the hot interstellar gas is the same in three different directions in the local interstellar medium, and is at least 8 times higher than derived for the local cloud surrounding our Sun. This provides new evidence for large thermal pressure imbalances in the local ISM, and directly contradicts the basic assumption of thermal pressure equilibrium used in almost all present models of the interstellar medium., 13 pages, 13 figures

The merger of two neutron stars leaves behind a rapidly spinning hypermassive object whose survival is believed to depend on the maximum mass supported by the nuclear equation of state, angular momentum redistribution by (magneto-)rotational instabilities, and spindown by gravitational waves. The high temperatures (~5-40 MeV) prevailing in the merger remnant may provide thermal pressure support that could increase its maximum mass and, thus, its life on a neutrino-cooling timescale. We investigate the role of thermal pressure support in hypermassive merger remnants by computing sequences of spherically-symmetric and axisymmetric uniformly and differentially rotating equilibrium solutions to the general-relativistic stellar structure equations. Using a set of finite-temperature nuclear equations of state, we find that hot maximum-mass critically spinning configurations generally do not support larger baryonic masses than their cold counterparts. However, subcritically spinning configurations with mean density of less than a few times nuclear saturation density yield a significantly thermally enhanced mass. Even without decreasing the maximum mass, cooling and other forms of energy loss can drive the remnant to an unstable state. We infer secular instability by identifying approximate energy turning points in equilibrium sequences of constant baryonic mass parametrized by maximum density. Energy loss carries the remnant along the direction of decreasing gravitational mass and higher density until instability triggers collapse. Since configurations with more thermal pressure support are less compact and thus begin their evolution at a lower maximum density, they remain stable for longer periods after merger., 20 pages, 12 figures. Accepted for publication in ApJ

While clusters of galaxies are regarded as one of the most important cosmological probes, the conventional spherical modeling of the intracluster medium (ICM) and the dark matter (DM), and the assumption of strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure) are very approximate at best. Extending our previous works, we developed further a method to reconstruct for the first time the full three-dimensional structure (triaxial shape and principal axis orientation) of both DM and intracluster (IC) gas, and the level of non-thermal pressure of the IC gas. We outline an application of our method to the galaxy cluster Abell 383, taken as part of the CLASH multi-cycle treasury program, presenting results of a joint analysis of X-ray and strong lensing measurements. We find that the intermediate-major and minor-major axis ratios of the DM are 0.71+/-0.10 and 0.55+/-0.06, respectively, and the major axis of the DM halo is inclined with respect to the line of sight of 21.1+/-10.1 deg. The level of non-thermal pressure has been evaluated to be about 10% of the total energy budget. We discuss the implications of our method for the viability of the CDM scenario, focusing on the concentration parameter C and the inner slope of the DM gamma, since the cuspiness of dark-matter density profiles in the central regions is one of the critical tests of the cold dark matter (CDM) paradigm for structure formation: we measure gamma=1.02+/-0.06 on scales down to 25 Kpc, and C=4.76+/-0.51, values which are close to the predictions of the standard model, and providing further evidences that support the CDM scenario. Our method allows us to recover the three-dimensional physical properties of clusters in a bias-free way, overcoming the limitations of the standard spherical modelling and enhancing the use of clusters as more precise cosmological probes., Comment: MNRAS in press

A phenomenon analogous to the conical refraction widely known in the crystalooptics and crystaloacoustics is discovered for the magnetohydrodynamical waves in the collisionless plasma with anisotropic thermal pressure. Angle of the conical refraction is calculated for the medium under study which is predicted to be $18^{\circ}26^{\prime}$. Possible experimental corroborating of the discovered phenomenon is discussed., 6 pages, REVTeX, Accepted in Physics of Plasmas

The localized deposition of the energy of a laser pulse, as it ablates a solid target, introduces high thermal pressure gradients in the plasma. The thermal expansion of this laser-heated plasma into the ambient medium (ionized residual gas) triggers the formation of non-linear structures in the collision-less plasma. Here an electron-proton plasma is modelled with a particle-in-cell (PIC) simulation to reproduce aspects of this plasma expansion. A jump is introduced in the thermal pressure of the plasma, across which the otherwise spatially uniform temperature and density change by the factor 100. The electrons from the hot plasma expand into the cool one and the charge imbalance drags a beam of cool electrons into the hot plasma. This double layer reduces the electron temperature gradient. The presence of the low-pressure plasma modifies the proton dynamics compared to the plasma expansion into a vacuum. The jump in the thermal pressure develops into a primary shock. The fast protons, which move from the hot into the cold plasma in form of a beam, give rise to the formation of phase space holes in the electron and proton distributions. The proton phase space holes develop into a secondary shock that thermalizes the beam., Comment: 17 pages, 9 Figures, accepted for publication in plasma physics and controlled fusion

We present a spherically symmetric model for the origin and evolution of the temperature profiles in the hot plasma filling galaxy groups and clusters. We find that the gas in clusters is generically not isothermal, and that the temperature declines with radius at large distances from the cluster center (outside the core- and scale radii). This temperature profile is determined by the accretion history of the halo, and is not quantitatively well-described by a polytropic model. We explain quantitatively how the large-scale temperature gradient persists in spite of thermal conduction and convection. These results are a consequence of the cosmological assembly of clusters and cannot be reproduced with non-cosmological simulations of isolated halos. We show that the variation in halo assembly histories produces a ~10% scatter in temperature at fixed mass. On top of this scatter, conduction decreases the temperature of the gas near the scale radius in massive clusters, which may bias hydrostatic mass estimates inferred from x-ray and SZ observations. As an example application of our model profiles, we use mixing-length theory to estimate the turbulent pressure support created by the magnetothermal instability (MTI): in agreement with our earlier MHD simulations, we find that the convection produced by the MTI can provide ~5% non-thermal pressure support near r_500. The magnitude of this turbulent pressure support is likely to be non-monotonic in halo mass, peaking in ~10^14.5 M_sun halos., accepted for publication in MNRAS

We compare the gravitational potential profiles of the elliptical galaxies NGC 4486 (M87) and NGC 1399 (the central galaxy in the Fornax cluster) derived from X-ray and optical data. This comparison suggests that the combined contribution of cosmic rays, magnetic fields and micro-turbulence to the pressure is ~10% of the gas thermal pressure in the cores of NGC 1399 and M87, although the uncertainties in our model assumptions (e.g., spherical symmetry) are sufficiently large that the contribution could be consistent with zero. In the absence of any other form of non-thermal pressure support, these upper bounds translate into upper limits on the magnetic field of ~10-20 muG at a distance of 1'-2' from the centers of NGC1399 and M87. We show that these results are consistent with the current paradigm of cool cluster cores, based on the assumption that AGN activity regulates the thermal state of the gas by injecting energy into the intra-cluster medium. The limit of ~10-20% on the energy density in the form of relativistic protons applies not only to the current state of the gas, but essentially to the entire history of the intra-cluster medium, provided that cosmic ray protons evolve adiabatically and that their spatial diffusion is suppressed., Comment: Accepted for MNRAS. 19 pages; 14 figures; expanded version in response to comments from the referee

We compare X-ray hydrostatic and weak-lensing mass estimates for a sample of 12 clusters that have been observed with both XMM-Newton and Subaru. At an over-density of Δ = 500, we obtain 1 − MX/MWL = 0.01 ± 0.07 for the whole sample. We also divided the sample into undisturbed and disturbed sub-samples based on quantitative X-ray morphologies using asymmetry and fluctuation parameters, obtaining 1 − MX/MWL = 0.09 ± 0.06 and −0.06 ± 0.12 for the undisturbed and disturbed clusters, respectively. In addition to non-thermal pressure support, there may be a competing effect associated with adiabatic compression and/or shock heating which leads to overestimate of X-ray hydrostatic masses for disturbed clusters, for example, in the famous merging cluster A1914. Despite the modest statistical significance of the mass discrepancy, on average, in the undisturbed clusters, we detect a clear trend of improving agreement between MX and MWL as a function of increasing overdensity, MX/MWL = (0.908 ± 0.004) + (0.187 ± 0.010) · log10(Δ/500). We also examine the gas mass fractions, fgas = Mgas/MWL, finding that they are an increasing function of cluster radius, with no dependence on dynamical state, in agreement with predictions from numerical simulations.Overall, our results demonstrate that XMM-Newton and Subaru are a powerful combination for calibrating systematic uncertainties in cluster mass measurements., We acknowledge support from KICP in Chicago for hospitality, and thank our LoCuSS collaborators, especially Masahiro Takada and Keiichi Umetsu, for helpful comments on the manuscript. Y.Y.Z. thanks Massimo Meneghetti and Gabriel Pratt for useful discussion. Y.Y.Z. acknowledges support by the DFG through Emmy Noether Research Grant RE 1462/ 2, through Schwerpunkt Program 1177, and through project B6 “Gravitational Lensing and X-ray Emission by Non-Linear Structures” of Transregional Collaborative Research Centre TRR 33 The Dark Universe, and support by the German BMBF through the Verbundforschung under grant 50 OR 0601. This work is supported by a Grant-in-Aid for the COE Program “Exploring New Science by Bridging Particle-Matter Hierarchy” and G-COE Program “Weaving Science Web beyond Particle-Matter Hierarchy” in Tohoku University, funded by theMinistry of Education, Science, Sports and Culture of Japan. This work is, in part, supported by a Grant-in-Aid for Science Research in a Priority Area “Probing the Dark Energy through an Extremely Wide and Deep Survey with Subaru Telescope” (18072001) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan. N.O. is, in part, supported by a Grant-in-Aid from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (20740099). A.F. acknowledges support from BMBF/DLR under grant 50 OR 0207 and MPG, and was partially supported by a NASA grant NNX08AX46G to UMBC. G.P.S. acknowledges support from the Royal Society and STFC. D.P.M. acknowledges support provided by NASA through Hubble Fellowship grant HF-51259.01 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS 5-26555.

The utility of large Sunyaev Zel'dovich (SZ) surveys for determining cosmological parameters from cluster abundances is limited by the theoretical uncertainties in the integrated SZ-flux-to-mass relation, Y-M. We explore how non-thermal pressure and the anisotropic shape of the gas distribution of the intracluster medium (ICM) impacts Y-M scaling using a suite of SPH simulations of the cosmic web. We contrast results for models with different treatments of entropy injection and transport, varying radiative cooling, star formation and accompanying supernova feedback, cosmic rays, and energetic feedback from active galactic nuclei (AGN). We find that the gas kinetic-to-thermal pressure ratio from internal bulk motions depends on the cluster mass, and increases in the outer-cluster due to enhanced substructure, as does the asphericity of the ICM gas. With only a ~5-10% correction to projected (observable) ellipticities, we can infer the 3D ellipticities. Our simulated Y-M-slope roughly follows the self-similar prediction, except for a steepening due to a deficit of gas in lower mass clusters at low redshift in our AGN-feedback simulations. AGN feedback enhances slightly the overall Y-M-scatter, from ~11% to ~13%, a reflection of accretion history variations due to cluster merging. If we split the cluster system into lower, middle and upper bands of both P_kin/P_th and long-to-short axis ratio, we find a ~10% effect on Y-M. Identifying observable second parameters related to internal bulk flows and anisotropy for cluster-selection to minimize Y-M scatter in a "fundamental plane" would allow tighter cosmological parameter constraints., Comment: To be submitted to ApJ, 22 pages, 22 figures

We show how the fully resummed thermal pressure is rendered ultraviolet finite by standard zero-temperature renormalisation. The analysis is developed in a 6-dimensional scalar model that mimics QED and has $N$ flavours. The $N\to\infty$ limit of the model can be calculated completely. At a critical temperature, one of the degrees of freedom has vanishing screening mass like the transverse gauge bosons in four-dimensional finite-temperature perturbation theory. The renormalised nonperturbative interaction pressure of this model is evaluated numerically., 27 pages, plain tex, with 10 figures embedded using epsf

The central mass distribution of clusters of galaxies can be inferred from gravitationally lensed arcs with known redshifts. For the Abell cluster 2218, this method yields a core mass which is larger by a factor of $2.5 \pm 0.5$ than the value deduced from X-ray observations, under the assumptions that the gas is supported by thermal pressure and that the cluster is spherical. We show that a non-thermal pressure support is the most plausible explanation for this discrepancy. Such a pressure can be naturally provided by strong turbulence and equipartition magnetic fields ($\sim 50��$G) that are tangled on small spatial scales (, 12 pages, Postcript file

(Abridged) The main purpose of this paper is to consider the contribution of all three non-thermal components to total mass measurements of galaxy clusters: cosmic rays, turbulence and magnetic pressures. To estimate the thermal pressure we used public XMM-\textit{Newton} archival data of 5 Abell clusters. To describe the magnetic pressure, we assume a radial distribution for the magnetic field, $B(r) \propto \rho_{g}^{\alpha}$, to seek generality we assume $\alpha$ within the range of 0.5 to 0.9, as indicated by observations and numerical simulations. For the turbulent component, we assumed an isotropic pressure, $P_{\rm turb} = {1/3}\rho_{\rm g}(\sigma_{r}^{2}+\sigma_{t}^{2})$. We also consider the contribution of cosmic ray pressure, $P_{cr}\propto r^{-0.5}$. It follows that a consistent description for the non-thermal component could yield variation in mass estimates that vary from 10% up to $\sim$30%. We verified that in the inner parts of cool-core clusters the cosmic ray component is comparable to the magnetic pressure, while in non cool-core cluster the cosmic ray component is dominant. For cool-core clusters the magnetic pressure is the dominant component, contributing with more than 50% of total mass variation due to non-thermal pressure components. However, for non cool-core clusters, the major influence comes from the cosmic ray pressure that accounts with more than 80% of total mass variation due to non-thermal pressure effects. For our sample, the maximum influence of the turbulent component to total mass variation can be almost 20%. We show that this analysis can be regarded as a starting point for a more detailed and refined exploration of the influence of non-thermal pressure in the intra-cluster medium (ICM)., Comment: 14 pages, 4 figures, accepted for publication in A&A

A phenomenon analogous to the conical refraction well-known in the crystalooptics and crystaloacoustics is considered for the magnetohydrodynamical waves in a collisionless plasma with anisotropic thermal pressure. Imposing the most general (generalization of Tsiklauri, 1996, Phys. Plasmas, 3, 800) condition for the existence of the phenomenon, angle of the conical refraction is calculated which appeared to be dependent on the ratio of the Alfven velocity and sound speed measured in the perpendicular direction in respect to the external magnetic field. Feasible ways of experimental demonstration of the phenomenon are discussed and a novelty brought by the general consideration is outlined., Comment: 6 pages, RevTeX, to appear in Phys. Plasmas, October 1996 issue

A simple four-mode Bose-Hubbard model with intrinsic time scale separation can be considered as a paradigm for mesoscopic quantum systems in thermal contact. In our previous work we showed that in addition to coherent particle exchange, a novel slow collective excitation can be identified by a series of Holstein-Primakoff transformations. This resonant energy exchange mode is not predicted by linear Bogoliubov theory, and its frequency is sensitive to interactions among Bogoliubov quasi-particles; it may be referred to as a second Josephson oscillation, in analogy to the second sound mode of liquid Helium II. In this paper we will explore this system beyond the Gross-Pitaevskii mean field regime. We directly compare the classical mean field dynamics to the exact full quantum many-particle dynamics and show good agreement over a large range of the system parameters. The second Josephson frequency becomes imaginary for stronger interactions, however, indicating dynamical instability of the symmetric state. By means of a generalized quantum phase model for the full four-mode system, we then show that, in this regime, high-energy Bogoliubov quasiparticles tend to accumulate in one pair of sites, while the actual particles preferentially occupy the opposite pair. We interpret this as a simple model for thermal pressure.

We present a modelization of the non-thermal pressure, $P_{NT}$, and we apply it to the X-ray (and Sunayev-Zel'dovich) derived radial profiles of the X-COP galaxy clusters. We relate the amount of non-thermal pressure support to the hydrostatic bias, $b$, and speculate on how we can interpret this $P_{NT}$ in terms of the expected levels of turbulent velocity and magnetic fields. Current upper limits on the turbulent velocity in the intracluster plasma are used to build a distribution $\mathcal{N}(, Comment: 6 pages, 3 figures. Accepted for publication in A&A Letters

We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar trends are also seen for the azimuthal temperature scatter and non-thermal pressure fraction, both of which correlate with density fluctuations, with these values also increasing with redshift. However, in agreement with recent work, we find only a weak correlation between the clumping, or its proxies, and the hydrostatic mass bias. To reduce the effect of clumping in the projected profiles, we compute the azimuthal median following recent observational studies, and find this reduces the scatter in the bias. We also attempt to correct the cluster masses by using a non-thermal pressure term and find over-corrected mass estimates ($1-b=0.86$ to $1-b=1.15$) from 3D gas profiles but improved mass estimates ($1-b=0.75$ to $1-b=0.85$) from projected gas profiles, with the caveat of systematically increased scatter. We conclude that the cluster-averaged mass bias is minimised from applying a non-thermal pressure correction ($1-b=0.85$) with more modest reductions from selecting clusters that have low clumping ($1-b=0.79$) or are dynamically relaxed ($1-b=0.80$). However, the latter selection is most effective at minimising the scatter for individual objects. Such results can be tested with next generation X-ray missions equipped with high-resolution spectrometers such as Athena., Comment: 13 pages, 10 figures, accepted by MNRAS

We analyze the physical properties of gas in the circumgalactic medium (CGM) of 132 Milky Way (MW)-like central galaxies at $z=0$ from the cosmological magneto-hydrodynamical simulation TNG50, part of the IllustrisTNG project. The properties and abundance of CGM gas across the sample are diverse, and the fractional budgets of different phases (cold, warm, and hot), as well as neutral HI mass and metal mass, vary considerably. Over our stellar mass range of $10^{10.5} < M_\star / \rm{M}_\odot < 10^{10.9}$, radial profiles of gas physical properties from $0.15 < R\rm{ / R_{\rm 200c}} < 1.0$ reveal great CGM structural complexity, with significant variations both at fixed distance around individual galaxies, and across different galaxies. CGM gas is multi-phase: the distributions of density, temperature and entropy are all multimodal, while metallicity and thermal pressure distributions are unimodal; all are broad. We present predictions for magnetic fields in MW-like halos: a median field strength of $|B|\sim\,1μ$G in the inner halo decreases rapidly at larger distance, while magnetic pressure dominates over thermal pressure only within $\sim0.2 \times \rm{R_{200c}}$. Virial temperature gas at $\sim 10^6\,$K coexists with a sub-dominant cool, $< 10^5\,$K component in approximate pressure equilibrium. Finally, the physical properties of the CGM are tightly connected to the galactic star formation rate, in turn dependent on feedback from supermassive black holes (SMBHs). In TNG50, we find that energy from SMBH-driven kinetic winds generates high-velocity outflows ($\gtrsim 500-2000$ km/s), heats gas to super-virial temperatures ($> 10^{6.5-7}$ K), and regulates the net balance of inflows versus outflows in otherwise quasi-static gaseous halos., Minor typos fixed. Additional visuals at www.tng-project.org/ramesh22!

We analyze the internal cross-section structure and dynamics of inter-halo cosmic-web filaments that connect to massive high-$z$ galaxies. Our analysis is based on a high-resolution \texttt{AREPO} cosmological simulation zooming in on a volume encompassing three filaments of $\sim 1 \, \Mpc$ length feeding three massive haloes of $\sim 10^{12}\,\msun$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a strong cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \kpc$, within which the gas is in virial equilibrium in the potential well set by an isothermal dark-matter filament. The cross-section of the filament can be divided into three radial zones. In the outer ``thermal'' (\textbf{T}) zone, $r \geq 0.65 \, r_{\rm shock}$, the inward gravity and ram-pressure forces from the inflowing gas are over-balanced by outward thermal pressure forces, causing the inflowing gas to decelerate and the shock to expand outwards. In the intermediate ``vortex'' (\textbf{V}) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outward force is provided mainly by centrifugal forces associated with these vortices and, to a lesser extent, by global rotation and thermal pressure, while shear and turbulent forces that are also associated with the vortex structure add to the force inward. The inner ``stream'' (\textbf{S}) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core defining the cold streams that feed the galaxies, with $T\sim 3 \times 10^4 \, \Kdegree$ and $n_{\rm H}\sim 0.01 \cmc$. The core is formed by an isobaric cooling flow that results from post-shock cooling on a free-fall time, and is associated with a decrease in the outward forces, though the gas within the core exhibits both inflows and outflows. [abridged], 30 pages, 20 figures, submitted to MNRAS

We present an analytic model for the cool, $T \approx 10^4$ K, circumgalactic medium (CGM), describing the gas distribution, thermal and ionization state. Our model assumes (total) pressure equilibrium with the ambient warm/hot CGM, photoionization by the metagalactic radiation field, and allows for non-thermal pressure support, parametrized by the ratio of thermal pressures, $\eta = P_{\rm hot,th}/P_{\rm cool,th}$. We apply the model to the COS-Halos data set and find that a nominal model with $\eta = 3$, gas distribution out to $r \approx 0.6 R_{\rm vir}$, and $M_{\rm cool} = 3 \times 10^9~{\rm M_{\odot}}$, corresponding to a volume filling fraction of $f_{\rm V,cool} \approx 1\%$, reproduces the mean measured column densities of HI and low/intermediate metal ions (CII, CIII, SiII, SiIII, MgII). Variation of $\pm 0.5$ dex in the non-thermal pressure or gas mass encompasses $\sim 2/3$ of the scatter between objects. Our nominal model underproduces the measured CIV and SiIV columns, and we show these can be reproduced with (i) a cool phase with $M_{\rm cool} \approx 10^{10}~{\rm M_{\odot}}$ and $\eta \approx 5$, or (ii) an additional component at intermediate temperatures, of cooling or mixing gas, with $M \approx 1.5 \times 10^{10}~{\rm M_{\odot}}$ and occupying $\sim 1/2$ of the total CGM volume. For cool gas with $f_{\rm V,cool} \approx 1\%$ we provide an upper limit on the cloud sizes, $R_{\rm cl} \lesssim 0.5$ kpc. Our results suggest that for the average galaxy CGM, the mass and non-thermal support in the cool phase are lower than estimated in previous works, and extreme scenarios for galactic feedback and non-thermal support may not be necessary. We estimate the rates of cool gas depletion and replenishment, and find accretion onto the galaxy can be entirely offset by condensation, outflows, and IGM accretion, allowing $\dot{M}_{\rm cool}\sim0$ over long timescales., Comment: Submitted to ApJ, 17 pages + appendices, comments welcome

In many published articles isobaric behavior in DAC chambers has been reported On the other hand, a clear isochoric behavior in the chambers were reported. Thus, Isochoric or isobaric condition remain an open question. In isochoric condition, upon increasing the temperature, at each pressure, the examined sample and the PTM both are subject to increases of their volumes which is suppressed by the chamber finite volume, provoking an increase in the thermal pressure over the whole system. For this reason, the pressure scale reported in DAC experiments do not represent the actual pressure experienced by the sample in the cell. The different response of the PTMs to P,T changes is the reason for the variety of melting curves reported in the literature. It is concluded that the pressure scale in DAC experiments need a correction taking in to account the thermal pressure shift. Trusting the first principals DFT- Z methodology it is shown that the melting curve derived by this method should serve as an anchor for the pressure scale correction., 14 pages, 4 figures. arXiv admin note: text overlap with arXiv:2209.09836

Water photodissociation in the 114 - 144 nm UV range forms excited OH which emits at mid-infrared wavelengths via highly excited rotational lines. These lines have only been detected with Spitzer in several proto-planetary disks and shocks. Previous studies have shown they are a unique diagnostic for water photodissociation. Thanks to its high sensitivity and angular resolution, the James Webb Space Telescope (JWST) could be able to detect them in other environments such as interstellar Photo-Dissociation Regions (PDRs). In order to predict the emerging spectrum of OH, we use the Meudon PDR Code to compute the thermal and chemical structure of PDRs. The influence of thermal pressure ($P_{\rm th}/k$ = $n_{\rm H} T_{\rm K}$) and UV field strength on the integrated intensities, as well as their detectability with the JWST are studied in details. OH mid-IR emission is predicted to originate very close to the H$^0$/H$_2$ transition and is directly proportional to the column density of water photodissociated in that layer. Because neutral gas-phase formation of water requires relatively high temperatures ($T_{\rm K} \gtrsim 300~$K), the resulting OH mid-IR lines are primarily correlated with the temperature at this position, and are therefore brighter in regions with high pressure. This implies that these lines are predicted to be only detectable in strongly irradiated PDRs ($G_0^{\rm incident}$ $>$ 10$^3$) with high thermal pressure ($P_{\rm th}/k$ $\gtrsim$ 5$\times$10$^7$ K cm$^{-3}$). In the latter case, OH mid-IR lines are less dependent on the strength of the incident UV field. The detection in PDRs like the Orion bar, which should be possible, is also investigated. To conclude, OH mid-IR lines observable by JWST are a promising diagnostics for dense and strongly irradiated PDRs., Accepted on the 17/10/2022 by Astronomy & Astrophysics (A&A)

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly-resolved $L^*$ galaxies from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM ($\Delta x \lesssim 1$ kpc $h^{-1}$ and $M \simeq 200M_\odot$). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of $\lesssim5$ kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond $\sim 0.25R_{200}$, or $\sim50$ kpc at $z=0$. Within $\sim0.25R_{200}$, turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the non-equilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies even when averaging out small-scale variations., Comment: 36 pages, 24 figures. Accepted to ApJ. See https://foggie.science/cgm_forces.html for animated Figures 6 and 7

The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$., 20 pages; 21 figures; Substantial Revision; MNRAS in press

Deep Chandra and VLA imaging reveals a clear correlation between X-ray and radio emission on scales $\sim 100$~kpc in the Spiderweb radio galaxy at z=2.16. The X-ray emission associated with the extended radio source is likely dominated by inverse Compton up-scattering of cosmic microwave background photons by the radio emitting relativistic electrons. For regions dominated by high surface brightness emission, such as hot spots and jet knots, the implied magnetic fields are $\sim 50~\mu$G to $70~\mu$G. The non-thermal pressure is these brighter regions is then $\sim 9\times 10^{-10}$ dyne cm$^{-2}$, or three times larger than the non-thermal pressure derived assuming minimum energy conditions, and an order of magnitude larger than the thermal pressure in the ambient cluster medium. Assuming ram pressure confinement implies an average advance speed for the radio source of $\sim 2400$ km s$^{-1}$, and a source age of $\sim 3\times 10^7$ years. Considering the lower surface brightness, diffuse radio emitting regions, we identify an evacuated cavity in the Ly$\alpha$ emission coincident with the tail of the eastern radio lobe. Making reasonable assumptions for the radio spectrum, we find that the relativistic electrons and fields in the lobe are plausibly in pressure equilibrium with the thermal gas, and close to a minimum energy configuration. The radio morphology suggests that the Spiderweb is a high-$z$ example of the rare class of hybrid morphology radio sources (or HyMoRS), which we attribute to interaction with the asymmetric gaseous environment indicated by the Ly$\alpha$ emission., Comment: 8 figures, 25 pages, to appear in the ApJ

We present the properties of magnetic fields (B-fields) in two clumps (clump 1 and clump 2), located at the waist of the bipolar H II region Sh2-201, based on JCMT SCUBA-2/POL-2 observations of 850 $\mu$m polarized dust emission. We find that B-fields in the direction of the clumps are bent and compressed, showing bow-like morphologies, which we attribute to the feedback effect of the H II region on the surface of the clumps. Using the modified Davis-Chandrasekhar-Fermi method we estimate B-fields strengths of 266 $\mu$G and 65 $\mu$G for clump 1 and clump 2, respectively. From virial analyses and critical mass ratio estimates, we argue that clump 1 is gravitationally bound and could be undergoing collapse, whereas clump 2 is unbound and stable. We hypothesize that the interplay between thermal pressure imparted by the H II region, B-field morphologies, and the various internal pressures of the clumps (such as magnetic, turbulent, and gas thermal pressure), has the following consequences: (a) formation of clumps at the waist of the H II region; (b) progressive compression and enhancement of the B-fields in the clumps; (c) stronger B-fields will shield the clumps from erosion by the H II region and cause pressure equilibrium between the clumps and the H II region, thereby allowing expanding I-fronts to blow away from the filament ridge, forming bipolar H II regions; and (d) stronger B-fields and turbulence will be able to stabilize the clumps. A study of a larger sample of bipolar H II regions would help to determine whether our hypotheses are widely applicable., Comment: 25 pages, 10 figures, and 4 tables. Accepted for publication in ApJ

We constrain the cosmic ray (CR) population in the circumgalactic medium (CGM) of Milky Way by comparing the observations of absorption lines of OVIII ion with predictions from analytical models of CGM : precipitation (PP) and isothermal (IT) model. For a CGM in hydrostatic equilibrium, the introduction of CR suppresses thermal pressure, and affects the OVIII ion abundance. We explore the allowances given to the ratio of CR pressure to thermal pressure ($\rm{P}_{\rm{CR}}/\rm{P}_{\rm{th}}=\eta$), with varying boundary conditions, CGM mass content, photoionization by extragalactic ultraviolet background and temperature fluctuations. We find that the allowed maximum values of $\eta$ are : $\eta\lesssim10$ in the PP model and $\eta\lesssim6$ in the IT model. We also explore the spatial variation of $\eta$ : rising ($\eta=Ax$) or declining ($\eta=A/x$) with radius, where A is the normalization of the profiles. In particular, the models with declining ratio of CR to thermal pressure fare better than those with rising ratio with suitable temperature fluctuation (larger $\sigma_{\rm ln T}$ for PP and lower for IT). The declining profiles allow $A\lesssim8$ and $A\lesssim10$ in the case of IT and PP models, respectively, thereby accommodating a large value of $\eta \,(\simeq 200)$ in the central region, but not in the outer regions. These limits, combined with the limits derived from $\gamma$-ray and radio background, can be useful for building models of Milky Way CGM including CR population. However, the larger amount of CR can be packed in cold phase which may be one way to circumvent these constraints., Comment: 11 pages, 5 figures, 1 table Accepted for publication in ApJ on Apr 22, 2022

We combine the Santa-Cruz Semi-Analytic Model (SAM) for galaxy formation and evolution with the circumgalactic medium (CGM) model presented in Faerman et al. (2020) to explore the CGM properties of $L^{*}$ galaxies. We use the SAM to generate a sample of galaxies with halo masses similar to the Milky Way (MW) halo, $M_{\rm vir} \approx 10^{12}~{\rm M_{sun}}$, and find that the CGM mass and mean metallicity in the sample are correlated. We use the CGM masses and metallicities of the SAM galaxies as inputs for the FSM20 model, and vary the amount of non-thermal support. The density profiles in our models can be approximated by power-law functions with slopes in the range of $0.75 < a_n < 1.25$, with higher non-thermal pressure resulting in flatter distributions. We explore how the gas pressure, dispersion measure, OVI-OVIII column densities, and cooling rates behave with the gas distribution and total mass. We show that for CGM masses below $\sim 3 \times 10^{10}~{\rm M_{sun}}$, photoionization has a significant effect on the column densities of OVI and OVIII. The combination of different MW CGM observations favors models with similar fractions in thermal pressure, magnetic fields/cosmic rays, and turbulent support, and with $M_{\rm gas} \sim 3-10 \times 10^{10}~{\rm M_{sun}}$. The MW OVI column requires $t_{\rm cool}/t_{\rm dyn} \sim 4$, independent of the gas distribution. The AGN jet-driven heating rates in the SAM are enough to offset the CGM cooling, although exact balance is not required in star-forming galaxies. We provide predictions for the columns densities of additional metal ions - NV, NeVIII, and MgX., Submitted to ApJ. 19 pages, 15 figures (+ appendix)

Observations of the ICM in the outskirts of the Virgo cluster with Suzaku have found the gas mass fraction in the northern direction to be significantly above the expected level, indicating that there may be a very high level of gas clumping on small scales in this direction. Here, we explore the XMM-Newton data in the outskirts of Virgo, dividing it into a Voronoi tessellation to separate the bulk ICM component from the clumped ICM component. As the nearest galaxy cluster, Virgo's large angular extent allows the spatial scale of the tessellation to be much smaller than has been achieved using the same technique on intermediate redshift clusters, allowing us to probe gas clumping on the scales down to 5$\times$5 kpc. We find that the level of gas clumping in the outskirts to the north is relatively mild, ($\sqrt{C} < 1.1$), suggesting that our point-source detection procedure may have excluded a significant fraction of clumps. While correcting for clumping brings the gas mass fraction at $r_{200}$ into agreement with the universal gas mass fraction, the values outside $r_{200}$ remain significantly above it. This may suggest that non-thermal pressure support in the outskirts to the north is significant, and we find that a non-thermal pressure support at the level of 20 per cent of the total pressure outside $r_{200}$ can explain the high gas mass fraction to the north, which is in agreement with the range expected from simulations., 10 pages, 12 figures, 1 table. Accepted for publication in MNRAS

A certain appeal to the alpha model for turbulence and related viscosity in accretion disks was that one scales the Reynolds stresses simply on the thermal pressure, assuming that turbulence driven by a certain mechanism will attain a characteristic Mach number in its velocity fluctuations. Besides the notion that there are different mechanism driving turbulence and angular momentum transport in a disk, we also find that within a single instability mechanism, here the Vertical Shear Instability, stresses do not linearly scale with thermal pressure. Here we demonstrate in numerical simulations the effect of the gas temperature gradient and the thermal relaxation time on the average stresses generated in the non-linear stage of the instability. We find that the stresses scale with the square of the exponent of the radial temperature profile at least for a range of $d \log T /d \log R = [-0.5, -1]$, beyond which the pressure scale height varies too much over the simulation domain, to provide clear results. Stresses are also dependent on thermal relaxation times, provided they are longer than $10^{-3}$ orbital periods. The strong dependence of viscous transport of angular momentum on the local conditions in the disk (especially temperature, temperature gradient, and surface density/optical depth) challenges the ideas of viscosity leading to smooth density distributions, opening a route for structure (ring) formation and time variable mass accretion., Comment: 9 pages, 14 figures, accepted for publication in MNRAS

Due to the lower ionospheric thermal pressure and existence of the crustal magnetism at Mars, the Martian ionopause is expected to behave differently from the ionopause at Venus. We study the solar wind interaction and pressure balance at the ionopause of Mars using both in situ and remote sounding measurements from the MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) instrument on the Mars Express orbiter. We show that the magnetic pressure usually dominates the thermal pressure to hold off the solar wind at the ionopause at Mars, with only 13% of the cases where the ionospheric thermal pressure plays a more important role in pressure balance. This percentage at Venus, however, is up to 65%. We also find that the ionopause altitude at Mars decreases as the normal component of the solar wind dynamic pressure increases, similar to the altitude variation of the ionopauses at Venus. Moreover, our results show that the ionopause thickness at Mars and Venus is mainly determined by the ion gyromotion and is equivalent to about 5 ion gyroradii.

We present the detection of excited fine-structure energy levels of singly-ionized silicon and neutral carbon associated with the proximate damped Lyman-$\alpha$ system at $z_{\rm abs}=2.811$ towards \qso. This absorber has an apparent relative velocity that is inconsistent with the Hubble flow indicating motion along the line-of-sight towards the quasar, i.e., $z_{\rm abs}>z_{\rm em}$. We measure the metallicity of the system to be ${\rm [Zn/H]}=-0.68\pm 0.02$. Using the relative populations of the fine-structure levels of SiII and CI, as well as the populations of H$_2$ rotational levels, we constrain the physical conditions of the gas. We derive hydrogen number densities of $n_{\rm H}=190^{+70}_{-50}$ cm$^{-3}$ and $260^{+30}_{-20}$ cm$^{-3}$ in two velocity components where both CI and H$_2$ are detected. Taking into account the kinetic temperature in each component, $\sim 150$K, we infer high values of thermal pressure in the cold neutral medium probed by the observations. The strengths of the UV field in Draine's unit are $I_{\rm UV} = 10^{+5}_{-3}$ and $14^{+3}_{-3}$ in each of these two components, respectively. Such enhanced UV fluxes and thermal pressure compared to intervening DLAs are likely due to the proximity of the quasar. The typical size of the absorber is $\sim 10^4$ a.u. Assuming the UV flux is dominated by the quasar, we constrain the distance between the quasar and the absorber to be $\sim 150-200$ kpc. This favours a scenario where the absorption occurs in a companion galaxy located in the group where the quasar-host galaxy resides. This is in line with studies in emission that revealed the presence of several galaxies around the quasar., Comment: 11 pages, 5 figures, accepted for publication in MNRAS

Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure., 22 pages, 16 figures, accepted in MNRAS

Martian bow shocks, the solar wind interacting with an unmagnetized planet, are studied. We theoretically investigated how solar parameters, such as the solar wind dynamic pressure and the solar extreme ultraviolet (EUV) flux, influence the bow shock location, which is still currently not well understood. We present the formula for the location of the bow shock nose of the unmagnetized planet. The bow shock location, the sum of the ionopause location and bow shock standoff distance, is calculated in the gasdynamics approach. The ionopause location is determined using thermal pressure continuity, i.e., the solar wind thermal pressure equal to the ionospheric pressure, according to tangential discontinuity. The analytical formula of the ionopause nose location and the ionopause profile around the nose are obtained. The standoff distance is calculated using the empirical model. Our derived formula shows that the shock nose location is a function of the scale height of ionosphere, the dynamic pressure of the solar wind and the peak ionospheric pressure. The theoretical model implies that the shock nose location is more sensitive to the solar EUV flux than solar wind dynamics pressure. Further, we theoretically show that the bow shock location is proportional to the solar wind dynamic pressure to the power of negative C, where C is about the ratio of the ionospheric scale height to the distance between bow shock nose and the planet center. This theory matches the gasdynamics simulation and is consistent with the spacecraft measurement result by Mars Express [Hall, et al. (2016) J. Geophys. Res. Space Physics, 121, 11,474-11,494].