Galaxies undergoing ram-pressure stripping: the influence of the bulge on morphology and star formation rate

We investigate the influence of stellar bulges on the star formation and morphology of disc galaxies that suffer from ram pressure. Several tree-SPH (smoothed particle hydrodynamics) simulations have been carried out to study the dependence of the star formation rate on the mass and size of a stellar bulge. In addition, different strengths of ram pressure and different alignments of the disc with respect to the intra-cluster medium (ICM) are applied. As claimed in previous works, when ram pressure is acting on a galaxy, the star formation rate (SFR) is enhanced and rises up to four times with increasing ICM density compared to galaxies that evolve in isolation. However, a bulge suppresses the SFR when the same ram pressure is applied. Consequently, fewer new stars are formed because the SFR can be lowered by up to 2 M_sun/yr. Furthermore, the denser the surrounding gas, the more inter-stellar medium (ISM) is stripped. While at an ICM density of 10^-28 g/cm^3 about 30% of the ISM is stripped, the galaxy is almost completely (more than 90%) stripped when an ICM density of 10^-27 g/cm^3 is applied. But again, a bulge prevents the stripping of the ISM and reduces the amount being stripped by up to 10%. Thereby, fewer stars are formed in the wake if the galaxy contains a bulge. The dependence of the SFR on the disc tilt angle is not very pronounced. Hereby a slight trend of decreasing star formation with increasing inclination angle can be determined. Furthermore, with increasing disc tilt angles, less gas is stripped and therefore fewer stars are formed in the wake. Reducing the disc gas mass fraction results in a lower SFR when the galaxies evolve in vacuum. On the other hand, the enhancement of the SFR in case of acting ram pressure is less pronounced with increasing gas mass fraction. Moreover, the fractional amount of stripped gas does not depend on the gas mass fraction.Comment: 11 pages, 18 figure

Disk Galaxy Formation in a LambdaCDM Universe

We describe hydrodynamical simulations of galaxy formation in a Lambda cold dark matter (CDM) cosmology performed using a subresolution model for star formation and feedback in a multiphase interstellar medium (ISM). In particular, we demonstrate the formation of a well-resolved disk galaxy. The surface brightness profile of the galaxy is exponential, with a B-band central surface brightness of 21.0 mag arcsec^-2 and a scale-length of R_d = 2.0 h^-1 kpc. We find no evidence for a significant bulge component. The simulated galaxy falls within the I-band Tully-Fisher relation, with an absolute magnitude of I = -21.2 and a peak stellar rotation velocity of V_rot=121.3 km s^-1. While the total specific angular momentum of the stars in the galaxy agrees with observations, the angular momentum in the inner regions appears to be low by a factor of ~2. The star formation rate of the galaxy peaks at ~7 M_sun yr^-1 between redshifts z=2-4, with the mean stellar age decreasing from \~10 Gyrs in the outer regions of the disk to ~7.5 Gyrs in the center, indicating that the disk did not simply form inside-out. The stars exhibit a metallicity gradient from 0.7 Z_sun at the edge of the disk to 1.3 Z_sun in the center. Using a suite of idealized galaxy formation simulations with different models for the ISM, we show that the effective pressure support provided by star formation and feedback in our multiphase model is instrumental in allowing the formation of large, stable disk galaxies. If ISM gas is instead modeled with an isothermal equation of state, or if star formation is suppressed entirely, growing gaseous disks quickly violate the Toomre stability criterion and undergo catastrophic fragmentation.Comment: 14 pages, 12 figures, LaTex (emulateapj.cls), submitted to ApJ, high resolution images available at http://www-cfa.harvard.edu/~brobertson/papers/galaxy

Exploring Galaxy Formation Models and Cosmologies with Galaxy Clustering

Using N-body simulations and galaxy formation models, we study the galaxy stellar mass correlation and the two-point auto-correlation. The simulations are run with cosmological parameters from the WMAP first, third and seven year results, which mainly differ in the perturbation amplitude of \sigma_{8}. The stellar mass of galaxies are determined using either a semi-analytical galaxy formation model or a simple empirical abundance matching method. Compared to the SDSS DR7 data at z=0 and the DEEP2 results at z=1, we find that the predicted galaxy clusterings from the semi-analytical model are higher than the data at small scales, regardless of the adopted cosmology. Conversely, the abundance matching method predicts good agreement with the data at both z=0 and z=1 for high \sigma_8 cosmologies (WMAP1 & WMAP7), but the predictions from a low \sigma_8 cosmology (WMAP3) are significantly lower than the data at z=0. We find that the excess clustering at small-scales in the semi-analytical model mainly arises from satellites in massive haloes, indicating that either the star formation is too efficient in low-mass haloes or tidal stripping is too inefficient at high redshift. Our results show that galaxy clustering is strongly affected by the models for galaxy formation, thus can be used to constrain the baryonic physics. The weak dependence of galaxy clustering on cosmological parameters makes it difficult to constrain the WMAP1 and WMAP7 cosmologies.Comment: 8 pages, 7 figures. Accepted to MNRA

Early Structure Formation and Reionization in a Cosmological Model with a Running Primordial Power Spectrum

(abridged) We study high redshift structure formation and reionization in a LCDM universe under the assumption that the spectral power index of primordial density fluctuations is a function of length scale. We adopt a particular formulation of the running spectral index (RSI) model as suggested by the recent WMAP data. While early structure forms hierarchically in the RSI model, the reduced power on small scales causes a considerable delay in the formation epoch of low mass (~ 10^6 Msun) ``mini-halos'' compared to the LCDM model. The extremely small number of gas clouds in the RSI model indicates that reionization is initiated later than z<15, generally resulting in a smaller total Thomson optical depth than in the LCDM model. By carrying out radiative transfer calculations, we also study reionization by stellar populations formed in galaxies. Even with a top-heavy intial mass function representing an early population of massive stars and/or an extraordinarily high photon emission rate from galaxies, the total optical depth can only be as large as tau ~ 0.1 for reasonable models of early star-formation. The RSI model is thus in conflict with the large Thomson optical depth inferred by the WMAP satellite.Comment: Version accepted by ApJ. Visualizations are shown at http://cfa-www.harvard.edu/cpac/Reion/stars.htm

Metal enrichment of the intra-cluster medium over a Hubble time for merging and relaxed galaxy clusters

We investigate the efficiency of galactic mass loss, triggered by ram-pressure stripping and galactic winds of cluster galaxies, on the chemical enrichment of the intra-cluster medium (ICM). We combine N-body and hydrodynamic simulations with a semi-numerical galaxy formation model. By including simultaneously different enrichment processes, namely ram-pressure stripping and galactic winds, in galaxy-cluster simulations, we are able to reproduce the observed metal distribution in the ICM. We find that the mass loss by galactic winds in the redshift regime z>2 is ~10% to 20% of the total galactic wind mass loss, whereas the mass loss by ram-pressure stripping in the same epoch is up to 5% of the total ram-pressure stripping mass loss over the whole simulation time. In the cluster formation epochs z<2 ram-pressure stripping becomes more dominant than galactic winds. We discuss the non-correlation between the evolution of the mean metallicity of galaxy clusters and the galactic mass losses. For comparison with observations we present two dimensional maps of the ICM quantities and radial metallicity profiles. The shape of the observed profiles is well reproduced by the simulations in the case of merging systems. In the case of cool-core clusters the slope of the observed profiles are reproduced by the simulation at radii below ~300 kpc, whereas at larger radii the observed profiles are shallower. We confirm the inhomogeneous metal distribution in the ICM found in observations. To study the robustness of our results, we investigate two different descriptions for the enrichment process interaction.Comment: 11 pages, 13 figures, accepted for publication in A&A, high resolution version can be found at <http://astro.uibk.ac.at/~wolfgang/kapferer.pdf

Origin of chemically distinct discs in the Auriga cosmological simulations

The stellar disc of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom-in simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the [α/Fe]–[Fe/H] plane. We find two main pathways to creating this dichotomy, which operate in different regions of the galaxies: (a) an early (z > 1) and intense high-[α/Fe] star formation phase in the inner region (R ≲ 5 kpc) induced by gas-rich mergers, followed by more quiescent low-[α/Fe] star formation; and (b) an early phase of high-[α/Fe] star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process (b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process (a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-[α/Fe] sequence stars. We discuss the implication that both processes occurred in the Milky Way

Metal enrichment of the intra-cluster medium by thermally and cosmic-ray driven galactic winds

We investigate the efficiency and time-dependence of thermally and cosmic ray driven galactic winds for the metal enrichment of the intra-cluster medium (ICM) using a new analytical approximation for the mass outflow. The spatial distribution of the metals are studied using radial metallicity profiles and 2D metallicity maps of the model clusters as they would be observed by X-ray telescopes like XMM-Newton. Analytical approximations for the mass loss by galactic winds driven by thermal and cosmic ray pressure are derived from the Bernoulli equation and implemented in combined N-body/hydrodynamic cosmological simulations with a semi-analytical galaxy formation model. Observable quantities like the mean metallicity, metallicity profiles, and 2D metal maps of the model clusters are derived from the simulations. We find that galactic winds alone cannot account for the observed metallicity of the ICM. At redshift $z=0$ the model clusters have metallicities originating from galactic winds which are almost a factor of 10 lower than the observed values. For massive, relaxed clusters we find, as in previous studies, a central drop in the metallicity due to a suppression of the galactic winds by the pressure of the ambient ICM. Combining ram-pressure stripping and galactic winds we find radial metallicity profiles of the model clusters which agree qualitatively with observed profiles. Only in the inner parts of massive clusters the observed profiles are steeper than in the simulations. Also the combination of galactic winds and ram-pressure stripping yields too low values for the ICM metallicities. The slope of the redshift evolution of the mean metallicity in the simulations agrees reasonably well with recent observations.Comment: 9 pages, 6 figures, accepted by A&

Simulating Black Hole White Dwarf Encounters

The existence of supermassive black holes lurking in the centers of galaxies and of stellar binary systems containing a black hole with a few solar masses has been established beyond reasonable doubt. The idea that black holes of intermediate masses ($\sim 1000$ \msun) may exist in globular star clusters has gained credence over recent years but no conclusive evidence has been established yet. An attractive feature of this hypothesis is the potential to not only disrupt solar-type stars but also compact white dwarf stars. In close encounters the white dwarfs can be sufficiently compressed to thermonuclearly explode. The detection of an underluminous thermonuclear explosion accompanied by a soft, transient X-ray signal would be compelling evidence for the presence of intermediate mass black holes in stellar clusters. In this paper we focus on the numerical techniques used to simulate the entire disruption process from the initial parabolic orbit, over the nuclear energy release during tidal compression, the subsequent ejection of freshly synthesized material and the formation process of an accretion disk around the black hole.Comment: 9 pages, 4 figures, Computer Physics Communications, accepted; movie can be found at http://www.faculty.iu-bremen.de/srosswog/; reference correcte

Substructure of the galaxy clusters in the REXCESS sample: observed statistics and comparison to numerical simulations

We study the substructure statistics of a representative sample of galaxy clusters by means of two currently popular substructure characterisation methods, power ratios and centroid shifts. We use the 31 clusters from the REXCESS sample, compiled from the southern ROSAT All-Sky cluster survey REFLEX with a morphologically unbiased selection in X-ray luminosity and redshift, all of which have been reobserved with XMM-Newton. We investigate the uncertainties of the substructure parameters and examine the dependence of the results on projection effects, finding that the uncertainties of the parameters can be quite substantial. Thus while the quantification of the dynamical state of individual clusters with these parameters should be treated with extreme caution, these substructure measures provide powerful statistical tools to characterise trends of properties in large cluster samples. The centre shift parameter, w, is found to be more sensitive in general. For the REXCESS sample neither the occurence of substructure nor the presence of cool cores depends on cluster mass. There is a significant anti-correlation between the existence of substantial substructure and cool cores. The simulated clusters show on average larger substructure parameters than the observed clusters, a trend that is traced to the fact that cool regions are more pronounced in the simulated clusters, leading to stronger substructure measures in merging clusters and clusters with offset cores. Moreover, the frequency of cool regions is higher in the simulations than in the observations, implying that the description of the physical processes shaping cluster formation in the simulations requires further improvement.Comment: Mauscript submitted to Astronomy and Astrophysics, 20 figure

Simulation techniques for cosmological simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.Comment: 42 pages, 16 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 12; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke