We investigate the relations between the properties of bars and their host galaxies in a sample of 77 nearby barred galaxies, spanning a wide range of morphological types and luminosities, with 34 SB0-SBa and 43 SBab-SBc galaxies. The sample includes all the galaxies with reliable direct measurement of their bar pattern speed based on long-slit or integral-field stellar spectroscopy using the Tremaine-Weinberg method. We limited our analysis to the galaxies with a relatively small relative error on the bar pattern speed (smaller than 50 per cent) and not hosting an ultrafast bar. For each galaxy, we collected the radius, strength, pattern speed, corotation radius, and rotation rate for the bar and we also collected the Hubble type and absolute SDSS r-band magnitude. We also used literature bulge-to-total luminosity ratio for a subsample of 53 galaxies with an available photometric decomposition. We confirmed earlier observational findings that longer bars rotate with lower bar pattern speeds, shorter bars are weaker, and bars with a small bar rotation rate rotate with higher bar pattern speeds and have smaller corotation radii. In addition, we found that stronger bars rotate with lower bar pattern speeds, as predicted from the interchange of angular momentum during bar evolution, which in turn may depend on different galaxy properties. Moreover, we report that brighter galaxies host longer bars, which rotate with lower bar pattern speeds and have larger corotation radii. This result is in agreement with a scenario of downsizing in bar formation, if more massive galaxies formed earlier and had sufficient time to slow down, grow in length, and push corotation outwards.

We recently developed an automated method, auto-GMM to decompose simulated galaxies. It extracts kinematic structures in an accurate, efficient, and unsupervised way. We use auto-GMM to study the stellar kinematic structures of disk galaxies from the TNG100 run of IllustrisTNG. We identify four to five structures that are commonly present among the diverse galaxy population. Structures having strong to moderate rotation are defined as cold and warm disks, respectively. Spheroidal structures dominated by random motions are classified as bulges or stellar halos, depending on how tightly bound they are. Disky bulges are structures that have moderate rotation but compact morphology. Across all disky galaxies and accounting for the stellar mass within 3 half-mass radii, the kinematic spheroidal structures, obtained by summing up stars of bulges and halos, contribute ~45% of the total stellar mass, while the disky structures constitute 55%. This study also provides important insights about the relationship between kinematically and morphologically derived galactic structures. Comparing the morphology of kinematic structures with that of traditional bulge+disk decomposition, we conclude: (1) the morphologically decomposed bulges are composite structures comprised of a slowly rotating bulge, an inner halo, and a disky bulge; (2) kinematically disky bulges, akin to what are commonly called pseudo bulges in observations, are compact disk-like components that have rotation similar to warm disks; (3) halos contribute almost 30% of the surface density of the outer part of morphological disks when viewed face-on; and (4) both cold and warm disks are often truncated in central regions.

The Milky Way's progenitor experienced several merger events which left their imprints on the stellar halo, including the Gaia-Sausage/Enceladus. Recently, it has been proposed that this event perturbed the proto-disk and gave rise to a metal rich ([Fe/H] $>-1$) low angular momentum ($v_{\phi} < 100$ km/s) stellar population. These stars have dynamical and chemical properties different from the accreted stellar halo, but are continuous with the canonical thick disk. In this letter, we use a hydrodynamical simulation of an isolated galaxy which develops clumps that produce a bimodal thin$+$thick disk chemistry to explore whether it forms such a population. We demonstrate clump scattering forms a metal-rich, low angular momentum population, without the need for a major merger. We show that, in the simulation, these stars have chemistry, kinematics and density distribution in good agreement with those in the Milky Way.

Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host's disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation.

A scenario for the formation of the bi-modality in the chemical space [$\alpha$/Fe] vs [Fe/H] of the Milky Way was recently proposed in which $\alpha$-enhanced stars are produced early and quickly in clumps. Besides accelerating the enrichment of the medium with $\alpha$-elements, these clumps scatter the old stars, converting in-plane to vertical motion, forming a geometric thick disc. In this paper, by means of a detailed analysis of the data from smooth particle hydrodynamical simulations, we investigate the geometric properties (in particular of the chemical thick disc) produced in this scenario. For mono-age populations we show that the surface radial density profiles of high-[$\alpha$/Fe] stars are well described by single exponentials, while that of low-[$\alpha$/Fe] stars require broken exponentials. This break is sharp for young populations and broadens for older ones. The position of the break does not depend significantly on age. The vertical density profiles of mono-age populations are characterized by single exponentials, which flare significantly for low-[$\alpha$/Fe] stars but only weakly (or not at all) for high-[$\alpha$/Fe] stars. For low-[$\alpha$/Fe] stars, the flaring level decreases with age, while for high-[$\alpha$/Fe] stars it weakly increases with age (although with large uncertainties). All these properties are in agreement with observational results recently reported for the Milky Way, making this a plausible scenario for the formation of the Galactic thick disc.

Vertically thickened bars, observed in the form of boxy/peanut (B/P) bulges, are found in the majority of massive barred disc galaxies in the local Universe, including our own. B/P bulges indicate that their host bars have suffered violent bending instabilities driven by anisotropic velocity distributions. We investigate for the first time how the frequency of B/P bulges in barred galaxies evolves from $z = 1$ to $z\approx 0$, using a large sample of non-edge-on galaxies with masses $M_{\star} > 10^{10}\:M_{\odot}$, selected from the HST COSMOS survey. We find the observed fraction increases from $0^{+3.6}_{-0.0}\%$ at $z = 1$ to $37.8^{+5.4}_{-5.1}\%$ at $z = 0.2$. We account for problems identifying B/P bulges in galaxies with low inclinations and unfavourable bar orientations, and due to redshift-dependent observational biases with the help of a sample from the Sloan Digital Sky Survey, matched in resolution, rest-frame band, signal-to-noise ratio and stellar mass and analysed in the same fashion. From this, we estimate that the true fraction of barred galaxies with B/P bulges increases from $\sim 10\%$ at $z \approx 1$ to $\sim 70\%$ at $z = 0$. In agreement with previous results for nearby galaxies, we find a strong dependence of the presence of a B/P bulge on galaxy stellar mass. This trend is observed in both local and high-redshift galaxies, indicating that it is an important indicator of vertical instabilities across a large fraction of the age of the Universe. We propose that galaxy formation processes regulate the thickness of galaxy discs, which in turn affect which galaxies experience violent bending instabilities of the bar.

Short (inner) bars of sub-kiloparsec radius have been hypothesized to be an important mechanism for driving gas inflows to small scales, thus feeding central black holes. Recent numerical simulations have shown that the growth of central black holes in galaxies can destroy short bars, when the black hole reaches a mass of $\sim 0.1\%$ of the total stellar mass of the galaxy. We study $N$-body simulations of galaxies with single and double bars to track the long-term evolution of the central stellar mass distribution. We find that the destruction of the short bar contributes significantly to the growth of the bulge. The final bulge mass is roughly equal to the sum of the masses of the initial pseudo bulge and short bar. The initially boxy/peanut-shaped bulge of S\'ersic index $n\lesssim1$ is transformed into a more massive, compact structure that bears many similarities to a classical bulge, in terms of its morphology ($n \approx 2$), kinematics (dispersion-dominated, isotropic), and location on standard scaling relations (Kormendy relation, mass-size relation, and correlations between black hole mass and bulge stellar mass and velocity dispersion). Our proposed channel for forming classical bulges relies solely on the destruction of short bars without any reliance on mergers. We suggest that some of the less massive, less compact classical bulges were formed in this manner.

We introduce the study of box/peanut (B/P) bulges in the action space of the initial axisymmetric system. We explore where populations with different actions end up once a bar forms and a B/P bulge develops. We find that the density bimodality due to the B/P bulge (the X-shape) is better traced by populations with low radial, JR,0, or vertical, Jz,0, actions, or high azimuthal action, J{\phi},0. Generally populations separated by JR,0 have a greater variation in bar strength and vertical heating than those separated by Jz,0. While the bar substantially weakens the initial vertical gradient of Jz,0, it also drives a strikingly monotonic vertical profile of JR,0. We then use these results to guide us in assigning metallicity to star particles in a pure N-body model. Because stellar metallicity in unbarred galaxies depends on age as well as radial and vertical positions, the initial actions are particularly well suited for assigning metallicities. We argue that assigning metallicities based on single actions, or on positions, results in metallicity distributions inconsistent with those observed in real galaxies. We therefore use all three actions to assign metallicity to an N-body model by comparing with the actions of a star-forming, unbarred simulation. The resulting metallicity distribution is pinched on the vertical axis, has a realistic vertical gradient and has a stronger X-shape in metal-rich populations, as found in real galaxies.

We present surface photometry and stellar kinematics of NGC 4264, a lenticular galaxy in the region of the Virgo Cluster undergoing a tidal interaction with its neighbour, NGC 4261. We measured the bar radius and strength from SDSS imaging and the bar pattern speed from MUSE integral-field spectroscopy. We find that NGC 4264 hosts a strong and large bar, which is rotating fast. The accurate measurement of the bar rotation rate allows us to exclude that the formation of the bar was triggered by the ongoing interaction.

Galaxies host a wide array of internal stellar components, which need to be decomposed accurately in order to understand their formation and evolution. While significant progress has been made with recent integral-field spectroscopic surveys of nearby galaxies, much can be learned from analyzing the large sets of realistic galaxies now available through state-of-the-art hydrodynamical cosmological simulations. We present an unsupervised machine learning algorithm, named auto-GMM, based on Gaussian mixture models, to isolate intrinsic structures in simulated galaxies based on their kinematic phase space. For each galaxy, the number of Gaussian components allowed by the data is determined through a modified Bayesian information criterion. We test our method by applying it to prototype galaxies selected from the cosmological simulation IllustrisTNG. Our method can effectively decompose most galactic structures. The intrinsic structures of simulated galaxies can be inferred statistically by non-human supervised identification of galaxy structures. We successfully identify four kinds of intrinsic structures: cold disks, warm disks, bulges, and halos. Our method fails for barred galaxies because of the complex kinematics of particles moving on bar orbits.

We aim at investigating the formation process of weak bars by measuring their properties in a sample of 29 nearby SAB galaxies, spanning a wide range of morphological types and luminosities. The sample galaxies were selected to have an intermediate inclination, a bar at an intermediate angle between the disc minor and major axes, and an undisturbed morphology and kinematics to allow the direct measurement of the bar pattern speed. Combining our analysis with previous studies, we compared the properties of weak and strong bars. We measured the bar radius and strength from the r-band images available in SDSS and bar pattern speed and corotation radius from the stellar kinematics obtained by CALIFA. We derived the bar rotation rate as the ratio between the corotation and bar radii. Thirteen out of 29 galaxies, which were morphologically classified as SABs from a visual inspection, do not actually host a bar component or their central elongated component is not in rigid rotation. We successfully derived the bar pattern speed in 16 objects. Two of them host an ultrafast bar. Using the bar strength to differentiate weak and strong bars, we found that the SABs host shorter bars with smaller corotation radii than their strongly barred counterparts. Weak and strong bars have similar bar pattern speeds and rotation rates, which are all consistent with being fast. We did not observe any difference between the bulge prominence in SAB and SB galaxies, whereas nearly all the weak bars reside in the disc inner parts, contrary to strong bars. We ruled out that the bar weakening is only related to the bulge prominence and that the formation of weak bars is triggered by the tidal interaction with a companion. Our observational results suggest that weak bars may be evolved systems exchanging less angular momentum with other galactic components than strong bars.

The direction of the spin vectors of disk galaxies change over time. We present the tilting rate of a sample of galaxies in the NIHAO suite of cosmological hydrodynamical simulations. Galaxies in our sample have been selected to be isolated and to have well determined spins. We compare the tilting rates to the predicted observing limit of Gaia, finding that our entire sample lies above the limit, in agreement with previous work. To test the role of dark matter and of gas we compare the weighted Pearson's correlation coefficients between the tilting rates and various properties. We find no correlation between the dark halo's tilting rate, shape, or misalignment with respect to the disc, and the tilting rate of the stellar disc. Therefore, we argue that, in the presence of gas, the dark halo plays a negligible role in the tilting of the stellar disc. On the other hand, we find a strong correlation between the tilting rate of the stellar disc and the misalignment of the cold gas warp. Adding the stellar mass fraction improves the correlation, while none of the dark matter's properties together with the cold gas misalignment improves the correlation to any significant extent. This implies that the gas cooling onto the disc is the principal driver of disc tilting.

We present surface photometry and stellar kinematics of NGC 4264, a barred lenticular galaxy in the region of the Virgo Cluster undergoing a tidal interaction with one of its neighbours, NGC 4261. We measured the bar radius (a_bar=3.2 +/-0.5 kpc) and strength (S_bar=0.31+/-0.04) of NGC 4264 from Sloan Digital Sky Survey imaging and its bar pattern speed (Omega_bar=71+/-4 km/s/kpc) using the Tremaine-Weinberg method with stellar-absorption integral-field spectroscopy performed with the Multi Unit Spectroscopic Explorer at the Very Large Telescope. We derived the circular velocity (V_circ=189+/-10 km/s) by correcting the stellar streaming velocity for asymmetric drift and calculated the corotation radius (R_cor=2.8+/-0.2 kpc) from the bar pattern speed. Finally, we estimated the bar rotation rate (R_cor/a_bar=0.88+/-0.23). We find that NGC 4264 hosts a strong and large bar extending out to the corotation radius. This means that the bar is rotating as fast as it can like nearly all the other bars measured so far even when the systematic error due to the uncertainty on the disc position angle is taken into account. The accurate measurement of the bar rotation rate allows us to infer that the formation of the bar of NGC 4264 was due to self-generated internal processes and not triggered by the ongoing interaction.

We present the late-time evolution of m12m, a cosmological simulation of a Milky Way-like galaxy from the FIRE project. The simulation forms a bar after redshift z = 0.2. We show that the evolution of the model exhibits behaviours typical of kinematic fractionation, with a bar weaker in older populations, an X-shape traced by the younger, metal-rich populations and a prominent X-shape in the edge-on mean metallicity map. Because of the late formation of the bar in m12m, stars forming after 10 Gyr (z = 0.34) significantly contaminate the bulge, at a level higher than is observed at high latitudes in the Milky Way, implying that its bar cannot have formed as late as in m12m. We also study the model's vertex deviation of the velocity ellipsoid as a function of stellar metallicity and age in the equivalent of Baade's Window. The formation of the bar leads to a non-zero vertex deviation. We find that metal-rich stars have a large vertex deviation (~ 40 degrees), which becomes negligible for metal-poor stars, a trend also found in the Milky Way, despite not matching in detail. We demonstrate that the vertex deviation also varies with stellar age and is large for stars as old as 9 Gyr, while 13 Gyr old stars have negligible vertex deviation. When we exclude stars that have been accreted, the vertex deviation is not significantly changed, demonstrating that the observed variation of vertex deviation with metallicity is not necessarily due to an accreted population.

Compact ellipticals (cEs) are outliers from the scaling relations of early-type galaxies, particularly the mass-metallicity relation which is an important outcome of feedback. The formation of such low-mass, but metal-rich and compact, objects is a long-standing puzzle. Using a pair of high-resolution N-body+gas simulations, we investigate the evolution of a gas-rich low-mass galaxy on a highly radial orbit around a massive host galaxy. As the infalling low-mass galaxy passes through the host's corona at supersonic speeds, its diffuse gas outskirts are stripped by ram pressure, as expected. However, the compactness increases rapidly because of bursty star formation in the gas tidally driven to the centre. The metal-rich gas produced by supernovae and stellar winds is confined by the ram pressure from the surrounding environment, leading to subsequent generations of stars being more metal-rich. After the gas is depleted, tidal interactions enhance the metallicity further via the stripping of weakly bound, old, and metal-poor stars, while the size of the satellite is changed only modestly. The outcome is a metal-rich cE that is consistent with observations. These results argue that classical cEs are neither the stripped remnants of much more massive galaxies nor the merger remnants of normal dwarfs. We present observable predictions that can be used to test our model.

We study the physical properties of the ionized gas in local disks using the sample of 38 nearby $\sim10^{8.5-11.2}$M$_\odot$ Star-Forming Main Sequence (SFMS) galaxies observed so far as part of the MUSE Atlas of Disks (MAD). Specifically, we use all strong emission lines in the MUSE wavelength range 4650-9300 \AA\ to investigate the resolved ionized gas properties on $\sim$100 pc scales. This spatial resolution enables us to disentangle HII regions from the Diffuse Ionized Gas (DIG) in the computation of gas metallicities and Star Formation Rates (SFRs) of star forming regions. The gas metallicities generally decrease with radius. The metallicity of the HII regions is on average $\sim$0.1 dex higher than that of the DIG, but the metallicity radial gradient in both components is similar. The mean metallicities within the inner galaxy cores correlate with the total stellar mass of the galaxies. On our <100 pc scales, we find two correlations previously reported at kpc scales: a spatially resolved Mass-Metallicity Relation (RMZR) and a spatially resolved SFMS (RSFMS). We find no secondary dependency of the RMZR with the SFR density. We find that both resolved relations have a local origin, as they do not depend on the total stellar mass. The observational results of this paper are consistent with the inside-out scenario for the growth of galactic disks.

The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies.

We recommend configuring the LSST coverage of the inner Galactic plane to allow the production of the definitive age/metallicity map of the Galactic bulge from LSST data, matched to external surveys where appropriate. This will allow the formation history of the Galactic bulge to be reconstructed, as well as furnishing a huge legacy dataset to support one of the key LSST science goals (Mapping the Milky Way). We recommend precise multi-color $grizy$~photometry as deeply as spatial crowding will allow (ideally completing early in the project), complemented by a single-filter survey spread over the entire 10-year time baseline. Both strands should cover as broad an area within the Bulge as possible, with image quality sufficient to reach at least the bulge main sequence turn-off in seeing-limited observations. We specify metrics and a figure of merit by which candidate observing strategies could be evaluated with respect to `static' bulge science (proper motions and photometry)

The main aim of this experiment is to provide a complete census of old (t > 10 Gyr, RR Lyrae, type II Cepheids, red horizontal branch), intermediate age (red clump, Miras) and young (classical Cepheids) stellar tracers across the Galactic Bulge. To fully exploit the unique photometric quality of LSST images, we plan to perform a Shallow minisurvey (ugrizy, -20 < l < 20 deg, -15 < b < 10 deg) and a Deep minisurvey (izy, -20 < l < 20 deg, -3 < b < 3 deg). The former one is aimed at constraining the 3D structure of the galactic Bulge across the four quadrants, and in particular, the transition between inner and outer Bulge. The u,g,r,i,z,y LSST bands provide fundamental diagnostics to constrain the evolutionary properties of low and intermediate-mass stars when moving from a metal-poor to a metal-rich regime. The deep minisurvey is aimed at tracing RR Lyrae, Red Clump stars, Miras and classical Cepheids in highly reddened regions of the Galactic center. These images will allow us to investigate the role that baryonic mass and dark matter played in the early formation and evolution of the MW.

Red clump stars are commonly used to map the reddening and morphology of the inner regions of the Milky Way. We use the new photometric catalogues of the VISTA Variables in the V\'ia L\'actea survey to achieve twice the spatial resolution of previous reddening maps for Galactic longitudes $-10^{\circ}<l<10^{\circ}$ and latitudes $-1.5^{\circ}<b<1.5^{\circ}$. We use these de-reddened catalogues to construct the $K_{s}$ luminosity function around the red clump in the Galactic plane. We show that the secondary peak (fainter than the red clump) detected in these regions does not correspond to the bulge red-giant branch bump alone, as previously interpreted. Instead, this fainter clump corresponds largely to the over-density of red clump stars tracing the spiral arm structure behind the Galactic bar. This result suggests that studies aiming to characterise the bulge red-giant branch bump should avoid low galactic latitudes ($|b|< 2^{\circ}$), where the background red clump population contributes significant contamination. It furthermore highlights the need to include this structural component in future modelling of the Galactic bar

We report the discovery from our IRSF/SIRIUS Near-Infrared (NIR) variability survey of five new classical Cepheids located in the inner Galactic Disk, at longitude $l\simeq -40^{\circ}$. The new Cepheids are unique in probing the kinematics and metallicity of young stars at the transition between the inner Disk and the minor axis of the central Bar, where they are expected to be less affected by its dynamical influence.This is also the first time that metallicity of Cepheids is estimated on the basis of medium-resolution ($R\sim3,000$) NIR spectra, and we validated our results with data in the literature, finding a minimal dependence on the adopted spectroscopic diagnostics. This result is very promising for using Cepheids as stellar proxy of the present-time chemical content of the obscured regions in the Disk. We found that the three Cepheids within 8--10 kpc from us have metallicities consistent with the mean radial metallicity gradient, and kinematics consistent with the Galactic rotation curve. Instead, the closest ($\sim$4 kpc)/farthest ($\sim$12 kpc) Cepheids have significant negative/positive residuals, both in velocity and in iron content. We discuss the possibility that such residuals are related to large-scale dynamical instabilities, induced by the bar/spiral-arm pattern, but the current sample is too limited to reach firm conclusions.

We present a detailed study of the nuclear star clusters (NSCs) and massive black holes (BHs) of four of the nearest low-mass early-type galaxies: M32, NGC205, NGC5012, and NGC5206. We measure dynamical masses of both the BHs and NSCs in these galaxies using Gemini/NIFS or VLT/SINFONI stellar kinematics, Hubble Space Telescope (HST) imaging, and Jeans Anisotropic Models. We detect massive BHs in M32, NGC5102, and NGC5206, while in NGC205, we find only an upper limit. These BH mass estimates are consistent with previous measurements in M32 and NGC205, while those in NGC5102 and NGC5206 are estimated for the first time, and both found to be $<$$10^6~M_{\odot}$. This adds to just a handful of galaxies with dynamically measured sub-million $M_{\odot}$ central BHs. Combining these BH detections with our recent work on NGC404's BH, we find that 80\% (4/5) of nearby, low-mass ($10^9-10^{10}~M_{\odot}$; $\sigma_{\star}\sim20-70$ km/s) early-type galaxies host BHs. Such a high occupation fraction suggests the BH seeds formed in the early epoch of cosmic assembly likely resulted in abundant seeds, favoring a low-mass seed mechanism of the remnants, most likely from the first generation of massive stars. We find dynamical masses of the NSCs ranging from $2-80\times10^6~M_{\odot}$ and compare these masses to scaling relations for NSCs based primarily on photometric mass estimates. Color gradients suggest younger stellar populations lie at the centers of the NSCs in three of the four galaxies (NGC205, NGC5102, and NGC5206), while the morphology of two are complex and are best-fit with multiple morphological components (NGC5102 and NGC5206). The NSC kinematics show they are rotating, especially in M32 and NGC5102 ($V/\sigma_{\star}\sim0.7$).

We report results from an exploratory study implementing a new probe of Galactic evolution using archival Hubble Space Telescope imaging observations. Precise proper motions are combined with photometric relative metallicity and temperature indices, to produce the proper motion rotation curves of the Galactic bulge separately for metal-poor and metal-rich Main Sequence samples. This provides a "pencil-beam" complement to large-scale wide-field surveys, which to-date have focused on the more traditional bright Giant Branch tracers. We find strong evidence that the Galactic bulge rotation curves drawn from "Metal-rich" and "Metal-poor" samples are indeed discrepant. The "Metal-rich" sample shows greater rotation amplitude and a steeper gradient against line of sight distance, as possibly a stronger central concentration along the line of sight. This may represent a new detection of differing orbital anisotropy between metal-rich and metal-poor bulge objects. We also investigate selection effects that would be implied for the longitudinal proper motion cut often used to isolate a "pure-bulge" sample. Extensive investigation of synthetic stellar populations suggest that instrumental and observational artifacts are unlikely to account for the observed rotation curve differences. Thus, proper motion-based rotation curves can be used to probe chemo-dynamical correlations for Main Sequence tracer stars, which are orders of magnitude more numerous in the Galactic Bulge than the bright Giant Branch tracers. We discuss briefly the prospect of using this new tool to constrain detailed models of Galactic formation and evolution.

The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.

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.

In Debattista et al. (2015), we proposed that a kiloparsec-scale nuclear disc is responsible for the high-velocity secondary peak in the stellar line-of-sight velocity distributions (LOSVDs) seen at positive longitudes in the bulge by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). Here, we make further qualitative but distinctive predictions of the kinematic properties of a nuclear disc, including for the LOSVDs at negative longitudes (which APOGEE-2 will observe) and examine the proper motions throughout the disc. Since a nuclear ring is also able to produce similar high-velocity LOSVD peaks, we present predictions for the proper motion signatures which distinguish between a nuclear disc and a nuclear ring. We also demonstrate that the stars in a nuclear disc, which would be on x2 orbits perpendicular to the bar, can remain on these orbits for a long time and can therefore be old. We show that such (old) nuclear discs of comparable size exist in external galaxies.

We present an extensively updated version of the purely ray-tracing 3D dust radiation transfer code DART-Ray. The new version includes five major upgrades : 1) a series of optimizations for the ray-angular density and the scattered radiation source function; 2) the implementation of several data and task parallelizations using hybrid MPI+OpenMP schemes; 3) the inclusion of dust self-heating; 4) the ability to produce surface brightness maps for observers within the models in HEALPix format; 5) the possibility to set the expected numerical accuracy already at the start of the calculation. We tested the updated code with benchmark models where the dust self-heating is not negligible. Furthermore, we performed a study of the extent of the source influence volumes, using galaxy models, which are critical in determining the efficiency of the DART-Ray algorithm. The new code is publicly available, documented for both users and developers, and accompanied by several programmes to create input grids for different model geometries and to import the results of N-body and SPH simulations. These programmes can be easily adapted to different input geometries, and for different dust models or stellar emission libraries.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

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 present tilting rates for galaxies comparable to the Milky Way (MW) in a $\Lambda$ cold dark matter cosmological hydrodynamical simulation, and compare these with the predicted tilting rate detection limit of the {\it Gaia} satellite $0.28\degrees$Gyr$^{-1}$. We first identify galaxies with mass comparable to the MW ($9 \times 10^{11} \le M_{200} \le 1.2 \times 10^{12} \Msun $) and consider the tilting rates between $z=0.3$ and $0$. This sample yields a tilting rate of $7.6\degrees \pm 4.5\degrees$Gyr$^{-1}$. We constrain our sample further to exclude any galaxies that have high stellar accretion during the same time. We still find significant tilting, with an average rate of $6.3\degrees$Gyr$^{-1}$. Both subsamples tilt with rates significantly above {\it Gaia}'s predicted detection limit. We show that our sample of galaxies covers a wide range of environments, including some similar to the MW's. We find galaxies in denser regions tilt with higher rates then galaxies in less dense regions. We also find correlations between the angular misalignment of the hot gas corona, and the tilting rate. {\it Gaia} is likely to be able to directly measure tilting in the MW. Such a detection will provide an important constraint on the environment of the MW, including the rate of gas cooling onto the disc, the shape and orientation of its dark matter halo, and the mass of the Large Magellanic Cloud. Conversely, failure to detect tilting may suggest the MW is in a very quiet configuration.

Using a cosmological $N$-body simulation, we investigate the origin and distribution of stars in the intracluster light (ICL) of a Fornax-like cluster. In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200 \simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling subhalos with models that include spheroid and disc components. As they fall into the cluster, the stars in some of these galaxies are stripped from their hosts, and form the ICL. We consider the separate contributions to the ICL from stars which originate in the haloes and the discs of the galaxies. We find that disc ICL stars are more centrally concentrated than halo ICL stars. The majority of the disc ICL stars are associated with one initially disc-dominated galaxy that falls to the centre of the cluster and is heavily disrupted, producing part of the cD galaxy. At radial distances greater than 200kpc, well beyond the stellar envelope of the cD galaxy, stars formerly from the stellar haloes of galaxies dominate the ICL. Therefore at large distances, the ICL population is dominated by older stars.

We present a novel interpretation of the previously puzzling different behaviours of stellar populations of the Milky Way's bulge. We first show, by means of pure N-body simulations, that initially co-spatial stellar populations with different in-plane random motions separate when a bar forms. The radially cooler populations form a strong bar, and are vertically thin and peanut-shaped, while the hotter populations form a weaker bar and become a vertically thicker box. We demonstrate that it is the radial, not the vertical, velocity dispersion that dominates this evolution. Assuming that early stellar discs heat rapidly as they form, then both the in-plane and vertical random motions correlate with stellar age and chemistry, leading to different density distributions for metal-rich and metal-poor stars. We then use a high-resolution simulation, in which all stars form out of gas, to demonstrate that this is what happens. When we apply these results to the Milky Way we show that a very broad range of observed trends for ages, densities, kinematics and chemistries, that have been presented as evidence for contradictory paths to the formation of the bulge, are in fact consistent with a bulge which formed from a continuum of disc stellar populations which were kinematically separated by the bar. For the first time we are able to account for the bulge's main trends via a model in which the bulge formed largely in situ. Since the model is generic, we also predict the general appearance of stellar population maps of external edge-on galaxies.

Stellar migration, whether due to trapping by transient spirals (churning), or to scattering by non-axisymmetric perturbations, has been proposed to explain the presence of stars in outer disks. After a review of the basic theory, we present compelling, but not yet conclusive, evidence that churning has been important in the outer disks of galaxies with type II (down-bending) profiles, while scattering has produced the outer disks of type III (up-bending) galaxies. In contrast, field galaxies with type I (pure exponential) profiles appear to not have experienced substantial migration. We conclude by suggesting work that would improve our understanding of the origin of outer disks.

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.

Spiral arms are the most singular features in disc galaxies. These structures can exhibit different patterns, namely grand design and flocculent arms, with easily distinguishable characteristics. However, their origin and the mechanisms shaping them are unclear. The overall role of spirals in the chemical evolution of disc galaxies is another unsolved question. In particular, it has not been fully explored if the \hii\,regions of spiral arms present different properties from those located in the interarm regions. Here we analyse the radial oxygen abundance gradient of the arm and interarm star forming regions of 63 face-on spiral galaxies using CALIFA Integral Field Spectroscopy data. We focus the analysis on three characteristic parameters of the profile: slope, zero-point, and scatter. The sample is morphologically separated into flocculent versus grand design spirals and barred versus unbarred galaxies. We find subtle but statistically significant differences between the arm and interarm distributions for flocculent galaxies, suggesting that the mechanisms generating the spiral structure in these galaxies may be different to those producing grand design systems, for which no significant differences are found. We also find small differences in barred galaxies, not observed in unbarred systems, hinting that bars may affect the chemical distribution of these galaxies but not strongly enough as to be reflected in the overall abundance distribution. In light of these results, we propose bars and flocculent structure as two distinct mechanisms inducing differences in the abundance distribution between arm and interarm star forming regions.

From a sample of 84 local barred, moderately inclined disc galaxies, we determine the fraction which host boxy or peanut-shaped (B/P) bulges (the vertically thickened inner parts of bars). We find that the frequency of B/P bulges in barred galaxies is a very strong function of stellar mass: 79% of the bars in galaxies with log (M_{star}/M_{sun}) >~ 10.4 have B/P bulges, while only 12% of those in lower-mass galaxies do. (We find a similar dependence in data published by Yoshino & Yamauchi 2015 for edge-on galaxies.) There are also strong trends with other galaxy parameters -- e.g., Hubble type: 77% of S0-Sbc bars, but only 15% of Sc-Sd bars, have B/P bulges -- but these appear to be side effects of the correlations of these parameters with stellar mass. In particular, despite indications from models that a high gas content can suppress bar buckling, we find no evidence that the (atomic) gas mass ratio M_{atomic}/M_{star} affects the presence of B/P bulges, once the stellar-mass dependence is controlled for. The semi-major axes of B/P bulges range from one-quarter to three-quarters of the full bar size, with a mean of R_{box}/L_{bar} = 0.42 +/- 0.09 and R_{box}/a_{max} = 0.53 +/- 0.12 (where R_{box} is the size of the B/P bulge and a_{max} and L_{bar} are lower and upper limits on the size of the bar).

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.

A minimum in stellar velocity dispersion is often observed in the central regions of disc galaxies. To investigate the origin of this feature, known as a {\sigma}-drop, we analyse the stellar kinematics of a high-resolution N-body + smooth particle hydrodynamical simulation, which models the secular evolution of an unbarred disc galaxy. We compared the intrinsic mass-weighted kinematics to the recovered luminosity-weighted ones. The latter were obtained by analysing synthetic spectra produced by a new code, SYNTRA, that generates synthetic spectra by assigning a stellar population synthesis model to each star particle based on its age and metallicity. The kinematics were derived from the synthetic spectra as in real spectra to mimic the kinematic analysis of real galaxies. We found that the recovered luminosity-weighted kinematics in the centre of the simulated galaxy are biased to higher rotation velocities and lower velocity dispersions due to the presence of young stars in a thin and kinematically cool disc, and are ultimately responsible for the {\sigma}-drop.

Bulges of edge-on galaxies are often boxy/peanut-shaped (B/PS), and unsharp masks reveal the presence of an X shape. Simulations show that these shapes can be produced by dynamical processes driven by a bar which vertically thickens the centre. In the Milky Way, which contains such a B/PS bulge, the X-shaped structure is traced by the metal-rich stars but not by the metal-poor ones. Recently Debattista et al. (2016) interpreted this property as a result of the varying effect of the bar on stellar populations with different starting kinematics. This kinematic fractionation model predicts that cooler populations at the time of bar formation go on to trace the X shape, whereas hotter populations are more uniformly distributed. As this prediction is not specific to the Milky Way, we test it with MUSE observations of the B/PS bulge in the nearby galaxy NGC 4710. We show that the metallicity map is more peanut-shaped than the density distribution itself, in good agreement with the prediction. This result indicates that the X-shaped structure in B/PS bulges is formed of relatively metal-rich stars that have been vertically redistributed by the bar, whereas the metal-poor stars have a more uniform, box-shaped distribution.

Using $N$-body+SPH simulations of galaxies falling into a cluster, we study the evolution of their radial density profiles. When evolved in isolation, galaxies develop a type~II (down-bending) profile. In the cluster, the evolution of the profile depends on the minimum cluster-centric radius the galaxy reaches, which controls the degree of ram pressure stripping. If the galaxy falls to $\sim 50\%$ of the virial radius, then the profile remains type~II, but if the galaxy reaches down to $\sim 20\%$ of the virial radius the break weakens and the profile becomes more type~I like. The velocity dispersions are only slightly increased in the cluster simulations compared with the isolated galaxy; random motion therefore cannot be responsible for redistributing material sufficiently to cause the change in the profile type. Instead we find that the joint action of radial migration driven by tidally-induced spirals and the outside-in quenching of star formation due to ram pressure stripping alters the density profile. As a result, this model predicts a flattening of the age profiles amongst cluster lenticulars with type~I profiles, which can be observationally tested. %

The majority of massive disk galaxies, including our own, have stellar bars with vertically thick inner regions -- so-called "boxy/peanut-shaped" (B/P) bulges. The most commonly suggested mechanism for the formation of B/P bulges is a violent vertical "buckling" instability in the bar, something that has been seen in N-body simulations for over twenty years, but never identified in real galaxies. Here, we present the first direct observational evidence for ongoing buckling in two nearby galaxies (NGC 3227 and NGC 4569), including characteristic asymmetric isophotes and (in NGC 4569) stellar-kinematic asymmetries that match buckling in simulations. This confirms that the buckling instability takes place and produces B/P bulges in real galaxies. A toy model of bar evolution yields a local fraction of buckling bars consistent with observations if the buckling phase lasts ~0.5--1 Gyr, in agreement with simulations.

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.

Globular clusters which exhibit chemical and dynamical complexity have been suggested to be the stripped nuclei of dwarf galaxies (e.g., M54, $\omega$ Cen). We use $N$-body simulations of nuclear star clusters forming via the mergers of star clusters to explore the persistence of substructure in the phase space. We find that the observed level of differentiation is difficult to reconcile with the observed if nuclear clusters form wholly out of the mergers of star clusters. Only the star clusters that merged most recently retain sufficiently distinct kinematics to be distinguishable from the rest of the nuclear cluster though the critical factor is the number of merger events not the elapsed time. In situ star formation must therefore be included to explain the observed properties of nuclear star clusters, in good agreement with previous results.

We used the new ESO VLT instrument MUSE to obtain spectral and imaging coverage of NGC 4710. The wide area and excellent sampling of the MUSE integral field spectrograph allows us to investigate the dynamical properties of the X-shaped bulge of NGC 4710 and compare it with the properties of the Milky Way's own X-shaped bulge. We measured the radial velocities, velocity dispersion, and stellar populations using a penalized pixel full spectral fitting technique adopting simple stellar populations models, on a 1' x 1' area centred on the bulge of NGC 4710. We have constructed the velocity maps of the bulge of NGC 4710 and we investigated the presence of vertical metallicity gradients. These properties were compared to those of the Milky Way bulge and as well as to a simulated galaxy with boxy/peanut bulge. We find the line-of-sight velocity maps and 1D rotation curves of the bulge of NGC 4710 to be remarkably similar to those of the Milky Way bulge. Some specific differences that were identified are in good agreement with the expectations from variations in the bar orientation angle. The bulge of NGC 4710 has a boxy-peanut morphology with a pronounced X-shape, showing no indication of any additional spheroidally distributed bulge population, in which we measure a vertical metallicity gradient of 0.35 dex/kpc. The general properties of NGC 4710 are very similar to those observed in the MW bulge. However, it has been suggested that the MW bulge has an additional component that is comprised of the oldest, most metal-poor stars and which is not part of the boxy-peanut bulge structure. Such a population is not observed in NGC 4710, but could be hidden in the integrated light we observed.

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

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.

The finding of a double red clump in the luminosity function of the Milky Way bulge has been interpreted as evidence for an X-shaped structure. Recently, an alternative explanation has been suggested, where the double red clump is an effect of multiple stellar populations in a classical spheroid. In this letter we provide an observational assessment of this scenario and show that it is not consistent with the behaviour of the red clump across different lines of sight, particularly at high distances from the Galactic plane. Instead, we confirm that the shape of the red clump magnitude distribution closely follows the distance distribution expected for an X-shaped bulge at critical Galactic latitudes. We also emphasize some key observational properties of the bulge red clump that should not be neglected in the search for alternative scenarios.

We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.

The Apache Point Observatory Galactic Evolution Experiment has measured the stellar velocities of red giant stars in the inner Milky Way. We confirm that the line of sight velocity distributions (LOSVDs) in the mid-plane exhibit a second peak at high velocities, whereas those at |b| = 2degrees do not. We use a high resolution simulation of a barred galaxy, which crucially includes gas and star formation, to guide our interpretation of the LOSVDs. We show that the data are fully consistent with the presence of a thin, rapidly rotating, nuclear disk extending to ~1 kpc. This nuclear disk is orientated perpendicular to the bar and is likely to be composed of stars on x2 orbits. The gas in the simulation is able to fall onto such orbits, leading to stars populating an orthogonal disk.

Large surveys have shown that red galaxies are preferentially aligned with their halos while blue galaxies have a more isotropic distribution. Since halos generally align with their filaments this introduces a bias in the measurement of the cosmic shear from weak lensing. It is therefore vitally important to understand why this difference arises. We explore the stability of different disc orientations within triaxial halos. We show that, in the absence of gas, the disc orientation is most stable when its spin is along the minor axis of the halo. Instead when gas cools onto a disc it is able to form in almost arbitrary orientation, including off the main planes of the halo (but avoiding an orientation perpendicular to the halo's intermediate axis). Substructure helps gasless galaxies reach alignment with the halo faster, but have less effect on galaxies when gas is cooling onto the disc. Our results provide a novel and natural interpretation for why red, gas poor galaxies are preferentially aligned with their halo, while blue, star-forming, galaxies have nearly random orientations, without requiring a connection between galaxies' current star formation rate and their merger history.

NGC 4395 is a bulgeless spiral galaxy, harboring one of the nearest known type 1 Seyfert nuclei. Although there is no consensus on the mass of its central engine, several estimates suggest it to be one of the lightest massive black holes (MBHs) known. We present the first direct dynamical measurement of the mass of this MBH from a combination of two-dimensional gas kinematic data, obtained with the adaptive optics assisted near infrared integral field spectrograph Gemini/NIFS, and high-resolution multiband photometric data from Hubble Space Telescope's Wide Field Camera 3 (HST/WFC3). We use the photometric data to model the shape and stellar mass-to-light ratio (M/L) of the nuclear star cluster. From the Gemini/NIFS observations, we derive the kinematics of warm molecular hydrogen gas as traced by emission through the H$_2$ 1--0 S(1) transition. These kinematics show a clear rotational signal, with a position angle orthogonal to NGC 4395's radio jet. Our best fitting tilted ring models of the kinematics of the molecular hydrogen gas contain a black hole with mass $M=4_{-3}^{+8}\times 10^5$ M$_\odot$ (3$\sigma$ uncertainties) embedded in a nuclear star cluster of mass $M=2 \times 10^6$ M$_\odot$. Our black hole mass measurement is in excellent agreement with the reverberation mapping mass estimate of Peterson et al. (2005), but shows some tension with other mass measurement methods based on accretion signals.

We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.

Nuclear star clusters are among the densest stellar systems known and are common in both early- and late-type galaxies. They exhibit scaling relations with their host galaxy which may be related to those of supermassive black holes. These may therefore help us to unravel the complex physical processes occurring at the centres of galaxies. The properties of nuclear stellar systems suggest that their formation requires both dissipational and dissipationless processes. They have stellar populations of different ages, from stars as old as their host galaxy to young stars formed in the last 100 Myr. Therefore star formation must be happening either directly in the nuclear star cluster or in its vicinity. The secular processes that fuel the formation of pseudobulges very likely also contributes to nuclear star cluster growth.

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 describe the calculation of the stochastically heated dust emission using the 3D ray-tracing dust radiative transfer code DART-Ray, which is designed to solve the dust radiative transfer problem for galaxies with arbitrary geometries. In order to reduce the time required to derive the non-equilibrium dust emission spectra from each volume element within a model, we implemented an adaptive SED library approach, which we tested for the case of axisymmetric galaxy geometries. To show the capabilities of the code, we applied DART-Ray to a high-resolution N-body+SPH galaxy simulation to predict the appearance of the simulated galaxy at a set of wavelengths from the UV to the sub-mm. We analyse the results to determine the effect of dust on the observed radial and vertical profiles of the stellar emission as well as on the attenuation and scattering of light from the constituent stellar populations. We also quantify the proportion of dust re-radiated stellar light powered by young and old stellar populations, both bolometrically and as a function of infrared wavelength.

The role of gas in the mass assembly at the nuclei of galaxies is still subject to some uncertainty. Stellar nuclear discs bridge the gap between the large-scale galaxy and the central massive objects that reside there. Using a high resolution simulation of a galaxy forming out of gas cooling and settling into a disc, we study the formation and properties of nuclear discs. Gas, driven to the centre by a bar, settles into a rotating star-forming nuclear disc (ND). This ND is thinner, younger, kinematically cooler, and more metal-rich than the surrounding bar. The ND is elliptical and orthogonal to the bar. The complex kinematics in the region of the ND are a result of the superposition of older stars streaming along the bar and younger stars circulating within the ND. The signature of the ND is therefore subtle in the kinematics. Instead the ND stands out clearly in metallicity and age maps. We compare the model to the density and kinematics of real galaxies with NDs finding qualitative similarities. Our results suggest that gas dissipation is very important for forming nuclear structures.

Athanassoula (2014) has claimed that measurements of the ratio of corotation radius to bar length in galaxies do not place a constraint on the disk to halo mass ratio. Her conclusion was based on a series of simulations published by Athanassoula et al. (2013). Here we show that these results are, in fact, consistent with previous work on the slow down of bars due to dynamical friction because gas inflow rearranges the disk mass before the bar settles. It therefore remains true that a bar rotating sufficiently fast that corotation is not far beyond the bar end requires a near maximum disk.

Using N-body simulations, we study the properties of a bar induced in a dwarf galaxy as a result of tidal interaction with the Milky Way. The dwarf is initially composed of a disk embedded in a dark matter halo and we follow its evolution on a typical orbit for 10 Gyr. It undergoes an evolution typical of tidally stirred dwarf galaxies: it loses mass, the stellar component transforms from a disk to a spheroid and the rotation of the stars is partially replaced by random motions. A tidally induced bar forms at the first pericentre passage and survives until the end of the evolution. Fourier decomposition of the distribution of stars reveals that only even modes are significant and preserve a hierarchy so that the bar mode is always the strongest. They show a characteristic profile with a maximum, similar to simulated bars forming in isolated galaxies and observed bars in real galaxies. We adopt the maximum of the bar mode as a measure of the bar strength and we estimate the bar length by comparing the density profiles along the bar and perpendicular to it. Both the bar strength and the bar length decrease with time, mainly at pericentres, as a result of tidal torques acting at those times and not to secular evolution. The pattern speed of the bar varies significantly on a time scale of 1 Gyr and seems to be fully controlled by the orientation of the tidal torque from the Milky Way. The bar is never tidally locked but we discover a hint of a 5/2 orbital resonance between the third and fourth pericentre passage. The speed of the bar expressed in terms of the ratio of the bar length to the corotation radius decreases in the long run so that the bar changes from initially rather fast to slow in the later stages. The boxy/peanut shape is clearly visible in the edge-on view between the second and third pericentre and its occurrence is preceded by a short period of buckling instability.

(Abridged) We study the incidence, as well as the nature, of composite bulges in a sample of 10 face-on barred galaxies to constrain the formation and evolutionary processes of the central regions of disk galaxies. We analyze the morphological, photometric, and kinematic properties of each bulge. Then, by using a case-by-case analysis we identify composite bulges and classify every component into a classical or pseudobulge. In addition, bar-related boxy/peanut (B/P) structures were also identified and characterised. We find only three galaxies hosting a single-component bulge (two pseudobulges and one classical bulge). We find evidence of composite bulges coming in two main types based on their formation: secular-built and merger- and secular-built. We call secular-built to composite bulges made of entirely by structures associated with secular processes such as pseudo bulges, central disks, or B/P bulges. We find four composite bulges of this kind in our sample. On the other hand, merger- and secular-built bulges are those where structures with different formation paths coexist within the same galaxy, i.e., a classical bulge coexisting with a secular-built structure (pseudobulge, central disk, or B/P). Three bulges of this kind were found in the sample. We remark on the importance of detecting kinematic structures such as sigma-drops to identify composite bulges. A large fraction (80%) of galaxies were found to host sigma-drops or sigma-plateaus in our sample revealing their high incidence in barred galaxies. The high frequency of composite bulges in barred galaxies points towards a complex formation and evolutionary scenario. Moreover, the evidence for coexisting merger- and secular-built bulges reinforce this idea. We discuss how the presence of different bulge types, with different formation histories and timescales, can constrain current models of bulge formation.

The redistribution of stars in galactic disks is an important aspect of disk galaxy evolution. Stars that efficiently migrate in such a way that does not also appreciably heat their orbits can drastically affect the stellar populations observed today and therefore influence constraints derived from such observations. Unfortunately, while the theoretical understanding of the migration process is becoming increasingly robust, there are currently few specific observable predictions. As a result, we do not yet have a clear handle on whether the process has been important for the Milky Way in the past or how to constrain it. I discuss some of the expected qualitative outcomes of migration as well as some current controversies.

We use an N-body simulation to study the 3-D density distribution of spirals, and the resulting stellar vertical velocities. Relative to the disc's rotation, the phase of the spiral's peak density away from the mid-plane trails that at the mid-plane. In addition, at fixed radius the density distribution is azimuthally skewed, having a shallower slope on the trailing side inside corotation and switching to shallower on the leading side beyond corotation. The spirals induce non-zero average vertical velocities, <V_z>, as large as <V_z> ~ 10-20 km/s, consistent with recent observations in the Milky Way. The vertical motions are compressive (towards the mid-plane) as stars enter the spiral, and expanding (away from the mid-plane) as they leave it. Since stars enter the spiral on the leading side outside corotation and on the trailing side within corotation, the relative phase of the expanding and compressive motions switches sides at corotation. Moreover, because stars always enter the spiral on the shallow density gradient side and exit on the steeper side, the expanding motions are larger than the compressing motions.

Using simulations of box/peanut- (B/P-) shaped bulges, we explore the nature of the X-shape of the Milky Way's bulge. An X-shape can be associated with a B/P-shaped bulge driven by a bar. By comparing in detail the simulations and the observations we show that the principal kinematic imprint of the X-shape is a minimum in the difference between the near and far side mean line-of-sight velocity along the minor axis. This minimum occurs at around |b| = 4{\deg}, which is close to the lower limit at which the X-shape can be detected. No coherent signature of an X-shape can be found in Galactocentric azimuthal velocities, vertical velocities, or any of the dispersions. After scaling our simulations, we find that a best fit to the BRAVA data leads to a bar angle of 15{\deg}. We also explore a purely geometric method for determining the distance to the Galactic Centre by tracing the arms of the X-shape. We find that we are able to determine this ill-known distance to an accuracy of about 5% per cent with sufficiently accurate distance measurements for the red clump stars in the arms.

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

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

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.

In some scenarios for the formation of the Milky Way bulge the stellar population at the edges of the boxy bulge may be younger than those on the minor axis, or close to the Galactic center. So far the only bulge region where deep color-magnitude diagrams have been obtained is indeed along the minor axis. To overcome this limitation, we aim at age-dating the bulge stellar populations far away from the bulge minor axis. Color-magnitude diagrams and luminosity functions have been obtained from deep near-IR VLT/HAWK-I images taken at the two Southern corners of the boxy bulge, i.e., near the opposite edges of the Galactic bar. The foreground disk contamination has been statistically removed using a pure disk field observed with the same instrument and located approximately at similar Galactic latitudes of the two bulge fields, and 30deg in longitude away from the Galactic center. For each bulge field, mean reddening and distance are determined using the position of red clump stars, and the metallicity distribution is derived photometrically using the color distribution of stars in the upper red giant branch. The resulting metallicity distribution function of both fields peaks around [Fe/H] -0.1 dex, with the bulk of the stellar population having a metallicity within the range: -1 dex < [Fe/H] <+0.4 dex, quite similar to that of other inner bulge fields. Like for the inner fields previously explored, the color-magnitude diagrams of the two bar fields are consistent with their stellar population being older than ~10 Gyr, with no obvious evidence for the presence of a younger population.

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.

Models of the Sagittarius Stream have consistently found that the Milky Way disc is oriented such that its short axis is along the intermediate axis of the triaxial dark matter halo. We attempt to build models of disc galaxies in such an `intermediate-axis orientation'. We do this with three models. In the first two cases we simply rigidly grow a disc in a triaxial halo such that the disc ends up perpendicular to the global intermediate axis. We also attempt to coax a disc to form in an intermediate-axis orientation by producing a gas+dark matter triaxial system with gas angular momentum about the intermediate axis. In all cases we fail to produce systems which remain with stellar angular momentum aligned with the halo's intermediate axis, even when the disc's potential flattens the inner halo such that the disc is everywhere perpendicular to the halo's local minor axis. For one of these unstable simulations we show that the potential is even rounder than the models of the Milky Way potential in the region probed by the Sagittarius Stream. We conclude that the Milky Way's disc is very unlikely to be in an intermediate axis orientation. However we find that a disc can persist off one of the principal planes of the potential. We propose that the disc of the Milky Way must be tilted relative to the principal axes of the dark matter halo. Direct confirmation of this prediction would constitute a critical test of Modified Newtonian Dynamics.

We present evidence that isolated growing discs, subject to internal spiral perturbations, thicken due to both heating \emph{and} radial migration. We show this by demonstrating that the thickness and vertical velocity dispersions of coeval stars depend on their age as well as the change in their radii. While the disc thickens due to internal processes, we find that this induces only a minor amount of flaring. We further demonstrate the consequences of such thickening on the structural properties of stellar populations and find that they qualitatively agree with recent studies of the Milky Way disc.

The nuclei of galaxies often host small stellar discs with scale-lengths of a few tens of parsecs and luminosities up to 10^7 Lsun. To investigate the formation and properties of nuclear stellar discs (NSDs), we look for their presence in a set of N-body simulations studying the dissipationless merging of multiple star clusters in galactic nuclei. A few tens of star clusters with sizes and masses comparable to those of globular clusters observed in the Milky Way are accreted onto a pre-existing nuclear stellar component: either a massive super star cluster or a rapidly rotating, compact disc with a scale-length of a few parsecs, mimicing the variety of observed nuclear structures. Images and kinematic maps of the simulation time-steps are then built and analysed as if they were real and at the distance of the Virgo cluster. We use the Scorza-Bender method to search for the presence of disc structures via photometric decomposition. In one case the merger remnant has all the observed photometric and kinematic properties of NSDs observed in real galaxies. This shows that current observations are consistent with most of the NSD mass being assembled from the migration and accretion of star clusters into the galactic centre. In the other simulation instead, we detect an elongated structure from the unsharp masked image, that does not develop the photometric or kinematic signature of a NSD. Thus, in the context of searches for a disc structure, the Scorza-Bender method is a robust and necessary tool.

It has recently been discovered that the Galactic bulge is X-shaped, with the two southern arms of the X both crossing the lines of sight at l=0 and |b|>4, hence producing a double red clump (RC) in the bulge CMD. Dynamical models predict the formation of X-shaped bulges, as extreme cases of boxy-peanut bulges. However, since X-shaped bulges were known to be present only in external galaxies, models have never been compared to 3D kinematical data for individual stars. We study the orbital motion of Galactic bulge stars, in the two arms of the X in the southern hemisphere. The goal is to provide observational constraints to bulge formation models that predict the formation of X-shapes through bar dynamical instabilities. Radial velocities have been obtained for a sample of 454 bulge giants, roughly equally distributed between the bright and faint RC, in a field at (l,b)=(0,-6). Proper motions were derived for all RC stars in the same field by combining images from two epochs obtained 11 years apart. The proper motions for the spectroscopic subsample are analyzed taking into account the radial velocities and metallicities measured from near-IR Calcium triplet lines, for a subsample of 352 stars. The radial velocity distribution of stars in the bright RC, tracing the closer overdensity of bulge stars, shows an excess of stars moving towards the Sun. Similarly, an excess of stars receding from the Sun is seen in the far overdensity, traced by faint RC stars. This is explained with the presence of stars on elongated orbits, most likely streaming along the arms of the X-shaped bulge. Proper motions for these stars are consistent with qualitative predictions of dynamical models of peanut-shaped bulges. Surprisingly, stars on elongated orbits have preferentially metal poor metallicities, while the metal rich ones, in both overdensities, are preferentially found in more axysimmetric orbits.

We show that direct detection and measurement of the vertically thickened parts of bars (so-called "boxy" or "peanut-shaped" bulges) is possible not only for edge-on galaxies but also for galaxies with moderate inclinations (i < 70 deg), and that examples are relatively common in the nearby universe. Analysis of a sample of 78 nearby, moderately inclined (i < 65 deg) early-type (S0--Sb) barred galaxies shows that the isophotal signature of the box/peanut can usually be detected for inclinations as low as i ~ 40 deg -- and in exceptional cases down to i ~ 30 deg. In agreement with the predictions from N-body simulations, the signature is most easily detectable when the bar's position angle is within ~ 50 deg of the galaxy major axis; in particular, galaxies where the bar lies very close to the minor axis do not show the signature clearly or at all. For galaxies with i = 40--65 deg and relative angles < 45 deg, we find evidence for the signature ~ 2/3 of the time; the true frequency of box/peanut structures in bars may be higher. Comparison with N-body models also allows us to link observed photometric morphology with 3D physical structures, and thus estimate the relative sizes of box/peanut structures and bars. For our local sample, we find that box/peanut structures range in radial size (measured along the bar major axis) from 0.4--3.8 kpc (mean = 1.5 +/- 0.9 kpc) and span 0.26--0.58 of the bar length (mean of 0.38 +/- 0.08). This is a clear observational confirmation that when bars thicken, it is not the entire bar which does so, but only the inner part. This technique can also be used to identify galaxies with bars which have not vertically thickened. We suggest that NGC 3049 and IC 676 may be particularly good examples, and that the fraction of S0--Sb bars which lack box/peanut structures is at least ~ 13%.

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

Recent Hubble Space Telescope (HST) observations have revealed that a majority of active galactic nuclei (AGN) at z ~ 1-3 are resident in isolated disk galaxies, contrary to the usual expectation that AGN are triggered by mergers. Here we develop a new test of the cosmic evolution of supermassive black holes (SMBHs) in disk galaxies by considering the local population of SMBHs. We show that substantial SMBH growth in spiral galaxies is required as disks assemble. SMBHs exhibit a tight relation between their mass and the velocity dispersion of the spheroid within which they reside, the M_BH-sigma relation. In disk galaxies the bulge is the spheroid of interest. We explore the evolution of the M_BH-sigma relation when bulges form together with SMBHs on the M_BH-sigma relation and then slowly reform a disk around them. The formation of the disk compresses the bulge raising its sigma. We present evidence for such compression in the form of larger velocity dispersion of classical bulges compared with elliptical galaxies at the same mass. This compression leads to an offset in the M_BH-sigma relation if it is not accompanied by an increased M_BH. We quantify the expected offset based on photometric data and show that, on average, SMBHs must grow by ~ 50-65% just to remain on the M_BH-sigma relation. We find no significant offset in the M_BH-sigma relations of classical bulges and of ellipticals, implying that SMBHs have been growing along with disks. Our simulations demonstrate that SMBH growth is necessary for the local population of disk galaxies to have remained on the M_BH-sigma relation.

Adaptive optics observations of the flattened nuclear star cluster in the nearby edge-on spiral galaxy NGC 4244 using the Gemini Near-Infrared Integral Field Spectrograph (NIFS) have revealed clear rotation. Using these kinematics plus 2MASS photometry we construct a series of axisymmetric two-component particle dynamical models with our improved version of NMAGIC, a flexible Chi^2-made-to-measure code. The models consist of a nuclear cluster disc embedded within a spheroidal particle population. We find a mass for the nuclear star cluster of M=1.6^+0.5_-0.2 x 10^7 M_sun within ~42.4 pc (2"). We also explore the presence of an intermediate mass black hole and show that models with a black hole as massive as M_bh = 5.0 x 10^5 M_sun are consistent with the available data. Regardless of whether a black hole is present or not, the nuclear cluster is vertically anisotropic (beta_z < 0), as was found with earlier two-integral models. We then use the models as initial conditions for N-body simulations. These simulations show that the nuclear star cluster is stable against non-axisymmetric perturbations. We also explore the effect of the nuclear cluster accreting star clusters at various inclinations. Accretion of a star cluster with mass 13% that of the nuclear cluster is already enough to destroy the vertical anisotropy, regardless of orbital inclination.

We seek to understand the origin of radial migration in spiral galaxies by analyzing in detail the structure and evolution of an idealized, isolated galactic disk. To understand the redistribution of stars, we characterize the time-evolution of properties of spirals that spontaneously form in the disk. Our models unambiguously show that in such disks, single spirals are unlikely, but that a number of transient patterns may coexist in the disk. However, we also show that while spirals are transient in amplitude, at any given time the disk favors patterns of certain pattern speeds. Using several runs with different numerical parameters we show that the properties of spirals that occur spontaneously in the disk do not sensitively depend on resolution. The existence of multiple transient patterns has large implications for the orbits of stars in the disk, and we therefore examine the resonant scattering mechanisms that profoundly alter angular momenta of individual stars. We confirm that the corotation scattering mechanism described by Sellwood & Binney (2002) is responsible for the largest angular momentum changes in our simulations.

We analyzed the radial surface brightness profile of the spiral galaxy NGC 7793 using HST/ACS images from the GHOSTS survey and a new HST/WFC3 image across the disk break. We used the photometry of resolved stars to select distinct populations covering a wide range of stellar ages. We found breaks in the radial profiles of all stellar populations at 280" (~5.1 kpc). Beyond this disk break, the profiles become steeper for younger populations. This same trend is seen in numerical simulations where the outer disk is formed almost entirely by radial migration. We also found that the older stars of NGC 7793 extend significantly farther than the underlying HI disk. They are thus unlikely to have formed entirely at their current radii, unless the gas disk was substantially larger in the past. These observations thus provide evidence for substantial stellar radial migration in late-type disks.

We use the Mitchell Spectrograph (formerly VIRUS-P) to observe 12 nearby disk galaxies. We successfully measure ages in the outer disk in six systems. In three cases (NGC 2684, NGC 6155, and NGC 7437), we find that a downward break in the disk surface brightness profile corresponds with a change in the dominant stellar population with the interior being dominated by active star formation and the exterior having older stellar populations that are best-fit with star formation histories that decline with time. The observed increase in average stellar ages beyond a profile break is similar to theoretical models that predict surface brightness breaks are caused by stellar migration, with the outer disk being populated from scattered old interior stars. In three more cases (IC 1132, NGC 4904, and NGC 6691), we find no significant change in the stellar population as one crosses the break radius. In these galaxies, both the inner and outer disks are dominated by active star formation and younger stellar populations. While radial migration can contribute to the stellar populations beyond the break, it appears more than one mechanism is required to explain all of our observed stellar profile breaks.

Double-barred galaxies are common in the local Universe, with approximately one third of barred spirals hosting an smaller, inner bar. Nested bars have been proposed as a mechanism for transporting gas to the very central regions of the galaxy, trigger star formation and contribute to the growth of the bulge. To test this idea, we perform for the first time a detailed analysis of the photometry, kinematics and stellar populations of a double-barred galaxy: NGC 357. We find that this galaxy is either hosting a pseudobulge or a classical bulge together with an inner disc. We compare the relative mean luminosity-weighted age, metallicity and alpha-enhancement between the (pseudo)bulge, inner bar and outer bar, finding that the three structures are nearly coeval and old. Moreover, the bulge and inner bar present the same metallicity and overabundance, whereas the outer bar tends to be less metal-rich and more alpha-enhanced. These results point out that, rather than the classical secular scenario in which gas and star formation play a major role, the redistribution of the existing stars is driving the formation of the inner structures.

The phase space coordinates of individual halo stars obtained by Galactic surveys enable the computation of their full 3-dimensional orbits. Spectral analysis of halo orbits can be used to construct "frequency maps" which provide a compact representation of the 6-dimensional phase space distribution function. Frequency maps identify important major orbit families, and the orbital abundances reflect the shape and orientation of the dark matter halo relative to the disk. We apply spectral analysis to halo orbits in a series of controlled simulations of disk galaxies. Although the shape of the simulated halo varies with radius, frequency maps of local samples of halo orbits confined to the inner halo contain most of the information about the global shape of the halo and its major orbit families. Quiescent or adiabatic disk formation results in significant trapping of halo orbits in resonant orbit families (i.e. orbits with commensurable frequencies). If a good estimate of the Galactic potential in the inner halo (within ~50kpc) is available, the appearance of strong, stable resonances in frequency maps of halo orbits will allow us to determine the degree of resonant trapping induced by the disk potential. The locations and strengths of these resonant families are determined both by the global shape of the halo and its distribution function. Identification of such resonances in the Milky Way's stellar halo would therefore provide evidence of an extended period of adiabatic disk growth. If the Galactic potential is not known exactly, a measure of the diffusion rate of large sample of 10^4 halo orbits can help distinguish between the true potential and an incorrect potential. The orbital spectral analysis methods described in this paper provide a strong complementarity to existing methods for constraining the potential of the Milky Way halo and its stellar distribution function (ABRIDGED).

We present observations and dynamical models of the stellar nuclear clusters (NCs) at the centres of NGC 4244 and M33. We then compare these to an extensive set of simulations testing the importance of purely stellar dynamical mergers on the formation and growth of NCs. Mergers of star clusters are able to produce a wide variety of observed properties, including densities, structural scaling relations, shapes (including the presence of young discs) and even rapid rotation. Nonetheless, difficulties remain, most notably that the second order kinematic moment V_rms = (V^2 + sigma^2)^(1/2) of the models is too centrally peaked to match observations. This can be remedied by the merger of star clusters onto a pre-existing nuclear disc, but the line-of-sight velocity V is still more slowly rising than in NGC 4244. Our results therefore suggest that purely stellar dynamical mergers cannot form NCs, and that gas dissipation is a necessary ingredient for at least ~50% of a NC's mass. The negative vertical anisotropy found in NGC 4244 however requires at least 10% of the mass to be accreted as stars, since gas dissipation and in situ star formation leads to positive vertical anisotropy.

Recent theoretical work suggests that it may be common for stars in the disks of spiral galaxies to migrate radially across significant distances in the disk. Such migrations are a result of resonant scattering with spiral arms and move the guiding centers of the stars while preserving the circularities of their orbits. Migration can therefore efficiently mix stars in all parts of the Galactic disk. We are rapidly approaching an important confluence of theory and observation, where we may soon be able to uncover signatures of such processes in our own Milky Way. The resolution and robustness of the physical modeling in simulations has improved drastically, while observational datasets are increasing in depth and astrometric accuracy. Here, we discuss the results from our idealized N-body/SPH simulations of disk formation and evolution, emphasizing specifically the observational consequences of stellar migration on the solar neighborhood and the vertical structure of the disk. We demonstrate that radial mixing of stars is a crucial dynamical process that we must try to understand if we are to draw significant conclusions about our Galactic history from the properties of stars in our vicinity.

The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies.

We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmological simulation to investigate the nature of the outer regions of discs and their relevance for the disc formation process. Specifically, we focus on the phenomenon of misaligned disc components and find that the outer disc warp is a consequence of the misalignment between the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks toward the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas reaches the central disc-forming region its angular momentum is completely aligned with the spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of the inner disc, a misaligned outer disc forms comprised of newly accreted material. The inner and outer components are misaligned with each other because they respond differently to infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous angular momentum losses, but small amounts of star formation in the warp itself form a low-metallicity thick disc. We show that observations of resolved stellar populations in warped galaxies in the local universe could provide evidence for the presence of these processes and therefore indirectly reveal ongoing gas accretion and the existence of hot gas halos.

Like the majority of spiral galaxies, NGC 6155 exhibits an exponential surface brightness profile that steepens significantly at large radii. Using the VIRUS-P IFU spectrograph, we have gathered spatially resolved spectra of the system. Modifying the GANDALF spectral fitting routine for use on the complex stellar populations found in spirals, we find that the average stellar ages increase significantly beyond the profile break radius. This result is in good agreement with recent simulations that predict the outskirts of disk galaxies are populated through stellar migration. With the ability to bin multiple fibers, we are able to measure stellar population ages down to mu_V~24 mag/sq arcsec.

We examine the nuclear morphology, kinematics, and stellar populations in nearby S0 galaxy NGC 404 using a combination of adaptive optics assisted near-IR integral-field spectroscopy, optical spectroscopy, and HST imaging. These observations enable study of the NGC 404 nucleus at a level of detail possible only in the nearest galaxies. The surface brightness profile suggests the presence of three components, a bulge, a nuclear star cluster, and a central light excess within the cluster at radii <3 pc. These components have distinct kinematics with modest rotation seen in the nuclear star cluster and counter-rotation seen in the central excess. Molecular hydrogen emission traces a disk with rotation nearly orthogonal to that of the stars. The stellar populations of the three components are also distinct, with half of the mass of the nuclear star cluster having ages of ~1 Gyr (perhaps resulting from a galaxy merger), while the bulge is dominated by much older stars. Dynamical modeling of the stellar kinematics gives a total nuclear star cluster mass of 1.1x10^7 Msol. Dynamical detection of a possible intermediate mass black hole is hindered by uncertainties in the central stellar mass profile. Assuming a constant mass-to-light ratio, the stellar dynamical modeling suggests a black hole mass of <1x10^5 Msol, while the molecular hydrogen gas kinematics are best fit by a black hole with mass of 4.5x10^5 Msol. Unresolved and possibly variable dust emission in the near-infrared and AGN-like molecular hydrogen emission line ratios do suggest the presence of an accreting black hole in this nearby LINER galaxy.

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.

We present new Hubble Space Telescope (HST) observations of NGC 300 taken as part of the ACS Nearby Galaxy Survey Treasury (ANGST). Individual stars are resolved in these images down to an absolute magnitude of M_F814W = 1.0 (below the red clump). We determine the star formation history of the galaxy in 6 radial bins by comparing our observed color-magnitude diagrams (CMDs) with synthetic CMDs based on theoretical isochrones. We find that the stellar disk out to 5.4 kpc is primarily old, in contrast with the outwardly similar galaxy M33. We determine the scale length as a function of age and find evidence for inside-out growth of the stellar disk: the scale length has increased from 1.1 +/- 0.1 kpc 10 Gyr ago to 1.3 +/- 0.1 kpc at present, indicating a buildup in the fraction of young stars at larger radii. As the scale length of M33 has recently been shown to have increased much more dramatically with time, our results demonstrate that two galaxies with similar sizes and morphologies can have very different histories. With an N-body simulation of a galaxy designed to be similar to NGC 300, we determine that the effects of radial migration should be minimal. We trace the metallicity gradient as a function of time and find a present day metallicity gradient consistent with that seen in previous studies. Consistent results are obtained from archival images covering the same radial extent but differing in placement and filter combination.

We investigate the effect of dust on the observed rotation rate of a stellar bar as measured by application of the Tremaine & Weinberg (TW) method. As the method requires that the tracer satisfies the continuity equation, it has been applied largely to early-type barred galaxies. We show using numerical simulations of barred galaxies that dust attenuation factors typically found in these systems change the observed bar pattern speed by 20-40%. We also address the effect of star formation on the TW method and find that it does not change the results significantly. This suggests that applications of the TW method can be extended to include barred galaxies covering the full range of Hubble type.

We present high resolution absorption-line spectroscopy of 3 face-on galaxies, NGC 98, NGC 600, and NGC 1703 with the aim of searching for box/peanut (B/P)-shaped bulges. These observations test and confirm the prediction of Debattista et al. (2005) that face-on B/P-shaped bulges can be recognized by a double minimum in the profile of the fourth-order Gauss-Hermite moment h_4. In NGC 1703, which is an unbarred control galaxy, we found no evidence of a B/P bulge. In NGC 98, a clear double minimum in h_4 is present along the major axis of the bar and before the end of the bar, as predicted. In contrast, in NGC 600, which is also a barred galaxy but lacks a substantial bulge, we do not find a significant B/P shape.

We summarize the recent results of our survey of the nearest nuclear star clusters. The purpose of the survey is to understand nuclear star cluster formation mechanisms and constrain the presence of black holes using adaptive optics assisted integral field spectroscopy, optical spectroscopy, and HST imaging in 13 galaxies within 5 Mpc. We discuss the formation history of the nuclear star cluster and possible detection of an intermediate mass BH in NGC 404, the nearest S0 galaxy.

We summarize recent developments in the use of spectral methods for analyzing large numbers of orbits in N-body simulations to obtain insights into the global phase space structure of dark matter halos. The fundamental frequencies of oscillation of orbits can be used to understand the physical mechanism by which the shapes of dark matter halos evolve in response to the growth of central baryonic components. Halos change shape primarily because individual orbits change their shapes adiabatically in response to the growth of a baryonic component, with those at small radii become preferentially rounder. Chaotic scattering of orbits occurs only when the central point mass is very compact and is equally effective for centrophobic long-axis tube orbits as it is for centrophilic box orbits.

Using spectral methods, we analyse the orbital structure of prolate/triaxial dark matter (DM) halos in N-body simulations to understand the processes that drive the evolution of shapes of DM halos and elliptical galaxies in which central masses are grown. A longstanding issue is whether the change in the shapes of DM halos is the result of chaotic scattering of box orbits, or whether they change shape adiabatically in response to the evolving galactic potential. We use orbital frequencies to classify orbits, to quantify orbital shapes, and to identify resonant orbits and chaotic orbits. The frequency-based method overcomes the limitations of Lyapunov exponents which are sensitive to numerical discreteness effects. Regardless of the distribution of the baryonic component, the shape of a DM halo changes primarily due to changes in the shapes of individual orbits within a given family. Orbits with small pericentric radii are more likely to change both their orbital type and shape than orbits with large pericentric radii. Whether the evolution is regular (and reversible) or chaotic (and irreversible), depends primarily on the radial distribution of the baryonic component. The growth of an extended baryonic component of any shape results in a regular rather than chaotic change in orbital populations. In contrast the growth of a massive and compact central component results in chaotic scattering of a significant fraction of both box and long-axis tube orbits. The growth of a disk causes a significant fraction of halo particles to become trapped by major global orbital resonances. Despite the fact that shape of a DM halo is always quite oblate following the growth of a central baryonic component, a significant fraction of its orbit population has characteristics of its triaxial or prolate progenitor (ABRIDGED).

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

We demonstrate that the chaotic nature of N-body systems can lead to macroscopic variations in the evolution of collisionless simulations containing rotationally supported discs. The unavoidable stochasticity that afflicts all simulations generally causes mild differences between the evolution of similar models but, in order to illustrate that this is not always true, we present a case that shows extreme bimodal divergence. The divergent behaviour occurs in two different types of code and is independent of all numerical parameters. We identify and give explicit illustrations of several sources of stochasticity, and also show that macroscopic variations in the evolution can originate from differences at the round-off error level. We obtain somewhat more consistent results from simulations in which the halo is set up with great care compared with those started from more approximate equilibria, but we have been unable to eliminate diverging behaviour entirely because the main sources of stochasticity are intrinsic to the disc. We show that the divergence is only temporary and that halo friction is merely delayed, for a substantial time in some cases. We argue that the delays are unlikely to arise in real galaxies, and that our results do not affect dynamical friction constraints on halo density. Stochastic variations in the evolution are inevitable in all simulations of disc-halo systems, irrespective of how they were created, although their effect is generally far less extreme than we find here. The possibility of divergent behaviour complicates comparison of results from different workers.

(Abridged) We perform dissipationless N-body simulations to elucidate the dynamical response of thin disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-sized CDM halos to derive the properties of substructure and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially-thin, fully-formed MW type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. We find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z~1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of observed antitruncated disks; (4) lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We conclude that substructure-disk encounters of the kind expected in the LCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.

Predicting the local flux of dark matter particles is vital for dark matter direct detection experiments. To date, such predictions have been based on simulations that model the dark matter alone. Here we include the influence of the baryonic matter for the first time. We use two different approaches. Firstly, we use dark matter only simulations to estimate the expected merger history for a Milky Way mass galaxy, and then add a thin stellar disc to measure its effect. Secondly, we use three cosmological hydrodynamic simulations of Milky Way mass galaxies. In both cases, we find that a stellar/gas disc at high redshift (z~1) causes merging satellites to be preferentially dragged towards the disc plane. This results in an accreted dark matter disc that contributes ~0.25 - 1 times the non-rotating halo density at the solar position. An associated thick stellar disc forms with the dark disc and shares a similar velocity distribution. If these accreted stars can be separated from those that formed in situ, future astronomical surveys will be able to infer the properties of the dark disc from these stars. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, increases the likelihood of direct detection at low recoil energy, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined (see contribution from T. Bruch this volume).

We present observations of the nuclear star cluster in the nearby edge-on spiral galaxy NGC 4244 using the Gemini Near-Infrared Integral Field Spectrograph (NIFS) with laser guide star adaptive optics. From a previous study of edge-on galaxies, this nuclear star cluster was found to be one of a sample of clusters that appear flattened along the plane of their host galaxies disks. Such clusters show evidence for multiple morphological components, with younger/bluer disk components and older/redder spheroidal components. Our new observations of NGC 4244 show clear rotation of 30 km/sec within the central 10 pc (0.5") of the cluster. The central velocity dispersion is found to be 28+/-2 km/sec. The multiple stellar populations inferred from the optical colors and spectra are seen as variations in the CO line strength in the NIFS spectra. The rotation is clearly detected even in the older, more spheroidal stellar component. We discuss evidence for similar structures and kinematics in the nuclear star clusters of other galaxies including M33 and the Milky Way. Our observations support two possible formation mechanisms: (1) episodic accretion of gas from the disk directly onto the nuclear star cluster, or (2) episodic accretion of young star clusters formed in the central part of the galaxy due to dynamical friction.

Stars in disks of spiral galaxies are usually assumed to remain roughly at their birth radii. This assumption is built into decades of modelling of the evolution of stellar populations in our own Galaxy and in external systems. We present results from self-consistent high-resolution $N$-body + Smooth Particle Hydrodynamics simulations of disk formation, in which stars migrate across significant galactocentric distances due to resonant scattering with transient spiral arms, while preserving their circular orbits. We investigate the implications of such migrations for observed stellar populations. Radial migration provides an explanation for the observed flatness and spread in the age-metallicity relation and the relative lack of metal poor stars in the solar neighborhood. The presence of radial migration also prompts rethinking of interpretations of extra-galactic stellar population data, especially for determinations of star formation histories.

The outer disks of galaxies present a unique laboratory for studying the process of disk formation. A considerable fraction of observed disks exhibit a break in their surface brightness profiles. The ubiquity of these features points to a crucial aspect of disk formation which must be explained. Recent theoretical work suggests that such breaks are related to significant amounts of radial migration. We discuss the current observational evidence which supports this picture.

In a LCDM cosmology, the Milky Way accretes satellites into the stellar disc. We use cosmological simulations to assess the frequency of near disc plane and higher inclination accretion events, and collisionless simulations of satellite mergers to quantify the final state of the accreted material and the effect on the thin disc. On average, a Milky Way-sized galaxy has 3 subhalos with vmax>80km/s; 7 with vmax>60km/s; and 15 with vmax>40km/s merge at redshift z>1. Assuming isotropic accretion, a third of these merge at an impact angle <20 degrees and are dragged into the disc plane by dynamical friction. Their accreted stars and dark matter settle into a thick disc. The stellar thick disc qualitatively reproduces the observed thick disc at the solar neighbourhood, but is less massive by a factor ~2-10. The dark matter disc contributes 0.25-1 times the halo density at the solar position. Although not likely to be dynamically interesting, the dark disc has important implications for the direct detection of dark matter because of its low velocity with respect to the Earth. Higher inclination encounters (>20 degrees) are twice as likely as low inclination ones. These lead to structures that closely resemble the recently discovered inner/outer stellar halos. They also do more damage to the Milky Way stellar disc creating a more pronounced flare, and warp; both long-lived and consistent with current observations. The most massive mergers (vmax>80km/s) heat the thin disc enough to produce a thick disc. These heated thin disc stars are essential for obtaining a thick disc as massive as that seen in the Milky Way; they likely comprise some ~50-90% of the thick disc stars. The Milky Way thin disc must reform from fresh gas after z=1 [abridged].

The presence of black holes (BHs) at the centers of dwarf elliptical galaxies (dEs) has been argued both theoretically and observationally. Using archival HST/WFPC2 data, we found the Virgo cluster dwarf elliptical galaxy VCC128 to harbor a binary nucleus, a feature that is usually interpreted as the observable signature of a stellar disk orbiting a central massive black hole. Debattista et al. 2006 estimated its mass M sim 6 10^6 - 5 10^7 Msun. One of the most robust means of verifying the existence of a BH is radio continuum and/or X-ray emission, however because of the deficiency of gas in dEs, radio continuum emission forms the best option here. We have tried to detect the X-band radio emission coming from the putative black hole in VCC128 when it accretes gas from the surrounding ISM. While we made a positive 4 sigma detection of a point source 4.63'' south-west of the binary nucleus, no statistically significant evidence for emission associated with the nuclei themselves was detected. This implies either that VCC128 has no massive central black hole, which makes the nature of the binary nucleus hard to explain, or, if it has a central black hole, that the physical conditions of the ISM (predominantly its density and temperature) and/or of the surrounding accretion disk do not allow for efficient gas accretion onto the black hole, making the quiescent black hole very hard to detect at radio wavelengths.

Cold dark matter cosmogony predicts triaxial dark matter halos, whereas observations find quite round halos. This is most likely due to the condensation of baryons leading to rounder halos. We examine the halo phase space distribution basis for such shape changes. Triaxial halos are supported by box orbits, which pass arbitrarily close to the density center. The decrease in triaxiality caused by baryons is thought to be due to the scattering of these orbits. We test this hypothesis with simulations of disks grown inside triaxial halos. After the disks are grown we check whether the phase space structure has changed by evaporating the disks and comparing the initial and final states. While the halos are substantially rounder when the disk is at full mass, their final shape after the disk is evaporated is not much different from the initial. Likewise, the halo becomes (more) radially anisotropic when the disk is grown, but the final anisotropy is consistent with the initial. Only if the baryons are unreasonably compact or massive does the halo change irreversibly. We show that the character of individual orbits is not generally changed by the growing mass. Thus the central condensation of baryons does not destroy enough box orbits to cause the shape change. Rather, box orbits merely become rounder along with the global potential. However, if angular momentum is transferred to the halo, either via satellites or via bars, a large irreversible change in the halo distribution occurs. The ability of satellites to alter the phase space distribution of the halo is of particular concern to galaxy formation simulations since halo triaxiality can profoundly influence the evolution of disks.

We present the results of high resolution absorption-line spectroscopy of 3 face-on galaxies, NGC 98, NGC 600, and NGC 1703 with the aim of searching for box/peanut (B/P)-shaped bulges. These observations test and confirm, for the first, time the prediction that face-on B/P-shaped bulges can be recognized by two minima in the profile along the bar's major axis of the fourth Gauss-Hermite moment h_4 of the line-of-sight velocity distribution (LOSVD). In NGC 98, a clear double minimum in h_4 is present along the major axis of the bar and before the end of the bar, as predicted. In contrast, in NGC 600, which is also a barred galaxy but lacks a substantial bulge, we do not find any significant kinematic signature for a B/P-shaped bulge. In NGC 1703, which is an unbarred control galaxy, we found no evidence of a B/P bulge. We also show directly that the LOSVD is broader at the location of the h_4 minimum in NGC 98 than elsewhere. This more direct method avoids possible artifacts associated with the degeneracy between the measurement of line-of-sight velocity dispersion and h_4.

Recent results from the Planetary Nebula Spectrograph (PN.S) survey have revealed a rapidly falling velocity dispersion profile in the nearby elliptical galaxy NGC 3379, casting doubts on whether this intermediate-luminosity galaxy has the kind of dark matter halo expected in LambdaCDM cosmology. We present a detailed dynamical study of this galaxy, combining long-slit spectroscopy, SAURON integral-field data, and PN.S velocities, reaching to more than seven effective radii (R_e). We construct spherical and axisymmetric dynamical models for these data with the flexible made-to-measure NMAGIC code, in a sequence of gravitational potentials with varying dark halo mass. We find that the data are consistent both with near-isotropic spherical systems dominated by the stellar mass, and with models in massive halos with strongly radially anisotropic outer parts (beta >~ 0.8 at 7R_e). Formal likelihood limits would exclude (at 1 sigma) the model with stars only, as well as halo models with v_circ(7R_e) >~ 250 km/s. A sequence of more realistic axisymmetric models of different inclinations and a small number of triaxial tests confirm the spherical results. All valid models fitting all the data are dynamically stable over Gyrs, including the most anisotropic ones. NGC 3379 may well have a dark matter halo as predicted by recent merger models within LambdaCDM cosmology, provided its outer envelope is strongly radially anisotropic. (abridged)

Lopsidedness is common in spiral galaxies. Often, there is no obvious external cause, such as an interaction with a nearby galaxy, for such features. Alternatively, the lopsidedness may have an internal cause, such as a dynamical instability. In order to explore this idea, we have developed a computer code that searches for self-consistent perturbations in razor-thin disc galaxies and performed a thorough mode-analysis of a suite of dynamical models for disc galaxies embedded in an inert dark-matter halo with varying amounts of rotation and radial anisotropy. Models with two equal-mass counter-rotating discs and fully rotating models both show growing lopsided modes. For the counter-rotating models, this is the well-known counter-rotating instability, becoming weaker as the net rotation increases. The m=1 mode of the maximally rotating models, on the other hand, becomes stronger with increasing net rotation. This rotating m=1 mode is reminiscent of the eccentricity instability in near-Keplerian discs. To unravel the physical origin of these two different m=1 instabilities, we studied the individual stellar orbits in the perturbed potential and found that the presence of the perturbation gives rise to a very rich orbital behaviour. In the linear regime, both instabilities are supported by aligned loop orbits. In the non-linear regime, other orbit families exist that can help support the modes. In terms of density waves, the counter-rotating m=1 mode is due to a purely growing Jeans-type instability. The rotating m=1 mode, on the other hand, grows as a result of the swing amplifier working inside the resonance cavity that extends from the disc center out to the radius where non-rotating waves are stabilized by the model's outwardly rising Q-profile.

We have performed a high mass and force resolution simulation of an idealized galaxy forming from dissipational collapse of gas embedded in a spherical dark matter halo. The simulation includes star formation and effects of stellar feedback. In our simulation a stellar disk forms with a surface density profile consisting of an inner exponential breaking to a steeper outer exponential. The break forms early on and persists throughout the evolution, moving outward as more gas is able to cool and add mass to the disk. The parameters of the break are in excellent agreement with observations. The break corresponds with a rapid drop in the star formation rate associated with a drop in the cooled gas surface density, but the outer exponential is populated by stars that were scattered outward on nearly circular orbits from the inner disk by spiral arms. The resulting profile and its associated break are therefore a consequence of the interplay between a radial star formation cutoff and redistribution of stellar mass by secular processes. A consequence of such evolution is a sharp change in the radial mean stellar age profile at the break radius.

We present a dynamical study of NGC 4697, an almost edge-on, intermediate-luminosity, E4 elliptical galaxy, combining new surface brightness photometry, new as well as published long-slit absorption line kinematic data, and published planetary nebulae (PNe) velocity data. The combined kinematic data set extends out to ~= 5' ~= 4.5 R_e and allows us to probe the galaxy's outer halo. For the first time, we model such a dataset with the new and flexible Chi^2-made-to-measure particle code NMAGIC. We extend NMAGIC to include seeing effects, introduce an efficient scheme to estimate the mass-to-light ratio, and incorporate a maximum likelihood technique to account for discrete velocity measurements. For modelling the PNe kinematics we use line-of-sight velocities and velocity dispersions computed on two different spatial grids, and we also use the individual velocity measurements with the likelihood method, in order to make sure that our results are not biased by the way we treat the PNe measurements. We generate axisymmetric self-consistent models as well models including various dark matter halos. These models fit all the mean velocity and velocity dispersion data with Chi^2/N<1, both in the case with only luminous matter and in potentials including quite massive halos. The likelihood analysis together with the velocity histograms suggest that models with low density halos such that the circular velocity v_c ~< 200 km/s at R_e are not consistent with the data. A range of massive halos with v_c ~= 250 km/s at 5R_e fit the PN data best. To derive stronger results would require PN velocities at even larger radii. The best fitting models are slightly radially anisotropic; the anisotropy parameter beta ~= 0.3 at the center, increasing to beta ~= 0.5 at radii >~ 2R_e.

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.

We present new estimates of ages and metallicities, based on FORS/VLT optical (4400-5500A) spectroscopy, of 16 dwarf elliptical galaxies (dE's) in the Fornax Cluster and in Southern Groups. These dE's are more metal-rich and younger than previous estimates based on narrow-band photometry and low-resolution spectroscopy. For our sample we find a mean metallicity [Z/H] = -0.33 dex and mean age 3.5 Gyr, consistent with similar samples of dE's in other environments (Local Group, Virgo). Three dE's in our sample show emission lines and very young ages. This suggests that some dE's formed stars until a very recent epoch and were self-enriched by a long star formation history. Previous observations of large near-infrared (~8500A) Ca II absorption strengths in these dE's are in good agreement with the new metallicity estimates, solving part of the so-called Calcium puzzle.

Red clump giants in the Galactic bulge are approximate standard candles and hence they can be used as distance indicators. We compute the proper motion dispersions of RCG stars in the Galactic bulge using the proper motion catalogue from the second phase of the Optical Gravitational Microlensing Experiment (OGLE-II, Sumi et al. 2004) for 45 fields. The proper motion dispersions are measured to a few per cent accuracy due to the large number of stars in the fields. The observational sample is comprised of 577736 stars. These observed data are compared to a state-of-the-art particle simulation of the Galactic bulge region. The predictions are in rough agreement with observations, but appear to be too anisotropic in the velocity ellipsoid. We note that there is significant field-to-field variation in the observed proper motion dispersions. This could either be a real feature, or due to some unknown systematic effect.

We investigate the effect of dust on the observed rotation rate of a stellar bar. The only direct way to measure this quantity relies on the Tremaine & Weinberg method which requires that the tracer satisfies the continuity equation. Thus it has been applied largely to early-type barred galaxies. We show using numerical simulations of barred galaxies that dust attenuation factors typically found in these systems change the observed bar pattern speed by 20-40 percent. We also address the effect of star formation on the TW method and find that it does not change the results significantly. The results presented here suggest that applications of the TW method can be extended to include barred galaxies covering the full range of Hubble type.

We describe a made-to-measure algorithm for constructing N-particle models of stellar systems from observational data (Chi-Squared-M2M), extending earlier ideas by Syer and Tremaine. The algorithm properly accounts for observational errors, is flexible, and can be applied to various systems and geometries. We implement this algorithm in a parallel code NMAGIC and carry out a sequence of tests to illustrate its power and performance: (i) We reconstruct an isotropic Hernquist model from density moments and projected kinematics and recover the correct differential energy distribution and intrinsic kinematics. (ii) We build a self-consistent oblate three-integral maximum rotator model and compare how the distribution function is recovered from integral field and slit kinematic data. (iii) We create a non-rotating and a figure rotating triaxial stellar particle model, reproduce the projected kinematics of the figure rotating system by a non-rotating system of the same intrinsic shape, and illustrate the signature of pattern rotation in this model. From these tests we comment on the dependence of the results from Chi-Squared-M2M on the initial model, the geometry, and the amount of available data.

We present surface photometry and stellar kinematics of NGC 4431, a barred dwarf galaxy in the Virgo cluster undergoing a tidal interaction with one of its neighbors, NGC 4436. We measured its bar pattern speed using the Tremaine-Weinberg method, and derived the ratio of the corotation radius, D_L, to the bar semi-major axis, a_B. We found D_L/a_B=0.6^{+1.2}_{-0.4} at 99% confidence level. Albeit with large uncertainty, the probability that the bar ends close to its corotation radius (i.e., 1.0 \leq D_L/a_B \leq 1.4) is about twice as likely as that the bar is much shorter than corotation radius (i.e., D_L/a_B > 1.4).

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.

{\it Hubble Space Telescope} (\hst) Wide Field Planetary Camera 2 ({\it WFPC2}) images of the Virgo Cluster dwarf elliptical galaxy VCC 128 reveal an apparently double nucleus. The two components, which are separated by $\sim 32$ pc in projection, have the same magnitude and color. We present a spectrum of this double nucleus and show that it is inconsistent with one or both components being emission-line background objects or foreground stars. The most likely interpretation is that, as suggested by \citet{lauer_etal_96} for the double nucleus in NGC 4486B, we are seeing a nuclear disk surrounding a supermassive black hole. This is only the second time an early-type dwarf (dE/dSph) galaxy has been suggested to host a SMBH.

We find 9 nuclear cluster candidates in a sample of 14 edge-on, late-type galaxies observed with HST/ACS. These clusters have magnitudes (M_I ~ -11) and sizes (r_eff ~ 3pc) similar to those found in previous studies of face-on, late-type spirals and dE galaxies. However, three of the nuclear clusters are significantly flattened and show evidence for multiple, coincident structural components. The elongations of these three clusters are aligned to within 10 degrees of the galaxies' major axes. Structurally, the flattened clusters are well fit by a combination of a spheroid and a disk or ring. The nuclear cluster disks/rings have F606W-F814W (~V-I) colors 0.3-0.6 magnitudes bluer than the spheroid components, suggesting that the stars in these components have ages < 1 Gyr. In NGC 4244, the nearest of the nuclear clusters, we further constrain the stellar populations and provide a lower limit on the dynamical mass via spectroscopy. We also present tentative evidence that another of the nuclear clusters (in NGC 4206) may also host a supermassive black hole. Based on our observational results we propose an in situ formation mechanism for nuclear clusters in which stars form episodically in compact nuclear disks, and then lose angular momentum or heat vertically to form an older spheroidal structure. We estimate the period between star formation episodes to be 0.5 Gyr and discuss possible mechanisms for tranforming the disk-like components into spheroids. We also note the connection between our objects and massive globular clusters (e.g. $\omega$ Cen), UCDs, and SMBHs. (Abridged)

We present a comprehensive series of $N$-body as well as $N$-body + SPH simulations to study the secular evolution of the structure of disk galaxies. Our simulations are organized in a hierarchy of increasing complexity, ranging from rigid-halo collisionless simulations to fully live simulations with gas and star formation. Comparisons between the different types of simulations allow us to isolate the role of various physical mechanisms. We focus on the evolution of systems expected in a LCDM universe. Our goal is to examine which structural properties of disk galaxies may result from secular evolution rather than from direct hierarchical assembly. In the vertical direction, we find that various mechanisms can lead to heating. The strongest heating occurs during the vertical buckling instability of a bar. Among the consequences of this instability is the formation of peanut-shaped bulges which produce clear kinematic signatures when observed face-on. We find that bars are robust structures that are not destroyed by buckling. They can be destroyed instead by a central mass concentration but we find that this mass needs to be a large fraction of the total mass of the disk. We then study the evolution of stellar surface density profiles showing how angular momentum redistribution leads to increasing central densities and disk scale lengths and to profile breaks at large radii. The breaks in these simulations are in excellent agreement with observed breaks, even when the evolution is purely collisionless. Disk scale-lengths increase even when the total disk angular momentum is conserved; thus mapping halo angular momenta to scale-lengths is non-trivial. [Abridged]

From observations of edge-on, late-type galaxies, we present morphological evidence that some nuclear star clusters have experienced in situ star formation. We find three nuclear clusters that, viewed from the edge-on perspective, have both a compact disk-like component and a spheroidal component. In each cluster, the disk components are closely aligned with the major axis of the host galaxy and have bluer colors than the spheroidal components. We spectroscopically verify that one of the observed multiple component clusters has multiple generations of stars. These observations lead us to suggest a formation mechanism for nuclear star clusters, in which stars episodically form in compact nuclear disks, and then lose angular momentum, eventually forming an older spheroid.

Despite the increasing number of studies of barred galaxies at intermediate and high redshifts, double-barred (S2B) systems have only been identified in the nearby (z<0.04) universe thus far. In this feasibility study we demonstrate that the detection and analysis of S2Bs is possible at intermediate redshifts (0.1 < z < 0.5) with the exquisite resolution of the Hubble Space Telescope Advanced Camera for Surveys (HST/ACS). We identify barred galaxies in the HST/ACS data of the Great Observatories Origins Deep Survey (GOODS) using a novel method. The radial profile of the Gini coefficient -- a model-independent structure parameter -- is able to detect bars in early-type galaxies that are large enough that they might host an inner bar of sufficient angular size. Using this method and subsequent examination with unsharp masks and ellipse fits we identified the two most distant S2Bs currently known (at redshifts z=0.103 and z=0.148). We investigate the underlying stellar populations of these two galaxies through a detailed colour analysis, in order to demonstrate the analysis that could be performed on a future sample of intermediate-redshift S2Bs. We also identify two S2Bs and five S2B candidates in the HST/ACS data of the Cosmic Evolution Survey (COSMOS). Our detections of distant S2Bs show that deep surveys like GOODS and COSMOS have the potential to push the limit for S2B detection and analysis out by a factor of ten in redshift and lookback time (z=0.5, t=5Gyr) compared to the previously known S2Bs. This in turn would provide new insight into the formation of these objects.

I review recent progress from $N$-body simulations in our understanding of the secular evolution of isolated disk galaxies. I describe some of the recent controversies in the field which have been commonly attributed to numerics. The numerical methods used are widely used in computational astronomy and the problems encountered are therefore of wider interest.

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

A bar rotating in a pressure-supported halo generally loses angular momentum and slows down due to dynamical friction. Valenzuela & Klypin report a counter-example of a bar that rotates in a dense halo with little friction for several Gyr, and argue that their result invalidates the claim by Debattista & Sellwood that fast bars in real galaxies require a low halo density. We show that it is possible for friction to cease for a while should the pattern speed of the bar fluctuate upward. The reduced friction is due to an anomalous gradient in the phase-space density of particles at the principal resonance created by the earlier evolution. The result obtained by Valenzuela & Klypin is probably an artifact of their adaptive mesh refinement method, but anyway could not persist in a real galaxy. The conclusion by Debattista & Sellwood still stands.

We present a kinematic diagnostic for peanut-shaped bulges in nearly face-on galaxies. The face-on view provides a novel perspective on peanuts which would allow study of their relation to bars and disks in greater detail than hitherto possible. The diagnostic is based on the fact that peanut shapes are associated with a flat density distribution in the vertical direction. We show that the kinematic signature corresponding to such a distribution is a minimum in the fourth-order Gauss-Hermite moment $s_4$. We demonstrate our method on $N$-body simulations of varying peanut strength, showing that strong peanuts can be recognized to inclinations $i \simeq 30\degrees$, regardless of the strength of the bar. We also consider compound systems in which a bulge is present in the initial conditions as may happen if bulges form at high redshift through mergers. We show that in this case, because the vertical structure of the bulge is not derived from that of the disk, that the signature of a peanut in $s_4$ is weakened. Thus the same kinematic signature of peanuts can be used to explore bulge formation mechanisms. The observational requirements for identifying peanuts with this method are challenging, but feasible.

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high-resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. The galaxy models are consistent with the LCDM paradigm of structure formation and the simulations include the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation of a close SMBH pair at the center of the remnant with separations limited solely by the adopted force resolution of ~ 100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can be entirely disrupted before the merger is completed leaving its SMBH wandering at a distance too far from the center of the remnant for the formation of a close pair. In contrast, we show that gas cooling facilitates the pairing process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that merging disk galaxies constructed to obey the M(BH)-sigma relation, move relative to it depending on whether they undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless simulations reveal that remnants tend to move away from the mean relation highlighting the role of gas-poor mergers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary fuel to the SMBHs and allows the merger remnants to satisfy the relation.

The recent discovery of hidden non-axysimmetric and disk-like structures in bright Virgo dwarf elliptical and lenticular galaxies (dE/dSph/dS0) indicates that they may have late-type progenitors. Using N-body simulations we follow the evolution of disk galaxies within a LambdaCDM cluster simulated with 10e7 particles, where the hierarchical growth and galaxy harassment are modeled self-consistently. Most of the galaxies undergo significant morphological transformation, even at the outskirts of the cluster, and move through the Hubble sequence from late type disks to dwarf spheroidals. None of the disks are completely destroyed therefore they can not be the progenitors of ultra compact dwarf galaxies (UCDs). The time evolution of the simulated galaxies is compared with unsharp-masked images obtained from VLT data and the projected kinematics of our models with the latest high resolution spectroscopic studies from the Keck and Palomar telescopes.

n many nucleated dwarf elliptical galaxies (dE,N's), the nucleus is offset by a significant fraction of the scale radius with respect to the center of the outer isophotes. Using a high-resolution N-body simulation, we demonstrate that the nucleus can be driven off-center by the m=1 counterstreaming instability, which is strong in flattened stellar systems with zero rotation. The model develops a nuclear offset on the order of 30% of the exponential scale length. We compare our numerical results with the photometry and kinematics of FCC 046, a Fornax Cluster dE,N with a nucleus offset by 1.2" we find good agreement between the model and FCC 046. We also discuss mechanisms that may cause counterrotation in dE,N's and conclude that the destruction of box orbits in an initially triaxial galaxy is the most promising.

We use high resolution collisionless $N$-body simulations to study the secular evolution of disk galaxies and in particular the final properties of disks that suffer a bar and perhaps a bar-buckling instability. Although we find that bars are not destroyed by the buckling instability, when we decompose the radial density profiles of the secularly-evolved disks into inner S\'ersic and outer exponential components, for favorable viewing angles, the resulting structural parameters, scaling relations and global kinematics of the bar components are in good agreement with those obtained for bulges of late-type galaxies. Round bulges may require a different formation channel or dissipational processes.

I show that Tremaine-Weinberg (TW) measurements of bar pattern speeds are sensitive to errors in the position angle of the disc, PA_disc. I use an $N$-body experiment to measure these errors; for typical random PA_disc errors, the resulting scatter in the measured values of the dimensionless bar speed parameter $R$ (defined as the ratio of the corotation radius to the bar semi-major axis) is of order the observational. I also consider how the systematic PA_disc errors produced by disc ellipticities affect TW measurements. The scatter produced by these errors may be significant, depending on the ellipticity distribution. Conversely, by using the sample of TW observations, I find that an upper limit of the typical disc (density) ellipticity is 0.07 at the 90 per cent confidence level, which is in good agreement with previous measurements. Taken together, the random and systematic scatter suggest that the intrinsic distribution of $R$ of gas-poor early-type barred galaxies may be as narrow as that of the gas-rich later-types.

We present Fabry-Perot absorption-line spectroscopy of the SB0 galaxy NGC 7079. This is the first use of Fabry-Perot techniques to measure the two-dimensional stellar kinematics of an early-type disk galaxy. We scan the infrared CaII line using the Rutgers Fabry-Perot (RFP), to obtain kinematic data extending to $I$-band surface brightness $\mu_I \simeq 21$ mag./arcsec^-2, in a field of radius $\sim 40\arcsec$. The kinematic data, consisting of line-of-sight velocities and velocity dispersions, are in good agreement with data obtained along the major axis of the disk with standard slit spectroscopy. Comparison of the exposure times required for slit and RFP spectroscopy to reach the same limiting magnitude shows that the RFP is significantly more efficient for mapping absorption-line galaxy kinematics. We use the velocity data, together with our own deep broad-band photometry,to measure the bar pattern speed, $\Omega_p$, of NGC 7079 with the model-independent Tremaine-Weinberg (TW) method. We find $\Omega_p = 8.4 \pm 0.2$ km/s/arcsec; this is the best-constrained pattern speed ever measured for a bar using the TW method. From the rotation curve, corrected for asymmetric drift, we calculate the co-rotation radius and find that the bar ends just inside this radius. The two-dimensional character of these data allow us to show that the TW method is sensitive to errors in the position angle (PA) of the disk. For example, a PA error of $2\degrees$ can give errors $\sim \pm 25%$ in $\Omega_p$.

We build a stellar-dynamical model of the Milky Way barred bulge and disk, using a newly implemented adaptive particle method. The underlying mass model has been previously shown to match the Galactic near-infrared surface brightness as well as gas-kinematic observations. Here we show that the new stellar-dynamical model also matches the observed stellar kinematics in several bulge fields, and that its distribution of microlensing event timescales reproduces the observed timescale distribution of the {\it MACHO} experiment with a reasonable stellar mass function. The model is therefore an excellent basis for further studies of the Milky Way. We also predict the observational consequences of this mass function for parallax shifted events.

We measured the bar pattern speed in a sample of 7 SB0 galaxies using the Tremaine-Weinberg method. This represents the largest sample of galaxies for which bar pattern speed has been measured this way. All the observed bars are as rapidly rotating as they can be. We compared this result with recent high-resolution N-body simulations of bars in cosmologically-motivated dark matter halos, and conclude that these bars are not located inside centrally concentrated halos.

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

We present surface photometry and stellar kinematics of a sample of 5 SB0 galaxies: ESO 139-G009, IC 874, NGC 1308, NGC 1440 and NGC 3412. We measured their bar pattern speed using the Tremaine-Weinberg method, and derived the ratio, R, of the corotation radius to the length of the bar semi-major axis. For all the galaxies, R is consistent with being in the range from 1.0 and 1.4, i.e. that they host fast bars. This represents the largest sample of galaxies for which R has been measured this way. Taking into account the measured distribution of R and our measurement uncertainties, we argue that this is probably the true distribution of R. If this is the case, then the Tremaine-Weinberg method finds a distribution of R which is in agreement with that obtained by hydrodynamical simulations. We compared this result with recent high-resolution N-body simulations of bars in cosmologically-motivated dark matter halos,and conclude that these bars are not located inside centrally concentrated dark matter halos.

We show how the continuity equation can be used to determine pattern speeds in the Milky Way Galaxy (MWG). This method, first discussed by Tremaine & Weinberg in the context of external galaxies, requires projected positions, $(l,b)$, and line-of-sight velocities for a spatially complete sample of relaxed tracers. If the local standard of rest (LSR) has a zero velocity in the radial direction ($u_{\rm LSR}$), then the quantity that is measured is $\Delta V \equiv \Omega_p R_0 - V_{\rm LSR}$, where $\Omega_p$ is the pattern speed of the non-axisymmetric feature, $R_0$ is the distance of the Sun from the Galactic centre and $V_{\rm LSR}$ is the tangential motion of the LSR, including the circular velocity. We use simple models to assess the reliability of the method for measuring a single, constant pattern speed of either a bar or spiral in the inner MWG. We then apply the method to the OH/IR stars in the ATCA/VLA OH 1612 MHz survey of Sevenster et al, finding $\Delta V = 252 \pm 41$ km/s, if $u_{\rm LSR} = 0$. Assuming further that $R_0 = 8$ kpc and $V_{\rm LSR} = 220$ \kms, this gives $\Omega_p = 59\pm 5$ km/s/kpc with a possible systematic error of perhaps 10 km/s/kpc. The non-axisymmetric feature for which we measure this pattern speed must be in the disc of the MWG.

We use the Tremaine-Weinberg method, adapted to the Milky Way, to measure a pattern speed in OH/IR stars located in the inner MW.

I present some new results, from simulations and observations, of the constraints bar pattern speeds place on the dark matter content of disk galaxies.

We measured the bar pattern speed, $\om$, of the SB0 galaxy NGC 1023 using the Tremaine-Weinberg (1984) method with stellar-absorption slit spectroscopy. The morphology and kinematics of the \hi gas outside NGC 1023 suggest it suffered a tidal interaction, sometime in the past, with one of its dwarf companions. At present, however, the optical disc is relaxed. If the disc had been stabilized by a massive dark matter halo and formed its bar in the interaction, then the bar would have to be slow. We found $\om = 5.0 \pm 1.8$ \kmsa, so that the bar ends near its co-rotation radius. It is therefore rotating rapidly and must have a maximum disc.

We have used 2D Fabry-Perot absorption-line spectroscopy of the SB0 galaxy NGC 7079 to measure its bar pattern speed, $\om$. As in all previous cases of bar pattern speed measurements, we find a fast bar. We estimate that NGC 7079 has been undisturbed for at least the past Gyr or roughly 8 bar rotations, long enough for the bar to have slowed down significantly through dynamical friction if the disk is sub-maximal.

We show that bars in galaxy models having halos of moderate density and a variety of velocity distributions all experience a strong drag from dynamical friction unless the halo has large angular momentum in the same sense as the disk. The frictional drag decreases the bar pattern speed, driving the corotation point out to distances well in excess of those estimated in barred galaxies. The halo angular momentum required to avoid strong braking is unrealistically large, even when rotation is confined to the inner halo only. We conclude, therefore, that bars are able to maintain their observed high pattern speeds only if the halo has a central density low enough for the disk to provide most of the central attraction in the inner galaxy. We present evidence that this conclusion holds for all bright galaxies.

We report on absorption-line Fabry-Perot spectroscopy of the SB0 galaxy NGC 7079. Our data allow us to determine velocity and dispersion maps, which compare well with the results of slit observations. Our goals are to obtain a bar pattern speed and to determine the velocity ellipsoid as a function of radius; this is work in progress.

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

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

We describe fully self-consistent N-body experiments of barred galaxies with massive halos. A rotating bar is braked through dynamical friction with the halo, which occurs on a short time scale when the central density of the halo is high. On the other hand, friction is weak in a model with a central halo density low enough that the disk dominates the rotation curve in the inner parts; this model supports a fast bar for a Hubble time. We conclude that real barred galaxies in which the bar rotates rapidly, which is believed to be the rule, must have close to maximum disks.

A strong bar rotating within a massive halo should lose angular momentum to the halo through dynamical friction, as predicted by Weinberg. We have conducted fully self-consistent, numerical simulations of barred galaxy models with a live halo population and find that bars are indeed braked very rapidly. Specifically, we find that the bar slows sufficiently within a few rotation periods that the distance from the centre to co-rotation is more than twice the semi-major axis of the bar. Observational evidence (meagre) for bar pattern speeds seems to suggest that this ratio typically lies between 1.2 to 1.5 in real galaxies. We consider, a number of possible explanations for this discrepancy between theoretical prediction and observation, and conclude that no conventional alternative seems able to account for it.