Feedback from central supermassive blackholes is often invoked to explain the low star formation rates in massive galaxies at the centers of galaxy clusters. However, the detailed physics of the coupling of the injected feedback energy with the intracluster medium is still unclear. Using high-resolution magnetohydrodynamic cosmological simulations of galaxy cluster formation, we investigate the role of anisotropic thermal conduction in shaping the thermodynamic structure of clusters, and, in particular, in modifying the impact of black hole feedback. Stratified anisotropically conducting plasmas are formally always unstable, and thus more prone to mixing, an expectation borne out by our results. The increased mixing efficiently isotropizes the injected feedback energy which in turn significantly improves the coupling between the feedback energy and the intracluster medium. This facilitates an earlier disruption of the cool core, reduces the star formation rate by more than an order of magnitude, and results in earlier quenching despite an overall lower amount of feedback energy injected into the cluster core. With conduction, the metallicity gradients and dispersions are lowered, aligning them better with observational constraints. These results highlight the important role of thermal conduction in establishing and maintaining quiescence of massive galaxies.
We analyze the mass, temperature, metal enrichment, and OVI abundance of the circumgalactic medium (CGM) around $z\sim 0.2$ galaxies of mass $10^9 M_\odot M_\bigstar < 10^{11.5} M_\odot$ in the Illustris simulation. Among star-forming galaxies, the mass, temperature, and metallicity of the CGM increase with stellar mass, driving an increase in the OVI column density profile of $\sim 0.5$ dex with each $0.5$ dex increase in stellar mass. Observed OVI column density profiles exhibit a weaker mass dependence than predicted: the simulated OVI abundance profiles are consistent with those observed for star-forming galaxies of mass $M_\bigstar = 10^{10.5-11.5} M_\odot$, but underpredict the observed OVI abundances by $\gtrsim 0.8$ dex for lower-mass galaxies. We suggest that this discrepancy may be alleviated with additional heating of the abundant cool gas in low-mass halos, or with increased numerical resolution capturing turbulent/conductive mixing layers between CGM phases. Quenched galaxies of mass $M_\bigstar = 10^{10.5-11.5} M_\odot$ are found to have 0.3-0.8 dex lower OVI column density profiles than star-forming galaxies of the same mass, in qualitative agreement with the observed OVI abundance bimodality. This offset is driven by AGN feedback, which quenches galaxies by heating the CGM and ejecting significant amounts of gas from the halo. Finally, we find that the inclusion of the central galaxy's radiation field may enhance the photoionization of the CGM within $\sim 50$ kpc, further increasing the predicted OVI abundance around star-forming galaxies.
In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport of cosmic rays in the form of diffusion or streaming not only affects the gas dynamics but also rearranges the magnetic fields themselves. The coupled dynamics of magnetic fields and cosmic rays can thus impact the formation and evolution of galaxies and the thermal evolution of galaxy clusters in critical ways. Numerically studying these effects requires solvers for anisotropic diffusion that are accurate, efficient, and robust, requirements that have proved difficult to be satisfied in practice. Here, we present an anisotropic diffusion solver on an unstructured moving mesh that is conservative, does not violate the entropy condition, allows for semi-implicit time integration with individual timesteps, and only requires solving a single linear system of equations per timestep. We apply our new scheme to a large number of test problems and show that it works as well or better than previous implementations. Finally, we demonstrate for a numerically demanding simulation of the formation of an isolated disc galaxy that our local time-stepping scheme reproduces the results obtained with global time-stepping at a fraction of the computational cost.
We study how X-rays from stellar binary systems and the hot intracluster medium (ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel implementation of gas cooling in the moving-mesh hydrodynamics code AREPO. X-rays from stellar binaries do not affect cooling at all as their emission spectrum is too hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple well with gas in the temperature range 10^4-10^6 K. Idealized simulations show that the hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters because of the high virial temperature ( > 10^7 K), making the interaction between the gas and incident photons very ineffective. Satellite galaxies in cluster environments, on the other hand, experience a high radiation flux due to the emission from the host halo. Low-mass satellites (< 10^12 M⊙) in particular have virial temperatures that are exactly in the regime where the effect of the radiation field is maximal. Idealized simulations of satellite galaxies including only the effect of host halo radiation (no ram pressure stripping or tidal effects) fields show a drastic reduction in the amount of cool gas formed (~40 per cent) on a short time-scale of about 0.5 Gyr. A galaxy merger simulation including all the other environmental quenching mechanisms, shows about 20 per cent reduction in the stellar mass of the satellite and about ~30 per cent reduction in star formation rate after 1 Gyr due to the host hot halo radiation field. These results indicate that the hot halo radiation fields potentially play an important role in quenching galaxies in cluster environments.
We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the AREPO code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limited such that the total flux is both locally conservative and extremum preserving. The approach makes use of harmonic averaging points and a simple, robust interpolation scheme that works well for strong heterogeneous and anisotropic diffusion problems. Moreover, the required discretization stencil is small. Efficient fully implicit and semi-implicit time integration schemes are also implemented. We perform several numerical tests that evaluate the stability and accuracy of the scheme, including applications such as point explosions with heat conduction and calculations of convective instabilities in conducting plasmas. The new implementation is suitable for studying important astrophysical phenomena, such as the conductive heat transport in galaxy clusters, the evolution of supernova remnants, or the distribution of heat from black hole-driven jets into the intracluster medium.
We present a dynamical classification system for galaxies based on the shapes of their circular velocity curves (CVCs). We derive the CVCs of 40 SAURON and 42 CALIFA galaxies across Hubble sequence via a full line-of-sight integration as provided by solutions of the axisymmetric Jeans equations. We use Principal Component Analysis (PCA) applied to the circular curve shapes to find characteristic features and use a k-means classifier to separate circular curves into classes. This objective classification method identifies four different classes, which we name Slow-Rising (SR), Flat (F), Sharp-Peaked (SP) and Round-Peaked (RP) circular curves. SR-CVCs are mostly represented by late-type spiral galaxies (Scd-Sd) with no prominent spheroids in the central parts and slowly rising velocities; F-CVCs span almost all morphological types (E,S0,Sab,Sb-Sbc) with flat velocity profiles at almost all radii; SP-CVCs are represented by early-type and early-type spiral galaxies (E,S0,Sb-Sbc) with prominent spheroids and sharp peaks in the central velocities. RP-CVCs are represented by only two morphological types (E,Sa-Sab) with prominent spheroids, but RP-CVCs have much rounder peaks in the central velocities than SP-CVCs. RP-CVCs are typical for high-mass galaxies, while SR-CVCs are found for low-mass galaxies. Intermediate-mass galaxies usually have F-CVCs and SP-CVCs. Circular curve classification presents an alternative to typical morphological classification and may be more tightly linked to galaxy evolution.
We study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high-resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas-rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localized to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation.
We use semi-analytic models and cosmological merger trees to provide the initial conditions for multimerger numerical hydrodynamic simulations, and exploit these simulations to explore the effect of galaxy interaction and merging on star formation (SF). We compute numerical realizations of 12 merger trees from z = 1.5 to 0. We include the effects of the large hot gaseous halo around all galaxies, following recent observations and predictions of galaxy formation models. We find that including the hot gaseous halo has a number of important effects. First, as expected, the star formation rate on long time-scales is increased due to cooling of the hot halo and refuelling of the cold gas reservoir. Secondly, we find that interactions do not always increase the SF in the long term. This is partially due to the orbiting galaxies transferring gravitational energy to the hot gaseous haloes and raising their temperature. Finally, we find that the relative size of the starburst, when including the hot halo, is much smaller than previous studies showed. Our simulations also show that the order and timing of interactions are important for the evolution of a galaxy. When multiple galaxies interact at the same time, the SF enhancement is less than when galaxies interact in series. All these effects show the importance of including hot gas and cosmologically motivated merger trees in galaxy evolution models.
The distribution of galaxy morphological types is a key test for models of galaxy formation and evolution, providing strong constraints on the relative contribution of different physical processes responsible for the growth of the spheroidal components. In this paper, we make use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupting galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., based on results from high-resolution hydrodynamical simulations, and we show that these new implementations reduce the efficiency of bulge formation through mergers. In addition, we compare our model results with a variety of observational measurements of the fraction of spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present uncertainties in the data represent an important limitation to our understanding of spheroid formation. Our results indicate that the main tension between theoretical models and observations does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather from the distribution of galaxies of different morphological types, as a function of their stellar mass.
We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high-redshift galaxies including the stellar mass-halo mass relation, the galaxy stellar mass function (GSMF) at low masses and the number density evolution of low-mass galaxies. The poor match of M⋆-Mh and the GSMF at high masses indicates that supernova feedback is insufficient to limit star formation in these haloes. At z = 0, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around milky way mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernova feedback, plays a major role in regulating the properties of low-mass galaxies at high redshift.
We present a first study of the effect of local photoionizing radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionizing radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm as used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionization feedback through the whole history of a galaxy's formation. The simulation of a Milky Way-like galaxy using the local photoionization model forms ~40 per cent less stars than a simulation that only includes a standard uniform background UV field. The local photoionization model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionizing sources is significant and should not be ignored in models of galaxy formation.
We investigate the idea that the interaction of dark matter (DM) sub-halos with the gaseous disks of galaxies can be the origin for the observed holes and shells found in their neutral hydrogen (HI) distributions. We use high-resolution hydrodynamic simulations to show that pure DM sub-halos impacting a galactic disk are not able to produce holes; on the contrary, they result in high-density regions in the disk. However, sub-halos containing a small amount of gas (a few percent of the total DM mass of the sub-halo) are able to displace the gas in the disk and form holes and shells. The sizes and lifetimes of these holes depend on the sub-halo gas mass, density, and impact velocity. A DM sub-halo, of mass 10^8 M⊙ and a gas mass fraction of ~3%, is able to create a kiloparsec-scale hole with a lifetime similar to those observed in nearby galaxies. We also register an increase in the star formation rate at the rim of the hole, again in agreement with observations. Even though the properties of these simulated structures resemble those found in observations, we find that the number of predicted holes (based on mass and orbital distributions of DM halos derived from cosmological N-body simulations) falls short compared to the observations. Only a handful of holes are produced per gigayear. This leads us to conclude that DM halo impact is not the major channel through which these holes are formed.
The heat flux in rotating Rayleigh-Bénard convection in a fluid of Prandtl number Pr=0.1 enclosed between free-slip top and bottom boundaries is investigated using direct numerical simulation in a wide range of Rayleigh numbers (104≤Ra≤108) and Taylor numbers (0≤Ta≤108). The Nusselt number Nu scales with the Rayleigh number Ra as Ra^β with β=2/7 for values of Nu greater than a critical value Nuc, which occurs for Ta/Ra∼1. The exponent β is not universal for Nu < Nuc (for Ta/Ra > 1) but a function of Ta showing a minimum for some intermediate value of Ta. The critical Nusselt number Nuc and the corresponding critical Rossby number Roc scale with Ta as Ta^0.277±0.001 and Ta^−0.015±0.003, respectively.
We calculate the effects of frame dragging on the Galactic-Center stars. Assuming the stars are only slightly relativistic, we derive an approximation to the Kerr metric, which turns out to be a weak-field Schwarzschild metric plus a frame-dragging term. By numerically integrating the resulting geodesic equations, we compute the effect on Keplerian elements and the kinematics. We find that the kinematic effect at pericenter passage is proportional to (a(1 - e^2))^-2. For known Galactic-center stars it is of order 10 m/s. If observed, this would provide a measurement of the spin of the black hole.