Resolved Stellar Populations in the Bulge and the Magellanic Clouds ASP Conference Series, Vol. eds. S. Points, A. Kunder c © Astronomical Society of the Pacific

## Formation Models of the Galactic Bulge

Ortwin Gerhard 1

1 Max Planck Institute for extraterrestrial Physics, PO Box 1312, Giessenbachstr., 85741 Garching, Germany

Abstract. The Galactic bulge is now considered to be the inner three-dimensional part of the Milky Way's bar. It has a peanut shape and is characterized by cylindrical rotation. In N-body simulations, box peanut bulges arise from disks through bar and / buckling instabilities. Models of this kind explain much of the structure and kinematics of the Galactic bulge and, in principle, also its vertical metallicity gradient. Cosmological disk galaxy formation models with high resolution and improved feedback models are now able to generate late-type disk galaxies with disk-like or barred bulges. These bulges often contain an early collapse stellar population and a population driven by later disk instabilities. Due to the inside-out disk formation, these bulges can be predominantly old, similar to the Milky Way bulge.

## 1. Introduction: the barred Galactic bulge

The Galactic bulge, subject of many studies here at CTIO observatory, is now considered to be the inner three-dimensional part of the Milky Way's bar. Its triaxial shape was first established by the COBE NIR photometry (Blitz &amp; Spergel 1991; Dwek et al. 1995; Binney et al. 1997; Bissantz &amp; Gerhard 2002). Star counts confirmed this but also pointed to a larger, in-plane bar, the so-called 'long bar' (Stanek et al. 1994; L´pezo Corredoira et al. 2005; Skrutskie et al. 2006; Benjamin et al. 2005; Cabrera-Lavers et al. 2007). Comparing the observed HI and CO lv-diagrams with hydrodynamic models showed that many of the observed features could be naturally, sometimes even quantitatively interpreted with gas flow models in barred potentials, modeled on the COBE data (Englmaier &amp; Gerhard 1999; Fux 1999; Bissantz et al. 2003; Rodriguez-Fernandez &amp;Combes 2008).

In recent years, Galactic bulge studies have been revitalized by new large photometric and spectroscopic surveys. With these new data, the density and metallicity structure of the Milky Way's bulge and its kinematics and dynamics have been mapped in much greater detail than possible before. Comparison of this information with models of bulge formation is promising to give us a better understanding of the origin of the Galactic bulge.

## 2. Box peanut bulge models from disk instability /

N-body simulations following the evolution of isolated disk galaxies have shown that stellar disks are often unstable to bar formation (Sellwood 1981; Athanassoula 2002; Debattista et al. 2006). Subsequently, the bar may quickly go through a second, so-

called buckling instability, resulting in the formation of an inner boxy bulge (Combes &amp; Sanders 1981; Combes et al. 1990; Raha et al. 1991). The bar may then evolve by losing angular momentum to the dark matter halo, grow in size, and eventually go through a second buckling instability. This results in a a strongly peanut-shaped bulge (MartinezValpuesta et al. 2006) much like those seen in barred galaxies (Athanassoula 2005), which is supported by three-dimensional 2:1:2 resonant orbit families that support a characteristic X-shape (Pfenniger &amp; Friedli 1991).

Recent star count and stellar kinematics data have shown that the predictions of these idealized models characterize the properties of the Milky Way bulge surprisingly well. The apparent magnitude plot for red clump stars in high-latitude pencil beams shows a clearly bimodal distance distribution, termed the 'split red clump'. This was first seen in 2MASS data (McWilliam &amp; Zoccali 2010; Saito et al. 2011) but has also been confirmed with OGLE data (Nataf et al. 2010) and the ARGOS spectroscopic bulge sample (Ness et al. 2012). The split red clump is a signature of the X-shaped stellar orbits in the bulge. Recent deconvolution of the VVV DR1 star count data has shown that the Galactic bulge has the shape of a highly elongated bar with a strong peanut shape (Wegg &amp; Gerhard 2013). The measured angle of the barred bulge to the line of sight is (27 ± 2) ◦ . Along the bar axes the density falls o ff roughly exponentially, with axis ratios (10 : 6.3 : 2.6) and exponential scale-lengths (0.70 : 0.44 : 0.18) kpc, but along the major axis the profile becomes shallower beyond 1.5 kpc.

N-body models with box peanut bulges grown from unstable disks have been used / to interpret not only the structure, but also the stellar kinematics and metallicity distribution of the Milky Way bulge. Gerhard &amp; Martinez-Valpuesta (2012) showed that the flattening in the longitude profiles seen for | l | &lt; 4 ◦ in low-latitude star counts (Nishiyama et al. 2005; Gonzalez et al. 2011) is naturally reproduced by the density distribution in a simulated box peanut bulge. / Li &amp; Shen (2012) showed that their Nbody model shows a similar split red clump as the Galactic bulge when viewed at an appropriate orientation. N-body bulges also explain the near-cylindrical rotation and the velocity dispersion profiles seen in the BRAVA (Kunder et al. 2012) and Argos (Ness et al. 2013a) spectroscopic surveys (Shen et al. 2010; Ness et al. 2013b).

The match of the models to the kinematic data is very good and Shen et al. (2010) concluded that it would be worsened even by a fairly small additional classical bulge. Such a classical bulge could have arisen from an early merger episode that preceded the growth and subsequent instability of the early disk. Saha et al. (2012); Saha &amp; Gerhard (2013) showed that angular momentum transfer from the bar and boxy bulge can spin up the preexisting bulge, so that it too develops cylindrical rotation. The radial dependence of the rotation is di ff erent between the two bulges though; thus some more work is needed on this issue.

The strongest argument for a classical bulge in the Milky Way has been the observation of a strong vertical metallicity gradient (Zoccali et al. 2008; Johnson et al. 2011; Gonzalez et al. 2013). However, a recent analysis of an N-body bulge showed that this argument is not compelling. If the Milky Way's bar and bulge formed rapidly from the disk at early times, then any preexisting metallicity gradients are preserved as vertical gradients in the final bulge, because of the incomplete violent relaxation during the instabilities (Martinez-Valpuesta &amp; Gerhard 2013). Therefore, at present there is no strong evidence for a classical bulge in the Milky Way. Further chemodynamical analysis of the metallicity distribution and kinematics will be needed to clarify whether the Milky Way has a small classical bulge underneath its dominant peanut-shaped bulge,

and if so how this is related to the bulge metallicity components identified by Ness et al. (2013a) in the ARGOS sample.

The connection between the Milky Way's peanut bulge and planar bar is much less explored, primarily because this region is observationally less accessible. Hammersley et al. (2000); Benjamin et al. (2005); Cabrera-Lavers et al. (2007) found evidence from starcounts for a 'long bar' at an angle of ∼ 45 , apparently tilted relative to and ◦ therefore distinct from the barred bulge. Martinez-Valpuesta &amp; Gerhard (2011) used an N-body model to argue that the long bar may nonetheless be the planar part of the same Galactic bar that contains the peanut bulge (see also Romero-G´mez et al. 2011). o In their model, the outer parts of the bar oscillate from trailing to leading, coupling with nearby spiral arms which rotate at lower pattern speed (Tagger et al. 1987). Future APOGEE kinematic observations in this region (e.g. Nidever et al. 2012) will test the model predictions and shed more light on the dynamical transition between bulge, two-dimensional bar, and inner disk.

## 3. Bulges in cosmological disk galaxy simulations

Galaxies in the hierarchical universe grow continuously by accretion of matter, and they are subjected to a variety of external perturbations, for example, from merging satellites. Therefore, eventually the idealized disk simulations will need to be superceded by more realistic disk galaxy formation simulations embedded in a cosmological setting. It has only become possible recently to follow realistic high-resolution, fully cosmological simulations from high redshift to now. These models are already stimulating our understanding of bulge formation, but they still have considerable uncertainties in how the star formation and feedback processes are modelled.

As a reference point, I start with a brief discussion of the dissipative collapse model by Samland &amp; Gerhard (2003). This may be viewed as a modern version of the dissipative collapse envisaged by Eggen et al. (1962). The model considered the formation of a large disk galaxy in a 2 × 10 12 M , spinning ( /circledot λ = 0 05) halo, with accre-. tion history taken from cosmological simulations. Star formation and feedback were modelled by a 2-phase fluid cloud fluid model calibrated on observed star formation / rates, and the enrichment of Fe and α -elements was tracked separately. Axisymmetric collapse was followed excluding all minor mergers.

In this dissipative collapse, the galaxy formed from inside out. The star formation history closely followed the gas accretion history. The disk grew from initially small, thick and dynamically hot to a large, thin and cold disk at redshift z = 0. The final bulge consisted of at least two stellar populations. One was formed in the early rapid collapse phase, the second was formed later in the disk which grew bar unstable. They di ff ered by their [ α / Fe]-ratios; the final bulge contained [ α / Fe]-enhanced stars for a range of [Fe H], and also some super-solar stars made from gas channeled inwards by the bar. / This illustrates the chemodynamical signatures which can be used to understand the origin of the Milky Way bulge from the new data.

Obreja et al. (2013) compared the properties of bulges made in some of the recent cosmological simulations, from Brook et al. (2012); Dom´ enech-Moral et al. (2012). Generally, these simulations employ strong feedback to prevent overproduction of stars at early times. As in the dissipative collapse model, the star formation rate (SFR) in these simulations closely follows the mass accretion rate, and the formation history can be divided into an early starburst-collapse phase and a later phase with lower SFR

driven by disk instabilities and minor mergers. The stellar populations show corresponding age and metallicity distributions, with the old population being more metal poor and α -enhanced. Most of the old population found in the final bulge formed in disjoint places along filaments at high redshift (as opposed to in the early collapsing halo in the dissipative collapse model).

One can associate the rapid phase with a classical bulge, and the late phase with a box peanut or disky bulge. / Such a classical bulge is clearly di ff erent from a bulge formed by a late major merger. In all of the five simulations analyzed by Obreja et al. (2013), both components were present, with di ff erent mass ratio. Guedes et al. (2013) and Okamoto (2013) followed the evolution of such 'pseudo-bulges' in very high resolution simulations. The structure of these bulges is characterized by a Sersic index around n = 1 ± 0 5, . possibly depending on the feedback model. One important result from these simulations is that, because of the inside-out formation, the inner bulge region can form early and consist of predominantly very old stars, whereas the disk further out has a much more uniform star formation history. In the Eris simulation of Guedes et al, the bulge is built from the early, low-angular momentum part of the disk. This model does not have a final box peanut bulge like in the Milky Way, but it illus-/ trates how the bulge stars could be very old as is observed in Baade's window in the Galactic bulge.

## References

Athanassoula, E. 2002, ApJ, 569, L83. arXiv:astro-ph/0203368

-2005, MNRAS, 358, 1477. arXiv:astro-ph/0502316

Benjamin, R. A., Churchwell, E., Babler, B. L., Indebetouw, R., &amp; et al. 2005, ApJ, 630, L149. arXiv:astro-ph/0508325

Binney, J., Gerhard, O., &amp; Spergel, D. 1997, MNRAS, 288, 365. arXiv:astro-ph/9609066

Bissantz, N., Englmaier, P., &amp; Gerhard, O. 2003, MNRAS, 340, 949. arXiv:astro-ph/ 0212516

Bissantz, N., &amp; Gerhard, O. 2002, MNRAS, 330, 591. arXiv:astro-ph/0110368

Blitz, L., &amp; Spergel, D. N. 1991, ApJ, 379, 631

Brook, C. B., Stinson, G., Gibson, B. K., Roˇ skar, R., Wadsley, J., &amp; Quinn, T. 2012, MNRAS, 419, 771. 1105.2562

Cabrera-Lavers, A., Hammersley, P. L., Gonz´ alez-Fern´ an dez, C., L´pez-Corredoira, M., &amp; et o al. 2007, A&amp;A, 465, 825. arXiv:astro-ph/0702109

Combes, F., Debbasch, F., Friedli, D., &amp; Pfenniger, D. 1990, A&amp;A, 233, 82

Combes, F., &amp; Sanders, R. H. 1981, A&amp;A, 96, 164

Debattista, V. P., Mayer, L., Carollo, C. M., Moore, B., Wadsley, J., &amp; Quinn, T. 2006, ApJ, 645, 209. arXiv:astro-ph/0509310

Dom´ enech-Moral, M., Mart´ ınez-Serrano, F. J., Dom´ ıngue z-Tenreiro, R., &amp; Serna, A. 2012,

MNRAS, 421, 2510. 1201.2641

Dwek, E., Arendt, R. G., Hauser, M. G., Kelsall, T., &amp; et al. 1995, ApJ, 445, 716

Eggen, O. J., Lynden-Bell, D., &amp; Sandage, A. R. 1962, ApJ, 136, 748

Englmaier, P., &amp; Gerhard, O. 1999, MNRAS, 304, 512. arXiv:astro-ph/9810208

Fux, R. 1999, A&amp;A, 345, 787. arXiv:astro-ph/9903154

Gerhard, O., &amp; Martinez-Valpuesta, I. 2012, ApJ, 744, L8. 1112.0179

Gonzalez, O. A., Rejkuba, M., Minniti, D., Zoccali, M., Valenti, E., &amp; Saito, R. K. 2011, A&amp;A, 534, L14. 1110.0925

Gonzalez, O. A., Rejkuba, M., Zoccali, M., Valent, E., Minniti, D., &amp; Tobar, R. 2013, A&amp;A, 552, A110. 1302.0243

Guedes, J., Mayer, L., Carollo, M., &amp; Madau, P. 2013, ApJ, 772, 36. 1211.1713

Hammersley, P. L., Garz´ on, F., Mahoney, T. J., L´pez-Corr edoira, M., &amp; Torres, M. A. P. 2000, o MNRAS, 317, L45. arXiv:astro-ph/0007232

Johnson, C. I., Rich, R. M., Fulbright, J. P., Valenti, E., &amp; McWilliam, A. 2011, ApJ, 732, 108. 1103.2143

Kunder, A., Koch, A., Rich, R. M., de Propris, R., Howard, C. D., Stubbs, S. A., Johnson, C. I., Shen, J., Wang, Y., Robin, A. C., Kormendy, J., Soto, M., Frinchaboy, P., Reitzel, D. B., Zhao, H., &amp; Origlia, L. 2012, AJ, 143, 57. 1112.1955

Li, Z.-Y., &amp; Shen, J. 2012, ApJ, 757, L7. 1206.3779

L´ opez-Corredoira, M., Cabrera-Lavers, A., &amp; Gerhard, O. E. 2005, A&amp;A, 439, 107. arXiv: astro-ph/0504608

Martinez-Valpuesta, I., &amp; Gerhard, O. 2011, ApJ, 734, L20 + . 1105.0928

-2013, ApJ, 766, L3. 1302.1613

Martinez-Valpuesta, I., Shlosman, I., &amp; Heller, C. 2006, ApJ, 637, 214. arXiv:astro-ph/ 0507219

McWilliam, A., &amp; Zoccali, M. 2010, ApJ, 724, 1491. 1008.0519

Nataf, D. M., Udalski, A., Gould, A., Fouqu´ e, P., &amp; Stanek, K. Z. 2010, ApJ, 721, L28. 1007. 5065

Ness, M., Freeman, K., Athanassoula, E., Wylie-de-Boer, E., Bland-Hawthorn, J., Asplund, M., Lewis, G. F., Yong, D., Lane, R. R., &amp; Kiss, L. L. 2013a, MNRAS, 430, 836. 1212.1540

Ness, M., Freeman, K., Athanassoula, E., Wylie-de-Boer, E., Bland-Hawthorn, J., Asplund, M., Lewis, G. F., Yong, D., Lane, R. R., Kiss, L. L., &amp; Ibata, R. 2013b, MNRAS, 432, 2092. 1303.6656

Ness, M., Freeman, K., Athanassoula, E., Wylie-De-Boer, E., Bland-Hawthorn, J., Lewis, G. F., Yong, D., Asplund, M., Lane, R. R., Kiss, L. L., &amp; Ibata, R. 2012, ApJ, 756, 22. 1207. 0888

Nidever, D. L., Zasowski, G., Majewski, S. R., Bird, J., &amp; et al. 2012, ApJ, 755, L25. 1207. 3797

Nishiyama, S., Nagata, T., Baba, D., Haba, Y., &amp; et al. 2005, ApJ, 621, L105. arXiv: astro-ph/0502058

Obreja, A., Dom´ ınguez-Tenreiro, R., Brook, C., Mart´ ınez -Serrano, F. J., Dom´nech-Moral, M., e

Serna, A., Moll´ a, M., &amp; Stinson, G. 2013, ApJ, 763, 26. 1211.3906

Okamoto, T. 2013, MNRAS, 428, 718. 1203.5372

Pfenniger, D., &amp; Friedli, D. 1991, A&amp;A, 252, 75

Raha, N., Sellwood, J. A., James, R. A., &amp; Kahn, F. D. 1991, Nat, 352, 411

Rodriguez-Fernandez, N. J., &amp; Combes, F. 2008, A&amp;A, 489, 115. 0806.4252

Romero-G´ omez, M., Athanassoula, E., Antoja, T., &amp; Figueras, F. 2011, MNRAS, 418, 1176. 1108.0660

Saha, K., &amp; Gerhard, O. 2013, MNRAS, 430, 2039. 1212.2317

Saha, K., Martinez-Valpuesta, I., &amp; Gerhard, O. 2012, MNRAS, 421, 333. 1105.5797

Saito, R. K., Zoccali, M., McWilliam, A., Minniti, D., Gonzalez, O. A., &amp; Hill, V. 2011, AJ, 142, 76. 1107.5360

Samland, M., &amp; Gerhard, O. E. 2003, A&amp;A, 399, 961. arXiv:astro-ph/0301499 Sellwood, J. A. 1981, A&amp;A, 99, 362

Shen, J., Rich, R. M., Kormendy, J., Howard, C. D., De Propris, R., &amp; Kunder, A. 2010, ApJ, 720, L72. 1005.0385

Skrutskie, M. F., Cutri, R. M., Stiening, R., Weinberg, M. D., &amp; et al. 2006, AJ, 131, 1163 Stanek, K. Z., Mateo, M., Udalski, A., Szymanski, M., Kaluzny, J., &amp; Kubiak, M. 1994, ApJ,

429, L73. arXiv:astro-ph/9404026

Tagger, M., Sygnet, J. F., Athanassoula, E., &amp; Pellat, R. 1987, ApJ, 318, L43

Wegg, C., &amp; Gerhard, O. 2013, ArXiv e-prints. 1308.0593

Zoccali, M., Hill, V ., Lecureur, A., Barbuy, B., Renzini, A., Minniti, D., G´mez, A., &amp; Ortolani, o S. 2008, A&amp;A, 486, 177. 0805.1218
Ortwin Gerhard
Formation Models of the Galactic Bulge
The Galactic bulge is now considered to be the inner three-dimensional part
of the Milky Way's bar. It has a peanut shape and is characterized by
cylindrical rotation. In N-body simulations, box/peanut bulges arise from disks
through bar and buckling instabilities. Models of this kind explain much of the
structure and kinematics of the Galactic bulge and, in principle, also its
vertical metallicity gradient. Cosmological disk galaxy formation models with
high resolution and improved feedback models are now able to generate late-type
disk galaxies with disk-like or barred bulges. These bulges often contain an
early collapse stellar population and a population driven by later disk
instabilities. Due to the inside-out disk formation, these bulges can be
predominantly old, similar to the Milky Way bulge.