YOUNG STARS IN AN OLD BULGE: A NATURAL OUTCOME OF INTERNAL EVOLUTION IN THE MILKY WAY

The center of our disk galaxy, the Milky Way, is dominated by a boxy/peanut-shaped bulge. Numerous studies of the bulge based on stellar photometry have concluded that the bulge stars are exclusively old. The perceived lack of young stars in the bulge strongly constrains its likely formation scenarios, providing evidence that the bulge is a unique population that formed early and separately from the disk. However, recent studies of individual bulge stars using the microlensing technique have reported that they span a range of ages, emphasizing that the bulge may not be a monolithic structure. In this Letter we demonstrate that the presence of young stars that are located predominantly nearer to the plane is expected for a bulge that has formed from the disk via dynamical instabilities. Using an N-body+ smoothed particle hydrodynamics simulation of a disk galaxy forming out of gas cooling inside a dark matter halo and forming stars, we find a qualitative agreement between our model and the observations of younger metal-rich stars in the bulge. We are also able to partially resolve the apparent contradiction in the literature between results that argue for a purely old bulge population and those that show a population comprised of a range in ages; the key is where to look

Forming double-barred galaxies from dynamically cool inner disks

About one-third of early-type barred galaxies host small-scale secondary bars. The formation and evolution of such double-barred (S2B) galaxies remain far from being well understood. In order to understand the formation of such systems, we explore a large parameter space of isolated pure-disk simulations. We show that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady secondary bar while the outer disk forms a large-scale primary bar. The independent bar instabilities of inner and outer disks result in long-lived double-barred structures whose dynamical properties are comparable to those in observations. This formation scenario indicates that the secondary bar might form from the general bar instability, the same as the primary bar. Under some circumstances, the interaction of the bars and the disk leads to the two bars aligning or single, nuclear, bars only. Simulations that are cool enough of the center to experience clump instabilities may also generate steady S2B galaxies. In this case, the secondary bars are “fast,” i.e., the bar length is close to the co-rotation radius. This is the first time that S2B galaxies containing a fast secondary bar are reported. Previous orbit-based studies had suggested that fast secondary bars were not dynamically possibl

Imprints of radial migration on the Milky Way’s metallicity distribution functions

Recent analysis of the SDSS-III/Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 12 stellar catalog has revealed that the Milky Way’s (MW) metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high-resolution, N-body+SPH simulation, we show that the changing skewness arises from radial migration—metal-rich stars form in the inner disk and subsequently migrate to the metal-poorer outer disk. These migrated stars represent a large fraction (>50%) of the stars in the outer disk; they populate the high-metallicity tail of the MDFs and are, in general, more metal-rich than the surrounding outer disk gas. The simulation also reproduces another surprising APOGEE result: the spatially invariant high-[α/Fe] MDFs. This arises in the simulation from the migration of a population formed within a narrow range of radii (3.2 ±1.2 kpc) and time (8.8 ± 0.6 Gyr ago), rather than from spatially extended star formation in a homogeneous medium at early times. These results point toward the crucial role radial migration has played in shaping our MW

Bar-Halo Friction in Galaxies II: Metastability

It is well-established that strong bars rotating in dense halos generally slow down as they lose angular momentum to the halo through dynamical friction. Angular momentum exchanges between the bar and halo particles take place at resonances. While some particles gain and others lose, friction arises when there is an excess of gainers over losers. This imbalance results from the generally decreasing numbers of particles with increasing angular momentum, and friction can therefore be avoided if there is no gradient in the density of particles across the major resonances. Here we show that anomalously weak friction can occur for this reason if the pattern speed of the bar fluctuates upwards. After such an event, the density of resonant halo particles has a local inflexion created by the earlier exchanges, and bar slowdown can be delayed for a long period; we describe this as a metastable state. We show that this behavior in purely collisionless N-body simulations is far more likely to occur in methods with adaptive resolution. We also show that the phenomenon could arise in nature, since bar-driven gas inflow could easily raise the bar pattern speed enough to reach the metastable state. Finally, we demonstrate that mild external, or internal, perturbations quickly restore the usual frictional drag, and it is unlikely therefore that a strong bar in a galaxy having a dense halo could rotate for a long period without friction.Comment: 13 pages, 11 figures, to appear in Ap

Direct Confirmation of Two Pattern Speeds in the Double Barred Galaxy NGC 2950

We present surface photometry and stellar kinematics of NGC 2950, which is a nearby and undisturbed SB0 galaxy hosting two nested stellar bars. We use the Tremaine-Weinberg method to measure the pattern speed of the primary bar. This also permits us to establish directly and for the first time that the two nested bars are rotating with different pattern speeds, and in particular that the rotation frequency of the secondary bar is higher than that of the primary one.Comment: 12 pages, 4 figures. To appear in ApJ Letter

A unified framework for the orbital structure of bars and triaxial ellipsoids

We examine a large random sample of orbits in two self-consistent simulations of N-body bars. Orbits in these bars are classified both visually and with a new automated orbit classification method based on frequency analysis. The well-known prograde x1 orbit family originates from the same parent orbit as the box orbits in stationary and rotating triaxial ellipsoids. However, only a small fraction of bar orbits (~4%) have predominately prograde motion like their periodic parent orbit. Most bar orbits arising from the x1 orbit have little net angular momentum in the bar frame, making them equivalent to box orbits in rotating triaxial potentials. In these simulations a small fraction of bar orbits (~7%) are long-axis tubes that behave exactly like those in triaxial ellipsoids: they are tipped about the intermediate axis owing to the Coriolis force, with the sense of tipping determined by the sign of their angular momentum about the long axis. No orbits parented by prograde periodic x2 orbits are found in the pure bar model, but a tiny population (~2%) of short-axis tube orbits parented by retrograde x4 orbits are found. When a central point mass representing a supermassive black hole (SMBH) is grown adiabatically at the center of the bar, those orbits that lie in the immediate vicinity of the SMBH are transformed into precessing Keplerian orbits that belong to the same major families (short-axis tubes, long-axis tubes and boxes) occupying the bar at larger radii. During the growth of an SMBH, the inflow of mass and outward transport of angular momentum transform some x1 and long-axis tube orbits into prograde short-axis tubes. This study has important implications for future attempts to constrain the masses of SMBHs in barred galaxies using orbit-based methods like the Schwarzschild orbit superposition scheme and for understanding the observed features in barred galaxies

Warped Galaxies From Misaligned Angular Momenta

A galaxy disk embedded in a rotating halo experiences a dynamical friction force which causes it to warp when the angular momentum axes of the disk and halo are misaligned. Our fully self-consistent simulations of this process induce long-lived warps in the disk which mimic Briggs's rules of warp behavior. They also demonstrate that random motion within the disk adds significantly to its stiffness. Moreover, warps generated in this way have no winding problem and are more pronounced in the extended \h1 disk. As emphasized by Binney and his co-workers, angular momentum misalignments, which are expected in hierarchical models of galaxy formation, can account for the high fraction of warped galaxies. Our simulations exemplify the role of misaligned spins in warp formation even when the halo density is not significantly flattened.Comment: 6 pages, 5 figures. Accepted for publication in Ap.J.

Dynamical Friction and the Distribution of Dark Matter in Barred Galaxies

We use fully self-consistent N-body simulations of barred galaxies to show that dynamical friction from a dense dark matter halo dramatically slows the rotation rate of bars. Our result supports previous theoretical predictions for a bar rotating within a massive halo. On the other hand, low density halos, such as those required for maximum disks, allow the bar to continue to rotate at a high rate. There is somewhat meager observational evidence indicating that bars in real galaxies do rotate rapidly and we use our result to argue that dark matter halos must have a low central density in all high surface brightness disk galaxies, including the Milky Way. Bars in galaxies that have larger fractions of dark matter should rotate slowly, and we suggest that a promising place to look for such candidate objects is among galaxies of intermediate surface brightness.Comment: 6 pages, Latex, 3 figures, Accepted by Ap.J.L., revised copy, includes an added paragrap

Disk galaxies with broken luminosity profiles from cosmological simulations

We present SPH cosmological simulations of the formation of three disk galaxies with a detailed treatment of chemical evolution and cooling. The resulting galaxies have properties compatible with observations: relatively high disk-to-total ratios, thin stellar disks and good agreement with the Tully-Fisher and the luminosity-size relations. They present a break in the luminosity profile at 3.0 +- 0.5 disk scale lengths, while showing an exponential mass profile without any apparent breaks, in line with recent observational results. Since the stellar mass profile is exponential, only differences in the stellar populations can be the cause of the luminosity break. Although we find a cutoff for the star formation rate imposed by a density threshold in our star formation model, it does not coincide with the luminosity break and is located at 4.3 +- 0.4 disk scale lengths, with star formation going on between both radii. The color profiles and the age profiles are "U-shaped", with the minimum for both profiles located approximately at the break radius. The SFR to stellar mass ratio increases until the break, explaining the coincidence of the break with the minimum of the age profile. Beyond the break we find a steep decline in the gas density and, consequently, a decline in the SFR and redder colors. We show that most stars (64-78%) in the outer disk originate in the inner disk and afterwards migrate there. Such stellar migrations are likely the main origin of the U-shaped age profile and, therefore, of the luminosity break.Comment: 8 pages, 4 figures. Accepted by ApJ