Boxy/peanut 'bulges': comparing the structure of galaxies with the underlying families of periodic orbits

The vertical profiles of disc galaxies are built by the material trapped around stable periodic orbits, which form their 'skeletons'. Therefore, knowledge of the stability of the main families of periodic orbits in appropriate 3D models enables one to predict possible morphologies for edge-on disc galaxies. In a pilot survey we compare the orbital structures that lead to the appearance of 'peanut'- and 'X'-like features with the edge-on profiles of three disc galaxies (IC 2531, NGC 4013 and UGC 2048). The subtraction from the images of a model representing the axisymmetric component of the galaxies reveals the contribution of the non-axisymmetric terms. We find a direct correspondence between the orbital profiles of 3D bars in models and the observed main morphological features of the residuals. We also apply a simple unsharp masking technique in order to study the sharpest features of the images. Our basic conclusion is that the morphology of the boxy 'bulges' of these galaxies can be explained by considering disc material trapped around stable 3D periodic orbits. In most models, these building-block periodic orbits are bifurcated from the planar central family of a non-axisymmetric component, usually a bar, at low-order vertical resonances. In such a case, the boxy 'bulges' are parts of bars seen edge-on. For the three galaxies we study, the families associated with the 'peanut' or 'X'-shape morphology are probably bifurcations at the vertical 2/1 or 4/1 resonance.     

Galactic bulges from Hubble Space Telescope NICMOS observations: ages and dust

We present optical and near-infrared colour maps of the central regions of bulges of S0 and spiral galaxies obtained with WFPC2 and NICMOS on the Hubble Space Telescope ( HST ). By combined use of HST and ground-based data, the colour information spans a region from a few tens of pc to a few kpc. In almost all galaxies, the colour profiles in the central 100–200 pc become more rapidly redder. We attribute the high central colour indices to a central concentration of dust. We infer an average extinction at the centre of A _V =0.6–1.0 mag. Several objects show central dust rings or discs at subkpc scales similar to those found by others in giant ellipticals. For galactic bulges of types S0 to Sb, the tightness of the B − I versus I − H relation suggests that the age spread among bulges of early-type galaxies is small, at most 2 Gyr. Colours at 1 R _eff, where we expect extinction to be negligible, are similar to those of elliptical galaxies in the Coma cluster, suggesting that these bulges formed at the same time as the bright galaxies in Coma. Furthermore, the galaxy ages are found to be independent of their environment. As it is likely that Coma was formed at redshift z >3, our bulges, which are in groups and in the field, must also have been formed at this epoch. Bulges of early-type spirals cannot be formed by secular evolution of bars at recent epochs, because such bulges would be much younger. There are three galaxies of type Sbc and later; their bulges are younger and could perhaps arise from secular evolution of transient bars. Our results are in good agreement with semi-analytic predictions, which also predict that bulges, in clusters and in the field, are as old as giant ellipticals in clusters.     

A SAURON look at galaxy bulges

Kinematic and population studies show that bulges are generally rotationally flattened systems similar to lowluminosity ellipticals. However, observations with state‐of‐the‐art integral field spectrographs, such as SAURON, indicate that the situation is much more complex, and allow us to investigate phenomena such as triaxiality, kinematic decoupling and population substructure, and to study their connection to current formation and evolution scenarios for bulges of early‐type galaxies. We present the examples of two S0 bulges from galaxies in our sample of nearby galaxies: one that shows all the properties expected from classical bulges (NGC5866), and another case that presents kinematic features appropriate for barred disk galaxies (NGC7332). (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Unravelling the mystery of the M31 bar

The inclination of M31 is too close to edge-on for a bar component to be easily recognized and is not sufficiently edge-on for a boxy/peanut bulge to protrude clearly out of the equatorial plane. Nevertheless, a sufficient number of clues allow us to argue that this galaxy is barred. We use fully self-consistent N -body simulations of barred galaxies and compare them with both photometric and kinematic observational data for M31. In particular, we rely on the near-infrared photometry presented in a companion paper. We compare isodensity contours to isophotal contours and the light profile along cuts parallel to the galaxy major axis and offset towards the north, or the south, to mass profiles along similar cuts on the model. All these comparisons, as well as position–velocity diagrams for the gaseous component, give us strong arguments that M31 is barred. We compare four fiducial N -body models to the data and thus set constraints on the parameters of the M31 bar, as its strength, length and orientation. Our 'best' models, although not meant to be exact models of M31, reproduce in a very satisfactory way the main relevant observations. We present arguments that M31 has both a classical and a boxy/peanut bulge. Its pseudo-ring-like structure at roughly 50 arcmin is near the outer Lindblad resonance of the bar and could thus be an outer ring, as often observed in barred galaxies. The shape of the isophotes also argues that the vertically thin part of the M31 bar extends considerably further out than its boxy bulge, that is, that the boxy bulge is only part of the bar, thus confirming predictions from orbital structure studies and from previous N -body simulations. It seems very likely that the backbone of M31's boxy bulge is families of periodic orbits, members of the x₁-tree and bifurcating from the x₁ family at its higher order vertical resonances, such as the x1v3 or x1v4 families.     

Dynamical Evolution of Barred Galaxies

Angular momentum redistribution within barred galaxies drives their dynamical evolution. Angular momentum is emitted mainly by near-resonant material in the bar region and absorbed by resonant material mainly in the outer disc and in the halo. This exchange determines the strength of the bar, the decrease of its pattern speed, as well as its morphology. If the galaxy has also a gaseous component and/or a companion or satellite, then these also take part in the angular momentum exchange. During the evolution a bar structure forms in the inner parts of the halo as well. This bar is shorter and fatter than the disc bar and stays so all through the simulation, although its length grows considerably with time. Viewed edge-on, the bar in the disc component acquires a boxy or peanut shape. I describe the families of periodic orbits that explain such structures and review the observations showing that boxy/peanut ‘bulges’ are in fact just bars seen edge-on.     

Constraining the star formation histories of spiral bulges

Stellar populations in spiral bulges are investigated using the Lick system of spectral indices. Long-slit spectroscopic observations of line strengths and kinematics made along the minor axes of four spiral bulges are reported. Comparisons are made between central line strengths in spiral bulges and those in other morphological types [elliptical, spheroidal (Sph) and S0]. The bulges investigated are found to have central line strengths comparable to those of single stellar populations of approximately solar abundance or above. Negative radial gradients are observed in line strengths, similar to those exhibited by elliptical galaxies. The bulge data are also consistent with correlations between Mg₂, Mg₂ gradient and central velocity dispersion observed in elliptical galaxies. In contrast to elliptical galaxies, central line strengths lie within the loci defining the range of 〈Fe〉 and Mg₂ achieved by Worthey's solar abundance ratio, single stellar populations (SSPs). The implication of solar abundance ratios indicates significant differences in the star formation histories of spiral bulges and elliptical galaxies. A 'single zone with infall' model of galactic chemical evolution, using Worthey's SSPs, is used to constrain the possible star formation histories of our sample. We show that the 〈Fe〉, Mg₂ and H β line strengths observed in these bulges cannot be reproduced using primordial collapse models of formation but can be reproduced by models with extended infall of gas and star formation (2–17 Gyr) in the region modelled. One galaxy (NGC 5689) shows a central population with a luminosity-weighted average age of ∼5 Gyr, supporting the idea of extended star formation. Kinematic substructure, possibly associated with a central spike in metallicity, is observed at the centre of the Sa galaxy NGC 3623.     

Structural properties of pseudo-bulges, classical bulges and elliptical galaxies: a Sloan Digital Sky Survey perspective

We have performed 2D bulge/bar/disc decompositions using g , r and i -band images of a representative sample of nearly 1000 galaxies from the Sloan Digital Sky Survey. We show that the Petrosian concentration index is a better proxy for the bulge-to-total ratio than the global Sérsic index. We show that pseudo-bulges can be distinguished from classical bulges as outliers in the Kormendy relation. We provide the structural parameters and distributions of stellar masses of ellipticals, classical bulges, pseudo-bulges, discs and bars, and find that 32 per cent of the total stellar mass in massive galaxies in the local universe is contained in ellipticals, 36 per cent in discs, 25 per cent in classical bulges, 3 per cent in pseudo-bulges and 4 per cent in bars. Pseudo-bulges are currently undergoing intense star formation activity and populate the blue cloud of the colour–magnitude diagram. Most (though not all) classical bulges are quiescent and populate the red sequence of the diagram. Classical bulges follow a correlation between the bulge Sérsic index and bulge-to-total ratio, while pseudo-bulges do not. In addition, for a fixed bulge-to-total ratio, pseudo-bulges are less concentrated than classical bulges. Pseudo-bulges follow a mass–size relation similar to that followed by bars, and different from that followed by classical bulges. In the fundamental plane, pseudo-bulges occupy the same locus as discs. While these results point out different formation processes for classical and pseudo-bulges, we also find a significant overlap in their properties, indicating that the different processes might happen concomitantly. Finally, classical bulges and ellipticals follow offset mass–size relations, suggesting that high-mass bulges might not be simply high-mass ellipticals surrounded by discs.     

Stellar populations of the bulges of four spiral galaxies

Key information to understand the formation and evolution of disk galaxies are imprinted in the stellar populations of their bulges. This paper has the purpose to make available new measurements of the stellar population properties of the bulges of four spiral galaxies. Both the central values and radial profiles of the line strength of some of the most common Lick indices are measured along the major‐ and minor‐ axis of the bulge‐dominated region of the sample galaxies. The corresponding age, metallicity, and α /Fe ratio are derived by using the simple stellar population synthesis model predictions. The central values and the gradients of the stellar population properties of ESO‐LV 1890070, ESO‐LV 4460170, and ESO‐LV 5140100 are consistent with previous findings for bulges of spiral galaxies. On the contrary, the bulge of ESO‐LV 4500200 shows peculiar chemical properties possibly due to the presence of a central kinematically‐decoupled component. The negative metallicity gradient found in our bulges sample indicates a relevant role for the dissipative collapse in bulge formation. However, the shallow gradients found for the age and α /Fe ratio suggests that merging can not be completely ruled out for the sample bulges. This is confirmed by the properties of ESO‐LV 4500200 which can hardly be explained without invoking the capture of external material. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase-space density constraints on the formation of bulges from discs

The colours of stellar bulges and of inner stellar discs are comparable, and consistent with rather similar mean metallicities and ages. Indeed, the mean chemical abundances of the Milky Way bulge and old disc are approximately equal. Further, the scalelengths of discs and bulges are correlated. These observations imply a close relationship between discs and bulges, and may support models in which stellar bulges form from stellar discs. The present paper discusses constraints on this scenario from the stellar phase-space density of bulges and of discs. Phase-space density cannot increase in the absence of collisional processes. We show here that the maximum phase-space density of galactic bulges is higher than that of inner discs, arguing that instabilities of purely stellar discs cannot form bulges. Rather, the high densities of bulges probably reflect gaseous dissipation. Gas inflow from the disc would complicate the interpretation of the similarities in stellar colours between discs and bulges. Gas inflow from the stellar halo, if one exists, may be favoured on angular momentum grounds, but this means of formation of the bulge would provide no explanation for the relationships between disc and bulge in any one galaxy. At least in the Milky Way, the metallicity distribution of the bulge is not consistent with the bulge being built up from the dense regions of accreted satellite galaxies and/or globular clusters.     

Galactic bulges

We model the evolution of the Galactic bulge and of the bulges of a selected sample of external spiral galaxies, via the multiphase multizone evolution model. We address a few questions concerning the role of the bulges within galactic evolution schemes and the properties of bulge stellar populations. We provide solutions to the problems of chemical abundances and spectral indices, the two main observational constraints to bulge structure.     

The black hole mass–stellar velocity dispersion correlation: bulges versus pseudo-bulges

We investigate the correlation between the supermassive black holes (SMBHs) mass ( M _bh) and the stellar velocity dispersion  (σ_*)  in two types of host galaxies: the early-type bulges (disc galaxies with classical bulges or elliptical galaxies) and pseudo-bulges. In the form  log ( M _bh/M_⊙) =α+β log (σ_*/200 km s⁻¹)  , the best-fitting results for the 39 early-type bulges are the slope  β= 4.06 ± 0.28  and the normalization  α= 8.28 ± 0.05  ; the best-fitting results for the nine pseudo-bulges are  β= 4.5 ± 1.3  and  α= 7.50 ± 0.18  . Both relations have intrinsic scatter in  log  M _bh  of ≲0.27 dex. The   M _bh–σ_*  relation for pseudo-bulges is different from the relation in the early-type bulges over the 3σ significance level. The contrasting relations indicate the formation and growth histories of SMBHs depend on their host type. The discrepancy between the slope of the   M _bh–σ_*  relations using different definition of velocity dispersion vanishes in our sample, a uniform slope will constrain the coevolution theories of the SMBHs and their host galaxies more effectively. We also find the slope for the 'core' elliptical galaxies at the high-mass range of the relation appears steeper  (β≃ 5–6)  , which may be the imprint of their origin of dissipationless mergers.     

The growth of supermassive black holes in pseudo-bulges, classical bulges and elliptical galaxies

Using results from structural analysis of a sample of nearly 1000 local galaxies from the Sloan Digital Sky Survey, we estimate how the mass in central black holes is distributed amongst elliptical galaxies, classical bulges and pseudo-bulges, and investigate the relation between their stellar masses and central stellar velocity dispersion σ. Assuming a single relation between elliptical galaxy/bulge mass, M _Bulge, and central black hole mass, M _BH, we find that  55⁺⁸₋₄  per cent of the mass in black holes in the local universe is in the centres of elliptical galaxies,  41⁺⁴₋₂  per cent in classical bulges and  4^+0.9_−0.4  per cent in pseudo-bulges. We find that ellipticals, classical bulges and pseudo-bulges follow different relations between their stellar masses and σ, and the most significant offset occurs for pseudo-bulges in barred galaxies. This structural dissimilarity leads to discrepant black hole masses if single   M _BH– M _Bulge  and   M _BH–σ  relations are used. Adopting relations from the literature, we find that the   M _BH–σ  relation yields an estimate of the total mass density in black holes that is roughly 55 per cent larger than if the   M _BH– M _Bulge  relation is used.     

Elliptical galaxies from mergers of discs

We analyse N -body galaxy merger experiments involving disc galaxies. Mergers of disc–bulge–halo models are compared to those of bulgeless, disc–halo models to quantify the effects of the central bulge on merger dynamics and the structure of the remnant. Our models explore galaxy mass ratios 1:1 through 3:1, and use higher bulge mass fractions than previous studies. A full comparison of the structural and dynamical properties with our observations is carried out. The presence of central bulges results in longer tidal tails, oblate final intrinsic shapes, surface brightness profiles with a higher Sérsic index, steeper rotation curves and oblate-rotator internal dynamics. Mergers of bulgeless galaxies do not generate long-lasting tidal tails, and their strong triaxiality seems inconsistent with observations; these remnants show shells, which we do not find in models including central bulges. Giant ellipticals with boxy isophotes and anisotropic dynamics cannot be produced by the mergers modelled here; they could be the result of mergers between lower luminosity ellipticals, themselves plausibly formed in disc-disc mergers.     

Morphology, photometry and kinematics of N-body bars – I. Three models with different halo central concentrations

We discuss the morphology, photometry and kinematics of the bars which have formed in three N -body simulations. These have initially the same disc and the same halo-to-disc mass ratio, but their haloes have very different central concentrations. The third model includes a bulge. The bar in the model with the centrally concentrated halo (model MH) is much stronger, longer and thinner than the bar in the model with the less centrally concentrated halo (model MD). Its shape, when viewed side-on, evolves from boxy to peanut and then to 'X'-shaped, as opposed to that of model MD, which stays boxy. The projected density profiles obtained from cuts along the bar major axis, for both the face-on and the edge-on views, show a flat part, as opposed to those of model MD which are falling rapidly. A Fourier analysis of the face-on density distribution of model MH shows very large  m=2  , 4, 6 and 8 components. Contrary to this, for model MD the components  m=6  and 8 are negligible. The velocity field of model MH shows strong deviations from axial symmetry, and in particular has wavy isovelocities near the end of the bar when viewed along the bar minor axis. When viewed edge-on, it shows cylindrical rotation, which the MD model does not. The properties of the bar of the model with a bulge and a non-centrally concentrated halo (MDB) are intermediate between those of the bars of the other two models. All three models exhibit a lot of inflow of the disc material during their evolution, so that by the end of the simulations the disc dominates over the halo in the inner parts, even for model MH, for which the halo and disc contributions were initially comparable in that region.     

Morphology and kinematics in thick box‐peanut bulge galaxies

Both, observational and theoretical results support the close connection between bars and box/peanut‐shaped bulges in spiral galaxies. However, for the sub‐class of thick box/peanut‐shaped bulges, an environmental process instead of an internal mechanism is suggested to account for the prominence of the box/peanut component. In this paper a multiwavelength study of the morphology and the kinematics of a sample of thick box/peanut‐shaped bulges is presented. The results suggest a scenario in which interaction events play a significant role in the history of thick box/peanut‐shaped bulge galaxies. Most likely these sources have experienced a mass infall. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stellar populations of classical and pseudo-bulges for a sample of isolated spiral galaxies

In this paper we present the stellar population synthesis results for a sample of 75 bulges in isolated spiral Sb-Sc galaxies, using the spectroscopic data from the Sloan Digital Sky Survey and the STARLIGHT code. We find that both pseudo-bulges and classical bulges in our sample are predominantly composed of old stellar populations, with mean mass-weighted stellar age around 10 Gyr. While the stellar population of pseudo-bulges is, in general, younger than that of classical bulges, the difference is not significant, which indicates that it is hard to distinguish pseudo-bulges from classical bulges, at least for these isolated galaxies, only based on their stellar populations. Pseudo-bulges have star formation activities with relatively longer timescale than classical bulges, indicating that secular evolution is more important in this kind of systems. Our results also show that pseudo-bulges have a lower stellar velocity dispersion than their classical counterparts, which suggests that classical bulges are more dispersion-supported than pseudo-bulges.     

A comparison of stellar populations in galaxy spheroids across a wide range of Hubble types

We present line-strengths and kinematics from the central regions of 32 galaxies with Hubble types ranging from E to Sbc. Spectral indices, based on the Lick system, are measured in the optical and near-infrared (NIR). The 24 indices measured, in conjunction with models of the effects of varying abundance ratios, permit the breaking of age/metallicity degeneracy, and allow estimation of enhancements in specific light elements (particularly C and Mg). The large range of Hubble types observed allows direct comparison of line-strengths in the centres of early-type galaxies (E and S0) with those in spiral bulges, free from systematic differences that have plagued comparisons of results from different studies. Our sample includes field and Virgo cluster galaxies. For early-type galaxies our data are consistent with previously reported trends of Mg₂ and Mgb with velocity dispersion. In spiral bulges we find trends in all indices with velocity dispersion. We estimate luminosity-weighted ages, metallicities and heavy-element abundance ratios (enhancements) from optical indices. These show that bulges are less enhanced in light ( α -capture) elements and have lower average age than early-type galaxies. Trends involving age and metallicity also differ sharply between early and late types. An anticorrelation exists between age and metallicity in early types, while, in bulges, metallicity is correlated with velocity dispersion. We consider the implications of these findings for models of the formation of these galaxies. We find that primordial collapse models of galaxy formation are ruled out by our observations, while several predictions of hierarchical clustering (merger) models are confirmed.     

The global structure of galactic discs

A statistical study of global galaxy parameters can help to improve our understanding of galaxy formation processes. In this paper we present the analysis of global galaxy parameters based on optical and near-infrared observations of a large sample of edge-on disc galaxies. We found a correlation between the ratio of the radial to vertical scale parameter and galaxy type: galaxies become systematically thinner when going from type S0 to type Sc, whereas the distribution seems to level off for later types. The observed scalelength ratios (and thus the radial colour gradients) largely represent the dust content of the galaxies. On average, the colour gradients indicated by the scalelength ratios increase from type Sa to at least type Sc. For galaxy types later than Sc, the average colour gradient seems to decrease again. The distribution of K -band (edge-on) disc central surface brightnesses is rather flat, although with a large scatter. However, the latest-type sample galaxies ( T  > 6) show an indication that their average disc central surface brightnesses may be fainter than those of the earlier types. This effect is probably not the result of dust extinction.     

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.     

Near-infrared surface photometry of bulges and disks of spiral galaxies. The data

We present optical and near-infrared (NIR) surface brightness and colour profiles, in bands ranging from U to K, for the disk and bulge components of a complete sample of 30 nearby S0 to Sbc galaxies with inclinations larger than 50 °. We describe in detail the observations and the determination of colour parameters. Calibrated monochromatic and real-colour images are presented, as well as colour index maps. This data set, tailored for the study of the population characteristics of galaxy bulges, provides useful information on the colours of inner disks as well. In related papers, we have used them to quantify colour gradients in bulges, and age differentials between bulge and inner disk.