Forming a large disc galaxy from a z < 1 major merger

Using high-resolution SPH simulations in a fully cosmological Λ cold dark matter context, we study the formation of a bright disc-dominated galaxy that originates from a 'wet' major merger at   z = 0.8  . The progenitors of the disc galaxy are themselves disc galaxies that formed from early major mergers between galaxies with blue colours. A substantial thin stellar disc grows rapidly following the last major merger and the present-day properties of the final remnant are typical of early-type spiral galaxies, with an i -band bulge-to-disc ratio ∼0.65, a disc scalelength of 7.2 kpc,   g − r = 0.5 mag  , an H  i linewidth ( W ₂₀/2) of 238 km s⁻¹ and total magnitude   i =−22.4  . The key ingredients for the formation of a dominant stellar disc component after a major merger are (i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, (ii) supernova feedback that is able to partially suppress star formation during mergers and (iii) relative fading of the spheroidal component. The gas fraction of the progenitors' discs does not exceed 25 per cent at   z < 3  , emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of discs and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disc galaxies by the high likelihood of interactions and mergers for galaxy-sized haloes at relatively low z .     

The pattern speed of the spiral galaxy NGC 1253

We present Very Large Array H  i observations of the gas-rich, interacting spiral galaxies, NGC 1253/1253A (Arp 279). The larger of the two galaxies, NGC 1253, has a very pronounced H  i ring and well-defined spiral structure. The velocity structure of the H  i data shows a sudden change at the position of the spiral arms; we identify this change as evidence of a strong spiral shock and hence proceed to estimate the pattern speed, Ω_p, of the spiral arms in NGC 1253. Assuming that the pattern speed is constant across the disc our derived value places the outer Lindblad resonance (OLR) at the position of the observed H  i ring. As an accumulation of gas is expected at the OLR when this falls within the disc of a galaxy this agreement provides independent support for the derived value of Ω_p.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Modelling gaseous and stellar kinematics in the disc galaxies NGC 772, 3898 and 7782

We present V -band surface photometry and major-axis kinematics of stars and ionized gas of three early-type spiral galaxies, namely NGC 772, 3898 and 7782. For each galaxy we present a self-consistent Jeans model for the stellar kinematics, adopting the light distribution of bulge and disc derived by means of a two-dimensional parametric photometric decomposition. This allows us to investigate the presence of non-circular gas motions, and derive the mass distribution of luminous and dark matter in these objects.
NGC 772 and 7782 have apparently normal kinematics with the ionized gas tracing the gravitational equilibrium circular speed. This is not true in the innermost region (| r |≲8 arcsec) of NGC 3898, where the ionized gas is rotating more slowly than the circular velocity predicted by dynamical modelling. This phenomenon is common in the bulge-dominated galaxies for which dynamical modelling enables us to make the direct comparison between the gas velocity and the circular speed, and it poses questions about the reliability of galaxy mass distributions derived by the direct decomposition of the observed ionized-gas rotation curve into the contributions of luminous and dark matter.     

Parameter tests within cosmological simulations of galaxy formation

Numerical simulations of galaxy formation require a number of parameters. Some of these are intrinsic to the numerical integration scheme (e.g., the time-step), while others describe the physical model (e.g., the gas metallicity). In this paper we present results of a systematic exploration of the effects of varying a subset of these parameters on simulations of galaxy formation. We use N -body and 'Smoothed Particle Hydrodynamics' techniques to follow the evolution of cold dark matter and gas in a small volume. We compare a fiducial model with 24 different simulations, in which one parameter at a time is varied, focusing on properties such as the relative fraction of hot and cold gas, and the abundance and masses of galaxies. We find that for reasonable choices of numerical values, many parameters have relatively little effect on the galaxies, with the notable exception of the parameters that control the resolution of the simulation and the efficiency with which gas cools.     

Formation of giant low surface brightness galaxies through disc instability

It is shown that the giant low surface brightness galaxies (GLSBs), characterized by a large but diffuse disc component, can result from ordinary spiral galaxies through dynamical evolution. Numerical simulations indicate that the formation of a bar in a gravitationally unstable disc with high surface density induces non-circular motions and radial mixing of disc matter, leading to the flattening of the disc density profile. The resulting decrease in the disc central surface brightness is ∼1.5 magnitude, while the disc scalelength is nearly doubled, transforming a typical high surface brightness galaxy to a GSLB. This scenario seems promising especially for the GSLBs possessing a significant bulge, which are difficult to incorporate into the traditional Hubble sequence. Namely, because this disc transmutation can operate even if a moderate bulge component exists, the GSLBs with a bulge are argued to have resulted from the high surface brightness galaxies which had already possessed a bulge. The current picture naturally explains other observed characteristics of the GSLBs as well, including the propensity for having grand-design spiral arms and a bar, a high incidence of active nuclei, and galaxy environments.     

How galaxies lose their angular momentum

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disc galaxies that are too compact with respect to the observations. High-resolution N -body/SPH simulations in a cosmological context are presented including cold gas and dark matter (DM). A halo with quiet merging activity since redshift   z ∼ 3.8  and with a high-spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disc. We show that the gas and the DM have similar specific angular momenta until a merger event occurs at   z ∼ 2  with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and dynamical friction processes falls quickly into the centre. In contrast, gas infall through small subclumps or accretion does not lead to catastrophic angular momentum loss. In fact, a new extended disc begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disc galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered starburst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected, late-type galaxies could form with a dominant extended disc component, resulting from late infall, a small bulge-to-disc ratio and a low baryon fraction, in agreement with observations.     

Supersonic motions of galaxies in clusters

We study motions of galaxies in galaxy clusters formed in the concordance Λ cold dark matter cosmology. We use high-resolution cosmological simulations that follow the dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing the motions of galaxies and the properties of intracluster gas in a sample of eight simulated clusters at z = 0, we study the velocity dispersion profiles of the dark matter, gas and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ≈1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move slightly faster than the dark matter particles. The magnitude of the velocity bias,   b _v ≈ 1.1  , is, however, smaller than the bias estimated for subhaloes in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ≈20–30 per cent of the gas sound speed. These random motions provide about 10 per cent of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.     

The angular momentum content of dwarf galaxies: new challenges for the theory of galaxy formation

We compute the specific angular momentum distributions for a sample of low-mass disc galaxies observed by Swaters. We compare these distributions to those of dark matter haloes obtained by Bullock et al. from high-resolution N -body simulations of structure formation in a ΛCDM universe. We find that although the disc mass fractions are significantly smaller than the universal baryon fraction, the total specific angular momenta of the discs are in good agreement with those of dark matter haloes. This suggests that discs form out of only a small fraction of the available baryons, but yet manage to draw most of the available angular momentum. In addition we find that the angular momentum distributions of discs are clearly distinct from those of the dark matter; discs lack predominantly both low and high specific angular momenta. Understanding these findings in terms of a coherent picture for disc formation is challenging. Cooling, feedback and stripping, which are the main mechanisms to explain the small disc mass fractions found, seem unable to simultaneously explain the angular momentum distributions of the discs. In fact, it seems that the baryons that make up the discs must have been born out of angular momentum distributions that are clearly distinct from those of ΛCDM haloes. However, the dark and baryonic mass components experience the same tidal forces, and it is therefore expected that they should have similar angular momentum distributions. Therefore, understanding the angular momentum content of disc galaxies remains an important challenge for our picture of galaxy formation.     

Maximum feedback and dark matter profiles of dwarf galaxies

The observed rotation curves of dark matter-dominated dwarf galaxies indicate low-density cores, contrary to the predictions of CDM models. A possible solution of this problem involves stellar feedback. A strong baryonic wind driven by vigorous star formation can remove a large fraction of the gas, causing the dark matter to expand. Using both numerical and analytical techniques, we explore the maximum effect of the feedback with an instantaneous removal of the gaseous disc. The energy input depends on the compactness of the disc, hence the specific angular momentum of the disc. For the plausible cosmological parameters and a wide range of the disc angular momenta, the feedback is insufficient to destroy the central halo cusp, while the inner density is lowered only by a modest factor of 2 to 6. Any realistic modelling of the feedback would have even lesser impact on dark matter. We find that no star formation effect can resolve the problems of CDM cusps.     

Dwarf galaxy rotation curves and the core problem of dark matter haloes

The standard cold dark matter (CDM) model has recently been challenged by the claim that dwarf galaxies have dark matter haloes with constant-density cores, whereas CDM predicts haloes with steeply cusped density distributions. Consequently, numerous alternative dark matter candidates have recently been proposed. In this paper we scrutinize the observational evidence for the incongruity between dwarf galaxies and the CDM model. To this end, we analyse the rotation curves of 20 late-type dwarf galaxies studied by Swaters. Taking the effects of beam smearing and adiabatic contraction into account, we fit mass models to these rotation curves with dark matter haloes with different cusp slopes, ranging from constant-density cores to r ⁻² cusps. Even though the effects of beam smearing are small for these data, the uncertainties in the stellar mass-to-light ratio and the limited spatial sampling of the halo's density distribution hamper a unique mass decomposition. Consequently, the rotation curves in our sample cannot be used to discriminate between dark haloes with constant-density cores and r ⁻¹ cusps. We show that the dwarf galaxies analysed here are consistent with CDM haloes in a ΛCDM cosmology, and that there is thus no need to abandon the idea that dark matter is cold and collisionless. However, the data are also consistent with any alternative dark matter model that produces dark matter haloes with central cusps less steep than r ^−1.5. In fact, we argue that based on existing H  i rotation curves alone, at best weak limits can be obtained on cosmological parameters and/or the nature of the dark matter. In order to make progress, rotation curves with higher spatial resolution and independent measurements of the mass-to-light ratio of the disc are required.     

Intrinsic correlation of galaxy shapes: implications for weak lensing measurements

Weak gravitational lensing is now established as a powerful method to measure mass fluctuations in the universe. It relies on the measurement of small coherent distortions of the images of background galaxies. Even low-level correlations in the intrinsic shapes of galaxies could however produce a significant spurious lensing signal. These correlations are also interesting in their own right, since their detection would constrain models of galaxy formation. Using     haloes found in N -body simulations, we compute the correlation functions of the intrinsic ellipticity of spiral galaxies assuming that the disc is perpendicular to the angular momentum of the dark matter halo. We also consider a simple model for elliptical galaxies, in which the shape of the dark matter halo is assumed to be the same as that of the light. For deep lensing surveys with median redshifts ∼1, we find that intrinsic correlations of ∼10⁻⁴ on angular scales     are generally below the expected lensing signal, and contribute only a small fraction of the excess signals reported on these scales. On larger scales we find limits to the intrinsic correlation function at a level ∼10⁻⁵, which gives a (model-dependent) range of separations for which the intrinsic signal is about an order of magnitude below the ellipticity correlation function expected from weak lensing. Intrinsic correlations are thus negligible on these scales for dedicated weak lensing surveys. For wider but shallower surveys such as SuperCOSMOS, APM and SDSS, we cannot exclude the possibility that intrinsic correlations could dominate the lensing signal. We discuss how such surveys could be used to calibrate the importance of this effect, as well as study spin–spin correlations of spiral galaxies.     

Dynamical response to supernova-induced gas removal in spiral galaxies with dark matter halo

We investigate the dynamical response, in terms of disc size and rotation velocity, to mass loss by supernovae in the evolution of spiral galaxies. A thin baryonic disc having the Kuzmin density profile embedded in a spherical dark matter halo having a density profile proposed by Navarro, Frenk & White is considered. For the purpose of comparison, we also consider the homogeneous and   r ⁻¹  profiles for dark matter in a truncated spherical halo. Assuming for simplicity that the dark matter distribution is not affected by mass-loss from discs and the change of baryonic disc matter distribution is homologous, we evaluate the effects of dynamical response in the resulting discs. We found that the dynamical response only for an adiabatic approximation of mass-loss can simultaneously account for the rotation velocity and disc size as observed particularly in dwarf spiral galaxies, thus reproducing the Tully–Fisher relation and the size versus magnitude relation over the full range of magnitude. Furthermore, we found that the mean specific angular momentum in discs after the mass-loss becomes larger than that before the mass-loss, suggesting that the mass-loss would most likely occur from the central disc region where the specific angular momentum is low.     

Numerical modelling of the vertical structure and dark halo parameters in disc galaxies

The non‐linear dynamics of bending instability and vertical structure of a galactic stellar disc embedded into a spherical halo are studied with N‐body numerical modelling. Development of the bending instability in stellar galactic disc is considered as the main factor that increases the disc thickness. Correlation between the disc vertical scale height and the halo‐to‐disc mass ratio is predicted from the simulations. The method of assessment of the spherical‐to‐disc mass ratio for edge‐on spiral galaxies with a small bulge is considered. Modelling of eight edge‐on galaxies: NGC 891, NGC 4738, NGC 5170, UGC 6080, UGC 7321, UGC 8286, UGC 9422 and UGC 9556 is performed. Parameters of stellar discs, dark haloes and bulges are estimated. The lower limit of the dark‐to‐luminous mass ratio in our galaxies is of the order of one within the limits of their stellar discs. The dark haloes dominate by mass in the galaxies with very thin stellar discs (NGC 5170, UGC 7321 and UGC 8286) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spurs and feathering in spiral galaxies

We present smoothed particle hydrodynamic (SPH) simulations of the response of gas discs to a spiral potential. These simulations show that the commonly observed spurs and feathering in spiral galaxies can be understood as being due to structures present in the spiral arms that are sheared by the divergent orbits in a spiral potential. Thus, dense molecular cloud-like structures generate the perpendicular spurs as they leave the spiral arms. Subsequent feathering occurs as spurs are further sheared into weaker parallel structures as they approach the next spiral passage. Self-gravity of the gas is not included in these simulations, stressing that these features are purely due to the hydrodynamics in spiral shocks. Instead, a necessary condition for this mechanism to work is that the gas need be relatively cold (1000 K or less) in order that the shock is sufficient to generate structure in the spiral arms, and such structure is not subsequently smoothed by the gas pressure.     

Dissipationless mergers of elliptical galaxies and the evolution of the fundamental plane

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter haloes and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter haloes [standard Navarro, Frenk & White (NFW) profile versus adiabatically contracted NFW profile], the global stellar to dark matter mass ratio and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy mass.     

The onset of warps in Spitzer observations of edge-on spiral galaxies

We analyse warps in the nearby edge-on spiral galaxies observed in the Spitzer /Infrared Array Camera (IRAC) 4.5-μm band. In our sample of 24 galaxies, we find evidence of warp in 14 galaxies. We estimate the observed onset radii for the warps in a subsample of 10 galaxies. The dark matter distribution in each of these galaxies are calculated using the mass distribution derived from the observed light distribution and the observed rotation curves. The theoretical predictions of the onset radii for the warps are then derived by applying a self-consistent linear response theory to the obtained mass models for six galaxies with rotation curves in the literature. By comparing the observed onset radii to the theoretical ones, we find that discs with constant thickness can not explain the observations; moderately flaring discs are needed. The required flaring is consistent with the observations. Our analysis shows that the onset of warp is not symmetric in our sample of galaxies. We define a new quantity called the onset-asymmetry index and study its dependence on galaxy properties. The onset asymmetries in warps tend to be larger in galaxies with smaller disc scalelengths. We also define and quantify the global asymmetry in the stellar light distribution, that we call the edge-on asymmetry in edge-on galaxies. It is shown that in most cases the onset asymmetry in warp is actually anticorrelated with the measured edge-on asymmetry in our sample of edge-on galaxies and this could plausibly indicate that the surrounding dark matter distribution is asymmetric.     

On the origin of warps and the role of the intergalactic medium

There is still no consensus as to what causes galactic discs to become warped. Successful models should account for the frequent occurrence of warps in quite isolated galaxies, their amplitude as well as the observed azimuthal and vertical distributions of the H  i layer. Intergalactic accretion flows and intergalactic magnetic fields may bend the outer parts of spiral galaxies. In this paper we consider the viability of these non-gravitational torques to take the gas off the plane. We show that magnetically generated warps are clearly flawed because they would wrap up into a spiral in less than two or three galactic rotations. The inclusion of any magnetic diffusivity to dilute the wrapping effect causes the amplitude of the warp to damp. We also consider the observational consequences of the accretion of an intergalactic plane-parallel flow at infinity. We have computed the amplitude and warp asymmetry in the accretion model, for a disc embedded in a flattened dark matter halo, including self-consistently the contribution of the modes with azimuthal wavenumbers   m = 0  and   m = 1  . Since the m = 0 component, giving a U-shaped profile, is not negligible compared to the m = 1 component, this model predicts quite asymmetric warps, maximum gas displacements on the two sides in the ratio 3 : 2 for the preferred Galactic parameters, and the presence of a fraction ∼3.5 per cent of U-shaped warps, at least. The azimuthal dependence of the moment transfer by the ram pressure would produce a strong asymmetry in the thickness of the H  i layer and asymmetric density distributions in z , in conflict with observational data for the warp in our Galaxy and in external galaxies. The amount of accretion that is required to explain the Galactic warp would give gas scaleheights in the far outer disc that are too small. We conclude that accretion of a flow with no net angular momentum cannot be the main and only cause of warps.     

The structure of spiral galaxies — II. Near-infrared properties of spiral arms

We have imaged a sample of 45 face-on spiral galaxies in the K band, to determine the morphology of the old stellar population, which dominates the mass in the disc. The K -band images of the spiral galaxies have been used to calculate different characteristics of the underlying density perturbation such as arm strengths, profiles and cross-sections, and spiral pitch angles. Contrary to expectations, no correlation was found between arm pitch angle and Hubble type, and combined with previous results this leads us to conclude that the morphology of the old stellar population bears little resemblance to the optical morphology used to classify galaxies. The arm properties of our galaxies seem inconsistent with predictions from the simplest density wave theories, and some observations, such as variations in pitch angle within galaxies, seem hard to reconcile even with more complex modal theories. Bars have no detectable effect on arm strengths for the present sample. We have also obtained B -band images of three of the galaxies. For these galaxies we have measured arm cross-sections and strengths, to investigate the effects of disc density perturbations on star formation in spiral discs. We find that B -band arms lead K -band arms and are narrower than K -band arms, apparently supporting predictions made by the large-scale shock scenario, although the effects of dust on B -band images may contribute towards these results.     

Simulations of spiral galaxies with an active potential: molecular cloud formation and gas dynamics

We describe simulations of the response of a gaseous disc to an active spiral potential. The potential is derived from an N -body calculation and leads to a multi-armed time-evolving pattern. The gas forms long spiral arms typical of grand-design galaxies, although the spiral pattern is asymmetric. The primary difference from a grand-design spiral galaxy, which has a consistent two-/four-armed pattern, is that instead of passing through the spiral arms, gas generally falls into a developing potential minimum and is released only when the local minimum dissolves. In this case, the densest gas is coincident with the spiral potential, rather than offset as in the grand-design spirals. We would therefore expect no offset between the spiral shock and star formation, and no obvious corotation radius. Spurs which occur in grand-design spirals when large clumps are sheared off leaving the spiral arms, are rare in the active, time-evolving spiral reported here. Instead, large branches are formed from spiral arms when the underlying spiral potential is dissolving due to the N -body dynamics. We find that the molecular cloud mass spectrum for the active potential is similar to that for clouds in grand-design calculations, depending primarily on the ambient pressure rather than the nature of the potential. The largest molecular clouds occur when spiral arms collide, rather than by agglomeration within a spiral arm.