Bulges are commonly believed to form in the dynamical violence of galaxy collisions and mergers. Here we model the stellar kinematics of the Bulge Radial Velocity Assay (BRAVA), and find no sign that the Milky Way contains a classical bulge formed by scrambling pre-existing disks of stars in major mergers. Rather, the bulge appears to be a bar, seen somewhat end-on, as hinted from its asymmetric boxy shape. We construct a simple but realistic N-body model of the Galaxy that self-consistently develops a bar. The bar immediately buckles and thickens in the vertical direction. As seen from the Sun, the result resembles the boxy bulge of our Galaxy. The model fits the BRAVA stellar kinematic data covering the whole bulge strikingly well with no need for a merger-made classical bulge. The bar in our best fit model has a half-length of ∼ 4 kpc and extends 20 • from the Sun-Galactic Center line. We use the new kinematic constraints to show that any classical bulge contribution cannot be larger than ∼ 8% of the disk mass. Thus the Galactic bulge is a part of the disk and not a separate component made in a prior merger. Giant, pure-disk galaxies like our own present a major challenge to the standard picture in which galaxy formation is dominated by hierarchical clustering and galaxy mergers.

More than two thirds of disk galaxies are barred to some degree. Many today harbor massive concentrations of gas in their centers, and some are known to possess supermassive black holes (SMBHs) and their associated stellar cusps. Previous theoretical work has suggested that a bar in a galaxy could be dissolved by the formation of a mass concentration in the center, although the precise mass and degree of central concentration required is not well-established. We report an extensive study of the effects of central masses on bars in high-quality N-body simulations of galaxies. We have varied the growth rate of the central mass, its final mass and degree of concentration to examine how these factors affect the evolution of the bar. Our main conclusions are: (1) Bars are more robust than previously thought. The central mass has to be as large as several percent of the disk mass to completely destroy the bar on a short timescale. (2) For a given mass, dense objects cause the greatest reduction in bar amplitude, while significantly more diffuse objects have a lesser effect. (3) The bar amplitude always decreases as the central mass is grown, and continues to decay thereafter on a cosmological time-scale. (4) The first phase of bar-weakening is due to the destruction by the CMC of lower-energy, bar-supporting orbits, while the second phase is a consequence of secular changes to the global potential which further diminish the number of bar-supporting orbits. We provide detailed phase-space and orbit analysis to support this suggestion. Thus current masses of SMBHs are probably too small, even when dressed with a stellar cusp, to affect the bar in their host galaxies. The molecular gas concentrations found in some barred galaxies are also too diffuse to affect the amplitude of the bar significantly.Comment: AASTeX v5.0 preprint; 44 pages, including 1 table and 16 figures. To appear in ApJ. High resolution version can be found at http://www.physics.rutgers.edu/~shen/bar_destruct/paper_high_res.pd

We apply the axisymmetric orbit superposition modeling to estimate the mass of the supermassive black hole and dark matter halo profile of NGC 4649. We have included data sets from the Hubble Space Telescope, stellar, and globular cluster observations. Our modeling gives M • = (4.5 ± 1.0) × 10 9 M ⊙ and M/L V, obs = 8.7 ± 1.0 (or M/L V = 8.0 ± 0.9 after foreground Galactic extinction is corrected). We confirm the presence of a dark matter halo, but the stellar mass dominates inside the effective radius. The parameters of the dark halo are less constrained due to the sparse globular cluster data at large radii. We find that in NGC 4649 the dynamical mass profile from our modeling is consistently larger than that derived from the X-ray data over most of the radial range by roughly 60% to 80%. It implies that either some forms of non-thermal pressure need to be included, the assumed hydrostatic equilibrium may not be a good approximation in the X-ray modelings of NGC 4649, or our assumptions used in the dynamical models are biased. Our new M • is about two times larger than the previous published value; the earlier model did not adequately sample the orbits required to match the large tangential anisotropy in the galaxy center. If we assume that there is no dark matter, the results on the black hole mass and M/L V, obs do not change significantly, which we attribute to the inclusion of HST spectra, the sparse globular cluster kinematics, and a diffuse dark matter halo. Without the HST data, the significance of the black hole detection is greatly reduced.

Recent ideas for the origin and persistence of the warps commonly observed in disc galaxies have focused on cosmic infall. We present N‐body simulations of an idealized form of cosmic infall on to a disc galaxy and obtain a warp that closely resembles those observed. The inner disc tilts remarkably rigidly, indicating strong cohesion due to self‐gravity. The line of nodes (LON) of the warp inside R26.5∼ 4.5Rd is straight, while that beyond R26.5 generally forms a loosely wound, leading spiral in agreement with Briggs's rules. We focus on the mechanism of the warp and show that the leading spiral arises from the torques from the misaligned inner disc and its associated inner oblate halo. The fact that the LON of most warps forms a leading spiral might imply that the disc mass is significant in the centre. If the LON can be traced to very large radii in future observations, it may reveal information on the mass distribution of the outer halo. The warp is not strongly damped by the halo because the precession rate of the inner disc is slow and the inner halo generally remains aligned with the inner disc. Thus, even after the imposed quadrupolar perturbation is removed, the warp persists for a few Gyr, by which time another infall event can be expected.

We present gas flow models for the Milky Way based on high-resolution grid-based hydrodynamical simulations. The basic galactic potential we use is from a N-body model constrained by the density of red clump giants in the Galactic bulge. We augment this potential with a nuclear bulge, two pairs of spiral arms and additional mass at the bar end to represent the long bar component. With this combined model we can reproduce many features in the observed (l, v) diagram with a bar pattern speed of 33 km s −1 kpc −1 and a spiral pattern speed of 23 km s −1 kpc −1 . The shape and kinematics of the nuclear ring, Bania's Clump 2, the Connecting arm, the Near and Far 3-kpc arms, the Molecular Ring, and the spiral arm tangent points in our simulations are comparable to those in the observations. Our results imply that a low pattern speed model for the bar in our Milky Way reproduces the observations for a suitable Galactic potential. Our best model gives a better match to the (l, v) diagram than previous high pattern speed hydrodynamical simulations.

A vertical X-shaped structure in the Galactic bulge was recently reported. Here we present evidence of a similar X-shaped structure in the Shen et al. (2010) bar/boxy bulge model that simultaneously matches the stellar kinematics successfully. The X-shaped structure is found in the central region of our bar/boxy bulge model, and is qualitatively consistent with the observed one in many aspects. End-to-end separations of the X-shaped structure in the radial and vertical directions are roughly 3 kpc and 1.8 kpc, respectively. The X-shaped structure contains about 7% of light in the boxy bulge region, but it is significant enough to be identified in observations. An X-shaped structure naturally arises in the formation of bar/boxy bulges, and is mainly associated with orbits trapped around the vertically-extended x 1 family. Like the bar in our model, the X-shaped structure tilts away from the Sun-Galactic center line by 20 • . The X-shaped structure becomes increasingly symmetric about the disk plane, so the observed symmetry may indicate that it formed at least a few billion years ago. The existence of the vertical X-shaped structure suggests that the formation of the Milky Way bulge is shaped mainly by internal disk dynamical instabilities.

Dust lanes, nuclear rings, and nuclear spirals are typical gas structures in the inner region of barred galaxies. Their shapes and properties are linked to the physical parameters of the host galaxy. We use high-resolution hydrodynamical simulations to study 2D gas flows in simple barred galaxy models. The nuclear rings formed in our simulations can be divided into two groups: one group is nearly round and the other is highly elongated. We find that roundish rings may not form when the bar pattern speed is too high or the bulge central density is too low. We also study the periodic orbits in our galaxy models, and find that the concept of inner Lindblad resonance (ILR) may be generalized by the extent of x 2 orbits. All roundish nuclear rings in our simulations settle in the range of x 2 orbits (or ILRs). However, knowing the resonances is insufficient to pin down the exact location of these nuclear rings. We suggest that the backbone of round nuclear rings is the x 2 orbital family, i.e., round nuclear rings are allowed only in the radial range of x 2 orbits. A round nuclear ring forms exactly at the radius where the residual angular momentum of infalling gas balances the centrifugal force, which can be described by a parameter f ring measured from the rotation curve. The gravitational torque on gas in high pattern speed models is larger, leading to a smaller ring size than in the low pattern speed models. Our result may have important implications for using nuclear rings to measure the parameters of real barred galaxies with 2D gas kinematics.

The anisotropy parameter β characterizes the extent to which orbits in stellar systems are predominantly radial or tangential, and is likely to constrain, for the stellar halo of the Milky Way, scenarios for its formation and evolution. We have measured the anisotropy β as a function of Galactocentric radius from 5 − 100 kpc for over 8600 metal poor ([Fe/H] < −1.3) halo K giants from the LAMOST catalog with line-of-sight velocities and distances, matched to proper motions from the second Gaia data release. We construct full 6-D positions and velocities for the K giants to directly measure the 3 components of the velocity dispersion (σ r , σ θ , σ φ ) (in spherical coordinates). We find that the orbits in the halo are radial over our full Galactocentric distance range reaching over 100 kpc. The anisotropy remains remarkably unchanged with Galactocentric radius from approximately 5 to 25 kpc, with an amplitude that depends on the metallicity of the stars, dropping from β ≈ 0.9 for −1.8 ≤ [Fe/H] < −1.3 (for the bulk of the stars) to β ≈ 0.6 for the lowest metallicities ([Fe/H] < −1.8). Considering our sample as a whole, β ≈ 0.8 and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to β < 0.3 past 100 kpc. Within 8 kpc, β < 0.8. The measurement of anisotropy is affected by substructure and streams, particularly beyond a Galactocentric distance of approximately 25 kpc, where the Sagittarius stream is prominent in the data. These results are complimentary to recent analysis of simulations by Loebman et al. and of SDSS/Gaia DR1 data by Belokurov et al.