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

Aims. The orbital structure of galaxies is strongly influenced by the accuracy of the force calculation during orbit integration. We explore the accuracy of force calculations for two expansion methods and determine which one is preferable for orbit integration. Methods. We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations.We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations. Results. For the generalized Dehnen density, the CylSP method is more accurate than the HO method for nearly all inner power-law indices and shapes at all radii. For $N$-body realizations of the Dehnen model, or snapshots of an $N$-body simulation, the CylSP method is more accurate than the HO method in the central region for the oblate, prolate, and triaxial Hernquist profiles if the particle number is more than $5\times 10^5$. For snapshots of the Hernquist models with spherical shape, the HO method is preferred. (Abridged)

Bars in disc-dominated galaxies are able to drive gas inflow inside the corotation radius, thus enhancing the central star formation rate (SFR). Previous work, however, has found that disc-dominated galaxies with centrally suppressed SFR frequently host a bar. Here we investigate possible causes for the suppression of central SFR in such cases. We compare physical properties of a sample of disc-dominated barred galaxies with high central SFR (HC galaxies) with those of a sample of disc-dominated barred galaxies with low central SFR (LC galaxies). We find that the two samples have on average similar HI content and bars of similar strength. But we also find that the HCs have bluer colors than LCs, and that outside the bar region they host stronger spiral arms than the LCs where closed rings are more often seen. We discuss and evaluate the possible causes for the suppression of the central SFR in the LC galaxies as opposed to its enhancement in the HC galaxies.

We show how wet major mergers can create disc galaxies in a cosmological context, using the Illustris simulation. We select a sample of 38 disc galaxies having experienced a major merger in their history with no subsequent significant minor merger, and study how they transform into discs after the merger. In agreement with what was previously found in controlled simulations of such mergers, we find that their disc is built gradually from young stars formed after the merger in the disc region, while the old stars born before the merger form an ellipsoidal component. Focusing on one fiducial case from our sample, we show how the gas was initially dispersed in the halo region right after the merger, but is then accreted onto a disc to form stars, and builds the disc component. We then select a sample of major mergers creating elliptical galaxies, to show that those cases correspond mainly to dry mergers, where the lack of star formation prevents the formation of a disc in the remnant galaxy. The amount of gas in the remnant galaxy after the merger is therefore essential to determine the final outcome of a major merger.

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.

Recently, many orbital studies in barred galaxy potentials have revealed the existence of orbits which are not trapped around x1 tree orbits, but could be potentially appropriate building blocks for bars. These findings question the uniqueness of the x1 family as the standard paradigm of orbital motion in galactic bars. The main goal of this paper is to investigate the role that such orbits could play in shaping the morphology of bars. We trace the morphological patterns appearing in the face-on and edge-on views of the non-periodic orbits presented in these studies and we show that they are introduced in the system by second type bifurcations of x1. For this purpose, we use a typical 3D Ferrers bar model and follow the radial and vertical bifurcations of the x1 family considered as being multi periodic, with multi = 2, 3, 5. The variation of the stability indices of x1 in the multi = 2, 3 cases give us also the four- and six- periodic orbits, respectively. We tabulate these orbits including all information necessary to assess their role as appropriate building blocks. We discuss their stability and their extent, as well as their size and morphological evolution, as a function of energy. We conclude that even the most important of the multi-periodic orbits presented in Tables 2 to 5 are less appropriate building blocks for bars than the families of the x1 tree at the same energy.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We investigate how the properties of spiral arms relate to other fundamental galaxy properties. To this end, we use previously published measurements of those properties, and our own measurements of arm-interarm luminosity contrasts for a large sample of galaxies, using 3.6$\mu$m images from the Spitzer Survey of Stellar Structure in Galaxies. Flocculent galaxies are clearly distinguished from other spiral arm classes, especially by their lower stellar mass and surface density. Multi-armed and grand-design galaxies are similar in most of their fundamental parameters, excluding some bar properties and the bulge-to-total luminosity ratio. Based on these results, we discuss dense, classical bulges as a necessary condition for standing spiral wave modes in grand-design galaxies. We further find a strong correlation between bulge-to-total ratio and bar contrast, and a weaker correlation between arm and bar contrasts.

We study the spatially resolved star formation history and molecular gas distribution of 58 nearby galaxies, using integral field spectroscopy from the CALIFA survey and CO $J=1\rightarrow 0$ intensity mapping from the CARMA EDGE survey. We use the 4000 \AA\ break (D$_n$4000), the equivalent width of the H$\delta$ absorption line (EW H$\delta_A$), and the equivalent width of the H$\alpha$ emission line (EW H$\alpha$) to measure the recent star formation history (SFH) of these galaxies. We measure radial profiles of the three SFH indicators and molecular gas mass surface density, from which we measure the level of centrally enhanced star formation and the molecular gas concentration. When we separate our galaxies into categories of barred (17 galaxies), unbarred (24 galaxies), and merging/paired (17 galaxies) we find that the galaxies which have centrally-enhanced star formation (19/58) are either barred (13/19) or in mergers/pairs (6/19) with relatively high molecular gas concentrations. A comparison between our barred galaxies and a snapshot of a hydrodynamic $N$-body simulation of a barred galaxy shows that the current theory of bar formation and evolution can qualitatively reproduce the main features of the observed galaxies in our sample, including both the sharp decrease of stellar age in the galactic center and the gradual decrease of age with increasing distance from center. These findings provide substantial evidence for a picture in which cold gas is transported inward by a bar or tidal interaction, which leads to the growth and rejuvenation of star formation in the central region.

Bars in galaxies may form not only through instability, but also as a result of an interaction with another galaxy. In particular, they may appear in disky dwarf galaxies interacting with Milky Way-like galaxies. Here we report the results of N-body/SPH simulations of such dwarfs orbiting in the static potential of a larger galaxy. We used several models of the dwarf galaxy, all of the same mass, but covering a large range of gas fractions: 0, 30 and 70%. We also tested the impact of subgrid star formation processes. In all cases bars of similar length formed in the stellar disk of the dwarfs at the first pericenter passage. However, unexpectedly, the gaseous component remained approximately axisymmetric and unaffected by the bar potential. The bar properties did not change significantly between two consecutive pericenters. The impact of the later encounters with the host depends strongly on the exact orientation of the bar at the pericenter. When the bar is spun up by the tidal force torque, it is also shortened. Conversely, if it is slowed down, it gets longer. In the models with a low gas fraction the bars were more pronounced and survived until the end of the simulations, while in the dwarfs with a high gas fraction the bars were destroyed after the second or third pericenter passage. In terms of the ratio of the corotation radius to the bar length, the bars are slow, and remain so independently of the encounters with the host.

The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, making this the largest sample studied so far in this way. Our sample is selected from the MaNGA first public data release as part of SDSS Data Release 13 according mainly to the axis ratio and position angle difference between the bar and disc, while kinematic data is from the later SDSS Data Release 14. We have used both the photometric position angle from the photometric image and the kinematic position angle from the stellar velocity map to derive the pattern speed. Combining three independent bar length measurements and the circular velocity from Jeans Anisotropic modelling (JAM), we also determine the dimensionless ratio $\cal{R}$ of the corotation radius to the bar length. We find that the galaxy's position angle is the main uncertainty in determining the bar pattern speed. The kinematic position angle leads to fewer ultrafast bars than the photometric position angle, and this could be due to the method of measuring the kinematic position angle. We study the dependence of $\cal{R}$ values on galaxy properties such as the dark matter fraction from JAM modelling and the stellar age and metallicity from stellar population synthesis (SPS). A positive correlation between the bar length and bar strength is found: the longer the bar, the stronger the bar. However, no other significant correlations are found. This may result from errors in deriving the $\cal{R}$ values or from the complex formation and slowdown processes of galactic bars.

Asymmetries in HI in galaxies are often caused by the interaction with close companions, making isolated galaxies an ideal framework to study secular evolution. The AMIGA project has demonstrated that isolated galaxies show the lowest level of asymmetry in their HI integrated profiles, yet some present significant asymmetries. CIG 96 (NGC 864) is a representative case reaching a 16% level. Our aim is to investigate the HI asymmetries of this spiral galaxy and what processes have triggered the star-forming regions observed in the XUV pseudoring. We performed deep optical observations at CAHA 1.23m, 2.2m and VST telescopes. We reach surface brightness (SB) limits of mu_2.2m = 27.5 mag arcsec-2 (Cous R) and mu_VST = 28.7mag arcsec-2 (r) that show the XUV pseudoring of the galaxy in detail. Additionally, a wavelet filtering of the HI data cube from our deep observations with E/VLA telescope allowed us to reach a column density of N_HI = 8.9x10^18 cm -2 (5sigma) (28"x28" beam), lower than in any isolated galaxy. We confirm that the HI extends farther than 4xr_25 in all directions. Furthermore, we detect for the first time two gaseous structures (10^6 Msol) in the outskirts. The g-r colour index image from 1.23m shows extremely blue colours in certain regions of the pseudoring where N_HI>8.5x10^20 cm-2 , whereas the rest show red colours. Galactic cirrus contaminate the field, setting an unavoidable detection limit at 28.5mag arcsec-2 (r). We detect no stellar link within 1degx1deg or gaseous link within 40'x40' between CIG 96 and any companion. The isolation criteria rule out interactions with other similar-sized galaxies for at least 2.7Gyr. Using existing stellar evolution models, the age of the pseudoring is estimated at 1Gyr or older. Undetected previously accreted companions and cold gas accretion remain as the main hypothesis to explain the optical pseudoring and HI features of CIG 96.

We present a hydrodynamical model of M33 and its recent interaction with M31. This scenario was previously proposed in the literature in order to explain the distorted gaseous and stellar disks of M33, as well as the increased star formation rate in both objects around 2 Gyr ago. We used an orbit integration scheme to find which estimate of the transverse velocity of M31 more favors the interaction scenario and then tried to reproduce it in our simulations. M33 was modeled as a stellar and gaseous disk embedded in a live dark-matter halo, while M31 was approximated only with a live dark halo. In the simulations the two galaxies passed each other with a pericenter distance of 37 kpc. Tides excited a two-armed spiral structure in the M33 disk, which is found to be the predominant spiral signal in the observed galaxy and has been long known as a feature easily induced by tidal interactions. We found that the gaseous warp produced by the interaction did not resemble enough the observed one and we performed an additional simulation including the hot gas halo of M31 to show that this feature can be properly reproduced by tidal forces and ram pressure stripping acting simultaneously on the gaseous disk. In addition to the spiral arms, tidal forces produced the stellar stream similar to the observed one and triggered a star formation burst at similar radii as it is observed.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

I present examples of how chemo-dynamical N-body simulations can help understanding the structure and evolution of the inner Galaxy. Such simulations reproduce the observed links between kinematics, morphology and chemistry in the bar/bulge region and explain them by the self-consistent cohabitation of a number of components. Galactic archaeology, applied to simulation snapshots, explains the sequence in which the stars of the various components were formed. The thick disc stars form earlier than those of the thin disc and in a much shorter time scale. The bar in the thick disc is horizontally thicker than that of the thin disc and has a different vertical morphology. The Galaxy's inner disc scalelength is much smaller than what is expected from nearby galaxies of similar stellar mass.

MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.

We use hydrodynamical simulations of a Cartwheel-like ring galaxy, modelled as a nearly head-on collision of a small companion with a larger disc galaxy, to probe the evolution of the gaseous structures and flows, and to explore the physical conditions setting the star formation activity. Star formation is first quenched by tides as the companion approaches, before being enhanced shortly after the collision. The ring ploughs the disc material as it radially extends, and almost simultaneously depletes its stellar and gaseous reservoir into the central region, through the spokes, and finally dissolve 200 Myr after the collision. Most of star formation first occurs in the ring before this activity is transferred to the spokes and then the nucleus. We thus propose that the location of star formation traces the dynamical stage of ring galaxies, and could help constrain their star formation histories. The ring hosts tidal compression associated with strong turbulence. This compression yields an azimuthal asymmetry, with maxima reached in the side furthest away from the nucleus, which matches the star formation activity distribution in our models and in observed ring systems. The interaction triggers the formation of star clusters significantly more massive than before the collision, but less numerous than in more classical galaxy interactions. The peculiar geometry of Cartwheel-like objects thus yields a star (cluster) formation activity comparable to other interacting objects, but with notable second order differences in the nature of turbulence, the enhancement of the star formation rate, and the number of massive clusters formed.

We investigate how the properties of spiral arms relate to other fundamental galaxy properties, including bars and disc breaks. We use previously published measurements of those properties, and our own measurements of arm and bar contrasts for a large sample of galaxies, using $3.6 \mathrm{\mu m}$ images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). Flocculent galaxies are clearly distinguished from other spiral arm classes, especially by their lower stellar mass and surface density. Multi-armed and grand-design galaxies are similar in most of their fundamental parameters, excluding some bar properties and the bulge-to-total ratio. Based on these results, we revisit the sequence of spiral arm classes, and discuss classical bulges as a necessary condition for standing spiral wave modes in grand-design galaxies. We find a strong correlation between bulge-to-total ratio and bar contrast, and a weaker correlation between arm and bar contrasts. Barred and unbarred galaxies exhibit similar arm contrasts, but the highest arm contrasts are found exclusively in barred galaxies. Interestingly, the bar contrast, and its increase from flocculent to grand-design galaxies, is systematically more significant than that of the arm contrast. We corroborate previous findings concerning a connection between bars and disc breaks. In particular, in grand-design galaxies the bar contrast correlates with the normalised disc break radius. This does not hold for other spiral arm classes or the arm contrast. Our measurements of arm and bar contrast and radial contrast profiles are publicly available.

Bars in galaxies may develop through a global instability or due to an interaction with another system. We study bar formation in disky dwarf galaxies orbiting a Milky Way-like galaxy. We employ $N$-body simulations to study the impact of initial orbital parameters: the size of the dwarf galaxy orbit and the inclination of its disc with respect to the orbital plane. In all cases a bar develops in the center of the dwarf during the first pericenter on its orbit around the host. Between subsequent pericenter passages the bars are stable, but at the pericenters they are usually weakened and shortened. The initial properties and details of the further evolution of the bars depend heavily on the orbital configuration. We find that for the exactly prograde orientation, the strongest bar is formed for the intermediate-size orbit. On the tighter orbit, the disc is too disturbed and stripped to form a strong bar. On the wider orbit, the tidal interaction is too weak. The dependence on the disc inclination is such that weaker bars form in more inclined discs. The bars experience either a very weak buckling or none at all. We do not observe any secular evolution, possibly because the dwarfs are perturbed at each pericenter passage. The rotation speed of the bars can be classified as slow ($R_\mathrm{CR}/l_\mathrm{bar}\sim2-3$). We attribute this to the loss of a significant fraction of the disc's rotation during the encounter with the host galaxy.

We study the effect of angular momentum on the surface density profiles of disc galaxies, using high resolution simulations of major mergers whose remnants have downbending radial density profiles (type II). As described in the previous papers of this series, in this scenario, most of the disc mass is acquired after the collision via accretion from a hot gaseous halo. We find that the inner and outer disc scalelengths, as well as the break radius, correlate with the total angular momentum of the initial merging system, and are larger for high angular momentum systems. We follow the angular momentum redistribution in our simulated galaxies, and find that, like the mass, the disc angular momentum is acquired via accretion, i.e. to the detriment of the gaseous halo. Furthermore, high angular momentum systems give more angular momentum to their discs, which affects directly their radial density profile. Adding simulations of isolated galaxies to our sample, we find that the correlations are valid also for disc galaxies evolved in isolation. We show that the outer part of the disc at the end of the simulation is populated mainly by inside-out stellar migration, and that in galaxies with higher angular momentum, stars travel radially further out. This, however, does not mean that outer disc stars (in type II discs) were mostly born in the inner disc. Indeed, generally the break radius increases over time, and not taking this into account leads to overestimating the number of stars born in the inner disc.

Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger simulations can produce a non-realistic central mass concentration, and we propose simple, parametric, AGN-like feedback as a solution to this problem. Our AGN-like feedback algorithm is very simple: at each time-step we take all particles whose local volume density is above a given threshold value and increase their temperature to a preset value. We also compare the GADGET3 and GIZMO codes, by applying both of them to the same initial conditions. Results: We show that the evolution of isolated protogalaxies resembles the evolution of disk galaxies, thus arguing that our protogalaxies are well suited for our merger simulations. We demonstrate that the problem with the unphysical central mass concentration in our merger simulations is further aggravated when we increase the resolution. We show that our AGN-like feedback removes this non-physical central mass concentration, and thus allows the formation of realistic bars. Note that our AGN-like feedback mainly affects the central region of a model, without significantly modifying the rest of the galaxy. We demonstrate that, in the context of our kind of simulation, GADGET3 gives results which are very similar to those obtained with the PSPH (density independent SPH) flavor of GIZMO. Moreover, in the examples we tried, the differences between the results of the two flavors of GIZMO, namely PSPH, and MFM (mesh-less algorithm) are similar to and, in some comparisons, larger than the differences between the results of GADGET3 and PSPH.

I review briefly some dynamical models of structures in the outer parts of disc galaxies, including models of polar rings, tidal tails and bridges. I then discuss the density distribution in the outer parts of discs. For this, I compare observations to results of a model in which the disc galaxy is in fact the remnant of a major merger, and find good agreement. This comparison includes radial profiles of the projected surface density and of stellar age, as well as time evolution of the break radius and of the inner and outer disc scale lengths. I also compare the radial projected surface density profiles of dynamically motivated mono-age populations and find that, compared to older populations, younger ones have flatter density profiles in the inner region and steeper in the outer one. The break radius, however, does not vary with stellar age, again in good agreement with observations.

We use N-body chemo-dynamic simulations to study the coupling between morphology, kinematics and metallicity of the bar/bulge region of our Galaxy. We make qualitative comparisons of our results with available observations and find very good agreement. We conclude that this region is complex, since it comprises several stellar components with different properties -- i.e. a boxy/peanut bulge, thin and thick disc components, and, to lesser extents, a disky pseudobulge, a stellar halo and a small classical bulge -- all cohabiting in dynamical equilibrium. Our models show strong links between kinematics and metallicity, or morphology and metallicity, as already suggested by a number of recent observations. We discuss and explain these links.

We present a pilot study on the nearby massive barred galaxy NGC 1291, in which we use dynamical modelling to constrain the disc mass-to-light ratio (M/L), thus breaking the degeneracy between the baryonic and dark matter in its central regions. We use the gas, specifically the morphology of the dust lanes on the leading side of the bar, as a tracer of the underlying gravitational potential. We run a large number of hydrodynamic gas response simulations, in potentials obtained directly from near-infrared images of the galaxy, which have three free parameters: the M/L, the bar pattern speed and the height function. We explore the three-dimensional parameter space, by comparing the morphology of the shocks created in the gas response simulations with those of the observed dust lanes, and find the best-fitting models; these suggest that the M/L of NGC 1291 agrees with that predicted by stellar population synthesis models in the near-infrared ($\approx$0.6\,$M_{\odot}/L_{\odot}$), which leads to a borderline maximum disc for this galaxy. Furthermore, we find that the bar rotates fast, with a corotation radius which is $\leq$ 1.4 times the bar length.

We have studied a representative sample of intermediate-mass galaxies at z~1, observed by the kinematic survey KMOS3D. We have re-estimated the kinematical parameters from the published kinematic maps and analysed photometric data from HST to measure optical disk inclinations and PAs. We find that only half of the z~1 galaxies show kinematic properties consistent with rotating disks, using the same classification scheme than that adopted by the KMOS3D team. Because merger orbital motions can also brought rotation, we have also analysed galaxy morphologies from the available HST imagery. Combining these results to those from kinematics, it leads to a full morpho-kinematic classification. To test the robustness of the latter for disentangling isolated disks from mergers, we confronted the results with an analysis of pairs from the open-grism redshift survey 3D-HST. All galaxies found in pairs are affected by either kinematic and/or morphological perturbations. Conversely, all galaxies classified as virialized spirals are found to be isolated. A significant fraction (one fourth) of rotating disks classified from kinematics by the KMOS3D team are found in pairs, which further supports the need for a morpho-kinematic classification. It results that only one third of z~1 galaxies are isolated and virialized spirals, while 58% of them are likely involved in a merger sequence, from first approach to disk rebuilding. The later fraction is in good agreement with the results of semi-empirical {\Lambda}CDM models, supporting a merger-dominated hierarchical scenario as being the main driver of galaxy formation at least during the last 8 billion years.

We study the surface brightness profiles of disc galaxies in the 3.6 micron images from the Spitzer Survey of Stellar Structure in Galaxies (S4G) and K_s-band images from the Near Infrared S0-Sa galaxy Survey (NIRS0S). We particularly connect properties of single exponential (type I), downbending double exponential (type II), and upbending double exponential (type III) disc profile types, to structural components of galaxies by using detailed morphological classifications, and size measurements of rings and lenses. We also study how the local environment of the galaxies affects the profile types by calculating parameters describing the environmental density and the tidal interaction strength. We find that in majority of type II profiles the break radius is connected with structural components such as rings, lenses, and spirals. The exponential disc sections of all three profile types, when considered separately, follow the disc scaling relations. However, the outer discs of type II, and the inner discs of type III, are similar in scalelength to the single exponential discs. Although the different profile types have similar mean environmental parameters, the scalelengths of the type III profiles show a positive correlation with the tidal interaction strength.

The dynamics and evolution of any galactic structure are strongly influenced by the properties of the orbits that constitute it. In this paper, we compare two orbit classification schemes, one by Laskar (NAFF) , and the other by Carpintero and Aguilar (CA), by applying both of them to orbits obtained by following individual particles in a numerical simulation of a barred galaxy. We find that, at least for our case and some provisos, the main frequencies calculated by the two methods are in good agreement: for $80\%$ of the orbits the difference between the results of the two methods is less than $5\%$ for all three main frequencies. However, it is difficult to evaluate the amount of regular or chaotic bar orbits in a given system. The fraction of regular orbits obtained by the NAFF method strongly depends on the critical frequency drift parameter, while in the CA method the number of fundamental frequencies strongly depends on the frequency difference parameter $L_{\rm r}$ and the maximum integer used for searching the linear independence of the fundamental frequencies. We also find that, for a given particle, in general the projection of its motion along the bar minor axis is more regular than the other two projections, while the projection along the intermediate axis is the least regular.

Magellanic spiral galaxies are a class of one-armed systems that often exhibit an offset stellar bar, and are rarely found around massive spiral galaxies. Using a set of N-body and hydrodynamic simulations we consider a dwarf-dwarf galaxy interaction as the driving mechanism for the formation of this peculiar class of systems. We investigate here the relation between the dynamical, stellar and gaseous disk center and the bar. In all our simulations the bar center always coincides with the dynamical center, while the stellar disk becomes highly asymmetric during the encounter causing the photometric center of the Magellanic galaxy disk to become mismatched with both the bar and the dynamical center. The disk asymmetries persist for almost 2 Gyrs, the time that it takes for the disk to be re-centered with the bar, and well after the companion has passed. This explains the nature of the offset bar found in many Magellanic-type galaxies, including the Large Magellanic Cloud (LMC) and NGC 3906. In particular, these results, once applied to the LMC, suggest that the dynamical center should reside in the bar center instead of the HI center as previously assumed, pointing to a variation in the current estimate of the north component of the LMC proper motion.

Orbits are the key building blocks of any density distribution and their study helps us understand the kinematical structure and the evolution of galaxies. Here we investigate orbits in a tidally induced bar of a dwarf galaxy, using an $N$-body simulation of an initially disky dwarf galaxy orbiting a Milky Way-like host. After the first pericenter passage, a tidally induced bar forms in the stellar component of the dwarf. The bar evolution is different than in isolated galaxies and our analysis focuses on the period before it buckles. We study the orbits in terms of their dominant frequencies, which we calculate in a Cartesian coordinate frame rotating with the bar. Apart from the well-known x$_1$ orbits we find many other types, mostly with boxy shapes of various degree of elongation. Some of them are also near-periodic, admitting frequency ratios of 4/3, 3/2 and 5/3. The box orbits have various degrees of vertical thickness but only a relatively small fraction of those have banana (i.e. smile/frown) or infinity-symbol shapes in the edge-on view. In the very center we also find orbits known from the potential of triaxial ellipsoids. The elongation of the orbits grows with distance from the center of the bar in agreement with the variation of the shape of the density distribution. Our classification of orbits leads to the conclusion that more than $80 \%$ of them have boxy shapes, while only $8 \%$ have shapes of classical x$_1$ orbits.

We present evidence of bar-induced secular evolution in galactic discs using 3.6 ${\mu m}$ images of nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find that among massive galaxies ($M_{star}/M_{sun}> 10^{10}$), longer bars tend to reside in inner discs having a flatter radial profile. Such galaxies show a light deficit in the disc surrounding the bar, within the bar radius and often show a $\Theta$-shaped morphology. We quantify this deficit and find that among all galaxies explored in this study (with $10^{9}<M_{star}/M_{sun}< 10^{11}$), galaxies with a stronger bar (i.e. longer and/or with a higher Bar/T) show a more pronounced deficit. We also examine simulation snapshots to confirm and extend results by Athanassoula and Misiriotis, showing that as bars evolve they become longer, while the light deficit in the disc becomes more pronounced. Theoretical studies have predicted that, as a barred galaxy evolves, the bar captures disc stars in its immediate neighbourhood so as to make the bar longer, stronger and thinner. Hence, we claim that the light deficit in the inner disc is produced by bars, which thus take part in shaping the mass distribution of their host galaxies.

In this letter we investigate the effect of boxy/peanut (b/p) bulges on bar-induced gas inflow to the central kiloparsec, which plays a crucial role on the evolution of disc galaxies. We carry out hydrodynamic gas response simulations in realistic barred galaxy potentials, including or not the geometry of a b/p bulge, to investigate the amount of gas inflow induced in the different models. We find that b/p bulges can reduce the gas inflow rate to the central kiloparsec by more than an order of magnitude, which leads to a reduction in the amount of gas available in the central regions. We also investigate the effect of the dark matter halo concentration on these results, and find that for maximal discs, the effect of b/p bulges on gas inflow remains significant. The reduced amount of gas reaching the central regions due to the presence of b/p bulges could have significant repercussions on the formation of discy- (pseudo-) bulges, on the amount of nuclear star formation and feedback, on the fuel reservoir for AGN activity, and on the overall secular evolution of galaxies.

Using N-body simulations we study the formation and evolution of tidally induced bars in disky galaxies in clusters. Our progenitor is a massive, late-type galaxy similar to the Milky Way, composed of an exponential disk and an NFW dark matter halo. We place the galaxy on four different orbits in a Virgo-like cluster and evolve it for 10 Gyr. As a reference case we also evolve the same model in isolation. Tidally induced bars form on all orbits soon after the first pericenter passage and survive until the end of the evolution. They appear earlier, are stronger, longer and have lower pattern speeds for tighter orbits. Only for the tightest orbit the properties of the bar are controlled by the orientation of the tidal torque from the cluster at pericenters. The mechanism behind the formation of the bars is the angular momentum transfer from the galaxy stellar component to its halo. All bars undergo extended periods of buckling instability that occur earlier and lead to more pronounced boxy/peanut shapes when the tidal forces are stronger. Using all simulation outputs of galaxies at different evolutionary stages we construct a toy model of the galaxy population in the cluster and measure the average bar strength and bar fraction as a function of clustercentric radius. Both are found to be mildly decreasing functions of radius. We conclude that tidal forces can trigger bar formation in cluster cores, but not in the outskirts, and thus cause larger concentrations of barred galaxies towards cluster center.

We have developed a computationally competitive N-body model of a previrialized aggregation of galaxies in a flat LambdaCDM universe to assess the role of the multiple mergers that take place during the formation stage of such systems in the configuration of the remnants assembled at their centres. An analysis of a suite of 48 simulations of low-mass forming groups (of about 1E13 solar masses) demonstrates that the gravitational dynamics involved in their hierarchical collapse is capable of creating realistic first-ranked galaxies without the aid of dissipative processes. Our simulations indicate that the brightest group galaxies (BGGs) constitute a distinct population from other group members, sketching a scenario in which the assembly path of these objects is dictated largely by the formation of their host system. We detect significant differences in the distribution of Sersic indices and total magnitudes, as well as a luminosity gap between BGGs and the next brightest galaxy that is positively correlated with the total luminosity of the parent group. Such gaps arise from both the grow of BGGs at the expense of lesser companions and the decrease in the relevance of second-ranked objects in equal measure. This results in a dearth of intermediate-mass galaxies which explains the characteristic central dip detected in their luminosity functions in dynamically young galaxy aggregations. The fact that the basic global properties of our BGGs define a thin mass fundamental plane strikingly similar to that followed giant early-type galaxies in the local universe reinforces confidence in the results obtained.

We present a pilot study on the nearby massive galaxy NGC 1291, in which we aim to constrain the dark matter in the inner regions, by obtaining a dynamical determination of the disc mass-to-light ratio (M/L). To this aim, we model the bar-induced dust lanes in the galaxy, using hydrodynamic gas response simulations. The models have three free parameters, the M/L of the disc, the bar pattern speed and the disc height function. We explore the parameter space to find the best fit models, i.e. those in which the morphology of the shocks in the gas simulations matches the observed dust lanes. The best-fit models suggest that the M/L of NGC 1291 agrees with that predicted by stellar population synthesis models in the near-infrared ($\approx$0.6\,$M_{\odot}/L_{\odot}$), which leads to a borderline maximum disc for this galaxy. The bar rotates fast, with corotation radius $\leq$ 1.4 times the bar length. Additionally, we find that the height function has a significant effect on the results, and can bias them towards lower or higher M/L.

Barred galaxies have interesting morphological features whose presence and properties set constraints on galactic evolution. Here we examine barlenses, i.e. lens-like components whose extent along the bar major axis is shorter than that of the bar and whose outline is oval or circular. We identify and analyse barlenses in $N$-body plus SPH simulations, compare them extensively with those from the NIRS0S (Near-IR S0 galaxy survey) and the S$^4$G samples (Spitzer Survey of Stellar Structure in Galaxies) and find very good agreement. We observe barlenses in our simulations from different viewing angles. This reveals that barlenses are the vertically thick part of the bar seen face-on, i.e. a barlens seen edge-on is a boxy/peanut/X bulge. In morphological studies, and in the absence of kinematics or photometry, a barlens, or part of it, may be mistaken for a classical bulge. Thus the true importance of classical bulges, both in numbers and mass, is smaller than currently assumed, which has implications for galaxy formation studies. Finally, using the shape of the isodensity curves, we propose a rule of thumb for measuring the barlens extent along the bar major axis of moderately inclined galaxies, thus providing an estimate of which part of the bar is thicker.

Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms subsequently and gradually in the remnant from the gas accreted mainly from the halo. It is vertically thin and well extended in its equatorial plane. A bar starts forming before the disk is fully in place, contrary to what is assumed in idealised simulations of isolated bar-forming galaxies. It has morphological features such as ansae and boxy/peanut bulges. Stars of different ages populate different parts of the box/peanut. A disky pseudobulge forms also, so that by the end of the simulation, all three types of bulges coexist. The oldest stars are found in the classical bulge, followed by those of the thick disk, then by those in the thin disk. The youngest stars are in the spiral arms and the disky pseudobulge. The disk surface density profiles are of type II (exponential with downbending), and the circular velocity curves are flat and show that the disks are submaximum in these examples: two clearly so and one near-borderline between maximum and submaximum. On average, only roughly between 10 and 20% of the stellar mass is in the classical bulge of the final models, i.e. much less than in previous simulations.

We present a detailed study of the shape of the innermost part of the rotation curves of a sample of 29 nearby spiral galaxies, based on high angular and spectral resolution kinematic Halpha Fabry-Perot observations. In particular, we quantify the steepness of the rotation curve by measuring its slope dRvc(0). We explore the relationship between the inner slope and several galaxy parameters, such as stellar mass, maximum rotational velocity, central surface brightness ({\mu}0), bar strength and bulge-to-total ratio. Even with our limited dynamical range, we find a trend for low-mass galaxies to exhibit shallower rotation curve inner slopes than high-mass galaxies, whereas steep inner slopes are found exclusively in high-mass galaxies. This trend may arise from the relationship between the total stellar mass and the mass of the bulge, which are correlated among them. We find a correlation between the inner slope of the rotation curve and the morphological T-type, complementary to the scaling relation between dRvc(0) and {\mu}0 previously reported in the literature. Although we find that the inner slope increases with the Fourier amplitude A2 and decreases with the bar torque Qb, this may arise from the presence of the bulge implicit in both A2 and Qb. As previously noted in the literature, the more compact the mass in the central parts of a galaxy (more concretely, the presence of a bulge), the steeper the inner slopes. We conclude that the baryonic matter dominates the dynamics in the central parts of our sample galaxies.

We present the stellar kinematics across the Galactic bulge and into the disk at positive longitudes from the SDSS-III APOGEE spectroscopic survey of the Milky Way. APOGEE includes extensive coverage of the stellar populations of the bulge along the mid-plane and near-plane regions. From these data, we have produced kinematic maps of 10,000 stars across longitudes 0 deg < l < 65 deg, and primarily across latitudes of |b| < 5 deg in the bulge region. The APOGEE data reveal that the bulge is cylindrically rotating across all latitudes and is kinematically hottest at the very centre of the bulge, with the smallest gradients in both kinematic and chemical space inside the inner-most region (l,|b|) < (5,5) deg. The results from APOGEE show good agreement with data from other surveys at higher latitudes and a remarkable similarity to the rotation and dispersion maps of barred galaxies viewed edge on. The thin bar that is reported to be present in the inner disk within a narrow latitude range of |b| < 2 deg appears to have a corresponding signature in [Fe/H] and [alpha/Fe]. Stars with [Fe/H] > -0.5 have dispersion and rotation profiles that are similar to that of N-body models of boxy/peanut bulges. There is a smooth kinematic transition from the thin bar and boxy bulge (l,|b|) < (15,12) deg out into the disk for stars with [Fe/H] > -1.0, and the chemodynamics across (l,b) suggests the stars in the inner Galaxy with [Fe/H] > -1.0 have an origin in the disk.

Using 3.6 and 4.5$\mu$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars.

Using numerical analysis of fundamental frequencies we study the orbital structure of a tidally induced bar formed in a simulated dwarf galaxy orbiting a Milky Way-like host. We find that only about 10% of stars have frequencies compatible with x1 orbits, the classical periodic orbits in a barred potential. The rest of the stars follows box orbits parallel to the bar, with varying degree of elongation.

The strong time-dependence of the dynamics of galactic bars yields a complex and rapidly evolving distribution of dense gas and star forming regions. Although bars mainly host regions void of any star formation activity, their extremities can gather the physical conditions for the formation of molecular complexes and mini-starbursts. Using a sub-parsec resolution hydrodynamical simulation of a Milky Way-like galaxy, we probe these conditions to explore how and where bar (hydro-)dynamics favours the formation or destruction of molecular clouds and stars. The interplay between the kpc-scale dynamics (gas flows, shear) and the parsec-scale (turbulence) is key to this problem. We find a strong dichotomy between the leading and trailing sides of the bar, in term of cloud fragmentation and in the age distribution of the young stars. After orbiting along the bar edge, these young structures slow down at the extremities of the bar, where orbital crowding increases the probability of cloud-cloud collision. We find that such events increase the Mach number of the cloud, leading to an enhanced star formation efficiency and finally the formation of massive stellar associations, in a fashion similar to galaxy-galaxy interactions. We highlight the role of bar dynamics in decoupling young stars from the clouds in which they form, and discuss the implications on the injection of feedback into the interstellar medium, in particular in the context of galaxy formation.

Resonance rings are the consequence of secular evolution processes that redistribute material and angular momentum in discs. We produced a Catalogue and an Atlas of the rings detected in the Spitzer Survey of Stellar Structure in Galaxies (S4G) and to conduct a statistical study of the data in the Catalogue. We traced the contours of rings previously identified by Buta et al. (in preparation) and fitted them with ellipses. We found the orientation of bars by studying the galaxy ellipse fits from S4G's Pipeline 4. We used the galaxy orientation data obtained by S4G's Pipeline 4 to obtain intrinsic ellipticities and orientations of rings and the bars. ARRAKIS contains data on 724 ringed galaxies in the S4G. The frequency of resonance rings in the S4G is of 16+-1% and 35+-1% for outer and inner features, respectively. Outer rings are mostly found in Hubble stages -1<=T<=4. Inner rings are found in a distribution that covers the range -1<=T<=7. We confirm that outer rings have two preferred orientations, parallel and perpendicular to the bar. We confirm a tendency for inner rings to be oriented parallel to the bar, but we find that a significant fraction (~50%) of them have random orientations with respect to the bar. These misaligned inner rings are mostly found in late-type galaxies (T>=4). This may be due to spiral modes decoupled from the bar dominating the Fourier amplitude spectrum at the radius of the inner ring. We find that the fraction of barred galaxies hosting outer (inner) rings is ~1.7 times (~1.3 times) that in unbarred galaxies. The fact that rings are only mildly favoured by bars suggests that those in unbarred galaxies either formed due to weak departures from the axisymmetry of the galactic potential or that they are born because of bars that have been destroyed after the ring formation.

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$\mu m$ IRAC band has been exploited to this effect, but such mass estimates can be biased by dust emission. We present our pipeline to reveal the old stellar flux at 3.6$\mu m$ and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G). This survey consists of images in two infrared bands (3.6 and 4.5$\mu m$), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. (2012) to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6$\mu m$ flux. We exclude from our ICA analysis galaxies with low signal-to-noise ratio (S/N < 10) and those with original [3.6]-[4.5] colors compatible with an old stellar population, indicative of little dust emission (mostly early Hubble types, which can directly provide good mass maps). For the remaining 1251 galaxies to which ICA was successfully applied, we find that as much as 10-30% of the total light at 3.6$\mu m$ typically originates from dust, and locally it can reach even higher values. This contamination fraction shows a correlation with specific star formation rates, confirming that the dust emission that we detect is related to star formation. Additionally, we have used our large sample of mass estimates to calibrate a relationship of effective mass-to-light ratio ($M/L$) as a function of observed [3.6]-[4.5] color: $\log(M/L)=-0.339 (\pm 0.057) \times ([3.6]-[4.5]) -0.336 (\pm 0.002)$. Our final pipeline products have been made public through IRSA, providing the astronomical community with an unprecedentedly large set of stellar mass maps ready to use for scientific applications.

We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.

We study the role of radial migration of stars on the chemical evolution of the Milky Way disk. In particular, we are interested in the impact of that process on the local properties of the disk (age-metallicity relation and its dispersion, metallicity distribution, evolution of abundance ratios) and on the morphological properties of the resulting thick and thin disks.We use a model with several new or up-dated ingredients: atomic and molecular gas phases, star formation depending on molecular gas, yields from the recent homogeneous grid provided by Nomoto et al. (2013), observationally inferred SNIa rates. We describe radial migration with parametrised time- and radius-dependent diffusion coefficients, based on the analysis of a N-body+SPH simulation. We also consider parametrised radial gas flows, induced by the action of the Galactic bar. Our model reproduces well the present day values of most of the main global observables of the MW disk and bulge, and also the observed "stacked" evolution of MW-type galaxies from van Dokkum et al. (2013). The azimuthally averaged radial velocity of gas inflow is constrained to less than a few tenths of km/s. Radial migration is constrained by the observed dispersion in the age-metallicity relation. Assuming that the thick disk is the oldest (>9 Gyr) part of the disk, we find that the adopted radial migration scheme can reproduce quantitatively the main local properties of the thin and thick disk. The thick disk extends up to ~11 kpc and has a scale length of 1.8 kpc, considerably shorter than the thin disk, because of the inside-out formation scheme. We also show how, in this framework, current and forthcoming spectroscopic observations can constrain the nucleosynthesis yields of massive stars for the metallicity range of 0.1 solar to 2-3 solar.

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is a volume, magnitude, and size-limited survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5um. In this paper we describe our surface photometry pipeline and showcase the associated data products that we have released to the community. We also identify the physical mechanisms leading to different levels of central stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local stellar mass-size relation at 3.6um for galaxies of different morphologies. Our radial profiles reach stellar mass surface densities below 1 Msun pc-2. Given the negligible impact of dust and the almost constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and shapes, and concentration indices. These and other data products from our various pipelines (science-ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly accessed at IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/).

We present an implementation of the Schwarzschild orbit superposition method which can be used for constructing self-consistent equilibrium models of barred or non-barred disc galaxies, or of elliptical galaxies with figure rotation. This is a further development of the publicly available code SMILE; its main improvements include a new efficient representation of an arbitrary gravitational potential using two-dimensional spline interpolation of Fourier coefficients in the meridional plane, as well as the ability to deal with rotation of the density profile and with multicomponent mass models. We compare several published methods for constructing composite axisymmetric disc--bulge--halo models and demonstrate that our code produces the models that are closest to equilibrium. We also apply it to create models of triaxial elliptical galaxies with cuspy density profiles and figure rotation, and find that such models can be found and are stable over many dynamical times in a wide range of pattern speeds and angular momenta, covering both slow- and fast-rotator classes. We then attempt to create models of strongly barred disc galaxies, using an analytic three-component potential, and find that it is not possible to make a stable dynamically self-consistent model for this density profile. Finally, we take snapshots of two N-body simulations of barred disc galaxies embedded in nearly-spherical haloes, and construct equilibrium models using only information on the density profile of the snapshots. We demonstrate that such reconstructed models are in near-stationary state, in contrast with the original N-body simulations, one of which displayed significant secular evolution.

Selecting centrally quiescent galaxies from the Sloan Digital Sky Survey (SDSS) to create high signal-to-noise (>100) stacked spectra with minimal emission line contamination, we accurately and precisely model the central stellar populations of barred and unbarred quiescent disk galaxies. By splitting our sample by redshift, we can use the fixed size of the SDSS fiber to model the stellar populations at different radii within galaxies. At 0.02<z<0.04, the SDSS fiber radius corresponds to ~1 kpc, which is the typical half-light radii of both classical bulges and disky pseudobulges. Assuming that the SDSS fiber primarily covers the bulges at these redshifts, our analysis shows that there are no significant differences in the stellar populations, i.e., stellar age, [Fe/H], [Mg/Fe], and [N/Fe], of the bulges of barred vs. unbarred quiescent disk galaxies. Modeling the stellar populations at different redshift intervals from z=0.020 to z=0.085 at fixed stellar masses produces an estimate of the stellar population gradients out to about half the typical effective radius of our sample, assuming null evolution over this ~1 Gyr epoch. We find that there are no noticeable differences in the slopes of the azimuthally averaged gradients of barred vs. unbarred quiescent disk galaxies. These results suggest that bars are not a strong influence on the chemical evolution of quiescent disk galaxies.

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.

We study the imprints of boxy/peanut structures on the 2D line-of-sight kinematics of simulated disk galaxies. The models under study belong to a family with varying initial gas fraction and halo triaxiality, plus few other control runs with different structural parameters; the kinematic information was extracted using the Voronoi-binning technique and parametrised up to the fourth order of a Gauss-Hermite series. Building on a previous work for the long-slit case, we investigate the 2D kinematic behaviour in the edge-on projection as a function of the boxy/peanut strength and position angle; we find that for the strongest structures the highest moments show characteristic features away from the midplane in a range of position angles. We also discuss the masking effect of a classical bulge and the ambiguity in discriminating kinematically this spherically-symmetric component from a boxy/peanut bulge seen end-on. Regarding the face-on case, we extend existing results to encompass the effect of a second buckling and find that this phenomenon spurs an additional set of even deeper minima in the fourth moment. Finally, we show how the results evolve when inclining the disk away from perfectly edge-on and face-on. The behaviour of stars born during the course of the simulations is discussed and confronted to that of the pre-existing disk. The general aim of our study is providing a handle to identify boxy/peanut structure and their properties in latest generation IFU observations of nearby disk galaxies.

Bars have a complex three-dimensional shape. In particular their inner part is vertically much thicker than the parts further out. Viewed edge-on, the thick part of the bar is what is commonly known as a boxy-, peanut- or X- bulge and viewed face-on it is referred to as a barlens. These components are due to disc and bar instabilities and are composed of disc material. I review here their formation, evolution and dynamics, using simulations, orbital structure theory and comparisons to observations.

We examine the effects that the modelling of a Boxy/Peanut (B/P) bulge will have on the estimates of the stellar gravitational potential, forces, orbital structure and bar strength of barred galaxies. We present a method for obtaining the potential of disc galaxies from surface density images, assuming a vertical density distribution (height function), which is let to vary with position, thus enabling it to represent the geometry of a B/P. We construct a B/P height function after the results from a high-resolution, N-body+SPH simulation of an isolated galaxy and compare the resulting dynamical model to those obtained with the commonly used, position-independent "flat" height functions. We show that methods that do not allow for a B/P can induce errors in the forces in the bar region of up to 40% and demonstrate that this has a significant impact on the orbital structure of the model, which in turn determines its kinematics and morphology. Furthermore, we show that the bar strength is reduced in the presence of a B/P. We conclude that neglecting the vertical extent of a B/P can introduce considerable errors in the dynamical modelling. We also examine the errors introduced in the model due to uncertainties in the parameters of the B/P and show that even for generous but realistic values of the uncertainties, the error will be noticeably less than that of not modelling a B/P bulge at all.

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is the largest available database of deep, homogeneous middle-infrared (mid-IR) images of galaxies of all types. The survey, which includes 2352 nearby galaxies, reveals galaxy morphology only minimally affected by interstellar extinction. This paper presents an atlas and classifications of S4G galaxies in the Comprehensive de Vaucouleurs revised Hubble-Sandage (CVRHS) system. The CVRHS system follows the precepts of classical de Vaucouleurs (1959) morphology, modified to include recognition of other features such as inner, outer, and nuclear lenses, nuclear rings, bars, and disks, spheroidal galaxies, X patterns and box/peanut structures, OLR subclass outer rings and pseudorings, bar ansae and barlenses, parallel sequence late-types, thick disks, and embedded disks in 3D early-type systems. We show that our CVRHS classifications are internally consistent, and that nearly half of the S4G sample consists of extreme late-type systems (mostly bulgeless, pure disk galaxies) in the range Scd-Im. The most common family classification for mid-IR types S0/a to Sc is SA while that for types Scd to Sm is SB. The bars in these two type domains are very different in mid-IR structure and morphology. This paper examines the bar, ring, and type classification fractions in the sample, and also includes several montages of images highlighting the various kinds of "stellar structures" seen in mid-IR galaxy morphology.

We present a new study investigating whether active galactic nuclei (AGN) beyond the local universe are preferentially fed via large-scale bars. Our investigation combines data from Chandra and Galaxy Zoo: Hubble (GZH) in the AEGIS, COSMOS, and GOODS-S surveys to create samples of face-on, disc galaxies at 0.2 < z < 1.0. We use a novel method to robustly compare a sample of 120 AGN host galaxies, defined to have 10^42 erg/s < L_X < 10^44 erg/s, with inactive control galaxies matched in stellar mass, rest-frame colour, size, Sersic index, and redshift. Using the GZH bar classifications of each sample, we demonstrate that AGN hosts show no statistically significant enhancement in bar fraction or average bar likelihood compared to closely-matched inactive galaxies. In detail, we find that the AGN bar fraction cannot be enhanced above the control bar fraction by more than a factor of two, at 99.7% confidence. We similarly find no significant difference in the AGN fraction among barred and non-barred galaxies. Thus we find no compelling evidence that large-scale bars directly fuel AGN at 0.2<z<1.0. This result, coupled with previous results at z=0, implies that moderate-luminosity AGN have not been preferentially fed by large-scale bars since z=1. Furthermore, given the low bar fractions at z>1, our findings suggest that large-scale bars have likely never directly been a dominant fueling mechanism for supermassive black hole growth.

We study the role of radial motions of stars and gas on the evolution of abundance profiles in the Milky Way disk. We investigate, in a parametrized way, the impact of radial flows of gas and radial migration of stars induced mainly by the Galactic bar and its iteraction with the spiral arms. We use a model with several new or up-dated ingredients (atomic and molecular gas phases, star formation depending on molecular gas, recent sets of metallicity-dependent stellar yields from H to Ni, observationally inferred SNIa rates), which reproduces well most global and local observables of the Milky Way. We obtain abundance profiles flattening both in the inner disk (because of radial flows) and in the outer disk (because of the adopted star formation law). The gas abundance profiles flatten with time, but the corresponding stellar profiles appear to be steeper for younger stars, because of radial migration. We find a correlation between the stellar abundance profiles and O/Fe, which is a proxy for stellar age. Our final abundance profiles are in overall agreement with observations, but slightly steeper (by 0.01-0.02 dex/kpc) for elements above S. We find an interesting "odd-even effect" in the behaviour of the abundance profiles (steeper slopes for odd elements) for all sets of stellar yields; however, this behaviour does not appear in observations, suggesting that the effect is, perhaps, overestimated in current stellar nucleosynthesis calculations.

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

We have measured the radial light profiles and global shapes of bars using two-dimensional 3.6 $\mu m $ image decompositions for 144 face-on barred galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The bar surface brightness profile is correlated with the stellar mass and bulge-to-total (B/T) ratio of their host galaxies. Bars in massive and bulge-dominated galaxies (B/T$>$0.2) show a flat profile, while bars in less massive, disk-dominated galaxies (B/T$\sim$0) show an exponential, disk-like profile with a wider spread in the radial profile than in the bulge-dominated galaxies. The global two-dimensional shapes of bars, however, are rectangular/boxy, independent of the bulge or disk properties. We speculate that because bars are formed out of disk, bars initially have an exponential (disk-like) profile which evolves over time, trapping more stars into the boxy bar orbits. This leads bars to become stronger and have flatter profiles. The narrow spread of bar radial profiles in more massive disks suggests that these bars formed earlier (z$>$1), while the disk-like profiles and a larger spread in the radial profile in less massive systems imply a later and more gradual evolution, consistent with the cosmological evolution of bars inferred from observational studies. Therefore, we expect that the flatness of the bar profile can be used as a dynamical age indicator of the bar to measure the time elapsed since the bar formation. We argue that cosmic gas accretion is required to explain our results on bar profile and the presence of gas within the bar region.

Using 3.6$\mu$m images of 97 early-type galaxies, we develop and verify methodology to measure globular cluster populations from the S$^4$G survey images. We find that 1) the ratio, T$_{\rm N}$, of the number of clusters, N$_{\rm CL}$, to parent galaxy stellar mass, M$_*$, rises weakly with M$_*$ for early-type galaxies with M$_* > 10^{10}$ M$_\odot$ when we calculate galaxy masses using a universal stellar initial mass function (IMF), but that the dependence of T$_{\rm N}$ on M$_*$ is removed entirely once we correct for the recently uncovered systematic variation of IMF with M$_*$, and 2) for M$_* < 10^{10}$ M$_\odot$ there is no trend between N$_{\rm CL}$ and M$_*$, the scatter in T$_{\rm N}$ is significantly larger (approaching 2 orders of magnitude), and there is evidence to support a previous, independent suggestion of two families of galaxies. The behavior of N$_{\rm CL}$ in the lower mass systems is more difficult to measure because these systems are inherently cluster poor, but our results may add to previous evidence that large variations in cluster formation and destruction efficiencies are to be found among low mass galaxies. The average fraction of stellar mass in clusters is $\sim$ 0.0014 for M$_* > 10^{10}$ M$_\odot$ and can be as large as $\sim 0.02$ for less massive galaxies. These are the first results from the S$^4$G sample of galaxies, and will be enhanced by the sample of early-type galaxies now being added to S$^4$G and complemented by the study of later type galaxies within S$^4$G.

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.

Boxy/Peanut/X-shaped (B/P/X) bulges are studied using the 3.6 mum images from the Spitzer Survey of Stellar Structure in Galaxies, and the K_s-band images from the Near-IR S0 galaxy Survey. They are compared with the properties of barlenses, defined as lens-like structures embedded in bars, with sizes of about 50% of bars and axial ratios of 0.6-0.9. Based on observations and recent simulation models we show evidence that barlenses are the more face-on counterparts of B/P/X-shaped bulges. Using unsharp masks 18 new X-shaped structures are identified, covering a large range of galaxy inclinations. The similar masses and red B-3.6 mum colors of the host galaxies, and the fact that the combined axial ratio distribution of the host galaxy disks is flat, support the interpretation that barlenses and X-shapes are physically the same phenomenon. In Hubble types -3<T<2 even half of the bars contain either a barlens or an X-shaped structure. Our detailed 2D multi-component decompositions for 29 galaxies, fitting the barlens/X-shape with a separate component, indicate very small or non-existent classical bulges. Taking into account that the structures we study have similar host galaxy masses as the Milky Way (MW), our results imply that MW mass galaxies with no significant classical bulges are common in the nearby Universe.

We present a catalogue and images of visually detected features, such as asymmetries, extensions, warps, shells, tidal tails, polar rings, and obvious signs of mergers or interactions, in the faint outer regions (at and outside of R_25) of nearby galaxies. This catalogue can be used in future quantitative studies that examine galaxy evolution due to internal and external factors. We are able to reliably detect outer region features down to a brightness level of 0.03 MJy/sr per pixel at 3.6 microns in the Spitzer Survey of Stellar Structure in Galaxies (S4G). We also tabulate companion galaxies. We find asymmetries in the outer isophotes in 22+/-1 per cent of the sample. The asymmetry fraction does not correlate with galaxy classification as an interacting galaxy or merger remnant, or with the presence of companions. We also compare the detected features to similar features in galaxies taken from cosmological zoom re-simulations. The simulated images have a higher fraction (33 per cent) of outer disc asymmetries, which may be due to selection effects and an uncertain star formation threshold in the models. The asymmetries may have either an internal (e.g., lopsidedness due to dark halo asymmetry) or external origin.

We demonstrate that failure to properly account for stellar evolution can bias results in determinations of the spatial morphology of Galactic bulge stars, focusing on the question of whether or not the X-shape is more pronounced among the more metal-rich stars than among the metal-poor stars. We argue that this trend, a result recently claimed by three separate groups, may have been overestimated as it is relatively easier to detect a bimodality in the distance distribution function at higher metallicities. This is due to three factors. First, the intrinsic colour of red clump and red giant stars vary with metallicity, at the level d(V-I)_{RC}/d\[M/H] ~ 0.25 mag dex^{-1}, and thus the ratio of red clump to red giant stars within a spectroscopic sample will depend on the photometric selection of any investigation. Second, the duration of ascent of the red giant branch goes down and the red clump lifetime goes up as metallicity increases, which has the effect of increasing the ratio of red clump to red giant stars by as much as ~33% over the range of the bulge metallicity-distribution function. Finally, over the same metallicity interval, the effective number of red giant branch bump stars is predicted to increase by ~200%, and their presence becomes degenerate with the observational parameters of the two red clumps, creating an illusory increase in signal-to-noise for a second peak in the distance modulus distribution.

HI kinematic asymmetries are common in late-type galaxies irrespective of environment, although the amplitudes are strikingly low in isolated galaxies. As part of our studies of the HI morphology and kinematics in isolated late-type galaxies we have chosen several very isolated galaxies from the AMIGA sample for HI mapping. Here we present GMRT 21-cm HI line mapping of CIG 340 which was selected because its integrated HI spectrum has a very symmetric profile, Aflux = 1.03. Optical images of the galaxy hinted at a warped disk in contrast to the symmetric integrated HI spectrum profile. Our aim is to determine the extent to which the optical asymmetry is reflected in the resolved HI morphology and kinematics. GMRT observations reveal significant HI morphological asymmetries in CIG 340 despite it's overall symmetric optical form and highly symmetric HI spectrum. The most notable HI features are: 1) a warp in the HI disk (with an optical counterpart), 2) the HI north/south flux ratio = 1.32 is much larger than expected from the integrated HI spectrum profile and 3) a ~ 45" (12 kpc) HI extension, containing ~ 6% of the detected HI mass on the northern side of the disk. We conclude that in isolated galaxies a highly symmetric HI spectrum can mask significant HI morphological asymmetries. The northern HI extension appears to be the result of a recent perturbation (10^8 yr), possibly by a satellite which is now disrupted or projected within the disk. This study provides an important step in our ongoing program to determine the predominant source of HI asymmetries in isolated galaxies. For CIG 340 the isolation from major companions, symmetric HI spectrum, optical morphology and interaction timescales have allowed us to narrow the possible causes the HI asymmetries and identify tests to further constrain the source of the asymmetries.

We combine data from the Spitzer Survey for Stellar Structure in Galaxies (S$^4$G), a recently calibrated empirical stellar mass estimator from Eskew et al., and an extensive database of HI spectral line profiles to examine the baryonic Tully-Fisher (BTF) relation. We find 1) that the BTF has lower scatter than the classic Tully-Fisher (TF) relation and is better described as a linear relationship, confirming similar previous results, 2) that the inclusion of a radial scale in the BTF decreases the scatter but only modestly, as seen previously for the TF relation, and 3) that the slope of the BTF, which we find to be $3.5\pm 0.2$ ($\Delta$ log $M_{baryon}/\Delta$ log $v_c$), implies that on average a nearly constant fraction ($\sim 0.4$) of all baryons expected to be in a halo are "condensed" onto the central region of rotationally supported galaxies. The condensed baryon fraction, $M_{baryon}/M_{total}$, is, to our measurement precision, nearly independent of galaxy circular velocity (our sample spans circular velocities, $v_c$, between 60 and 250 km s$^{-1}$, but is extended to $v_c\sim 10$ km s$^{-1}$ using data from the literature). The observed galaxy-to-galaxy scatter in this fraction is generally $\le$ a factor of 2 despite fairly liberal selection criteria. These results imply that cooling and heating processes, such as cold vs. hot accretion, mass loss due to stellar winds, and AGN driven feedback, to the degree that they affect the global galactic properties involved in the BTF, are independent of halo mass for galaxies with $10 < v_c < 250$ km/s and typically introduce no more than a factor of two range in the resulting $M_{baryon}/M_{total}$. Recent simulations by Aumer et al. of a small sample of disk galaxies are in excellent agreement with our data, suggesting that current simulations are capable of reproducing the global properties of individual disk galaxies.

We present a new approach for estimating the 3.6 micron stellar mass-to-light ratio in terms of the [3.6]-[4.5] colors of old stellar populations. Our approach avoids several of the largest sources of uncertainty in existing techniques. By focusing on mid-IR wavelengths, we gain a virtually dust extinction-free tracer of the old stars, avoiding the need to adopt a dust model to correctly interpret optical or optical/NIR colors normally leveraged to assign M/L. By calibrating a new relation between NIR and mid-IR colors of GLIMPSE giant stars we also avoid discrepancies in model predictions for the [3.6]-[4.5] colors of old stellar populations due to uncertainties in molecular line opacities. We find that the [3.6]-[4.5] color, which is driven primarily by metallicity, provides a tight constraint on M/L_3.6, which varies intrinsically less than at optical wavelengths. The uncertainty on M/L_3.6 of ~0.07 dex due to unconstrained age variations marks a significant improvement on existing techniques for estimating the stellar M/L with shorter wavelength data. A single M/L_3.6=0.6 (assuming a Chabrier IMF), independent of [3.6]-[4.5] color, is also feasible as it can be applied simultaneously to old, metal-rich and young, metal-poor populations, and still with comparable (or better) accuracy (~0.1 dex) as alternatives. We expect our M/L_3.6 to be optimal for mapping the stellar mass distributions in S4G galaxies, for which we have developed an Independent Component Analysis technique to first isolate the old stellar light at 3.6 micron from non-stellar emission (e.g. hot dust and the 3.3 PAH feature). Our estimate can also be used to determine the fractional contribution of non-stellar emission to global (rest-frame) 3.6 micron fluxes, e.g. in WISE imaging, and establishes a reliable basis for exploring variations in the stellar IMF.

We have performed two-dimensional multicomponent decomposition of 144 local barred spiral galaxies using 3.6 $\mu {\rm m}$ images from the Spitzer Survey of Stellar Structure in Galaxies. Our model fit includes up to four components (bulge, disk, bar, and a point source) and, most importantly, takes into account disk breaks. We find that ignoring the disk break and using a single disk scale length in the model fit for Type II (down-bending) disk galaxies can lead to differences of 40% in the disk scale length, 10% in bulge-to-total luminosity ratio (B/T), and 25% in bar-to-total luminosity ratios. We find that for galaxies with B/T $\geq$ 0.1, the break radius to bar radius, $r_{\rm br}/R_{\rm bar}$, varies between 1 and 3, but as a function of B/T the ratio remains roughly constant. This suggests that in bulge-dominated galaxies the disk break is likely related to the outer Lindblad Resonance (OLR) of the bar, and thus moves outwards as the bar grows. For galaxies with small bulges, B/T $<$ 0.1, $r_{\rm br}/R_{\rm bar}$ spans a wide range from 1 to 6. This suggests that the mechanism that produces the break in these galaxies may be different from that in galaxies with more massive bulges. Consistent with previous studies, we conclude that disk breaks in galaxies with small bulges may originate from bar resonances that may be also coupled with the spiral arms, or be related to star formation thresholds.

`Conspiracy' between the dark and the baryonic mater prohibits an unambiguous decomposition of disc galaxy rotation curves into the corresponding components. Several methods have been proposed to counter this difficulty, but their results are widely discrepant. In this paper, I revisit one of these methods, which relies on the relation between the halo density and the decrease of the bar pattern speed. The latter is routinely characterised by the ratio ${\cal R}$ of the corotation radius $R_{CR}$ to the bar length $L_b$, ${\cal R}=R_{CR}/L_b$. I use a set of $N$-body+SPH simulations, including sub-grid physics, whose initial conditions cover a range of gas fractions and halo shapes. The models, by construction, have roughly the same azimuthally averaged circular velocity curve and halo density and they are all submaximal, i.e. according to previous works they are expected to have all roughly the same ${\cal R}$ value, well outside the fast bar range (1.2 $\pm$ 0.2). Contrary to these expectations, however, these simulations end up having widely different ${\cal R}$ values, either within the fast bar range, or well outside it. This shows that the ${\cal R}$ value can not constrain the halo density, nor determine whether galactic discs are maximal or submaximal. I argue that this is true even for early type discs (S0s and Sas).

In this paper, we use the Galaxy Zoo 2 dataset to study the behavior of bars in disk galaxies as a function of specific star formation rate (SSFR), and bulge prominence. Our sample consists of 13,295 disk galaxies, with an overall (strong) bar fraction of $23.6\pm 0.4\%$, of which 1,154 barred galaxies also have bar length measurements. These samples are the largest ever used to study the role of bars in galaxy evolution. We find that the likelihood of a galaxy hosting a bar is anti-correlated with SSFR, regardless of stellar mass or bulge prominence. We find that the trends of bar likelihood and bar length with bulge prominence are bimodal with SSFR. We interpret these observations using state-of-the-art simulations of bar evolution which include live halos and the effects of gas and star formation. We suggest our observed trends of bar likelihood with SSFR are driven by the gas fraction of the disks; a factor demonstrated to significantly retard both bar formation and evolution in models. We interpret the bimodal relationship between bulge prominence and bar properties as due to the complicated effects of classical bulges and central mass concentrations on bar evolution, and also to the growth of disky pseudobulges by bar evolution. These results represent empirical evidence for secular evolution driven by bars in disk galaxies. This work suggests that bars are not stagnant structures within disk galaxies, but are a critical evolutionary driver of their host galaxies in the local universe ($z<1$).

Star-forming regions that are visible at 3.6 microns and Halpha but not in the u,g,r,i,z bands of the Sloan Digital Sky survey (SDSS), are measured in five nearby spiral galaxies to find extinctions averaging ~3.8 mag and stellar masses averaging ~5x10^4 Msun. These regions are apparently young star complexes embedded in dark filamentary shock fronts connected with spiral arms. The associated cloud masses are ~10^7 Msun. The conditions required to make such complexes are explored, including gravitational instabilities in spiral shocked gas and compression of incident clouds. We find that instabilities are too slow for a complete collapse of the observed spiral filaments, but they could lead to star formation in the denser parts. Compression of incident clouds can produce a faster collapse but has difficulty explaining the semi-regular spacing of some regions along the arms. If gravitational instabilities are involved, then the condensations have the local Jeans mass. Also in this case, the near-simultaneous appearance of equally spaced complexes suggests that the dust lanes, and perhaps the arms too, are relatively young.

We study radial migration and chemical evolution in a bar-dominated disk galaxy, by analyzing the results of a fully self-consistent, high resolution N-body+SPH simulation. We find different behaviours for gas and star particles. Gas within corotation is driven in the central regions by the bar, where it forms a pseudo-bulge (disky-bulge), but it undergoes negligible radial displacement outside the bar region. Stars undergo substantial radial migration at all times, caused first by transient spiral arms and later by the bar. Despite the important amount of radial migration occurring in our model, its impact on the chemical properties is limited. The reason is the relatively flat abundance profile, due to the rapid early evolution of the whole disk. We show that the implications of radial migration on chemical evolution can be studied to a good accuracy by post-processing the results of the N-body+SPH calculation with a simple chemical evolution model having detailed chemistry and a parametrized description of radial migration. We find that radial migration impacts on chemical evolution both directly (by moving around the long-lived agents of nucleosynthesis, like e.g. SNIa or AGB stars, and thus altering the abundance profiles of the gas) and indirectly (by moving around the long-lived tracers of chemical evolution and thus affecting stellar metallicity profiles, local age-metallicity relations and metallicity distributions of stars, etc.).

The morphology of galaxies can be quantified to some degree using a set of scale-invariant parameters. Concentration (C), Asymmetry (A), Smoothness (S), the Gini index (G), relative contribution of the brightest pixels to the second order moment of the flux (M20), ellipticity (E), and the Gini index of the second order moment (GM) have all been applied to morphologically classify galaxies at various wavelengths. Here we present a catalog of these parameters for the Spitzer Survey of Stellar Structure in Galaxies (S4G), a volume-limited near-infrared imaging survey of nearby galaxies using the 3.6 and 4.5 micron channels of the IRAC camera. Our goal is to provide a reference catalog of near-infrared quantified morphology for high-redshift studies and galaxy evolution models with enough detail to resolve stellar mass morphology. We explore where normal, non-interacting galaxies -those typically found on the Hubble tuning fork- lie in this parameter space and show that there is a tight relation between Concentration and M20 for normal galaxies. M20 can be used to classify galaxies into earlier and later types (e.g., to separate spirals from irregulars). Several criteria using these parameters exist to select systems with a disturbed morphology, i.e., those that appear to be undergoing a tidal interaction. We examine the applicability of these criteria to Spitzer near-infrared imaging. We find that four relations, based on the parameters A&S, G&M20, GM, and C&M20, respectively, select outliers in the morphological parameter space, but each selects different subsets of galaxies. Two criteria (GM > 0.6, G > -0.115 x M20 + 0.384) seem most appropriate to identify possible mergers and the merger fraction in near-infrared surveys. We find no strong relation between lopsidedness and most of these morphological parameters, except for a weak dependence of lopsidedness on Concentration and M20.

Dynamical simulations are a fundamental tool for studying the secular evolution of disc galaxies. Even at their maximum resolution, they still follow a limited number of particles and typically resolve scales of the order of a few tens of parsecs. Generally, the spatial resolution is defined by (some multiple of) the softening length, whose value is set as a compromise between the desired resolution and the need for limiting small-scale noise. Several works have studied the question whether a softening scale fixed in space and time provides a good enough modelling of an astrophysical system. Here we address this question within the context of dynamical simulations and disc instabilities. We first follow the evolution of a galaxy-like object in isolation and then set up a simulation of an idealised merger event. Alongside a run using the standard fixed-softening approach, we performed simulations where the softening lengths were let to vary from particle to particle according to the evolution of the local density field in space and time. Even though during the most violent phases of the merging the fixed-softening simulation tends to underestimate the resulting matter densities, as far as the evolution of the disc component is concerned we found no significant differences among the runs. We conclude that using an appropriate fixed softening scale is a safe approach to the problem of modelling an N-body, non-cosmological disc galaxy system.

Stellar bars can lead to gas inflow toward the center of a galaxy and stimulate nuclear star formation. However, there is no compelling evidence on whether they also feed a central supermassive black hole: by measuring the fractions of barred active and inactive galaxies, previous studies have yielded conflicting results. In this paper, we aim to understand the lack of observational evidence for bar-driven active galactic nucleus (AGN) activity by studying a sample of 41 nearby (d < 35 Mpc) barred galaxies from the Spitzer Survey for Stellar Structure in Galaxies. We use Chandra observations to measure nuclear 2--10 keV X-ray luminosities and estimate Eddington ratios, together with Spitzer 3.6um imaging to quantify the strength of the stellar bar in two independent ways: (1) from its structure, as traced by its ellipticity and boxiness, and (2) from its gravitational torque Q_b, taken as the maximum ratio of the tangential force to the mean background radial force. In this way, rather than discretizing the presence of both stellar bars and nuclear activity, we are able to account for the continuum of bar strengths and degrees of AGN activity. We find nuclear X-ray sources in 31 out of 41 galaxies with median X-ray luminosity and Eddington ratio of L_X=4.3x10^{38} erg/s and L_bol/L_Edd=6.9x10^{-6} respectively, consistent with low-luminosity AGN activity. Including upper limits for those galaxies without nuclear detections, we find no significant correlation between any of the bar strength indicators and the degree of nuclear activity, irrespective of galaxy luminosity, stellar mass, Hubble type, or bulge size. Strong bars do not favor brighter or more efficient nuclear activity, implying that at least for the low-luminosity regime, supermassive black hole fueling is not closely connected to large scale features.

We study the m=1 distortions (lopsidedness) in the stellar components of 167 nearby galaxies that span a wide range of morphologies and luminosities. We confirm the previous findings of 1) a high incidence of lopsidedness in the stellar distributions, 2) increasing lopsidedness as a function of radius out to at least 3.5 exponential scale lengths, and 3) greater lopsidedness, over these radii, for galaxies of later type and lower surface brightness. Additionally, the magnitude of the lopsidedness 1) correlates with the character of the spiral arms (stronger arm patterns occur in galaxies with less lopsidedness), 2) is not correlated with the presence or absence of a bar, or the strength of the bar when one is present, 3) is inversely correlated to the stellar mass fraction, f_*, within one radial scale length, and 4) correlates directly with f_* measured within the radial range over which we measure lopsidedness. We interpret these findings to mean that lopsidedness is a generic feature of galaxies and does not, generally, depend on a rare event, such as a direct accretion of a satellite galaxy onto the disk of the parent galaxy. While lopsidedness may be caused by several phenomena, moderate lopsidedness (<A_1>_i + <A_1>_o)/2 < 0.3) is likely to reflect halo asymmetries to which the disk responds or a gravitationally self-generated mode . We hypothesize that the magnitude of the stellar response depends both on how centrally concentrated the stars are with respect to the dark matter and whether there are enough stars in the region of the lopsidedness that self-gravity is dynamically important.

We have analyzed the radial distribution of old stars in a sample of 218 nearby face-on disks, using deep 3.6um images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). In particular, we have studied the structural properties of those disks with a broken or down-bending profile. We find that, on average, disks with a genuine single exponential profile have a scale-length and a central surface brightness which are intermediate to those of the inner and outer components of a down-bending disk with the same total stellar mass. In the case of barred galaxies, the ratio between the break and the bar radii (Rbr/Rbar) depends strongly on the total stellar mass of the galaxy. For galaxies more massive than 10^10 Msun, the distribution is bimodal, peaking at Rbr/Rbar~2 and ~3.5. The first peak, which is the most populated one, is linked to the Outer Lindblad Resonance of the bar, whereas the second one is consistent with a dynamical coupling between the bar and the spiral pattern. For galaxies below 10^10 Msun, breaks are found up to ~10 Rbar, but we show that they could still be caused by resonances given the rising nature of rotation curves in these low-mass disks. While not ruling out star formation thresholds, our results imply that radial stellar migration induced by non-axysymmetric features can be responsible not only for those breaks at 2 Rbar, but also for many of those found at larger radii.

We present the kinematic results from our ARGOS spectroscopic survey of the Galactic bulge of the Milky Way. Our aim is to understand the formation of the Galactic bulge. We examine the kinematics of about 17,400 stars in the bulge located within 3.5 kpc of the Galactic centre, identified from the 28,000 star ARGOS survey. We aim to determine if the formation of the bulge has been internally driven from disk instabilities as suggested by its boxy shape, or if mergers have played a significant role as expected from Lambda CDM simulations. From our velocity measurements across latitudes b = -5 deg, -7.5 deg and -10 deg we find the bulge to be a cylindrically rotating system that transitions smoothly out into the disk. Within the bulge, we find a kinematically distinct metal-poor population ([Fe/H] < -1.0) that is not rotating cylindrically. The 5% of our stars with [Fe/H] < -1.0 are a slowly rotating spheroidal population, which we believe are stars of the metal weak thick disk and halo which presently lie in the inner Galaxy. The kinematics of the two bulge components that we identified in ARGOS paper III (mean [Fe/H] = -0.25 and [Fe/H] = +0.15, respectively) demonstrate that they are likely to share a common formation origin and are distinct from the more metal poor populations of the thick disk and halo which are colocated inside the bulge. We do not exclude an underlying merger generated bulge component but our results favour bulge formation from instabilities in the early thin disk.

The fractions and dimension of bars, rings and lenses are studied in the Near-IR S0 galaxy Survey (NIRS0S). We find evidence that multiple lenses in some barred S0s are related to bar resonances in a similar manner as the inner and outer rings, for which the outer/inner length ratio 2. Inner lenses in the non-barred galaxies normalized to galaxy diameter are clearly smaller than those in the barred systems. Interestingly, these small lenses in the non-barred galaxies have similar sizes as barlenses (lens-like structures embedded in a bar), and therefore might actually be barlenses in former barred galaxies, in which the outer, more elongated bar component, has been destroyed. We also find that fully developed inner lenses are on average a factor 1.3 larger than bars, whereas inner rings have similar sizes as bars. The fraction of inner lenses is found to be constant in all family classes (A, AB, B). Nuclear bars appear most frequently among the weakly barred (AB) galaxies, which is consistent with the theoretical models by Maciejewski & Athanassoula (2008). Similar sized bars as the nuclear bars were detected in seven 'non-barred' S0s. Galaxy luminosity does not uniquely define the sizes of bars or bar-related structures, neither is there any upper limit in galaxy luminosity for bar formation. Although all the family classes cover the same range of galaxy luminosity, the non-barred (A) galaxies are on average 0.6 mag brighter than the strongly barred (B) systems. Overall, our results are consistent with the idea that bars play an important role in the formation of the structure components of galaxies. The fact that multiple lenses are common in S0s, and that at least the inner lenses can have very old stellar populations, implies that the last destructive merger, or major gas accretion event, must have taken place at a fairly high redshift.

We present a study of the kinematics of the isolated spiral galaxy NGC 864, using H{\alpha} Fabry-Perot data obtained with the GH{\alpha}FaS instrument at the William Herschel Telescope in La Palma, complemented with images at 3.6 {\mu}m, in the R band and in H{\alpha} filter, and integral field spectroscopic data. The resulting data cubes and velocity maps allow the study of the kinematics of the galaxy, including in-depth investigations of the rotation curve, velocity moment maps, velocity residual maps, gradient maps and position-velocity diagrams. We find asymmetries in the velocity field in the bar zone, caused by non-circular motions, probably in response to the potential of the bar. We also find a flat-profile bar, in agreement with the strong bar, with the grand design spiral pattern, and with the gap between the ends of the bar and the start of the spiral arms. We quantify the rate of massive star formation, which is concentrated in the two spiral arms.

We present the metallicity results from the ARGOS spectroscopic survey of the Galactic bulge. Our aim is to understand the formation of the Galactic bulge: did it form via mergers, as expected from Lambda CDM theory, or from disk instabilities, as suggested by its boxy/peanut shape, or both? We have obtained spectra for 28,000 stars at a spectral resolution of R = 11,000. From these spectra, we have determined stellar parameters and distances to an accuracy of < 1.5 kpc. The stars in the inner Galaxy span a large range in [Fe/H], -2.8 < [Fe/H] < +0.6. From the spatial distribution of the red clump stars as a function of [Fe/H] (Ness et al. 2012a), we propose that the stars with [Fe/H] > -0.5 are part of the boxy/peanut bar/bulge. We associate the lower metallicity stars ([Fe/H] < -0.5) with the thick disk, which may be puffed up in the inner region, and with the inner regions of the metal-weak thick disk and inner halo. For the bulge stars with [Fe/H] > -0.5, we find two discrete populations; (i) stars with [Fe/H] ~ -0.25 which provide a roughly constant fraction of the stars in the latitude interval b = -5 deg to -10 deg, and (ii) a kinematically colder, more metal-rich population with mean [Fe/H] ~ +0.15 which is more prominent closer to the plane. The changing ratio of these components with latitude appears as a vertical abundance gradient of the bulge. We attribute both of these bulge components to instability-driven bar/bulge formation from the thin disk. We do not exclude a weak underlying classical merger-generated bulge component, but see no obvious kinematic association of any of our bulge stars with such a classical bulge component. [abridged]

We describe the motivation, field locations and stellar selection for the ARGOS spectroscopic survey of 28,000 stars in the bulge and inner disk of the Milky Way galaxy across latitudes of b = -5 deg to -10 deg. The primary goal of this survey is to constrain the formation processes of the bulge and establish whether it is predominantly a merger or instability remnant. From the spectra (R = 11,000), we have measured radial velocities and determined stellar parameters, including metallicities and [alpha/Fe] ratios. Distances were estimated from the derived stellar parameters and about 14,000 stars are red giants within 3.5 kpc of the Galactic centre. In this paper, we present the observations and analysis methods. Subsequent papers (III and IV) will discuss the stellar metallicity distribution and kinematics of the Galactic bulge and inner disk, and the implications for the formation of the bulge.

Bars play a major role in driving the evolution of disk galaxies and in shaping their present properties. They cause angular momentum to be redistributed within the galaxy, emitted mainly from (near-)resonant material at the inner Lindblad resonance of the bar, and absorbed mainly by (near-)resonant material in the spheroid (i.e., the halo and, whenever relevant, the bulge) and in the outer disk. Spheroids delay and slow down the initial growth of the bar they host, but, at the later stages of the evolution, they strengthen the bar by absorbing angular momentum. Increased velocity dispersion in the (near-)resonant regions delays bar formation and leads to less strong bars. When bars form they are vertically thin, but soon their inner parts puff up and form what is commonly known as the boxy/peanut bulge. This gives a complex and interesting shape to the bar which explains a number of observations and also argues that the COBE/DIRBE bar and the Long bar in our Galaxy are, respectively, the thin and the thick part of a single bar. The value of the bar pattern speed may be set by optimising the balance between emitters and absorbers, so that a maximum amount of angular momentum is redistributed. As they evolve, bars grow stronger and rotate slower. Bars also redistribute matter within the galaxy, create a disky bulge (pseudo-bulge), increase the disk scale-length and extent and drive substructures such as spirals and rings. They also affect the shape of the inner part of the spheroid, which can evolve from spherical to triaxial.

We follow the formation and evolution of bars in N-body simulations of disc galaxies with gas and/or a triaxial halo. We find that both the relative gas fraction and the halo shape play a major role in the formation and evolution of the bar. In gas-rich simulations, the disc stays near-axisymmetric much longer than in gas-poor ones, and, when the bar starts growing, it does so at a much slower rate. Due to these two effects combined, large-scale bars form much later in gas-rich than in gas-poor discs. This can explain the observation that bars are in place earlier in massive red disc galaxies than in blue spirals. We also find that the morphological characteristics in the bar region are strongly influenced by the gas fraction. In particular, the bar at the end of the simulation is much weaker in gas-rich cases. In no case did we witness bar destruction. Halo triaxiality has a dual influence on bar strength. In the very early stages of the simulation it induces bar formation to start earlier. On the other hand, during the later, secular evolution phase, triaxial haloes lead to considerably less increase of the bar strength than spherical ones. The shape of the halo evolves considerably with time. The inner halo parts may become more elongated, or more spherical, depending on the bar strength. The main body of initially triaxial haloes evolves towards sphericity, but in initially strongly triaxial cases it stops well short of becoming spherical. Part of the angular momentum absorbed by the halo generates considerable rotation of the halo particles that stay located relatively near the disc for long periods of time. Another part generates halo bulk rotation, which, contrary to that of the bar, increases with time but stays small.

The mechanism causing breaks in the radial surface brightness distribution of spiral galaxies is not yet well known. Despite theoretical efforts, there is not a unique explanation for these features and the observational results are not conclusive. In an attempt to address this problem, we have selected a sample of 34 highly inclined spiral galaxies present both in the Sloan Digital Sky Survey and in the Spitzer Survey of Stellar Structure in Galaxies. We have measured the surface brightness profiles in the five Sloan optical bands and in the 3.6$\mu m$ Spitzer band. We have also calculated the color and stellar surface mass density profiles using the available photometric information, finding two differentiated features: an innermost break radius at distances of $\sim 8 \pm 1$ kpc [$0.77 \pm 0.06$ $R_{25}$] and a second characteristic radius, or truncation radius, close to the outermost optical extent ($\sim 14 \pm 2$ kpc [$1.09 \pm 0.05$ $R_{25}$]) of the galaxy. We propose in this work that the breaks might be a phenomena related to a threshold in the star formation, while truncations are more likely a real drop in the stellar mass density of the disk associated with the maximum angular momentum of the stars.

Breaks in the radial luminosity profiles of galaxies have been until now mostly studied averaged over discs. Here we study separately breaks in thin and thick discs in 70 edge-on galaxies using imaging from the Spitzer Survey of Stellar Structure in Galaxies. We built luminosity profiles of the thin and the thick discs parallel to midplanes and we found that thin discs often truncate (77%). Thick discs truncate less often (31%), but when they do, their break radius is comparable with that in the thin disc. This suggests either two different truncation mechanisms - one of dynamical origin affecting both discs simultaneously and another one only affecting the thin disc - or a single mechanism that creates a truncation in one disc or in both depending on some galaxy property. Thin discs apparently antitruncate in around 40% of galaxies. However, in many cases, these antitruncations are an artifact caused by the superposition of a thin disc and a thick disc with the latter having a longer scale length. We estimate the real thin disc antitruncation fraction to be less than 15%. We found that the ratio of the thick and thin stellar disc mass is roughly constant (0.2<M_T/M_t<0.7) for circular velocities v_c>120 km/s, but becomes much larger at smaller velocities. We hypothesize that this is due to a combination of a high efficiency of supernova feedback and a slower dynamical evolution in lower-mass galaxies causing stellar thin discs to be younger and less massive than in higher-mass galaxies.

We present the results from the Giant Metrewave Radio Telescope (GMRT) interferometric HI and 20 cm radio continuum observations of CIG 85, an isolated asymmetric galaxy from the AMIGA (Analysis of the Interstellar Medium of Isolated GAlaxies) sample. Despite being an isolated galaxy, CIG 85 showed an appreciable optical and HI spectral asymmetry and therefore was an excellent candidate for resolved HI studies to understand the reasons giving rise to asymmetries in isolated galaxies. The galaxy was imaged in HI and 20 cm radio continuum using the GMRT. For a detailed discussion of the results we also made use of multi-wavelength data from archival SDSS, GALEX and Halpha imaging. We find the HI in CIG 85 to have a clumpy, asymmetric distribution which in the NW part is correlated with optical tail like features, but the HI velocity field displays a relatively regular rotation pattern. Evaluating all the observational evidence, we come to a conclusion that CIG 85 is most likely a case of a disturbed spiral galaxy which now appears to have the morphology of an irregular galaxy. Although it is currently isolated from major companions, the observational evidence is consistent with HI asymmetries, a highly disturbed optical disk and recent increase in star formation having been caused by a minor merger, remnants of which are now projected in front of the optical disk. If this is correct, the companion will be fully accreted by CIG 85 in the near future.

We present observational evidence for the inhibition of bar formation in dispersion-dominated (dynamically hot) galaxies by studying the relationship between galactic structure and host galaxy kinematics in a sample of 257 galaxies between 0.1 $<$ z $\leq$ 0.84 from the All-Wavelength Extended Groth Strip International Survey (AEGIS) and the Deep Extragalactic Evolutionary Probe 2 (DEEP2) survey. We find that bars are preferentially found in galaxies that are massive and dynamically cold (rotation-dominated) and on the stellar Tully-Fisher relationship, as is the case for barred spirals in the local Universe. The data provide at least one explanation for the steep ($\times$3) decline in the overall bar fraction from z=0 to z=0.84 in L$^*$ and brighter disks seen in previous studies. The decline in the bar fraction at high redshift is almost exclusively in the lower mass (10 $<$ log M$_{*}$(\Msun)$<$ 11), later-type and bluer galaxies. A proposed explanation for this "downsizing" of the bar formation / stellar structure formation is that the lower mass galaxies may not form bars because they could be dynamically hotter than more massive systems from the increased turbulence of accreting gas, elevated star formation, and/or increased interaction/merger rate at higher redshifts. The evidence presented here provides observational support for this hypothesis. However, the data also show that not every disk galaxy that is massive and cold has a stellar bar, suggesting that mass and dynamic coldness of a disk are necessary but not sufficient conditions for bar formation -- a secondary process, perhaps the interaction history between the dark matter halo and the baryonic matter, may play an important role in bar formation.

The baryonic discs of galaxies are believed to alter the shapes of the dark matter haloes in which they reside. We perform a set of hydrodynamical N-body simulations of disc galaxies with triaxial dark matter haloes, using elliptical discs with a gaseous component as initial conditions. We explore models of different halo triaxiality and also of different initial gas fractions, which allows us to evaluate how each affects the formations of the bar. Due to star formation, models of all halo shapes and of all initial gas fractions reach approximately the same gas content at the end of the simulation. Nevertheless, we find that the presence of gas in the early phases has important effects on the subsequent evolution. Bars are generally weaker for larger initial gas content and for larger halo triaxiality. The presence of gas, however, is a more efficient factor in inhibiting the formation of a strong bar than halo triaxiality is.

We discuss the properties of spiral arms in a N-body simulation of a barred galaxy and present evidence that these are manifold-driven. The strongest evidence comes from following the trajectories of individual particles. Indeed, these move along the arms while spreading out a little. In the neighbourhood of the Lagrangian points they follow a variety of paths, as expected by manifold-driven trajectories. Further evidence comes from the properties of the arms themselves, such as their shape and growth pattern. The shape of the manifold arms changes considerably with time, as expected from the changes in the bar strength and pattern speed. In particular, the radial extent of the arms increases with time, thus bringing about a considerable increase of the disc size, by as much as ~50% in about a Gyr.

Near the minor axis of the Galactic bulge, at latitudes b < -5 degrees, the red giant clump stars are split into two components along the line of sight. We investigate this split using the three fields from the ARGOS survey that lie on the minor axis at (l,b) = (0,-5), (0,-7.5), (0,-10) degrees. The separation is evident for stars with [Fe/H] > -0.5 in the two higher-latitude fields, but not in the field at b = -5 degrees. Stars with [Fe/H] < -0.5 do not show the split. We compare the spatial distribution and kinematics of the clump stars with predictions from an evolutionary N-body model of a bulge that grew from a disk via bar-related instabilities. The density distribution of the peanut-shaped model is depressed near its minor axis. This produces a bimodal distribution of stars along the line of sight through the bulge near its minor axis, very much as seen in our observations. The observed and modelled kinematics of the two groups of stars are also similar. We conclude that the split red clump of the bulge is probably a generic feature of boxy/peanut bulges that grew from disks, and that the disk from which the bulge grew had relatively few stars with [Fe/H] < -0.5

We study the properties of two bars formed in fully cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies. In one case, the bar formed in a system with disc, bulge and halo components and is relatively strong and long, as could be expected for a system where the spheroid strongly influences the evolution. The second bar is less strong, shorter, and formed in a galaxy with no significant bulge component. We study the strength and length of the bars, the stellar density profiles along and across the bars and the velocity fields in the bar region. We compare them with the results of dynamical (idealised) simulations and with observations, and find, in general, a good agreement, although we detect some important differences as well. Our results show that more or less realistic bars can form naturally in a $\Lambda$CDM cosmology, and open up the possibility to study the bar formation process in a more consistent way than previously done, since the host galaxies grow, accrete matter and significantly evolve during the formation and evolution of the bar.

About two-thirds of present-day, large galaxies are spirals such as the Milky Way or Andromeda, but the way their thin rotating disks formed remains uncertain. Observations have revealed that half of their progenitors, six billion years ago, had peculiar morphologies and/or kinematics, which exclude them from the Hubble sequence. Major mergers, i.e., fusions between galaxies of similar mass, are found to be the likeliest driver for such strong peculiarities. However, thin disks are fragile and easily destroyed by such violent collisions, which creates a critical tension between the observed fraction of thin disks and their survival within the L-CDM paradigm. Here we show that the observed high occurrence of mergers amongst their progenitors is only apparent and is resolved when using morpho-kinematic observations which are sensitive to all the phases of the merging process. This provides an original way of narrowing down observational estimates of the galaxy merger rate and leads to a perfect match with predictions by state-of-the-art L-CDM semi-empirical models with no particular fine-tuning needed. These results imply that half of local thin disks do not survive but are actually rebuilt after a gas-rich major merger occurring in the past nine billion years, i.e., two-thirds of the lifetime of the Universe. This emphasizes the need to study how thin disks can form in halos with a more active merger history than previously considered, and to investigate what is the origin of the gas reservoir from which local disks would reform.

Tidal debris around galaxies can yield important clues on their evolution. We have identified tidal debris in 11 early type galaxies (T \leq 0) from a sample of 65 early types drawn from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The tidal debris includes features such as shells, ripples and tidal tails. A variety of techniques, including two-dimensional decomposition of galactic structures, was used to quantify the residual tidal features. The tidal debris contributes ~3 - 10% to the total 3.6 {\mu}m luminosity of the host galaxy. Structural parameters of the galaxies were estimated using two-dimensional profile fitting. We investigate the locations of galaxies with tidal debris in the Fundamental Plane and Kormendy relation. We find that galaxies with tidal debris lie within the scatter of early type galaxies without tidal features. Assuming that the tidal debris is indicative of recent gravitational interaction or merger, this suggests that these galaxies have either undergone minor merging events so that the overall structural properties of the galaxies are not significantly altered, or they have undergone a major merging events but already have experienced sufficient relaxation and phase-mixing so that their structural properties become similar to those of the non-interacting early type galaxies.

We examine the effect of circumstellar dust extinction on the near-IR contribution of asymptotic giant branch (AGB) stars in intermediate-age clusters throughout the disk of M100. For our sample of 17 AGB-dominated clusters we extract optical-to-mid-IR SEDs and find that NIR brightness is coupled to the mid-IR dust emission in such a way that a significant reduction of AGB light, of up to 1 mag in K-band, follows from extinction by the dust shell formed during this stage. Since the dust optical depth varies with AGB chemistry (C-rich or O-rich), our results suggest that the contribution of AGB stars to the flux from their host clusters will be closely linked to the metallicity and the progenitor mass of the AGB star, to which dust chemistry and mass-loss rate are sensitive. Our sample of clusters--each the analogue of a ~1 Gyr old post-starburst galaxy--has implications within the context of mass and age estimation via SED modelling at high z: we find that the average ~0.5 mag extinction estimated here may be sufficient to reduce the AGB contribution in (rest-frame) K-band from ~70%, as predicted in the latest generation of synthesis models, to ~35%. Our technique for selecting AGB-dominated clusters in nearby galaxies promises to be effective for discriminating the uncertainties associated with AGB stars in intermediate-age populations that plague age and mass estimation in high-z galaxies.

Many haloes of nearby disc galaxies contain faint and extended features, including loops, which are often interpreted as relics of satellite infall in the main galaxy's potential well. In most cases, however, the residual nucleus of the satellite is not seen, although it is predicted by numerical simulations. We test whether such faint and extended features can be associated to gas-rich, major mergers, which may also lead to disc rebuilding and thus be a corner stone for the formation of spiral galaxies. Using the TreeSPH code GADGET-2, we model the formation of an almost bulge-less galaxy similar to NGC 5907 (B/T $\le$ 0.2) after a gas-rich major merger. We indeed find that 3:1 major mergers can form features similar to the loops found in many galactic haloes, including in NGC 5907, and can reproduce an extended thin disc, a bulge, as well as the pronounced warp of the gaseous disc. Even though it remains difficult to fully cover the large volume of free parameters, the present modelling of the loops in NGC 5907 proves that they could well be the result of a major merger. It has many advantages over the satellite infall scenario; e.g., it solves the problem of the visibility of the satellite remnant, and it may explain some additional features in the NGC 5907 halo, as well as some gas properties of this system. For orbital parameters derived from cosmological simulations, the loops in NGC 5907 can be reproduced by major mergers (3:1 to 5:1) and possibly by intermediate mergers (5:1 to 12:1). The major merger scenario thus challenges the minor merger one and could explain many properties that haloes of spiral galaxies have in common, including their red colours and the presence of faint extended features.

I review some of the work on bars which is closely linked to the bar/bulge system in our Galaxy. Several independent studies, using totally independent methods, come to the same results about the 3D structure of a bar, i.e., that a bar is composed of a vertically thick inner part and a vertically thin outer part. I give examples of this from simulations and substantiate the discussion with input from orbital structure analysis and from observations. The thick part has a considerably shorter radial extent than the thin part. I then see how this applies to our Galaxy, where two bars have been reported, the COBE/DIRBE bar and the Long bar. Comparing their extents and making the reasonable and necessary assumption that our Galaxy has properties similar to those of other galaxies of similar type, leads to the conclusion that these two bars can not form a standard double bar system. I then discuss arguments in favour of the two bars being simply different parts of the same bar, the COBE/DIRBE bar being the thick inner part and the Long bar being the thin outer part of this bar. I also very briefly discuss some related new results. I first consider bar formation and evolution in disc galaxies with a gaseous component - including star formation, feedback and evolution - and a triaxial halo. Then I consider bar formation in a fully cosmological context using hydrodynamical LCDM simulations, where the host galaxies grow, accrete matter and significantly evolve during the formation and evolution of the bar.

We use both N-body simulations and integration in fixed potentials to explore the stability and the long-term secular evolution of self-consistent, equilibrium, non-rotating, triaxial spheroidal galactic models. More specifically, we consider Dehnen models built with the Schwarzschild method. We show that short-term stability depends on the degree of velocity anisotropy (radially anisotropic models are subject to rapid development of radial-orbit instability). Long-term stability, on the other hand, depends mainly on the properties of the potential, and in particular, on whether it admits a substantial fraction of strongly chaotic orbits. We show that in the case of a weak density cusp (gamma=1 Dehnen model) the N-body model is remarkably stable, while the strong-cusp (gamma=2) model exhibits substantial evolution of shape away from triaxiality, which we attribute to the effect of chaotic diffusion of orbits. The different behaviour of these two cases originates from the different phase space structure of the potential; in the weak-cusp case there exist numerous resonant orbit families that impede chaotic diffusion. We also find that it is hardly possible to affect the rate of this evolution by altering the fraction of chaotic orbits in the Schwarzschild model, which is explained by the fact that the chaotic properties of an orbit are not preserved by the N-body evolution. There are, however, parameters in Schwarzschild modelling that do affect the stability of an N-body model, so we discuss the recipes how to build a `good' Schwarzschild model.

The dynamical evolution of stellar clusters is driven to a large extent by their environment. Several studies so far have considered the effect of tidal fields and their variations, such as, e.g., from giant molecular clouds, galactic discs, or spiral arms. In this paper we will concentrate on a tidal field whose effects on star clusters have not yet been studied, namely that of bars. We present a set of direct N-body simulations of star clusters moving in an analytic potential representing a barred galaxy. We compare the evolution of the clusters moving both on different planar periodic orbits in the barred potential and on circular orbits in a potential obtained by axisymmetrising its mass distribution. We show that both the shape of the underlying orbit and its stability have strong impact on the cluster evolution as well as the morphology and orientation of the tidal tails and the sub-structures therein. We find that the dissolution time-scale of the cluster in our simulations is mainly determined by the tidal forcing along the orbit and, for a given tidal forcing, only very little by the exact shape of the gravitational potential in which the cluster is moving.

With the aim of constructing accurate 2D maps of the stellar mass distribution in nearby galaxies from S4G 3.6 and 4.5 micron images, we report on the separation of the light from old stars from the emission contributed by contaminants (e.g. hot dust and the 3.3 micron PAH feature). Results for a small sample of six disk galaxies (NGC 1566, NGC 2976, NGC 3031, NGC 3184, NGC 4321, and NGC 5194) with a range of morphological properties, dust contents and star formation histories are presented to demonstrate our approach. We use an Independent Component Analysis (ICA) technique designed to separate statistically independent source distributions, maximizing the distinction in the [3.6]-[4.5] colors of the sources. The technique also removes emission from intermediate-age evolved red objects with a low mass-to-light ratio, such as asymptotic giant branch (AGB) and red supergiant (RSG) stars, revealing maps of the underlying old distribution of light with [3.6]-[4.5] colors consistent with the colors of K and M giants. Contaminants are identified via comparison to the non-stellar emission imaged at 8 microns, which is dominated by the broad PAH feature. Using the measured 3.6/8 micron ratio to select the individual contaminants, we find that hot dust and PAH together contribute between ~5-15% to the integrated light at 3.6 microns, while light from regions dominated by intermediate-age stars accounts for only 1-5%. Locally, however, the contribution from either contaminant can reach much higher levels; dust contributes on average 22% to the emission in star-forming regions throughout the sample, while intermediate age-stars contribute upwards of 50% in localized knots. The removal of these contaminants with ICA leaves maps of the old stellar disk that retain a high degree of structural information and are ideally suited for tracing the stellar mass, as will be the focus in a companion paper.

Using the deepest and the most complete set of observations of distant galaxies, we investigate how extended disks could have formed. Observations include spatially-resolved kinematics, detailed morphologies and photometry from UV to mid-IR. Six billion years ago, half of the present-day spiral progenitors had anomalous kinematics and morphologies, as well as relatively high gas fractions. We argue that gas-rich major mergers, i.e., fusions between gas-rich disk galaxies of similar mass, can be the likeliest driver for such strong peculiarities. This suggests a new channel of disk formation, e.g. many disks could be reformed after gas-rich mergers. This is found to be in perfect agreement with predictions from the state-of-the-art LCDM semi-empirical models: due to our sensitivity in detecting mergers at all phases, from pairs to relaxed post-mergers, we find a more accurate merger rate. The scenario can be finally confronted to properties of nearby galaxies, including M31 and galaxies showing ultra-faint, gigantic structures in their haloes.

Most, if not all, disk galaxies have a thin (classical) disk and a thick disk. In most models thick disks are thought to be a necessary consequence of the disk formation and/or evolution of the galaxy. We present the results of a study of the thick disk properties in a sample of carefully selected edge-on galaxies with types ranging from T=3 to T=8. We fitted one-dimensional luminosity profiles with physically motivated functions - the solutions of two stellar and one gaseous isothermal coupled disks in equilibrium - which are likely to yield more accurate results than other functions used in previous studies. The images used for the fits come from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We found that thick disks are on average more massive than previously reported, mostly due to the selected fitting function. Typically, the thin and the thick disk have similar masses. We also found that thick disks do not flare significantly within the observed range in galactocentric radii and that the ratio of thick to thin disk scaleheights is higher for galaxies of earlier types. Our results tend to favor an in situ origin for most of the stars in the thick disk. In addition the thick disk may contain a significant amount of stars coming from satellites accreted after the initial build-up of the galaxy and an extra fraction of stars coming from the secular heating of the thin disk by its own overdensities. Assigning thick disk light to the thin disk component may lead to an underestimate of the overall stellar mass in galaxies, because of different mass to light ratios in the two disk components. On the basis of our new results, we estimate that disk stellar masses are between 10% and 50% higher than previously thought and we suggest that thick disks are a reservoir of "local missing baryons".

NGC 4013 is a nearby Sb edge-on galaxy known for its "prodigious" HI warp and its "giant" tidal stream. Previous work on this unusual object shows that it cannot be fitted satisfactorily by a canonical thin+thick disk structure. We have produced a new decomposition of NGC 4013, considering three stellar flattened components (thin+thick disk plus an extra and more extended component) and one gaseous disk. All four components are considered to be gravitationally coupled and isothermal. To do so, we have used the 3.6micron images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find evidence for NGC4013 indeed having a thin and a thick disk and an extra flattened component. This smooth and extended component (scaleheight z_EC kpc) could be interpreted as a thick disk or as a squashed ellipsoidal halo and contains ~20% of the total mass of all three stellar components. We argue it is unlikely to be related to the ongoing merger or due to the off-plane stars from a warp in the other two disk components. Instead, we favor a scenario in which the thick disk and the extended component were formed in a two-stage process, in which an initially thick disk has been dynamically heated by a merger soon enough in the galaxy history to have a new thick disk formed within it.

We study the bar-driven dynamics in the inner part of the Milky Way by using invariant manifolds. This theory has been successfully applied to describe the morphology and kinematics of rings and spirals in external galaxies, and now, for the first time, we apply it to the Milky Way. We compute the orbits confined by the invariant manifolds of the unstable periodic orbits located at the ends of the bar. We discuss whether the COBE/DIRBE bar and the Long bar compose a single bar or two independent bars and perform a number of comparisons which, taken together, argue strongly in favour of the former. More specifically, we favour the possibility that the so-called COBE/DIRBE bar is the boxy/peanut bulge of a bar whose outer thin parts are the so-called Long bar. This possibility is in good agreement both with observations of external galaxies, with orbital structure theory and with simulations. We then analyse in detail the morphology and kinematics given by five representative Galactic potentials. Two have a Ferrers bar, two have a quadrupole bar and the last one a composite bar. We first consider only the COBE/DIRBE bar and then extend it to include the effect of the Long bar. We find that the large-scale structure given by the manifolds describes an inner ring, whose size is similar to the near and far 3-kpc arm, and an outer ring, whose properties resemble those of the Galactic Molecular Ring. We also analyse the kinematics of these two structures, under the different galactic potentials, and find they reproduce the main over-densities found in the galactic longitude-velocity CO diagram. Finally, we consider for what model parameters, the global morphology of the manifolds may reproduce the two outer spiral arms. We conclude that this would necessitate either more massive and more rapidly rotating bars, or including in the potential an extra component describing the spiral arms.

We study the Kennicutt-Schmidt star formation law and efficiency in the gaseous disk of the isolated galaxy CIG 96 (NGC 864), with special emphasis on its unusually large atomic gas (HI) disk (r_HI/r_25 = 3.5, r_25 = 1.'85). We present deep GALEX near and far ultraviolet observations, used as a recent star formation tracer, and we compare them with new, high resolution (16", or 1.6 kpc) VLA HI observations. The UV and HI maps show good spatial correlation outside the inner 1', where the HI phase dominates over H_2. Star-forming regions in the extended gaseous disk are mainly located along the enhanced HI emission within two (relatively) symmetric giant gaseous spiral arm-like features, which emulate a HI pseudo-ring at a r \simeq 3' . Inside such structure, two smaller gaseous spiral arms extend from the NE and SW of the optical disk and connect to the previously mentioned HI pseudo-ring. Interestingly, we find that the (atomic) Kennicutt-Schmidt power law index systematically decreases with radius, from N \simeq 3.0 +- 0.3 in the inner disk (0.'8 - 1.'7) to N = 1.6 +- 0.5 in the outskirts of the gaseous disk (3.'3 - 4.'2). Although the star formation efficiency (SFE), the star formation rate per unit of gas, decreases with radius where the HI component dominates as is common in galaxies, we find that there is a break of the correlation at r = 1.5 r_25. At radii 1.5 r_25 < r < 3.5 r_25, mostly within the HI pseudo-ring structure, there exist regions whose SFE remains nearly constant, SFE \simeq 10^-11 yr^-1. We discuss about possible mechanisms that might be triggering the star formation in the outskirts of this galaxy, and we suggest that the constant SFE for such large radii r > 2 r_25 and at such low surface densities might be a common characteristic in extended UV disk galaxies.

We present Expanded Very Large Array (EVLA) D-array observations in the 21-cm line of neutral hydrogen (HI) of CIG 292, an isolated SA(s)b galaxy at a distance of ~ 24.3Mpc. From previous HI single dish observations the galaxy was known to have a mildly asymmetric HI profile (A_flux =1.23 {\pm}0.3). Our EVLA observations show there is ~ 12% more HI projected South of the optical centre (approaching velocities) than in the North (receding velocities), despite the HI extending ~ 16% further to the North than the South. The HI projected within the optical disk must have been perturbed within the HI relaxation time (~ 10^8 yr) which implies that this can not have been caused by any of the three nearest companions, as their distance (~ 0.5Mpc) is too large. Neither HI-rich companions nor tidal tails were found within our field of view and velocity range covered. Our kinematical data suggest that the inner part harbors an oval distortion whereas the outer regions show signs of a modest warp. The mild asymmetry in the HI global profile thus actually masks stronger asymmetries in the two-dimensional distributions of gas and star forming regions in this galaxy. Since the galaxy is isolated, this must predominantly be due to processes related to its formation and secular evolution.

Aims: We aim to develop an algorithm for constructing equilibrium initial conditions for simulations of disk galaxies with a triaxial halo and/or a gaseous component. This will pave the way for N-body simulations of realistic disk galaxies. Methods: We use the iterative method, which we presented in a previous article. The idea of this method is very simple. It relies on constrained evolution. Results: We develop an algorithm for constructing equilibrium models of disk galaxies including a gaseous disk and a triaxial or axisymmetric halo. We discuss two test models. The first model consists of a spherical halo, a stellar disk, and an isothermal gaseous disk. The second model consists of a triaxial halo, a stellar disk, and a star-forming gaseous disk. We demonstrate that both test models are very close to equilibrium, as we had intended.

The distinction between chaotic and regular behavior of orbits in galactic models is an important issue and can help our understanding of galactic dynamical evolution. In this paper, we deal with this issue by applying the techniques of the Smaller (and Generalized) ALingment Indices, SALI (and GALI), to extensive samples of orbits obtained by integrating numerically the equations of motion in a barred galaxy potential. We estimate first the fraction of chaotic and regular orbits for the two-degree-of-freedom (DOF) case (where the galaxy extends only in the (x,y)-space) and show that it is a non-monotonic function of the energy. For the three DOF extension of this model (in the z-direction), we give similar estimates, both by exploring different sets of initial conditions and by varying the model parameters, like the mass, size and pattern speed of the bar. We find that regular motion is more abundant at small radial distances from the center of the galaxy, where the relative non-axisymmetric forcing is relatively weak, and at small distances from the equatorial plane, where trapping around the stable periodic orbits is important. We also find that the variation of the bar pattern speed, within a realistic range of values, does not affect much the phase space's fraction of regular and chaotic motions. Using different sets of initial conditions, we show that chaotic motion is dominant in galaxy models whose bar component is more massive, while models with a fatter or thicker bar present generally more regular behavior. Finally, we find that the fraction of orbits that are chaotic correlates strongly with the bar strength.

If thick disks are ubiquitous and a natural product of disk galaxy formation and/or evolution processes, all undisturbed galaxies which have evolved during a significant fraction of a Hubble time should have a thick disk. The late-type spiral galaxy NGC 4244 has been reported as the only nearby edge-on galaxy without a confirmed thick disk. Using data from the Spitzer Survey of Stellar Structure in Galaxies (S4G) we have identified signs of two disk components in this galaxy. The asymmetries between the light profiles on both sides of the mid-plane of NGC 4244 can be explained by a combination of the galaxy not being perfectly edge-on and a certain degree of opacity of the thin disk. We argue that the subtlety of the thick disk is a consequence of either a limited secular evolution in NGC 4244, a small fraction of stellar material in the fragments which built the galaxy, or a high amount of gaseous accretion after the formation of the galaxy.

[Abridged] We present an updated and improved M_bh-sigma diagram containing 64 galaxies for which M_bh measurements (not just upper limits) are available. Due to new and increased black hole masses at the high-mass end, and a better representation of barred galaxies at the low-mass end, the "classical" (all morphological type) M_bh-sigma relation for predicting black hole masses is log(M_bh/M_Sun) = (5.13+/-0.34)log[sigma/200] + (8.13+/-0.05), with an rms scatter of 0.43 dex. Modifying the regression analysis to correct for a hitherto over-looked sample bias in which black holes with masses <10^6 M_Sun are not (yet) detectable, the relation steepens further to give log(M_bh/M_Sun) = (5.95+/-0.44)log[sigma/200] + (8.15+/-0.06). We have also updated the "barless" and "elliptical-only" M_bh-sigma relations introduced by Graham and Hu in 2008 due to the offset nature of barred/disc galaxies. These relations have a total scatter as low as 0.34 dex and currently define the upper envelope of points in the M_bh-sigma diagram. These relations also have a slope consistent with the value 5, in agreement with the prediction by Silk & Rees based on feedback from massive black holes in bulges built by monolithic-collapse. Using updated virial products and velocity dispersions from 28 active galactic nuclei, we determine that the optimal scaling factor f - which brings their virial products in line with the 64 directly measured black hole masses - is 2.8^{+0.7}_{-0.5}. This is roughly half the value reported by Onken et al. and Woo et al., and consequently halves the mass estimates of most high-redshift quasars. We have explored the results after separating the samples into barred and non-barred galaxies, and we have also developed a preliminary corrective term to the velocity dispersion based on bar dynamics.

We analyze our accurate kinematical data for the old clusters in the inner regions of M31. These velocities are based on high S/N Hectospec data (Caldwell et al 2010). The data are well suited for analysis of M31's inner regions because we took particular care to correct for contamination by unresolved field stars from the disk and bulge in the fibers. The metal poor clusters show kinematics which are compatible with a pressure-supported spheroid. The kinematics of metal-rich clusters, however, argue for a disk population. In particular the innermost region (inside 2 kpc) shows the kinematics of the x2 family of bar periodic orbits, arguing for the existence of an inner Lindblad resonance in M31.

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.

Spitzer Space Telescope Infrared Array Camera (IRAC) imaging provides an opportunity to study all known morphological types of galaxies in the mid-IR at a depth significantly better than ground-based near-infrared and optical images. The goal of this study is to examine the imprint of the de Vaucouleurs classification volume in the 3.6 micron band, which is the best Spitzer waveband for galactic stellar mass morphology owing to its depth and its reddening-free sensitivity mainly to older stars. For this purpose, we have prepared classification images for 207 galaxies from the Spitzer archive, most of which are formally part of the Spitzer Survey of Stellar Structure in Galaxies (S^4G), a Spitzer post-cryogenic ("warm") mission Exploration Science Legacy Program survey of 2,331 galaxies closer than 40 Mpc. For the purposes of morphology, the galaxies are interpreted as if the images are {\it blue light}, the historical waveband for classical galaxy classification studies. We find that 3.6 micron classifications are well-correlated with blue-light classifications, to the point where the essential features of many galaxies look very similar in the two very different wavelength regimes. Drastic differences are found only for the most dusty galaxies. Consistent with a previous study by Eskridge et al. (2002), the main difference between blue light and mid-IR types is an approximately 1 stage interval difference for S0/a to Sbc or Sc galaxies, which tend to appear "earlier" in type at 3.6 microns due to the slightly increased prominence of the bulge, the reduced effects of extinction, and the reduced (but not completely eliminated) effect of the extreme population I stellar component. We present an atlas of all of the 207 galaxies analyzed here, and bring attention to special features or galaxy types that are particularly distinctive in the mid-IR.

We present dynamical models of NGC 4494, which we built using our iterative method presented in a previous paper. These models are live N-body models consisting of equal mass particles, and they are steady state as confirmed by a fully self-consistent evolution. Our goals were twofold. The first one -- namely to test whether our iterative method could indeed be used to construct galactic models following given observational constraints, both photometric and kinematic -- was fully achieved. Our method allowed us to go beyond a simple spherical model and to make full sets of rotating, axisymmetric models without any limitations to the velocity distribution. Our second goal was to understand better the structure of NGC 4494, and more specifically to set constraints on its halo mass. For this we tried three families of models: without halo, with a light halo and with a heavy halo, respectively. Our models reproduce well the photometry and the kinematics, the latter except specific regions where some non-equilibrium or non-axisymmetric structure could be present in the galaxy (e.g. the kinematically decoupled core). However, the lower order moments of the velocity distribution (up to and including the second order) do not allow us to discriminate between the three halos. On the other hand, when we extend the comparison to the higher order moments of the velocity distribution obtained from the long-slit data, we find that our light halo model fits the data better than the no halo, or the heavy halo models. They also reproduce the shape of the angular dependence of the PNe velocity dispersion in the outermost parts of the galaxy, but not the amplitude of its azimuthal variation. This may imply that a yet more general class of models, such as triaxial, may be necessary for a yet better fit.

In a series of papers, we propose a theory to explain the formation and properties of rings and spirals in barred galaxies. The building blocks of these structures are orbits guided by the manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In this paper, the last of the series, we present more comparisons of our theoretical results to observations and also give new predictions for further comparisons. Our theory provides the right building blocks for the rectangular-like bar outline and for ansae. We consider how our results can be used to give estimates for the pattern speed values, as well as their effect on abundance gradients in barred galaxies. We present the kinematics along the manifold loci, to allow comparisons with the observed kinematics along the ring and spiral loci. We consider gaseous arms and their relations to stellar ones. We discuss several theoretical aspects and stress that the orbits that constitute the building blocks of the spirals and rings are chaotic. They are, nevertheless, spatially well confined by the manifolds and are thus able to outline the relevant structures. Such chaos can be termed `confined chaos' and can play a very important role in understanding the formation and evolution of galaxy structures and in galactic dynamics in general. This work, in agreement with several others, argues convincingly that galactic dynamic studies should not be limited to the study of regular motions and orbits.

Cosmological N-body simulations indicate that the dark matter haloes of galaxies should be generally triaxial. Yet, the presence of a baryonic disc is believed to alter the shape of the haloes. Here we aim to study how bar formation is affected by halo triaxiality and how, in turn, the presence of the bar influences the shape of the halo. We perform a set of collisionless N-body simulations of disc galaxies with triaxial dark matter haloes, using elliptical discs as initial conditions. We study models of different halo triaxialities and, to investigate the behaviour of the halo shape in the absence of bar formation, we run models with different disc masses, halo concentrations, disc velocity dispersions and also models where the disc shape is kept artificially axisymmetric. We find that the introduction of a massive disc causes the halo triaxiality to be partially diluted. Once the disc is fully grown, a strong stellar bar develops within the halo that is still non-axisymmetric, causing it to lose its remaining non-axisymmetry. In triaxial haloes in which the initial conditions are such that a bar does not form, the halo is able to remain triaxial and the circularisation of its shape on the plane of the disc is limited to the period of disc growth. We conclude that part of the circularisation of the halo is due to disc growth, but part must be attributed to the formation of a bar. We find that initially circular discs respond excessively to the triaxial potential and become highly elongated. They also lose more angular momentum than the initially elliptical discs and thus form stronger bars. Because of that, the circularisation that their bars induce on their haloes is also more rapid. We also analyse halo vertical shapes and observe that their vertical flattenings remain considerable, meaning that the haloes become approximately oblate by the end of the simulations. [abridged]

We extract at redshift z=0 a Milky Way sized object including gas, stars and dark matter (DM) from a recent, high-resolution cosmological N-body simulation with baryons. Its resolution is sufficient to witness the formation of a rotating disk and bulge at the center of the halo potential. The phase-space structure of the central galactic halo reveals the presence of a dark disk component, that is co-rotating with the stellar disk. At the Earth's location, it contributes to around 25% of the total DM local density, whose value is rho_DM ~ 0.37 GeV/cm^3. The velocity distributions also show strong deviations from pure Gaussian and Maxwellian distributions, with a sharper drop of the high velocity tail. We give a detailed study of the impact of these features on the predictions for DM signals in direct detection experiments. In particular, the question of whether the modulation signal observed by DAMA is or is not excluded by limits set by other experiments (CDMS, XENON and CRESST...) is re-analyzed and compared to the case of a standard Maxwellian halo, in both the elastic and the inelastic scattering scenarios. We find that the compatibility between DAMA and the other experiments is improved. In the elastic scenario, the DAMA modulation signal is slightly enhanced in the so-called channeling region, as a result of several effects. For the inelastic scenario, the improvement of the fit is mainly attributable to the departure from a Maxwellian distribution at high velocity.

Galaxies, in particular disc galaxies, contain a number of structures and substructures with well defined morphological, photometric and kinematic properties. Considerable theoretical effort has been put into explaining their formation and evolution, both analytically and with numerical simulations. In some theories, structures form during the natural evolution of the galaxy, i.e. they are a result of nature. For others, it is the interaction with other galaxies, or with the intergalactic medium -- i.e. nurture -- that accounts for a structure. Either way, the existence and properties of these structures reveal important information on the underlying potential of the galaxy, i.e. on the amount and distribution of matter -- including the dark matter -- in it, and on the evolutionary history of the galaxy. Here, I will briefly review the various formation scenarios and the respective role of nature and nurture in the formation, evolution and properties of the main structures and substructures.

We compare the results of the photometrical analysis of barred galaxies with those of a similar analysis from N-body simulations. The photometry is for a sample of nine barred galaxies observed in the J and Ks bands with the CANICA near infrared (NIR) camera at the 2.1-m telescope of the Observatorio Astrofisico Guillermo Haro (OAGH) in Cananea, Sonora, Mexico. The comparison includes radial ellipticity profiles and surface brightness (density for the N-body galaxies) profiles along the bar major and minor axes. We find very good agreement, arguing that the exchange of angular momentum within the galaxy plays a determinant role in the evolution of barred galaxies.

In this series of papers, we propose a theory to explain the formation and properties of rings and spirals in barred galaxies. The building blocks of these structures are orbits guided by the manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In this paper we focus on a comparison of the morphology of observed and of theoretical spirals and rings and we also give some predictions for further comparisons. Our theory can account for spirals as well as both inner and outer rings. The model outer rings have the observed $R_1$, $R_1'$, $R_2$, $R_2'$ and $R_1R_2$ morphologies, including the dimples near the direction of the bar major axis. We explain why the vast majority of spirals in barred galaxies are two armed and trailing, and discuss what it would take for higher multiplicity arms to form. We show that the shapes of observed and theoretical spirals agree and we predict that stronger non-axisymmetric forcings at and somewhat beyond corotation will drive more open spirals. We compare the ratio of ring diameters in theory and in observations and predict that more elliptical rings will correspond to stronger forcings. We find that the model potential may influence strongly the numerical values of these ratios.

Studying the atomic gas (HI) properties of the most isolated galaxies is essential to quantify the effect that the environment exerts on this sensitive component of the interstellar medium. We observed and compiled HI data for a well defined sample of ~ 800 galaxies in the Catalog of Isolated Galaxies, as part of the AMIGA project (Analysis of the ISM in Isolated GAlaxies, http://amiga.iaa.es), which enlarges considerably previous samples used to quantify the HI deficiency in galaxies located in denser environments. By studying the shape of 182 HI profiles, we revisited the usually accepted result that, independently of the environment, more than half of the galaxies present a perturbed HI disk. In isolated galaxies this would certainly be a striking result if these are supposed to be the most relaxed systems, and has implications in the relaxation time scales of HI disks and the nature of the most frequent perturbing mechanisms in galaxies. Our sample likely exhibits the lowest HI asymmetry level in the local Universe. We found that other field samples present an excess of ~ 20% more asymmetric HI profiles than that in CIG. Still a small percentage of galaxies in our sample present large asymmetries. Follow-up high resolution VLA maps give insight into the origin of such asymmetries.

Abr: We investigate whether the Hubble sequence can be reproduced by the relics of merger events. We verify that, at zmed=0.65, the abundant population of anomalous starbursts is mainly linked to the local spirals. Their morphologies are dominated by young stars and are related to their ionised-gas kinematics. We show that both morphologies and kinematics can be reproduced by using gas modelling from Barnes' (2002) study of major mergers. Using our modelling to estimate the gas-to-stars transformation during a merger, we identify the gas fraction in the progenitors to be generally above 50%. All distant and massive starbursts can be distributed along a temporal sequence from the first passage to the nuclei fusion and then to the disk rebuilding phase. It confirms that the rebuilding spiral disk scenario is possibly an important channel for the formation of present-day disks in spirals. Because half of the present-day spirals had peculiar morphologies and anomalous kinematics at zmed=0.65, they could indeed be in major mergers phases 6 Gyrs ago, and almost all at z~1. It is time now to study in detail the formation of spiral disks and of their substructures, including bulge, disks, arms, bars, rings that may mainly originate from instabilities created during the last major merger.

Two-armed, grand design spirals and inner and outer rings in barred galaxies can be due to orbits guided by the manifolds emanating from the vicinity of the L1 and L2 Lagrangian points, located at the ends of the bar. We first summarise the necessary theoretical background and in particular we describe the dynamics around the unstable equilibrium points in barred galaxy models, and the corresponding homoclinic and heteroclinic orbits. We then discuss two specific morphologies and the circulation of material within the corresponding manifolds. We also discuss the case where mass concentrations at the end of the bar can stabilise the L1 and L2 and the relevance of this work to the gas concentrations in spirals and rings.

In this paper we present building blocks which can explain the formation and properties both of spirals and of inner and outer rings in barred galaxies. We first briefly summarise the main results of the full theoretical description we have given elsewhere, presenting them in a more physical way, aimed to an understanding without the requirement of extended knowledge of dynamical systems or of orbital structure. We introduce in this manner the notion of manifolds, which can be thought of as tubes guiding the orbits. The dynamics of these manifolds can govern the properties of spirals and of inner and outer rings in barred galaxies. We find that the bar strength affects how unstable the L1 and L2 Lagrangian points are, the motion within the 5A5A5Amanifold tubes and the time necessary for particles in a manifold to make a complete turn around the galactic centre. We also show that the strength of the bar, or, to be more precise, of the non-axisymmetric forcing at and somewhat beyond the corotation region, determines the resulting morphology. Thus, less strong bars give rise to R1 rings or pseudorings, while stronger bars drive R2, R1R2 and spiral morphologies. We examine the morphology as a function of the main parameters of the bar and present descriptive two dimensional plots to that avail. We also derive how the manifold morphologies and properties are modified if the L1 and L2 Lagrangian points become stable. Finally, we discuss how dissipation affects the manifold properties and compare the manifolds in gas-like and in stellar cases. Comparison with observations, as well as clear predictions to be tested by observations will be given in an accompanying paper.

Abreg: By combining HST/UDF imagery with kinematics from VLT/GIRAFFE we derive a physical model of distant galaxy J033245.11-274724.0 in a way similar to what can be done in the nearby Universe. Here we study the properties of a distant compact LIRGs galaxy. Given the photometric and spectro photometric accuracies, we can decompose the galaxy in sub components and correct them for reddening. The galaxy is dominated by a dust enshrouded disk revealed by UDF imagery. The disk radius is half that of the Milky Way and the galaxy have a SFR=20Mo/yr. Morphology and kinematics show that gas and stars together spiral inwards rapidly to feed the disk and the central regions. A combined system of a bar and two non rotating spiral arms regulates the material accretion, induces large sigma, with sigma larger than 100 km/s and redistributes the angular momentum (AM). The detailed physical properties resemble to the expectations from modeling a merger of two equal mass, gaseous rich galaxies, 0.5 Gyr after the merger. In its later evolution, this galaxy could become a late type galaxy which falls on the T-F relation, with an AM mostly induced by the orbital AM of the merger.

We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters and/or constraints. This method is very powerful, to a large extent due to its simplicity. It can be used for mass distributions with an arbitrary geometry and a large variety of kinematical constraints. We present several examples illustrating it. Applications of this method include the creation of initial conditions for N-body simulations and the modelling of galaxies from their photometric and kinematic observations.

The Efstathiou, Lake and Negroponte (1982) criterion can not distinguish bar stable from bar unstable discs and thus should not be used in semi-analytic galaxy formation simulations. I discuss the reasons for this, illustrate it with examples and point out shortcomings in the recipes used for spheroid formation. I propose an alternative, although much less straightforward, possibility.

We study the role of the unstable equilibrium points in the transfer of matter in a galaxy using the potential of a rotating triaxial system. In particular, we study the neighbourhood of these points for energy levels and for main model parameters where the zero velocity curves just open and form a bottleneck in the region. For these energies, the transfer of matter from the inner to the outer parts and vice versa starts being possible. We study how the dynamics around the unstable equilibrium points is driven, by performing a partial normal form scheme and by computing the invariant manifolds of periodic orbits and quasi-periodic orbits using the reduced Hamiltonian. In particular, we compute some homoclinic and heteroclinic orbits playing a crucial role. Our results also show that in slow rotating and/or axisymmetric systems the hyperbolic character of the equilibrium points is cancelled, so that no transfer of matter is possible through the bottleneck.

In this paper, we estimate the gamma-ray and neutrino fluxes coming from dark matter annihilation in a Milky Way framework provided by a recent N-BODY HORIZON simulation. We first study the characteristics of the simulation and highlight the mass distribution within the galactic halo. The general dark matter density has a typical $r^{-3}$ power law for large radii, but the inner behaviour is poorly constrained below the resolution of the simulation ($\sim 200$ pc). We identify clumps and subclumps and analyze their distribution, as well as their internal structure. Inside the clumps, the power law is rather universal, $r^{-2.5}$ in the outer part with again strong uncertainties for smaller radii, especially for light clumps. We show a full-sky map of the astrophysical contribution to the gamma-ray or neutrino fluxes in this N-body framework. Using quite model independent and general assumptions for the high energy physics part, we evaluate the possible absolute fluxes and show some benchmark regions for the experiments GLAST, EGRET, and a km3 size extension of ANTARES like the KM3NeT project. While individual clumps seem to be beyond detection reach, the galactic center region is promising and GLAST could be sensitive to the geometry and the structure of its dark matter distribution. The detection by a km3 version of ANTARES is, however, more challenging due to a higher energy threshold. We also point out that the lack of resolution leaves the inner structure of subhalos poorly constrained. Using the same clump spectrum and mass fraction, a clump luminosity boost of order ten can be achieved with a steeper profile in the inner part of the sub-halos.

We show that stable double-frequency orbits form the backbone of double bars, because they trap around themselves regular orbits, as stable closed periodic orbits do in single bars, and in both cases the trapped orbits occupy similar volume of phase-space. We perform a global search for such stable double-frequency orbits in a model of double bars by constructing maps of trajectories with initial conditions well sampled over the available phase-space. We use the width of a ring sufficient to enclose a given map as the indicator of how tightly the trajectory is trapped around a double-frequency orbit. We construct histograms of these ring widths in order to determine the fraction of phase-space occupied by ordered motions. We build 22 further models of double bars, and we construct histograms showing the fraction of the phase-space occupied by regular orbits in each model. Our models indicate that resonant coupling between the bars may not be the dominant factor reducing chaos in the system.

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.

The class `bulges' contains objects with very different formation and evolution paths and very different properties. I review two types of `bulges', the boxy/peanut bulges (B/Ps) and the discy bulges. The former are parts of bars seen edge-on, have their origin in vertical instabilities of the disc and are somewhat shorter in extent than bars. Their stellar population is similar to that of the inner part of the disc from which they formed. Discy bulges have a disc-like outline, i.e., seen face-on they are circular or oval and seen edge-on they are thin. Their extent is of the order of 5 times smaller than that of the boxy/peanut bulges. They form from the inflow of mainly gaseous material to the centre of the galaxy and from subsequent star formation. They thus contain a lot of young stars and gas. Bulges of different types often coexist in the same galaxy. I review the main known results on these two types of bulges and present new simulation results. B/Ps form about 1Gyr after the bar, via a vertical buckling. At that time the bar strength decreases, its inner part becomes thicker -- forming the peanut or boxy shape -- and the ratio $\sigma_z^2/\sigma_r^2$ increases. A second buckling episode is seen in simulations with strong bars, also accompanied by a thickening of the peanut and a weakening of the bar. The properties of the B/Ps correlate strongly with those of the bar: stronger bars have stronger peanuts, a more flat-topped vertical density distribution and have experienced more bucklings. I also present simulations of disc galaxy formation, which include the formation of a discy bulge. Decomposition of their radial density profile into an exponential disc and a Sersic bulge gives realistic values for the disc and bulge scale-lengths and mass ratios, and a Sersic shape index of the order of 1.

I discuss results from the COSMOS survey, showing that the fraction of disc galaxies that is barred decreases considerably with look-back time from z ~ 0.2 to z ~ 0.8. This decrease is more important for small mass and low luminosity spirals. Classical bar formation theory provides a promising framework for understanding these results. I also discuss the formation of discy bulges using N-body simulations reproducing well the properties of observed discy bulges. Thus, these simulated discy bulges have the shape of a disc, they have Sersic profiles with small values of the shape index and their size is of the order of a kpc. They are formed by radial inflow of material driven by the bar and are thus composed of both gas and stars and have a considerable fraction of young stars. They can harbour spiral structure, or an inner bar.

We study the dynamics of 2D and 3D barred galaxy analytical models, focusing on the distinction between regular and chaotic orbits with the help of the Smaller ALigment Index (SALI), a very powerful tool for this kind of problems. We present briefly the method and we calculate the fraction of chaotic and regular orbits in several cases. In the 2D model, taking initial conditions on a Poincar\'{e} $(y,p_y)$ surface of section, we determine the fraction of regular and chaotic orbits. In the 3D model, choosing initial conditions on a cartesian grid in a region of the $(x, z, p_y)$ space, which in coordinate space covers the inner disc, we find how the fraction of regular orbits changes as a function of the Jacobi constant. Finally, we outline that regions near the $(x,y)$ plane are populated mainly by regular orbits. The same is true for regions that lie either near to the galactic center, or at larger relatively distances from it.

We propose a theory to explain the formation of spiral arms and of all types of outer rings in barred galaxies, extending and applying the technique used in celestial mechanics to compute transfer orbits. Thus, our theory is based on the chaotic orbital motion driven by the invariant manifolds associated to the periodic orbits around the hyperbolic equilibrium points. In particular, spiral arms and outer rings are related to the presence of heteroclinic or homoclinic orbits. Thus, R1 rings are associated to the presence of heteroclinic orbits, while R1R2 rings are associated to the presence of homoclinic orbits. Spiral arms and R2 rings, however, appear when there exist neither heteroclinic nor homoclinic orbits. We examine the parameter space of three realistic, yet simple, barred galaxy models and discuss the formation of the different morphologies according to the properties of the galaxy model. The different morphologies arise from differences in the dynamical parameters of the galaxy.

We propose a new theory to explain the formation of spiral arms and of all types of outer rings in barred galaxies. We have extended and applied the technique used in celestial mechanics to compute transfer orbits. Thus, our theory is based on the chaotic orbital motion driven by the invariant manifolds associated to the periodic orbits around the hyperbolic equilibrium points. In particular, spiral arms and outer rings are related to the presence of heteroclinic or homoclinic orbits. Thus, R1 rings are associated to the presence of heteroclinic orbits, while R1R2 rings are associated to the presence of homoclinic orbits. Spiral arms and R2 rings, however, appear when there exist neither heteroclinic nor homoclinic orbits. We examine the parameter space of three realistic, yet simple, barred galaxy models and discuss the formation of the different morphologies according to the properties of the galaxy model. The different morphologies arise from differences in the dynamical parameters of the galaxy.

We have analyzed the redshift-dependent fraction of galactic bars over 0.2<z<0.84 in 2,157 luminous face-on spiral galaxies from the COSMOS 2-square degree field. Our sample is an order of magnitude larger than that used in any previous investigation, and is based on substantially deeper imaging data than that available from earlier wide-area studies of high-redshift galaxy morphology. We find that the fraction of barred spirals declines rapidly with redshift. Whereas in the local Universe about 65% of luminous spiral galaxies contain bars (SB+SAB), at z ~0.84 this fraction drops to about 20%. Over this redshift range the fraction of strong (SB) bars drops from about 30% to under 10%. It is clear that when the Universe was half its present age, the census of galaxies on the Hubble sequence was fundamentally different from that of the present day. A major clue to understanding this phenomenon has also emerged from our analysis, which shows that the bar fraction in spiral galaxies is a strong function of stellar mass, integrated color and bulge prominence. The bar fraction in very massive, luminous spirals is about constant out to z ~ 0.84 whereas for the low mass, blue spirals it declines significantly with redshift beyond z=0.3. There is also a slight preference for bars in bulge dominated systems at high redshifts which may be an important clue towards the co-evolution of bars, bulges and black holes. Our results thus have important ramifications for the processes responsible for galactic downsizing, suggesting that massive galaxies matured early in a dynamical sense, and not just as a result of the regulation of their star formation rate.

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.

We study the effect of spiral arm passages on the evolution of star clusters on planar and circular orbits around the centres of galaxies. Individual passages with different relative velocity (V_drift) and arm width are studied using N-body simulations. When the ratio of the time it takes the cluster to cross the density wave to the crossing time of stars in the cluster is much smaller than one, the energy gain of stars can be predicted accurately in the impulsive approximation. When this ratio is much larger than one, the cluster is heated adiabatically and the net effect of heating is largely damped. For a given duration of the perturbation, this ratio is smaller for stars in the outer parts of the cluster compared to stars in the inner part. The cluster energy gain due to perturbations of various duration as obtained from our N-body simulations is in good agreement with theoretical predictions taking into account the effect of adiabatic damping. Perturbations by the broad stellar component of the spiral arms on a cluster are in the adiabatic regime and, therefore, hardly contribute to the energy gain and mass loss of the cluster. We consider the effect of crossings through the high density shocked gas in the spiral arms, which result in a more impulsive compression of the cluster. The time scale of disruption is shortest at ~0.8-0.9 R_CR since there V_drift is low. This location can be applicable to the solar neighbourhood. In addition, the four-armed spiral pattern of the Milky Way makes spiral arms contribute more to the disruption of clusters than in a similar but two-armed galaxy. Still, the disruption time due to spiral arm perturbations there is about an order of magnitude higher than what is observed for the solar neighbourhood.[ABRIDGED]

In this and in a previous paper (Romero-Gomez et al. 2006) we propose a theory to explain the formation of both spirals and rings in barred galaxies using a common dynamical framework. It is based on the orbital motion driven by the unstable equilibrium points of the rotating bar potential. Thus, spirals, rings and pseudo-rings are related to the invariant manifolds associated to the periodic orbits around these equilibrium points. We examine the parameter space of three barred galaxy models and discuss the formation of the different morphological structures according to the properties of the bar model. We also study the influence of the shape of the rotation curve in the outer parts, by making families of models with rising, flat, or falling rotation curves in the outer parts. The differences between spiral and ringed structures arise from differences in the dynamical parameters of the host galaxies. The results presented here will be discussed and compared with observations in a forthcoming paper.

We present a purely theoretical study of the morphological evolution of self-gravitating systems formed through the dissipationless collapse of N-point sources. We explore the effects of resolution in mass and length on the growth of triaxial structures formed by an instability triggered by an excess of radial orbits. We point out that as resolution increases, the equilibria shift, from mildly prolate, to oblate. A number of particles N ~= 100000 or larger is required for convergence of axial aspect ratios. An upper bound for the softening, e ~ 1/256, is also identified. We then study the properties of a set of equilibria formed from scale-free cold initial mass distributions, ro ~ r^-g with 0 <= g <= 2. Oblateness is enhanced for initially more peaked structures (larger values of g). We map the run of density in space and find no evidence for a power-law inner structure when g <= 3/2 down to a mass fraction <~0.1 per cent of the total. However, when 3/2 < g <= 2, the mass profile in equilibrium is well matched by a power law of index ~g out to a mass fraction ~ 10 per cent. We interpret this in terms of less-effective violent relaxation for more peaked profiles when more phase mixing takes place at the centre. We map out the velocity field of the equilibria and note that at small radii the velocity coarse-grained distribution function (DF) is Maxwellian to a very good approximation.

Bars in galaxies are mainly supported by particles trapped around stable periodic orbits. These orbits represent oscillatory motion with only one frequency, which is the bar driving frequency, and miss free oscillations. We show that a similar situation takes place in double bars: particles get trapped around parent orbits, which in this case represent oscillatory motion with two frequencies of driving by the two bars, and which also lack free oscillations. Thus the parent orbits, which constitute the backbone of an oscillating potential of two independently rotating bars, are the double-frequency orbits. These orbits do not close in any reference frame, but they map onto loops, first introduced by Maciejewski & Sparke (1997). Trajectories trapped around the parent double-frequency orbit map onto a set of points confined within a ring surrounding the loop.

We perform collisionless N-body simulations of 1:1 galaxy mergers, using models which include a galaxy halo, disc and bulge, focusing on the behaviour of the halo component. The galaxy models are constructed without recourse to a Maxwellian approximation. We investigate the effect of varying the galaxies' orientation, their mutual orbit, and the initial velocity anisotropy or cusp strength of the haloes upon the remnant halo density profiles and shape, as well as on the kinematics. We observe that the halo density profile (determined as a spherical average, an approximation we find appropriate) is exceptionally robust in mergers, and that the velocity anisotropy of our remnant haloes is nearly independent of the orbits or initial anisotropy of the haloes. The remnants follow the halo anisotropy - local density slope (\beta-\gamma) relation suggested by Hansen & Moore (2006} in the inner parts of the halo, but \beta is systematically lower than this relation predicts in the outer parts. Remnant halo axis ratios are strongly dependent on the initial parameters of the haloes and on their orbits. We also find that the remnant haloes are significantly less spherical than those described in studies of simulations which include gas cooling.

We use the Smaller ALignment Index (SALI) method of chaos detection, to study the global dynamics of conservative dynamical systems described by differential or difference equations. In particular, we consider the well--known 2 and 4--dimensional symplectic standard map, as well as an autonomous Hamiltonian system of 2 and 3 degrees of freedom describing the motion of a star in models of barred galaxies. The application of SALI helped us to compute rapidly and accurately the percentage of regular and chaotic motion for particular values of the parameters of these systems. We were also able to perform a computationally efficient determination of the dependence of these percentages on the variation of several parameters of the studied models.

N-body simulations argue that the inner haloes of barred galaxies should not be spherical, nor even axisymmetric, but triaxial. The departure from sphericity is strongest near the centre and decreases outwards; typical axial ratios for the innermost parts are of the order of 0.8. The halo shape is prolate-like in the inner parts up to a certain radius and then turns to oblate-like. I call this inner halo structure the `halo bar' and analyse here in depth its structure and kinematics in a representative model. It is always considerably shorter than the disc bar. It lags the disc bar by only a few degrees at all radii and the difference between the two bar phases increases with distance from the centre. The two bars turn with roughly the same pattern speed. This means that the halo bar is a slow bar, since its corotation radius is much larger than its length. The bisymmetric component in the halo continues well outside the halo bar in the form of an open spiral, trailing behind the disc bar. The inner parts of the halo display some mean rotation in the same sense as the disc rotation. This is more important for particles nearer to the equatorial plane and decreases with increasing distance from it, but is always much smaller than the disc rotation.

We acquired radial velocities of a significant number of globular clusters (GCs) on wide fields between galaxies in the nearby Fornax cluster of galaxies, in order to derive their velocity dispersion radial profile and to probe the dynamics of the cluster. We used FLAMES on the VLT to obtain accurate velocities for 149 GCs, within a ~500x150 kpc strip centered on NGC 1399, the Fornax central galaxy. These objects are at the very bright tail (M_V < -9.5) of the GC luminosity function, overlapping the so-called ``ultra-compact dwarfs'' magnitude range. Eight of the brightest FLAMES-confirmed members indeed show hints of resolution in the subarcsecond pre-imaging data we used for selecting the ~500 targets for FLAMES spectroscopy. Ignoring the GCs around galaxies by applying 3d_25 diameter masks, we find 61 GCs of 20.0 < V < 22.2 lying in the intra-cluster (IC) medium. The velocity dispersion of the population of ICGCs is 200 km/s at ~150 kpc from the central NGC 1399 and rises to nearly 400 km/s at 200 kpc, a value which compares with the velocity dispersion of the population of dwarf galaxies, thought to be infalling from the surroundings of the cluster.

I discuss several aspects of secular evolution linked to bars and to boxy/peanut bulges, based on a very large number of high resolution, fully self-consistent $N$-body simulations. When the bar forms, it is as thin as the disc. Its three-dimensional shape, however, evolves, so that, at later times, it has a thick inner part and a thin, more extended outer part. The former, when viewed edge-on, is called a boxy/peanut bulge, because of its shape. The strength of the box/peanut correlates with the bar strength, the strongest cases having formed after two buckling episodes. The extent of the box/peanut is considerably shorter than the bar length, in good agreement with orbital structure studies and with observations. Viewed at an angle near to, but not quite edge-on, barred galaxies show specific isodensity/isophotal shapes, which are different in the thick and in the thin part of the bar. The isophotes of M31 also have such shapes. This, taken together with radial photometric profiles and kinematics, argue that M31 is a barred galaxy. Thus, the pseudo-ring seen at roughly 50' could be an outer ring formed at the outer Lindblad resonance of the bar.

We investigate encounters between giant molecular clouds (GMCs) and star clusters. We propose a single expression for the energy gain of a cluster due to an encounter with a GMC, valid for all encounter distances and GMC properties. This relation is verified with N-body simulations of cluster-GMC encounters and excellent agreement is found. The fractional mass loss from the cluster is 0.25 times the fractional energy gain. This is because 75% of the injected energy goes to the velocities of escaping stars, that are higher than the escape velocity. We derive an expression for the cluster disruption time (t_dis) based on the mass loss from the simulations, taking into account the effect of gravitational focusing by the GMC. The disruption time depends on the cluster mass (M_c) and half-mass radius (r_h) as t_dis=2.0 S (M_c/10^4 M_sun)(3.75 pc/r_h)^3 Gyr, with S=1 for the solar neighbourhood and inversely proportional to the GMC density. The observed shallow relation between cluster radius and mass gives t_dis a power-law dependence on the mass with index 0.7, similar to that found from observations and from simulations of clusters dissolving in tidal fields (0.62). The constant of 2.0 Gyr is about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving under the combined effect of galactic tidal field and stellar evolution. It is somewhat higher than the observationally determined value of 1.3 Gyr. It suggests, however, that the combined effect of tidal field and encounters with GMCs can explain the lack of old open clusters in the solar neighbourhood. GMC encounters can also explain the (very) short disruption time that was observed for star clusters in the central region of M51, since there rho_n is an order of magnitude higher than in the solar neighbourhood.

Kn-band images, unsharp-masked images, as well as major-axis and vertically- summed surface brightness profiles are presented for 30 edge-on spiral galaxies, most with a boxy or peanut-shaped (B/PS) bulge. Such galaxies have more complex morphologies than galaxies of other bulge types, more often showing (off-)centered X structures, secondary major-axis maxima and spiral-like structures. Those features are also observed in N-body simulations of barred discs and may trace the main bar orbit families. The surface brightness profiles of galaxies with a B/PS bulge are also more complex, with typically 3 or more clearly separated regions, including a flat intermediate region (Freeman Type II profiles). Those radial breaks offer evidence for bar-driven transfer of angular momentum and radial redistribution of material. The profiles also suggest a rapid variation of the scaleheight of the disc material, contrary to conventional wisdom but again as expected from vertical resonances and instabilities in barred discs. The steep inner region of the surface brightness profiles is often shorter than the isophotally thick part of the galaxies, itself always shorter than the flat region of the profiles. Contrary to the standard `bulge + disc' model, we thus propose that galaxies with a B/PS bulge are composed of a thin concentrated disc (a disc-like bulge) contained within a partially thick bar (the B/PS bulge) and a thin outer disc. The inner disc likely formed secularly through bar-driven processes and is responsible for the steep inner region of the surface brightness profiles, while the bar is responsible for the flat region and the thick complex morphological structures observed. Those components are strongly coupled dynamically and are formed mostly of the same (disc) material. [Abridged]

The inclination of M31 is too close to edge-on for a bar component to be easily recognised and is not sufficiently edge-on for a boxy/peanut bulge to protrude clearly out of the equatorial plane. Nevertheless, a sufficient number of clues allow us to argue that this galaxy is barred. We use fully self-consistent N-body simulations of barred galaxies and compare them with both photometric and kinematic observational data for M31. In particular, we rely on the near infrared photometry presented in a companion paper. We compare isodensity contours to isophotal contours and the light profile along cuts parallel to the galaxy major axis and offset towards the North, or the South, to mass profiles along similar cuts on the model. All these comparisons, as well as position velocity diagrams for the gaseous component, give us strong arguments that M31 is barred. We compare four fiducial N-body models to the data and thus set constraints on the parameters of the M31 bar, as its strength, length and orientation. Our `best' models, although not meant to be exact models of M31, reproduce in a very satisfactory way the main relevant observations. We present arguments that M31 has both a classical and a boxy/peanut bulge. Its pseudo-ring-like structure at roughly 50' is near the outer Lindblad resonance of the bar and could thus be an outer ring, as often observed in barred galaxies. The shape of the isophotes also argues that the vertically thin part of the M31 bar extends considerably further out than its boxy bulge, i.e. that the boxy bulge is only part of the bar, thus confirming predictions from orbital structure studies and from previous N-body simulations.

A new, 2.8 deg^2 J,H,K_s infrared survey from the 2MASS 6x program across the extent of the optical disk of the Andromeda (M31) galaxy provides a clear view of the M31 center almost completely unfettered by dust extinction, and reveals a high contrast bulge with extremely boxy isophotes dominating the NIR light to a semi-major axis of ~700''(2.6 kpc). The inner bulge (<~50'') is relatively circular, but shows some isophotal twisting. Beyond this, the ellipticity and boxiness of the bulge increase with radius -- achieving a boxiness that rivals that of any other known disk galaxy observed in the near infrared -- and the position angle is constant at ~50 deg, which is about 10 deg higher than the position angle of the M31 disk. Boxy bulges in highly inclined disks have been shown to be the vertical structure of bars, and self-consistent, N-body modeling specific to the NIR images presented here can reproduce the observed NIR M31 features with a combination of a classical bulge and a boxy bulge/bar. Beyond the boxy bulge region and nearly along the 40 deg position angle of the disk a narrow ridge of infrared flux, which can be identified with the thin part of the bar, more or less symmetrically extends into the inner disk at semi-major axis radii of 700'' to 1200'' or more. Little variation in the morphology or relative brightnesses of these various M31 structures is seen across the NIR bands (i.e., no color gradients are seen). These new data verify that M31 is a barred spiral galaxy like the Milky Way.

We study the encounters between stars clusters and giant molecular clouds (GMCs). The effect of these encounters has previously been studied analytically for two cases: 1) head-on encounters, for which the cluster moves through the centre of the GMC and 2) distant encounters, where the encounter distance p > 3*R_n, with p the encounter parameter and R_n the radius of the GMC. We introduce an expression for the energy gain of the cluster due to GMC encounters valid for all values of p and R_n. This analytical result is confronted with results from N-body simulations and excellent agreement is found. From the simulations we find that the fractional mass loss is only 25% of the fractional energy gain. This is because stars escape with velocities much higher than the escape velocity. Based on the mass loss, we derive a disruption time for star clusters due to encounters with GMCs of the form t_dis [Gyr] = 2.0*S*(M_c/10^4 M_sun)^gamma, with S=1 for the solar neighbourhood and inversely proportional with the global GMC density and gamma=1-3lambda, with lambda the index that relates the cluster half-mass radius to the cluster mass (r_h ~ M_c^lambda). The observed shallow relation between cluster radius and mass (e.g. lambda=0.1), makes the index (gamma=0.7) similar to the index found both from observations and from simulations of clusters dissolving in tidal fields (gamma=0.62). The constant of 2.0 Gyr, which is the disruption time of a 10^4 M_sun cluster in the solar neighbourhood, is close to the value of 1.3 Gyr which was empirically determined from the age distribution of open clusters. This suggests that the combined effect of GMC encounters, stellar evolution and galactic tidal field can explain the lack of old open clusters in the solar neighbourhood.

We propose a new theory for the formation of rR_1 ring structures, i.e. for ring structures with both an inner and an outer ring, the latter having the form of ``8''. We propose that these rings are formed by material from the stable and unstable invariant manifolds associated with the Lyapunov orbits around the equilibrium points of a barred galaxy. We discuss the shape and velocity structure of the rings thus formed and argue that they are in agreement with the observed properties of rR_1 structures.

The chaotic or ordered character of orbits in galactic models is an important issue, since it can influence dynamical evolution. This distinction can be achieved with the help of the Smaller Alingment Index - (SALI). We describe here briefly this method and its advantages. Then we apply it to that case of 2D and 3D barred galaxy potentials. In particular, we find the fraction of chaotic and ordered orbits in such potentials and present how this fraction changes when the main parameters of the model are varied. For this, we consider models with different bar mass, bar thickness or pattern speed. Varying only one parameter at a time, we find that bars that are more massive, or thinner, or faster, have a larger fraction of chaotic orbits.

A very important issue in the area of galactic dynamics is the detection of chaotic and ordered motion inside galaxies. In order to achieve this target, we use the Smaller ALignment Index (SALI) method, which is a very suitable tool for this kind of problems. Here, we apply this index to 3D barred galaxy potentials and we present some results on the chaotic behavior of the model when its main parameters vary.

The dynamical evolution of barred galaxies depends crucially on the fraction and their spacial distribution of chaotic orbits in them. In order to distinguish between the two kinds of orbits, we use the Smaller Alignment Index (SALI) method, a very powerful method which can be applied to problems of galactic dynamics. Using model potentials, and taking into account the full 3D distribution of matter, we discuss how the distribution of chaotic orbits depends on the main model parameters, like the mass of the various components and the bar axial ratio.

From this vast subject, I will pick out and review three specific topics, namely the formation and evolution of bars, the formation of bulges, and the evolution during multiple major mergers. Bars form naturally in galactic discs. Their evolution is driven by the exchange of angular momentum within the galaxy. This is emitted mainly by near-resonant material in the inner disc (bar), and is absorbed by near-resonant material in the outer disc and in the halo. As a result of this, the bar becomes stronger and rotates slower. Bulges are not a homogeneous class of objects. Based on their formation history, one can distinguish three types. Classical bulges are mainly formed before the actual disc component, from collapses or mergers and the corresponding dissipative processes. Boxy/peanut bulges are parts of bars seen edge-on. Finally, disc-like bulges are formed by the inflow of material to the center due to bar torques. Major mergers bring strong and fast evolution and can turn discs into ellipticals. I present results from simulations of multiple mergers in groups of disc or of elliptical galaxies and discuss the orbital anisotropy in the merger remnant.

The halo plays a crucial role in the evolution of barred galaxies. Its near-resonant material absorbs angular momentum emitted from some of the disc particles and helps the bar become stronger. As a result, a bar (oval) forms in the inner parts of the halo of strongly barred disc galaxies. It is thinner in the inner parts (but still considerably fatter than the disc bar) and tends to spherical at larger radii. Its length increases with time, while always staying shorter than the disc bar. It is roughly aligned with the disc bar, which it trails only slightly, and it turns with roughly the same pattern speed. The bi-symmetric component of the halo density continues well outside the halo bar, where it clearly trails behind the disc bar. The length and strength of the disc and halo bars correlate; the former being always much stronger than the latter. If the halo is composed of weakly interacting massive particles, then the formation of the halo bar, by redistributing the matter in the halo and changing its shape, could influence the expected annihilation signal. This is indeed found to be the case if the halo has a core, but not if it has a steep cusp. The formation and evolution of the bar strongly affect the halo orbits. A fraction of them becomes near-resonant, similar to the disc near-resonant orbits at the same resonance, while another fraction becomes chaotic. Finally, a massive and responsive halo makes it harder for a central mass concentration to destroy the disc bar.

We study the impact of the halo shape and geometry on the expected weakly interacting massive particle (WIMP) dark matter annihilation signal from the galactic center. As the halo profile in the innermost region is still poorly constrained, we consider different density behaviors like flat cores, cusps and spikes, as well as geometrical distortions. We show that asphericity has a strong impact on the annihilation signal when the halo profile near the galactic center is flat, but becomes gradually less significant for cuspy profiles, and negligible in the presence of a central spike. However, the astrophysical factor is strongly dependent on the WIMP mass and annihilation cross-section in the latter case.

We have extended previous analytical studies of the interaction of dark halos with galactic disks by introducing for the halo particles anisotropic distribution functions in phase space. For this purpose we have employed the shearing sheet model of a patch of a galactic disk embedded in a homogeneous halo. We find that velocity anisotropy increases considerably the maximum growth factor of perturbations in the disk.

Kn-band imaging of a sample of 30 edge-on spiral galaxies with a boxy or peanut-shaped (B/PS) bulge is discussed. Galaxies with a B/PS bulge tend to have a more complex morphology than galaxies with other bulge types, unsharp-masked images revealing structures that trace the major orbit families of three-dimensional bars. Their surface brightness profiles are also more complex, typically containing 3 or more clearly separated regions, including a shallow or flat intermediate region (Freeman Type II profiles), suggestive of bar-driven transfer of angular momentum and radial redistribution of material. The data also suggest abrupt variations of the discs' scaleheights, as expected from the vertical resonances and instabilities present in barred discs but contrary to conventional wisdom. Counter to the standard `bulge + disc' model, we thus propose that galaxies with a B/PS bulge are composed of a thin concentrated disc (a disc-like bulge) contained within a partially thick bar (the B/PS bulge), itself contained within a thin outer disc. The inner disc most likely formed through bar-driven processes while the thick bar arises from buckling instabilities. Both are strongly coupled dynamically and are formed mostly of the same (disc) material.

We study the effect of a central mass concentration (CMC) on the secular evolution of a barred disc galaxy. Unlike previous studies, we use fully self-consistent 3D $N$-body simulations with live haloes, which are known to be important for bar evolution. The CMC is introduced gradually, to avoid transients. In all cases where the mass of the CMC is of the order of, or more than, a few per cent of the mass of the disc, the strength of the bar decreases noticeably. The amount of this decrease depends strongly on the bar type. For the same CMC, bars with exponential surface-density profile, which formed in a disk-dominated galaxy (MD-type bars), can be totally destroyed, while strong bars with a flat surface-density profile, whose evolution is largely due to the halo (MH-type bars), witness only a decrease of their strength. This decrease occurs simultaneously from both the innermost and outermost parts of the bar. The CMC has a stronger effect on the Fourier components of higher azimuthal wave number $m$, leading to fatter and/or less rectangular bars. Furthermore, the CMC changes the side-on outline from peanut-shaped to boxy or, for massive CMCs, to elliptical. Similarly, side-on initially boxy outlines can be destroyed. The CMC also influences the velocity dispersion profiles. Most of the decrease of the bar strength occurs while the mass of the CMC increases and it is accompanied by an increase of the pattern speed. In all our simulations, the mass of the CMC necessary in order to destroy the bar is at least several per cent of the mass of the disc. This argues that observed super-massive black holes are not likely to destroy pre-existing bars.

We present a HI synthesis imaging study of NGC 864 (CIG 96), a spiral galaxy well isolated from similarly sized companions, yet presenting an intriguing asymmetry in its integral HI spectrum. The asymmetry in the HI profile is associated with a strong kinematical perturbation in the gaseous envelope of the galaxy, where at one side the decay of the rotation curve is faster than keplerian. We detect a small (M(HI) = 5 x 10^6 Msun) galaxy with a faint optical counterpart at around 80 kpc projected distance from NGC 864. This galaxy is probably not massive enough to have caused the perturbations in NGC 864. We discuss alternatives, such as the accretion of a gaseous companion at a radial velocity lower than the maximum.

A live halo plays an active role in the formation and evolution of bars by participating in the angular momentum redistribution which drives the dynamical evolution. Angular momentum is emitted mainly by near-resonant material in the bar region and is absorbed mainly by near-resonant material in the halo and in the outer disc. This exchange determines the strength of the bar, the decrease of its pattern speed, as well as its morphology. Thus, contrary to previous beliefs, a halo can help the bar grow, so that bars growing in galaxies with responsive massive haloes can become stronger than bars growing in disc dominated galaxies. During the evolution the halo does not stay axisymmetric. It forms a bar which is shorter and fatter than the disc bar and stays so all through the simulation, although its length grows considerably with time. I discuss the orbital structure in the disc and the halo and compare it with periodic orbits in analytical barred galaxy potentials. A central mass concentration (e.g. a central black hole, or a central disc) weakens a bar and increases its pattern speed. The effect of the central mass concentration depends strongly on the model, being less strong in models with a massive concentrated halo and a strong bar.

Present in over 45% of local spirals, boxy and peanut-shaped bulges are generally interpreted as edge-on bars and may represent a key phase in the evolution of bulges. Aiming to test such claims, the kinematic properties of self-consistent 3D N-body simulations of bar-unstable disks are studied. Using Gauss-Hermite polynomials to describe the stellar kinematics, a number of characteristic bar signatures are identified in edge-on disks: 1) a major-axis light profile with a quasi-exponential central peak and a plateau at moderate radii (Freeman Type II profile); 2) a ``double-hump'' rotation curve; 3) a sometime flat central velocity dispersion peak with a plateau at moderate radii and occasional local central minimum and secondary peak; 4) an h3-V correlation over the projected bar length. All those kinematic features are spatially correlated and can easily be understood from the orbital structure of barred disks. They thus provide a reliable and easy-to-use tool to identify edge-on bars. Interestingly, they are all produced without dissipation and are increasingly realized to be common in spirals, lending support to bar-driven evolution scenarios for bulge formation. So called ``figure-of-eight'' position-velocity diagrams are never observed, as expected for realistic orbital configurations. Although not uniquely related to triaxiality, line-of-sight velocity distributions with a high velocity tail (i.e. an h3-V correlation) appear as particularly promising tracers of bars. The stellar kinematic features identified grow in strength as the bar evolves and vary little for small inclination variations. Many can be used to trace the bar length. Comparisons with observations are encouraging and support the view that boxy and peanut-shaped bulges are simply thick bars viewed edge-on.

Objects designated as bulges in disc galaxies do not form a homogeneous class. I distinguish three types. The classical bulges, whose properties are similar to those of ellipticals and which form by collapse or merging. Boxy and peanut bulges, which are seen in near edge-on galaxies and which are in fact just a part of the bar seen edge-on. Finally disc-like bulges, which result from the inflow of (mainly) gas to the center-most parts, and subsequent star formation. I make a detailed comparison of the properties of boxy and peanut bulges with those of $N$-body bars seen edge-on and answer previously voiced objections about the links between the two. I also present and analyse simulations where a boxy /peanut feature is present at the same time as a classical spheroidal bulge and compare them with observations. Finally, I propose a nomenclature that can help distinguish between the three types of bulges and avoid considerable confusion.

Angular momentum redistribution within barred galaxies drives their dynamical evolution. Angular momentum is emitted mainly by near-resonant material in the bar region and absorbed by resonant material mainly in the outer disc and in the halo. This exchange determines the strength of the bar, the decrease of its pattern speed, as well as its morphology. If the galaxy has also a gaseous component and/or a companion or satellite, then these also take part in the angular momentum exchange. During the evolution a bar structure forms in the inner parts of the halo as well. This bar is shorter and fatter than the disc bar and stays so all through the simulation, although its length grows considerably with time. Viewed edge-on, the bar in the disc component acquires a boxy or peanut shape. I describe the families of periodic orbits that explain such structures and review the observations showing that boxy/peanut `bulges' are in fact just bars seen edge-on.

A multi-faceted approach is described to constrain the importance of bar-driven evolution in disk galaxies, particularly bulge formation. N-body simulations are used to construct stellar kinematic bar diagnostics for edge-on systems and to quantify the expected vertical structure of bars, and they are compared to observations of 30 edge-on spirals, most with a boxy bulge. Long-slit spectra of the galaxies show characteristic double-hump rotation curves, dispersion profiles with secondary peaks and/or flat maxima, and correlated h3 and V profiles, indicating that most of them harbor edge-on bars. The presence of cold, quasi-axisymmetric central stellar disks is also suggested, presumably formed through bar-driven gaseous inflow and star formation. K-band imaging of the same galaxies spectacularly highlights radial variations of the bars' scaleheights, as expected from vertical disk instabilities. The light profiles also vary radially in shape but never approach a classic de Vaucouleurs law. Filtering of the images further isolates the specific orbit families at the origin of the boxy structure, which can be directly related to periodic orbit calculations in 3D barred potentials. Bars are thus shown to contribute substantially to the formation of both large-scale triaxial bulges and embedded central disks. Relevant results from the SAURON survey of the stellar/ionized-gas kinematics and stellar populations of spheroids are also described. Examples are used to illustrate the potential of coupling stellar kinematics and linestrengths (age and metallicity), here specifically to unravel the dynamical evolution and related chemical enrichment history of bars and bulges. [Abridged]

Isolated barred galaxies evolve by redistributing their internal angular momentum, which is emitted mainly at the inner disc resonances and absorbed mainly at the resonances in the outer disc and the halo. This causes the bar to grow stronger and its pattern speed to decrease with time. A massive, responsive halo enhances this process. I show correlations and trends between the angular momentum absorbed by the halo and the bar strength, pattern speed and morphology. It is thus possible to explain why some disc galaxies are strongly barred, while others have no bar, or only a short bar or an oval. In some cases, a bar is found also in the halo component. This ``halo bar'' is triaxial, but more prolate-like, is shorter than the disc bar and rotates with roughly the same pattern speed. I finally discuss whether bars can modify the density cusps found in cosmological CDM simulations of dark matter haloes.

We study the regeneration of stellar bars triggered by a tidal interaction, using numerical simulations of either purely stellar or stellar+gas disc galaxies. We find that interactions which are sufficiently strong to regenerate the bar in the purely stellar models do not lead to a regeneration in the dissipative models, owing to the induced gas inflow in those models. In models in which the bar can be regenerated, we find a tight correlation between the strength and the pattern speed of the induced bar. This relation can be explained by a significant radial redistribution of angular momentum in the disc due to the interaction, similar to the processes and correlations found for isolated barred spirals. We furthermore show that the regenerated bars show the same dynamical properties as their isolated counterparts.

We argue that inner rings in barred spiral galaxies are associated with specific 2D and 3D families of periodic orbits located just beyond the end of the bar. These are families located between the inner radial ultraharmonic 4:1 resonance and corotation. They are found in the upper part of a type-2 gap of the x1 characteristic, and can account for the observed ring morphologies without any help from families of the x1-tree. Due to the evolution of the stability of all these families, the ring shapes that are favored are mainly ovals, as well as polygons with `corners' on the minor axis, on the sides of the bar. On the other hand pentagonal rings, or rings of the NGC 7020 type hexagon, should be less probable. The orbits that make the rings belong in their vast majority to 3D families of periodic orbits and orbits trapped around them.

Isolated barred galaxies evolve by redistributing their angular momentum, which, emitted by material in the inner disc at resonance with the bar, can be absorbed by resonant material in the outer disc, or in the halo. The amount of angular momentum that can be emitted/absorbed at a given resonance depends on the distribution function of the emitting/absorbing material. It thus depends not only on the amount of material on resonant orbits, but also on the velocity dispersion of that material. As it loses angular momentum, the bar becomes stronger and it also rotates slower. Thus the strength of the bar and the decrease of its pattern speed with time are set by the amount of angular momentum exchanged within the galaxy, which, in turn, is regulated by the mass distribution and the velocity dispersion of the material in the disc and spheroidal components. Correlations between the pattern speed of the bar, its strength and the angular momentum absorbed by the spheroid (halo plus bulge) argue strongly that it is the amount of angular momentum exchanged that determines the strength and the slowdown rate of the bar. The decrease of the bar pattern speed with time should not be used to set constraints on the halo-to-disc mass ratio, since it depends also on the velocity dispersion of the halo and disc material.

Angular momentum exchange is a driving process for the evolution of barred galaxies. Material at resonance in the bar region loses angular momentum which is taken by material at resonance in the outer disc and/or the halo. By losing angular momentum, the bar grows stronger and slows down. This evolution scenario is backed by both analytical calculations and by $N$-body simulations. The morphology of the bar also depends on the amount of angular momentum exchanged.

We study the conditions that favour boxiness of isodensities in the face-on views of orbital 3D models for barred galaxies. Using orbital weighted profiles we show that boxiness is in general a composite effect that appears when one considers stable orbits belonging to several families of periodic orbits. 3D orbits that are introduced due to vertical instabilities, play a crucial role in the face-on profiles and enhance their rectangularity. This happens because at the 4:1 radial resonance region we have several orbits with boxy face-on projections, instead of few rectangular-like x1 orbits, which, in a fair fraction of the models studied so far, are unstable at this region. Massive bars are characterized by rectangular-like orbits. However, we find that it is the pattern speed that affects most the elongation of the boxy feature, in the sense that fast bars are more elongated than slow ones. Boxiness in intermediate distances between the center of the model and the end of the bar can be attributed to x1v1 orbits, or to a combination of families related to the radial 3:1 resonance.

We investigate the dynamical effects of an interaction between an initially barred galaxy and a small spherical companion using an N-body/SPH algorithm. In the models described here the small companion passes through the disc of the larger galaxy near-perpendicular to its plane. The impact positions and times are varied with respect to the phase of the bar and the dynamical evolution of the disc. The interactions produce expanding ring structures, offset bars, spokes,and other asymmetries in the stars and gas. These characteristic signatures of the interaction are present in the disc for about 1 Gyr. We find that in some cases it is possible to destroy the bar while keeping the disc structure. In general the central impacts cause larger damage to the bar and the disc than the peripheral ones. The interaction tends to accelerate the transition from a strongly barred galaxy to a weakly or non-barred galaxy. The final disc morphology is determined more by the impact position relative to the bar rather than the impact time.

We present a near-infrared K_n-band photometric study of edge-on galaxies with a box/peanut-shaped `bulge'. The morphology of the galaxies is analysed using unsharp masking and fits to the vertical surface brightness profiles, and the results are compared to N-body simulations and orbital calculations of barred galaxies. Both theoretical approaches reproduce the main structures observed.

I follow a bar from its formation, via its evolution, to its destruction and, perhaps, regeneration. I discuss the main features at each stage and particularly the role of the halo. Bars can form even in sub-maximum discs. In fact, such bars can be stronger than bars which have grown in maximum discs. This is due to the response of the halo and, in particular, to the exchange of energy and angular momentum between the disc particles constituting the bar and the halo particles at resonance with it. The bar slowdown depends on the initial central concentration of the halo and the initial value of the disc Q. Contrary to the halo mass distribution, the disc changes its radial density profile considerably during the evolution. Applying the Sackett criterion, I thus find that discs become maximum in many simulations in which they have started off as sub-maximum. I briefly discuss the evolution if a gaseous component is present, as well as the destruction and regeneration of bars.

We present families, and sets of families, of periodic orbits that provide building blocks for boxy and peanut (hereafter b/p) edge-on profiles. We find cases where the b/p profile is confined to the central parts of the model and cases where a major fraction of the bar participates in this morphology. A b/p feature can be built either by 3D families associated with 3D bifurcations of the x1 family, or, in some models, even by families related with the z-axis orbits and existing over large energy intervals. The {\sf `X'} feature observed inside the boxy bulges of several edge-on galaxies can be attributed to the peaks of successive x1v1 orbits (Skokos et al. 2002a, hereafter paper I), provided their stability allows it. However in general, the x1v1 family has to overcome the obstacle of a S\ar\D\ar S transition in order to support the structure of a b/p feature. Other families that can be the backbones of b/p features are x1v4 and z3.1s. The morphology and the size of the boxy or peanut-shaped structures we find in our models is determined by the presence and stability of the families that support b/p features. The present study favours the idea that the observed edge-on profiles are the imprints of families of periodic orbits that can be found in appropriately chosen Hamiltonian systems, describing the potential of the bar.

Boxy edge-on profiles can be accounted for not only in models of barred galaxies, but also in models of normal (non-barred) galaxies. Thus, the presence of a bar is not a sine qua non condition for the appearance of this feature, as often assumed. We show that a `boxy' or a `peanut' structure in the central parts of a model is due to the presence of vertical resonances at which stable families of periodic orbits bifurcate from the planar x1 family. The orbits of these families reach in their projections on the equatorial plane a maximum distance from the center, beyond which they increase their mean radii by increasing only their deviations from the equatorial plane. The resulting orbital profiles are `stair-type' and constitute the backbone for the observed boxy structures in edge-on views of $N$-body models and, we believe, in edge-on views of disc galaxies. Since the existence of vertical resonances is independent of barred or spiral perturbations in the disc, `boxy' profiles may appear also in almost axisymmetric cases.

NGC 3642 was classified as a spiral galaxy with three rings and no bar. We have performed an HI and optical study of this nearly face-on galaxy. We find that the nuclear ring might in fact be part of an inner one-armed spiral, that could be driving nuclear accretion and feeding the central activity in the inner kpc. The inner ring is faint, and the outer ring is a rather ill-defined pseudoring. Furthermore, the size ratio of the rings is such that they cannot be due to a single pattern speed linking them together. The outer pseudoring is peculiar, since it lies in the faint outer parts of the disk, where star formation is still going on at 1.4 times the optical radius. Higher HI column densities are associated with these regions and the atomic gas layer is warped. These perturbations affect only the outer disk, since the kinematics within the main body conforms well to an ordinary differentially rotating disk. We propose here that both nuclear activity and star formation in the warped outer parts might be linked to the fact that NGC 3642 is located in a rich environment, where its close neighbors show clear signs of merging. Our suggestion is that NGC 3642 has captured recently a low-mass, gas-rich dwarf, and star formation was triggered in this infalling external gas that produced also a pronounced warp in the gaseous disk.

The purpose of this article is to show that when dynamically cold, dissipationless self-gravitating systems collapse, their evolution is a strong function of the symmetry in the initial distribution. We explore with a set of pressure-less homogeneous fluids the time-evolution of ellipsoidal distributions and map the depth of potential achieved during relaxation as function of initial ellipsoid axis ratios. We then perform a series of $N$-body numerical simulations and contrast their evolution with the fluid solutions. We verify an analytic relation between collapse factor ${\cal C}$ and particle number $N$ in spherical symmetry, such that ${\cal C} \propto N^{1/3}$. We sought a similar relation for axisymmetric configurations, and found an empirical scaling relation such that ${\cal C} \propto N^{1/6}$ in these cases. We then show that when mass distributions do not respect spherical- or axial-symmetry, the ensuing gravitational collapse deepens with increasing particle number $N$ but only slowly: 86% of triaxial configurations may collapse by a factor of no more than 40 as $N\to\infty$. For $N\approx 10^5$ and larger, violent relaxation develops fully under the Lin-Mestel-Shu instability such that numerical $N$-body solutions now resolve the different initial morphologies adequately.

We have compiled a new sample of galaxies with published catalogs of HII region coordinates. This sample, together with the former catalog of Garcia-Gomez and Athanassoula (1991) will form the basis for subsequent studies of the spiral structure in disc galaxies. In this paper we address the problem of the deprojection of the galaxy images. For this purpose we use two deprojection methods based on the HII region distribution and compare the results with the values found in the literature using other deprojection methods. Taking into account the results of all the methods, we propose optimum values for the position and inclination angles of all the galaxies in our sample.

In this series of papers we investigate the orbital structure of 3D models representing barred galaxies. In the present introductory paper we use a fiducial case to describe all families of periodic orbits that may play a role in the morphology of three-dimensional bars. We show that, in a 3D bar, the backbone of the orbital structure is not just the x1 family, as in 2D models, but a tree of 2D and 3D families bifurcating from x1. Besides the main tree we have also found another group of families of lesser importance around the radial 3:1 resonance. The families of this group bifurcate from x1 and influence the dynamics of the system only locally. We also find that 3D orbits elongated along the bar minor axis can be formed by bifurcations of the planar x2 family. They can support 3D bar-like structures along the minor axis of the main bar. Banana-like orbits around the stable Lagrangian points build a forest of 2D and 3D families as well. The importance of the 3D x1-tree families at the outer parts of the bar depends critically on whether they are introduced in the system as bifurcations in $z$ or in $\dot{z}$.

We investigate the orbital structure in a class of 3D models of barred galaxies. We consider different values of the pattern speed, of the strength of the bar and of the parameters of the central bulge of the galactic model. The morphology of the stable orbits in the bar region is associated with the degree of folding of the x1-characteristic. This folding is larger for lower values of the pattern speed. The elongation of rectangular-like orbits belonging to x1 and to x1-originated families depends mainly on the pattern speed. The detailed investigation of the trees of bifurcating families in the various models shows that major building blocks of 3D bars can be supplied by families initially introduced as unstable in the system, but becoming stable at another energy interval. In some models without radial and vertical 2:1 resonances we find, except for the x1 and x1-originated families, also families related to the z-axis orbits, which support the bar. Bifurcations of the x2 family can build a secondary 3D bar along the minor axis of the main bar. This is favoured in the slow rotating bar case.

I show that strong bars can grow in galactic discs, even when the latter are immersed in haloes whose mass within the disc radius is comparable to, or larger than, the mass of the disc. I argue that this is due to the response of the halo and in particular to the destabilising influence of the halo resonant stars. Via this instability mechanism the halo can stimulate, rather than restrain, the growth of the bar.

We discuss the evolution of a disc galaxy due to the formation of a bar and, subsequently, a peanut. After the formation stage there is still considerable evolution, albeit slower. In purely stellar cases the pattern speed of the bar decreases with time, while its amplitude grows. However, if a considerable gaseous component is present in the disc, the pattern speed may increase with time, while the bar strength may decrease. In some cases the gas can be brought sufficiently close to the center to create a strong central concentration, which, in turn, may modify the properties of the bar. More violent evolution can take place during interactions, so that some disc substructures can be either formed or destroyed in a time scale which is small compared to a Hubble time. These include spirals, bars, bridges, tails, rings, thick discs and bulges. In some cases interactions may lead to mergings. We briefly review comparisons of the properties of merger remnants with those of elliptical galaxies, both for the case of pairwise mergings and the case of multiple mergings.

In this paper I argue that, far from necessarily hindering bar formation in disc galaxies, inner haloes may stimulate it. This constitutes a new instability mechanism by which bars can grow. To show this I use a number of N-body simulations whose initial conditions have identical discs and more or less concentrated haloes. They show that the bar that grows in the more halo-dominated environment is considerably stronger than the bar that grows in the more disc-dominated environment. This result is obtained from simulations with live haloes, i.e. composed of particles which respond to the disc and take part in the evolution. On the other hand, if the halo is rigid, it hinders or quenches bar formation, as expected. Comparison of two simulations which are identical in everything, except that the halo is live in the first one and rigid in the second one, leads me to suggest that the halo response can help the bar grow. Following the orbits of the stars in the halo, I find that a considerable fraction of the halo particles are in resonance with the bar. The halo may thus take angular momentum from the bar and stimulate its growth. I finally discuss whether and how the results of the N-body simulations can be applied to real galaxies.

I present N-body simulations of isolated and interacting galaxies, made on GRAPE machines. In particular I discuss the formation and evolution of N-body bars and compare their properties to those of bars in early-type and late-type galactic discs. I argue that the halo can help the bar grow, contrary to previous beliefs, by taking positive angular momentum from it via its resonant stars. I then focus on the interaction and subsequent merging of a barred disc galaxy with a spheroidal satellite. The evolution depends strongly on the mass (density) of the satellite and may lead to its destruction or to the destruction of the bar.

We discuss the morphology, photometry and kinematics of the bars which have formed in three N-body simulations. These have initially the same disc and the same halo-to-disc mass ratio, but their haloes have very different central concentrations. The third model includes a bulge. The bar in the model with the centrally concentrated halo (model MH) is much stronger, longer and thinner than the bar in the model with the less centrally concentrated halo (model MD). Its shape, when viewed side-on, evolves from boxy to peanut and then to X-shaped, as opposed to that of model MD, which stays boxy. The projected density profiles obtained from cuts along the bar major axis, both for the face-on and the edge-on views, show a flat part, as opposed to those of model MD which are falling rapidly. A Fourier analysis of the face-on density distribution of model MH shows very large m = 2, 4, 6 and 8 components. Contrary to this for model MD the components m = 6 and 8 are negligible. The velocity field of model MH shows strong deviations from axial symmetry, and in particular has wavy isovelocities near the end of the bar when viewed along the bar minor axis. When viewed edge-on, it shows cylindrical rotation, which the MD model does not. The properties of the bar of the model with a bulge and a non-centrally concentrated halo (MDB) are intermediate between those of the bar of the other two models. All three models exhibit a lot of inflow of the disc material during their evolution, so that by the end of the simulations the disc dominates over the halo in the inner parts, even for model MH, for which the halo and disc contributions were initially comparable in that region.

In this paper we measure the two-body relaxation time from the angular deflection of test particles launched in a rigid configuration of field particles. We find that centrally concentrated configurations have relaxation times that can be shorter than those of the corresponding homogeneous distributions by an order of magnitude or more. For homogeneous distributions we confirm that the relaxation time is proportional to the number of particles. On the other hand centrally concentrated configurations have a much shallower dependence, particularly for small values of the softening. The relaxation time increases with the inter-particle velocities and with softening. The latter dependence is not very strong, of the order of a factor of two when the softening is increased by an order of magnitude. Finally we show that relaxation times are the same on GRAPE-3 and GRAPE-4, dedicated computer boards with limited and high precision respectively.

We present N-body simulations of galaxy groups embedded in a common halo of matter. We study the influence of the different initial conditions upon the evolution of the group and show that denser configurations evolve faster, as expected. We then concentrate on the influence of the initial radial density profile of the common halo and of the galaxy distribution. We select two kinds of density distributions, a singular profile (modeled by a Hernquist distribution) and a profile with a flat core (modeled by a Plummer sphere). In all cases we witness the formation of a central massive object due to mergings of individual galaxies and to accretion of stripped material, but both its formation history and its properties depend heavily on the initial distribution. In Hernquist models the formation is due to a ``burst'' of mergings in the inner parts, due to the large initial concentration of galaxies in the center. The merging rate is much slower in the initial phases of the evolution of a Plummer distribution, where the contribution of accretion to the formation of the central object is much more important. The central objects formed within Plummer distributions have projected density profiles which are not in agreement with the radial profiles of observed brightest cluster members, unless the percentage of mass in the common halo is small. On the contrary the central object formed in initially cusped models has projected radial profiles in very good agreement with those of brightest cluster members, sometimes also showing luminosity excess over the $r^{1/4}$ law in the outer parts, as is observed in cD galaxies.

I briefly review the properties of the gas flow in and around the region of the bar in a disc galaxy and discuss the corresponding inflow and the loci of star formation. I then review the flow of gas in barred galaxies which have an additional secondary bar. Finally I discuss the signatures of bars in edge-on galaxies.

In N-body simulations the force calculated between particles representing a given mass distribution is usually softened, to diminish the effect of graininess. In this paper we study the effect of such a smoothing, with the aim of finding an optimal value of the softening parameter. As already shown by Merritt (1996), for too small a softening the estimates of the forces will be too noisy, while for too large a softening the force estimates are systematically misrepresented. In between there is an optimal softening, for which the forces in the configuration approach best the true forces. The value of this optimal softening depends both on the mass distribution and on the number of particles used to represent it. For higher number of particles the optimal softening is smaller. More concentrated mass distributions necessitate smaller softening, but the softened forces are never as good an approximation of the true forces as for not centrally concentrated configurations. We give good estimates of the optimal softening for homogeneous spheres, Plummer spheres, and Dehnen spheres. We also give a rough estimate of this quantity for other mass distributions, based on the harmonic mean distance to the $k$th neighbour ($k$ = 1, .., 12), the mean being taken over all particles in the configuration. Comparing homogeneous Ferrers ellipsoids of different shapes we show that the axial ratios do not influence the value of the optimal softening. Finally we compare two different types of softening, a spline softening (Hernquist & Katz 1989) and a generalisation of the standard Plummer softening to higher values of the exponent. We find that the spline softening fares roughly as well as the higher powers of the power-law softening and both give a better representation of the forces than the standard Plummer softening.

We use numerical simulations to investigate the cusp at the centre of elliptical galaxies, due to the slow growth of a super-massive black hole. We study this problem for axisymmetric models of galaxies, with or without rotation. The numerical simulations are based on the `Perturbation Particles' method, and use GRAPEs to compute the force due to the cusp. We study how the density cusp is affected by the initial flattening of the model, as well as the role played by initial rotation. The logarithmic slope of the density cusp is found to be very much insensitive to flattening; as a consequence, we deduce that tangential velocity anisotropy -which supports the flattening- is also of little influence on the final cusp. We investigate via two different kinds of rotating models the efficiency with which a rotation velocity component builds within the cusp. A cusp in rotation develops only for models where a net rotation component is initially present at high energy levels. The eventual observation of a central rotational velocity peak in E galaxies has therefore some implications for the galaxy dynamical history.

In this paper I focus on three topics related to the dynamical evolution of small galaxy groups, for which the input of N-body simulations has been decisive. These are the merging rates in compact groups, the properties of remnants of multiple mergers, and the evolution of disc galaxies surrounded by one or more satellites. The short dynamical times of compact groups make it difficult to understand why such groups are observed at all. N-body simulations have pointed out two possible classes of solutions to this problem. The first one proposes that there is on-going formation of compact groups, or that the longevity of the group is due to secondary infall. For the second class of solutions the longevity of compact groups is due either to their specific initial conditions, or to a massive common halo, encompassing the whole group. I discuss here these alternatives, together with their respective advantages and disadvantages. I then turn to the structure of remnants of multiple mergers and compare the results of N-body simulations with the properties of observed elliptical galaxies. Finally I discuss the dynamical evolution of a disc galaxy surrounded by one or more spherical satellites.

I first comment on some recent advances in computing equipment for CPU-intensive numerical simulations, and on possible developments in the near future. I then discuss some particularly important and yet unsolved problems in dynamics and evolution of galaxies, on which analytical, numerical and observational effort should be focused.

We develop diagnostics based on gas kinematics to identify the presence of a bar in an edge-on spiral galaxy and determine its orientation. We use position-velocity diagrams (PVDs) obtained by projecting edge-on two-dimensional hydrodynamical simulations of the gas flow in a barred galaxy potential. We show that when a nuclear spiral is formed, the presence of a gap in the PVDs, between the signature of the nuclear spiral and that of the outer parts of the disk, reliably indicates the presence of a bar. This gap is due to the presence of shocks and inflows in the simulations, leading to a depletion of the gas in the outer bar region. If no nuclear spiral signature is present in a PVD, only indirect arguments can be used to argue for the presence of a bar. The shape of the signature of the nuclear spiral, and to a lesser extent that of the outer bar region, allows to determine the orientation of the bar with respect to the line-of-sight. The presence of dust can also help to discriminate between viewing angles on either side of the bar. Simulations covering a large fraction of parameter space constrain the bar properties and mass distribution of observed galaxies. The strongest constraint comes from the presence or absence of the signature of a nuclear spiral in the PVD.

Using an N-body+SPH code we have performed numerical simulations to investigate the dynamical effects of an interaction between an initially barred galaxy and a small spherical companion. In the models described here the small companion passes through the disk of the larger galaxy perpendicular to its plane. The impact positions and times are varied with respect to the evolutionary phase of the bar and disc. The interactions produce expanding ring structures, offset bars, spokes, and other asymmetries in the stars and gas. They also affect the strength and pattern speed of the bar.

We develop diagnostics to detect the presence and orientation of a bar in an edge-on disk, using its kinematical signature in the position-velocity diagram (PVD) of a spiral galaxy observed edge-on. Using a well-studied barred spiral galaxy mass model, we briefly review the orbital properties of two-dimensional non-axisymmetric disks and identify the main families of periodic orbits. We use those families as building blocks to model real galaxies and calculate the PVDs obtained for various realistic combinations of periodic orbit families and for a number of viewing angles with respect to the bar. We show that the global structure of the PVD is a reliable bar diagnostic in edge-on disks. Specifically, the presence of a gap between the signatures of the families of periodic orbits in the PVD follows directly from the non-homogeneous distribution of the orbits in a barred galaxy. Similarly, material in the two so-called forbidden quadrants of the PVD results from the elongated shape of the orbits. We show how the shape of the signatures of the dominant x1 and x2 families of periodic orbits in the PVD can be used efficiently to determine the viewing angle with respect to the bar and, to a lesser extent, to constrain the mass distribution of an observed galaxy. We also address the limitations of the models when interpreting observational data.

Prompted by work on the buckling instability in barred spiral galaxies, much effort has been devoted lately to the study of boxy/peanut-shaped (B/PS) bulges. Here, we present new bar diagnostics for edge-on spiral galaxies based on periodic orbits calculations and hydrodynamical simulations. Both approaches provide reliable ways to identify bars and their orientations in edge-on systems. We also present the results of an observational search for bars in a large sample of edge-on spirals with and without B/PS bulges. We show that most B/PS bulges are due to the presence of a thick bar viewed edge-on while only a few may be due to accretion. This strongly supports the bar-buckling mechanism for the formation of B/PS bulges.

Disc galaxies can be substantially modified by close encounters and mergers, since their discs are very responsive components. Close interactions can be held responsible for the formation of bridges and tails, as well as for the formation of some bars, asymmetries and grand design spirals. Bound clumps can form in the tails, due to self-gravity, and could evolve to dwarf galaxies. Off-centerings and asymmetries in the central parts of barred galaxies can be made by off-centered and/or oblique impacts of sufficiently massive and compact companions. Similar impacts, but preferably centered, on non-barred galaxies can form ring galaxies. Companions on initially near-circular orbits can also cause changes to the target disc as they spiral gradually inwards. Low density companions are disrupted before reaching the center of the target and their debris form a thick disc. On the other hand most of the mass of the high density companions reaches the center, where it may form a bulge, thus entailing evolution along the Hubble sequence. Such companions thicken and expand the target disc and may also destroy bars in it. If their initial orbital plane is at an angle to the plane of the disc of the target, they can cause the latter to tilt substantially, depending on their mass and initial inclination.

I review the Marseille GRAPE systems and the N-body simulations done with them. I first describe briefly the available hardware and software, their possibilities and their limitations. I then describe work done on interacting galaxies and groups of galaxies. This includes simulations of the formation of ring galaxies, simulations of bar destruction by massive compact satellites, of merging in compact groups and of the formation of brightest members in clusters of galaxies.

In this paper we propose an improved method for calculating Henon's stability parameter, which is based on the differential of the Poincare map using the first variational equation. We show that this method is very accurate and give some examples where it gives correct results, while the previous method could not cope.

The performance and accuracy of a GRAPE-3 system for collisionless N-body simulations is discussed. After a description of the hardware configurations available to us at Marseille, and the usefulness of on-line analysis, we concentrate on the performance and accuracy of the direct summation and tree code software. For the former we discuss the sources of round-off errors. The standard Barnes-Hut tree code has been modified for GRAPE-3, by dividing particles into blocks and performing the tree traversal for the whole block. The forces are then calculated by direct summation over the whole interaction list. The performance of the tree code depends on the number of particles in the block, the optimum number depending on the speed of the front end and the number of boards. The time per step scales with the number of particles, and decreases only weakly as the tolerance increases. We discuss the effect of the front end and compare performances with those of general purpose computers. The accuracy of both direct summation and the tree code is discussed as function of the number of particles and the softening. The increased role of direct summation in the force calculation enhances the accuracy compared to the standard tree code. Finally, we follow the evolution of an isolated barred galaxy using different hardware and software in order to assess the reliability of our results. We find excellent agreement between the pattern speed of the bar in direct summation simulations run on the high precision GRAPE-4 machines and on our GRAPE-3 system. The agreement with the tree code is very good for tolerance values smaller than about 1.0. We conclude that GRAPE-3 systems are well suited for collisionless simulations, in particular those of galaxies.

We present N-body simulations of groups of galaxies with a number of very different initial conditions. These include spherical isotropic, nonspherical anisotropic collapses and virialised spherical systems. In all cases but one the merging instability leads to the formation of a giant central galaxy in the center of the group. The initial conditions of the exception are such that no galaxies are present in the central part of the group. Thus some central seed of material is necessary to trigger the formation of a giant central galaxy. We concentrate on the properties of these giant central galaxies. Spherical virialised systems give rise to relatively round and isotropic systems, while aspherical initial conditions give rise to triaxial objects with anisotropic velocity dispersion tensors. In the latter cases the orientation of the resulting central galaxy is well correlated with that of the initial cluster. We compare the projected properties of the objects formed with the properties of real brightest cluster member galaxies. The surface density profiles are in good agreement with the observed surface brightness profiles. In the case of extended virialised groups the projected properties of the giant central galaxy are the same as the properties of cD galaxies. These include a halo of luminous material and a nearly flat velocity dispersion profile.

We analyze the morphology and dynamics of the galaxies of the Hickson 96 compact group by means of deep CCD images in the B, V and R photometric bands and long-slit spectroscopy. The four galaxies of this spiral rich group show signs of gravitational interaction. Two long tails come out from the contact region of the close pair H96ac. It is also there where both galaxies, bi-symmetrical in their inner parts, loose one of their arms. Moreover, both the photometric and kinematical center of H96a are displaced relative to the center of the disk. H96b, a giant elliptical galaxy, shows significant deviations from a $r^{1/4}$ law at the inner parts and small rotation along its major axis. We find at the centre of the galaxy an elongated component which is kinematically decoupled. A wide faint plume seems to emerge from this galaxy. H96d, the smallest galaxy of the group, seems to be influenced by the bigger members, as indicated by three prominent knots of recent bursts of star formation in its blue disk, and by hints of optical bridges joining H96d with H96a and b. Finally the group has a low velocity dispersion (160 km/s) and is well isolated - no galaxies with comparable magnitude to H96a or H96b are found in its neighborhood. All these results lead us to conclude that Hickson 96 constitutes a real physical system.

We use N-body simulations to study the formation of rings in a disc galaxy by the impact of a small spherical companion. Both barred and nonbarred target discs are considered. We discuss the effect of the properties of the target disc (distribution of mass in the disc, velocity dispersion, etc.) as well as of the mass and orbit of the companion on the properties of the rings, such as their amplitude, width, shape, expansion velocity and lifetime. In particular the amplitude, width, lifetime and expansion velocity of the first ring increase considerably with companion mass, and so does the expansion velocity of the particles in it and the total extent of the disc after the interaction. We also discuss the formation and evolution of spokes and show that they can be caused by companions of relatively small mass. In our three examples of oblique impacts on barred target galaxies we note important transient displacements of the bar, as well as changes of its pattern speed and size. An asymmetric pseudoring is formed in each case, and during the first stages of its evolution the bar forms part of it.

We discuss the merging rates in compact groups of 5 identical elliptical galaxies. All groups have the same mass and binding energy. We consider both cases with individual halos and cases where the halo is common to all galaxies and enveloping the whole group. In the latter situation the merging rate is slower if the halo is more massive. The mass of individual halos has little influence on the merging rates, due to the fact that all galaxies in our simulations have the same mass, and so the more extended ones have a smaller velocity dispersion. Groups with individual halos merge faster than groups with common halos if the configuration is centrally concentrated, like a King distribution of index 10. On the other hand for less concentrated configurations the merging is initially faster for individual halo cases, and slower after part of the group has merged. In cases with common halo, centrally concentrated configurations merge faster for high halo-to-total mass ratios and slower for low halo-to-total mass ratios. Groups whose virial ratio is initially less than one merge faster, while groups that have initially cylindrical rotation merge slower than groups starting in virial equilibrium. In order to test how long a virialised group can survive before merging we followed the evolution of a group with a high halo-to-total mass ratio and a density distribution with very little central concentration. We find that the first merging occurred only after a large number of crossing times, which with areasonable calibration should be larger than a Hubble time. Hence, at least for appropriate initial conditions, the longevity of compact groups is not necessarily a problem, which is an alternative explanation to why we observe so many compact groups despite the fact that their lifetimes seem short.

In this paper we test the performance of explicit simulations of groups of galaxies, (i.e. simulations in which each galaxy is treated as a mass point and the physics of the interactions is modelled by specific analand the physics of the interactions is modelled by specific analytical prescriptions for merging conditions), by comparing them with fully self-consistent simulations starting from identical initial conditions. The quality of the explicit simulations is very unequal. For somThe quality of the explicit simulations is very unequal. For some prescriptions the results are complete disagreement with the self-consistent simulations. The inclusion of other dynamical effects like dynamical friction gives, in some cases, better agreeement.like dynamical friction gives, in some cases, better agreeement. We also propose a new merging criterion, which, combined with dynamical friction, gives much better agreement with self-consistent simulations in a variety of initial conditions, but even this criterion has a limited ranof initial conditions, but even this criterion has a limited range of applicability.