The formation of the Hubble sequence of disc galaxies: the effects of early viscous evolution

We investigate a model of disc galaxies whereby viscous evolution of the gaseous disc drives material inwards to form a protobulge. We start from the standard picture of disc formation through the settling of gas into a dark halo potential well, with the disc initially coming into centrifugal equilibrium with detailed conservation of angular momentum. We derive generic analytic solutions for the disc–halo system after adiabatic compression of the dark halo, with free choice of the input virialized dark halo density profile and of the specific angular momentum distribution. We derive limits on the final density profile of the halo in the central regions. Subsequent viscous evolution of the disc is modelled by a variation of the specific angular momentum distribution of the disc, providing analytic solutions to the final disc structure. The assumption that the viscous evolution time-scale and the star formation time-scale are similar leads to predictions of the properties of the stellar components. Focusing on small 'exponential' bulges, i.e., ones that may be formed through a disc instability, we investigate the relationship between the assumed initial conditions, such as halo 'formation', or assembly, redshift z _f, spin parameter λ , baryonic fraction F , and final disc properties such as global star formation time-scale, gas fraction, and bulge-to-disc ratio. We find that the present properties of discs, such as the scalelength, are compatible with a higher initial formation redshift if the redistribution by viscous evolution is included than if it is ignored. We also quantify the dependence of final disc properties on the ratio F λ , thus including the possibility that the baryonic fraction varies from galaxy to galaxy, as perhaps may be inferred from the observations.     

Stability of a finite disc under the influence of a spherical halo

We have studied the stability of finite gaseous discs, against large-scale perturbations, under the influence of spherical, massive haloes. A surface-density distribution consistent with the observed spiral-tracer profiles in disc galaxies is considered for the disc. We find that growing eigenmodes with both ‘trailing’ and ‘leading’ spirals exist in ‘cold’ discs for a wide range of values of the halo mass and its radius. The amplification rates of the unstable modes reduce as the ratio of the mass of the halo to the mass of the disc is increased. A uniform halo is not very effective towards stabilizing the disc against these modes. The results from the present study are consideredvis-a-vis previous studies on the global modes of self-gravitating discs.     

The Monitor project: JW 380 – a 0.26-, 0.15-M⊙, pre-main-sequence eclipsing binary in the Orion nebula cluster

We present a general recipe for constructing N -body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N -body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.     

On the mass of M31

Recent work by several groups has established the properties of the dwarf satellites to M31. We reexamine the reported kinematics of this group employing a fresh technique we have developed previously. By calculating the distribution of a χ statistic (which we define in the paper) for the M31 system, we conclude that the total mass (disc plus halo) of the primary is unlikely to be as great as that of our own Milky Way. In fact the χ distribution for M31 indicates that, like NGC 3992, it does not have a massive halo. In contrast, the analysis of the satellites of NGC 1961 and NGC 5084 provides strong evidence for massive haloes surrounding both 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.     

The edge-on spiral gravitational lens B1600+434

We present new observations of the gravitational lens (GL) system B1600+434, strongly suggesting that the lens is an edge-on spiral galaxy. These observations are used to constrain the mass model of the system, in particular the oblateness and velocity dispersion of the dark matter halo around the lensing galaxy. From an analytical model we find a lower limit on the halo oblateness q _halo=( c/a )_ρ≳0.4; more detailed numerical models give a lower limit of q _halo≳0.5. We determine an average halo velocity dispersion of σ_halo=190±15 km s⁻¹ over all non-singular isothermal elliptical (NIE) halo models. Constraining the models to larger and more massive discs decreases this average by only 10 km s⁻¹. A lower limit of σ_halo≳150 km s⁻¹ is found, even for disc masses larger than the mass inside the Einstein radius. This lower limit indicates the need for a massive dark matter halo, contributing at least half of the mass inside the Einstein radius. Time-delay calculations give (54±3)/ h ₅₀ d for the NIE halo model and (70±4)/ h ₅₀ d for the modified Hubble profile (MHP) halo model. Although the time delay for both NIE and MHP halo models is well constrained on our parameter grid, it strongly depends on the halo surface density profile. We furthermore find that the presence of a flat luminous mass distribution can severely alter the statistical properties of the lens.     

Period changes in the pulsating extreme helium stars V652 Her and BX Cir

A new statement of the eigenvalue problem of studying small perturbations in arbitrary integrable self-gravitating systems is presented. An example of such a system, a 2D stellar disc, is considered in detail. The theory, based on the general equation for disc eigenmodes, reveals mechanisms for the formation and growth of global galactic structures. This new point of view specifies the limits of the unified theory of bar-like and spiral modes that was based on the assumption that global galactic structures could be understood in terms of low-frequency disc modes.     

Dynamics of an interacting luminous disc, dark halo and satellite companion

This paper describes a method, based on linear perturbation theory, to determine the dynamical interaction between extended halo and spheroid components and an environmental disturbance. One finds that resonant interaction between a galaxy and passing interlopers or satellite companions can carry the disturbance inward, deep inside the halo, where it can perturb the disc.    Applied to the Milky Way for example, the LMC and SMC appear to be sufficient to cause the observed Galactic warp and possibly seed other asymmetries. This is a multi-scale interaction in which the halo wake has a feature at roughly half the satellite orbital radius owing to a 2:1 orbital resonance. The rotating disturbance then excites an m  = 1 vertical disc mode which has the classic integral-sign morphology. A polar satellite orbit produces the largest warp and therefore the inferred LMC orbit is nearly optimal for maximum warp production.   Both the magnitude and morphology of the response depend on the details of the disc and halo models. Most critically, a change in the halo profile will shift the resonant frequencies and response location and consequently alter the coupling to the bending disc. Increasing the halo support relative to the disc, a sub-maximal disc model, decreases the warp amplitude.   Finally, the results and prognosis for N -body simulations are discussed. Discreteness noise in the halo, similar to that arising from a population of 10⁶-M⊙ black holes, can produce observable warping.     

Time evolution of galactic warps in prolate haloes

A recent observation with the Hipparcos satellite and some numerical simulations imply that the interaction between an oblate halo and a disc is inappropriate for the persistence of galactic warps. Following on from this , we have compared the time evolution of galactic warps in a prolate halo with that in an oblate halo. The haloes were approximated as fixed potentials, while the discs were represented by N -body particles. We have found that the warping in the oblate halo continues to wind up, and finally disappears. On the other hand, for the prolate halo model, the precession rate of the outer disc increases when the precession of the outer disc recedes from that of the inner disc, and vice versa. Consequently, the warping in the prolate halo persisted to the end of the simulation by retaining the alignment of the line of nodes of the warped disc. Therefore, our results suggest that prolate haloes could sustain galactic warps. The physical mechanism of the persistence of warp is discussed on the basis of the torque between a halo and a disc and that between the inner and outer regions of the disc.     

Bar-induced evolution of dark matter cusps

The evolution of a stellar bar transforms not only the galactic disc, but also the host dark matter halo. We present high-resolution, fully self-consistent N -body simulations that clearly demonstrate that dark matter halo central density cusps flatten as the bar torques the halo. This effect is independent of the bar formation mode and occurs even for rather short bars. The halo and bar evolution is mediated by resonant interactions between orbits in the halo and the bar pattern speed, as predicted by linear Hamiltonian perturbation theory. The bar lengthens and slows as it loses angular momentum, a process that occurs even in rather warm discs. We demonstrate that the bar and halo response can be critically underestimated for experiments that are unable to resolve the relevant resonant dynamics; this occurs when the phase space in the resonant region is undersampled or plagued by noise.     

Contrasting the chemical evolution of the Milky Way and Andromeda

The chemical evolution history of a galaxy hides clues about how it formed and has been changing through time. We have studied the chemical evolution history of the Milky Way (MW) and Andromeda (M31) to find which are common features in the chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a semi-analytic multizone chemical evolution model. Such models have succeeded in explaining the mean trends of the observed chemical properties in these two Local Group spiral galaxies with similar mass and morphology. Our results suggest that while the evolution of the MW and M31 shares general similarities, differences in the formation history are required to explain the observations in detail. In particular, we found that the observed higher metallicity in the M31 halo can be explained by either (i) a higher halo star formation efficiency (SFE), or (ii) a larger reservoir of infalling halo gas with a longer halo formation phase. These two different pictures would lead to (i) a higher [O/Fe] at low metallicities, or (ii) younger stellar populations in the M31 halo, respectively. Both pictures result in a more massive stellar halo in M31, which suggests a possible correlation between the halo metallicity and its stellar mass.     

Self-consistent response of a galactic disc to vertical perturbations

We study the self-consistent, linear response of a galactic disc to vertical perturbations, as induced, say, by a tidal interaction. We calculate the self-gravitational potential corresponding to a non-axisymmetric, self-consistent density response of the disc using the Green's function method. The response potential is shown to oppose the perturbation potential because the self-gravity of the disc resists the imposed potential, and this resistance is stronger in the inner parts of a galactic disc. For the   m = 1  azimuthal wavenumber, the disc response opposes the imposed perturbation up to a radius that spans a range of 4–6 disc scalelengths, so that the disc shows a net warp only beyond this region. This physically explains the well known but so far unexplained observation that warps typically set in beyond this range of radii. We show that the inclusion of a dark matter halo in the calculation only marginally changes (by ∼10 per cent) the radius for the onset of warps. For perturbations with higher azimuthal wavenumbers, the net signature of the vertical perturbations can only be seen at larger radii – for example, beyond 7 exponential disc scalelengths for   m = 10  . Also, for the high- m cases, the magnitude of the negative disc response due to the disc self-gravity is much smaller. This is shown to result in corrugations of the mid-plane density, which explains the puzzling scalloping with   m = 10  detected in H  i in the outermost regions ∼30 kpc in the Galaxy.     

Dark matter halo response to the disc growth

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter haloes in the Λ cold dark matter (ΛCDM) paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though haloes and discs may assemble simultaneously. We further demonstrate the validity of the simple approximation for ΛCDM haloes with isotropic velocity distributions using three-dimensional N -body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response of a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low angular momentum orbits in the vicinity of the disc and the growth time-scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modelling the response of dark matter haloes to the growth of a disc.     

Tidal tails in cold dark matter cosmologies

We study the formation of tidal tails in pairs of merging disc galaxies with structural properties motivated by current theories of cold dark matter (CDM) cosmologies. In a recent study, Dubinski, Mihos & Hernquist showed that the formation of prominent tidal tails can be strongly suppressed by massive and extended dark haloes. For the large halo-to-disc mass ratio expected in CDM cosmologies their sequence of models failed to produce strong tails like those observed in many well-known pairs of interacting galaxies. In order to test whether this effect can constrain the viability of CDM cosmologies, we construct N ‐body models of disc galaxies with structural properties derived in analogy to the recent analytical work of Mo, Mao & White. With a series of self-consistent collisionless simulations of galaxy–galaxy mergers we demonstrate that even the discs of very massive dark haloes have no problems developing long tidal tails, provided the halo spin parameter is large enough. For our class of models, the halo-to-disc mass ratio is not a good indicator of the ability to produce tails. Instead, the relative size of disc and halo or, alternatively, the ratio of circular velocity to local escape speed at the half mass radius of the disc is a more useful criterion. This result holds in all CDM models. While tidal tails can provide useful information on the structure of galaxies, it thus appears unlikely that they are able to constrain the values of the cosmological parameters within these models.     

The kinematics of lopsided galaxies

Lopsidedness is a common feature in galaxies, both in the distribution of light and in the kinematics. We investigate the kinematics of a model for lopsided galaxies that consists of a disc lying off-centre in a dark halo, and circling around the halo centre. We search for families of stable, closed, non-crossing orbits, and assume that gas in our galaxies moves on these orbits. Several of our models show strong lopsided gas kinematics, especially those in which the disc spins around its axis in a retrograde sense compared with its motion around the halo centre. We are able to reproduce the H  i velocity map of the kinematically lopsided galaxy NGC 4395.
The lopsidedness in our models is most pronounced in the models where the halo provides a relatively large fraction of the total mass at small radii. This may explain why the gas shows lopsidedness more frequently in late-type galaxies, which are dominated by dark matter. Surfaces of section show large regions of irregular orbits in the models where the halo density is low. This may indicate that these models are unstable.     

The environment of formation as a second parameter for globular cluster classification

We perform an evolutionary multivariate analysis of a sample of 54 Galactic globular clusters with high-quality colour–magnitude diagrams and well-determined ages. The four parameters adopted for the analysis are: metallicity, age, maximum temperature on the horizontal branch and absolute V magnitude. Our cladistic analysis breaks the sample into three novel groups. An a posteriori kinematical analysis puts groups 1 and 2 in the halo, and group 3 in the thick disc. The halo and disc clusters separately follow a luminosity–metallicity relation of much weaker slope than galaxies. This property is used to propose a new criterion for distinguishing halo and disc clusters. A comparison of the distinct properties of the two halo groups with those of Galactic halo field stars indicates that the clusters of group 1 originated in the inner halo, while those of group 2 formed in the outer halo of the Galaxy. The inner halo clusters were presumably initially the most massive one, which allowed the formation of more strongly helium-enriched second generation stars, thus explaining the presence of Cepheids and of very hot horizontal-branch stars exclusively in this group. We thus conclude that the 'second parameter' is linked to the environment in which globular clusters form, the inner halo favouring the formation of the most massive clusters which subsequently become more strongly self-enriched than their counterparts of the galactic outer halo and disc.     

Confirmation of a faint red halo around NGC 5907

We report the detection of extended infrared emission 5.2 kpc above the plane of the edge-on spiral galaxy NGC 5907 in the J K bands. This provides confirmation at two other wavelengths of the R -band result of Sackett et al. The halo J − K colour is found to be 1.3 ± 0.3. If we make the assumption that the halo emission is of the form necessary to give a flat rotation curve at large radii, with a power-law index in projection of −1.0, the K surface brightness, M/L_K ratio and R − K colour are found to be 23.0 mag arcsec⁻², 100 and 3.5 respectively, at 95-arcsec radius. Adopting a halo index of −2.5 in projection, typical of known stellar haloes but inconsistent with the halo measured by Sackett et al., the K surface brightness at 95 arcsec drops to 24.0 mag arcsec⁻² and the J − K colour to 2.5, while the M/L_K ratio is undetermined since such a halo would be unrelated to the dark matter component. If the emission is interpreted as arising from the massive halo, then it is broadly consistent with a population of stars at the very lowest masses, possibly extending into the brown dwarf regime.     

Do high-velocity clouds form by thermal instability?

We examine the proposal that the H  i 'high-velocity' clouds (HVCs) surrounding the Milky Way and other disc galaxies form by condensation of the hot galactic corona via thermal instability. Under the assumption that the galactic corona is well represented by a non-rotating, stratified atmosphere, we find that for this formation mechanism to work the corona must have an almost perfectly flat entropy profile. In all other cases, the growth of thermal perturbations is suppressed by a combination of buoyancy and thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes smaller than  10 kpc  could form in the corona of the Milky Way only at radii larger than  100 kpc  , in contradiction with the determined distances of the largest HVC complexes. Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. We conclude that unless even slow rotation qualitatively changes the dynamics of a corona, thermal instability is unlikely to be a viable mechanism for formation of cold clouds around disc galaxies.     

Chemo-spectrophotometric evolution of spiral galaxies – II. Main properties of present-day disc galaxies

We study the chemical and spectrophotometric evolution of galactic discs with detailed models calibrated on the Milky Way and using simple scaling relations, based on currently popular semi-analytic models of galaxy formation. We compare our results with a large body of observational data on present-day galactic discs, including disc sizes and central surface brightness, Tully–Fisher relations in various wavelength bands, colour–colour and colour–magnitude relations, gas fractions versus magnitudes and colours and abundances versus local and integrated properties, as well as spectra for different galactic rotational velocities. Despite the extremely simple nature of our models, we find satisfactory agreement with all those observables, provided that the time-scale for star formation in low-mass discs is longer than for more massive ones. This assumption is apparently in contradiction with the standard picture of hierarchical cosmology. We find, however, that it is extremely successful in reproducing major features of present-day discs, like the change in the slope of the Tully–Fisher relation with wavelength, the fact that more massive galaxies are on average 'redder' than low-mass ones (a generic problem of standard hierarchical models) and the metallicity–luminosity relation for spirals. It is concluded that, on a purely empirical basis, this new picture is at least as successful as the standard one. Observations at high redshifts could help to distinguish between the two possibilities.