Source author record

Juntai Shen

Juntai Shen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2022arXiv

3D intrinsic shapes of quiescent galaxies in observations and simulations

We study the intrinsic 3D shapes of quiescent galaxies over the last half of cosmic history based on their axial ratio distribution. To this end, we construct a sample of unprecedented size, exploiting multi-wavelength $u$-to-$K_s$ photometry from the deep wide area surveys KiDS+VIKING paired with high-quality $i$-band imaging from HSC-SSP. Dependencies of the shapes on mass, redshift, photometric bulge prominence and environment are considered. For comparison, the intrinsic shapes of quenched galaxies in the IllustrisTNG simulations are analyzed and contrasted to their formation history. We find that over the full $0<z<0.9$ range, and in both simulations and observations, spheroidal 3D shapes become more abundant at $M_* > 10^{11}\ M_{\odot}$, with the effect being most pronounced at lower redshifts. In TNG, the most massive galaxies feature the highest ex-situ stellar mass fractions, pointing to violent relaxation via mergers as the mechanism responsible for their 3D shape transformation. Larger differences between observed and simulated shapes are found at low to intermediate masses. At any mass, the most spheroidal quiescent galaxies in TNG feature the highest bulge mass fractions, and conversely observed quiescent galaxies with the highest bulge-to-total ratios are found to be intrinsically the roundest. Finally, we detect an environmental influence on galaxy shape, at least at the highest masses, such that at fixed mass and redshift quiescent galaxies tend to be rounder in denser environments.

preprint2022arXiv

Blanco DECam Bulge Survey (BDBS) IV: Metallicity Distributions and Bulge Structure from 2.6 Million Red Clump Stars

We present photometric metallicity measurements for a sample of 2.6 million bulge red clump stars extracted from the Blanco DECam Bulge Survey (BDBS). Similar to previous studies, we find that the bulge exhibits a strong vertical metallicity gradient, and that at least two peaks in the metallicity distribution functions appear at b < -5. We can discern a metal-poor ([Fe/H] ~ -0.3) and metal-rich ([Fe/H] ~ +0.2) abundance distribution that each show clear systematic trends with latitude, and may be best understood by changes in the bulge's star formation/enrichment processes. Both groups exhibit asymmetric tails, and as a result we argue that the proximity of a star to either peak in [Fe/H] space is not necessarily an affirmation of group membership. The metal-poor peak shifts to lower [Fe/H] values at larger distances from the plane while the metal-rich tail truncates. Close to the plane, the metal-rich tail appears broader along the minor axis than in off-axis fields. We also posit that the bulge has two metal-poor populations -- one that belongs to the metal-poor tail of the low latitude and predominantly metal-rich group, and another belonging to the metal-poor group that dominates in the outer bulge. We detect the X-shape structure in fields with |Z| > 0.7 kpc and for stars with [Fe/H] > -0.5. Stars with [Fe/H] < -0.5 may form a spheroidal or "thick bar" distribution while those with [Fe/H] > -0.1 are strongly concentrated near the plane.

preprint2022arXiv

Large-scale Hydrodynamical Shocks as the Smoking Gun Evidence for a Bar in M31

The formation and evolutionary history of M31 are closely related to its dynamical structures, which remain unclear due to its high inclination. Gas kinematics could provide crucial evidence for the existence of a rotating bar in M31. Using the position-velocity diagram of [OIII] and HI, we are able to identify clear sharp velocity jump (shock) features with a typical amplitude over 100 km/s in the central region of M31 (4.6 kpc X 2.3 kpc, or 20 arcmin X 10 arcmin). We also simulate gas morphology and kinematics in barred M31 potentials and find that the bar-induced shocks can produce velocity jumps similar to those in [OIII]. The identified shock features in both [OIII] and HI are broadly consistent, and they are found mainly on the leading sides of the bar/bulge, following a hallmark pattern expected from the bar-driven gas inflow. Shock features on the far side of the disk are clearer than those on the near side, possibly due to limited data coverage on the near side, as well as obscuration by the warped gas and dust layers. Further hydrodynamical simulations with more sophisticated physics are desired to fully understand the observed gas features and to better constrain the parameters of the bar in M31.

preprint2022arXiv

Milky Way Mass with K Giants and BHB Stars Using LAMOST, SDSS/SEGUE, and Gaia: 3D Spherical Jeans Equation and Tracer Mass Estimator

We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it Gaia} DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to $\sim70$ and $\sim50$ kpc, the Jeans method yields total enclosed masses of $4.3\pm0.95$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ and $4.1\pm1.2$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ for the KG and BHB stars, respectively. For the KG and BHB samples we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) $\pm0.083$ (systematic) $\times10^{12}$ M$_\odot$ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) $\pm0.15$ (systematic) $\times10^{12}$ M$_\odot$, respectively.

preprint2022arXiv

The North/South Asymmetry of the Galaxy: Possible Connection to the Vertical Phase Space Snail

The Galaxy is found to be in disequilibrium based on recent findings of the North/South (N/S) asymmetry and the phase mixing signatures, such as a phase spiral (snail) structure in the vertical phase space ($z-V_{z}$). We show that the N/S asymmetry in a tracer population of dwarfs may be quantitatively modeled with a simple phase snail model superimposed on a smooth equilibrium background. As the phase snail intersects with the $z$ axis, the number density is enhanced, and the velocity dispersion ($σ_{z}$) is decreased relative to the other side of the Galactic plane. Fitting only to the observed asymmetric N/S $σ_{z}$ profiles, we obtain reasonable parameters for the phase space snail and the potential utilized in modeling the background, despite the complex dependence of the model on the potential parameters and the significant selection effects of the data. Both the snail shape and the N/S number density difference given by our best-fit model are consistent with previous observations. The equilibrium background implies a local dark matter density of $0.0151^{+0.0050}_{-0.0051}$ ${\rm M}_{\odot}\,{\rm pc}^{-3}$. The vertical bulk motion of our model is similar to the observation, but with a $\sim$1.2 $\rm km\,s^{-1}$ shift. Our work demonstrates the strong correlation between the phase space snail and the N/S asymmetry. Future observational constraints will facilitate more comprehensive snail models to unravel the Milky Way potential and the perturbation history encoded in the snail feature.

preprint2021arXiv

Gas Dynamics in the Galaxy: Total Mass Distribution and the Bar Pattern Speed

Gas morphology and kinematics in the Milky Way contain key information for understanding the formation and evolution of our Galaxy. We present a high resolution hydrodynamical simulation based on a realistic barred Milky Way potential constrained by recent observations. Our model can reproduce most features in the observed longitude-velocity diagram, including the Central Molecular Zone, the Near and Far 3-kpc arms, the Molecular Ring, and the spiral arm tangents. It can also explain the non-circular motions of masers obtained by the recent BeSSeL2 survey. The central gas kinematics are consistent with a mass of $6.9\times10^8\; {\rm M}_{\odot}$ in the Nuclear Stellar Disk. Our model predicts the formation of an elliptical gaseous ring surrounding the bar, which is composed of the 3-kpc arms, Norma arm, and the bar-spiral interfaces. This ring is similar to those "inner" rings in some Milky Way analogs with a boxy/peanut-shaped bulge. The kinematics of gas near the solar neighbourhood are governed by the Local arm, which is induced by the four major stellar spiral arms. The bar pattern speed constrained by our gas model is $37.5-40\; {\rm km}\;{\rm s}^{-1}\;{\rm kpc}^{-1}$, corresponding to a corotation radius of $R_{\rm CR}=6.0-6.4\;{\rm kpc}$. The rotation curve of our model rises gently within the central $\sim5\;{\rm kpc}$, which is significantly less steep than those predicted by modern zoom-in cosmological simulations such as Auriga.

preprint2020arXiv

Anisotropy of the Milky Way's stellar halo using K giants from LAMOST and $Gaia$

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 $β\approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $β\approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $β\approx0.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.

preprint2020arXiv

Dissecting the Phase Space Snail Shell

The on-going vertical phase mixing, manifesting itself as a snail shell in the $Z-V_{Z}$ phase space, has been discovered with the Gaia DR2 data. To better understand the origin and properties of the phase mixing process, we study the vertical phase-mixing signatures in arches (including the classical ``moving groups'') of the $V_{R}-V_ϕ$ phase space near the Solar circle. Interestingly, the phase space snail shell exists only in the arches with $|V_ϕ - V_{\rm LSR}| \lesssim 30$ km/s, i.e., stars on dynamically ``colder'' orbits. The snail shell becomes much weaker and eventually disappears for increasingly larger radial action ($J_{R}$), quantifying the ``hotness'' of orbits. Thus one should pay closer attention to the colder orbits in future phase mixing studies. We also confirm that the Hercules stream has two branches (at fast and slow $V_ϕ$), which may not be explained by a single mechanism, since only the fast branch shows the prominent snail shell feature. The hotter orbits may have phase-wrapped away already due to the much larger dynamical range in radial variation to facilitate faster phase mixing. To explain the lack of a well-defined snail shell in the hotter orbits, the disk should have been perturbed at least $500$ Myr ago. Our results offer more support to the recent satellite-disk encounter scenario than the internal bar buckling perturbation scenario as the origin of the phase space mixing. Origin of the more prominent snail shell in the $V_ϕ$ color-coded phase space is also discussed.

preprint2020arXiv

Mapping the tilt of the Milky Way bulge velocity ellipsoids with ARGOS and $Gaia$ DR2

Until the recent advent of $Gaia$ Data Release 2 (DR2) and deep multi-object spectroscopy, it has been difficult to obtain 6-D phase space information for large numbers of stars beyond 4 kpc, in particular towards the Galactic centre, where dust and crowding effects are significant. In this study we combine line-of-sight velocities from the Abundances and Radial velocity Galactic Origins Survey (ARGOS) spectroscopic survey with proper motions from $Gaia$ DR2, to obtain a sample of $\sim$ 7,000 red clump stars with 3-D velocities. We perform a large scale stellar kinematics study of the Milky Way (MW) bulge to characterize the bulge velocity ellipsoids. We measure the tilt $l_{v}$ of the major-axis of the velocity ellipsoid in the radial-longitudinal velocity plane in 20 fields across the bulge. The tilt or vertex deviation, is characteristic of non-axisymmetric systems and a significant tilt is a robust indicator of non-axisymmetry or bar presence. We compare the observations to the predicted kinematics of an N-body boxy-bulge model formed from dynamical instabilities. In the model, the $l_{v}$ values are strongly correlated with the angle ($α$) between the bulge major-axis and the Sun-Galactic centre line-of-sight. We use a maximum likelihood method to obtain an independent measurement of $α$, from bulge stellar kinematics alone. The most likely value of $α$ given our model is $α= (29 \pm 3)^{\circ}$. In the Baade's window, the metal-rich stars display a larger vertex deviation ($l_{v} = -40^{\circ}$) than the metal-poor stars ($l_{v} = 10^{\circ}$) but we do not detect significant $l_{v}-$metallicity trends in the other fields.

preprint2020arXiv

Testing the Prediction of Fuzzy Dark Matter Theory in the Milky Way Center

The fuzzy dark matter model (FDM, also known as quantum wave dark matter model) argues that light bosons with a mass of $\sim10^{-22}{\;{\rm eV}}$ are a possible candidate for dark matter in the Universe. One of the most important predictions of FDM is the formation of a soliton core instead of a density cusp at the center of galaxies. If FDM is the correct theory of dark matter, then the predicted soliton core can help to form the Central Molecular Zone (CMZ) in the Milky Way. We present high-resolution hydrodynamical simulations of gas flow patterns to constrain the properties of the soliton core based on a realistic Milky Way potential. We find that a dense center is required to form a reasonable CMZ. The size and kinematics of the CMZ offer a relatively strong constraint on the inner enclosed mass profile of the Galaxy. If a soliton core is not considered, a compact nuclear bulge alone with a radially varying mass-to-light ratio can match the observed size and kinematics of the CMZ. A soliton core model with a mass of $\approx4.0\times10^8{\; {\rm M}_{\odot}}$ and a core radius of $\approx0.05{\;{\rm kpc}}$, together with a less massive nuclear bulge with a constant mass-to-light ratio, also agrees nicely with the current data. Such a FDM soliton core corresponds to a boson mass of $\sim2-7\times10^{-22}{\;{\rm eV}}$, which could be further constrained by the improved determination of the mass-to-light ratio in the Galactic center.

preprint2020arXiv

Understanding the velocity distribution of the Galactic Bulge with APOGEE and Gaia

We revisit the stellar velocity distribution in the Galactic bulge/bar region with APOGEE DR16 and {\it Gaia} DR2, focusing in particular on the possible high-velocity (HV) peaks and their physical origin. We fit the velocity distributions with two different models, namely with Gauss-Hermite polynomial and Gaussian mixture model (GMM). The result of the fit using Gauss-Hermite polynomials reveals a positive correlation between the mean velocity ($\bar{V}$) and the "skewness" ($h_{3}$) of the velocity distribution, possibly caused by the Galactic bar. The $n=2$ GMM fitting reveals a symmetric longitudinal trend of $|μ_{2}|$ and $σ_{2}$ (the mean velocity and the standard deviation of the secondary component), which is inconsistent to the $x_{2}$ orbital family predictions. Cold secondary peaks could be seen at $|l|\sim6^\circ$. However, with the additional tangential information from {\it Gaia}, we find that the HV stars in the bulge show similar patterns in the radial-tangential velocity distribution ($V_{\rm R}-V_{\rm T}$), regardless of the existence of a distinct cold HV peak. The observed $V_{\rm R}-V_{\rm T}$ (or $V_{\rm GSR}-μ_{l}$) distributions are consistent with the predictions of a simple MW bar model. The chemical abundances and ages inferred from ASPCAP and CANNON suggest that the HV stars in the bulge/bar are generally as old as, if not older than, the other stars in the bulge/bar region.

preprint2016arXiv

Gas dynamics in the Milky Way: a low pattern speed model

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\; {\rm km}\;{\rm s}^{-1}\;{\rm kpc}^{-1}$ and a spiral pattern speed of $23\; {\rm km}\;{\rm s}^{-1}\;{\rm 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.

preprint2016arXiv

Kinematic properties of double-barred galaxies: simulations vs. integral-field observations

Using high resolution $N$-body simulations, we recently reported that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady double-barred (S2B) structure. Here we study the kinematics of these S2B simulations, and compare them to integral-field observations from ATLAS$^{3D}$ and SAURON. We show that S2B galaxies exhibit several distinct kinematic features, namely: (1) significantly distorted isovelocity contours at the transition region between the two bars, (2) peaks in $σ_\mathrm{LOS}$ along the minor axis of inner bars, which we term "$σ$-humps", that are often accompanied by ring/spiral-like features of increased $σ_\mathrm{LOS}$, (3) $h_3-\bar{v}$ anti-correlations in the region of the inner bar for certain orientations, and (4) rings of positive $h_4$ when viewed at low inclinations. The most impressive of these features are the $σ$-humps, these evolve with the inner bar, oscillating in strength just as the inner bar does as it rotates relative to the outer bar. We show that, in cylindrical coordinates, the inner bar has similar streaming motions and velocity dispersion properties as normal large-scale bars, except for $σ_z$, which exhibits peaks on the minor axis, i.e., humps. These $σ_z$ humps are responsible for producing the $σ$-humps. For three well-resolved early-type S2Bs (NGC 2859, NGC 2950, and NGC 3941) and a potential S2B candidate (NGC 3384), the S2B model qualitatively matches the integral-field data well, including the "$σ$-hollows" previously identified. We also discuss the kinematic effect of a nuclear disk in S2Bs.

preprint2015arXiv

A high-velocity bulge RR Lyrae variable on a halo-like orbit

We report on the RR Lyrae variable star, MACHO 176.18833.411, located toward the Galactic bulge and observed within the data from the ongoing Bulge RR Lyrae Radial Velocity Assay (BRAVA-RR), which has the unusual radial velocity of -372 +- 8 km/s and true space velocity of -482 +- 22 km/s relative to the Galactic rest frame. Located less than 1 kpc from the Galactic center and toward a field at (l,b)=(3,-2.5), this pulsating star has properties suggesting it belongs to the bulge RR Lyrae star population yet a velocity indicating it is abnormal, at least with respect to bulge giants and red clump stars. We show that this star is most likely a halo interloper and therefore suggest that halo contamination is not insignificant when studying metal-poor stars found within the bulge area, even for stars within 1 kpc of the Galactic center. We discuss the possibility that MACHO 176.18833.411 is on the extreme edge of the bulge RR Lyrae radial velocity distribution, and also consider a more exotic scenario in which it is a runaway star moving through the Galaxy.

preprint2015arXiv

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

We examine a large random sample of orbits in self-consistent simulations of N-body bars. Orbits in the 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. 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 due 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 (PKOs) which belong the same major families (short axis tubes, long axis tubes and boxes) occupying the bar at larger radii. During the growth of a SMBH the inflow of mass and outward transport of angular momentum transforms 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.

preprint2015arXiv

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 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 with 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 double-barred 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 double-barred galaxies containing a fast secondary bar are reported. Previous orbit-based studies had suggested that fast secondary bars are not dynamically possible.

preprint2015arXiv

Hydrodynamical Simulations of Nuclear Rings in Barred Galaxies

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_{\rm 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.

preprint2015arXiv

Kinematics of the X-shaped Milky Way Bulge: Expectations from a Self-consistent N-body Model

We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage.

preprint2015arXiv

Mapping the Three-Dimensional "X-Shaped Structure" in Models of the Galactic Bulge

Numerical simulations have shown that the X-shaped structure in the Milky Way bulge can naturally arise from the bar instability and buckling instability. To understand the influence of the buckling amplitude on the morphology of the X-shape, we analyze three self-consistent numerical simulations of barred galaxies with different buckling amplitudes (strong, intermediate and weak). We derive the three-dimensional density with an adaptive kernel smoothing technique. The face-on iso-density surfaces are all elliptical, while in the edge-on view, the morphology of buckled bars transitions with increasing radius, from a central boxy core to a peanut bulge and then to an extended thin bar. Based on these iso-density surfaces at different density levels, we find no clear evidence for a well-defined structure shaped like a letter X. The X-shaped structure is more peanut-like, whose visual perception is probably enhanced by the pinched inner concave iso-density contours. The peanut bulge can reproduce qualitatively the observed bimodal distributions which were used as evidence for the discovery of the X-shape. The central boxy core is shaped like an oblong tablet, extending to $\sim$ 500 pc above and below the Galactic plane ($|b| \sim 4^\circ$). From the solar perspective, lines of sight passing through the central boxy core do not show bimodal distributions. This generally agrees with the observations that the double peaks merge at $|b| \sim 4^\circ - 5^\circ$ from the Galactic plane, indicating the presence of a possibly similar structure in the Galactic bulge.

preprint2015arXiv

Resonant Clumping and Substructure in Galactic Discs

We describe a method to extract resonant orbits from N-body simulations exploiting the fact that they close in a frame rotating with a constant pattern speed. Our method is applied to the N-body simulation of the Milky Way by Shen et al. (2010). This simulation hosts a massive bar, which drives strong resonances and persistent angular momentum exchange. Resonant orbits are found throughout the disc, both close to the bar itself and out to the very edges of the disc. Using Fourier spectrograms, we demonstrate that the bar is driving kinematic substructure even in the very outer parts of the disc. We identify two major orbit families in the outskirts of the disc that make significant contributions to the kinematic landscape, namely the m:l = 3:-2 and 1:-1 families resonating with the pattern speed of the bar. A mechanism is described that produces bimodal distributions of Galactocentric radial velocities at selected azimuths in the outer disc. It occurs as a result of the temporal coherence of particles on the 3:-2 resonant orbits, which causes them to arrive simultaneously at pericentre or apocentre. This resonant clumping, due to the in-phase motion of the particles through their epicycle, leads to both inward and outward moving groups which belong to the same orbital family and consequently produce bimodal radial velocity distributions. This is a possible explanation of the bimodal velocity distributions observed towards the Galactic anti-Centre by Liu et al. (2012). Another consequence is that transient overdensities appear and dissipate (in a symmetric fashion) on timescales equal to the their epicyclic period resulting in a periodic pulsing of the disc's surface density.

preprint2015arXiv

Resonant Orbits and the High Velocity Peaks Towards the Bulge

We extract the resonant orbits from an N-body bar that is a good representation of the Milky Way, using the method recently introduced by Molloy et al. (2015). By decomposing the bar into its constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards the Galactic centre. The resonant orbits are shown to have distinct kinematic features and may be used to explain the cold velocity peak seen in the APOGEE commissioning data (Nidever at al., 2012). We show that high velocity peaks are a natural consequence of the motions of stars in the 2:1 orbit family and that stars on other higher order resonances can contribute to the peaks. The locations of the peaks vary with bar angle and, with the tacit assumption that the observed peaks are due to the 2:1 family, we find that the locations of the high velocity peaks correspond to bar angles in the range 10 < theta_bar < 25 (deg). However, some important questions about the nature of the peaks remain, such as their apparent absence in other surveys of the Bulge and the deviations from symmetry between equivalent fields in the north and south. We show that the absence of a peak in surveys at higher latitudes is likely due to the combination of a less prominent peak and a lower number density of bar supporting orbits at these latitudes.

preprint2015arXiv

Theoretical Models of the Galactic Bulge

Near infrared images from the COBE satellite presented the first clear evidence that our Milky Way galaxy contains a boxy shaped bulge. Recent years have witnessed a gradual paradigm shift in the formation and evolution of the Galactic bulge. Bulges were commonly believed to form in the dynamical violence of galaxy mergers. However, it has become increasingly clear that the main body of the Milky Way bulge is not a classical bulge made by previous major mergers, instead it appears to be a bar seen somewhat end-on. The Milky Way bar can form naturally from a precursor disk and thicken vertically by the internal firehose/buckling instability, giving rise to the boxy appearance. This picture is supported by many lines of evidence, including the asymmetric parallelogram shape, the strong cylindrical rotation (i.e., nearly constant rotation regardless of the height above the disk plane), the existence of an intriguing X-shaped structure in the bulge, and perhaps the metallicity gradients. We review the major theoretical models and techniques to understand the Milky Way bulge. Despite the progresses in recent theoretical attempts, a complete bulge formation model that explains the full kinematics and metallicity distribution is still not fully understood. Upcoming large surveys are expected to shed new light on the formation history of the Galactic bulge.

preprint2014arXiv

Are high velocity peaks in the Milky Way bulge due to the bar?

Recently the commissioning APOGEE observations of the Galactic bulge reported that a significant fraction of stars ($\sim10%$) are in a cold ($σ_{\rm V} \approx 30$ km/s) high velocity peak (Galactocentric radial velocity $\approx 200$ km/s). These stars are speculated to reflect the stellar orbits in the Galactic bar. In this study, we use two $N$-body models of the Milky Way-like disk galaxy with different bar strengths to critically examine this possibility. The general trends of the Galactocentric radial velocity distribution in observations and simulations are similar, but neither our models nor the BRAVA data reveal a statistically significant cold high velocity peak. A Monte Carlo test further suggests that it is possible for a spurious high velocity peak to appear if there are only a limited number of stars observed. Thus, the reported cold high velocity peak, even if it is real, is unlikely due to stars on the bar-supporting orbits. Our models do predict an excess of stars with high radial velocity, but not in a distinct peak. In the distance--velocity diagram, the high velocity particles in different fields exist at a similar distance $\sim8.5 \pm 1$ kpc away from the Sun. This result may be explained with geometric intersections between the line-of-sight and the particle orbits; high velocity stars naturally exist approximately at the tangent point, without constituting a distinct peak. We further demonstrate that even without the presence of a bar structure, particle motions in an axisymmetric disk can also exhibit an excess of high velocity stars.

preprint2014arXiv

The X-Shaped Milky Way Bulge in OGLE-III Photometry and in N-Body Models

We model the split red clump of the Galactic bulge in OGLE-III photometry, and compare the results to predictions from two N-body models. Our analysis yields precise maps of the brightness of the two red clumps, the fraction of stars in the more distant peak, and their combined surface density. We compare the observations to predictions from two N-body models previously used in the literature. Both models correctly predict several features as long as one assumes an angle $α_{\rm{Bar}} \approx 30^{\circ}$ between the Galactic bar's major axis and the line of sight to the Galactic centre. In particular that the fraction of stars in the faint red clump should decrease with increasing longitude. The biggest discrepancies between models and data are in the rate of decline of the combined surface density of red clump stars toward negative longitudes and of the brightness difference between the two red clumps toward positive longitudes, with neither discrepancy exceeding $\sim$25% in amplitude. Our analysis of the red giant luminosity function also yields an estimate of the red giant branch bump parameters toward these high-latitude fields, and evidence for a high rate ($\sim$25%) of disk contamination in the bulge at the colour and magnitude of the red clump, with the disk contamination rate increasing toward sightlines further distant from the plane.

preprint2014arXiv

Uncertainties in the Deprojection of the Observed Bar Properties

In observations, it is important to deproject the two fundamental quantities characterizing a bar, i.e., its length ($a$) and ellipticity ($e$), to face-on values before any careful analyses. However, systematic estimation on the uncertainties of the commonly used deprojection methods is still lacking. Simulated galaxies are well suited in this study. We project two simulated barred galaxies onto a 2D plane with different bar orientations and disk inclination angles ($i$). Bar properties are measured and deprojected with the popular deprojection methods in the literature. Generally speaking, deprojection uncertainties increase with increasing $i$. All the deprojection methods behave badly when $i$ is larger than $60^\circ$, due to vertical thickness of the bar. Thus, future statistical studies of barred galaxies should exclude galaxies more inclined than $60^\circ$. At moderate inclination angles ($i\leq60^\circ$), 2D deprojection methods (analytical and image stretching) and Fourier-based methods (Fourier decomposition and bar-interbar contrast) perform reasonably well with uncertainties $\sim10\%$ in both the bar length and ellipticity. Whereas the uncertainties of the 1D analytical deprojection can be as high as $100\%$ in certain extreme case. We find that different bar measurement methods show systematic differences in the deprojection uncertainties. We further discuss the deprojection uncertainty factors with the emphasis on the most important one, i.e., the 3D structure of the bar itself. We construct two triaxial toy bar models that can qualitatively reproduce the results of the 1D and 2D analytical deprojections; they confirm that the vertical thickness of the bar is the main source of uncertainties.

preprint2013arXiv

A Schwarzschild model of the Galactic bar with initial density from N-body simulations

Using the potential from N-body simulations, we construct the Galactic bar models with the Schwarzschild method. By varying the pattern speed and the position angle of the bar, we find that the best-fit bar model has pattern speed $Ω_{\rm p}=40\ \rm{km\ s^{-1}\ kpc^{-1}}$, and bar angle $θ_{\rm bar}=45^{\circ}$. $N$-body simulations show that the best-fit model is stable for more than 1.5 Gyrs. Combined with the results in Wang et al. (2012), we find that the bar angle and/or the pattern speed are not well constrained by BRAVA data in our Schwarzschild models. The proper motions predicted from our model are slightly larger than those observed in four fields. In the future, more kinematic data from the ground and space-based observations will enable us to refine our model of the Milky Way bar.

preprint2013arXiv

On the Offset of Barred Galaxies From the Black Hole $M_{BH}$-$σ$ Relationship

We use collisionless $N$-body simulations to determine how the growth of a supermassive black hole (SMBH) influences the nuclear kinematics in both barred and unbarred galaxies. In the presence of a bar, the increase in the velocity dispersion $σ$ (within the effective radius) due to the growth of an SMBH is on average $\lesssim 10%$, whereas the increase is only $\lesssim 4%$ in an unbarred galaxy. In a barred galaxy, the increase results from a combination of three separate factors (a) orientation and inclination effects; (b) angular momentum transport by the bar that results in an increase in the central mass density; (c) an increase in the vertical and radial velocity anisotropy of stars in the vicinity of the SMBH. In contrast the growth of the SMBH in an unbarred galaxy causes the velocity distribution in the inner part of the nucleus to become less radially anisotropic. The increase in $σ$ following the growth of the SMBH is insensitive to a variation of a factor of 10 in the final mass of the SMBH, showing that it is the growth process rather than the actual SMBH mass that alters bar evolution in a way that increases $σ$. We argue that using an axisymmetric stellar dynamical modeling code to measure SMBH masses in barred galaxies could result in a slight overestimate of the derived $M_{BH}$, especially if a constant M/L ratio is assumed. We conclude that the growth of a black hole in the presence of a bar could result in an increase in $σ$ which is roughly of $4-8%$ larger than the increase that occurs in an axisymmetric system. While the increase in $σ$ due to SMBH growth in a barred galaxy might partially account for the claimed offset of barred galaxies and pseudo bulges from the $M_{BH}-σ$ relation obtained for elliptical galaxies and classical bulges in unbarred galaxies, it is inadequate to account for all of the offset.

preprint2013arXiv

The effect of bars on the M-sigma relation: offset, scatter and residuals correlations

We analyse a set of collisionless disc galaxy simulations to study the consequences of bar formation and evolution on the M-sigma relation of supermassive black holes. The redistribution of angular momentum driven by bars leads to a mass increase within the central region, raising the velocity dispersion of the bulge, sigma, on average by ~12% and as much as ~20%. If a disc galaxy with a SMBH satisfying the M-sigma relation forms a bar, and the SMBH does not grow in the process, then the increase in sigma moves the galaxy off the M-sigma relation. We explore various effects that can affect this result including contamination from the disc and anisotropy. The displacement from the M-sigma relation for individual model barred galaxies correlates with both M(B)/M(B+D) and beta_phi(B+D) measured within the effective radius of the bulge. Overall, this process leads to an M-sigma for barred galaxies offset from that of unbarred galaxies, as well as an increase in its scatter. We assemble samples of unbarred and barred galaxies with classical bulges and find tentative hints of an offset between the two consistent with the predicted. Including all barred galaxies, rather than just those with a classical bulge, leads to a significantly larger offset.

preprint2012arXiv

Made-to-measure galaxy models - III Modelling with Milky Way observations

We demonstrate how the Syer & Tremaine made-to-measure method of stellar dynamical modelling can be adapted to model a rotating galactic bar. We validate our made-to-measure changes using observations constructed from the existing Shen et al. (2010) N-body model of the Milky Way bar, together with kinematic observations of the Milky Way bulge and bar taken by the Bulge Radial Velocity Assay (BRAVA). Our results for a combined determination of the bar angle and bar pattern speed (~30 degrees and ~40 km/s/kpc) are consistent with those determined by the N-body model. Whilst the made-to-measure techniques we have developed are applied using a particular N-body model and observational data set, they are in fact general and could be applied to other Milky Way modelling scenarios utilising different N-body models and data sets. Additionally, we use the exercise as a vehicle for illustrating how N-body and made-to-measure methods might be combined into a more effective method.

preprint2012arXiv

The Vertical X-shaped Structure in the Milky Way: Evidence from a Simple Boxy Bulge Model

A vertical X-shaped structure was recently reported in the Galactic bulge. 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 degrees. 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.

preprint2011arXiv

Orbit-Based Dynamical Models of the Sombrero Galaxy (NGC 4594)

We present axisymmetric, orbit-based models to study the central black hole, stellar mass-to-light ratio, and dark matter halo of NGC 4594 (M104, the Sombrero Galaxy). For stellar kinematics, we use published high-resolution kinematics of the central region taken with the Hubble Space Telescope, newly obtained Gemini long-slit spectra of the major axis, and integral field kinematics from the SAURON instrument. At large radii, we use globular cluster kinematics to trace the mass profile and apply extra leverage to recovering the dark matter halo parameters. We find a black hole of mass M_{\bullet}=(6.6 +/- 0.4) x 10^8 M_{\odot}, and determine the stellar M/L_I=3.4 +/- 0.05 (uncertainties are the 68% confidence band marginalized over the other parameters). Our best fit dark matter halo is a cored logarithmic model with asymptotic circular speed V_c=376 +/- 12 km/s and core radius r_c= 4.7 +/- 0.6 kpc. The fraction of dark to total mass contained within the half-light radius is 0.52. Taking the bulge and disk components into account in our calculation of σ_e puts NGC 4594 squarely on the M-σrelation. We also determine that NGC 4594 lies directly on the M-L relation.

preprint2011arXiv

The Bulge Radial Velocity Assay (BRAVA): II. Complete Sample and Data Release

We present new radial velocity measurements from the Bulge Radial Velocity Assay (BRAVA), a large scale spectroscopic survey of M-type giants in the Galactic bulge/bar region. The sample of ~4500 new radial velocities, mostly in the region -10 deg < l < +10 deg and b ~ -6 deg more than doubles the existent published data set. Our new data extend our rotation curve and velocity dispersion profile to +20 deg, which is ~2.8 kpc from the Galactic Center. The new data confirm the cylindrical rotation observed at -6 deg and -8 deg, and are an excellent fit to the Shen et al. (2010) N-body bar model. We measure the strength of the TiO molecular band as a first step towards a metallicity ranking of the stellar sample, from which we confirm the presence of a vertical abundance gradient. Our survey finds no strong evidence of previously unknown kinematic streams. We also publish our complete catalog of radial velocities, photometry, TiO band strengths, and spectra, which is available at the IRSA archive: http://irsa.ipac.caltech.edu/ as well as at UCLA: http://brava.astro.ucla.edu/.

preprint2010arXiv

Long-lived double-barred galaxies in N-body simulations

Many barred galaxies harbor small-scale secondary bars in the center. The evolution of such double-barred galaxies is still not well understood, partly because of a lack of realistic N-body models with which to study them. Here we report the generation of such systems in the presence of rotating pseudobulges. We demonstrate with high mass and force resolution collisionless N-body simulations that long-lived secondary bars can form spontaneously without requiring gas, contrary to previous claims. We find that secondary bars rotate faster than primary ones. The rotation is not rigid: the secondary bars pulsate, with their amplitude and pattern speed oscillating as they rotate through the primary bars. This self-consistent study supports previous work based on orbital analysis in the potential of two rigidly rotating bars. We also characterize the density and kinematics of the N-body simulations of the double-barred galaxies, compare with observations to achieve a better understanding of such galaxies. The pulsating nature of secondary bars may have important implications for understanding the central region of double-barred galaxies.

preprint2010arXiv

Our Milky Way as a Pure-Disk Galaxy -- A Challenge for Galaxy Formation

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 ~ 4kpc and extends 20 degrees 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.

preprint2010arXiv

The Supermassive Black Hole and Dark Matter Halo of NGC 4649 (M60)

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 the black hole mass = 4.5 \pm 1.0 10^9 \Msun and M/L = 8.7 \pm 1.0 (or 8.0 \pm 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 black hole mass 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 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.

preprint2010arXiv

The VIRUS-P Exploration of Nearby Galaxies (VENGA): Survey Design and First Results

VENGA is a large-scale extragalactic IFU survey, which maps the bulges, bars and large parts of the outer disks of 32 nearby normal spiral galaxies. The targets are chosen to span a wide range in Hubble types, star formation activities, morphologies, and inclinations, at the same time of having vast available multi-wavelength coverage from the far-UV to the mid-IR, and available CO and 21cm mapping. The VENGA dataset will provide 2D maps of the SFR, stellar and gas kinematics, chemical abundances, ISM density and ionization states, dust extinction and stellar populations for these 32 galaxies. The uniqueness of the VIRUS-P large field of view permits these large-scale mappings to be performed. VENGA will allow us to correlate all these important quantities throughout the different environments present in galactic disks, allowing the conduction of a large number of studies in star formation, structure assembly, galactic feedback and ISM in galaxies.