Our location in the Milky Way provides an exceptional opportunity to gain insight on the galactic evolution processes, and complement the information inferred from observations of external galaxies. Since the Milky Way is a barred galaxy, the study of motions of individual stars in the bulge and disc is useful to understand the role of the bar. The Gaia mission enables such study by providing the most precise parallaxes and proper motions to date. In this theoretical work, we explore the effects of the bar on the distribution of higher-order moments --the skewness and kurtosis-- of the proper motions by confronting two simulated galaxies, one with a bar and one nearly axisymmetric, with observations from the latest Gaia data release (GaiaDR2). We introduce the code ASGAIA to account for observational errors of Gaia in the kinematical structures predicted by the numerical models. As a result, we find clear imprints of the bar in the skewness distribution of the longitudinal proper motion $\mu_\ell$ in GaiaDR2, as well as other features predicted for the next Gaia data releases.

Stellar feedback plays a significant role in modulating star formation, redistributing metals, and shaping the baryonic and dark structure of galaxies -- however, the efficiency of its energy deposition to the interstellar medium is challenging to constrain observationally. Here we leverage HST and ALMA imaging of a molecular gas and dust shell ($M_{H2} \sim 2\times 10^{5} ~{\rm M}_{\odot}$) in an outflow from the nuclear star forming ring of the galaxy NGC 3351, to serve as a boundary condition for a dynamical and energetic analysis of the outflowing ionised gas seen in our MUSE TIMER survey. We use \texttt{STARBURST99} models and prescriptions for feedback from simulations to demonstrate that the observed star formation energetics can reproduce the ionised and molecular gas dynamics -- provided a dominant component of the momentum injection comes from direct photon pressure from young stars, on top of supernovae, photoionisation heating and stellar winds. The mechanical energy budget from these sources is comparable to low luminosity AGN, suggesting that stellar feedback can be a relevant driver of bulk gas motions in galaxy centres - although here $\lesssim 10^{-3}$ of the ionized gas mass is escaping the galaxy. We test several scenarios for the survival/formation of the cold gas in the outflow, including in-situ condensation and cooling. Interestingly, the geometry of the molecular gas shell, observed magnetic field strengths and emission line diagnostics are consistent with a scenario where magnetic field lines aided survival of the dusty ISM as it was initially launched (with mass loading factor $\lesssim 1$) from the ring by stellar feedback. This system's unique feedback driven morphology can hopefully serve as a useful litmus test for feedback prescriptions in magnetohydrodynamical galaxy simulations.

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

Bars are common galactic structures in the local universe that play an important role in the secular evolution of galaxies, including the Milky Way. In particular, the velocity distribution of individual stars in our galaxy is useful to shed light on stellar dynamics, and provides information complementary to that inferred from the integrated light of external galaxies. However, since a wide variety of models reproduce the distribution of velocity and the velocity dispersion observed in the Milky Way, we look for signatures of the bar on higher-order moments of the line-of-sight velocity ($V_{los}$) distribution. We make use of two different numerical simulations --one that has developed a bar and one that remains nearly axisymmetric-- to compare them with observations in the latest APOGEE data release (SDSS DR14). This comparison reveals three interesting structures that support the notion that the Milky Way is a barred galaxy. A high skewness region found at positive longitudes constrains the orientation angle of the bar, and is incompatible with the orientation of the bar at $\ell=0^\circ$ proposed in previous studies. We also analyse the $V_{los}$ distributions in three regions, and introduce the Hellinger distance to quantify the differences among them. Our results show a strong non-Gaussian distribution both in the data and in the barred model, confirming the qualitative conclusions drawn from the velocity maps. In contrast to earlier work, we conclude it is possible to infer the presence of the bar from the kurtosis distribution.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) provides spectroscopic information of regions of the inner Milky Way inaccessible to optical surveys. We present the first large study of the metallicity distribution of the innermost Galactic regions based on homogeneous measurements from the SDSS Data Release 12 for 7545 red giant stars within 4.5 kpc of the Galactic center, with the goal to shed light on the structure and origin of the Galactic bulge. Stellar metallicities are found, through multiple-Gaussian decompositions, to be distributed in several components indicative of the presence of various stellar populations such as the bar, or the thin and the thick disk. A super-solar ([Fe/H]=+0.32) and a solar ([Fe/H]=+0.00) metallicity components, tentatively associated with the thin disk and the Galactic bar, respectively, seem to be the major contributors near the midplane. The solar-metallicity component extends outwards in the midplane but is not observed in the innermost regions. The central regions (within 3 kpc of the Galactic center) reveal, on the other hand, the presence of a significant metal-poor population ([Fe/H] =-0.46), tentatively associated with the thick disk, and which becomes the dominant component far from the midplane ($|Z| > +0.75 kpc). Varying contributions from these different components produce a transition region at +0.5 kpc < |Z| < +1.0 kpc characterized by a significant vertical metallicity gradient.

Boxy/peanut bulges are considered to be part of the same stellar structure as bars and both could be linked through the buckling instability. The Milky Way is our closest example. The goal of this letter is determining if the mass assembly of the different components leaves an imprint in their stellar populations allowing to estimate the time of bar formation and its evolution. To this aim we use integral field spectroscopy to derive the stellar age distributions, SADs, along the bar and disc of NGC 6032. The analysis shows clearly different SADs for the different bar areas. There is an underlying old (>=12 Gyr) stellar population for the whole galaxy. The bulge shows star formation happening at all times. The inner bar structure shows stars of ages older than 6 Gyrs with a deficit of younger populations. The outer bar region presents a SAD similar to that of the disc. To interpret our results, we use a generic numerical simulation of a barred galaxy. Thus, we constrain, for the first time, the epoch of bar formation, the buckling instability period and the posterior growth from disc material. We establish that the bar of NGC 6032 is old, formed around 10 Gyr ago while the buckling phase possibly happened around 8 Gyr ago. All these results point towards bars being long-lasting even in the presence of gas.

For several decades it has been known that stellar bars in disc galaxies can be triggered by interactions, or by internal processes such as dynamical instabilities. In this work, we explore the differences between these two mechanisms using numerical simulations. We perform two groups of simulations based on isolated galaxies, one group in which a bar develops naturally, and another group in which the bar could not develop in isolation. The rest of the simulations recreate 1:1 coplanar fly-by interactions computed with the impulse approximation. The orbits we use for the interactions represent the fly-bys in groups or clusters of different masses accordingly to the velocity of the encounter. In the analysis we focus on bars' amplitude, size, pattern speed and their rotation parameter, ${\cal R}=R_{CR}/R_{bar}$. The latter is used to define fast (${\cal R}<1.4$) and slow rotation (${\cal R}>1.4$). Compared with equivalent isolated galaxies we find that bars affected or triggered by interactions: (i) remain in the slow regime for longer; (ii) are more boxy in face-on views; (iii) they host kinematically hotter discs. Within this set of simulations we do not see strong differences between retrograde or prograde fly-bys. We also show that slow interactions can trigger bar formation.

Much of the inner Milky Way's (MW) global rotation and velocity dispersion patterns can be reproduced by models of secularly-evolved, bar-dominated bulges. More sophisticated constraints, including the higher moments of the line-of-sight velocity distributions (LOSVDs) and limits on the chemodynamical substructure, are critical for interpreting observations of the unresolved inner regions of extragalactic systems and for placing the MW in context with other galaxies. Here, we use SDSS-APOGEE data to develop these constraints, by presenting the first maps of the LOSVD skewness and kurtosis of metal-rich and metal-poor inner MW stars (divided at [Fe/H] = -0.4), and comparing the observed patterns to those that are seen both in N-body models and in extragalactic bars. Despite closely matching the mean velocity and dispersion, the models do not reproduce the observed LOSVD skewness patterns in different ways, which demonstrates that our understanding of the detailed orbital structure of the inner MW remains an important regime for improvement. We find evidence in the MW of the skewness-velocity correlation that is used as a diagnostic of extragalactic bar/bulges. This correlation appears in metal-rich stars only, providing further evidence for different evolutionary histories of chemically differentiated populations. We connect these skewness measurements to previous work on high-velocity "peaks" in the inner Galaxy, confirming the presence of that phenomenon, and we quantify the cylindrical rotation of the inner Galaxy, finding that the latitude-independent rotation vanishes outside of lon ~ 7 deg. Finally, we evaluate the MW data in light of select extragalactic bar diagnostics and discuss progress and challenges of using the MW as a resolved analog of unresolved stellar populations.

Massive high-redshift quiescent compact galaxies (nicknamed red nuggets) have been traditionally connected to present-day elliptical galaxies, often overlooking the relationships that they may have with other galaxy types. We use large bulge-disk decomposition catalogues based on the Sloan Digital Sky Survey (SDSS) to check the hypothesis that red nuggets have survived as compact cores embedded inside the haloes or disks of present-day massive galaxies. In this study, we designate a "compact core" as the bulge component that satisfies a prescribed compactness criterion. Photometric and dynamic mass-size and mass-density relations are used to show that, in the inner regions of galaxies at z ~ 0.1, there are "abundant" compact cores matching the peculiar properties of the red nuggets, an abundance comparable to that of red nuggets at z ~ 1.5. Furthermore, the morphology distribution of the present-day galaxies hosting compact cores is used to demonstrate that, in addition to the standard channel connecting red nuggets with elliptical galaxies, a comparable fraction of red nuggets might have ended up embedded in disks. This result generalises the inside-out formation scenario; present-day massive galaxies can begin as dense spheroidal cores (red nuggets), around which either a spheroidal halo or a disk are formed later.

Classical bulges in spiral galaxies are known to rotate but the origin of this observed rotational motion is not well understood. It has been shown recently that a low-mass classical bulge (ClB) in a barred galaxy can acquire rotation from absorbing a significant fraction of the angular momentum emitted by the bar. Our aim here is to investigate whether bars can spin up also more massive ClBs during the secular evolution of the bar, and to study the kinematics and dynamics of these ClBs. We use a set of self-consistent N-body simulations to study the interaction of ClBs with a bar that forms self-consistently in the disk. We use orbital spectral analysis to investigate the angular momentum gain by the classical bulge stars. We show that the ClBs gain, on average, about 2 - 6% of the disk's initial angular momentum within the bar region. Most of this angular momentum gain occurs via low-order resonances, particularly 5:2 resonant orbits. A density wake forms in the ClB which corotates and aligns with the bar at the end of the evolution. The spin-up process creates a characteristic linear rotation profile and mild tangential anisotropy in the ClB. The induced rotation is small in the centre but significant beyond $\sim2$ bulge half mass radii, where it leads to mass-weighted $V/\sigma \sim 0.2$, and reaches a local $V_{max}/\sigma_{in} \sim 0.5$ at around the scale of the bar. The resulting $V/\sigma$ is tightly correlated with the ratio of the bulge size to the bar size. In all models, a box/peanut bulge forms suggesting that composite bulges may be common. Bar-bulge resonant interaction in barred galaxies can provide some spin up of massive ClBs, but the process appears to be less efficient than for low-mass ClBs. Further angular momentum transfer due to nuclear bars or gas inflow would be required to explain the observed rotation if it is not primordial.

We present a kinematical study of 29 spiral galaxies included in the Spitzer Survey of Stellar Structure in Galaxies, using Halpha Fabry-Perot data obtained with the Galaxy Halpha Fabry-Perot System instrument at the William Herschel Telescope in La Palma, complemented with images in the R-band and in Halpha. The primary goal is to study the evolution and properties of the main structural components of galaxies through the kinematical analysis of the FP data, complemented with studies of morphology, star formation and mass distribution. In this paper we describe how the FP data have been obtained, processed and analysed. We present the resulting moment maps, rotation curves, velocity model maps and residual maps. Images are available in FITS format through the NASA/IPAC Extragalactic Database and the Centre de Donn\'ees Stellaires. With these data products we study the non-circular motions, in particular those found along the bars and spiral arms. The data indicate that the amplitude of the non-circular motions created by the bar does not correlate with the bar strength indicators. The amplitude of those non-circular motions in the spiral arms does not correlate with either arm class or star formation rate along the spiral arms. This implies that the presence and the magnitude of the streaming motions in the arms is a local phenomenon.

To provide a quantitative cosmological context to ongoing observational work on the formation histories and location of compact massive galaxies, we locate and study a sample of exceptionally compact systems in the Bolshoi simulation, using the dark matter structural parameters from a real, compact massive galaxy (NGC 1277) as a basis for our working criteria. We find that over 80% of objects in this nominal compact category are substructures of more massive groups or clusters, and that the probability of a given massive substructure being this compact increases significantly with the mass of the host structure; rising to ~30% for the most massive clusters in the simulation. Tracking the main progenitors of this subsample back to z=2, we find them all to be distinct structures with scale radii and densities representative of the population as a whole at this epoch. What does characterise their histories, in addition to mostly becoming substructures, is that they have almost all experienced below-average mass accretion since z=2; a third of them barely retaining, or even losing mass during the intervening 10 Gyr.

We construct dynamical models of the Milky Way's Box/Peanut (B/P) bulge, using the recently measured 3D density of Red Clump Giants (RCGs) as well as kinematic data from the BRAVA survey. We match these data using the NMAGIC Made-to-Measure method, starting with N-body models for barred discs in different dark matter haloes. We determine the total mass in the bulge volume of the RCGs measurement (+-2.2 x +- 1.4 x +- 1.2 kpc) with unprecedented accuracy and robustness to be 1.84 +- 0.07 x10^10 Msun. The stellar mass in this volume varies between 1.25-1.6 x10^10 Msun, depending on the amount of dark matter in the bulge. We evaluate the mass-to-light and mass-to-clump ratios in the bulge and compare them to theoretical predictions from population synthesis models. We find a mass-to-light ratio in the K-band in the range 0.8-1.1. The models are consistent with a Kroupa or Chabrier IMF, but a Salpeter IMF is ruled out for stellar ages of 10 Gyr. To match predictions from the Zoccali IMF derived from the bulge stellar luminosity function requires about 40% or 0.7 x10^10 Msun dark matter in the bulge region. The BRAVA data together with the RCGs 3D density imply a low pattern speed for the Galactic B/P bulge of 25-30 km.s-1.kpc-1. This would place the Galaxy among the slow rotators (R >= 1.5). Finally, we show that the Milky Way's B/P bulge has an off-centred X structure, and that the stellar mass involved in the peanut shape accounts for at least 20% of the stellar mass of the bulge, significantly larger than previously thought.

The Galactic bulge is a massive, old component of the Milky Way. It is known to host a bar, and it has recently been demonstrated to have a pronounced boxy/peanut structure in its outer region. Several independent studies suggest the presence of more than one stellar populations in the bulge, with different origins and a relative fraction changing across the bulge area. This is the first of a series of papers presenting the results of the Giraffe Inner Bulge Survey, carried out at the ESO-VLT with the multifibre spectrograph FLAMES. Spectra of ~5000 red clump giants in 24 bulge fields have been obtained at resolution R=6500, in the infrared Calcium triplet wavelength region at 8500 {\AA}. They are used to derive radial velocities and metallicities, based on new calibration specifically devised for this project. Radial velocities for another ~1200 bulge red clump giants, obtained from similar archive data, have been added to the sample. Higher resolution spectra have been obtained for 450 additional stars at latitude b=-3.5, with the aim of investigating chemical abundance patterns variations with longitude, across the inner bulge. In total we present here radial velocities for 6392 RC stars. We derive a radial velocity, and velocity dispersion map of the Milky Way bulge, useful to be compared with similar maps of external bulges, and to infer the expected velocities and dispersion at any line of sight. The K-type giants kinematics is consistent with the cylindrical rotation pattern of M-giants from the BRAVA survey. Our sample enables to extend this result to latitude b=-2, closer to the Galactic plane than probed by previous surveys. Finally, we find strong evidence for a velocity dispersion peak at (0,-1) and (0,-2), possibly indicative of a high density peak in the central 250 pc of the bulge

We study the dark and luminous mass distributions, circular velocity curves (CVC), line-of-sight kinematics, and angular momenta for a sample of 42 cosmological zoom simulations of massive galaxies. Using a temporal smoothing technique, we are able to reach large radii. We find that: (i)The dark matter halo density profiles outside a few kpc follow simple power-law models, with flat dark matter CVCs for lower-mass systems, and rising CVCs for high-mass haloes. The projected stellar density distributions at large radii can be fitted by Sersic functions with n>10, larger than for typical ETGs. (ii)The massive systems have nearly flat total CVCs at large radii, while the less massive systems have mildly decreasing CVCs. The slope of the CVC at large radii correlates with v_circ itself. (iii)The dark matter fractions within Re are in the range 15-30% and increase to 40-65% at 5Re. Larger and more massive galaxies have higher dark matter fractions. (iv)The short axes of simulated galaxies and their host dark matter haloes are well aligned and their short-to-long axis ratios are correlated. (v)The stellar vrms(R) profiles are slowly declining, in agreement with planetary nebulae observations in the outer haloes of most ETGs. (vi)The line-of-sight velocity fields v show that rotation properties at small and large radii are correlated. Most radial profiles for the cumulative specific angular momentum parameter lambda(R) are nearly flat or slightly rising, with values in [0.06,0.75] from 2Re to 5Re. (vii)Stellar mass, ellipticity at 5Re, and lambda(5Re) are correlated: the more massive systems have less angular momentum and are rounder, as for observed ETGs. (viii)More massive galaxies with a large fraction of accreted stars have radially anisotropic velocity distributions outside Re. Tangential anisotropy is seen only for galaxies with high fraction of in-situ stars. (Full abstract in PDF)

NGC 4494 is one of several intermediate-luminosity elliptical galaxies inferred to have an unusually diffuse dark matter halo. We use the chi^2-made-to-measure particle code NMAGIC to construct axisymmetric models of NGC 4494 from photometric and various kinematic data. The extended kinematics include light spectra in multiple slitlets out to 3.5 R_e, and hundreds of planetary nebulae velocities out to ~7 R_e, thus allowing us to probe the dark matter content and orbital structure in the halo. We use Monte Carlo simulations to estimate confidence boundaries for the halo parameters, given our data and modelling set-up. We find that the true potential of the dark matter halo is recovered within Delta G (merit function)<26 (Delta chi^2<59) at 70% confidence level (C.L.), and within Delta G<32 (Delta chi^2<70) at 90% C.L.. These numbers are much larger than the usually assumed Delta chi^2=2.3 (4.6) for 70% (90%) C.L. for two free parameters, perhaps case-dependent, but calling into question the general validity of the standard assumptions used for halo and black hole mass determinations. The best-fitting models for NGC 4494 have a dark matter fraction of about 0.6\pm0.1 at 5R_e (70% C.L.), and are embedded in a dark matter halo with circular velocity ~200 km/s. The total circular velocity curve (CVC) is approximately flat at v_c=220 km/s outside ~0.5R_e. The orbital anisotropy of the stars is moderately radial. These results are independent of the assumed inclination of the galaxy, and edge-on models are preferred. Comparing with the halos of NGC 3379 and NGC 4697, whose velocity dispersion profiles also decrease rapidly from the center outwards, the outer CVCs and dark matter halos are quite similar. NGC 4494 shows a particularly high dark matter fraction inside ~3R_e, and a strong concentration of baryons in the center.

Observations show a clear vertical metallicity gradient in the Galactic bulge, which is often taken as a signature of dissipative processes in the formation of a classical bulge. Various evidence shows, however, that the Milky Way is a barred galaxy with a boxy bulge representing the inner three-dimensional part of the bar. Here we show with a secular evolution N-body model that a boxy bulge formed through bar and buckling instabilities can show vertical metallicity gradients similar to the observed gradient, if the initial axisymmetric disk had a comparable radial metallicity gradient. In this framework the range of metallicities in bulge fields constrains the chemical structure of the Galactic disk at early times, before bar formation. Our secular evolution model was previously shown to reproduce inner Galaxy star counts and we show here that it also has cylindrical rotation. We use it to predict a full mean metallicity map across the Galactic bulge from a simple metallicity model for the initial disk. This map shows a general outward gradient on the sky as well as longitudinal perspective asymmetries. We also briefly comment on interpreting metallicity gradient observations in external boxy bulges.

Commissioning observations with the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III, have produced radial velocities (RVs) for ~4700 K/M-giant stars in the Milky Way bulge. These high-resolution (R \sim 22,500), high-S/N (>100 per resolution element), near-infrared (1.51-1.70 um; NIR) spectra provide accurate RVs (epsilon_v~0.2 km/s) for the sample of stars in 18 Galactic bulge fields spanning -1<l<20 deg, |b|<20 deg, and dec>-32 deg. This represents the largest NIR high-resolution spectroscopic sample of giant stars ever assembled in this region of the Galaxy. A cold (sigma_v~30 km/s), high-velocity peak (V_GSR \sim +200 km/s) is found to comprise a significant fraction (~10%) of stars in many of these fields. These high RVs have not been detected in previous MW surveys and are not expected for a simple, circularly rotating disk. Preliminary distance estimates rule out an origin from the background Sagittarius tidal stream or a new stream in the MW disk. Comparison to various Galactic models suggests that these high RVs are best explained by stars in orbits of the Galactic bar potential, although some observational features remain unexplained.

Secular evolution is one of the key routes through which galaxies evolve along the Hubble sequence. Not only the disk undergoes morphological and kinematic changes, but also a preexisting classical bulge may be dynamically changed by the secular processes driven primarily by the bar. We study the influence of a growing bar on the dynamical evolution of a low mass classical bulge such as might be present in galaxies like the Milky Way. Using self-consistent high resolution {\it N}-body simulations, we study how an initially isotropic non-rotating small classical bulge absorbs angular momentum emitted by the bar. The basic mechanism of this angular momentum exchange is through resonances and a considerable fraction of the angular momentum is channeled through Lagrange point (-1:1) and ILR (2:1) orbits. In the phase of rapid dynamical growth, also retrograde non-resonant orbits absorb significant angular momentum. As a result of this angular momentum gain, the initially non-rotating classical bulge transforms into a fast rotating, radially anisotropic and triaxial object, embedded in the similarly fast rotating boxy bulge formed from the disk. Towards the end of the evolution, the classical bulge develops cylindrical rotation. By that time, its inner regions host a "classical bulge-bar" whose distinct kinematics could serve as direct observational evidence for the secular evolution in the galaxy. Some implications of these results are discussed briefly.

Recent data from the VVV survey have strengthened evidence for a structural change in the Galactic bulge inwards of |l|<=4 deg. Here we show with an N-body barred galaxy simulation that a boxy bulge formed through the bar and buckling instabilities effortlessly matches measured bulge longitude profiles for red clump stars. The same simulation snapshot was earlier used to clarify the apparent boxy bulge - long bar dichotomy, for the same orientation and scaling. The change in the slope of the model longitude profiles in the inner few degrees is caused by a transition from highly elongated to more nearly axisymmetric isodensity contours in the inner boxy bulge. This transition is confined to a few degrees from the Galactic plane, thus the change of slope is predicted to disappear at higher Galactic latitudes. We also show that the nuclear star count map derived from this simulation snapshot displays a longitudinal asymmetry similar to that observed in the 2MASS data, but is less flattened to the Galactic plane than the 2MASS map. These results support the interpretation that the Galactic bulge originated from disk evolution, and question the evidence advanced from star count data for the existence of a secondary nuclear bar in the Milky Way.

The detailed dynamical structure of the bulge in the Milky Way is currently under debate. Although kinematics of the bulge stars can be well reproduced by a boxy-bulge, the possible existence of a small embedded classical bulge can not be ruled out. We study the dynamical evolution of a small classical bulge in a model of the Milky Way using a self-consistent high resolution N-body simulation. Detailed kinematics and dynamical properties of such a bulge are presented.

It has been known for some time that the Milky Way is a barred disk galaxy. More recently several studies inferred from starcount observations that the Galaxy must contain a separate, new, flat long bar component, twisted relative to the barred bulge. Here we use a simulation with a boxy bulge and bar to suggest that these observations can be reproduced with a single structure. In this simulation a stellar bar evolved from the disk, and the boxy bulge originated from it through secular evolution and the buckling instability. We calculate starcount distributions for this model at different longitudes and latitudes, in a similar way as observers have done for resolved starcounts. Good agreement between the simulation and the observations can be achieved for a suitable snapshot, even though the simulation has a single boxy bulge and bar} structure. In this model, part of the long bar signature is due to a volume effect in the starcounts, and another part is due to chosing a snapshot in which the planar part of the boxy bulge and bar has developed leading ends through interaction with the adjacent spiral arm heads. We also provide predictions from this model for the line-of-sight velocity distributions at the longitudes with the long bar signature, for comparison with upcoming surveys.

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A.

We present the most complete atlas of nuclear rings to date. We include 113 rings found in 107 galaxies, six of which are elliptical galaxies, five are highly inclined disc galaxies, 18 are unbarred disc galaxies, and 78 are barred disc galaxies. Star-forming nuclear rings occur in 20% of disc galaxies with types between T=-3 and T=7. We aim to explore possible relationships between the size and morphology of the rings and various galactic parameters. We produce colour index and structure maps, as well as Halpha and Paalpha continuum-subtracted images from HST archival data. We derive ellipticity profiles from H-band 2MASS images in order to detect bars and find their metric parameters. We measure the non-axisymmetric torque parameter, Qg, and search for correlations between bar, ring metric parameters, and Qg. Our atlas of nuclear rings includes star-forming and dust rings. Nuclear rings span a range from a few tens of parsecs to a few kiloparsecs in radius. Star-forming nuclear rings can be found in a wide range of morphological types, from S0 to Sd, with a peak in the distribution between Sab and Sb, and without strong preference for barred galaxies. Dust nuclear rings are found in elliptical and S0 galaxies. For barred galaxies, the maximum radius that a nuclear ring can reach is a quarter of the bar radius. We found a nearly random distribution of PA offsets between nuclear rings and bars. There is some evidence that nuclear ring ellipticity is limited by bar ellipticity. We confirm that the maximum relative size of a star-forming nuclear ring is inversely proportional to the non-axisymmetric torque parameter, Qg, and that the origin of nuclear rings, even the ones in non-barred hosts, are closely linked to the existence of dynamical resonances.

We test the theoretical prediction that the straightest dust lanes in bars are found in strongly barred galaxies, or more specifically, that the degree of curvature of the dust lanes is inversely proportional to the strength of the bar. The test used archival images of barred galaxies for which a reliable non-axisymmetric torque parameter (Qb) and the radius at which Qb has been measured (r(Qb)) have been published in the literature. Our results confirm the theoretical prediction but show a large spread that cannot be accounted for by measurement errors. We simulate 238 galaxies with different bar and bulge parameters in order to investigate the origin of the spread in the dust lane curvature versus Qb relation. From these simulations, we conclude that the spread is greatly reduced when describing the bar strength as a linear combination of the bar parameters Qb and the quotient of the major and minor axis of the bar, a/b. Thus we conclude that the dust lane curvature is predominantly determined by the parameters of the bar.

We have applied the Tremaine-Weinberg method to 10 late-type barred spiral galaxies using data cubes, in H-alpha emission, from the GHAFAS and FANTOMM Fabry-Perot spectrometers. We have combined the derived bar (and/or spiral) pattern speeds with angular frequency plots to measure the corotation radii for the bars in these galaxies. We base our results on a combination of this method with a morphological analysis designed to estimate the corotation radius to bar-length ratio using two independent techniques on archival near infrared images, and although we are aware of the limitation of the application of the Tremaine-Weinberg method using ionised gas observations, we find consistently excellent agreement between bar and spiral arm parameters derived using different methods. In general, the corotation radius, measured using the Tremaine-Weinberg method, is closely related to the bar length, measured independently from photometry and consistent with previous studies. Our corotation/bar-length ratios and pattern speed values are in good agreement with general results from numerical simulations of bars. In systems with identified secondary bars, we measure higher H-alpha velocity dispersion in the circumnuclear regions, whereas in all the other galaxies, we detect flat velocity dispersion profiles. The excellent agreement between the Tremaine-Weinberg method results and the morphological analysis and bar parameters in numerical simulations, suggests that although the H-alpha emitting gas does not obey the continuity equation, it can be used to derive the bar pattern speed.

We have conducted a deep multi-wavelength analysis (0.15-160 mum) to study the outer region of the nearby galaxy M94. We show that the non-optical data support the idea that the outskirts of this galaxy is not formed by a closed stellar ring (as traditionally claimed in the literature) but by a spiral arm structure. In this sense, M94 is a good example of a Type III (anti-truncated) disk galaxy having a very bright outer disk. The outer disk of this galaxy contains ~23% of the total stellar mass budget of the galaxy and contributes ~10% of the new stars created showing that this region of the galaxy is active. In fact, the specific star formation rate of the outer disk (~0.012 Gyr^{-1}) is a factor of ~2 larger (i.e. the star formation is more efficient per unit stellar mass) than in the inner disk. We have explored different scenarios to explain the enhanced star formation in the outer disk. We find that the inner disk (if considered as an oval distortion) can dynamically create a spiral arm structure in the outer disk which triggers the observed relatively high star formation rate as well as an inner ring similar to what is found in this galaxy.

We discuss the formation and evolution of boxy/peanut bulges (B/Ps) and present new simulations results. Orbital structure studies show that B/Ps are parts of bars seen edge-on, they have their origin in vertical instabilities of the disc material and they are somewhat shorter in extent than bars. When the bar forms it is vertically thin, but after a time of the order of a Gyr it experiences a vertical instability and buckles. At that time the strength of the bar decreases, its inner part becomes thicker, so that, seen edge-on, it acquires a peanut or boxy shape. A second buckling episode is seen in simulations with strong bars, accompanied by a further thickening of the B/P and a weakening of the bar. Quantitatively, this evolution depends considerably on the properties of the halo and particularly on the extent of its core. This influences the amount of angular momentum exchanged within the galaxy, emitted by near-resonant material in the bar region and absorbed by near-resonant material in the halo and in the outer disc. Haloes with small cores generally harbour stronger bars and B/Ps and they often witness double buckling.

Boxy/peanut bulges in disk galaxies have been associated to stellar bars. In this talk, we discuss the different properties of such bulges and their relation with the corresponding bar, using a very large sample of a few hundred numerical N-body simulations. We present and inter-compare various methods of measuring the boxy/peanut bulge properties, namely its strength, shape and possible asymmetry. Some of these methods can be applied to both simulations and observations. Our final goal is to get correlations that will allow us to obtain information on the boxy/peanut bulge for a galaxy viewed face-on as well as information on the bars of galaxies viewed edge-on.

After presenting three ways of defining a bulge component in disc galaxies, we introduce the various types of bulges, namely the classical bulges, the boxy/peanut bulges and the disc-like bulges. We then discuss three specific topics linked to bulge formation and evolution, namely the coupled time evolution of the bar, buckling and peanut strengths; the effect of velocity anisotropy on peanut formation; and bulge formation via bar destruction.

Many barred galaxies show a set of symmetric enhancements at the ends of the stellar bar, called {\it ansae}, or the ``handles'' of the bar. The ansa bars have been in the literature for some decades, but their origin has still not been specifically addressed, although, they could be related to the growth process of bars. But even though ansae have been known for a long time, no statistical analysis of their relative frequency of occurrence has been performed yet. Similarly, there has been no study of the varieties in morphology of ansae even though significant morphological variations are known to characterise the features. In this paper, we make a quantitative analysis of the occurrence of ansae in barred galaxies, making use of {\it The de Vaucouleurs Atlas of Galaxies} by Buta and coworkers. We find that $\sim 40%$ of SB0's show ansae in their bars, thus confirming that ansae are common features in barred lenticulars. The ansa frequency decreases dramatically with later types, and hardly any ansae are found in galaxies of type Sb or later. The bars in galaxies with ansae are stronger in the median than those in galaxies without ansae, but the presence of inner and outer rings is not related to the presence of ansae. Implications of these results and theories for the possible origin of ansae are discussed briefly.

We analyze evolution of live disk-halo systems in the presence of various gas fractions, f_gas less than 8% in the disk. We addressed the issue of angular momentum (J) transfer from the gas to the bar and its effect on the bar evolution. We find that the weakening of the bar, reported in the literature, is not related to the J-exchange with the gas, but is caused by the vertical buckling instability in the gas-poor disks and by a steep heating of a stellar velocity dispersion by the central mass concentration (CMC) in the gas-rich disks. The gas has a profound effect on the onset of the buckling -- larger f_gas brings it forth due to the more massive CMCs. The former process leads to the well-known formation of the peanut-shaped bulges, while the latter results in the formation of progressively more elliptical bulges, for larger f_gas. The subsequent (secular) evolution of the bar differs -- the gas-poor models exhibit a growing bar while gas-rich models show a declining bar whose vertical swelling is driven by a secular resonance heating. The border line between the gas-poor and -rich models lies at f_gas ~ 3% in our models, but is model-dependent and will be affected by additional processes, like star formation and feedback from stellar evolution. The overall effect of the gas on the evolution of the bar is not in a direct J transfer to the stars, but in the loss of J by the gas and its influx to the center that increases the CMC. The more massive CMC damps the vertical buckling instability and depopulates orbits responsible for the appearance of peanut-shaped bulges. The action of resonant and non-resonant processes in gas-poor and gas-rich disks leads to a converging evolution in the vertical extent of the bar and its stellar dispersion velocities, and to a diverging evolution in the bulge properties.

Evolution of stellar bars in disk galaxies is accompanied by dynamical instabilities and secular changes. Following the vertical buckling instability, the bars are known to weaken dramatically and develop a pronounced boxy/peanut shape when observed edge-on. Using high-resolution N-body simulations of stellar disks embedded in live axisymmetric dark matter halos, we have investigated the long-term changes in the bar morphology, specifically the evolution of the bar size, its vertical structure and exchange of angular momentum. We find that following the initial buckling, the bar resumes its growth from deep inside the corotation radius and follows the Ultra-Harmonic resonance thereafter. We also find that this secular bar growth triggers a spectacular secondary vertical buckling instability which leads to the appearance of characteristic boxy/peanut/X-shaped bulges. The secular bar growth is crucial for the recurrent buckling to develop. Furthermore, the secondary buckling is milder, persists over a substantial period of time, ~3 Gyr, and can have observational counterparts. Overall, the stellar bars show recurrent behavior in their properties and evolve by increasing their linear and vertical extents, both dynamically and secularly. We also demonstrate explicitly that the prolonged growth of the bar is mediated by continuous angular momentum transfer from the disk to the surrounding halo, and that this angular momentum redistribution is resonant in nature -- a large number of lower resonances trap disk and halo particles and this trapping is robust, in a broad agreement with the earlier results in the literature.

Young stellar bars in disk galaxies experience a vertical buckling instability which terminates their growth and thickens them, resulting in a characteristic peanut/boxy shape when viewed edge on. Using N-body simulations of galactic disks embedded in live halos, we have analyzed the bar structure throughout this instability and found that the outer third of the bar dissolves completely while the inner part (within the vertical inner Lindblad resonance) becomes less oval. The bar acquires the frequently observed peanut/boxy-shaped isophotes. We also find that the bar buckling is responsible for a mass injection above the plane, which is subsequently trapped by specific 3-D families of periodic orbits of particular shapes explaining the observed isophotes, in line with previous work. Using a 3-D orbit analysis and surfaces of sections, we infer that the outer part of the bar is dissolved by a rapidly widening stochastic region around its corotation radius -- a process related to the bar growth. This leads to a dramatic decrease in the bar size, decrease in the overall bar strength and a mild increase in its pattern speed, but is not expected to lead to a complete bar dissolution. The buckling instability appears primarily responsible for shortening the secular diffusion timescale to a dynamical one when building the boxy isophotes. The sufficiently long timescale of described evolution, ~1 Gyr, can affect the observed bar fraction in local universe and at higher redshifts, both through reduced bar strength and the absence of dust offset lanes in the bar.