We analyze the kinematics and spatial distribution of fundamental-mode 15,599 RR Lyrae (RRL) stars in the Galactic bulge by combining OGLE-IV photometric data and Gaia DR2 proper motions. We show that the longitudinal proper motion and the line-of-sight velocity give similar results for the rotation in the Galactic central regions. The angular velocity of bulge RRLs is found to be around $35$ km s$^{-1}$ kpc$^{-1}$, significantly smaller than that for the majority of bulge stars ($50-60$ km s$^{-1}$ kpc$^{-1}$); bulge RRLs have larger velocity dispersion (120$-$140 km s$^{-1}$) than younger stars. The dependence of the kinematics of the bulge RRLs on their metallicities is shown by their rotation curves and spatial distributions. Metal-poor RRLs ([Fe/H]<$-1$) show a smaller bar angle than metal-rich ones. We also find clues suggesting that RRLs in the bulge are not dominated by halo stars. These results might explain some previous conflicting results over bulge RRLs and help understand the chemodynamical evolution of the Galactic bulge.

GRB 200219A is a short gamma-ray burst (GRB) with an extended emission (EE) lasting $\sim 90$s. By analyzing data observed with the {\em Swift}/BAT and {\em Fermi}/GBM, we find that a cutoff power law (CPL) model can adequately well fit the spectra of the initial short pulse with $\rm E_{p}=1387^{+232}_{-134}$ keV. More interestingly, together with the EE component and early X-ray data, it exhibits a plateau emission smoothly connected with a $\sim t^{-1}$ segment, and followed by an extremely steep decay. The short GRB composed of those three segments is unique in the {\em Swift} era, and is very difficult to explain with the standard internal/external shock model of a black hole central engine, but could be consistent with the prediction of a magnetar central engine from the merger of an NS binary. We suggest that the plateau emission followed by a $\sim t^{-1}$ decay phase is powered by the spin-down of millisecond magnetar, which loses its rotation energy via GW quadrupole radiation. Then, the abrupt drop decay is caused by magnetar collapsing into black hole before switching to the EM-dominated emission. This is the first short GRB for which the X-ray emission has such an intriguing feature powered by magnetar via GW-dominated radiation. If this is the case, one can estimate the physical parameters of a magnetar, the GW signal powered by magnetar and the merger-nova emission are also discussed.

We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST, and Gaia, and we use these data to construct a full 7D set consisting of positions, space motions, and metallicity. We remove substructure from our samples using integrals of motion based on the method of Xue et al. We find radially dominated kinematic profiles with nearly constant anisotropy within 20 kpc, beyond which the anisotropy profile gently declines although remains radially dominated to the furthest extents of our sample. Independent of star type or substructure removal, the anisotropy depends on metallicity, such that the orbits of the stars become less radial with decreasing metallicity. For $-1.7<$ [Fe/H] $<-1$, the smooth, diffuse halo anisotropy profile begins to decline at Galactocentric distances $\sim20$ kpc, from $\beta\sim0.9$ to 0.7 for K giants and from $\beta\sim0.8$ to 0.1 for blue horizontal branch stars. For [Fe/H] $<-1.7$, the smooth, diffuse halo anisotropy remains constant along all distances with $0.2<\beta<0.7$ depending on the metallicity range probed, although independent on star type. These samples are ideal for estimating the total Galactic mass as they represent the virialized stellar halo system.

To understand the history and formation mechanisms of galaxies it is crucial to determine their current multidimensional structure. Here we focus on stellar population properties, such as metallicity and [$\alpha$/Fe] enhancement. We devise a new technique to recover the distribution of these parameters using spatially resolved, line-of-sight averaged data. Our chemodynamical method is based on the made-to-measure (M2M) framework and results in an $N$-body model for the abundance distribution. We test our method on a mock dataset and find that the radial and azimuthal profiles are well-recovered, however only the overall shape of the vertical profile matches the true profile. We apply our procedure to spatially resolved maps of mean [Z/H] and [$\alpha$/Fe] for the Andromeda Galaxy, using an earlier barred dynamical model of M31. We found that the metallicity is enhanced along the bar, with possible maxima at the ansae. In the edge-on view the [Z/H] distribution has an X shape due to the boxy/peanut bulge; the average vertical metallicity gradient is equal to $-0.132\pm0.006$ dex/kpc. We identify a metallicity-enhanced ring around the bar which also has relatively lower [$\alpha$/Fe]. The highest [$\alpha$/Fe] is found in the centre, due to the classical bulge. Away from the centre, the $\alpha$-overabundance in the bar region increases with height, which could be an indication of a thick disc. We argue that the galaxy assembly resulted in a sharp peak of metallicity in the central few hundred parsecs and a more gentle negative gradient in the remaining disc, but no [$\alpha$/Fe] gradient. The formation of the bar lead to the re-arrangement of the [Z/H] distribution, causing a flat gradient along the bar. Subsequent star formation close to the bar ends may have produced the metallicity enhancements at the ansae and the [Z/H] enhanced lower-$\alpha$ ring.

We examined the solar gravitational lens (SGL) as the means to produce direct high-resolution, multipixel images of exoplanets. The properties of the SGL are remarkable: it offers maximum light amplification of ~1e11 and angular resolution of ~1e-10 arcsec. A probe with a 1-m telescope in the SGL focal region can image an exoplanet at 30 pc with 10-kilometer resolution on its surface, sufficient to observe seasonal changes, oceans, continents, surface topography. We reached and exceeded all objectives set for our study: We developed a new wave-optical approach to study the imaging of exoplanets while treating them as extended, resolved, faint sources at large but finite distances. We properly accounted for the solar corona brightness. We developed deconvolution algorithms and demonstrated the feasibility of high-quality image reconstruction under realistic conditions. We have proven that multipixel imaging and spectroscopy of exoplanets with the SGL are feasible. We have developed a new mission concept that delivers an array of optical telescopes to the SGL focal region relying on three innovations: i) a new way to enable direct exoplanet imaging, ii) use of smallsats solar sails fast transit through the solar system and beyond, iii) an open architecture to take advantage of swarm technology. This approach enables entirely new missions, providing a great leap in capabilities for NASA and the greater aerospace community. Our results are encouraging as they lead to a realistic design for a mission that will be able to make direct resolved images of exoplanets in our stellar neighborhood. It could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies. The architecture and mission concepts for a mission to the strong interference region of the SGL are promising and should be explored further.

Short GRBs with extended emission (EE) that are composed initially of a short-hard spike and followed by a long-lasting EE, are thought to be classified as a subsection of short GRBs. The narrow energy band available during the {\em Swift} era combined with a lack of spectral information prevented discovery of the intrinsic properties of those events. In this paper, we performed a systematic search of short GRBs with EE by using all available {\em Fermi}/GBM data. The search identified 26 GBM-detected short GRBs with EE that are similar to GRB 060614 observed by {\em Swift}/BAT. We focus on investigating the spectral and temporal properties for both the hard spike and the EE components of all 26 GRBs, and explore differences and possible correlations between them. We find that while the peak energy ($E_{\rm p}$) of the hard spikes is a little bit harder than that of the EE, but their fluences are comparable. The harder $E_{\rm p}$ seems to correspond to a larger fluence and peak flux with a large scatter for both the hard spike and EE components. Moreover, the $E_{\rm p}$ of both the hard spikes and EE are compared to other short GRBs. Finally, we also compare the properties of GRB 170817A with those short GRBs with EE and find no significant statistical differences between them. We find that GRB 170817A has the lowest $E_{\rm p}$, likely because it was off-axis.

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_{\phi}$ phase space near the Solar circle. Interestingly, the phase space snail shell exists only in the arches with $|V_{\phi} - 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_{\phi}$), 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_{\phi}$ color-coded phase space is also discussed.

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.

In the present study, we developed algorithms that are capable of measuring apparent magnitudes and the effective temperature of stars using raw images shot with amateur telescopes. The regularized Radial Basis Function (RBF) network, one of the machine learning algorithms, was employed to measure the effective temperature, and the simple function fitting method was adopted to measure the apparent magnitude. The achieved results are satisfying. After the white balance and noise cancellation process was simply calibrated, it was demonstrated that the measurements of the effective temperature had mean fraction errors at around 9%, and the measurements of the magnitudes had absolute error at nearly 0.1.

Tremaine and Weinberg (TW) proposed a conceptually simple procedure relying on long-slit spectroscopy to measure the pattern speeds of bars ($\Omega_{\rm p}$) in disk galaxies. Using a simulated galaxy, we investigate the potential biases and uncertainties of TW measurements using increasingly popular integral-field spectrographs (IFSs), for which multiple pseudo-slits (and thus independent measurements) can be constructed with a single observation. Most important, to establish the spatial coverage required and ensure the validity of the measurements, the inferred $\Omega_{\rm p}$ must asymptotically converge as the (half-)length of each pseudo-slit used is increased. The requirement for our simulation is to reach $\approx1.3$ times the half-light radius, but this may vary from galaxy to galaxy. Only those slits located within the bar region yield accurate measurements. We confirm that the position angle of the disk is the dominant source of systematic error in TW $\Omega_{\rm p}$ measurements, leading to under/overestimates of tens of percents for inaccuracies of even a few degrees. Recasting the data so that the data grid aligns with the disk major axis leads to slightly reduced uncertainties. Accurate measurements are obtained only for well-defined ranges of the bar angle (relative to the galaxy major axis) $\phi_{\rm bar}$ and the inclination angle $i$, here $10\lesssim\phi_{\rm bar}\lesssim75^{\circ}$ and $105\lesssim\phi_{\rm bar}\lesssim170^{\circ}$ and $15\lesssim i\lesssim70^{\circ}$. The adopted (pseudo-)slit widths, spatial resolution, and (unless extremely aggressive) spatial binning of IFS data have no significant impact on the measurements. Our results thus provide useful guidelines for reliable and accurate direct $\Omega_{\rm p}$ measurements with IFS observations.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

The Zwicky Transit Factory (ZTF) is a powerful time domain survey facility with a large field of view. We apply the synthetic tracking technique to integrate a ZTF's long-dwell dataset, which consists of 133 nominal 30-second exposure frames spanning about 1.5 hours, to search for slowly moving asteroids down to approximately 23rd magnitude. We found more than one thousand objects from searching 40 CCD-quadrant subfields, each of which covers a field size of $\sim$0.73 deg$^2$. While most of the objects are main belt asteroids, there are asteroids belonging to families of Trojan, Hilda, Hungaria, Phocaea, and near-Earth-asteroids. Such an approach is effective and productive. Here we report the data process and results.

The inner disc of the local group galaxy M33 appears to be in settled rotational balance, and near IR images reveal a mild, large-scale, two-arm spiral pattern with no strong bar. We have constructed N-body models that match all the extensive observational data on the kinematics and surface density of stars and gas in the inner part of M33. We find that currently favoured models are unstable to the formation of a strong bar of semi-major axis 2 < a_B < 3 kpc within 1 Gyr, which changes the dynamical properties of the models to become inconsistent with the current, apparently well-settled, state. The formation of a bar is unaffected by how the gas component is modelled, by increasing the mass of the nuclear star cluster, or by making the dark matter halo counter-rotate, but it can be prevented by either reducing the mass-to-light ratio of the stars to Upsilon_V ~ 0.6 or Upsilon_K ~ 0.23 in solar units or by increasing the random motions of the stars. Also a shorter and weaker bar results when the halo is rigid and unresponsive. However, all three near-stable models support multi-arm spirals, and not the observed large-scale bi-symmetric spiral. A two-arm spiral pattern could perhaps be tidally induced, but such a model would require an implausibly low mass disc to avoid a bar and there is no visible culprit. Thus the survival of the current state of this exceptionally well-studied galaxy is not yet understood. We also suspect that many other unbarred galaxies present a similar puzzle.

The event GW170817/GRB 170817A, discovered via the successful joint observation of its gravitational wave radiation and its multi-wavelength electromagnetic counterparts, was the first definite "smoking-gun" from the merger of two neutron stars (NSs). However, the remnant of the merger remains unknown. Piro et al. recently claimed that a low-significance X-ray variability in GRB 170817A. By systematically comparing the properties of variability in the afterglow of GRB 170817A and X-ray flares in GRB afterglows, we find that this X-ray variability seems to share similar statistical correlations with X-ray flares in GRB afterglows. We further investigate several possible merger product scenarios to see whether they can produce the observed X-ray variability in GRB 170817A. The first scenario invokes a stable magnetar as the central engine producing the later X-ray variability via differential rotation or fall-back accretion onto the NS. The second scenario invokes a black hole as the central engine with a fall-back accretion process. The final scenario is a central engine with a long-lived supra-massive NS. We find that the first two scenarios have difficulty producing the later X-ray variability, which requires either an impractical NS magnetic field or an extraordinarily large stellar envelope and an extremely long accretion timescale. However, the third scenario seems to be consistent with observations, and the later X-ray variability can be produced by the magnetosphere which is expelled following the collapse of the NS with a $B_p\in(3.6, 13.5)\times10^{13}$ G.

The accelerating expansion of universe can be described by the non-zero cosmological constant or the dark energy. However, the origin of the dark energy remains a mystery of modern physics. The local Lorentz invariance is the most exact symmetry of the Nature on the one hand, but all quantum gravity theories predict Lorentz violation on the other hand. The local Lorentz violation induced by the quantum gravity at the very early universe may be transformed into large scale by inflation. Combining the low-$l$ anomalies of the CMB spectrum, we propose that the local Lorentz invariance may be broken at the large scale. We construct the effective gravity at the cosmic scale with a local $SO(3)$ symmetry. The theory exhibits non-trivial contortion distribution even with scalar matter source. The FRW like solution of the theory is analyzed and the contortion distribution contributes a dark energy like effect which is responsible for the accelerating expansion of the universe. It reveals that the dark energy may be the remnants of quantum gravity in this sense.

The favored progenitor model for Gamma-ray Bursts (GRBs) with Supernova (SN) association is the core collapse of massive stars. One possible outcome of such a collapse is a rapidly spinning, strongly magnetized neutron star ("magnetar"). We systematically analyze the multi-wavelength data of GRB/SN associations detected by several instruments before 2017 June. Twenty GRB/SN systems have been confirmed via direct spectroscopic evidence or a clear light curve bump, as well as some spectroscopic evidence resembling a GRB-SN. We derive/collect the basic physical parameters of the GRBs and the SNe, and look for correlations among these parameters. We find that the peak brightness, $\rm ^{56}Ni$ mass, and explosion energy of SNe associated with GRBs are statistically higher than other Type Ib/c SNe. A statistically significant relation between the peak energy of GRBs and the peak brightness of their associated SNe is confirmed. No significant correlations are found between the GRB energies (either isotropic or beaming-corrected) and the supernova energy. We investigate the energy partition within these systems and find that the beaming-corrected GRB energy of most systems is smaller than the SN energy, with less than 30\% of the total energy distributed in the relativistic jet. The total energy of the systems is typically smaller than the maximum available energy of a millisecond magnetar ($2\times 10^{52}$ erg), especially if aspherical SN explosions are considered. The data are consistent with-though not a proof of-the hypothesis that most, but not all, GRB/SN systems are powered by millisecond magnetars.

The anisotropy parameter $\beta$ 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 $\beta$ 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 $(\sigma_r, \sigma_\theta, \sigma_\phi)$ (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 $\beta \approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $\beta \approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $\beta\approx0.8$ and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to $\beta<0.3$ past 100 kpc. Within 8 kpc, $\beta<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.

The Solar Gravitational Lens (SGL) allows for major brightness amplification ($\sim 10^{11}$ at wavelength of $1~\mu$m) and extreme angular resolution ($\sim10^{-10}$ arcsec) within a narrow field of view. A meter-class telescope, with a modest coronagraph to block solar light with 1e-6 suppression placed in the focal area of the SGL, can image an exoplanet at a distance of 30 parsec with few kilometer-scale resolution on its surface. Notably, spectroscopic broadband SNR is $\sim 10^{-6}$ in two weeks of integration time, providing this instrument with incredible remote sensing capabilities. A mission capable of exploiting the remarkable optical properties of the SGL allows for direct high-resolution imaging/spectroscopy of a potentially habitable exoplanet. Such missions could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies.

The remarkable optical properties of the solar gravitational lens (SGL) include major brightness amplification (~1e11 at wavelength of 1 um) and extreme angular resolution (~1e-10 arcsec) in a narrow field of view. A mission to the SGL carrying a modest telescope and coronagraph opens up a possibility for direct megapixel imaging and high-resolution spectroscopy of a habitable Earth-like exoplanet at a distance of up to 100 light years. The entire image of such a planet is compressed by the SGL into a region with a diameter of ~1.3 km in the vicinity of the focal line. The telescope, acting as a single pixel detector while traversing this region, can build an image of the exoplanet with kilometer-scale resolution of its surface, enough to see its surface features and signs of habitability. We report here on the results of our initial study of a mission to the deep outer regions of our solar system, with the primary mission objective of conducting direct megapixel high-resolution imaging and spectroscopy of a potentially habitable exoplanet by exploiting the remarkable optical properties of the SGL. Our main goal was to investigate what it takes to operate spacecraft at such enormous distances with the needed precision. Specifically, we studied i) how a space mission to the focal region of the SGL may be used to obtain high-resolution direct imaging and spectroscopy of an exoplanet by detecting, tracking, and studying the Einstein ring around the Sun, and ii) how such information could be used to detect signs of life on another planet. Our results indicate that a mission to the SGL with an objective of direct imaging and spectroscopy of a distant exoplanet is challenging, but possible. We composed a list of recommendations on the mission architectures with risk and return tradeoffs and discuss an enabling technology development program.

The shape of LOSVDs (line-of-sight velocity distributions) carries important information about the internal dynamics of galaxies. The skewness of LOSVDs represents their asymmetric deviation from a Gaussian profile. Correlations between the skewness parameter ($h_3$) and the mean velocity (\vm) of a Gauss-Hermite series reflect the underlying stellar orbital configurations of different morphological components. Using two self-consistent $N$-body simulations of disk galaxies with different bar strengths, we investigate \hv\ correlations at different inclination angles. Similar to previous studies, we find anticorrelations in the disk area, and positive correlations in the bar area when viewed edge-on. However, at intermediate inclinations, the outer parts of bars exhibit anticorrelations, while the core areas dominated by the boxy/peanut-shaped (B/PS) bulges still maintain weak positive correlations. When viewed edge-on, particles in the foreground/background disk (the wing region) in the bar area constitute the main velocity peak, whereas the particles in the bar contribute to the high-velocity tail, generating the \hv\ correlation. If we remove the wing particles, the LOSVDs of the particles in the outer part of the bar only exhibit a low-velocity tail, resulting in a negative \hv\ correlation, whereas the core areas in the central region still show weakly positive correlations. We discuss implications for IFU observations on bars, and show that the variation of the \hv\ correlation in the disk galaxy may be used as a kinematic indicator of the bar and the B/PS bulge.

We present high-quality fluid dynamical simulations of isothermal gas flows in a rotating barred potential. We show that a large quantity of gas is driven right into the nucleus of a model galaxy when the potential lacks a central mass concentration, but the inflow stalls at a nuclear ring in comparison simulations that include a central massive object. The radius of the nuclear gas ring increases linearly with the mass of the central object. We argue that bars drive gas right into the nucleus in the early stages of disk galaxy formation, where a nuclear star cluster and perhaps a massive black hole could be created. The process is self-limiting, however, because inflow stalls at a nuclear ring once the mass of gas and stars in the nucleus exceeds ~1% of the disk mass, which shuts off rapid growth of the black hole. We briefly discuss the relevance of these results to the seeding of massive black holes in galaxies, the merger model for quasar evolution, and the existence of massive black holes in disk galaxies that lack a significant classical bulge.

A cold high-velocity (HV, $\sim$ 200 km/s) peak was first reported in several Galactic bulge fields based on the APOGEE commissioning observations. Both the existence and the nature of the high-velocity peak are still under debate. Here we revisit this feature with the latest APOGEE DR13 data. We find that most of the low latitude bulge fields display a skewed Gaussian distribution with a HV shoulder. However, only 3 out of 53 fields show distinct high-velocity peaks around 200 km/s. The velocity distribution can be well described by Gauss-Hermite polynomials, except the three fields showing clear HV peaks. We find that the correlation between the skewness parameter ($h_{3}$) and the mean velocity ($\bar{v}$), instead of a distinctive HV peak, is a strong indicator of the bar. It was recently suggested that the HV peak is composed of preferentially young stars. We choose three fields showing clear HV peaks to test this hypothesis using the metallicity, [$\alpha$/M] and [C/N] as age proxies. We find that both young and old stars show HV features. The similarity between the chemical abundances of stars in the HV peaks and the main component indicates that they are not systematically different in terms of chemical abundance or age. In contrast, there are clear differences in chemical space between stars in the Sagittarius dwarf and the bulge stars. The strong HV peaks off-plane are still to be explained properly, and could be different in nature.

Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.

The growth of black holes (BHs) in disk galaxies lacking classical bulges, which implies an absence of significant mergers, appears to be driven by secular processes. Short bars of sub-kiloparsec radius have been hypothesized to be an important mechanism for driving gas inflows to small scale, feeding central BHs. In order to quantify the maximum BH mass allowed by this mechanism, we examine the robustness of short bars to the dynamical influence of BHs. Large-scale bars are expected to be robust, long-lived structures; extremely massive BHs, which are rare, are needed to completely destroy such bars. However, we find that short bars, which are generally embedded in large-scale outer bars, can be destroyed quickly when BHs of mass $M_{\rm bh}\sim0.05-0.2\%$ of the total stellar mass ($M_\star$) are present. In agreement with this prediction, all galaxies observed to host short bars have BHs with a mass fraction less than $0.2\%M_\star$. Thus the dissolution of short inner bars is possible, perhaps even frequent, in the universe. An important implication of this result is that inner bar-driven gas inflows may be terminated when BHs grow to $\sim0.1\%M_\star$. We predict that $0.2\%M_\star$ is the maximum mass of BHs allowed if they are fed predominately via inner bars. This value matches well the maximum ratio of BH-to-host galaxy stellar mass observed in galaxies with pseudo-bulges and most nearby AGN host galaxies. This hypothesis provides a novel explanation for the lower $M_{\rm bh}/M_\star$ in galaxies which have avoided significant mergers compared with galaxies with classical bulges.

We study the structure of the inner Milky Way using the latest data release of the Vista Variables in Via Lactea (VVV) survey. The VVV is a deep near-infrared, multi-colour photometric survey with a coverage of 300 square degrees towards the Bulge/Bar. We use Red Clump (RC) stars to produce a high-resolution dust map of the VVV's field of view. From de-reddened colour-magnitude diagrams we select Red Giant Branch stars to investigate their 3D density distribution within the central 4 kpc. We demonstrate that our best-fit parametric model of the Bulge density provides a good description of the VVV data, with a median percentage residual of 5$\%$ over the fitted region. The strongest of the otherwise low-level residuals are overdensities associated with a low-latitude structure as well as the so-called X-shape previously identified using the split RC. These additional components contribute only $\sim5\%$ and $\sim7\%$ respectively to the Bulge mass budget. The best-fit Bulge is `boxy' with an axis ratio of [1:0.44:0.31] and is rotated with respect to the Sun-Galactic Centre line by at least $20^{\circ}$. We provide an estimate of the total, full sky, mass of the Bulge of $M_\mathrm{Bulge}^{\mathrm{Chabrier}} = 2.36 \times 10^{10} M_{\odot}$ for a Chabrier initial mass function. We show there exists a strong degeneracy between the viewing angle and the dispersion of the Red Clump absolute magnitude distribution. The value of the latter is strongly dependent on the assumptions made about the intrinsic luminosity function of the Bulge.

The Milky Way bulge shows a box/peanut or X-shaped bulge (hereafter BP/X) when viewed in infrared or microwave bands. We examine orbits in an N-body model of a barred disk galaxy that is scaled to match the kinematics of the Milky Way (MW) bulge. We generate maps of projected stellar surface density, unsharp masked images, 3D excess-mass distributions (showing mass outside ellipsoids), line-of-sight number count distributions, and 2D line-of-sight kinematics for the simulation as well as co-added orbit families, in order to identify the orbits primarily responsible for the BP/X shape. We estimate that between 19-23\% of the mass of the bar is associated with the BP/X shape and that most bar orbits contribute to this shape which is clearly seen in projected surface density maps and 3D excess mass for non-resonant box orbits, "banana" orbits, "fish/pretzel" orbits and "brezel" orbits. {We find that nearly all bar orbit families contribute some mass to the 3D BP/X-shape. All co-added orbit families show a bifurcation in stellar number count distribution with heliocentric distance that resembles the bifurcation observed in red clump stars in the MW. However, only the box orbit family shows an increasing separation of peaks with increasing galactic latitude $|b|$, similar to that observed.} Our analysis shows that no single orbit family fully explains all the observed features associated with the MW's BP/X shaped bulge, but collectively the non-resonant boxes and various resonant boxlet orbits contribute at different distances from the center to produce this feature. We propose that since box orbits have three incommensurable orbital fundamental frequencies, their 3-dimensional shapes are highly flexible and, like Lissajous figures, this family of orbits is most easily able to adapt to evolution in the shape of the underlying potential.

In barred galaxies, the contours of stellar velocity dispersions ($\sigma$) are generally expected to be oval and aligned with the orientation of bars. However, many double-barred (S2B) galaxies exhibit distinct $\sigma$ peaks on the minor axis of inner bars, which we termed "$\sigma$-humps," while two local $\sigma$ minima are present close to the ends of inner bars, i.e., "$\sigma$-hollows." Analysis of numerical simulations shows that $\sigma_z$-humps or hollows should play an important role in generating the observed $\sigma$-humps+hollows in low-inclination galaxies. In order to systematically investigate the properties of $\sigma_z$ in barred galaxies, we apply the vertical Jeans equation to a group of well-designed three-dimensional bar+disk(+bulge) models. A vertically thin bar can lower $\sigma_z$ along the bar and enhance it perpendicular to the bar, thus generating $\sigma_z$-humps+hollows. Such a result suggests that $\sigma_z$-humps+hollows can be generated by the purely dynamical response of stars in the presence of a, sufficiently massive, vertically thin bar, even without an outer bar. Using self-consistent $N$-body simulations, we verify the existence of vertically thin bars in the nuclear-barred and S2B models which generate prominent $\sigma$-humps+hollows. Thus the ubiquitous presence of $\sigma$-humps+hollows in S2Bs implies that inner bars are vertically thin. The addition of a bulge makes the $\sigma_z$-humps more ambiguous and thus tends to somewhat hide the $\sigma_z$-humps+hollows. We show that $\sigma_z$ may be used as a kinematic diagnostic of stellar components that have different thickness, providing a direct perspective on the morphology and thickness of nearly face-on bars and bulges with integral field unit spectroscopy.

Using the LAMOST-Gaia common stars, we demonstrate that the in-plane velocity field for the nearby young stars are significantly different from that for the old ones. For the young stars, the probably perturbed velocities similar to the old population are mostly removed from the velocity maps in the $X$--$Y$ plane. The residual velocity field shows that the young stars consistently move along $Y$ with faster $v_\phi$ at the trailing side of the local arm, while at the leading side, they move slower in azimuth direction. At both sides, the young stars averagely move inward with $v_R$ of $-5\sim-3$ km s$^{-1}$. The divergence of the velocity in $Y$ direction implies that the young stars are associated with a density wave nearby the local arm. We therefore suggest that the young stars may reflect the formation of the local spiral arm by correlating themselves with a density wave. The range of the age for the young stars is around 2 Gyr, which is sensible since the transient spiral arm can sustain that long. We also point out that alternative explanations of the peculiar velocity field for the young population cannot be ruled out only from this observed data.

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 $\sigma_\mathrm{LOS}$ along the minor axis of inner bars, which we term "$\sigma$-humps", that are often accompanied by ring/spiral-like features of increased $\sigma_\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 $\sigma$-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 $\sigma_z$, which exhibits peaks on the minor axis, i.e., humps. These $\sigma_z$ humps are responsible for producing the $\sigma$-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 "$\sigma$-hollows" previously identified. We also discuss the kinematic effect of a nuclear disk in S2Bs.

We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy-peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal rich stars also exhibit peaky features in their line-of-sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy-peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The 'banana' or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data.

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.

We present 947 radial velocities of RR Lyrae variable stars in four fields located toward the Galactic bulge, observed within the data from the ongoing Bulge RR Lyrae Radial Velocity Assay (BRAVA-RR). We show that these RR Lyrae stars exhibit hot kinematics and null or negligible rotation and are therefore members of a separate population from the bar/pseudobulge that currently dominates the mass and luminosity of the inner Galaxy. Our RR Lyrae stars predate these structures, and have metallicities, kinematics, and spatial distribution that are consistent with a "classical" bulge, although we cannot yet completely rule out the possibility that they are the metal-poor tail of a more metal rich ([Fe/H] ~ -1 dex) halo-bulge population. The complete catalog of radial velocities for the BRAVA-RR stars is also published electronically.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 $\alpha_{\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.

We investigate the matter power spectra in the power law and exponential types of viable $f(R)$ theories along with massive neutrinos. The enhancement of the matter power spectrum is found to be a generic feature in these models. In particular, we show that in the former type, such as the Starobinsky model, the spectrum is magnified much larger than the latter one, such as the exponential model. A greater scale of the total neutrino mass, $\Sigma m_{\nu}$, is allowed in the viable $f(R)$ models than that in the $\Lambda$CDM one. We obtain the constraints on the neutrino masses by using the CosmoMC package with the modified MGCAMB. Explicitly, we get $\Sigma m_{\nu} < 0.451 \ (0.214)\ \mathrm{eV}$ at 95% C.L. in the Starobinsky (exponential) model, while the corresponding one for the $\Lambda$CDM model is $\Sigma m_{\nu} < 0.200\ \mathrm{eV}$. Furthermore, by treating the effective number of neutrino species $N_{\mathrm{eff}}$ as a free parameter along with $\Sigma m_{\nu}$, we find that $N_{\mathrm{eff}} = 3.78^{+0.64}_{-0.84} (3.47^{+0.74}_{-0.60})$ and $\Sigma m_{\nu} = 0.533^{+0.254}_{-0.411}$ ($< 0.386) \ \mathrm{eV}$ at 95% C.L. in the Starobinsky (exponential) model.

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.

Recently the commissioning APOGEE observations of the Galactic bulge reported that a significant fraction of stars ($\sim10%$) are in a cold ($\sigma_{\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.

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 $\sigma$ (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 $\sigma$ 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 $\sigma$. 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 $\sigma$ which is roughly of $4-8%$ larger than the increase that occurs in an axisymmetric system. While the increase in $\sigma$ 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}-\sigma$ relation obtained for elliptical galaxies and classical bulges in unbarred galaxies, it is inadequate to account for all of the offset.

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.

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 $\Omega_{\rm p}=40\ \rm{km\ s^{-1}\ kpc^{-1}}$, and bar angle $\theta_{\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.

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.

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.

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/.

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 \sigma_e puts NGC 4594 squarely on the M-\sigma relation. We also determine that NGC 4594 lies directly on the M-L relation.

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.

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.

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.

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.

Galactic bars are the most important driver of secular evolution in galaxies. They can efficiently drive gas into the central kiloparsec of galaxies, thus feed circumnuclear starbursts, and possibly help to fuel AGN. The connection between bars and AGN activities has been actively debated in the past two decades. Previous work used fairly small samples and often lacked a proper control sample. They reported conflicting results on the correlation between bars and AGN activity. Here we revisit the bar-AGN and bar-starburst connections using the analysis of bars in a large sample of about 2000 SDSS disk galaxies (Barazza, Jogee, & Marinova 2008). We find that AGN and star-forming galaxies have similar optical bar fractions, 47% and 50%, respectively. Both bar fractions are higher than that in inactive galaxies (29%). We discuss the implications of the study on the relationship between host galaxies and their central activities.

This Letter is to investigate the physics of a newly discovered phenomenon -- contracting flare loops in the early phase of solar flares. In classical flare models, which were constructed based on the phenomenon of expansion of flare loops, an energy releasing site is put above flare loops. These models can predict that there is a vertical temperature gradient in the top of flare loops due to heat conduction and cooling effects. Therefore, the centroid of an X-ray looptop source at higher energy bands will be higher in altitude, for which we can define as normal temperature distribution. With observations made by {\it RHESSI}, we analyzed 10 M- or X-class flares (9 limb flares). For all these flares, the movement of looptop sources shows an obvious U-shaped trajectory, which we take as the signature of contraction-to-expansion of flare loops. We find that, for all these flares, normal temperature distribution does exist, but only along the path of expansion. The temperature distribution along the path of contraction is abnormal, showing no spatial order at all. The result suggests that magnetic reconnection processes in the contraction and expansion phases of these solar flares are different.

We have measured the stellar velocity dispersions (\sigma_*) and estimated the central black hole (BH) masses for over 900 broad-line active galactic nuclei (AGNs) observed with the Sloan Digital Sky Survey. The sample includes objects which have redshifts up to z=0.452, high quality spectra, and host galaxy spectra dominated by an early-type (bulge) component. The AGN and host galaxy spectral components were decomposed using an eigenspectrum technique. The BH masses (M_BH) were estimated from the AGN broad-line widths, and the velocity dispersions were measured from the stellar absorption spectra of the host galaxies. The range of black hole masses covered by the sample is approximately 10^6 < M_BH < 10^9 M_Sun. The host galaxy luminosity-velocity dispersion relationship follows the well-known Faber-Jackson relation for early-type galaxies, with a power-law slope 4.33+-0.21. The estimated BH masses are correlated with both the host luminosities (L_{H}) and the stellar velocity dispersions (\sigma_*), similar to the relationships found for low-redshift, bulge-dominated galaxies. The intrinsic scatter in the correlations are large (~0.4 dex), but the very large sample size allows tight constraints to be placed on the mean relationships: M_BH ~ L_H^{0.73+-0.05} and M_BH ~ \sigma_*^{3.34+-0.24}. The amplitude of the M_BH-\sigma_* relation depends on the estimated Eddington ratio, such that objects with larger Eddington ratios have smaller black hole masses than expected at a given velocity dispersion.

At the intersection of galactic dynamics, evolution and global structure, issues such as the relation between bars and spirals and the persistence of spiral patterns can be addressed through the characterization of the angular speeds of the patterns and their possible radial variation. The Radial Tremaine-Weinberg (TWR) Method, a generalized version of the Tremaine-Weinberg method for observationally determining a single, constant pattern speed, allows the pattern speed to vary arbitrarily with radius. Here, we perform tests of the TWR method with regularization on several simulated galaxy data sets. The regularization is employed as a means of smoothing intrinsically noisy solutions, as well as for testing model solutions of different radial dependence (e.g. constant, linear or quadratic). We test these facilities in studies of individual simulations, and demonstrate successful measurement of both bar and spiral pattern speeds in a single disk, secondary bar pattern speeds, and spiral winding (in the first application of a TW calculation to a spiral simulation). We also explore the major sources of error in the calculation and find uncertainty in the major axis position angle most dominant. In all cases, the method is able to extract pattern speed solutions where discernible patterns exist to within 20% of the known values, suggesting that the TWR method should be a valuable tool in the area of galactic dynamics. For utility, we also discuss the caveats in, and compile a prescription for, applications to real galaxies.

Although at least one quarter of early-type barred galaxies host secondary stellar bars embedded in their large-scale primary counterparts, the dynamics of such double barred galaxies are still not well understood. Recently we reported success at simulating such systems in a repeatable way in collisionless systems. In order to further our understanding of double-barred galaxies, here we characterize the density and kinematics of the N-body simulations of these galaxies. This will facilitate comparison with observations and lead to a better understanding of the observed double-barred galaxies. We find the shape and size of our simulated secondary bars are quite reasonable compared to the observed ones. We demonstrate that an authentic decoupled secondary bar may produce only a weak twist of the kinematic minor axis in the stellar velocity field, due to the relatively large random motion of stars in the central region. We also find that the edge-on nuclear bars are probably not related to boxy peanut-shaped bulges which are most likely to be edge-on primary large-scale bars. Finally we demonstrate that the non-rigid rotation of the secondary bar causes its pattern speed not to be derived with great accuracy using the Tremaine-Weinberg method. We also compare with observations of NGC 2950, a prototypical double-barred early-type galaxy, which suggest that the nuclear bar may be rotating in the opposite sense as the primary.

A large fraction of barred galaxies host secondary bars that are embedded in their large-scale primary counterparts. These are common also in gas poor early-type barred galaxies. 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 a new mechanism for generating such systems, namely the presence of rotating pseudo-bulges. We demonstate 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, however, 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. The pulsating nature of secondary bars may have important implications for understanding the central region of double-barred galaxies.

We have studied the wave dynamics and the Hawking radiation for the scalar field as well as the brane-localized gravitational field in the background of the braneworld black hole with tidal charge containing information of the extra dimension. Comparing with the four-dimensional black holes, we have observed the signature of the tidal charge which presents the signals of the extra dimension both in the wave dynamics and the Hawking radiation.

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 onto 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 of the warp inside $R_{26.5} ~ 4.5 \Rd$ is straight, while that beyond $R_{26.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 line of nodes of most warps forms a leading spiral might imply that the disc mass is significant in the centre. If the line of nodes 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 Gyrs, by which time another infall event can be expected.

We discuss an analytic approach for modeling structure formation in sheets, filaments and knots. This is accomplished by combining models of triaxial collapse with the excursion set approach: sheets are defined as objects which have collapsed along only one axis, filaments have collapsed along two axes, and halos are objects in which triaxial collapse is complete. In the simplest version of this approach, which we develop here, large scale structure shows a clear hierarchy of morphologies: the mass in large-scale sheets is partitioned up among lower mass filaments, which themselves are made-up of still lower mass halos. Our approach provides analytic estimates of the mass fraction in sheets, filaments and halos, and its evolution, for any background cosmological model and any initial fluctuation spectrum. In the currently popular $\Lambda$CDM model, our analysis suggests that more than 99% of the cosmic mass is in sheets, and 72% in filaments, with mass larger than $10^{10} M_{\odot}$ at the present time. For halos, this number is only 46%. Our approach also provides analytic estimates of how halo abundances at any given time correlate with the morphology of the surrounding large-scale structure, and how halo evolution correlates with the morphology of large scale structure.

Using an eigenspectrum decomposition technique, we separate the host galaxy from the broad line active galactic nucleus (AGN) in a set of 4666 spectra from the Sloan Digital Sky Survey (SDSS), from redshifts near zero up to about 0.75. The decomposition technique uses separate sets of galaxy and quasar eigenspectra to efficiently and reliably separate the AGN and host spectroscopic components. The technique accurately reproduces the host galaxy spectrum, its contributing fraction, and its classification. We show how the accuracy of the decomposition depends upon S/N, host galaxy fraction, and the galaxy class. Based on the eigencoefficients, the sample of SDSS broad-line AGN host galaxies spans a wide range of spectral types, but the distribution differs significantly from inactive galaxies. In particular, post-starburst activity appears to be much more common among AGN host galaxies. The luminosities of the hosts are much higher than expected for normal early-type galaxies, and their colors become increasingly bluer than early-type galaxies with increasing host luminosity. Most of the AGNs with detected hosts are emitting at between 1% and 10% of their estimated Eddington luminosities, but the sensitivity of the technique usually does not extend to the Eddington limit. There are mild correlations among the AGN and host galaxy eigencoefficients, possibly indicating a link between recent star formation and the onset of AGN activity. The catalog of spectral reconstruction parameters is available as an electronic table.

Employing the Quark Mass Denisity- and temperature- dependent model and the Hartle's method, We have studied the slowly rotating strange star with uniform angular velocity. The mass-radius relation, the moment of inertia and the frame dragging for different frequencies are given. We found that we cannot use the strange star to solve the challenges of Stella and Vietri for the horizontal branch oscillations and the moment of inertia $I_{45}/(M/M_ s)>2.3$. Furthermore, we extended the Hartle's method to study the differential rotating strange star and found that the differential rotation is an effective way to get massive strange star.

Using the recently obtained holographic cosmic duality, we reached a reasonable quantitative agreement between predictions of the Cosmic Microwave Background Radiation at small l and the WMAP observations, showing the power of the holographic idea. We also got constraints on the dark energy and its behaviour as a function of the redshift upon relating it to the small l CMB spectrum. For a redshift independent dark energy, our constraint is consistent with the supernova results, which again shows the correctness of the cosmic duality prescription. We have also extended our study to the redshift dependence of the dark energy.

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.

Any connection between central activity and the large-scale dynamics of disk galaxies requires an efficient mechanism to remove angular momentum from the orbiting material. The only viable means of achieving inflow from kiloparsec scales is through gravitational stresses created by bars and/or mergers. The inflow of gas in bars today appears to stall at a radius of few hundred parsec, however, forming a nuclear ring. Here we suggest that bars in the early Universe may have avoided this problem, and propose that the progenitors of central supermassive black holes (SMBHs) are created by gas that is driven deep into the centers of galaxies by bars in the early stages of disk formation. The coincidence of the QSO epoch with galaxy formation, the short lifetimes of QSOs, and the existence of SMBHs in the centers of most bright galaxies are all naturally accounted for by disk dynamics in this model. The progenitor SMBHs are the seeds for brighter QSO flares during galaxy mergers. We present a new study of bar weakening by central mass concentrations, which shows that bars are less easily destroyed than previously thought. An extremely massive and compact central mass can, however, dissolve the bar, which creates a pseudo-bulge component in the center of the disk.

Many barred galaxies today harbor massive concentrations of gas in their centers, and some are known to possess supermassive black holes (SBHs) 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. Here we report an extensive study of the effects of central masses on bars in high-quality N-body simulations of galaxies. We study models containing both strong and weak bars, and quantify the change in bar amplitude as a mass concentration is grown in their centers. We find that bars are more robust than previously thought. We have varied the growth rate of the central mass, its final mass and the degree of concentration of the mass. Our main conclusions are: (1) the central mass has to be as large as several percent of the disk mass to completely destroy the bar; (2) for a given mass, objects whose scale-length is a few pc or less cause the greatest reduction in bar amplitude, while significantly more diffuse objects have a lesser effect; and (3) the bar amplitude always decreases as the central mass is grown, and continues to decay thereafter on a cosmological timescale. Thus the masses of SBHs are probably too small, even when dressed with a stellar cusp, to affect the bar amplitude significantly.